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°,与浊流主体流向基本一致。还认为中三叠世时右江盆地沉积物源来自桂东和桂东南蚀源区,海水由其南的越南境内进入。     

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.     

Timing of UHP metamorphism in the Hong’an area,western Dabie Mountains,China: evidence from zircon U–Pb age,trace element and Hf isotope composition

The Hong’an area (western Dabie Mountains) is the westernmost terrane in the Qinling-Dabie-Sulu orogen that preserves UHP eclogites. The ages of the UHP metamorphism have not been well constrained, and thus hinder our understanding of the tectonic evolution of this area. LA-ICPMS U–Pb age, trace element and Hf isotope compositions of zircons of a granitic gneiss and an eclogite from the Xinxian UHP unit in the Hong’an area were analyzed to constrain the age of the UHP metamorphism. Most zircons are unzoned or show sector zoning. They have low trace element concentrations, without significant negative Eu anomalies. These metamorphic zircons can be further subdivided into two groups according to their U–Pb ages, and trace element and Lu–Hf isotope compositions. One group with an average age of 239 ± 2 Ma show relatively high and variable HREE contents (527 ≥ Lu_N ≥ 14) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.00008–0.000931), indicating their growth prior to a great deal of garnet growth in the late stage of continental subduction. The other group yields an average age of 227 ± 2 Ma, and shows consistent low HREE contents and ¹⁷⁶Lu/¹⁷⁷Hf ratios, suggesting their growth with concurrent garnet crystallization and/or recrystallization. These two groups of age are taken as recording the time of prograde HP to UHP and retrograde UHP–HP stages, respectively. A few cores have high Th/U ratios, high trace element contents, and a clear negative Eu anomaly. These features support a magmatic origin of these zircon cores. The upper intercept ages of 771 ± 86 and 752 ± 70 Ma for the granitic gneiss and eclogite, respectively, indicate that their protoliths probably formed as a bimodal suite in rifting zones in the northern margin of the Yangtze Block. Young Hf model ages (T _DM1) of magmatic cores indicate juvenile (mantle-derived) materials were involved in their protolith formation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

北方泥炭沼泽中酚类物质释放的热敏感性研究:以英国威尔士地区雨养和矿养泥炭沼泽为例

利用高精度多梯度热培养系统,对英国威尔士地区两类典型的北方泥炭沼泽中酚类物质和溶解有机碳释放的热敏感性进行了为期1年的调查研究.研究结果表明,雨养泥炭沼泽(Bog)酚类物质释放的热敏感性Q10(总酚)月变化为0.92～1.57,而矿养泥炭沼泽(Fen)的月变化范围为0.93～1.30.酚类释放的热敏感性与土壤温度大致呈正相关关系.此外,土壤温度与溶解有机碳释放的热敏感性Q10(DOC)也呈正相关关系.总体上,雨养泥炭沼泽的热敏感性比矿养泥炭沼泽略高.除温度外,水文条件和植被等多种环境因子对酚类物质和溶解有机碳释放的热敏感性也有较大影响.研究结果还表明酚碳热敏感性比值Q10o(总酚)/Q10(DOC)相对稳定,土壤温度的变化对其影响有限.     

Comparative analyses of East Texas forest cover maps generated from Landsat and AVHRR data

Comparing satellite data derived map products are affected by differences in data characteristics, image acquisition dates, processing techniques, and classification schemes used for assigning pixels to a thematic class. By comparing two forest maps generated from Landsat Enhanced Thematic Mapper Plus (ETM+) and Advanced Very High Resolution Radiometer (AVHRR) images acquired on the same day, and processed using identical classification scheme and methods these differences were minimized. The ETM+ derived map had higher classification accuracy values and more precise area estimates than the AVHRR derived map. In the ETM+ derived map, 87 of the 599 verification data were misclassified, whereas in the AVHRR derived map, 155 of the 469 verification data were misclassified. Detailed error analyses by land cover class revealed that a land use based definition of forest accounted for 74% (64 out of 87) and 57% (89 out of 155) of the classification errors in ETM+ and AVHRR derived maps, respectively.     

圆柱螺线法测斜计算中的数值方法

圆柱螺线法是井眼轨迹计算中最常使用的方法之一,在实钻轨迹监控、中靶分析和预测、井身质量评价等实际问题中都有重要的应用.但是当井段2个端点的井斜角之差或者方位角之差等于零时,圆柱螺线法坐标计算公式需要进行特殊处理才能正常计算;在井斜角之差或者方位角之差很小时,坐标计算会出现较大误差.对圆柱螺线法坐标计算公式进行了改进,用一组计算公式统一处理井斜角之差或者方位角之差等于零和不等于零的情况,提高了坐标的数值计算的精度和稳定性.对圆柱螺线型井眼轨迹的平均井眼曲率的数值计算进行了研究,提出了使用高斯型数值积分公式计     

基于AGA的SVM需水预测模型研究

需水预测是一个由城市人口、工业水平、社会经济水平共同作用的多因素、多层次的复杂非线性系统.其结果将直接影响受区域水资源承载力约束的产业结构、布局形态等决策.作为一种集中参数预报方法,支持向量机方法具有对未来样本的较好的泛化性能,对于这类资料缺乏、系统结构尚欠清晰的问题可以取得较好的模拟和预测结果.基于此,本文将支持向量机方法引入需水预测领域,建立了需水预测支持向量机模型.同时,本文将加速遗传算法和支持向量机方法耦合起来,构造了支持向量机模型参数的自适应优化算法.模型在珠海市的应用实例表明:与简单遗传算法比较,AGA的模型参数寻优效率更高;与BP神经网络模型相比,SVM模型较好地解决了小样本、经验性等问题,并取得了较高的预测精度.