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.

Basaltic explosive volcanism: Constraints from deposits and models

Basaltic pyroclastic volcanism takes place over a range of scales and styles, from weak discrete Strombolian explosions (～10²–10³ kg s^?1) to Plinian eruptions of moderate intensity (10⁷–10⁸ kg s^?1). Recent well-documented historical eruptions from Etna, Kīlauea and Stromboli typify this diversity. Etna is Europe's largest and most voluminously productive volcano with an extraordinary level and diversity of Strombolian to subplinian activity since 1990. Kīlauea, the reference volcano for Hawaiian fountaining, has four recent eruptions with high fountaining (>400 m) activity in 1959, 1960, 1969 (–1974) and 1983–1986 (–2008); other summit (1971, 1974, 1982) and flank eruptions have been characterized by low fountaining activity. Stromboli is the type location for mildly explosive Strombolian eruptions, and from 1999 to 2008 these persisted at a rate of ca. 9 per hour, briefly interrupted in 2003 and 2007 by vigorous paroxysmal eruptions. Several properties of basaltic pyroclastic deposits described here, such as bed geometry, grain size, clast morphology and vesicularity, and crystal content are keys to understand the dynamics of the parent eruptions.The lack of clear correlations between eruption rate and style, as well as observed rapid fluctuations in eruptive behavior, point to the likelihood of eruption style being moderated by differences in the fluid dynamics of magma and gas ascent and the mechanism by which the erupting magma fragments. In all cases, the erupting magma consists of a mixture of melt and gaseous bubbles. The depth and rate of degassing, melt rheology, bubble rise and coalescence rates, and extent of syn-eruptive microlite growth define complex feedbacks that permit reversible shifts between fragmentation mechanisms and in eruption style and intensity. However, many basaltic explosive eruptions end after an irreversible shift to open-system outgassing and microlite crystallization in melt within the conduit.Clearer understanding of the factors promoting this diversity of basaltic pyroclastic eruptions is of fundamental importance in order to improve understanding of the range of behaviors of these volcanoes and assess hazards of future explosive events at basaltic volcanoes. The three volcanoes used for this review are the sites of large and growing volcano-tourism operations and there is a public need both for better knowledge of the volcanoes’ behavior and improved forecasting of the likely course of future eruptions.     

Low-pressure,water-assisted anatexis of basic dykes in a contact metamorphic aureole,Fuerteventura (Canary Islands): oxygen isotope evidence for a meteoric fluid origin

Migmatites produced by low-pressure anatexis of basic dykes are found in a contact metamorphic aureole around a pyroxenite–gabbro intrusion (PX2), on Fuerteventura. Dykes outside and inside the aureole record interaction with meteoric water, with low or negative δ¹⁸O whole-rock values (+0.2 to −3.4‰), decreasing towards the contact. Recrystallised plagioclase, diopside, biotite and oxides, from within the aureole, show a similar evolution with lowest δ¹⁸O values (−2.8, −4.2, −4.4 and −7.6‰, respectively) in the migmatite zone, close to the intrusion. Relict clinopyroxene phenocrysts preserved in all dykes, retain typically magmatic δ¹⁸O values up to the anatectic zone, where the values are lower and more heterogeneous. Low δ¹⁸O values, decreasing towards the intrusion, can be ascribed to the advection of meteoric water during magma emplacement, with increasing fluid/rock ratios (higher dyke intensities towards the intrusion acting as fluid-pathways) and higher temperatures promoting increasing exchange during recrystallisation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Temperature effects on non-bridging oxygen and aluminum coordination number in calcium aluminosilicate glasses and melts

Configurational changes with temperature are important for the thermodynamic and transport properties of most aluminosilicate melts, but in general are not well understood. Here, we present high-resolution ²⁷Al and ¹⁷O NMR data on several calcium aluminosilicate glasses prepared with varying quench rates and thus with fictive temperatures that span ranges up to about 200 K. In all compositions the content of five-coordinated aluminum increases with fictive temperature, in agreement with recent high temperature NMR data on melts. In a glass of CaAl₂Si₂O₈ (“anorthite”) composition, the content of non-bridging oxygens also increases with temperature; however this effect was not observed in a sample with a much higher CaO/Al₂O₃ ratio. We present a consistent notation for reactions among structural species in these systems that clarify why in some cases, high-coordinated network cations may appear on the same side of the reaction, while in others they occur on the opposite sides: the key difference is in accounting for all coordination changes for oxygens. Mixing of non-bridging oxygens and of high-coordinated aluminum make significant contributions to the overall configurational entropy and heat capacity of the melts, as does the mixing of various bridging oxygens and of tetrahedral network cations. Other, less well known, types of increase in disorder with temperature may be important as well.     

