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.     

Application of precise 142Nd/144Nd analysis of small samples to inclusions in diamonds (Finsch,South Africa) and Hadean Zircons (Jack Hills,Western Australia)

¹⁴⁶Sm–¹⁴²Nd and ¹⁴⁷Sm–¹⁴³Nd systematics were investigated in garnet inclusions in diamonds from Finsch (S. Africa) and Hadean zircons from Jack Hills (W. Australia) to assess the potential of these systems as recorders of early Earth evolution. The study of Finsch inclusions was conducted on a composite sample of 50 peridotitic pyropes with a Nd model age of 3.3 Ga. Analysis of the Jack Hills zircons was performed on 790 grains with ion microprobe ²⁰⁷Pb/²⁰⁶Pb spot ages from 3.95 to 4.19 Ga. Finsch pyropes yield 100 × ?¹⁴²Nd = ? 6 ± 12 ppm, ?¹⁴³Nd = ? 32.5, and ¹⁴⁷Sm/¹⁴⁴Nd = 0.1150. These results do not confirm previous claims for a 30 ppm ¹⁴²Nd excess in South African cratonic mantle. The lack of a ¹⁴²Nd anomaly in these inclusions suggests that isotopic heterogeneities created by early mantle differentiation were remixed at a very fine scale prior to isolation of the South African lithosphere. Alternatively, this result may indicate that only a fraction of the mantle experienced depletion during the first 400 Myr of its history. Analysis of the Jack Hills zircon composite yielded 100 × ?¹⁴²Nd = 8 ± 10 ppm, ?¹⁴³Nd = 45 ± 1, and ¹⁴⁷Sm/¹⁴⁴Nd = 0.5891. Back-calculation of this present-day ?¹⁴³Nd yields an unrealistic estimate for the initial ?¹⁴³Nd of ? 160 ?-units, clearly indicating post-crystallization disturbance of the ¹⁴⁷Sm–¹⁴³Nd system. Examination of ^146,147Sm–^142,143Nd data reveals that the Nd budget of the Jack Hills sample is dominated by non-radiogenic Nd, possibly contained in recrystallized zircon rims or secondary subsurface minerals. This secondary material is characterized by highly discordant U–Pb ages. Although the mass fraction of altered zircon is unlikely to exceed 5–10% of total sample, its high LREE content precludes a reliable evaluation of ¹⁴⁶Sm–¹⁴²Nd systematics in Jack Hills zircons.     

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.     

Geomorphological hazards related to deep dissolution phenomena in the Western Italian Alps: Distribution, assessment and interaction with human activities

