Diamond in metasedimentary crustal rocks from Pohorje,Eastern Alps: a window to deep continental subduction

We report the first finding of diamond and moissanite in metasedimentary crustal rocks of Pohorje Mountains (Slovenia) in the Austroalpine ultrahigh‐pressure (UHP) metamorphic terrane of the Eastern Alps. Microscopic observations and Raman spectroscopy show that diamond occurs in situ as inclusions in garnet, being heterogeneously distributed. Under the optical microscope, diamond‐bearing inclusions are of cuboidal to rounded shape and of pinkish, yellow to brownish colour. The Raman spectra of the investigated diamond show a sharp, first order peak of sp³‐bonded carbon, in most cases centred between 1332 and 1330 cm^?1, with a full width at half maximum between 3 and 5 cm^?1. Several spectra show Raman bands typical for disordered graphitic (sp²‐bonded) carbon. Detailed observations show that diamond occurs either as a monomineralic, single‐crystal inclusion or it is associated with SiC (moissanite), CO₂ and CH₄ in polyphase inclusions_. This rare record of diamond occurring with moissanite as fluid‐inclusion daughter minerals implies the crystallization of diamond and moissanite from a supercritical fluid at reducing conditions. Thermodynamic modelling suggests that diamond‐bearing gneisses attained P–T conditions of ≥3.5 GPa and 800–850 °C, similar to eclogites and garnet peridotites. We argue that diamond formed when carbonaceous sediment underwent UHP metamorphism at mantle depth exceeding 100 km during continental subduction in the Late Cretaceous (c. 95–92 Ma). The finding of diamond confirms UHP metamorphism in the Pohorje Mountains, the most deeply subducted part of Austroalpine units.     

Evaluation of NLDAS‐2 evapotranspiration against tower flux site observations

The North American Land Data Assimilation System project phase 2 (NLDAS‐2) has run four land surface models for a 30‐year (1979–2008) retrospective period. Land surface evapotranspiration (ET) is one of the most important model outputs from NLDAS‐2 for investigating land–atmosphere interaction or to monitor agricultural drought. Here, we evaluate hourly ET using in situ observations over the Southern Great Plains (Atmospheric Radiation Measurement/Cloud and Radiation Testbed network) for 1 January 1997–30 September 1999 and daily ET u‐sing in situ observations at the AmeriFlux network over the conterminous USA for an 8‐year period (2000–2007). The NLDAS‐2 models compare well against observations, with the National Centers for Environmental Prediction's Noah land surface model performing best, followed, in order, by the Variable Infiltration Capacity, Sacramento Soil Moisture Accounting, and Mosaic models. Daily evaluation across the AmeriFlux network shows that for all models, performance depends on season and vegetation type; they do better in spring and fall than in winter or summer and better for deciduous broadleaf forest and grasslands than for croplands or evergreen needleleaf forest. Copyright © 2014 John Wiley & Sons, Ltd.     

Insights into the complexity of crustal differentiation: K2O‐poor leucosomes within metasedimentary migmatites from the Southern Marginal Zone of the Limpopo Belt,South Africa

Differentiation of the continental crust is the result of complex interactions between a large number of processes, which govern partial melting of the deep crust, magma formation and segregation, and magma ascent to significantly higher crustal levels. The anatectic metasedimentary rocks exposed in the Southern Marginal Zone of the Limpopo Belt represent an unusually well‐exposed natural laboratory where the portion of these processes that operate in the deep crust can be directly investigated in the field. The formation of these migmatites occurred via absent incongruent melting reactions involving biotite, which produced cm‐ to m‐scale, K₂O‐poor garnet‐bearing stromatic leucosomes, with high Ca/Na ratios relative to their source rocks. Field investigation combined with geochemical analyses, and phase equilibrium modelling designed to investigate some aspects of disequilibrium partial melting show that the outcrop features and compositions of the leucosomes suggest several steps in their evolution: (1) Melting of a portion of the source, with restricted plagioclase availability due to kinetic controls, to produce a magma (melt + entrained peritectic minerals in variable proportions relative to melt); (2) Segregation of the magma at near peak metamorphic conditions into melt accumulation sites (MAS), also known as future leucosome; (3a) Re‐equilibration of the magma with a portion of the bounding mafic residuum via chemical diffusion (H₂O, K₂O), which triggers the co‐precipitation of quartz and plagioclase in the MAS; (3b) Extraction of melt‐dominated magma to higher crustal levels, leaving peritectic minerals entrained from the site of the melting reaction, and the minerals precipitated in the MASs to form the leucosome in the source. The key mechanism controlling this behaviour is the kinetically induced restriction of the amount of plagioclase available to the melting reaction. This results in elevated melt H₂O and K₂O and chemical potential gradient for these components across the leucosome/mafic residuum contact. The combination of all of these processes accurately explains the composition of the K₂O‐poor leucosomes. These findings have important implications for our understanding of melt segregation in the lower crust and minimum melt residency time which, according to the chemical modelling, is <5 years. We demonstrate that in some migmatitic granulites, the leucosomes constitute a type of felsic refractory residuum, rather than evidence of failed magma extraction. This provides a new insight into the ways that source heterogeneity may control anatexis.     

