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.

An experimental study of the grain-scale processes of peridotite melting: implications for major and trace element distribution during equilibrium and disequilibrium melting

The grain-scale processes of peridotite melting were examined at 1,340°C and 1.5 GPa using reaction couples formed by juxtaposing pre-synthesized clinopyroxenite against pre-synthesized orthopyroxenite or harzburgite in graphite and platinum-lined molybdenum capsules. Reaction between the clinopyroxene and orthopyroxene-rich aggregates produces a melt-enriched, orthopyroxene-free, olivine + clinopyroxene reactive boundary layer. Major and trace element abundance in clinopyroxene vary systematically across the reactive boundary layer with compositional trends similar to the published clinopyroxene core-to-rim compositional variations in the bulk lherzolite partial melting studies conducted at similar P–T conditions. The growth of the reactive boundary layer takes place at the expense of the orthopyroxenite or harzburgite and is consistent with grain-scale processes that involve dissolution, precipitation, reprecipitation, and diffusive exchange between the interstitial melt and surrounding crystals. An important consequence of dissolution–reprecipitation during crystal-melt interaction is the dramatic decrease in diffusive reequilibration time between coexisting minerals and melt. This effect is especially important for high charged, slow diffusing cations during peridotite melting and melt-rock reaction. Apparent clinopyroxene-melt partition coefficients for REE, Sr, Y, Ti, and Zr, measured from reprecipitated clinopyroxene and coexisting melt in the reactive boundary layer, approach their equilibrium values reported in the literature. Disequilibrium melting models based on volume diffusion in solid limited mechanism are likely to significantly underestimate the rates at which major and trace elements in residual minerals reequilibrate with their surrounding melt. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

Geochemical and mineralogical evidence of tectonic and sedimentary factors controlling the origin of ferromanganese crusts associated to stratigraphic discontinuities (Betic Cordilleras, SE of Spain)

Ferromanganese crusts were found in carbonates of tectonostratigraphic units located in the northern and southern areas of the eastern External Subbetic of the Betic Cordilleras (SE Spain). The crusts are associated with four stratigraphic discontinuities of the Jurassic pelagic swells sequences: D1 (Late Carixian-Early Domerian), D2 (Middle Toarcian-Early Bajocian), D3 (Middle Bathonian-Middle Oxfordian), D4 (Early Tithonian-Late Albian). Two main textural types of crusts are distinguished. Type I crusts are thin and characterized by the presence of goethite, quartz, albite and phyllosilicates. Moreover, they show Si, Al, Mg, Na, Ti and K contents close to the European Shale Composite contents and Fe/Mn ratios (>350) higher than type II crusts. Type II crusts occur as thicker banded laminae and/or macrooncoids. They consist mainly of goethite and Mn-oxyhydroxides, which are enriched in REE, Co, Ni and Cu and show a strong Ce positive anomaly. After stratigraphical, mineralogical and geochemical data, the crust formation would be produced by the exposition of bottom sediments during long periods to a thin layer of oxidizing sea and porewater enriched in metallic elements. The textural and compositional variations between crusts can be explained by taking into account the bathymetric conditions. In shallower swells, the precipitation of a thick layer of banded type II crusts and in deeper areas, thin type I crusts were formed. Organic influence was only important in crusts from D3 of the northern area where textural evidence indicates the existence of seasonal periodically alternation between organism accretion and fine sedimentation. These were preceded and followed by phases in which the inorganic precipitation of oxides prevailed together with the fine sedimentation.     

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.     

Coupled hydrogeological and reactive transport modelling of the Simpevarp area (Sweden)

The Simpevarp area is one of the alternative sites being considered for the deep geological disposal of high level radioactive waste in Sweden. In this paper, a coupled regional groundwater flow and reactive solute transport model of the Simpevarp area is presented that integrates current hydrogeological and hydrochemical data of the area. The model simulates the current hydrochemical pattern of the groundwater system in the area. To that aim, a conceptual hydrochemical model was developed in order to represent the dominant chemical processes. Groundwater flow conditions were reproduced by taking into account fluid-density-dependent groundwater flow and regional hydrogeologic boundary conditions. Reactive solute transport calculations were performed on the basis of the velocity field so obtained. The model was calibrated and sensitivity analyses were carried out in order to investigate the effects of heterogeneities of hydraulic conductivity in the subsurface medium. Results provided by the reactive transport model are in good agreement with much of the measured hydrochemical data. This paper emphasizes the appropriateness of the use of reactive solute transport models when water-rock interaction reactions are involved, and demonstrates what powerful tools they are for the interpretation of hydrogeological and hydrochemical data from site geological repository characterization programs, by providing a qualitative framework for data analysis and testing of conceptual assumptions in a process-oriented approach.     

Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.     

Drought and institutional adaptation in the Great Plains of Alberta and Saskatchewan, 1914–1939

Agriculture in the southern Great Plains of Canada has been particularly vulnerable to prolonged episodes of drought. Using climate data and a precipitation minus potential evapotranspiration index, the extent of the region’s exposure to drought is examined. Between 1914 and 1917, the Dry Belt was particularly vulnerable to drought, whereas after 1928, a much larger region known as the Palliser Triangle covering most of southern Alberta and Saskatchewan was much more exposed to drought. These droughts provoked major institutional adaptation, in particular the establishment of the Special Areas Board by the Government of Alberta, and the creation of the Prairie Farm Rehabilitation Administration by the Government of Canada. Both organizations have proved to be relatively permanent public adaptations to the natural hazard of drought in the region. Moreover, these earlier experiences with prolonged drought as well as institution-building may be of value in helping the residents of the Palliser Triangle adapt to predicted climate changes in the future as well as anticipate some of the barriers to effective institutional adaptation.