Silicon isotope fractionation in bamboo and its significance to the biogeochemical cycle of silicon

A systematic investigation on silica contents and silicon isotope compositions of bamboos was undertaken. Seven bamboo plants and related soils were collected from seven locations in China. The roots, stem, branch and leaves for each plant were sampled and their silica contents and silicon isotope compositions were determined. The silica contents and silicon isotope compositions of bulk and water-soluble fraction of soils were also measured. The silica contents of studied bamboo organs vary from 0.30% to 9.95%. Within bamboo plant the silica contents show an increasing trend from stem, through branch, to leaves. In bamboo roots the silica is exclusively in the endodermis cells, but in stem, branch and leaves, the silica is accumulated mainly in epidermal cells. The silicon isotope compositions of bamboos exhibit significant variation, from −2.3‰ to 1.8‰, and large and systematic silicon isotope fractionation was observed within each bamboo. The δ³⁰Si values decrease from roots to stem, but then increase from stem, through branch, to leaves. The ranges of δ³⁰Si values within each bamboo vary from 1.0‰ to 3.3‰. Considering the total range of silicon isotope composition in terrestrial samples is only 7‰, the observed silicon isotope variation in single bamboo is significant and remarkable. This kind of silicon isotope variation might be caused by isotope fractionation in a Rayleigh process when SiO₂ precipitated in stem, branches and leaves gradually from plant fluid. In this process the Si isotope fractionation factor between dissolved Si and precipitated Si in bamboo (α_pre-sol) is estimated to be 0.9981. However, other factors should be considered to explain the decrease of δ³⁰Si value from roots to stem, including larger ratio of dissolved H₄SiO₄ to precipitated SiO₂ in roots than in stem. There is a positive correlation between the δ³⁰Si values of water-soluble fractions in soils and those of bulk bamboos, indicating that the dissolved silicon in pore water and phytoliths in soil is the direct sources of silicon taken up by bamboo roots. A biochemical silicon isotope fractionation exists in process of silicon uptake by bamboo roots. Its silicon isotope fractionation factor (α_bam-wa) is estimated to be 0.9988. Considering the distribution patterns of SiO₂ contents and δ³⁰Si values among different bamboo organs, evapotranspiration may be the driving force for an upward flow of a silicon-bearing fluid and silica precipitation. Passive silicon uptake and transportation may be important for bamboo, although the role of active uptake of silicic acid by roots may not be neglected. The samples with relatively high δ³⁰Si values all grew in soils showing high content of organic materials. In contrast, the samples with relatively low δ³⁰Si values all grew in soil showing low content of organic materials. The silicon isotope composition of bamboo may reflect the local soil type and growth conditions. Our study suggests that bamboos may play an important role in global silicon cycle.     

Framing in geographical analysis of environmental conflicts: Theory, methodology and three case studies

When applied to environmental concerns, framing offers a rigorous conceptual and analytic approach with potentially practical significance for dealing with complex issues relevant to geography. The purpose of this article is to introduce framing concepts, typology and modes of analysis to address issues which geographers commonly treat: disputes over land ownership and uses, competition for water resources, cultural clashes over control of territory and resources and the impact of spatial patterns and structures when siting noxious sites. Framing analysis is applied to three cases of environmental conflict in Israel and, more generally, for understanding conflicts revolving around the management of human habitats in relation to the physical environment.     

Models of natural and human dynamics in forest landscapes: Cross-site and cross-cultural synthesis

We synthesize the study of coupled natural and human systems across sites and cultures through a process of simplification and abstraction based on multiple dimensions of human-nature connectedness: satisfaction of basic needs, psycho-cultural connectedness and regulation of use of natural resources. We thus provide both a place-based and general understanding of value-driven anthropogenic environmental change and response. Two questions guide this research: what are the crucial stakeholder values that drive land use decisions and thus land cover change? And how can knowledge of these values be used to make decisions and policies that sustain both the human and natural systems in a place? To explore these questions we build simulation models of four study sites, two in the State of Texas, United States, and two in Venezuela. All include protected areas, though they differ in the specifics of vegetation and land use. In the Texas sites, relatively affluent individuals are legally converting forests to residential, commercial, and industrial uses, while in Venezuela landless settlers are extra-legally converting forests for purposes of subsistence agriculture. Contemporary modeling techniques now facilitate simulations of stakeholder and ecosystem dynamics revealing emergent patterns. Such coupled human and natural systems are currently recognized as a form of biocomplexity. Our modeling framework is flexible enough to allow adaptation to each of the study sites, capturing the essential features of the respective natural and anthropogenic land use changes and stakeholder reactions. The interactions between human stakeholders are simulated using multi-agent models that act on forest landscape models, and receive feedback of the effects of these actions on ecological habitats and hydrological response. The multi-agent models employ a formal logic-based method for the Venezuelan sites and a decision analysis approach using multi-attribute utility functions for the Texas sites, differing more in style and emphasis than in substance. Our natural-systems models are generic and can be tailored according to site-specific conditions. Similar models of tree growth and patch transitions are used for all the study sites and the differing responses to environmental variables are specified for each local species and terrain conditions.     

