Jorge A. Vazquez

Months between rejuvenation and volcanic eruption at Yellowstone caldera, Wyoming

Rejuvenation of previously intruded silicic magma is an important process leading to effusive rhyolite, which is the most common product of volcanism at calderas with protracted histories of eruption and unrest such as Yellowstone caldera (Wyoming), Long Valley caldera (California), and Valles caldera (New Mexico) in the United States. Although orders of magnitude smaller in volume than rare caldera-forming supereruptions, these relatively frequent effusions of rhyolite are comparable to the largest eruptions of the 20th century, and pose a considerable volcanic hazard. However, the physical pathway from rejuvenation to eruption of silicic magma is unclear, particularly because the time between reheating of a subvolcanic intrusion and eruption is poorly quantified. This study uses nanometer-scale trace element diffusion in sanidine crystals to reveal that rejuvenation of a near-solidus or subsolidus silicic intrusion occurred in ∼10 mo or less following a protracted period (220 k.y.) of volcanic repose, and resulted in effusion of ∼3 km3 of high-silica rhyolite lava at the onset of Yellowstone’s last volcanic interval. The future renewal of effusive silicic volcanism at Yellowstone will likely require a comparable energetic intrusion of magma that remelts the shallow subvolcanic reservoir and generates eruptible rhyolite on month to annual time scales.

Complex empirical orthogonal functions analysis of ERS‐1 and TOPEX/POSEIDON combined altimetric data in the region of the Algerian current

Maps of sea level anomalies (SLA) relative to the 1993 annual mean sea level combine the data from the two altimetric missions, ERS-1 and TOPEX/POSEIDON, during the overlap period (October 1992 to December 1993). These regular maps in space and time of residual sea level every 10 days on a 0.2° regular grid are used in the region of the Algerian current where the mesoscale eddies are of primary importance to the circulation of all the Mediterranean water masses. They are first compared with ERS-1 along-track scanning radiometer sea surface temperature images to get information on two anticyclonic eddies produced by instabilities of the Algerian current and visible in both infrared and altimetric data sets. Then, an analysis of complex empirical orthogonal functions (CEOFs) is performed on the SLA data set to see the correlation of the different dynamic features of the observed variability. The CEOF analysis is applied to the complex time series formed from the original SLA time series and their Hubert transforms to separate the variability into spatially coherent modes. The spatially correlated signal in the study area (0–15°E and 35°–40°N) was found to be dominated by the first two CEOFs. These first two modes explain nearly 85% of the variability, with 80% of the total variance for the first one and 5% of the total variance for the second one. The temporal phase of the first mode indicates that a constant frequency of one cycle per year is clearly dominant, corresponding to the seasonal signal. The strongest amplitude is obtained in the southern part of the channel of Sardinia and south of the Strait of Sicily. The temporal amplitude and the temporal phase of the second mode show a periodicity of about 6 months which appears to be associated with the variability of the Algerian current as the phase isolines are parallel to the mean current path along the Algerian coast. The strongest amplitude of the second mode is located near the African coast at ∼4°E and 8°E. These two points of high variability could correspond to eddy detachments from the main current.

Peninsular terrane basement ages recorded by Paleozoic and Paleoproterozoic zircon in gabbro xenoliths and andesite from Redoubt volcano, Alaska

Historically Sactive Redoubt volcano is an Aleutian arc basalt-to-dacite cone constructed upon the Jurassic–Early Tertiary Alaska–Aleutian Range batholith. The batholith intrudes the Peninsular tectonostratigraphic terrane, which is considered to have developed on oceanic basement and to have accreted to North America, possibly in Late Jurassic time. Xenoliths in Redoubt magmas have been thought to be modern cumulate gabbros and fragments of the batholith. However, new sensitive high-resolution ion microprobe (SHRIMP) U-Pb ages for zircon from gabbro xenoliths from a late Pleistocene pyroclastic deposit are dominated by much older, ca. 310 Ma Pennsylvanian and ca. 1865 Ma Paleoproterozoic grains. Zircon age distributions and trace-element concentrations indicate that the ca. 310 Ma zircons date gabbroic intrusive rocks, and the ca. 1865 Ma zircons also are likely from igneous rocks in or beneath Peninsular terrane basement. The trace-element data imply that four of five Cretaceous–Paleocene zircons, and Pennsylvanian low-U, low-Th zircons in one sample, grew from metamorphic or hydrothermal fluids. Textural evidence of xenocrysts and a dominant population of ca. 1865 Ma zircon in juvenile crystal-rich andesite from the same pyroclastic deposit show that this basement has been assimilated by Redoubt magma. Equilibration temperatures and oxygen fugacities indicated by Fe-Ti–oxide minerals in the gabbros and crystal-rich andesite suggest sources near the margins of the Redoubt magmatic system, most likely in the magma accumulation and storage region currently outlined by seismicity and magma petrology at ∼4–10 km below sea level. Additionally, a partially melted gabbro from the 1990 eruption contains zircon with U-Pb ages between ca. 620 Ma and ca. 1705 Ma, as well as one zircon with a U-Th disequilibrium model age of 0 ka. The zircon ages demonstrate that Pennsylvanian, and probably Paleoproterozoic, igneous rocks exist in, or possibly beneath, Peninsular terrane basement. Discovery of Pennsylvanian gabbro similar in age to Skolai arc plutons 500 km to the northeast indicates that the Peninsular terrane, along with the Wrangellia and Alexander terranes, has been part of the Wrangellia composite terrane since at least Pennsylvanian time. Moreover, the zircon data suggest that a Paleoproterozoic continental fragment may be present in the mid-to-upper crust in southern Alaska.

