Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Metamorphic phase relations in orthopyroxene-bearing granitoids: implication for high-pressure metamorphism and prograde melting in the continental crust

In this work, the factors controlling the formation and preservation of high-pressure mineral assemblages in the metamorphosed orthopyroxene-bearing metagranitoids of the Sandmata Complex, Aravalli-Delhi Mobile Belt (ADMB), northwestern India have been modelled. The rocks range in composition from farsundite through quartz mangerite to opdalite, and with varying K₂O, Ca/(Ca + Na)_rock and FeO_tot + MgO contents. A two stage metamorphic evolution has been recorded in these rocks.
An early hydration event stabilized biotite with or without epidote at the expense of magmatic orthopyroxene and plagioclase. Subsequent high-pressure granulite facies metamorphism (∼15 kbar, ∼800 °C) of these hydrated rocks produced two rock types with contrasting mineralogy and textures. In the non-migmatitic metagranitoids, spectacular garnet ± K-feldspar ± quartz corona was formed around reacting biotite, plagioclase, quartz and/or pyroxene. In contrast, biotite ± epidote melting produced migmatites, containing porphyroblastic garnet incongruent solids and leucosomes.
Applying NCKFMASHTO T–M (H₂O) and P–T pseudosection modelling techniques, it is demonstrated that the differential response of these magmatic rocks to high-pressure metamorphism is primarily controlled by the scale of initial hydration. Rocks, which were pervasively hydrated, produced garnetiferous migmatites, while for limited hydration, the same metamorphism formed sub-solidus garnet-bearing coronae. Based on the sequence of mineral assemblage evolution and the mineral compositional zoning features in the two metagranitoids, a clockwise metamorphic P–T path is constrained for the high-pressure metamorphic event. The finding has major implications in formulating geodynamic model of crustal amalgamation in the ADMB.     

Reliability assessment of granular filters in embankment dams

Empirical criteria have been used successfully to design filters of most embankment large dam projects throughout the world. However, these empirical rules are only applicable to a particular range of soils tested in laboratory and do not take into account the variability of the base material and filter particle sizes. In addition, it is widely accepted that the safety of fill dams is mainly dependent on the reliability of their filter performance. The work herein presented consists in a new general method for assessing the probability of fulfilling any empirical filter design criteria accounting for base and filter heterogeneity by means of first‐order reliability methods (FORM), so that reliability indexes and probabilities of fulfilling any particular criteria are obtained. This method will allow engineers to estimate the safety of existing filters in terms of probability of fulfilling their design criteria and might also be used as a decision tool on sampling needs and material size tolerances during construction. In addition, sensitivity analysis makes possible to analyse how reliabilities are influenced by different sources of input data. Finally, in case of a portfolio risk assessment, this method will allow engineers to compare the safety of several existing dams in order to prioritize safety investments and it is expected to be a very useful tool to evaluate probabilities of failure due to internal erosion. 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.     

Late Quaternary vegetation changes around Lake Rutundu,Mount Kenya,East Africa: evidence from grass cuticles,pollen and stable carbon isotopes

Woody, subalpine shrubs and grasses currently surround Lake Rutundu, Mount Kenya. Multiple proxies, including carbon isotopes, pollen and grass cuticles, from a 755‐cm‐long core were used to reconstruct the vegetation over the past 38 300 calendar years. Stable carbon‐isotope ratios of total organic carbon and terrestrial biomarkers from the lake sediments imply that the proportion of terrestrial plants using the C₄ photosynthetic pathway was greater during the Late Pleistocene than in the Holocene. Pollen data show that grasses were a major constituent of the vegetation throughout the Late Pleistocene and Holocene. The proportion of grass pollen relative to the pollen from other plants was greatest at the last glacial maximum (LGM). Grass cuticles confirm evidence that C₄ grass taxa were present at the LGM and that the majority followed the cold‐tolerant NADP‐MEC₄ subpathway. Copyright © 2003 John Wiley & Sons, Ltd.     

