Arc Magmas from Slab to Eruption: The Case of Kliuchevskoy Volcano

Arc magmas are generated by a number of mantle and crustal processes. Our multidisciplinary, long-term research is aimed at deciphering these processes for a single arc volcano, Kliuchevskoy volcano in Kamchatka. Some key results of the study follow: 1) Modeling of trace element and H2O contents in melt inclusions suggests that the primary magmas originate via hydrous flux-melting of the mantle wedge at temperatures close to the dry peridotite solidus. The role of decompression melting is minor or absent at Kliuchevskoy and other arc volcanoes built on relatively thick crust. 2) Geochemistry of high-Mg olivine suggests that primary Kliuchevskoy magmas have substantial contribution from olivine-free pyroxenite (up to 30 %), which could be formed by reaction of slab melts (or supercritical fluids) with mantle wedge peridotite. 3) Parental Kliuchevskoy melts start to crystallize as deep as the Moho boundary, and the erupted magmas reflect multistage and complex processes of crystallization, magma mixing and crustal assimilation. None of the Kliuchevskoy rocks analyzed thus far represent true primary melt compositions. 4) The Kliuchevskoy Holocene eruptive history is not steady-state in terms of eruption rate and geochemistry. There are two millenial cycles with major and trace element and OSr- Nd-Pb and U-series isotope compositions of the magmas changing gradually from more to less affected by crustal (?) assimilation. The onset of the cycles correlates with periods of enhanced volcanic activity in Kamchatka, suggesting that the extent of magma-crust interaction is inversely related to magma production rate and thus magma flux from the mantle

Rare sulfur and triple oxygen isotope geochemistry of volcanogenic sulfate aerosols

We present analyses of stable isotopic ratios ^(17)O/^(16)O, ^(18)O/^(16)O, ^(34)S/^(32)S, and ^(33)S/^(32)S, ^(36)S/^(32)S in sulfate leached from volcanic ash of a series of well known, large and small volcanic eruptions. We consider eruptions of Mt. St. Helens (Washington, 1980, ∼1 km^3), Mt. Spurr (Alaska, 1953, <1 km3), Gjalp (Iceland, 1996, 1998, <1 km^3), Pinatubo (Phillipines, 1991, 10 km^3), Bishop tuff (Long Valley, California, 0.76 Ma, 750 km^3), Lower Bandelier tuff (Toledo Caldera, New Mexico, 1.61 Ma, 600 km^3), and Lava Creek and Huckleberry Ridge tuffs (Yellowstone, Wyoming, 0.64 Ma, 1000 km^3 and 2.04 Ma 2500 km^3, respectively). This list covers much of the diversity of sizes and the character of silicic volcanic eruptions. Particular emphasis is paid to the Lava Creek tuff for which we present wide geographic sample coverage. This global dataset spans a significant range in δ^(34)S, δ^(18)O, and Δ^(17)O of sulfate (29‰, 30‰, and 3.3‰, respectively) with oxygen isotopes recording mass-independent (Δ^(17)O > 0.2‰) and sulfur isotopes exhibiting mass-dependent behavior. Products of large eruptions account for most of‘ these isotopic ranges. Sulfate with Δ^(17)O > 0.2‰ is present as 1–10 μm gypsum crystals on distal ash particles and records the isotopic signature of stratospheric photochemical reactions. Sediments that embed ash layers do not contain sulfate or contain little sulfate with Δ^(17)O near 0‰, suggesting that the observed sulfate in ash is of volcanic origin. Mass-dependent fractionation of sulfur isotopic ratios suggests that sulfate-forming reactions did not involve photolysis of SO2, like that inferred for pre-2.3 Ga sulfates from Archean sediments or Antarctic ice-core sulfate associated with few dated eruptions. Even though the sulfate sulfur isotopic compositions reflect mass-dependent processes, the products of caldera-forming eruptions display a large δ^(34)S range and exhibit fractionation relationships that do not follow the expected equilibrium slopes of 0.515 and 1.90 for ^(33)S/^(32)S vs. ^(34)S/^(32)S and ^(36)S/^(32)S vs. ^(34)S/^(32)S, respectively. The data presented here are consistent with modification of a chemical mass-dependent fractionation of sulfur isotopes in the volcanic plume by either a kinetic gas phase reaction of volcanic SO_2 with OH and/or a Rayleigh processes involving a residual Rayleigh reactant—volcanic SO_2 gas, rather than a Rayleigh product. These results may also imply at least two removal pathways for SO_2 in volcanic plumes. Above-zero Δ^(17)O values and their positive correlation with δ^(18)O in sulfate can be explained by oxidation by high-δ^(18)O and high-Δ^(17)O compounds such as ozone and radicals such as OH that result from ozone break down. Large caldera-forming eruptions have the highest Δ^(17)O values, and the largest range of δ^(18)O, which can be explained by stratospheric reaction with ozone-derived OH radicals. These results suggest that massive eruptions are capable of causing a temporary depletion of the ozone layer. Such depletion may be many times that of the measured 3–8% depletion following 1991 Pinatubo eruption, if the amount of sulfur dioxide released scales with the amount of ozone depletion