Episodic Holocene eruption of the Salton Buttes rhyolites, California, from paleomagnetic, U‐Th, and Ar/Ar dating

In the Salton Trough, CA, five rhyolite domes form the Salton Buttes: Mullet Island, Obsidian Butte, Rock Hill, North and South Red Hill, from oldest to youngest. Results presented here include 40Ar/39Ar anorthoclase ages, 238U-230Th zircon crystallization ages, and comparison of remanent paleomagnetic directions with the secular variation curve, which indicate that all domes are Holocene. 238U-230Th zircon crystallization ages are more precise than but within uncertainty of 40Ar/39Ar anorthoclase ages, suggesting that zircon crystallization proceeded until shortly before eruption in all cases except one. Remanent paleomagnetic directions require three eruption periods: (1) Mullet Island, (2) Obsidian Butte, and (3) Rock Hill, North Red Hill, and South Red Hill. Borehole cuttings logs document up to two shallow tephra layers. North and South Red Hills likely erupted within 100 years of each other, with a combined 238U-230Th zircon isochron age of: 2.83 ± 0.60 ka (2 sigma); paleomagnetic evidence suggests this age predates eruption by hundreds of years (1800 cal BP). Rock Hill erupted closely in time to these eruptions. The Obsidian Butte 238U-230Th isochron age (2.86 ± 0.96 ka) is nearly identical to the combined Red Hill age, but its Virtual Geomagnetic Pole position suggests a slightly older age. The age of aphyric Mullet Island dome is the least well constrained: zircon crystals are resorbed and the paleomagnetic direction is most distinct; possible Mullet Island ages include ca. 2300, 5900, 6900, and 7700 cal BP. Our results constrain the duration of Salton Buttes volcanism to between ca. 5900 and 500 years.

The role of mantle‐derived magmas in the isotopic evolution of Yellowstone's magmatic system

Injection of mantle-derived magmas into the Earths crust provides the heat necessary to develop and maintain large silicic magmatic systems. However, the role of mantle-derived magmas in controlling the compositional evolution of large silicic systems remains poorly understood. Here we examine the role of mantle-derived magmas in the post-caldera magmatic system at Yellowstone Plateau, the youngest magmatism associated with the Yellowstone hotspot. Using microbeam techniques we characterize the age and Hf isotope composition of single zircon crystals hosted in rhyolites from the most recent eruptive episode at Yellowstone Plateau, which produced the Central Plateau Member rhyolites. We place these zircon data into context by comparing them to new solution Hf isotope data for the Central Plateau Member glasses, Yellowstone basalts, and potential local crustal sources. Zircons in the Central Plateau Member rhyolites record a wide range of Hf isotope compositions relative to their host melts and extend from values similar to previously erupted Yellowstone rhyolites to values similar to Yellowstone basalts. Most zircons (∼90%) are in isotopic equilibrium with their host melt, but a significant proportion show eHf values higher than their host melt, thus providing the direct evidence that silicic derivatives of mantle-derived basalts have recharged Yellowstones magmatic system. Mixing models confirm that the isotopic characteristics of the youngest Yellowstone rhyolites can be explained by recharge of Yellowstones magma reservoir with silicic derivatives of underplating, mantle-derived basalts (∼5% to 10% material added by mass). This process helps drive the long-term isotopic evolution of Yellowstones magmatic system. This article is protected by copyright. All rights reserved.

