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.     

Analytical solutions of one‐dimensional infiltration before and after ponding

This paper is a continuation of previous research, which obtained a convenient solution for arbitrary surface fluxes before ponding. By means of Fourier Transformation this has been extended to derive analytical solutions of a linearized Richards' equation for arbitrary input fluxes after surface saturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Effects of Different Diets on the Dietary Attractability and Selectivity of Chinese Shrimp, Fenneropenaeus chinensis

Attractabilities of different diets and dietary selectivity of Chinese shrimp, Fenneropenaeus chinensis were studied through behavior observation and feeding experiment, respectively. The five diets used in the experiment are: Fish Flesh (FF), Shrimp Flesh (SF), Clam Foot (CF), Polychaete Worm (PW), and Formulated Diet (FD). No significant differences of attractability exist between any two different diets when every two natural diets or all five diets are provided simultaneously. On the other hand, significant differences of attractability exist between FD and every single natural diet when they are provided simultaneously. Results of behavioral observation indicate that natural diets are more attractive than FD. In feeding experiment, Chinese shrimp has distinct selectivity on different diets. It positively selects CF and PW, negatively selects FF and SF, and excludes FD absolutely. The results of the present studies indicate that the dietary selectivity of shrimp was based not only on the attractabilities of the diets, but also on the responses such as growth and food conversion.     

Peak metamorphic temperatures in type 6 ordinary chondrites: An evaluation of pyroxene and plagioclase geothermometry

Abstract— Quantifying the peak temperatures achieved during metamorphism is critical for understanding the thermal histories of ordinary chondrite parent bodies. Various geothermometers have been used to estimate equilibration temperatures for chondrites of the highest metamorphic grade (type 6), but results are inconsistent and span hundreds of degrees. Because different geothermometers and calibration models were used with different meteorites, it is unclear whether variations in peak temperatures represent actual ranges of metamorphic conditions within type 6 chondrites or differences in model calibrations. We addressed this problem by performing twopyroxene geothermometry, using QUILF95, on the same type 6 chondrites for which peak temperatures were estimated using the plagioclase geothermometer (Nakamuta and Motomura 1999). We also calculated temperatures for published pyroxene analyses from other type 6 H, L, and LL chondrites to determine the most representative peak metamorphic temperatures for ordinary chondrites. Pyroxenes record a narrow, overlapping range of temperatures in H6 (865–926 °C), L6 (812–934 °C), and LL6 (874–945 °C) chondrites. Plagioclase temperature estimates are 96–179 °C lower than pyroxenes in the same type 6 meteorites. Plagioclase estimates may not reflect peak metamorphic temperatures because chondrule glass probably recrystallized to plagioclase prior to reaching the metamorphic peak. The average temperature for H, L, and LL chondrites (～900 °C), which agrees with previously published oxygen isotope geothermometry, is at least 50 °C lower than the peak temperatures used in current asteroid thermal evolution models. This difference may require minor adjustments to thermal model calculations.     

中国东部第四纪自然环境演化及庐山“冰碛层”堆积时期的自然环境研究

根据第四纪地层内孢粉组合和动物群等新资料恢复了我国东部各大区域第四纪各个时期的自然环境特征,其温暖期与干冷期的交替与深海沉积O~(18)温度变化趋势大致对应。庐山三次“冰期冰碛层”内孢粉组合表明均为温暖湿润针阔叶混交林景观,而不是寒冷的冰川环境。     

Numerical simulation of January 28 cold air outbreak during GALE Part I: The model and sensitivity tests of turbulence closures

A mesoscale planetary boundary layer (PBL) numerical model has been developed to study airflow over complex topography. Turbulence closures using the turbulent kinetic energy (TKE) and dissipation () equations are investigated in combination with the level 2.5 scheme of Mellor and Yamada (1982) to determine eddy diffusivities for momentum and heat. This modified E- closure is simpler than the level 3 one which requires more prognostic equations for moist turbulent transport.One-dimensional (1-D) model results show that the PBL mean flows under various stability conditions are not significantly sensitive to the modified Blackadar and Kolmogorov eddy mixing-length formulations used in this E- model, although the latter yields excessively large mixing lengths in the entrainment region of the upper PBL. Eddy mixing lengths in the Kolmogorov-type formulation can be better defined by introducing background dissipation. Using the same prognostic TKE equation, the 1-D model results are not significantly affected by different diagnostic formulations in the closures. The simulated results compare well with large-eddy simulations and those obtained using higher-order closure schemes including the level 3 one. The results are found to be insensitive to eddy Prandtl number, in contrast to the 2-D model results (see Part II).     

TL dating of loess in China

