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.     

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

Relationship between footwall composition,crustal contamination,and fluid–rock interaction in the Platreef,Bushveld Complex,South Africa

In the northern limb of the 2.06-Ga Bushveld Complex, the Platreef is a platinum group elements (PGE)-, Cu-, and Ni-mineralized zone of pyroxenite that developed at the intrusion margin. From north to south, the footwall rocks of the Platreef change from Archaean granite to dolomite, hornfels, and quartzite. Where the footwall is granite, the Sr-isotope system is more strongly perturbed than where the footwall is Sr-poor dolomite, in which samples show an approximate isochron relationship. The Nd-isotope system for samples of pyroxenite and hanging wall norite shows an approximate isochron relationship with an implied age of 2.17 ± 0.2 Ga and initial Nd-isotope ratio of 0.5095. Assuming an age of 2.06 Ga, the ɛNd values range from −6.2 to −9.6 (ave. −7.8, n = 17) and on average are slightly more negative than the Main Zone of the Bushveld. These data are consistent with local contamination of an already contaminated magma of Main Zone composition. The similarity in isotope composition between the Platreef pyroxenites and the hanging wall norites suggests a common origin. Where the country rock is dolomite, the Platreef has generally higher plagioclase and pyroxene δ ¹⁸O values, and this indicates assimilation of the immediate footwall. Throughout the Platreef, there is considerable petrographic evidence for sub-solidus interaction with fluids, and the Δ _{plagioclase–pyroxene} values range from −2 to +6, which indicates interaction at both high and low temperatures. Whole-rock and mineral δD values suggest that the Platreef interacted with both magmatic and meteoric water, and the lack of disturbance to the Sr-isotope system suggests that fluid–rock interaction took place soon after emplacement. Where the footwall is granite, less negative δD values suggest a greater involvement of meteoric water. Consistently higher values of Δ _{plagioclase–pyroxene} in the Platreef pyroxenites and hanging wall norites in contact with dolomite suggest prolonged interaction with CO₂-rich fluid derived from decarbonation of the footwall rocks. The overprint of post crystallization fluid–rock interaction is the probable cause of the previously documented lack of correlation between PGE and sulfide content on the small scale. The Platreef in contact with dolomite is the focus of the highest PGE grades, and this suggests that dolomite contamination played a role in PGE concentration and deposition, but the exact link remains obscure. It is a possibility that the CO₂ produced by decarbonation of assimilated dolomite enhanced the process of PGE scavenging by sulfide precipitation.     

Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

In situ strength measurements on natural upper-mantle minerals

Using in situ strength measurements at pressures up to 10 GPa and at room temperature, 400, 600, and 700°C, we examined rheological properties of olivine, orthopyroxene, and chromian-spinel contained in a mantle-derived xenolith. Mineral strengths were estimated using widths of X-ray diffraction peaks as a function of pressure, temperature, and time. Differential stresses of all minerals increase with increasing pressure, but they decrease with increasing temperature because of elastic strain on compression and stress relaxation during heating. During compression at room temperature, all minerals deform plastically at differential stress of 4–6 GPa. During subsequent heating, thermally induced yielding is observed in olivine at 600°C. Neither orthopyroxene nor spinel shows complete stress relaxation, but both retain some stress even at 700°C. The strength of the minerals decreases in the order of chromian-spinel ≈ orthopyroxene > olivine for these conditions. This order of strength is consistent with the residual pressure of fluid inclusions in mantle xenoliths.     

Studies on the distribution and antibacterial properties of Streptomyces in the intertidal zone of Southern Fujian

The oceans have rich resources of microbes. In recent years, searching for antibiotic of new structure and special biological activity from marine microbes has received great attention. This paper reports the results of the studies on the distribution and the antibacterial properties of Streptomyces in the intertidal zone of Southern Fujian. These studies have important significance to further investigate     

Constraining the origin and evolution of confined turbidite systems: southern Cretan margin,Eastern Mediterranean Sea (34°30–36°N)

Bathymetric, 9.5-kHz long-range sidescan sonar (OKEAN), seismic reflection and sediment-core data are used in the analysis of two tectonic troughs south of Crete, Eastern Mediterranean Sea. Here, up to 1.2 s two-way travel time (TWTT) of strata have accumulated since the Middle Miocene in association with extension in the South Aegean region. The study area comprises >100-km- long by >25-km-wide basins filled by sediments subdivided into two seismic units: (1) an upper Unit 1 deposited in sub-basins which follow the present-day configuration of the southern Cretan margin; (2) a basal Unit 2, more than 500 ms (TWTT) thick, accumulated in deeper half-graben/grabens distinct from the present-day depocentres. Both units overlap a locally stratified Unit 3 comprising the pre-Neogene core complex of Crete and Gavdos. In this work, the interpreted seismic units are correlated with the onshore stratigraphy, demonstrating that denudation processes occurring on Crete and Gavdos in response to major tectonic events have been responsible for high sedimentation rates along the proximal southern Cretan margin. Consequently, topographically confined sedimentary units have been deposited south of Crete in the last 12 Ma, including turbidites and other mass-flow deposits fed by evolving transverse and axial channel systems. Surface processes controlling facies distribution include the direct inflow of sediment from alluvial-fan systems and incising mountain rivers onto the Cretan slope, where significant sediment instability processes occur at present. In this setting, seismic profiles reveal eight different types of stratigraphic contacts on basin-margin highs, and basinal areas show evidence of halokinesis and/or fluid escape. The acquired data also show that significant changes to the margin’s configuration occurred in association with the post-Alpine tectonic and eustatic episodes affecting the Eastern Mediterranean.     

Geomorphological evidence for upslope canyon-forming processes on the northern KwaZulu-Natal shelf,SW Indian Ocean,South Africa

A geomorphological and statistical analysis of slope canyons from the northern KwaZulu-Natal continental margin is documented and compared with submarine canyons from the Atlantic margin of the USA. The northern KwaZulu-Natal margin is characterized by increasing upslope relief, concave slope-gradient profiles and features related to upslope growth of the canyon forms. Discounting slope-gradient profile, this morphology is strikingly similar to canyon systems of the New Jersey slope. Several phases of canyon incision indicate that downslope erosion is also an important factor in the evolution of the northern KwaZulu-Natal canyon systems. Despite the strong similarities between the northern KwaZulu-Natal and New Jersey slope-canyon systems, key differences are evident: (1) the concavity of the northern KwaZulu-Natal slope, contrasting with the ∼linear New Jersey slope; (2) the relative isolation of the northern KwaZulu-Natal canyons, rather than the dense clustering of the New Jersey canyons; and (3) the absence of strongly shelf-breaching canyons along the northern KwaZulu-Natal margin. In comparison with the New Jersey margin, we surmise a more youthful stage of canyon evolution, a result of either the canyons themselves being younger or the formative processes being less active. Less complicated patterns of erosion resulting from reduced sediment availability have developed in northern KwaZulu-Natal. The reduction in slope concavity on the New Jersey margin may be the result of grading of the upper slope by intensive headward erosion, a process more subdued—or less evident—on the KwaZulu-Natal margin.