Observations on the long-period variability of the Gulf Stream downstream of Cape Hatteras

To examine the long-period variability of the Gulf Stream, sea level residuals relative to a 2-year mean sea level in the Gulf Stream downstream of Cape Hatteras (between 75°W and 60°W longitude) are used. Residuals, as derived from Geosat altimetry between November 1986 and December 1988, were gridded in space and time at a temporal resolution of 10 days and spatial resolution of 1/4°. Complex empirical orthogonal function (CEOF) analysis was applied to the data set to extract the spatially correlated signal with the original data sub sampled to 1/2°. In addition to determining the space-time scales and propagation characteristics of the different modes, wavenumber-frequency spectral techniques were used to separate the variability into propagating and stationary components. The CEOF technique applied to the data set indicated that the first four CEOF modes accounted for 60% of the variability and were found to be above the noise level 99% of the time. CEOF 1 was associated with westward propagation at 5 km/d at a wavelength of 2000 km and eastward propagation at 1–2 km/d centered at a 500-km wavelength. This first CEOF is in good agreement with thin-jet equivalent barotropic models which predict westward propagation for wavelengths greater than 1130 km. A deflection of the wavelike pattern at 65°W also indicates a possible topographic effect. A simple scaling of the effect of topography indicates that for length scales longer than the internal Rossby radius of deformation, the topographic term is at least of the same order of magnitude as the beta effect. The second CEOF was more broadbanded in wavenumber space, with eastward propagation occurring in a wavenumber-frequency band between 300 and 1400 km and 0.5 and 2.0 cycles/yr. The third CEOF is similar in structure to the first, but with less energy. CEOF 4 was clearly identifiable with higher frequencies than the first three with westward propagation at 4 km/d. The spatial location of this mode along with the westward propagation indicates possible influences from eddy-stream interactions. Thus topography, Rossby wave dynamics, and eddy-stream interactions all appear to have a significant role in determining the space-time scales and propagation properties of the long-period response of sea level in the Gulf Stream.

Fitful and protracted magma assembly leading to a giant eruption, Youngest Toba Tuff, Indonesia

The paroxysmal eruption of the 74 ka Youngest Toba Tuff (YTT) of northern Sumatra produced an extraordinary 2800 km3 of nonwelded to densely welded ignimbrite and coignimbrite ashfall. We report insights into the duration of YTT magma assembly obtained from ion microprobe U-Th and U-Pb dates, including continuous age spectra over >50% of final zircon growth, for pumices and a welded tuff spanning the compositional range of the YTT. A relatively large subpopulation of zircon crystals nucleated before the penultimate caldera-related eruption at 501 ka, but most zircons yielded interior dates 100–300 ka thereafter. Zircon nucleation and growth was likely episodic and from diverse conditions over protracted time intervals of >100 to >500 ka. Final zircon growth is evident as thin rim plateaus that are in Th/U chemical equilibrium with hosts, and that give crystallization ages within tens of ka of eruption. The longevity and chemical characteristics of the YTT zircons, as well as evidence for intermittent zircon isolation and remobilization associated with magma recharge, is especially favored at the cool and wet eutectoid conditions that characterize at least half of the YTT, wherein heat fluxes could dissolve major phases but have only a minor effect on larger zircon crystals. Repeated magma recharge may have contributed to the development of compositional zoning in the YTT but, considered together with limited allanite, quartz, and other mineral dating and geospeedometry, regular perturbations to the magma reservoir over >400 ka did not lead to eruption until 74 ka ago.

The role of crystallization-driven exsolution on the sulfur mass balance in volcanic arc magmas†

