Temporal evolution and compositional signatures of two supervolcanic systems recorded in zircons from Mangakino volcanic centre,New Zealand

Mangakino, the oldest rhyolitic caldera centre delineated in the Taupo Volcanic Zone of New Zealand, generated two very large (super-sized) ignimbrite eruptions, the 1.21 ± 0.04 Ma >500 km³ Ongatiti and ~1.0 Ma ~1,200 km³ Kidnappers events, the latter of which was followed after a short period of erosion by the ~200 km³ Rocky Hill eruption. We present U/Pb ages and trace-element analyses on zircons from pumice clasts from these three eruptions by Secondary Ion Mass Spectrometry (SIMS) using SHRIMP-RG instruments to illustrate the evolution of the respective magmatic systems. U–Pb age spectra from the Ongatiti imply growth of the magmatic system over ~250 kyr, with a peak of crystallisation around 1.32 Ma, ~100 kyr prior to eruption. The zircons are inferred to have then remained stable in a mush with little crystallisation and/or dissolution before later rejuvenation of the system at the lead-in to eruption. The paired Kidnappers and Rocky Hill eruptions have U–Pb zircon ages and geochemical signatures that suggest they were products of a common system grown over ~200 kyr. The Kidnappers and Rocky Hill samples show similar weakly bimodal age spectra, with peaks at 1.1 and 1.0 Ma, suggesting that an inherited antecrystic population was augmented by crystals grown at ages within uncertainty of the eruption age. In the Kidnappers, this younger age peak is dominantly seen in needle-shaped low U grains with aspect ratios of up to 18. In all three deposits, zircon cores show larger ranges and higher absolute concentrations of trace elements than zircon rims, consistent with zircon crystallisation from evolving melts undergoing crystal fractionation involving plagioclase and amphibole. Abundances and ratios of many trace elements frequently show variations between different sectors within single grains, even where there is no visible sector zoning in cathodoluminescence (CL) imaging. Substitution mechanisms, as reflected in the molar (Sc + Y + REE³⁺)/P ratio, differ in the same growth zone between the sides (along a-axis and b-axis: values approaching 1.0) and tips (c-axis: values between 1.5 and 5.0) of single crystals. These observations have implications for the use of zircons for tracking magmatic processes, particularly in techniques where CL zonation within crystals is not assessed and small analytical spot sizes cannot be achieved. These observations also limit applicability of the widely used Ti-in-zircon thermometer. The age spectra for the Ongatiti and Kidnappers/Rocky Hill samples indicate that both magmatic systems were newly built in the time-breaks after respective previous large eruptions from Mangakino. Trace element variations defining three-component mixing suggest that zircons, sourced from multiple melts, contributed to the population in each system.     

Building zoned ignimbrites by recycling silicic cumulates: insight from the 1,000 km<Superscript>3</Superscript> Carpenter Ridge Tuff,CO

The ~1,000 km³ Carpenter Ridge Tuff (CRT), erupted at 27.55 Ma during the mid-tertiary ignimbrite flare-up in the western USA, is among the largest known strongly zoned ash-flow tuffs. It consists primarily of densely welded crystal-poor rhyolite with a pronounced, highly evolved chemical signature (high Rb/Sr, low Ba, Zr, Eu), but thickly ponded intracaldera CRT is capped by a more crystal-rich, less silicic facies. In the outflow ignimbrite, this upper zone is defined mainly by densely welded crystal-rich juvenile clasts of trachydacite composition, with higher Fe–Ti oxide temperatures, and is characterized by extremely high Ba (to 7,500 ppm), Zr, Sr, and positive Eu anomalies. Rare mafic clasts (51–53 wt% SiO₂) with Ba contents to 4,000–5,000 ppm and positive Eu anomalies are also present. Much of the major and trace-element variations in the CRT juvenile clasts can be reproduced via in situ differentiation by interstitial melt extraction from a crystal-rich, upper-crustal mush zone, with the trachydacite, crystal-rich clasts representing the remobilized crystal cumulate left behind by the melt extraction process. Late recharge events, represented by the rare mafic clasts and high-Al amphiboles in some samples, mixed in with parts of the crystal cumulate and generated additional scatter in the whole-rock data. Recharge was important in thermally remobilizing the silicic crystal cumulate by partially melting the near-solidus phases, as supported by: (1) ubiquitous wormy/sieve textures and reverse zoning patterns in feldspars and biotites, (2) absence of quartz in this very silicic unit stored at depths of >4–5 km, and (3) heterogeneous melt compositions in the trachydacite fiamme and mafic clasts, particularly in Ba, indicating local enrichment of this element due mostly to sanidine and biotite melting. The injection of hot, juvenile magma into the upper-crustal cumulate also imparted the observed thermal gradient to the deposits and the mixing overprint that partly masks the in situ differentiation process. The CRT provides a particularly clear perspective on processes of in situ crystal-liquid separation into a lower crystal-rich zone and an upper eruptible cap, which appears common in incrementally built upper-crustal magma reservoirs of high-flux magmatic provinces.     

