CrystalMom: a new model for the evolution of crystal size distributions in magmas with the quadrature-based method of moments

Nucleation and growth of crystals, and the resulting crystal size distribution, play a fundamental role in controlling the physical properties of magmas and consequently the dynamics of the eruptions. In the past decades, laboratory experiments demonstrated that size and shape of crystals strongly control the physical properties of magma and lava. Additionally, natural and experimental samples are usually characterized in terms of their crystal size distribution to link it with physical processes that are not directly observable, such as cooling or decompression mechanisms. In this paper, we present CrystalMoM, a new predictive model, based on the quadrature-based method of moments, developed for studying the kinetic of crystallization in volcanic systems. The quadrature-based method of moments, well established in the field of chemical engineering, represents a mesoscale modelling approach that rigorously simulates the space–time evolution of a distribution of particles, by considering its moments. The method is applied here, for the first time, for studying the equilibrium/disequilibrium crystallization in magma, modelling the temporal evolution of the moments of a crystal size distribution. The model, verified against numerical and experimental data, represents a valuable tool to infer the cooling and decompression rates from the crystal size distribution observed in natural samples.     

Abyssal and hydrated mantle wedge serpentinised peridotites: a comparison of the 15$$^\circ$$20$$'$$′N fracture zone and New Caledonia serpentinites

Subduction of serpentinised mantle transfers oxidised and hydrated mantle lithosphere into the Earth, with consequences for the oxidation state of sub-arc mantle and the genesis of arc-related ore deposits. The role of subducted serpentinised mantle lithosphere in earth system processes is uncertain because subduction fluxes are poorly constrained. Most subducted serpentinised mantle is serpentinised on the ocean floor settings. Yet this material is poorly represented in the literature because it is difficult to access. Large volumes of accessible serpentinite are available in ophiolite complexes, and most interpretations of subduction fluxes associated with ultramafic rocks are based on ophiolite studies. Seafloor and ophiolite serpentinisation can occur under different conditions, so it is necessary to assess if ophiolite serpentinites are a good proxy for seafloor serpentinites. Serpentinites sampled during ODP cruise 209 were compared with serpentinites from New Caledonia. The ODP209 serpentinites were serpentinised by modified seawater in a shallow hydrothermal seafloor setting. The New Caledonia serpentinites were serpentinised in a mantle wedge setting by slab-derived fluids, with possible contributions from oceanic serpentinisation and post-obduction serpentinisation. Petrological, whole rock and mineralogical analyses were combined to compare the two sample sets. Petrologically, the evolution of serpentinisation was close to identical in the two environments. However, more oxidised iron, Cl, S and C is present in serpentine from the ODP209 serpentinites relative to the New Caledonia serpentinites. Given these observations, the use of serpentinites from different geodynamic settings as a proxy for abyssal serpentinites from spreading settings must be undertaken with caution.     

Monazite trumps zircon: applying SHRIMP U–Pb geochronology to systematically evaluate emplacement ages of leucocratic,low-temperature granites in a complex Precambrian orogen

Although zircon is the most widely used geochronometer to determine the crystallisation ages of granites, it can be unreliable for low-temperature melts because they may not crystallise new zircon. For leucocratic granites U–Pb zircon dates, therefore, may reflect the ages of the source rocks rather than the igneous crystallisation age. In the Proterozoic Capricorn Orogen of Western Australia, leucocratic granites are associated with several pulses of intracontinental magmatism spanning ~800 million years. In several instances, SHRIMP U–Pb zircon dating of these leucocratic granites either yielded ages that were inconclusive (e.g., multiple concordant ages) or incompatible with other geochronological data. To overcome this we used SHRIMP U–Th–Pb monazite geochronology to obtain igneous crystallisation ages that are consistent with the geological and geochronological framework of the orogen. The U–Th–Pb monazite geochronology has resolved the time interval over which two granitic supersuites were emplaced; a Paleoproterozoic supersuite thought to span ~80 million years was emplaced in less than half that time (1688–1659 Ma) and a small Meso- to Neoproterozoic supersuite considered to have been intruded over ~70 million years was instead assembled over ~130 million years and outlasted associated regional metamorphism by ~100 million years. Both findings have consequences for the duration of associated orogenic events and any estimates for magma generation rates. The monazite geochronology has contributed to a more reliable tectonic history for a complex, long-lived orogen. Our results emphasise the benefit of monazite as a geochronometer for leucocratic granites derived by low-temperature crustal melting and are relevant to other orogens worldwide.     

Stability of phlogopite in ultrapotassic kimberlite-like systems at 5.5–7.5 GPa