The release of large amounts of sulfur to the stratosphere during explosive eruptions affects the radiative balance in the atmosphere and consequentially impacts climate for up to several years after the event. Providing quantitative estimates for the processes that control the mass balance of sulfur between melt, crystals and vapor bubbles is needed to better understand the potential sulfur yield of individual eruption events and the conditions that favor large sulfur outputs to the atmosphere.The processes that control sulfur partitioning in magmas are (1) exsolution of volatiles (dominantly H2O) during decompression (first boiling) and during isobaric crystallization (second boiling), (2) the crystallization and breakdown of sulfide or sulfate phases in the magma and (3) the transport of sulfur-rich vapor transport (gas influx) from deeper unerupted regions of the magma reservoir. Vapor exsolution and the formation/breakdown of sulfur-rich phases can all be considered as closed system process where mass balance arguments are generally easier to constrain, whereas the contribution of sulfur by vapor transport (open system process) is more difficult to quantify. The ubiquitous “Excess Sulfur”, which refers to the much higher sulfur mass released during eruptions than what can be accounted for by the melt inclusion data (petrologic estimate), reflects the challenges in closing the sulfur mass balance between crystals, melt and vapor before and during a volcanic eruption. In this work, we try to quantify the relative importance of closed and open system processes for silicic arc volcanoes using kinetic models of sulfur partitioning during exsolution. Our calculations show that crystallization-induced exsolution (second boiling) can generate a significant fraction of the “Excess Sulfur” observed in crystal-rich arc magmas. This result does not preclude vapor migration to play an important role in the sulfur mass balance, but rather points out that second boiling (in-situ exsolution) can provide the necessary yield to drive the excess sulfur to the levels observed for these eruptions. In contrast, recharges of magma releasing sulfur-rich bubbles are necessary and most likely the primary contributor to the sulfur mass balance in silicic crystal-poor units. Finally, we apply our model to account for the effect of sulfur partitioning during second boiling and its impact on sulfur released during the Cerro Galan super-eruption in Argentina (2.08 Ma), and show the importance of second boiling in releasing a large amount of sulfur to the atmosphere during the eruption of large crystal-rich ignimbrites.

Cogenetic late Pleistocene rhyolite and cumulate diorites from Augustine Volcano revealed by SIMS 238U‐230Th dating of zircon, and implications for silicic magma generation by extraction from mush

Augustine Volcano, a frequently active andesitic island stratocone, erupted a late Pleistocene rhyolite pumice fall that is temporally linked through zircon geochronology to cumulate dioritic blocks brought to the surface in Augustines 2006 eruption. Zircon from the rhyolite yield a 238U-230Th age of ∼25 ka for their unpolished rims, and their interiors yield a bimodal age populations at ∼26 ka and a minority at ∼41 ka. Zircon from dioritic blocks, ripped from Augustines shallow magmatic plumbing system and ejected during the 2006 eruption, have interiors defining a ∼26 ka age population that is indistinguishable from that for the rhyolite; unpolished rims on the dioritic zircon are dominantly younger (≤12 ka) indicating subsequent crystallization. Zircon from rhyolite and diorite overlap in U, Hf, Ti, and REE concentrations although diorites also contain a second population of high-U, high temperature grains. Andesites that brought dioritic blocks to the surface in 2006 contain zircon with young (≤9 ka) rims and a scattering of older ages, but few zircon that crystallized during the 26 ka interval. Both the Pleistocene-age rhyolite and the 2006 dioritic inclusions plot along a whole-rock compositional trend distinct from mid-Holocene–present andesites and dacites, and the diorites, rhyolite, and two early Holocene dacites define linear unmixing trends often oblique to the main andesite array and consistent with melt (rhyolite) extraction from a mush (dacites), leaving behind a cumulate amphibole-bearing residue (diorites). Rare zircon antecrysts up to ∼300 ka from all rock types indicate that a Quaternary center has been present longer than preserved surficial deposits.

Late Pleistocene ages for the most recent volcanism and glacial‐pluvial deposits at Big Pine volcanic field, California, USA, from cosmogenic 36Cl dating

The Big Pine volcanic field is one of several Quaternary volcanic fields that poses a potential volcanic hazard along the tectonically active Owens Valley of east-central California, and whose lavas are interbedded with deposits from Pleistocene glaciations in the Sierra Nevada Range. Previous geochronology indicates an ∼1.2 Ma history of volcanism, but the eruption ages and distribution of volcanic products associated with the most-recent eruptions have been poorly resolved. To delimit the timing and products of the youngest volcanism, we combine field mapping and cosmogenic 36Cl dating of basaltic lava flows in the area where lavas with youthful morphology and well-preserved flow structures are concentrated. Field mapping and petrology reveal approximately 15 vents and 6 principal flow units with variable geochemical composition and mineralogy. Cosmogenic 36Cl exposure ages for lava flow units from the top, middle, and bottom of the volcanic stratigraphy indicate eruptions at ∼17, 27, and 40 ka, revealing several different and previously unrecognized episodes of late Pleistocene volcanism. Olivine to plagioclase-pyroxene phyric basalt erupted from several vents during the most recent episode of volcanism at ∼17 ka, and produced a lava flow field covering ∼35 km2. The late Pleistocene 36Cl exposure ages indicate that moraine and pluvial shoreline deposits that overlie or modify the youngest Big Pine lavas reflect Tioga stage glaciation in the Sierra Nevada and the shore of paleo-Owens Lake during the last glacial cycle.