Constraints on the timing and conditions of high‐grade metamorphism,charnockite formation and fluid–rock interaction in the Trivandrum Block,southern India

Incipient charnockites have been widely used as evidence for the infiltration of CO₂‐rich fluids driving dehydration of the lower crust. Rocks exposed at Kakkod quarry in the Trivandrum Block of southern India allow for a thorough investigation of the metamorphic evolution by preserving not only orthopyroxene‐bearing charnockite patches in a host garnet–biotite felsic gneiss, but also layers of garnet–sillimanite metapelite gneiss. Thermodynamic phase equilibria modelling of all three bulk compositions indicates consistent peak‐metamorphic conditions of 830–925 °C and 6–9 kbar with retrograde evolution involving suprasolidus decompression at high temperature. These models suggest that orthopyroxene was most likely stabilized close to the metamorphic peak as a result of small compositional heterogeneities in the host garnet–biotite gneiss. There is insufficient evidence to determine whether the heterogeneities were inherited from the protolith or introduced during syn‐metamorphic fluid flow. U–Pb geochronology of monazite and zircon from all three rock types constrains the peak of metamorphism and orthopyroxene growth to have occurred between the onset of high‐grade metamorphism at c. 590 Ma and the onset of melt crystallization at c. 540 Ma. The majority of metamorphic zircon growth occurred during protracted melt crystallization between c. 540 and 510 Ma. Melt crystallization was followed by the influx of aqueous, alkali‐rich fluids likely derived from melts crystallizing at depth. This late fluid flow led to retrogression of orthopyroxene, the observed outcrop pattern and to the textural and isotopic modification of monazite grains at c. 525–490 Ma.     

Water uptake of riparian plants in the lower Lhasa River Basin,South Tibetan Plateau using stable water isotopes

Riparian plants can adapt their water uptake strategies based on climatic and hydrological conditions within a river basin. The response of cold-alpine riparian trees to changes in water availability is poorly understood. The Lhasa River is a representative cold-alpine river in South Tibet and an under-studied environment. Therefore, a 96 km section of the lower Lhasa River was selected for a study on the water-use patterns of riparian plants. Plant water, soil water, groundwater and river water were measured at three sites for δ¹⁸O and δ²H values during the warm-wet and cold-dry periods in 2018. Soil profiles differed in isotope values between seasons and with the distance along the river. During the cold-dry period, the upper parts of the soil profiles were significantly affected by evaporation. During the warm-wet period, the soil profile was influenced by precipitation infiltration in the upper reaches of the study area and by various water sources in the lower reaches. Calculations using the IsoSource model indicated that the mature salix and birch trees (Salix cheilophila Schneid. and Betula platyphylla Suk.) accessed water from multiple sources during the cold-dry period, whereas they sourced more than 70% of their requirement from the upper 60–80 cm of the soil profile during the warm-wet period. The model indicated that the immature rose willow tree (Tamarix ramosissima Ledeb) accessed 66% of its water from the surface soil during the cold-dry period, but used the deeper layers during the warm-wet period. The plant type was not the dominant factor driving water uptake patterns in mature plants. Our findings can contribute to strategies for the sustainable development of cold-alpine riparian ecosystems. It is recommended that reducing plantation density and collocating plants with different rooting depths would be conducive to optimal plant growth in this environment.     

Interpreting zirconium‐in‐rutile thermometric results

At first sight, experimental results and observations on rocks suggest that the Zr content in rutile, where equilibrated with quartz and zircon, should be a useful thermometer for metamorphic rocks. However, diffusion data for Zr in rutile imply that thermometry should not, for plausible rates of cooling, give the high temperatures commonly observed in high‐grade metamorphic rocks. It is suggested here that such observations can be accounted for by high‐T diffusive closure of Si in rutile, causing the interior of rutile grains to become insensitive to the thermometer equilibrium well above the temperature of Zr diffusive closure. Paired with comparatively slow grain boundary diffusion and problematic zircon nucleation, this allows for cases of Zr retention in rutile through temperatures where Zr is still diffusively mobile within rutile grains. Other observations that may be accounted for in this context are large inter‐grain ranges of rutile Zr contents uncorrelated with rutile grain size, and flat Zr profiles across individual rutile grains, counter to what would be expected from diffusive closure. A consequence is that it is unlikely that Zr‐in‐rutile thermometry will be useful for estimating rock cooling rates.     

Dissolved inorganic carbon evolution in neutral discharge from mine tailings piles