In situ dating and investigation of remarkably depleted –27.3‰ SMOW “Slushball” Earth zircons

Paleoproterozoic amphibolites and gneisses - that are remarkably depleted in ^(18)O are found in the Belomorian Belt in Karelia, Russia [1,2]. We mapped their extent to exceed 200x20km and affect metamorphosed mafic intrusions (est. ~2.4 Ga intrusion age) and host 2.6Ga gneisses found in this 1.9 Ga collisional belt. δ^(18)O values of –7 to –27.3‰ characterize minerals and rocks from several of these localities; some of these rocks are also remarkably depleted with respect to δD (-212 to –235‰ amphiboles). All have typical terrestrial Δ^(17)O values of 0‰. Based on previous paleogeographic reconstructions, we attribute the origin of these exotic O and H isotope compositions to the hydrothermal alteration associated with subglacial rifting during the Paleoproterozoic panglobal ice ages, but discuss additional possibilities: extremely low-δ^(18)O Paleo- proterozoic sea water, and excursion of Karelia to polar latitudes. Given that at high-T hydrothermal exchange equilibrium Δ^(18)O(rockwater) is close to zero, but water-rock interaction is rarely 100% efficient, the lowest measured δ^(18)O value in silicates likely gives the upper δ^(18)O bound for the altering meteoric fluid; we thus continues our quest to find the lowest δ^(18)O material such as a mineral assemblage or a tiny zircon fragment that would provide record of δ^(18)Owater

Multi-Cyclic and Isotopically Diverse Silicic Magma Generation in an Arc Volcano: Gorely Eruptive Center, Kamchatka, Russia

