Enrichments of the mantle sources beneath the Southern Volcanic Zone (Andes) by fluids and melts derived from abraded upper continental crust

Mafic basaltic-andesitic volcanic rocks from the Andean Southern Volcanic Zone (SVZ) exhibit a northward increase in crustal components in primitive arc magmas from the Central through the Transitional and Northern SVZ segments. New elemental and Sr–Nd-high-precision Pb isotope data from the Quaternary arc volcanic centres of Maipo (NSVZ) and Infernillo and Laguna del Maule (TSVZ) are argued to reflect mainly their mantle source and its melting. For the C-T-NSVZ, we identify two types of source enrichment: one, represented by Antuco in CSVZ, but also present northward along the arc, was dominated by fluids which enriched a pre-metasomatic South Atlantic depleted MORB mantle type asthenosphere. The second enrichment was by melts having the characteristics of upper continental crust (UCC), distinctly different from Chile trench sediments. We suggest that granitic rocks entered the source mantle by means of subduction erosion in response to the northward increasingly strong coupling of the converging plates. Both types of enrichment had the same Pb isotope composition in the TSVZ with no significant component derived from the subducting oceanic crust. Pb–Sr–Nd isotopes indicate a major crustal compositional change at the southern end of the NSVZ. Modelling suggests addition of around 2 % UCC for Infernillo and 5 % for Maipo.     

Amphibole as an archivist of magmatic crystallization conditions: problems,potential, and implications for inferring magma storage prior to the paroxysmal 2010 eruption of Mount Merapi,Indonesia

Amphibole is widely employed to calculate crystallization temperature and pressure, although its potential as a geobarometer has always been debated. Recently, Ridolfi et al. (Contrib Mineral Petrol 160:45–66, 2010) and Ridolfi and Renzulli (Contrib Mineral Petrol 163:877–895, 2012) have presented calibrations for calculating temperature, pressure, fO₂, melt H₂O, and melt major and minor oxide composition from amphibole with a large compositional range. Using their calibrations, we have (i) calculated crystallization conditions for amphibole from eleven published experimental studies to examine the problems and the potential of the new calibrations; and (ii) calculated crystallization conditions for amphibole from basaltic–andesitic pyroclasts erupted during the paroxysmal 2010 eruption of Mount Merapi in Java, Indonesia, to infer pre-eruptive conditions. Our comparison of experimental and calculated values shows that calculated crystallization temperatures are reasonable estimates. Calculated fO₂ and melt SiO₂ content yields potentially useful estimates at moderately reduced to moderately oxidized conditions and intermediate to felsic melt compositions. However, calculated crystallization pressure and melt H₂O content are untenable estimates that largely reflect compositional variation in the crystallizing magmas and crystallization temperature and not the calculated parameters. Amphibole from Merapi’s pyroclasts yields calculated conditions of ~200–800 MPa, ~900–1,050 °C, ~NNO + 0.3–NNO + 1.1, ~3.7–7.2 wt% melt H₂O, and ~58–71 wt% melt SiO₂. We interpret the variations in calculated temperature, fO₂, and melt SiO₂ content as reasonable estimates, but conclude that the large calculated pressure variation for amphibole from Merapi and many other arc volcanoes is evidence for thorough mixing of mafic to felsic magmas and not necessarily evidence for crystallization over a large depth range. In contrast, bimodal pressure estimates obtained for other arc magmas reflect amphibole crystallization from mafic and more evolved magmas, respectively, and should not necessarily be taken as evidence for crystallization in two reservoirs at variable depth.     

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.     

A 2000-yr reconstruction of air temperature in the Great Basin of the United States with specific reference to the Medieval Climatic Anomaly

A sediment core representing the past two millennia was recovered from Stella Lake in the Snake Range of the central Great Basin in Nevada. The core was analyzed for sub-fossil chironomids and sediment organic content. A quantitative reconstruction of mean July air temperature (MJAT) was developed using a regional training set and a chironomid-based WA-PLS inference model (r²_jack = 0.55, RMSEP = 0.9°C). The chironomid-based MJAT reconstruction suggests that the interval between AD 900 and AD 1300, corresponding to the Medieval Climate Anomaly (MCA), was characterized by MJAT elevated 1.0°C above the subsequent Little Ice Age (LIA), but likely not as warm as recent conditions. Comparison of the Stella Lake temperature reconstruction to previously published paleoclimate records from this region indicates that the temperature fluctuations inferred to have occurred at Stella Lake between AD 900 and AD 1300 correspond to regional records documenting hydroclimate variability during the MCA interval. The Stella Lake record provides evidence that elevated summer temperature contributed to the increased aridity that characterized the western United States during the MCA.     