We measured the concentrations of dissolved inorganic carbon (DIC) and major ions and the stable carbon isotope ratios of DIC (δ¹³C_DIC) in two creeks discharging from carbonate‐rich sulphide‐containing mine tailings piles. Our aim was to assess downstream carbon evolution of the tailings discharge as it interacted with the atmosphere. The discharge had pH of 6.5–8.1 and was saturated with respect to carbonates. Over the reach of one creek, the DIC concentrations decreased by 1.1 mmol C/l and δ¹³C_DIC increased by ~4.0‰ 200 m from the seep source. The decrease in the DIC concentrations was concomitant with decreases in the partial pressure of CO_2(aq) because of the loss of excess CO_2(aq) from the discharge. The corresponding enrichment in the δ¹³C_DIC is because of kinetic isotope fractionation accompanying the loss of CO_2(g). Over the reach of the other creek, there was no significant decrease in the DIC concentrations or notable changes in the δ¹³C_DIC. The insignificant change in the DIC concentrations and the δ¹³C_DIC is because the first water sample was collected 160 m away from the discharge seep, not accessible during this research. In this case, most of the excess CO_2(aq) was lost before our first sampling station. Our results indicate that neutral discharges from tailings piles quickly lose excess CO_2(aq) to the atmosphere and the DIC becomes enrich in ¹³C. We suggest that a significant amount of carbon cycling in neutral discharges from tailings piles occur close to the locations where the discharge seeps to the surface. Copyright © 2015 John Wiley & Sons, Ltd.     

Oxygen isotope speedometry in granulite facies garnet recording fluid/melt–rock interaction (Sør Rondane Mountains,East Antarctica)

In situ analysis of a garnet porphyroblast from a granulite facies gneiss from Sør Rondane Mountains, East Antarctica, reveals discontinuous step‐wise zoning in phosphorus and large δ¹⁸O variations from the phosphorus‐rich core to the phosphorus‐poor rim. The gradually decreasing profile of oxygen isotope from the core (δ¹⁸O = ~15‰) to the rim (δ¹⁸O = ~11‰) suggests that the ¹⁸O/¹⁶O zoning was originally step‐wise, and modified by diffusion after the garnet rim formation at ~800°C and 0.8 GPa. Fitting of the ¹⁸O/¹⁶O data to the diffusion equation constrains a duration of the high‐T event (~800°C) to c. 0.5–40 Ma after the garnet rim formation. The low δ¹⁸O value of the garnet rim, together with the previously reported low δ¹⁸O values in metacarbonates, indicates regional infiltration, probably along a detachment fault, of low δ¹⁸O fluid/melt possibly derived from meta‐mafic to ultramafic rocks.     

Modelling the flow of self‐compacting concrete

A Lagrangian particle‐based method, smooth particle hydrodynamics (SPH), is used in this paper to model the flow of self‐compacting concretes (SCC) with or without short steel fibres. An incompressible SPH method is presented to simulate the flow of such non‐Newtonian fluids whose behaviour is described by a Bingham‐type model, in which the kink in the shear stress vs shear strain rate diagram is first appropriately smoothed out. The viscosity of the SCC is predicted from the measured viscosity of the paste using micromechanical models in which the second phase aggregates are treated as rigid spheres and the short steel fibres as slender rigid bodies. The basic equations solved in the SPH are the incompressible mass conservation and Navier–Stokes equations. The solution procedure uses prediction–correction fractional steps with the temporal velocity field integrated forward in time without enforcing incompressibility in the prediction step. The resulting temporal velocity field is then implicitly projected on to a divergence‐free space to satisfy incompressibility through a pressure Poisson equation derived from an approximate pressure projection. The results of the numerical simulation are benchmarked against actual slump tests carried out in the laboratory. The numerical results are in excellent agreement with test results, thus demonstrating the capability of SPH and a proper rheological model to predict SCC flow and mould‐filling behaviour. Copyright © 2010 John Wiley & Sons, Ltd.     

Influence of heterogeneity on unsaturated hydraulic properties: (1) local heterogeneity and scale effect

Spatial heterogeneity is ubiquitous in nature, which may significantly affect the soil hydraulic property curves. The models of a closed‐form functional relationship of soil hydraulic property curves (e.g. VG model or exponential model) are valid at point or local scale based on a point‐scale hydrological process, but how do scale effects of heterogeneity have an influence on the parameters of these models when the models are used in a larger scale process? This paper uses a two‐dimensional variably saturated flow and solute transport finite element model (VSAFT2) to simulate variations of pressure and moisture content in the soil flume under a constant head boundary condition. By changing different numerical simulation block sizes, a quantitative evaluation of parameter variations in the VG model, resulting from the scale effects, is presented. Results show that the parameters of soil hydraulic properties are independent of scale in homogeneous media. Parameters of α and n in homogeneous media, which are estimated by using the unsaturated hydraulic conductivity curve (UHC) or the soil water retention curve (WRC), are identical. Variations of local heterogeneities strongly affect the soil hydraulic properties, and the scale affects the results of the parameter estimations when numerical experiments are conducted. Furthermore, the discrepancy of each curve becomes considerable when moisture content becomes closer to a dry situation. Parameters estimated by UHC are totally different from the ones estimated by WRC. Copyright © 2012 John Wiley & Sons, Ltd.