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.     

多接收器电感耦合等离子体质谱法测定镉标准溶液中镉的同位素组成

建立了多接收器电感耦合等离子体质谱仪对Cd标准溶液同位素组成进行分析的实验方法。仪器的质量分馏校正采用Standard-Sample Bracketing法。实验结果用δ114/110Cd来表达。在此研究基础上,以SPEX Cd标准溶液为参考标准,对国外4种Cd标准溶液进行了测定。结果表明,实验测定的精度在0.07‰~0.13‰(δ114/110Cd),与目前文献报道的结果具有相似的精度。以最新SPEX Cd标准样品(δ114/110Cd=0)为基准,计算的δ114/110CdJMC、δ114/110CdSPEX-1、δ114/110CdBAM1012和δ114/110CdM櫣nster的值分别为0.55‰、0.56‰、-0.65‰和5.14‰,说明不同批次SPEX标准溶液的Cd同位素组成是明显不同的,最新的SPEX Cd标准溶液与SPEX-1 Cd的δ114/110Cd值存在着0.56‰的差别。将以SPEX Cd为参考标准的δ114/110CdBAM1012和δ114/110CdM櫣nster转化为以SPEX-1或JMC为参考标准后,得到的结果与文献报道的结果在误差范围内一致。     

Zircon U–Pb and Hf isotopic study of gneissic rocks from the Chinese Altai: Progressive accretionary history in the early to middle Palaeozoic

Gneissic rocks in the Chinese Altai Mountains have been interpreted as either Paleozoic metasedimentary rocks or Precambrian basement. This study reports geochemical and geochronological data for banded paragneisses and associated gneissic granitoids collected along a NE–SW traverse in the northwestern Chinese Altai. Petrological and geochemical data suggest that the protoliths of the banded gneisses were possibly immature sediments with significant volcanic input and that the gneissic granitoids were derived from I-type granites formed in a subduction environment. Three types of morphological features can be recognized in zircons from the banded gneisses and are interpreted to correlate with different sources. Zircons from five samples of banded paragneiss cluster predominantly between 466 and 528 Ma, some give Neoproterozoic ages, and a few yield discordant Paleoproterozoic to Archean ages. Zircon Hf isotopic compositions indicate that both juvenile/mantle and crust materials were involved in the generation of the source rocks from which these zircons were derived. In contrast, zircons occur ubiquitously as elongated euhedral prismatic crystals in the four samples of the gneissic granitoids, and define single populations for each sample with mean ages between 380 and 453 Ma. The general absence of Precambrian inheritance and positive zircon ?Hf values for these granitoids suggest insignificant crustal contribution to the generation of the precursor magmas. Our data can be interpreted in terms of a progressive accretionary history in early to middle Palaeozoic times, and the Chinese Altai may possibly represent a magmatic arc built on a continental margin dominated by Neoproterozoic rocks.     

Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Biaxial test simulations using a packing of polygonal particles

The mechanical response of cohesionless granular materials under monotonic loading is studied by performing molecular dynamic simulations. The diversity of shapes of soil grains is modelled by using randomly generated convex polygons as granular particles. Results of the biaxial test obtained for dense and loose media show that samples achieve the same void ratio at large strains independent of their initial density state. This limit state resembles the so‐called critical state of soil mechanics, except for some stress fluctuations, which remain for large deformations. These fluctuations are studied at the micro‐mechanical level, by following the evolution of the co‐ordination number, force chains and the fraction of the sliding contacts of the sample. Copyright © 2007 John Wiley & Sons, Ltd.