Hydrous K-rich kimberlite-like systems are studied experimentally at 5.5–7.5 GPa and 1200–1450?°C in terms of phase relations and conditions for formation and stability of phlogopite. The starting samples are phlogopite–carbonatite–phlogopite sandwiches and harzburgite–carbonatite mixtures consisting of Ol?+?Grt?+?Cpx?+?L (±Opx), according to the previous experimental results obtained at the same P–T parameters but in water-free systems. Carbonatite is represented by a K- and Ca-rich composition that may form at the top of a slab. In the presence of carbonatitic melt, phlogopite can partly melt in a peritectic reaction at 5.5 GPa and 1200–1350?°C, as well as at 6.3–7.0 GPa and 1200?°C: 2Phl?+?CaCO₃ (L)?Cpx?+?Ol?+?Grt?+?K₂CO₃ (L)?+?2H₂O (L). Synthesis of phlogopite at 5.5 GPa and 1200–1350?°C, with an initial mixture of H₂O-bearing harzburgite and carbonatite, demonstrates experimentally that equilibrium in this reaction can be shifted from right to left. Therefore, phlogopite can equilibrate with ultrapotassic carbonate–silicate melts in a?≥?150?°C region between 1200 and 1350?°C at 5.5 GPa. On the other hand, it can exist but cannot nucleate spontaneously and crystallize in the presence of such melts in quite a large pressure range in experiments at 6.3–7.0 GPa and 1200?°C. Thus, phlogopite can result from metasomatism of peridotite at the base of continental lithospheric mantle (CLM) by ultrapotassic carbonatite agents at depths shallower than 180–195 km, which creates a mechanism of water retaining in CLM. Kimberlite formation can begin at 5.5 GPa and 1350?°C in a phlogopite-bearing peridotite source generating a hydrous carbonate–silicate melt with 10–15 wt% SiO₂, Ca# from 45 to 60, and high K enrichment. Upon further heating to 1450?°C due to the effect of a mantle plume at the CLM base, phlogopite disappears and a kimberlite-like melt forms with SiO₂ to 20 wt% and Ca#?=?35–40.     

Crystallization conditions of peraluminous charnockites: constraints from mineral thermometry and thermodynamic modelling

Most igneous charnockites are interpreted to have crystallized at hot and dry conditions, i.e. at >800?°C and <3 wt.% H₂O and with an important CO₂ component in the system. These charnockites are metaluminous to weakly peraluminous and their formation involves a significant mantle-derived component. This study, in contrast, investigates the crystallization conditions of strongly peraluminous, metasediment-sourced charnockites from the Qinzhou Bay Granitic Complex, South China. To constrain the temperature-melt H₂O crystallization paths for the studied peraluminous charnockites, petrographic characterization was combined with fluid inclusion compositional data, mineral thermometry, and thermodynamic modelling. The uncertainties of the thermodynamic modelling in reconstructing the crystallization conditions of the granitic magmas have been evaluated by comparison between modelled and experimental phase relations for a moderately evolved, peraluminous granite (~70 wt.% SiO₂). The comparison suggests that the modelling reproduces the experimentally derived phase saturation boundaries with uncertainties of 20–60?°C and 0.5–1 wt.% H₂O for systems with ≤1–2 wt.% initial melt H₂O at ~0.2 GPa. For the investigated natural systems, the thermometric estimates and modelling indicate that orthopyroxene crystallized at relatively low temperature (750–790?±?30?°C) and moderately high to high melt H₂O content (3.5–5.6?±?0.5 wt.%). The charnockites finally solidified at relatively “cold” and “wet” conditions. This suggests that thermodynamic modelling affords a possible approach to constrain charnockite crystallization as tested here for peraluminous, moderately low pressure (≤0.3 GPa), and overall H₂O-poor systems (≤1–2 wt.% H₂O total), but yields results with increasing uncertainty for high-pressure or H₂O-rich granitic systems.     

Discussion: “Xenoliths in ultrapotassic volcanic rocks in the Lhasa block: direct evidence for crust–mantle mixing and metamorphism in the deep crust” by Wang et al. 2016 (Contributions to Mineralogy and Petrology) 171:62

Wang et al. (Contrib Mineral Petrol 171:62, 2016a) present data on composition of xenolith from Southern Tibet and conclude that ulrapotassic melts from the region formed by melting mantle, and complex interaction with a crustal component. In this discussion we demonstrate that numerous observations presented by Wang et al. (2016a) can be explained by partial melting of crust followed by interaction between that melt and the mantle. We show that this model can explain the variability of magmas in such suits without evoking occurrence of coincidental, unrelated events. Moreover we demonstrate that our model of a crustal origin of the proto-shoshonite melts is now supported by independent lines of evidence such as geochemistry of restites after high- and ultrahigh- pressure melting and melt inclusion studies.     

A geochemical approach to constraining the formation of glassy fallout debris from nuclear tests