Vanadium and niobium behavior in rutile as a function of oxygen fugacity: evidence from natural samples

Vanadium occurs in multiple valence states in nature, whereas Nb is exclusively pentavalent. Both are compatible in rutile, but the relationship of V–Nb partitioning and dependence on oxygen fugacity (expressed as fO₂) has not yet been systematically investigated. We acquired trace-element concentrations on rutile grains (n = 86) in nine eclogitic samples from the Dabie-Sulu orogenic belt by laser ablation inductively coupled plasma mass spectrometry (LA–ICP–MS) and combined them with published results in order to assess the direct and indirect effects of oxygen fugacity on the partitioning of V and Nb into rutile. A well-defined negative correlation between Nb (7–1,200 ppm) and V concentrations (50–3,200 ppm) was found, documenting a competitive relationship in the rutile crystal that does not appear to be controlled by bulk rock or mineral compositions. Based on the published relationship of ^RtD_V and V valence with ?QFM, we suggest that the priority order of V incorporation into rutile is V⁴⁺ > V³⁺ > V⁵⁺. The inferred Nb–V competitive relationship in rutile from the Dabie-Sulu orogenic belt could be explained by decreasing fO₂ due to dehydration reactions involving loss of oxidizing fluids during continental subduction: The increased proportion of V³⁺ (expressed as V³⁺/∑V) and attendant decrease in ^RtD_V is suggested to lead to an increase in rutile lattice sites available for Nb⁵⁺. A similar effect may be observed under more oxidizing conditions. When V⁵⁺/∑V increases, ^RtD_V shows a dramatic decline and Nb concentration increases considerably. This is possibly documented by rutile in highly metasomatized and oxidized MARID-type (MARID: mica–amphibole–rutile–ilmenite–diopside) mantle xenoliths from the Kaapvaal craton, which also show a negative V–Nb covariation. In addition, their Nb/Ta covaries with V concentrations: For V concentrations <1,250 ppm, Nb/Ta ranges between 35 and 45, whereas for V > 1,250 ppm, Nb/Ta is considerably lower (5–15). This relationship is mainly controlled by a change in Nb concentrations, suggesting that the indirect dependence of ^RtD_Nb on fO₂, which is not mirrored in ^RtD_Ta, can exert considerable influence on rutile Nb–Ta fractionation.     

Jurassic plume-origin ophiolites in Japan: accreted fragments of oceanic plateaus

The Mikabu and Sorachi–Yezo belts comprise Jurassic ophiolitic complexes in Japan, where abundant basaltic to picritic rocks occur as lavas and hyaloclastite blocks. In the studied northern Hamamatsu and Dodaira areas of the Mikabu belt, these rocks are divided into two geochemical types, namely depleted (D-) and enriched (E-) types. In addition, highly enriched (HE-) type has been reported from other areas in literature. The D-type picrites contain highly magnesian relic olivine phenocrysts up to Fo_93.5, and their Fo–NiO trend indicates fractional crystallization from a high-MgO primary magma. The MgO content is calculated as high as 25 wt%, indicating mantle melting at unusually high potential temperature (T _p) up to 1,650 °C. The E-type rocks represent the enrichment in Fe and LREE and the depletion in Mg, Al and HREE relative to the D-type rocks. These chemical characteristics are in good accordance with those of melts from garnet pyroxenite melting. Volcanics in the Sorachi–Yezo belts can be divided into the same types as the Mikabu belt, and the D-type picrites with magnesian olivines also show lines of evidence for production from high T _p mantle. Evidence for the high T _p mantle and geochemical similarities with high-Mg picrites and komatiites from oceanic and continental large igneous provinces (LIPs) indicate that the Mikabu and Sorachi–Yezo belts are accreted oceanic LIPs that were formed from hot large mantle plumes in the Late Jurassic Pacific Ocean. The E- and D-type rocks were formed as magmas generated by garnet pyroxenite melting at an early stage of LIP magmatism and by depleted peridotite melting at the later stage, respectively. The Mikabu belt characteristically bears abundant ultramafic cumulates, which could have been formed by crystal accumulation from a primary magma generated from Fe-rich peridotite mantle source, and the HE-type magma were produced by low degrees partial melting of garnet pyroxenite source. They should have been formed later and in lower temperatures than the E- and D-type rocks. The Mikabu and Sorachi Plateaus were formed in a low-latitude region of the Late Jurassic Pacific Ocean possibly near a subduction zone, partially experienced high P/T metamorphism during subduction, and then uplifted in association with (or without, in case of Mikabu) the supra-subduction zone ophiolite. The Mikabu and Sorachi Plateaus may be the Late Jurassic oceanic LIPs that could have been formed in brotherhood with the Shatsky Rise.     

