THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

A transdisciplinary analysis of water problems in the mountainous karst areas of Morocco

Water management and engineering in the karstic High Atlas of Morocco are difficult tasks under the prevailing geological, hydrogeological, geomorphological, vegetational and climatic conditions. It is important to be able to understand and predict the characteristics and availability of water for future water planning in the region under changing climatic and agricultural conditions. An interdisciplinary analysis of problems and adequate hydrological modelling tools developed by geologists, hydrologists and biologists are necessary. The karst areas of the High Atlas Mountains are characterised by impermeable triassic basalt underlying substantial subsurface reservoirs with high potential discharge rates. The karst groundwater aquifers are extensive but largely unknown in dimension, probably with a hierarchical network of groundwater flow paths. It is estimated that approximately 70% of the surface water is directly lost to groundwater. Steep landslide- and debris flow prone slopes exist next to coarse-grained, highly porous river beds. Infrequent, high intensity rainfall or snowmelt causes a particularly high flood risk to these karst areas. In addition, agriculture and land use changes have degraded the karst areas. The most important driving forces for degradation include permanent overgrazing even during droughts and the use of firewood by a continually growing population. Large scale degradation of vegetation has occurred in the oro-mediterranean (mountainous Mediterranean) zone, between 2600 and 3400 m which coincides with the most important zone for karstic groundwater creation. The combination of high amounts of groundwater flow and rapid surface flow due to sparse vegetation has increased the problems of flood flow.     

Water and sanitation in the Buenos Aires metropolitan region: Fragmented markets, splintering effects?

This paper discusses the provision of water and sanitation services and the related urban impacts in Buenos Aires circa 2005. The first part of the paper focuses on the metropolitan region (BAMR) which is larger than the area served by the Aguas Argentinas S.A. (AASA) concession (i.e. the capital and its conurbano). It highlights a form of institutional fragmentation that is reflected in the diverse management practices of the numerous service providers and results in differentiated levels of and access to services. The paper discusses the technical, economic and regulatory factors that have historically resulted in this institutional fragmentation. The second part of the article focuses on AASA’s water concession agreement. It demonstrates that, although the agreement had an in-built splintering logic, its implementation led to a more “integrated approach” based on the “reality on the ground”. On the whole, the conclusion regarding the application of the splintering urbanism concept to Buenos Aires is a mixed one and depends on the scale of observation. Although at regional level (BAMR), a form of “institutional fragmentation” would appear to be obvious, it is the result of a long and complex historical process and, while recent neo-liberal reforms in service provisions have “cemented” this fragmentation, they did not actually cause it. If we take the AASA service area in isolation, the trend in the decade following the contracting out of the service has been one of improved service provision and access—although this was not achieved through increased integration of the supply system. Finally, in the case of the Buenos Aires metropolitan region, the splintering urbanism thesis defended by Marvin and Graham needs to be nuanced in several significant ways.     

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.     

Geostatistical estimation of chemical contamination in stream sediments: The case study of Vale das Gatas mine (northern Portugal)

Based on an environmental geochemistry case study carried out in the neighbourhood of a W–Sn abandoned mine, the pollution in stream sediments was modelled through a Global Contamination Index. Such an index permits one to summarize the combination of deleterious elements in a single variable, obtained by the projection of samples onto the first axis of a PCASD (Principal Components Analysis of Standardized Data) applied to the entire n × p matrix containing the available concentrations of p = 16 elements in the set of n = 220 collected samples.In order to provide a sound basis for a coherent planning of the remediation process which will be put in operation in the affected area, it is necessary to balance the costs of reclaiming with the probabilities of exceeding the upper limits accepted for concentrations of environmentally harmful elements in sediments. Given these limits, they are back-transformed in the index values, providing a practical threshold between ‘clean’ and ‘contaminated’ samples. On the other hand, the minimum dimension of the cell to be reclaimed is restrained by the selected remediation process to be applied in the affected area. Hence, to meet the constraints of such a remediation process, it is required to estimate the probabilities of exceeding the index threshold in technologically meaningful sub-areas. For this end, the Indicator Block Kriging technique was applied, producing a series of maps where sub-areas to be reclaimed can be spotted for different probability levels. These maps, on which the decision making remediation agency can rely for its cost-benefit analysis, take into account both the spatial structure of ‘clean’ vs. ‘contaminated’ samples and the constraints of the reclaiming process.     