The Kamchatka Peninsula is home to some of the most frequent and prolific subduction-related volcanic activity in the world, with the largest number of caldera-forming eruptions per length of the volcanic arc. Among those, Gorely volcano has a topographically prominent Late Pleistocene caldera (13 km × 12 km, estimated to have produced >100 km3 of magma), which is now almost completely filled by a central cone. We report new 40Ar/39Ar ages and geochemical and isotopic data for newly recognized Mid-Pleistocene ignimbrite units of large but unknown volume sourced from the Gorely eruptive center, most of which were deposited in marginal glacial conditions. These ignimbrites have crystallinities of 9-24% and most are quartz-, amphibole-, and zircon-undersaturated. Additionally, we studied 32 eruptive units, including stratigraphically constrained Holocene tephra, and pre- and post-caldera lava sequences, to understand the petrogenetic and temporal evolution of this long-lived, multi-cyclic, arc volcano. Material erupted prior to the formation of the modern Gorely edifice, including the voluminous ignimbrites and eruptions of the ‘pra-Gorely' stage, consist primarily of dacite and andesite, whereas sequences of the modern Gorely edifice are represented by basalt to basaltic andesite. MELTS and EC-AFC modeling shows that it is possible to obtain silicic compositions near those of the evolved ignimbrite compositions through 60-75% fractional crystallization at 1 kbar and nickel-nickel oxide (NNO) oxygen fugacity. However, our newly compiled major and trace element data for Gorely yield two separate bimodal peaks in a SiO2-frequency diagram, showing a prominent Daly Gap, with a deficiency in andesite. Trace element concentrations define two separate trends, one for more silicic and another for more mafic sequences. Additionally, δ18Omelt values reconstructed from coexisting plagioclase and clinopyroxene phenocrysts range from a low value of 4·85‰ to a normal value of 6·22‰. The low δ18O values range throughout the known lifespan of Gorely, with the lowest value being from the first known ignimbrite to erupt, indicating episodic but temporally decreasing crustal assimilation of previously hydrothermally altered material. 87Sr/86Sr and 143Nd/144Nd ratios show wide ranges from 0·70328 to 0·70351 and from 0·51303 to 0·51309 respectively, also suggesting incorporation of surrounding crust, although there are less clear trends throughout the lifespan of Gorely. The combination of light and diverse δ18O values with elevated 87Sr/86Sr and low 143Nd/144Nd ratios suggests contamination by older and isotopically diverse, low-δ18O country-rocks, such as the neighboring 11 Ma Akhomten granitic massif, which shows ranges in δ18O, 87Sr/86Sr, and 144Nd/143Nd values overlapping with the Gorely magmas. In addition, the presence of glomerocrysts and mafic enclaves in the majority of Gorely dacites indicates a period of crystal settling and subsequent intrusion of hot, primitive basalt that probably triggered eruption. Finally, elevated Nb concentrations relative to other Kamchatkan volcanoes suggest that Gorely magmas may involve an enriched component, probably caused by delamination of a lower crustal root. Our results argue for an incremental view of silicic magma generation at so-called ‘long-term eruptive centers', in Kamchatka and worldwide, consisting of alternating episodes of magmatic and hydrothermal activity, and glacial advances and retreats. We demonstrate that large-volume, isotopically distinct, silicic magma can be generated rapidly between cone-building phases of volcanic activity through a combination of fractional crystallization, assimilation of older country rocks, and shallow assimilation of hydrothermally altered but otherwise petrochemically similar older intracaldera tuffs and intrusions. These transient shallow silicic magma chambers empty nearly completely in ignimbrite-forming eruptions after 103-105 years of assembly, partially triggered by glacial surface dynamic

Changing Mantle Sources and the Effects of Crustal Passage on the Steens Basalt, SE Oregon: Chemical and Isotopic Constraints

Continental flood basalts are more prone to compositional modification from passage through thicker and (or) more felsic crust in comparison to their oceanic counterparts. The Steens Basalt in southeast Oregon (~17 Ma) is among the oldest and most mafic members of the Columbia River Basalt Group and provides a record of the early stages of flood basalt volcanism. We evaluate the balance of mantle sources in time during the onset of Columbia River Basalt Group magmatism and assess the effect of crustal passage using stratigraphically controlled Sr, Nd, Pb, Hf, Os, and O isotopic compositions, as well as whole rock major and trace element data. Mixing models indicate that depleted and enriched mantle sources identified by previous workers contribute in varying proportions during the life of the magmatic system, with the greatest contribution by depleted mantle when eruption rate and presumed intrusion rate increase. During waxing, enrichment of δ18O in some flows signals cryptic deep fractionation of abundant clinopyroxene followed by shallow fractionation of olivine ± clinopyroxene ± plagioclase. Os concentrations are among the highest worldwide at a given MgO (0.29–0.86 ppb at 6.0 to 10.9 wt.%). We argue that high Os results from scavenging of sulfides by recharging magmas passing through earlier crystallized magmas. Elevated 87Sr/86Sr in the latest stage supports modest assimilation of partial melts from mafic accreted terranes, facilitated by thermal priming of crust by persistent magmatism. This work provides a more detailed schematic view of the Steens Basalt magmatic system, from mantle origin through crustal staging

