Experimental and theoretical study of stability of dense hydrous magnesium silicates in the deep upper mantle

Fluid-undersaturated experiments were conducted to determine the phase relations in the simplified peridotite system MgO-SiO₂-H₂O (MSH) at 11.0-14.5 GPa and 800-1400 °C. Stability relations of dense hydrous magnesium silicates (DHMSs) under fluid-undersaturated conditions were experimentally examined. From the fluid-absent experimental results, we retrieved thermodynamic data for clinohumite, phase A, phase E, and hydrous wadsleyite, consistent with the published data set for dry mantle minerals. With this new data set, we have calculated phase equilibria in the MSH system including dehydration reactions. The dehydration reactions calculated with lower water activities of 0.68-0.60 match the fluid-present experiments of this study above 11.0 GPa and 1000 °C, indicating that considerable amounts of silicate component were dissolved into the fluid phase. The calculated phase equilibria illustrate the stability of the post-antigorite phase A-bearing assemblages. In the cold subducting slab peridotite, phase A + enstatite assemblage survives into the transition zone, whereas phase A + forsterite + enstatite assemblage forms hydrous wadsleyite at a much shallower depth of about 360-km. The slab is subducted with no dehydration reactions occurring when entering the transition zone. The phase equilibria also show the high temperature stability of phase E. Phase E is stable up to 1200 °C at 13.5 GPa, a plausible condition in the mantle of relatively low temperature, i.e., beneath subduction zones. Phase E is a possible water reservoir in the mantle as well as wadsleyite and ringwoodite.     

长沙地震台北东方向深震sP震相讨论

长沙地震台北东方向深源地震的记录中,发生于日本海边缘和中朝俄接界区域的深源地震存在sP震相,它超前于S波,由于震相能量较大,容易与S波相混淆。本文讨论了如何对它进行辨别。     

Analysis of Short-Period Rayleigh Waves Recorded in the Bohemian Massif Area During Celebration 2000 Experiment

The aim of this paper is to show the application of short-period surface waves recorded during deep seismic sounding experiment for constraining shallow velocity structure of the crust. Phase velocity of fundamental mode Rayleigh waves, observed along the CELEBRATION 2000 experiment profile CEL09, were obtained by a p-ω method and has been subsequently inverted for one-dimensional shear velocity models for the top 2 km. Multiple filter technique applied to one shot gather was used to carry out a joint inversion of phase and group velocity data and to provide γ_R data to be used for Q_β inversion. Validity of obtained V_S and Q_β models was confirmed by the reflectivity method. Noticeably, no clear dispersive wawes were observed in the Tepla-Barrandian Unit. Quasi-2D model based on the individual 1D V_S models is well correlated with the surface geology. Lower V_S are observed in the Saxothuringian Zone in comparison to the Moldanubian Zone. In the vicinity of the Central Bohemian and Moldanubian Plutons, the near-surface V_S values are relatively low, but below 1 km depth, they are higher than in surrounding areas. We interpret it as the result of the weathering and cracks within the granitoid rocks.     

Correcting the Smoothing Effect of Ordinary Kriging Estimates

The smoothing effect of ordinary kriging is a well-known dangerous effect associated with this estimation technique. Consequently kriging estimates do not reproduce both histogram and semivariogram model of sample data. A four-step procedure for correcting the smoothing effect of ordinary kriging estimates is shown to be efficient for the reproduction of histogram and semivariogram without loss of local accuracy. Furthermore, this procedure provides a unique map sharing both local and global accuracies. Ordinary kriging with a proper correction for smoothing effect can be revitalized as a reliable estimation method that allows a better use of the available information.     

指数平滑技术在斜坡位移预测中的应用

应用趋势型数据指数平滑模型预测了链子崖危岩体监测点GA的位移量.根据其观测数据呈线性趋势的特点,选取趋势型二次指数平滑的线性预测公式和合适的平滑常数值进行了计算预测.预测结果表明位移预测值与实际观测值之间的误差很小,说明该模型可很好地应用于斜坡变形位移的预测.     

Polymorphic phase transition in Superhydrous Phase B

We synthesized superhydrous phase B (shy-B) at 22 GPa and two different temperatures: 1200°C (LT) and 1400°C (HT) using a multi-anvil apparatus. The samples were investigated by transmission electron microscopy (TEM), single crystal X-ray diffraction, Raman and IR spectroscopy. The IR spectra were collected on polycrystalline thin-films and single crystals using synchrotron radiation, as well as a conventional IR source at ambient conditions and in situ at various pressures (up to 15 GPa) and temperatures (down to −180°C). Our studies show that shy-B exists in two polymorphic forms. As expected from crystal chemistry, the LT polymorph crystallizes in a lower symmetry space group (Pnn2), whereas the HT polymorph assumes a higher symmetry space group (Pnnm). TEM shows that both modifications consist of nearly perfect crystals with almost no lattice defects or inclusions of additional phases. IR spectra taken on polycrystalline thin films exhibit just one symmetric OH band and 29 lattice modes for the HT polymorph in contrast to two intense but asymmetric OH stretching bands and at least 48 lattice modes for the LT sample. The IR spectra differ not only in the number of bands, but also in the response of the bands to changes in pressure. The pressure derivatives for the IR bands are higher for the HT polymorph indicating that the high symmetry form is more compressible than the low symmetry form. Polarized, low-temperature single-crystal IR spectra indicate that in the LT-polymorph extensive ordering occurs not only at the Mg sites but also at the hydrogen sites.     

Tetrahedral Plots of the Phase Relations for Basalts

