Origins and evolution of rhyolitic magmas in the central Snake River Plain: insights from coupled high-precision geochronology,oxygen isotope,and hafnium isotope analyses of zircon

We present new high-precision CA-ID-TIMS and in situ U–Pb ages together with Hf and O isotopic analyses (analyses performed all on the same grains) from four tuffs from the 15?10 Ma Bruneau–Jarbidge center of the Snake River Plain and from three rhyolitic units from the Kimberly borehole in the neighboring 10?6 Ma Twin Falls volcanic center. We find significant intrasample diversity in zircon ages (ranges of up to 3 Myr) and in δ¹⁸O (ranges of up to 6‰) and ε_Hf (ranges of up to 24 ε units) values. Zircon rims are also more homogeneous than the associated cores, and we show that zircon rim growth occurs faster than the resolution of in situ dating techniques. CA-ID-TIMS dating of a subset of zircon grains from the Twin Falls samples reveals complex crystallization histories spanning 10⁴–10⁶ years prior to some eruptions, suggesting that magma genesis was characterized by the cyclic remelting of buried volcanic rocks and intrusions associated with previous magmatic episodes. Age-dependent trends in zircon isotopic compositions show that rhyolite production in the Yellowstone hotspot track is driven by the mixing of mantle-derived melts (normal δ¹⁸O and ε_Hf) and a combination of Precambrian basement rock (normal δ¹⁸O and ε_Hf down to ??60) and shallow Mesozoic and Cenozoic age rocks, some of which are hydrothermally altered (to low δ¹⁸O values) by earlier stages of Snake River Plain magmatism. These crustal melts hybridize with juvenile basalts and rhyolites to produce the erupted rhyolites. We also observe that the Precambrian basement rock is only an important component in the erupted magmas in the first eruption at each caldera center, suggesting that the accumulation of new intrusions quickly builds an upper crustal intrusive body which is isolated from the Precambrian basement and evolves towards more isotopically juvenile and lower-δ¹⁸O compositions over time.     

Apatite-hosted melt inclusions from the Panzhihua gabbroic-layered intrusion associated with a giant Fe–Ti oxide deposit in SW China: insights for magma unmixing within a crystal mush

The ～?2-km-thick Panzhihua gabbroic-layered intrusion in SW China is unusual because it hosts a giant Fe–Ti oxide deposit in its lower zone. The deposit consists of laterally extensive net-textured and massive Fe–Ti oxide ore layers, the thickest of which is ～?60 m. To examine the magmatic processes that resulted in the Fe enrichment of parental high-Ti basaltic magma and the formation of thick, Fe–Ti oxide ore layers, we carried out a detailed study of melt inclusions in apatite from a ～?500-m-thick profile of apatite-bearing leucogabbro in the middle zone of the intrusion. The apatite-hosted melt inclusions are light to dark brown in color and appear as polygonal, rounded, oval and negative crystal shapes, which range from ～?5 to ～?50 µm in width and from ～?5 to ～?100 µm in length. They have highly variable compositions and show a large and continuous range of SiO₂ and FeO_t with contrasting end-members; one end-member being Fe-rich and Si-poor (40.2 wt% FeO_t and 17.7 wt% SiO₂) and the other being Si-rich and Fe-poor (74.0 wt% SiO₂ and 1.20 wt% FeO_t). This range in composition may be attributed to entrapment of the melt inclusions over a range of temperature and may reflect the presence of µm-scale and immiscible Fe-rich and Si-rich components in different proportions. Simulating results for the motion of Si-rich droplets within a crystal mush indicate that Si-rich droplets would be separated from Fe-rich melt and migrate upward due to density differences in the interstitial liquid when the magma unmixed. Migration of the Si-rich, immiscible liquid component from the interstitial liquid caused the remaining Fe-rich melt in the lower part to react with plagioclase primocrysts (An_59–60), as evidenced by fine-grained lamellar intergrowth of An-rich plagioclase (An_79–84)?+?clinopyroxene in the oxide gabbro of the lower zone. Therefore, magma unmixing within a crystal mush, combined with gravitationally driven loss of the Si-rich component, resulted in the formation of Fe-rich, melagabbro and major Fe–Ti oxide ores in the lower part and Si-rich, leucogabbro in the upper part of the intrusion.     