Alteration of volcaniclastic deposits at Minna Bluff : geochemical insights on mineralizing environment and climate during the Late Miocene in Antarctica

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry, Geophysics, Geosystems 15 (2014): 3258–3280, doi:10.1002/2014GC005422.Secondary minerals in volcaniclastic deposits at Minna Bluff, a 45 km long peninsula in the Ross Sea, are used to infer processes of alteration and environmental conditions in the Late Miocene. Glassy volcaniclastic deposits are altered and contain phillipsite and chabazite, low to high-Mg carbonates, chalcedony, and clay. The δ18O of carbonates and chalcedony is variable, ranging from −0.50 to 21.53‰ and 0.68 to 10.37‰, respectively, and δD for chalcedony is light (−187.8 to −220.6‰), corresponding to Antarctic meteoric water. A mean carbonate 87Sr/86Sr ratio of 0.70327 ± 0.0009 (1σ, n = 12) is comparable to lava and suggests freshwater, as opposed to seawater, caused the alteration. Minerals were precipitated at elevated temperatures (91 and 104°C) based on quartz-calcite equilibrium, carbonate 13C-18C thermometry (Δ47 derived temperature = 5° to 43°C) and stability of zeolites in geothermal systems (>10 to ∼100°C). The alteration was a result of isolated, ephemeral events involving the exchange between heated meteoric water and glass during or soon after the formation of each deposit. Near-surface evaporative distillation can explain 18O-enriched compositions for some Mg-rich carbonates and chalcedony. The δ18Owater calculated for carbonates (−15.8 to −22.9‰) reveals a broad change, becoming heavier between ∼12 and ∼7 Ma, consistent with a warming climate. These findings are independently corroborated by the interpretation of Late Miocene sedimentary sequences recovered from nearby sediment cores. However, in contrast to a cold-based thermal regime proposed for ice flow at core sites, wet-based conditions prevailed at Minna Bluff; a likely consequence of high heat flow associated with an active magma system.This research was funded by a collaborative grant NSF OPP 05-38033. It also was supported by UNED/NSF 250550001146, NSF grants EAR-0949191, ARC-1215551, EAR-1325054, EAR-1352212, EAR-1049351, ACS grant 51182-DNI2, DOE grants DE-FG02-13ER16402, and DE-SC0010288, a Hellman Fellowship, and a Katzner grant (BGSU).2015-02-1

Across-arc geochemical variations in the Southern Volcanic Zone, Chile (34.5- 38.0°S): Constraints on Mantle Wedge and Input Compositions

Crustal assimilation (e.g. Hildreth and Moorbath, 1988) and/or subduction erosion (e.g. Stern, 1991; Kay et al., 2005) are believed to control the geochemical variations along the northern portion of the Chilean Southern Volcanic Zone. In order to evaluate these hypotheses, we present a comprehensive geochemical data set (major and trace elements and O-Sr-Nd-Hf-Pb isotopes) from Holocene primarily olivine-bearing volcanic rocks across the arc between 34.5-38.0°S, including volcanic front centers from Tinguiririca to Callaqui, the rear arc centers of Infernillo Volcanic Field, Laguna del Maule and Copahue, and extending 300 km into the backarc. We also present an equivalent data set for Chile Trench sediments outboard of this profile. The volcanic arc (including volcanic front and rear arc) samples primarily range from basalt to andesite/trachyandesite, whereas the backarc rocks are low-silica alkali basalts and trachybasalts. All samples show some characteristic subduction zone trace element enrichments and depletions, but the backarc samples show the least. Backarc basalts have higher Ce/Pb, Nb/U, Nb/Zr, and Ta/Hf, and lower Ba/Nb and Ba/La, consistent with less of a slab-derived component in the backarc and, consequently, lower degrees of mantle melting. The mantle-like δ18O in olivine and plagioclase phenocrysts (volcanic arc = 4.9-5.6 and backarc = 5.0-5.4 per mil) and lack of correlation between δ18O and indices of differentiation and other isotope ratios, argue against significant crustal assimilation. Volcanic arc and backarc samples almost completely overlap in Sr and Nd isotopic composition. High precision (double-spike) Pb isotope ratios are tightly correlated, precluding significant assimilation of older sialic crust but indicating mixing between a South Atlantic Mid Ocean-Ridge Basalt (MORB) source and a slab component derived from subducted sediments and altered oceanic crust. Hf-Nd isotope ratios define separate linear arrays for the volcanic arc and backarc, neither of which trend toward subducting sediment, possibly reflecting a primarily asthenospheric mantle array for the volcanic arc and involvement of enriched Proterozoic lithospheric mantle in the backarc. We propose a quantitative mixing model between a mixed-source, slab-derived melt and a heterogeneous mantle beneath the volcanic arc. The model is consistent with local geodynamic parameters, assuming water-saturated conditions within the slab