The phase relations of quaternary systems are generally represented by projections onto ternary compositional planes. Such projections often obscure relationships that would only be evident in a three-dimensional tetrahedral plot. The tetrahedral plot requires that compositions of the minerals and melts be transformed into Cartesian coordinates. It is shown here how this transformation is carried out. The application is demonstrated by tetrahedral plots of experimental melt compositions of partially molten lherzolite. Furthermore, the plot can be used to evaluate whether or not a particular basaltic composition represents a primary melt. The methods are applicable to any four-component system.     

Fluid generation and evolution during exhumation of deeply subducted UHP continental crust: Petrogenesis of composite granite–quartz veins in the Sulu belt,China

Composite granite–quartz veins occur in retrogressed ultrahigh pressure (UHP) eclogite enclosed in gneiss at General's Hill in the central Sulu belt, eastern China. The granite in the veins has a high‐pressure (HP) mineral assemblage of dominantly quartz+phengite+allanite/epidote+garnet that yields pressures of 2.5–2.1 GPa (Si‐in‐phengite barometry) and temperatures of 850–780°C (Ti‐in‐zircon thermometry) at 2.5 GPa (~20°C lower at 2.1 GPa). Zircon overgrowths on inherited cores and new grains of zircon from both components of the composite veins crystallized at c. 221 Ma. This age overlaps the timing of HP retrograde recrystallization dated at 225–215 Ma from multiple localities in the Sulu belt, consistent with the HP conditions retrieved from the granite. The ε_Hf(t) values of new zircon from both components of the composite veins and the Sr–Nd isotope compositions of the granite consistently lie between values for gneiss and eclogite, whereas δ¹⁸O values of new zircon are similar in the veins and the crustal rocks. These data are consistent with zircon growth from a blended fluid generated internally within the gneiss and the eclogite, without any ingress of fluid from an external source. However, at the peak metamorphic pressure, which could have reached 7 GPa, the rocks were likely fluid absent. During initial exhumation under UHP conditions, exsolution of H₂O from nominally anhydrous minerals generated a grain boundary supercritical fluid in both gneiss and eclogite. As exhumation progressed, the volume of fluid increased allowing it to migrate by diffusing porous flow from grain boundaries into channels and drain from the dominant gneiss through the subordinate eclogite. This produced a blended fluid intermediate in its isotope composition between the two end‐members, as recorded by the composite veins. During exhumation from UHP (coesite) eclogite to HP (quartz) eclogite facies conditions, the supercritical fluid evolved by dissolution of the silicate mineral matrix, becoming increasingly solute‐rich, more ‘granitic’ and more viscous until it became trapped. As crystallization began by diffusive loss of H₂O to the host eclogite concomitant with ongoing exhumation of the crust, the trapped supercritical fluid intersected the solvus for the granite–H₂O system, allowing phase separation and formation of the composite granite–quartz veins. Subsequently, during the transition from HP eclogite to amphibolite facies conditions, minor phengite breakdown melting is recorded in both the granite and the gneiss by K‐feldspar+plagioclase+biotite aggregates located around phengite and by K‐feldspar veinlets along grain boundaries. Phase equilibria modelling of the granite indicates that this late‐stage melting records P–T conditions towards the end of the exhumation, with the subsolidus assemblage yielding 0.7–1.1 GPa at <670°C. Thus, the composite granite–quartz veins represent a rare example of a natural system recording how the fluid phase evolved during exhumation of continental crust. The successive availability of different fluid phases attending retrograde metamorphism from UHP eclogite to amphibolite facies conditions will affect the transport of trace elements through the continental crust and the role of these fluids as metasomatic agents interacting with the mantle wedge in the subduction channel.     

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.     

Low‐P and high‐T metamorphism of basalts: Insights from the Sudbury impact melt sheet aureole and thermodynamic modelling

Low‐pressure and high‐temperature (LP–HT) metamorphism of basaltic rocks, which occurs globally and throughout geological time, is rarely constrained by forward phase equilibrium modelling, yet such calculations provide valuable supplementary thermometric information and constraints on anatexis that are not possible to obtain from conventional thermometry. Metabasalts along the southern margin of the Sudbury Igneous Complex (SIC) record evidence of high‐grade contact metamorphism involving partial melting and melt segregation. Peak metamorphic temperatures reached at least ~925°C at ~1–3 kbar near the SIC contact. Preservation of the peak mineral assemblage indicates that most of the generated melt escaped from these rocks leaving a residuum characterized by a plagioclase–orthopyroxene–clinopyroxene–ilmenite‐magnetite±melt assemblage. Peak temperatures reached ~875°C up to 500 m from the SIC lower contact, which marks the transition to metabasalts that only experienced incipient partial melting without melt loss. Metabasalts ~500 to 750 m from the SIC contact are characterized by a similar two‐pyroxene mineral assemblage, but typically contain abundant hornblende that overgrew clino‐ and orthopyroxene along an isobaric cooling path. Metabasalts ~750 to 1,000 m from the SIC contact are characterized by a hornblende–plagioclase–quartz–ilmenite assemblage indicating temperatures up to ~680°C. Mass balance and phase equilibria calculations indicate that anatexis resulted in 10–20% melt generation in the inner ~500 m of the aureole, with even higher degrees of melting towards the contact. Comparison of multiple models, experiments, and natural samples indicates that modelling in the Na₂O–CaO–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (NCFMASHTO) system results in the most reliable predictions for the temperature of the solidus. Incorporation of K₂O in the most recent amphibole solution model now successfully predicts dehydration melting by the coexistence of high‐Ca amphibole and silicate melt at relatively low pressures (~1.5 kbar). However, inclusion of K₂O as a system component results in prediction of the solidus at too low a temperature. Although there are discrepancies between modelling predictions and experimental results, this study demonstrates that the pseudosection approach to mafic rocks is an invaluable tool to constrain metamorphic processes at LP–HT conditions.