Warm‐temperate,marine, carbonate sedimentation in an Early Miocene,tide‐influenced,incised valley; Provence,south‐east France

The Saumane‐Venasque compound palaeovalley succession accumulated in a strongly tide‐influenced embayment or estuary. Warm‐temperate normal marine to brackish conditions led to deposition of extensive cross‐bedded biofragmental calcarenites. Echinoids, bryozoans, coralline algae, barnacles and benthic foraminifera were produced in seagrass meadows, on rocky substrates colonized by macroalgae and within subaqueous dune fields. There are two sequences, S1 and S2, the first of which contains three high‐frequency sequences (S1a, S1b and S1c). Sequence 1 is largely confined to the palaeovalley with its upper part covering interfluves. Each of these has a similar upward succession of deposits that includes: (i) a basal erosional surface that is bored and glauconitized; (ii) a discontinuous lagoonal lime mudstone or wackestone; (iii) a thin conglomerate generated by tidal ravinement; (iv) a transgressive systems tract series of cross‐bedded calcarenites; (v) a maximum flooding interval of argillaceous, muddy quartzose, open‐marine limestones; and (vi) a thin highstand systems tract of fine‐grained calcarenite. Tidal currents during stages S1a, S1b and S1c were accentuated by the constricted valley topography, whereas basin‐scale factors enhanced tidal currents during the deposition of S2. The upper part of the succession in all but S1c has been removed by later erosion. There is an overall upward temporal change with quartz, barnacles, encrusting corallines and epifaunal echinoids decreasing but bryozoans, articulated corallines and infaunal echinoids increasing. This trend is interpreted to be the result of changing oceanographic conditions as the valley was filled, bathymetric relief was reduced, rocky substrates were replaced as carbonate factories by seagrass meadows and subaqueous dunes, and the setting became progressively less confined and more open marine. These limestones are characteristic of a suite of similar cool‐water calcareous sand bodies in environments with little siliciclastic or fresh water input during times of high‐amplitude sea‐level change wherein complex inboard antecedent topography was flooded by a rising ocean.     

Scales and causes of heterogeneity in bars in a large multi‐channel river: Río Paraná, Argentina

To date, published studies of alluvial bar architecture in large rivers have been restricted mostly to case studies of individual bars and single locations. Relatively little is known about how the depositional processes and sedimentary architecture of kilometre‐scale bars vary within a multi‐kilometre reach or over several hundreds of kilometres downstream. This study presents Ground Penetrating Radar and core data from 11, kilometre‐scale bars from the Río Paraná, Argentina. The investigated bars are located between 30 km upstream and 540 km downstream of the Río Paraná – Río Paraguay confluence, where a significant volume of fine‐grained suspended sediment is introduced into the network. Bar‐scale cross‐stratified sets, with lengths and widths up to 600 m and thicknesses up to 12 m, enable the distinction of large river deposits from stacked deposits of smaller rivers, but are only present in half the surface area of the bars. Up to 90% of bar‐scale sets are found on top of finer‐grained ripple‐laminated bar‐trough deposits. Bar‐scale sets make up as much as 58% of the volume of the deposits in small, incipient mid‐channel bars, but this proportion decreases significantly with increasing age and size of the bars. Contrary to what might be expected, a significant proportion of the sedimentary structures found in the Río Paraná is similar in scale to those found in much smaller rivers. In other words, large river deposits are not always characterized by big structures that allow a simple interpretation of river scale. However, the large scale of the depositional units in big rivers causes small‐scale structures, such as ripple sets, to be grouped into thicker cosets, which indicate river scale even when no obvious large‐scale sets are present. The results also show that the composition of bars differs between the studied reaches upstream and downstream of the confluence with the Río Paraguay. Relative to other controls on downstream fining, the tributary input of fine‐grained suspended material from the Río Paraguay causes a marked change in the composition of the bar deposits. Compared to the upstream reaches, the sedimentary architecture of the downstream reaches in the top ca 5 m of mid‐channel bars shows: (i) an increase in the abundance and thickness (up to metre‐scale) of laterally extensive (hundreds of metres) fine‐grained layers; (ii) an increase in the percentage of deposits comprised of ripple sets (to >40% in the upper bar deposits); and (iii) an increase in bar‐trough deposits and a corresponding decrease in bar‐scale cross‐strata (<10%). The thalweg deposits of the Río Paraná are composed of dune sets, even directly downstream from the Río Paraguay where the upper channel deposits are dominantly fine‐grained. Thus, the change in sedimentary facies due to a tributary point‐source of fine‐grained sediment is primarily expressed in the composition of the upper bar deposits.