Deep dissolution affects great part of soluble rocks (e.g. gypsum and anhydrite) of the Western Italian Alps. The related superficial phenomena (sinkholes, gravity-induced processes and a local worsening of geomechanical rock properties) are not limited to typical karsts landscape and cause slope instability also affecting populated sites and infrastructures. The paper aims to describe general characteristic of dissolution phenomena, to interpret their conditioning factors and evolutionary stages and to assess possible hazards due to their superficial effects.The search for evidences of deep dissolution leads to the selection of representative sites in the central part of the Western Italian Alps (Piemonte and Valle d'Aosta Region). Detailed geological and geomorphological studies have been used to classify the selected sites by type, size and variable state of activity. Very different evolutionary stages of dissolution phenomena have been interpreted by comparison of case-studies: some are early “embryonic”; others are more evolved, up to typical sinkholes, or even remodelled by other phenomena. Some cases show an extreme complexity in the interactions between corrosion phenomena and other geomorphic processes: slope deformations, from one side, and karst, fluvial and glacial phenomena, to the other. A wide range of movement rates on slope instabilities induced by deep dissolution have been estimated by topographic and geomorphic data. Geochemical data on removed rocks by dissolution indicate 0.4 mm/year values for local subsidence. Historical and technical data indicate low frequency of major dissolution-induced collapses, but highlight widespread damages to tunnels, roads and buildings, especially along slopes.     

Development and morphometry of sinkholes in coastal plains of Apulia, southern Italy. Preliminary sinkhole susceptibility assessment

Evolution of coastlines in karst areas may be strongly controlled by dissolution processes which favour the development of surface and subsurface landforms. The generation of caves in these environments is commonly favoured by the mixing between fresh and brackish waters. The sinkholes resulting from the upward propagation of the caves may interfere with the anthropogenic environment and cause damage to human elements (property and activities). To highlight the often underestimated importance of karst phenomena in coastal areas, we have analyzed a coastal stretch of Apulia, in southern Italy. The study area, covering an extension of about 6 km², is situated in the Ionian coast, and presents several interesting karst landforms that are generally connected to caves. Tens of sinkholes were mapped through field surveys, multi-year aerial-photographs (dating back to the 1940s) and archival research. We have performed a morphometric analysis of the sinkholes. The analysis describes the main parameters of the sinkholes (area, length, width, and depth), and the control exerted by the main discontinuity systems in the area. The detrimental effects derived from interaction between human environment and these karst landforms is also under consideration. A sinkhole susceptibility map, which may provide useful information for planners, developers and the insurance industry has eventually been produced through application of a decision tree model.     

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.     

Oxygen isotopic compositions of Allende Type C CAIs: Evidence for isotopic exchange during nebular melting and asteroidal metamorphism

In situ oxygen isotopic measurements of primary and secondary minerals in Type C CAIs from the Allende CV3 chondrite reveal that the pattern of relative enrichments and depletions of ¹⁶O in the primary minerals within each individual CAI are similar to the patterns observed in Types A and B CAIs from the same meteorite. Spinel is consistently the most ¹⁶O-rich (Δ¹⁷O = −25‰ to −15‰), followed by Al,Ti-dioside (Δ¹⁷O = −20‰ to −5‰) and anorthite (Δ¹⁷O = −15‰ to 0‰). Melilite is the most ¹⁶O-depleted primary mineral (Δ¹⁷O = −5‰ to −3‰). We conclude that the original melting event that formed Type C CAIs occurred in a ¹⁶O-rich (Δ¹⁷O  −20‰) nebular gas and they subsequently experienced oxygen isotopic exchange in a ¹⁶O-poor reservoir. At least three of these (ABC, TS26F1 and 93) experienced remelting at the time and place where chondrules were forming, trapping and partially assimilating ¹⁶O-poor chondrule fragments. The observation that the pyroxene is ¹⁶O-rich relative to the feldspar, even though the feldspar preceded it in the igneous crystallization sequence, disproves the class of CAI isotopic exchange models in which partial melting of a ¹⁶O-rich solid in a ¹⁶O-poor gas is followed by slow crystallization in that gas. For the typical (not associated with chondrule materials) Type C CAIs as well for as the Types A and B CAIs, the exchange that produced internal isotopic heterogeneity within each CAI must have occurred largely in the solid state. The secondary phases grossular, monticellite and forsterite commonly have similar oxygen isotopic compositions to the melilite and anorthite they replace, but in one case (CAI 160) grossular is ¹⁶O-enriched (Δ¹⁷O = −10‰ to −6‰) relative to melilite (Δ¹⁷O = −5‰ to −3‰), meaning that the melilite and anorthite must have exchanged its oxygen subsequent to secondary alteration. This isotopic exchange in melilite and anorthite likely occurred on the CV parent asteroid, possibly during fluid-assisted thermal metamorphism.     

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.     

Kinetics of Cd sorption, desorption and fixation by calcite: A long-term radiotracer study

Time-dependent sorption and desorption of Cd on calcite was studied over 210 days utilizing ¹⁰⁹Cd as a tracer to distinguish between ‘labile’ and ‘non-labile’ forms of sorbed Cd. Stabilizing the calcite suspensions for 12 months under atmospheric P_CO2 and controlled temperature was necessary to reliably follow Cd dynamics following initial sorption. Results revealed time-dependant Cd sorption and marked desorption hysteresis by calcite under environmentally relevant conditions. Data obtained were fitted to a first-order kinetic model and a concentric shell diffusion model. Both models described the progressive transfer of Cd²⁺ to a less reactive form within calcite and subsequent desorption of Cd subject to different initial contact times. The kinetic model provided a better fit to the combined sorption and desorption data (R² = 0.992). It differentiates between two ‘pools’ of sorbed Cd²⁺ on calcite, ‘labile’ and ‘non-labile’, in which labile sorbed Cd is in immediate equilibrium with the free Cd²⁺ ion activity in solution whereas non-labile Cd is kinetically restricted. For the diffusion model (R² = 0.959), the rate constants describing Cd dynamics in calcite produced a half-life for Cd desorption of ∼175 d, for release to a ‘zero-sink’ solution. Results from this study allow comment on the likely mechanisms occurring at the calcite surface following long-term Cd sorption.