Glassy nuclear fallout debris from near-surface nuclear tests is fundamentally reprocessed earth material. A geochemical approach to analysis of glassy fallout is uniquely suited to determine the means of reprocessing and shed light on the mechanisms of fallout formation. An improved understanding of fallout formation is of interest both for its potential to guide post-detonation nuclear forensic investigations and in the context of possible affinities between glassy debris and other glasses generated by high-energy natural events, such as meteorite impacts and lightning strikes. This study presents a large major-element compositional dataset for glasses within aerodynamic fallout from the Trinity nuclear test (“trinitite”) and a geochemically based analysis of the glass compositional trends. Silica-rich and alkali-rich trinitite glasses show compositions and textures consistent with formation through melting of individual mineral grains—quartz and alkali feldspar, respectively—from the test-site sediment. The volumetrically dominant glass phase—called the CaMgFe glass—shows extreme major-element compositional variability. Compositional trends in the CaMgFe glass are most consistent with formation through volatility-controlled condensation from compositionally heterogeneous plasma. Radioactivity occurs only in CaMgFe glass, indicating that co-condensation of evaporated bulk ground material and trace device material was the main mechanism of radioisotope incorporation into trinitite. CaMgFe trinitite glasses overlap compositionally with basalts, rhyolites, fulgurites, tektites, and microtektites but display greater compositional diversity than all of these naturally formed glasses. Indeed, the most refractory CaMgFe glasses compositionally resemble early solar system condensates—specifically, CAIs.     

Closed Circuit Television (CCTV) Earthquake Behaviour Coding Methodology: analysis of Christchurch Public Hospital video data from the 22 February Christchurch earthquake event

Earthquake epidemiological research indicates that the behavioural response influences the nature and severity of injuries sustained. However, there is no observational evidence of the actions individuals engage in during and immediately following earthquake shaking, and the context in which earthquake injuries and deaths are caused. Closed Circuit Television Earthquake Behaviour Coding Methodology has been developed as a tool to classify human behaviour during and immediately following earthquake shaking using real event video data. The coding methodology was applied to security video data captured during the 22 February 2011 M_w6.3 ‘Christchurch’ earthquake event from the Christchurch Public Hospital which experienced shaking intensity of MM9 lasting approximately 12–15 s. We applied this coding methodology to determine: demography, trans-event behavioural responses, post-event behavioural responses, influence of social context on behaviour, and influence of behaviour on injuries. A total of 213 individuals from 31 different camera views were analysed. Sixty-six per cent of the individuals were adult-aged females. The primary trans-event responses were to hold (26%) onto furniture, walls, and/or other people close to them and to look around (30%). No individuals were observed to perform all ‘Drop, Cover, Hold’ actions, the recommended action during strong earthquake shaking in New Zealand. Post-event behaviour included: running, walking, providing assistance, moving towards others, visual communication, and some individuals gave instructions. Social contextual behaviour varied depending on the role of the adult. There were no serious injuries linked to behaviour. The results of this initial study indicate the coding methodology can record the distribution of and variation in human behaviours. Therefore, objective observation of earthquake video data can provide a useful quantitative measure of human behaviour. Significantly, the process will enable researchers to look more closely at behaviours, as well as the social and physical contexts associated with injury risk during and immediately following earthquake shaking.     

Magnetostratigraphic correlation of the Oxfordian–Kimmeridgian boundary

A magnetic polarity pattern for Boreal and Sub-Boreal ammonite zones of the Upper Oxfordian to Lower Kimmeridgian was established and confirmed in four British sections, including the proposed Global Boundary Stratotype Section and Point (GSSP) on the Isle of Skye (Scotland) to define the base of the international Kimmeridgian Stage. A coeval pattern for Sub-Mediterranean ammonite zones was compiled from seven sections in Poland, one German section and multi-section composites from France and Spain. The mean paleopole for the European Craton (excluding Spain) at the Oxfordian–Kimmeridgian boundary is 74.2°N, 181.3°E (Α₉₅ = 3.8°). The common magnetic polarity scale enables inter-correlation of ammonite subzones among these three faunal provinces and to the marine magnetic-anomaly M-Sequence. The proposed GSSP at the base of the Pictonia baylei Zone is near the base of an extended interval dominated by reversed polarity, which is interpreted to be Chron M26r. This GSSP level projects to the lower to middle part of the Epipeltoceras bimammatum Subzone, which is the middle subzone of this E. bimammatum Zone in the Sub-Mediterranean standard zonation. In contrast, the traditional placement of the Oxfordian–Kimmeridgian boundary in that Sub-Mediterranean standard zonation (base of Sutneria platynota Zone) is at the base of Chron M25r, or nearly 1 million years younger.