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.     

Determination of fluid/melt partition coefficients by LA-ICPMS analysis of co-existing fluid and silicate melt inclusions: Controls on element partitioning

Analyses of co-existing silicate melt and fluid inclusions, entrapped in quartz crystals in volatile saturated magmatic systems, allowed direct quantitative determination of fluid/melt partition coefficients. Investigations of various granitic systems (peralkaline to peraluminous in composition, log fO₂ = NNO−1.7 to NNO+4.5) exsolving fluids with various chlorinities (1-14 mol/kg) allowed us to assess the effect of these variables on the fluid/melt partition coefficients (D). Partition coefficients for Pb, Zn, Ag and Fe show a nearly linear increase with the chlorinity of these fluid (D_Pb ∼ 6 ∗ m_Cl, D_Zn ∼ 8 ∗ m_Cl, D_Ag ∼ 4 ∗ m_Cl, D_Fe ∼ 1.4 ∗ m_Cl, where m_Cl is the molinity of Cl). This suggests that these metals are dissolved primarily as Cl-complexes and neither oxygen fugacity nor the composition of the melt affects significantly their fluid/melt partitioning. By contrast, partition coefficients for Mo, B, As, Sb and Bi are highest in low salinity (1-2 mol/kg Cl) fluids with maximum values of D_Mo ∼ 20, D_B ∼ 15, D_As ∼ 13, D_Sb ∼ 8, D_Bi ∼ 15 indicating dissolution as non-chloride (e.g., hydroxy) complexes. Fluid/melt partition coefficients of copper are highly variable, but highest between vapor like fluids and silicate melt (D_Cu ? 2700), indicating an important role for ligands other than Cl. Partition coefficients for W generally increase with increasing chlorinity, but are exceptionally low in some of the studied brines which may indicate an effect of other parameters. Fluid/melt partition coefficients of Sn show a high variability but likely increase with the chlorinity of the fluid (D_Sn = 0.3-42, D_W = 0.8-60), and decrease with decreasing oxygen fugacity or melt peraluminosity.     

Sedimentary petrology and geochemistry of siliciclastic rocks from the upper Jurassic Tordillo Formation (Neuquén Basin, western Argentina): Implications for provenance and tectonic setting

The Upper Jurassic Tordillo Formation is exposed along the western edge of the Neuquén Basin (west central Argentina) and consists of fluvial strata deposited under arid/semiarid conditions. The pebble composition of conglomerates, mineralogical composition of sandstones and pelitic rocks, and major- and trace-element geochemistry of sandstones, mudstones, and primary pyroclastic deposits are evaluated to determine the provenance and tectonic setting of the sedimentary basin. Conglomerates and sandstones derived almost exclusively from volcanic sources. The stratigraphic sections to the south show a clast population of conglomerates dominated by silicic volcanic fragments and a predominance of feldspathic litharenites. This framework composition records erosion of Triassic–Jurassic synrift volcaniclastic rocks and basement rocks from the Huincul arch, which was exhumed as a result of Late Jurassic inversion. In the northwestern part of the study area, conglomerates show a large proportion of mafic and acidic volcanic rock fragments, and sandstones are characterised by a high content of mafic volcanic rock fragments and plagioclase. These data suggest that the source of the sandstones and conglomerates was primarily the Andean magmatic arc, located west of the Neuquén Basin. The clay mineral assemblage is interpreted as the result of a complex set of factors, including source rock, climate, transport, and diagenesis. Postdepositional processes produced significant variations in the original compositions, especially the fine-grained deposits. The Tordillo sediments are characterised by moderate SiO₂ contents, variable abundances of K₂O and Na₂O, and a relatively high proportion of ferromagnesian elements. The degree of chemical weathering in the source area, expressed as the chemical index of alteration, is low to moderate. The major element geochemistry and Th/Sc, K/Rb, Co/Th, La/Sc, and Cr/Th values point to a significant input of detrital volcanic material of calcalkaline felsic and intermediate composition. However, major element geochemistry is not useful for interpreting the tectonic setting. Discrimination plots based on immobile trace elements, such as Ti, Zr, La, Sc, and Th, show that most data lie in the active continental margin field. Geochemical information is not sufficiently sensitive to differentiate the two different source areas recognized by petrographic and modal analyses of conglomerates and sandstones.     

The passivation of calcite by acid mine water. Column experiments with ferric sulfate and ferric chloride solutions at pH 2

Column experiments, simulating the behavior of passive treatment systems for acid mine drainage, have been performed. Acid solutions (HCl or H₂SO₄, pH 2), with initial concentrations of Fe(III) ranging from 250 to 1500 mg L⁻¹, were injected into column reactors packed with calcite grains at a constant flow rate. The composition of the solutions was monitored during the experiments. At the end of the experiments (passivation of the columns), the composition and structure of the solids were measured. The dissolution of calcite in the columns caused an increase in pH and the release of Ca into the solution, leading to the precipitation of gypsum and Fe–oxyhydroxysulfates (Fe(III)–SO₄–H⁺ solutions) or Fe–oxyhydroxychlorides (Fe(III)–Cl–H⁺ solutions). The columns worked as an efficient barrier for some time, increasing the pH of the circulating solutions from 2 to 6–7 and removing its metal content. However, after some time (several weeks, depending on the conditions), the columns became chemically inert. The results showed that passivation time increased with decreasing anion and metal content of the solutions. Gypsum was the phase responsible for the passivation of calcite in the experiments with Fe(III)–SO₄–H⁺ solutions. Schwertmannite and goethite appeared as the Fe(III) secondary phases in those experiments. Akaganeite was the phase responsible for the passivation of the system in the experiments with Fe(III)–Cl–H⁺ solutions.     

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.     

Tidal inundation mapping under enhanced land subsidence in Semarang,Central Java Indonesia

Tidal inundation by high tide under enhanced land subsidence is a damaging phenomenon and a major threat to the Semarang urban area in Indonesia. It impacts on economic activities, as well as the cost of an emergency program and causes interruption of pubic services, danger of infectious diseases and injury to human lives. This study examines a spatial analysis tool on the GIS-raster system for the tidal inundation mapping based on the subsidence-benchmark data and modified detail digital elevation model. Neighborhood operation and iteration model as a spatial analysis tool have been applied in order to calculate the encroachment of the tidal inundation on the coastal area. The resulting map shows that the tidal flood spreads to the lowland area and causes the inundation of coastal settlement, infrastructure, as well as productive agricultural land, i.e., the fish-pond area. The monitoring of the vulnerable area due to the tidal inundation under the scenario of extended land subsidence plays an important role in long-term coastal zone management in Semarang.