Magma evolution beneath Bequia,Lesser Antilles,deduced from petrology of lavas and plutonic xenoliths

Extrusive and intrusive igneous rocks represent different parts of a magmatic system and ultimately provide complementary information about the processes operating beneath volcanoes. To shed light on such processes, we have examined and quantified the textures and mineral compositions of plutonic and cumulate xenoliths and lavas from Bequia, Lesser Antilles arc. Both suites contain assemblages of iddingsitized olivine, plagioclase, clinopyroxene and spinel with rare orthopyroxene and ilmenite. Mineral zoning is widespread, but more protracted in lavas than xenoliths. Plagioclase cores and olivine have high anorthite (An?≤?98) and low forsterite (Fo?≤?84) compositions respectively, implying crystallisation from a hydrous mafic melt that was already fractionated. Xenolith textures range from adcumulate to orthocumulate with variable mineral crystallisation sequences. Textural criteria are used to organize the xenoliths into six groups. Amphibole, notably absent from lavas, is a common feature of xenoliths, together with minor biotite and apatite. Bulk compositions of xenoliths deviate from the liquid line of descent of lavas supporting a cumulate origin with varying degrees of reactive infiltration by evolved hydrous melts, preserved as melt inclusions in xenolith crystals. Volatile saturation pressures in melt inclusions indicate cumulate crystallization over a 162–571 MPa pressure range under conditions of high dissolved water contents (up to 7.8 wt% H₂O), consistent with a variety of other thermobarometric estimates. Phase assemblages of xenoliths are consistent with published experimental data on volatile-saturated low-magnesium and high-alumina basalts and basaltic andesite from the Lesser Antilles at pressures of 200–1000 MPa, temperatures of 950–1050 °C and dissolved H₂O contents of 4–7 wt%. Once extracted from mid-crustal mushes, residual melts ascend to higher levels and undergo H₂O-saturated crystallization in shallow, pre-eruptive reservoirs to form phenocrysts and glomerocrysts. The absence of amphibole from lavas reflects instability at low pressures, whereas its abundance in xenoliths testifies to its importance in mid-crustal differentiation processes. A complex, vertically extensive (6 to at least 21 km depth) magmatic system is inferred beneath Bequia. Xenoliths represent fragments of the mush incorporated into ascending magmas. The widespread occurrence of evolved melts in the mush, but the absence of erupted evolved magmas, in contrast to islands in the northern Lesser Antilles, may reflect the relative immaturity of the Bequia magmatic system.     

Phase-equilibrium geobarometers for silicic rocks based on rhyolite-MELTS. Part 4: Plagioclase,orthopyroxene, clinopyroxene,glass geobarometer,and application to Mt. Ruapehu,New Zealand

A new phase equilibria geobarometer determines magmatic storage and crystallization conditions, including pressure, temperature, oxygen fugacity (\({f_{{{\text{o}}_2}}}\)), and the presence of a fluid phase for glass-bearing rocks containing the assemblage plagioclase?+?pyroxene(s). This newly developed geobarometer can better constrain crystallization conditions of shallow (<?500 MPa; <~?20 km), glass-bearing andesites to dacites. The geobarometer utilizes rhyolite-MELTS to determine crystallization conditions in natural pumice and scoria samples. The validity of the geobarometer is tested by comparing it to results from experiments. Uncertainties are assessed using Monte Carlo simulations. We apply the geobarometer to the plag?+?opx?+?cpx-bearing system of Mt. Ruapehu, in the southern Taupo Volcanic Zone, New Zealand. The samples from Mt. Ruapehu are tested from ~?5 to ~?400 MPa and from super-liquidus to 90% crystalline (~ 1200 to ~ 700 °C). Mt. Ruapehu serves as a methodological testing ground for the geobarometer, and results from our geobarometer agree with recent Mt. Ruapehu studies. Results show a distribution of crystallization pressures ranging from 50 to 150 MPa (~?2.0 to 5.9 km) for different eruptions, with modes of 110 MPa (~ 4.3 km) and 130 MPa (~ 5.1 km). These are consistent with field interpretations of different eruptive styles based on juvenile clast textures and previous knowledge of the magma plumbing system. Mt. Ruapehu magmas are fluid saturated, with \({f_{{{\text{o}}_2}}}\) of ΔQFM ~ + 1 (NNO).     

