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.     

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.     

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.     

Association of dissolved aluminum with silica: Connecting molecular structure to surface reactivity using NMR

We studied uptake mechanisms for dissolved Al on amorphous silica by combining bulk-solution chemistry experiments with solid-state Nuclear Magnetic Resonance techniques (²⁷Al magic-angle spinning (MAS) NMR, ²⁷Al{¹H} cross-polarization (CP) MAS NMR and ²⁹Si{¹H} CP-MAS NMR). We find that reaction of Al (1 mM) with amorphous silica consists of at least three reaction pathways; (1) adsorption of Al to surface silanol sites, (2) surface-enhanced precipitation of an aluminum hydroxide, and (3) bulk precipitation of an aluminosilicate phase. From the NMR speciation and water chemistry data, we calculate that 0.20 (±0.04) tetrahedral Al atoms nm⁻² sorb to the silica surface. Once the surface has sorbed roughly half of the total dissolved Al (∼8% site coverage), aluminum hydroxides and aluminosilicates precipitate from solution. These precipitation reactions are dependent upon solution pH and total dissolved silica concentration. We find that the Si:Al stoichiometry of the aluminosilicate precipitate is roughly 1:1 and suggest a chemical formula of NaAlSiO₄ in which Na⁺ acts as the charge compensating cation. For the adsorption of Al, we propose a surface-controlled reaction mechanism where Al sorbs as an inner-sphere coordination complex at the silica surface. Analogous to the hydrolysis of , we suggest that rapid deprotonation by surface hydroxyls followed by dehydration of ligated waters results in four-coordinate (>SiOH)₂Al(OH)₂ sites at the surface of amorphous silica.     

Isotopic studies of leaf water. Part 1: A physically based two-dimensional model for pine needles

Understanding of isotopic variations in leaf water is important for reconstruction of paleoclimate and assessment of global biochemical processes. We report here a study of isotopic distributions within a single needle of two pine species, Pinus resinosa Ait and Pinus strobes L., with the objective of understanding how isotopic compositions of leaf water are controlled by environmental and physiological variables. A 2D model was developed to simulate along-leaf isotopic variations and bulk leaf water isotopic compositions. In addition to variables common to all leaf water isotopic models, this 2D model also takes into account the specific geometry and dimensions of pine needles and the isotopic transport in xylem and mesophyll. The model can successfully simulate oxygen isotopic variations along a single needle and averaged over a leaf (bulk leaf water). The simulations suggest that isotopic composition of the bulk leaf water does not always depend only upon the average transpiration rate, which in turn raises questions about using leaf water isotopic values to estimate transpiration rates. An unsuccessful attempt to simulate along-needle hydrogen isotopic variations suggests that certain unknown biological process(es) may not have been incorporated into our 2D model, and if so, it calls for a reevaluation of all other models for hydrogen isotopic simulations of leaf water since they too lack these processes.Existing leaf water isotopic models are reviewed in this work. In particular, we evaluate the most frequently used model, the stomatal boundary layer model (also referred to as the Craig-Gordon model). We point out that discrepancy between the boundary layer model and the measured bulk leaf water seems to depend upon relative humidity. Using our 2D model, we show that this humidity dependency is a result of an interplay between environmental and physiological conditions: if the transpiration rate of plant leaves decreases with increasing relative humidity, our 2D model can reproduce the pattern of isotopic discrepancy between boundary layer model predictions and observations, enabling us to understand better the reason behind this discrepancy.     

The Fe-C system at 5 GPa and implications for Earth’s core

Earth’s core may contain C, and it has been suggested that C in the core could stabilize the formation of a solid inner core composed of Fe₃C. We experimentally examined the Fe-C system at a pressure of 5 GPa and determined the Fe-C phase diagram at this pressure. In addition, we measured solid metal/liquid metal partition coefficients for 17 trace elements and examined the partitioning behavior between Fe₃C and liquid metal for 14 trace elements. Solid metal/liquid metal partition coefficients are similar to those found in one atmosphere studies, indicating that the effect of pressure to 5 GPa is negligible. All measured Fe₃C/liquid metal partition coefficients investigated are less than one, such that all trace elements prefer the C-rich liquid to Fe₃C. Fe₃C/liquid metal partition coefficients tend to decrease with decreasing atomic radii within a given period. Of particular interest, our 5 GPa Fe-C phase diagram does not show any evidence that the Fe-Fe₃C eutectic composition shifts to lower C contents with increasing pressure, which is central to the previous reasoning that the inner core may be composed of Fe₃C.     

Limitations of the ferrozine method for quantitative assay of mineral systems for ferrous and total iron