Emplacement ages and sources of kimberlites and related rocks in southern Africa: U–Pb ages and Sr–Nd isotopes of groundmass perovskite

Groundmass perovskite has been dated by LA-ICPMS in 135 kimberlites and related rocks from 110 localities across southern Africa. Sr and/or Nd isotopes have been analysed by LA-MC-ICPMS in a subset of these and integrated with published data. The age distribution shows peaks at 1,600–1,800, 1,000–1,200, 500–800 and 50–130 Ma. The major “bloom” of Group I kimberlites at ca 90 ± 10 Ma was preceded by a slow build-up in magmatic activity from ca 180 Ma. The main pulse of Group II kimberlites at 120–130 Ma was a distinct episode within this build-up. Comparison of the isotopic data with seismic tomography images suggests that metasomatized subcontinental lithospheric mantle (SCLM) with very low ε _Nd and high ⁸⁷Sr/⁸⁶Sr, (the isotopic signature of Group II kimberlites) was focused in low-Vs zones along translithospheric structures. Such metasomatized zones existed as early as 1,800 Ma, but were only sporadically tapped until the magmatic build-up began at ca 180 Ma, and contributed little to the kimberlitic magmas after ca 110 Ma. We suggest that these metasomatized volumes resided in the deep SCLM and that their low-melting point components were “burned off” by rising temperatures, presumably during an asthenospheric upwelling that led to SCLM thinning and a rise in the ambient geotherm between 120 and 90 Ma. The younger Group I kimberlites therefore rarely interacted with such SCLM, but had improved access to shallower volumes of differently metasomatized, ancient SCLM with low ⁸⁷Sr/⁸⁶Sr and intermediate ε _Nd (0–5). The kimberlite compositions therefore reflect the evolution of the SCLM of southern Africa, with metasomatic-enrichment events from as early as 1.8 Ga, through a major thermal and compositional change at ca 110 Ma, and the major kimberlite “bloom” around 90 Ma.     

The tripole: A new coherent vortex structure of incompressible two-dimensional flows

Abstract

Using a contour dynamical algorithm, we have found rotating tripolar V-state solutions for the inviscid Euler equations in two-dimensions. We have studied their geometry as a function of their physical parameters. Their stability was investigated with the aid of contour surgery, and most of the states were found to be stable. Under finite-amplitude perturbations, tripoles are shown to either fission into two asymmetric dipoles or to evolve into a shielded axisymmetric vortex, demonstrating the existence of two new ‘‘reversible transitions'’ between topologically distinct coherent vortex structures. These dynamical results are confirmed by pseudo-spectral simulations, with which we also show how continuous tripolar long-lived coherent vortex structures can be generated in a variety of ways.     

Numerical Studies of Fluid Flow in Microgravity Conditions for Confined Crystal Growth

We study the convective flow induced by residual accelerations in microgravity conditions for different geometric arrangements which are relevant to crystal growth experiments. We consider both constantand oscillating acceleration and focus mostly on the transient relaxation dynamics. Results are relevant to estimate impact of more realistic residual accelerations in crystal growth experiments.     

Continual long-term monitoring of methane in wells above the Utica Shale using total dissolved gas pressure probes

Hydrogeology Journal - Monitoring of dissolved methane concentrations in groundwater is required to identify impacts from oil and gas development and to understand temporal variability under...     

New Validity Conditions for the Multivariate Matérn Coregionalization Model,with an Application to Exploration Geochemistry

Mathematical Geosciences - This paper addresses the problem of finding parametric constraints that ensure the validity of the multivariate Matérn covariance for modeling the spatial...     

Construction of Structural Geological Model Using Monte Carlo Simulation

Geotechnical and Geological Engineering - To optimize the prediction of structural geological conditions in the underground as of data collected at the surface, due to the usual great uncertainties...