Crystal nucleation and growth produced by continuous decompression of Pinatubo magma

High-temperature decompression experiments demonstrate that crystal textures preserve a record of the style and rate of magmatic ascent. To reinforce this link, we performed a suite of isothermal decompression experiments using starting material from the climactic 1991 Pinatubo eruption. We decompressed experiments from 220 MPa to final, quench pressures of 75 or 30 MPa using continuous decompression rates of 100, 30, 10, 3, 1, and 0.3 MPa h^?1. Amphibole, clinopyroxene, and plagioclase crystallized during the experiments, with plagioclase microlites dominating the assemblage. Total microlite number densities range from 10^7.6±0.4 up to 10^8.2±0.2 cm^?3, with plagioclase accounting for up to 65% of the total number. Plagioclase microlite area increased systematically from 19?±?8 to 937?±?487 µm² with increasing experiment duration. Our textures provide time-integrated records of crystal kinetics. Average nucleation and areal growth rates of plagioclase are highest in the fastest decompressions (~?10^7.5 cm^?3 h^?1 and 10.1?±?4.1 µm² h^?1, respectively) and more than an order of magnitude lower in the slowest experiments (~?10^5.5 cm^?3 h^?1 and 0.8?±?0.2 µm² h^?1, respectively). Both nucleation and growth rates are highest at high degrees of disequilibrium. We find that peak supersaturation-dependent instantaneous rates are generally more than an order of magnitude faster than average rates. We use those instantaneous nucleation and growth rates to introduce an iterative model to evaluate the effects of different decompression rates, decompression paths (continuous, single-step or multistep), and the presence of phenocrysts on final crystallinity and microlite size distribution.     

Zoning and exsolution in alkali feldspars from Laacher See volcano (Western Germany): constraints on temperature history prior to eruption

Compositional zoning and exsolution patterns of alkali feldspars in carbonatite-bearing cognate syenites from the 6.3 km³ (D.R.E) phonolitic Laacher See Tephra (LST) deposit in western Germany (12.9 ka) are reported. These rocks represent the cooler outer portion and crystal-rich products of a cooling magma reservoir at upper crustal levels. Major and trace-element difference between cores and rims in sanidine crystals represent two generations of crystal growth separated by unmixing of a carbonate melt. Trace-element differences measured by LA–ICP–MS are in accordance with silicate–carbonate unmixing. Across the core–rim boundary, we extracted gray-scale profiles from multiple accumulations of back-scattered electron images. Gray scales directly represent K/Na ratios owing to low concentrations of Ba and Sr (<?30 ppm). Diffusion gradients are modeled to solve for temperature using known pre-eruptive U–Th zircon ages (0–20 ky) of each sample (Schmitt et al., J Petrol 51:1053–1085, 2010). Estimated temperatures range from 630 °C to 670 °C. For the exsolution boundaries, a diffusive homogenization model is constrained by the solvus temperature of ~ 712^_725 °C and gives short time scales of only 15–50 days. Based on our results, we present a model for the temperature–time history of these rocks. The model also constrains the thermal variation across the cooling crystal-rich carapace of the magma reservoir over 20 ka and suggests a thermal reactivation of cumulates, the cooling carapace, and probably the entire system only a few years prior to the explosive eruption of the remaining molten core of the phonolitic magma reservoir.     

