Numerical Calculation of Equivalent Cell Permeability Tensors for General Quadrilateral Control Volumes

A new method for upscaling fine scale permeability fields to general quadrilateral-shaped coarse cells is presented. The procedure, referred to as the conforming scale up method, applies a triangle-based finite element technique, capable of accurately resolving both the coarse cell geometry and the subgrid heterogeneity, to the solution of the local fine scale problem. An appropriate averaging of this solution provides the equivalent permeability tensor for the coarse scale quadrilateral cell. The general level of accuracy of the technique is demonstrated through application to a number of flow problems. The real strength of the conforming scale up method is demonstrated when the method is applied in conjunction with a flow-based gridding technique. In this case, the approach is shown to provide results that are significantly more accurate than those obtained using standard techniques.     

Timing of deformation in the Norseman-Wiluna Belt, Yilgarn Craton, Western Australia

Establishing relative and absolute time frameworks for the sedimentary, magmatic, tectonic and gold mineralisation events in the Norseman-Wiluna Belt of the Archean Yilgarn Craton of Western Australia, has long been the main aim of research efforts. Recently published constraints on the timing of sedimentation and absolute granite ages have emphasized the shortcomings of the established rationale used for interpreting the timing of deformation events. In this paper the assumptions underlying this rationale are scrutinized, and it is shown that they are the source of significant misinterpretations. A revised time chart for the deformation events of the belt is established. The first shortening phase to affect the belt, D₁, was preceded by an extensional event D_1e and accompanied by a change from volcanic-dominated to plutonic-dominated magmatism at approximately 2685–2675 Ma. Later extension (D_2e) controlled deposition of the ca 2655 Ma Kurrawang Sequence and was followed by D₂, a major shortening event, which folded this sequence. D₂ must therefore have started after 2655 Ma—at least 20 Ma later than previously thought and after the voluminous 2670–2655 Ma high-Ca granite intrusion. Younger transcurrent deformation, D₃–D₄, waned at around 2630 Ma, suggesting that the crustal shortening deformation cycle D₂–D₄ lasted approximately 20–30 Ma, contemporaneous with low-volume 2650–2630 Ma low-Ca granites and alkaline intrusions. Time constraints on gold deposits suggest a late mineralisation event between 2640–2630 Ma. Thus, D₂–D₄ deformation cycle and late felsic magmatism define a 20–30 Ma long tectonothermal event, which culminated with gold mineralisation. The finding that D₂ folding took place after voluminous high-Ca granite intrusion led to research into the role of competent bodies during folding by means of numerical models. Results suggest that buoyancy-driven doming of pre-tectonic competent bodies trigger growth of antiforms, whereas non-buoyant, competent granite bodies trigger growth of synforms. The conspicuous presence of pre-folding granites in the cores of anticlines may be a result from active buoyancy doming during folding.     

Rivers, chemical weathering and Earth's climate

We detail the results of recent studies describing and quantifying the large-scale chemical weathering of the main types of continental silicate rocks: granites and basalts. These studies aim at establishing chemical weathering laws for these two lithologies, describing the dependence of chemical weathering on environmental parameters, such as climate and mechanical erosion. As shown within this contribution, such mathematical laws are of primary importance for numerical models calculating the evolution of the partial pressure of atmospheric CO₂ and the Earth climate at geological timescales. The major results can be summarized as follow: (1) weathering of continental basaltic lithologies accounts for about 30% of the total consumption of atmospheric CO₂ through weathering of continental silicate rocks. This is related to their high weatherability (about eight times greater than the granite weatherability); (2) a simple weathering law has been established for basaltic lithologies, giving the consumption of atmospheric CO₂ as a function of regional continental runoff, and mean annual regional temperature; (3) no such simple weathering law can be proposed for granitic lithologies, since the effect of temperature can only be identified for regions displaying high continental runoff; (4) a general law relating mechanical erosion and chemical weathering has been validated on small and large catchments. The consequences of these major advances on the climatic evolution of the Earth are discussed. Particularly, the impacts of the onset of the Deccan trapps and the Himalayan orogeny on the global carbon cycle are reinvestigated. To cite this article: B. Dupré et al., C. R. Geoscience 335 (2003).     

Re-Os Isotopic Measurements at the Femtomole Level in Natural Samples

Here we describe a new chemical separation method for Os and an improved mass spectrometric procedure for Re and Os. This technique is based on the selective extraction of OsO₄ from aqueous solution in liquid bromine. Among other advantages, this procedure avoids the cumbersome distillation procedure for Os and uses only commercial "off the shelf" PFA teflon labware. Blank levels for 0.5 g sample sizes are: 0.06 pg and 0.5 pg for Os and Re respectively. Samples containing as little as 1 pg g⁻¹ Os can be analysed reliably with this method. Adaptation to other dissolution methods is also discussed.     

Transport in a 2-D saturated porous medium: A new method for particle tracking