Pressures of Crystallization of Icelandic Magmas

Iceland lies astride the Mid-Atlantic Ridge and was createdby seafloor spreading that began about 55 Ma. The crust is anomalouslythick (20–40 km), indicating higher melt productivityin the underlying mantle compared with normal ridge segmentsas a result of the presence of a mantle plume or upwelling centeredbeneath the northwestern edge of the Vatnajökull ice sheet.Seismic and volcanic activity is concentrated in 50 km wideneovolcanic or rift zones, which mark the subaerial Mid-AtlanticRidge, and in three flank zones. Geodetic and geophysical studiesprovide evidence for magma chambers located over a range ofdepths (1·5–21 km) in the crust, with shallow magmachambers beneath some volcanic centers (Katla, Grimsvötn,Eyjafjallajökull), and both shallow and deep chambers beneathothers (e.g. Krafla and Askja). We have compiled analyses ofbasalt glass with geochemical characteristics indicating crystallizationof ol–plag–cpx from 28 volcanic centers in the Western,Northern and Eastern rift zones as well as from the SouthernFlank Zone. Pressures of crystallization were calculated forthese glasses, and confirm that Icelandic magmas crystallizeover a wide range of pressures (0·001 to 1 GPa), equivalentto depths of 0–35 km. This range partly reflects crystallizationof melts en route to the surface, probably in dikes and conduits,after they leave intracrustal chambers. We find no evidencefor a shallow chamber beneath Katla, which probably indicatesthat the shallow chamber identified in other studies containssilica-rich magma rather than basalt. There is reasonably goodcorrelation between the depths of deep chambers (> 17 km)and geophysical estimates of Moho depth, indicating that magmaponds at the crust–mantle boundary. Shallow chambers (<7·1 km) are located in the upper crust, and probablyform at a level of neutral buoyancy. There are also discretechambers at intermediate depths (11 km beneath the rift zones),and there is strong evidence for cooling and crystallizing magmabodies or pockets throughout the middle and lower crust thatmight resemble a crystal mush. The results suggest that themiddle and lower crust is relatively hot and porous. It is suggestedthat crustal accretion occurs over a range of depths similarto those in recent models for accretionary processes at mid-oceanridges. The presence of multiple stacked chambers and hot, porouscrust suggests that magma evolution is complex and involvespolybaric crystallization, magma mixing, and assimilation. KEY WORDS: Iceland rift zones; cotectic crystallization; pressure; depth; magma chamber; volcanic glass     

Characterization and cycling of atmospheric mercury along the central US Gulf Coast

Concentrations of atmospheric Hg species, elemental Hg (Hg^°), reactive gaseous Hg (RGM), and fine particulate Hg (Hg-PM_2.5) were measured at a coastal site near Weeks Bay, Alabama from April to August, 2005 and January to May, 2006. Mean concentrations of the species were 1.6 ± 0.3 ng m⁻³, 4.0 ± 7.5 pg m⁻³ and 2.7 ± 3.4 pg m⁻³, respectively. A strong diel pattern was observed for RGM (midday maximum concentrations were up to 92.7 pg m⁻³), but not for Hg^° or Hg-PM_2.5. Elevated RGM concentrations (>25 pg m⁻³) in April and May of 2005 correlated with elevated average daytime O₃ concentrations (>55 ppbv) and high light intensity (>500 W m⁻²). These conditions generally corresponded with mixed continental-Gulf and exclusively continental air mass trajectories. Generally lower, but still elevated, RGM peaks observed in August, 2005 and January–March, 2006 correlated significantly (p < 0.05) with peaks in SO₂ concentration and corresponded to periods of high light intensity and lower average daytime O₃ concentrations. During these times air masses were dominated by trajectories that originated over the continent. Elevated RGM concentrations likely resulted from photochemical oxidation of Hg^° by atmospheric oxidants. This process may have been enhanced in and by the near-shore environment relative to inland sites. The marine boundary layer itself was not found to be a significant source of RGM.