The accuracy of standard enthalpies and entropies for phases of petrological interest derived from density-functional calculations

The internal energies and entropies of 21 well-known minerals were calculated using the density functional theory (DFT), viz. kyanite, sillimanite, andalusite, albite, microcline, forsterite, fayalite, diopside, jadeite, hedenbergite, pyrope, grossular, talc, pyrophyllite, phlogopite, annite, muscovite, brucite, portlandite, tremolite, and CaTiO₃–perovskite. These thermodynamic quantities were then transformed into standard enthalpies of formation from the elements and standard entropies enabling a direct comparison with tabulated values. The deviations from reference enthalpy and entropy values are in the order of several kJ/mol and several J/mol/K, respectively, from which the former is more relevant. In the case of phase transitions, the DFT-computed thermodynamic data of involved phases turned out to be accurate and using them in phase diagram calculations yields reasonable results. This is shown for the Al₂SiO₅ polymorphs. The DFT-based phase boundaries are comparable to those derived from internally consistent thermodynamic data sets. They even suggest an improvement, because they agree with petrological observations concerning the coexistence of kyanite?+?quartz?+?corundum in high-grade metamorphic rocks, which are not reproduced correctly using internally consistent data sets. The DFT-derived thermodynamic data are also accurate enough for computing the P–T positions of reactions that are characterized by relatively large reaction enthalpies (>?100 kJ/mol), i.e., dehydration reactions. For reactions with small reaction enthalpies (a few kJ/mol), the DFT errors are too large. They, however, are still far better than enthalpy and entropy values obtained from estimation methods.     

高效液相色谱-电感耦合等离子体质谱联用检测土壤中的无机硒形态

土壤样品中亚硒酸盐Se(Ⅳ)和硒酸盐Se(Ⅵ)的形态分析中,提取剂的选择和检测方法是技术的关键。以往的提取剂容易导致硒形态发生转变或无法同时提取Se(Ⅳ)和Se(Ⅵ),常用的氢化物发生原子荧光光谱法无法直接测定Se(Ⅵ),而是通过差减法得出Se(Ⅵ)含量。本文对比了不同提取剂的提取能力,确定使用0. 1 mol/L氢氧化钠溶液作为提取剂,在55℃超声萃取土壤样品30 min,提取液经高效液相色谱分离,电感耦合等离子体质谱检测,建立了土壤中Se(Ⅳ)和Se(Ⅵ)的形态分析方法。采用Hamilton PRP X-100色谱柱,以6 mmol/L柠檬酸为流动相,pH=5. 5,在8 min内可完全分离Se(Ⅳ)和Se(Ⅵ),两者的检出限分别为0. 15μg/L、0. 16μg/L,线性相关系数(r~2)均大于0. 999。以土壤为基体进行加标回收试验,Se(Ⅳ)和Se(Ⅵ)的回收率在84. 2%～95. 8%之间,相对标准偏差为1. 4%～5. 3%(n=6)。该方法简单快速,具有良好的精密度和准确度,适用于土壤中无机硒的形态分析。     

丹阳天王寺边坡稳定性评价及防治措施

针对苏南地区典型露采边坡的地层岩性及地质构造特点,以极限平衡理论为基础,运用GeoStudio软件中的SLOPE/W模块和VADOSE/W模块计算分析江苏丹阳天王寺边坡在天然工况、暴雨工况及地震工况下的稳定性。结果表明:天然工况下,潜在滑动面分布范围最小,稳定性最高,存在滑动的危险性;暴雨工况下,潜在滑动面分布范围有所增大,稳定性降低,边坡滑动的可能性较大;地震工况下,潜在滑动面分布范围最广,稳定性最低,边坡滑动的可能性很大。基于上述结果,提出防治措施。