Flexural isostatic response of the Alps to increased Quaternary erosion recorded by foreland basin remnants,SE France

We test the hypothesis that flexural isostatic compensation of the mass removed by enhanced Quaternary erosion is responsible for uplift of the Western European Alps and their forelands. We use two well‐preserved and well‐dated (1.8 Ma) abandonment surfaces of foreland basin remnants in SE France (the Chambaran and Valensole plateaux) as passive benchmarks for tilting of the foreland. Estimating their initial slope from morphometric scaling relationships, we determine bulk post‐depositional tilting of 0.5–0.8% for these surfaces. The calculated isostatic response of the Alpine lithosphere to erosional unloading, using the method recently proposed by Champagnac et al. [Geology 35 (2007) 195–198] , yields a predicted tilting of 0.3–0.4% in the considered areas, explaining approximately half of the determined post‐depositional tilting. Such long‐term deformation being insensitive to cyclic loading/unloading because of glaciations, we suspect the other half to be related to as yet undetermined long‐wavelength and long‐lived tectonic process(es).     

Sur la présence de nickel dans un gisement de bauxite près de Mégare

Résumé L'exploitation récente d'un gisement de bauxite près de Mégare (Grèce) a fait apparaître sur les calcaires de base une concentration verte d'oxyde de nickel relativement élevée. La présence dans le gisement d'argiles rouges, renfermant de fortes proportions d'oxyde de fer avec des teneurs relativement faibles en alumine et des traces d'oxyde de nickel, induit à penser qu'il s'agit d'un produit latéritique dérivé d'une roche ultrabasique et transporté ultérieurement sur la région calcaire karstifiée avec ses remplissages bauxitiques. L'étude se termine par une hypothèse sur la genèse de ce gîte singulier.

---

The recent mining of a bauxite deposit near Megara (Greece) has shown on the footwall limestone a relatively high concentration of a green nickel oxide. The occurrence in the same deposit of red clays with a rather high percentage of iron oxide and a low content of alumina, together with some traces of nickel, induces to think that this material is composed of a lateritic soil derived from an ultrabasic rock and subsequently transported on karstified limestone with its bauxitic fillings. — The study ends with an hypothesis on the genesis of this unusual deposit.

     

Longitudinal strength of a high-speed ferry

The longitudinal strength of the high-speed ferry was investigated by subjecting the ship's hull girder to long-term loads obtained from a frequency-domain panel code. Prior to the statistical analysis, linearly computed transfer functions were corrected for nonlinear effects, yielding two sets of transfer functions valid for different wave amplitudes. One set corresponded to the hogging condition; the other set, to the sagging condition. Two regular equivalent design waves were specified that resulted in loads representing the most severe global design load conditions. The still-water loading condition, yielding a still-water vertical bending moment in hogging, was superimposed on the wave-induced loads to obtain the total (design) loads in hogging. For the sagging condition only, additional impact-related loads were superimposed to obtain the total (design) loads in sagging. A finite element model of the ship's structure was subjected to pressure distributions according to the two regular design waves. For comparison with classification society rule values, a simple beam theory strength analysis of the ship's midship section was performed first, and then another finite element analysis was carried out, whereby the imposed loads were tuned to the rule values of vertical bending moments. Rule-based magnitudes of nominal maximum longitudinal stress deviated significantly (25–39%) from comparable stresses obtained by the panel code based finite element analysis. However, stresses obtained from the rule-based finite element analysis agreed more favorably, especially in hogging. In the uppermost deck, for example, the panel code based compressive stress was only 9% larger than the comparable stress from the rule-based finite element analysis.