A new method for solving the transport equation based on the management of a large numbe of particles in a discretized 2-D domain is presented. The method uses numerical variables to represent the number of particles in a given mesh and is more complex than the 1-D problem. The first part of the paper focuses on the specific management of particles in a 2-D problem. The method also would be valid for three dimensions as long as the medium can be modeled similar to a layered system. As the particles are no longer tracked individually, the algorithm is fast and does not depend on the number of particles present. The numerical tests show that the method is nearly numerical dispersion free and permits accurate calculations even for simulations of low-concentration transport. Because each mesh is considered as a closed system between two successive time steps, it is easy to add adsorption phenomenon without any problem of numerical stability. The model is tested under conditions that are extremely demanding for its operating mode and gives a good fit to analytical solutions. The conditions in which it can be used to best advantage are discussed.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Influence of water level fluctuation on the spatial dynamic of the water chemistry at Lake Fertõ/Neusiedler See

Lake Fertõ/Neusiedler See is a steppe lake situated on the Hungarian-Austrian border with a surface area of 309 km². It is an extremely shallow lake with regulated outflow. 54% of the whole lake and 85% of the Hungarian part (75 km²) is covered by reed. There are numerous reedless areas (inner ponds) of variable size within the reed belt, which is enmeshed with canals connecting the inner ponds with the open water areas. Fundamental differences were recorded in the water chemistry of the three water type areas, open water, canals and inner ponds. The degree of these differences is basically determined by the distance from the open water areas and the water level of the actual years. An extremely high, an average and an extremely low water level year was examined, in order to present the spatial and temporal differences in the water chemistry of the investigated water bodies. Differences due to the water level fluctuation between years in the three different water types are discussed on the basis of electrical conductivity, cation, and anion concentration changes of the water. Our results support the concept of high instability of Lake Fertõ and describe habitats with remarkable temporal and spatial variability. The high influence of water level on habitat conditions at different sites of the lake recommend a carefull water level regulation strategy.     

Organic matter in a low productivity anoxic intraplatform basin in the Triassic Tethys

This paper presents isotopic, bulk geochemical and biomarker data measured on organic matter accumulated in a narrow extensional basin developed at the oceanward margin of the huge Triassic carbonate platform in the Alps–Appenines domain. The integration of isotope signatures, organic petrographical and biomarker evidence together with the composition of kerogen pyrolysates suggests immature organic matter predominantly of algal origin with a minor, but not negligible, higher plant derived and moderate bacterial contribution for the entire sequence. The mineral sources are dominated by platform-derived subtidal Dachstein limestone with a minor palaeosol input and a moderate contribution of autochtonous quartz. Nevertheless, parallel variations observed in the mineral content, as well as in the amount and the quality of the organic matter reflect variations in the palaeoenvironment. The increased humidity, existing in the period of the accumulation of the upper section of the sequence, led to the restriction of dolomitization. A slightly greater higher plant derived contribution, in this section, is evidenced by the composition of bitumen and the results on GC and GC/MS analyses on the non-aromatic hydrocarbon fraction of bitumen. Moreover, the climate-induced weathering enchanced the primary productivity and resulted in a pronounced increase in the TOC content. The average estimated value of the planktonic productivity is about three times higher for the calcite-rich sequence than the dolomite-rich one, being 44.2 and 15.3 tC_org/m²/Ma, respectively. The low to moderate planktonic productivity shows that anoxic conditions, observed for the entire succession, are a consequence of the stagnant water stratification rather than high planktonic productivity. Depth trends in the data measured on kerogens (HI, OI, δ¹³Corg values, composition of pyrolysate) together with the δ¹⁸O excursions and Δδ¹³C values appear to be controlled by sea-level fluctuations. Consistent with the high abundance of alkyltiophenes in kerogen pyrolysates, the high S_org/C ratios (ranging between 0.05 and 0.10) suggest the importance of natural sulfurization in the formation of the sulfur-rich type-II-S kerogen occurring in all of the samples.     

Transfer of Taiwanese ideas and technology to The Gambia, West Africa: a viable approach to rural development?

The development of Taiwan's agricultural sector contributed significantly to the country's transformation from virtual poverty after World War II to a modern industrial nation. Success resulted from the use of appropriate technologies, developed through an exceptionally close working relationship between government, research and development institutions, extension services and farmers in Taiwan. As part of their overseas development aid programme, the Republic of China on Taiwan has established Technical Missions in those developing countries where it receives political recognition. Among these is The Gambia, West Africa, where one focus of Taiwan's current Technical Mission is the development of women's market gardens. Although the gardens have proved successful, there is doubt as to whether their achievements will be sustained once Taiwan's Technical Mission in The Gambia has ended. Field research in Taiwan and in The Gambia, based largely on participatory techniques, has revealed that the low-tech methods introduced from Taiwan are appropriate for the Gambian environment. However, much depends on whether women participating in the schemes can commit sufficient time and energy to the gardens, as at present this is a limiting factor. The potential of the gardens may also be constrained by a lack of appropriate institutional infrastructure, particularly the scope for profitable disposal of garden produce. A major conclusion of the article is that successful transfer of appropriate technology to The Gambia is likely to be impeded by poorly developed coordination between groups that proved to be key 'players' in the development of agriculture in Taiwan: government, research and development institutions, extension services and producers.