TL ages of 56 loess samples collected from Xinjiang and the middle reaches of the Yellow River were determined by using fine-grained quartz and man-made light source bleaching techniques.The results indicate that :1.The well-preserved loess sections in Xinjiang began depositing 400000 a ago.The loess/paleosol series may be comparable with that observed in the middle reaches of the Yellow River but the depositing time was probably delayed by 30000-50000a.2.The bottom boundary age of the Malan loess varies from one place to another.The TL dates can be divided into two group:70000-130000 a and 30000-50000 a.The TL data are con-sistent with the appearance of two warmer climate periods since the Pleistocene,indicating obvious regional differences in stratigraphic development of loess,i.e.,paleosol occurred in the warm climate period in the same location,but loess was deposited in other locations.The initial age of loess deposition tends to become older from west to east and from south to north .In contrast ,the cold cli-matic belt showed an opposite trend.The TL ages of loess provide strong evidence for the Quaternary stratigraphy throughout geological times and the shift tendency of climatic belt.3.The depositing rate of the Malan loess along the middle reaches of the Yellow River varies from place to place but the TL age shows a linear relationship with the depositing thickness.The de-positing rate of loess in the north is higher than that in the south,suggesting that loess materials are derived from desert areas of northern China.     

Neutrino mass,luminosity variation,and spectrum of SN1987A

In this paper, the neutrino mass has been determined from SN1987a observation in a manner that the simultaneity of neutrino emission is not regarded as the starting point, but is itself defined through the analysis by Monte-Carlo simulation. The result is that the neutrino mass lies in 3–4 eV, possiblym _v 3.6 eV. Neutrino luminosity variation and neutrino spectrum are also obtained. Comparison with theories gives further support to the mass determination, and also predicts the mass of progenitor star of SN1987a to be in the range of 12–25M .The project supported by the National Natural Science Foundation of China.     

Palaeoenvironmental significance of Late Permian palaeosols in the South‐Eastern Iberian Ranges,Spain

The Late Permian (Wuchiapingian) Alcotas Formation in the SE Iberian Ranges consists of one red alluvial succession where abundant soil profiles developed. Detailed petrographical and sedimentological studies in seven sections of the Alcotas Formation allow six different types of palaeosols, with distinctive characteristics and different palaeogeographical distribution, to be distinguished throughout the South‐eastern Iberian Basin. These characteristics are, in turn, related to topographic, climatic and tectonic controls. The vertical distribution of the palaeosols is used to differentiate the formation in three parts from bottom to top showing both drastic and gradual vertical upwards palaeoenvironmental changes in the sections. Reconstruction of palaeoenvironmental conditions based on palaeosols provides evidence for understanding the events that occurred during the Late Permian, some few millions of years before the well‐known Permian‐Triassic global crisis.     

Methane‐related microbial gypsum calcitization in stromatolites of a marine evaporative setting (Münder Formation,Upper Jurassic,Hils Syncline,north Germany)

Fossil stromatolites may reveal information about their hydrochemical palaeoenvironment, provided that assignment to a specific microbial community and a corresponding biogeochemical mechanism of formation can be made. Tithonian stromatolites of the Münder Formation at Thüste, north Germany, have traditionally been considered as formed by intertidal cyanobacterial communities. However, thin sections of the stromatolites show elongated angular traces of former gypsum crystals in a dense arrangement, but no algal or cyanobacterial filament traces. Moreover, high Fe²⁺ and Mn²⁺ contents, oxygen‐isotope and sulphur‐isotope ratios of carbonate‐bound sulphates, and sulphurized hydrocarbon biomarkers of the stromatolitic carbonate indicate that CaCO₃ precipitation occurred near the oxic–anoxic interface as a result of intensive bacterial sulphur cycling rather than photosynthetic activity. Furthermore, anaerobic oxidation of methane by Archaea may have driven CaCO₃ precipitation in deeper parts of the biofilm community, as reflected by high concentrations of squalane with a strongly negative δ¹³C in conjunction with evaporite pseudomorphs showing extremely low δ¹³C_Carb ratios. Consequently, the Thüste stromatolites are now interpreted as having initially formed by gypsum impregnation of biofilms. Subsequently, early Mg‐calcitic calcitization within the biofilms occurred because of combined bacterial iron, manganese and sulphate reduction, with an increasing contribution of anaerobic oxidation of methane with depth. This model plausibly explains the prominent preservation of signals derived from oxygen‐independent metabolic pathways, whereas virtually no geochemical record exists for an aerobic community that may, nevertheless, have prevailed at the stromatolite surface. Photic‐zone stromatolites with a prominent signal of anaerobic oxidation of methane may be common in, and indicative of, oxygen‐depleted sulphate‐bearing environments with high rates of methane production, conditions that possibly were fulfilled at the Archaean to Proterozoic transition.