The quantitative assay of clay minerals, soils, and sediments for Fe(II) and total Fe is fundamental to understanding biogeochemical cycles occurring therein. The commonly used ferrozine method was originally designed to assay extracted forms of Fe(II) from non-silicate aqueous systems. It is becoming, however, increasingly the method of choice to report the total reduced state of Fe in soils and sediments. Because Fe in soils and sediments commonly exists in the structural framework of silicates, extraction by HCl, as used in the ferrozine method, fails to dissolve all of the Fe. The phenanthroline (phen) method, on the other hand, was designed to assay silicate minerals for Fe(II) and total Fe and has been proven to be highly reliable. In the present study potential sources of error in the ferrozine method were evaluated by comparing its results to those obtained by the phen method. Both methods were used to analyze clay mineral and soil samples for Fe(II) and total Fe. Results revealed that the conventional ferrozine method under reports total Fe in samples containing Fe in silicates and gives erratic results for Fe(II). The sources of error in the ferrozine method are: (1) HCl fails to dissolve silicates and (2) if the analyte solution contains Fe³⁺, the analysis for Fe²⁺ will be photosensitive, and reported Fe(II) values will likely be greater than the actual amount in solution. Another difficulty with the ferrozine method is that it is tedious and much more labor intensive than the phen method. For these reasons, the phen method is preferred and recommended. Its procedure is simpler, takes less time, and avoids the errors found in the ferrozine method.     

Sulfur in peridotites and gabbros at Lost City (30°N, MAR): Implications for hydrothermal alteration and microbial activity during serpentinization

Variations in sulfur mineralogy and chemistry of serpentinized peridotites and gabbros beneath the Lost City Hydrothermal Field at the southern face of the Atlantis Massif (Mid-Atlantic Ridge, 30°N) were examined to better understand serpentinization and alteration processes and to study fluid fluxes, redox conditions, and the influence of microbial activity in this active, peridotite-hosted hydrothermal system. The serpentinized peridotites are characterized by low total sulfur contents and high bulk δ³⁴S values close to seawater composition. Low concentrations of ³⁴S-enriched sulfide phases and the predominance of sulfate with seawater-like δ³⁴S values indicate oxidation, loss of sulfide minerals and incorporation of seawater sulfate into the serpentinites. The predominance of pyrite in both serpentinites and gabbros indicates relatively high fO₂ conditions during progressive serpentinization and alteration, which likely result from high fluid fluxes during hydrothermal circulation and evolution of the Lost City system from temperatures of ∼250 to 150 °C. Sulfate and sulfide minerals in samples from near the base of hydrothermal carbonate towers at Lost City show δ³⁴S values that reflect the influence of microbial activity. Our study highlights the variations in sulfur chemistry of serpentinized peridotites in different marine environments and the influence of long-lived, moderate temperature peridotite-hosted hydrothermal system and high seawater fluxes on the global sulfur cycle.     

Lateral structural variability in zone of eocene island-arc-continent collision,Kamchatka

The lateral variability of structural elements in the collision zone of the Cretaceous-Paleocene Achaivayam-Valagin island arc with the northeastern Asian margin is considered. The similarity and difference of Eocene collision structural elements in the north and the south of Kamchatka are shown. In northern Kamchatka, the continent-arc boundary is traced along the Lesnaya-Vatyn Thrust Fault, which completed its evolution about 45 Ma ago. The thin, near-horizontal allochthon of this thrust, composed of island-arc rocks, overlies the deformed but unmetamorphosed terrigeneous sequences of the Asian margin. The general structure of this suture in the Kamchatka Isthmus and southern Koryakia is comparable with the uppermost subduction zone, where a thin lithospheric wedge overlaps intensely deformed sediments detached from the plunging plate. In southern Kamchatka (Malka Uplift of the Sredinny Range), the arc-continent collision started 55–53 Ma ago with thrusting of island-arc complexes over terrigenous rocks of continental margin. However, the thickness of the allochthon was much greater than in the north. Immediately after this event, both the autochthon and lower part of allochthon were deformed and subsided to a significant depth. This subsidence gave rise to metamorphism of both the autochthon (Kolpakov and Kamchatka groups, Kheivan Formation) and lower allochthon (Andrianovka and Khimka formations). The anomalously fast heating of the crust was most likely related to the ascent of asthenospheric masses due to slab breakoff, when the Eurasian Plate was plunging beneath the Achaivayam-Valagin arc.     

Late Precambrian rifting in the geological history of the western slope of the South Urals

The rift-related geodynamic setting of the Late Precambrian geological evolution on the western slope of the South Urals is reconstructed on the basis of localization of lithotectonic complexes of this age, their formation conditions, and the geochemistry of rocks. The Early Riphean stage comprises accumulation of coarse-clastic rocks intercalating with alkaline volcanic rocks of the Navysh Complex, which is a constituent of the Ai Formation, and emplacement of doleritic and picritic intrusions of the Shuida Complex and melanocratic dolerite and gabbrodolerite of the Yusha Complex. The Middle Riphean stage is characterized by wide-spread coarse-clastic terrigenous rocks of the Mashak Formation that intercalate with volcanic rocks of the bimodal basalt-rhyolite association, the Berdyaush pluton of rapakivi granite, the Kusa-Kopan layered intrusive complex, the Lapyshta Complex of dolerites and picrites, and numerous occurrences of gabbrodolerites. The terrigenous rocks of the Vendian stage include conglomerate, gravelstone, and sandstone of the Asha Group, while igneous rocks comprise alkaline volcanics of the Arsha Complex, alkali gabbroids of the Miseli Complex, and melanocratic syenite of the Avashla Complex. The geological evolution of the region is distinguished by local (failed or aborted) rifting. The occurrence of lithotectonic complexes is controlled by dynamic conditions of rifting. A certain inheritance in the evolution may be traced for the Early and Middle Riphean and partly for the Late Riphean and Vendian.