Palaeotemperature reconstructions of the European permafrost zone during marine oxygen isotope Stage 3 compared with climate model results

A palaeotemperature reconstruction based on periglacial phenomena in Europe north of approximately 51 °N, is compared with high‐resolution regional climate model simulations of the marine oxygen isotope Stage 3 (Stage 3) palaeoclimate. The experiments represent Stage 3 warm (interstadial), Stage 3 cold (stadial) and Last Glacial Maximum climatic conditions. The palaeotemperature reconstruction deviates considerably for the Stage 3 cold climate experiments, with mismatches up to 11 °C for the mean annual air temperature and up to 15 °C for the winter temperature. However, in this reconstruction various factors linking climate and permafrost have not been taken into account. In particular a relatively thin snow cover and high climatic variability of the glacial climate could have influenced temperature limits for ice‐wedge growth. Based on modelling the 0 °C mean annual ground temperature proves to be an appropriate upper temperature limit. Using this limit, mismatches with the Stage 3 cold climate experiments have been reduced but still remain. We therefore assume that the Stage 3 ice wedges were generated during short (decadal time‐scale) intervals of extreme cold climate, below the mean temperatures indicated by the Stage 3 cold climate model simulations. Copyright © 2003 John Wiley & Sons, Ltd.     

The importance of oxygen isotope provenance in relation to solute content of bulk meltwaters at Imersuaq Glacier,West Greenland

Stable oxygen isotope analysis and measurement of several dissolved cations and anions of bulk meltwater samples have provided information about the hydrochemical environment of the glacial hydrological system at Imersuaq Glacier, an outlet tongue from the Greenland ice‐sheet, West Greenland. The samples were collected at frequent intervals during the period 20–28 July 2000 in a small (<20 L s^?1) englacial meltwater outlet at the glacier margin. The results document the following findings: (i) a marked diurnal variation of δ¹⁸O is related to the composition of oxygen isotope provenances, mainly near‐marginal local superimposed ice and basal up‐sheared ice further up‐glacier; (ii) a relationship is seen between all base cations (Na⁺, K⁺, Ca²⁺, Mg²⁺), SO₄^2? and δ¹⁸O, indicating that solute acquisition is provided by solid–solution contact with the up‐sheared ice—as the relationship with Cl^? is weak the influence of seasalt‐derived solutes is small in the area; (iii) when the melt rate is high, two diurnal maxima of δ¹⁸O values and solute concentrations are measured, and it is suggested that a snow meltwater component is responsible for the second maximum of δ¹⁸O—a short residence time leads to a delayed decrease in ion concentrations. Copyright © 2003 John Wiley & Sons, Ltd.     

Morphological analysis of the drainage system in the Eastern Alps

We study the morphology of the major rivers draining the Eastern Alps to test whether the active tectonics of this part of the orogen is reflected in the shape of channel profiles of the river network. In our approach we compare channel profiles measured from digital elevation models with numerically modelled channel profiles using a stream power approach. It is shown that regions of high stream power coincide largely with regions of highest topography and largest uplift rates, while the forelands and the Pannonian Basin are characterised by a significantly lower stream power. From stream power modelling we conclude that there is young uplift at the very east of the Eastern Alps, in the Bohemian Massif and in the Pohorje Range. The impact of the Pleistocene glaciations is explored by comparing properties of rivers that drain in proximal and distal positions relative to the ice sheet during the last glacial maximum. Our analysis shows that most knick points, wind gaps and other non-equilibrium features of catchments covered by ice during the last glaciations (Salzach, Enns) can be correlated with glacial processes. In contrast the ice free catchments of the Mur and Drava are characterized by channels in morphological equilibrium at the first approximation and are showing only weak evidence of the strong tectonic activity within these catchments. Finally, the channel profiles of the Adige and the divide between the upper Rhine and Danube catchments differ significantly from the other catchments. We relate this to the fact that the Adige and the Rhine respond to different base levels from the remainder of the Eastern Alps: The Adige may preserve a record from the Messininan base level change and the Rhine is subject to the base level lowering in the Rhine Graben.