Large sulfur isotope fractionations in Martian sediments at Gale crater

Variability in the sulfur isotopic composition in sediments can reflect atmospheric, geologic and biological processes. Evidence for ancient fluvio-lacustrine environments at Gale crater on Mars and a lack of efficient crustal recycling mechanisms on the planet suggests a surface environment that was once warm enough to allow the presence of liquid water, at least for discrete periods of time, and implies a greenhouse effect that may have been influenced by sulfur-bearing volcanic gases. Here we report in situ analyses of the sulfur isotopic compositions of SO2 volatilized from ten sediment samples acquired by NASA’s Curiosity rover along a 13 km traverse of Gale crater. We find large variations in sulfur isotopic composition that exceed those measured for Martian meteorites and show both depletion and enrichment in 34S. Measured values of δ34S range from −47 ± 14‰ to 28 ± 7‰, similar to the range typical of terrestrial environments. Although limited geochronological constraints on the stratigraphy traversed by Curiosity are available, we propose that the observed sulfur isotopic signatures at Gale crater can be explained by equilibrium fractionation between sulfate and sulfide in an impact-driven hydrothermal system and atmospheric processing of sulfur-bearing gases during transient warm periods

Geochemical variations in the Central Southern Volcanic Zone, Chile (38-43°S): The role of fluids in generating arc magmas

We present new Sr-Nd-Pb-Hf-O isotope data from the volcanic arc (VA, volcanic front and rear arc) in Chile and the backarc (BA) in Argentina of the Central Southern Volcanic Zone in Chile (CSVZ; 38-43°S). Compared to the Transitional (T) SVZ (34.5-38°S; Jacques et al., 2013), the CSVZ VA has erupted greater volumes over shorter time intervals (Völker et al., 2011) and produced more tholeiitic melts. Although the CSVZ VA monogenetic cones are similar to the TSVZ VA samples, the CSVZ VA stratovolcanoes have higher ratios of highly fluid-mobile to less fluid-mobile trace elements (e.g. U/Th, Pb/Ce, Ba/Nb) and lower more- to less-incompatible fluid-immobile element ratios (e.g. La/Yb, La/Sm, Th/Yb, Nb/Yb), consistent with an overall higher fluid flux and greater degree of flux melting beneath the CSVZ stratovolcanoes compared to the CSVZ monogenetic centers and the TSVZ VA. The CSVZ monogenetic centers overlap the TSVZ in Sr and Nd isotopes, but the stratovolcanoes are shifted to higher Sr and/or Nd isotope ratios. The Pb isotopic composition of the CSVZ overlaps the TSVZ, which is clearly dominated by the composition of the trench sediments, but the CSVZ monogenetic samples extend to less radiogenic Pb isotope ratios. δ18Omelt from the CSVZ stratovolcano samples are below the MORB range, whereas the CSVZ monogenetic and the TSVZ samples fall within and slightly above the MORB range. The Nd and Hf isotopic ratios of the CSVZ VA extend to more radiogenic compositions than found in the TSVZ VA, indicating a greater contribution from a more depleted source. These correlations are interpreted to reflect derivation of fluids from hydrothermally altered oceanic crust and/or serpentinized upper mantle of the subducting plate. CSVZ BA basalts largely overlap TSVZ BA basalts, displaying less or no subduction influence compared to the VA, but some CSVZ BA basalts tap more enriched mantle, possibly subcontinental lithosphere, with distinctively lower Nd and Hf and elevated 207Pb/204Pb and 208Pb/204Pb isotope ratios