A high-precision age estimate of the Holocene Plinian eruption of Mount Mazama, Oregon, USA

© The Author(s) 2015 The climactic eruption of Mount Mazama in Oregon, North America, resulted in the deposition of the most widespread Holocene tephra deposit in the conterminous United States and south-western Canada. The tephra forms an isochronous marker horizon for palaeoenvironmental, sedimentary and archaeological reconstructions, despite the current lack of a precise age estimate for the source eruption. Previous radiocarbon age estimates for the eruption have varied, and Greenland ice-core ages are in disagreement. For the Mazama tephra to be fully utilised in tephrochronology and palaeoenvironmental research, a refined (precise and accurate) age for the eruption is required. Here, we apply a meta-analysis of all previously published radiocarbon age estimations (n = 81), and perform Bayesian statistical modelling to this data set, to provide a refined age of 7682–7584 cal. yr BP (95.4% probability range). Although the depositional histories of the published ages vary, this estimate is consistent with that estimated from the GISP2 ice-core of 7627 ± 150 yr BP (Zdanowicz et al., 1999)

Petrology and Geochemistry of the Galápagos Islands : portrait of a pathological mantle plume

We report new major element, trace element, isotope ratio, and geochronological data on the Galapagos Archipelago. Magmas erupted from the large western volcanos are generally moderately fractionated tholeiites of uniform composition; those erupted on other islands are compositionally diverse, ranging from tholeiites to picritic basanitoids. While these volcanoes do not form a strictly linear age progressive chain, the ages of the oldest dated flows on any given volcano do form a reasonable progression from youngest in the west to oldest in the east, consistent with motion of the Nazca plate with respect to the fixed hotspot reference frame. Isotope ratios in the Galapagos display a considerable range, from values typical of mid-ocean ridge basalt on Genovesa (⁸⁷Sr/⁸⁶Sr: 0.70259, εNd: +9.4, ²⁰⁶Pb/²⁰⁴Pb: 18.44), to typical oceanic island values on Floreana (⁸⁷Sr/⁸⁶Sr: 0.70366, εNd: +5.2, ²⁰⁶Pb/²⁰⁴Pb: 20.0). La/SmN ranges from 0.45 to 6.7; other incompatible element abundances and ratios show comparable ranges. Isotope and incompatible element ratios define a horseshoe pattern with the most depleted signatures in the center of the Galapagos Archipelago and the more enriched signatures on the eastern, northern, and southern periphery. These isotope and incompatible element patterns appear to reflect thermal entrainment of asthenosphere by the Galapagos plume as it experiences velocity shear in the uppermost asthenosphere. Both north-south heterogeneity within the plume itself and regional variations in degree and depth of melting also affect magma compositions. Rare earth systematics indicate that melting beneath the Galapagos begins in the garnet peridotite stability field, except beneath the southern islands, where melting may occur entirely in the spinel peridotite stability field. The greatest degree of melting occurs beneath the central western volcanos and decreases both to the east and to the north and south. Si₈.₀, Fe₈.₀, and Na₈.₀, values are generally consistent with these inferences. This suggests that interaction between the plume and surrounding asthenosphere results in significant cooling of the plume. Superimposed on this thermal pattern produced by plume-asthenosphere interaction is a tendency for melting to be less extensive and to occur at shallower depths to the south, presumably reflecting a decrease in ambient asthenospheric temperatures away from the Galapagos Spreading Center.Copyrighted by American Geophysical Union

Rheology of magmas with bimodal crystal size and shape distributions: insights from analog experiments

Magmas in volcanic conduits commonly contain microlites in association with preexisting phenocrysts, as often indicated by volcanic rock textures. In this study, we present two different experiments that inves- tigate the flow behavior of these bidisperse systems. In the first experiments, rotational rheometric methods are used to determine the rheology of monodisperse and polydisperse suspensions consisting of smaller, prolate particles (microlites) and larger, equant particles (phenocrysts) in a bubble‐free Newtonian liquid (silicate melt). Our data show that increasing the relative proportion of prolate microlites to equant pheno- crysts in a magma at constant total particle content can increase the relative viscosity by up to three orders of magnitude. Consequently, the rheological effect of particles in magmas cannot be modeled by assuming a monodisperse population of particles. We propose a new model that uses interpolated parameters based on the relative proportions of small and large particles and produces a considerably improved fit to the data than earlier models. In a second series of experiments we investigate the textures produced by shearing bimodal suspensions in gradually solidifying epoxy resin in a concentric cylinder setup. The resulting textures show the prolate particles are aligned with the flow lines and spherical particles are found in well‐organized strings, with sphere‐depleted shear bands in high‐shear regions. These observations may explain the measured variation in the shear thinning and yield stress behavior with increasing solid fraction and particle aspect ratio. The implications for magma flow are discussed, and rheological results and tex- tural observations are compared with observations on natural samples

Shear-induced pressure changes and seepage phenomena in a deforming porous layer-I

We present a model for flow and seepage in a deforming, shear-dilatant sensitive porous layer that enables estimates of the excess pore fluid pressures and flow rates in both the melt and solid phase to be captured simultaneously as a function of stress rate. Calculations are relevant to crystallizing magma in the solidosity range 0.5–0.8 (50–20 per cent melt), corresponding to a dense region within the solidification front of a crystallizing magma chamber. Composition is expressed only through the viscosity of the fluid phase, making the model generally applicable to a wide range of magma types. A natural scaling emerges that allows results to be presented in non-dimensional form. We show that all length-scales can be expressed as fractions of the layer height H, timescales as fractions of H2(nβ'θ+ 1)/(θk) and pressures as fractions of . Taking as an example the permeability k in the mush of the order of magnitude 1015 m2 Pa1 s1, a layer thickness of tens of metres and a mush strength (θ) in the range 108–1012 Pa, an estimate of the consolidation time for near-incompressible fluids is of the order of 105–109 s. Using mush permeability as a proxy, we show that the greatest maximum excess pore pressures develop consistently in rhyolitic (high-viscosity) magmas at high rates of shear ( , implying that during deformation, the mechanical behaviour of basaltic and rhyolitic magmas will differ. Transport parameters of the granular framework including tortuosity and the ratio of grain size to layer thickness (a/H) will also exert a strong effect on the mechanical behaviour of the layer at a given rate of strain. For dilatant materials under shear, flow of melt into the granular layer is implied. Reduction in excess pore pressure sucks melt into the solidification front at a velocity proportional to the strain rate. For tectonic rates (generally 1014 s1), melt upwelling (or downwelling, if the layer is on the floor of the chamber) is of the order of cm yr1. At higher rates of loading comparable with emplacement of some magmatic intrusions (1010 s1), melt velocities may exceed effects due to instabilities resulting from local changes in density and composition. Such a flow carries particulates with it, and we speculate that these may become trapped in the granular layer depending on their sizes. If on further solidification the segregated grain size distribution of the particulates is frozen in the granular layer, structure formation including layering and grading may result. Finally, as the process settles down to a steady state, the pressure does not continue to decrease. We find no evidence for critical rheological thresholds, and the process is stable until so much shear has been applied that the granular medium fails, but there is no hydraulic failure

Late Silurian plutons in Yucatan

This is the published version. Copyright 1996 American Geophysical Union. All Rights Reserved.U-Pb measurements of zircons from two composite plutons in the Maya Mountains of the Yucatan Block (Belize) give Late Silurian ages. Zircons from one of the five compositional phases of the Mountain Pine Ridge pluton yield an age of 418±3.6 Ma. A second compositional phase gives a minimum age of 404 Ma, and zircons from a third phase, although plagued with high common Pb, yield ages consistent with the other two. Zircons from one compositional phase of the Hummingbird-Mullins River pluton indicate an age of about 410–420 Ma. These data demonstrate that two of the three Maya Mountains plutons residing among the strata of the Late Pennsylvanian through Permian Santa Rosa Group are older than that sedimentation. Although the third pluton was not dated, both the similarity of sedimentary facies patterns adjacent to it to those adjacent to one of the plutons dated as Late Silurian and a published single Rb-Sr age of 428 ± 41 Ma suggest this third pluton also was emergent during Santa Rosa deposition. Thus the new U/Pb dates and other data suggest that all three Maya Mountains plutons pre-date Late Carboniferous sedimentation and that none intrude the Santa Rosa Group. Although very uniform ages of about 230 Ma amongst all plutons, derived from abundant earlier dating by the K-Ar system, led to the conclusion that intrusion mostly had occurred in the Late Triassic, the U-Pb ages (obtained from the same sites as the K-Ar dates) demonstrate that the K-Ar ages do not derive from a Late Triassic intrusive episode. The K-Ar dates probably are a signature of the rifting associated with Pangean breakup and formation of the Gulf of Mexico. In a reconstructed Pangea, the position of the Maya Mountains Late Silurian plutons suggests that the Late Silurian Acadian-Caledonian orogen of eastern North America extended through the region of the future Gulf of Mexico. Finally, the U-Pb ages of the Maya Mountains plutons are the same as those of a group of shocked zircons found in the Chicxulub impact structure and its fallout layer. The presence of these ages in both locations suggests that the Maya Mountains exposures may be representative of the basement of the Yucatan Block, hence of the basement impacted by the Chicxulub object

A criterion for the fragmentation of bubbly magma based on brittle failure theory

The fragmentation of bubbly magma is a defining point in a volcanic eruption-before fragmentation the magma flows relatively slowly, during fragmentation the bubbles break up to release compressed gas and, afterwards, the eruption becomes a violent gas flow carrying suspended magma particles. Seemingly benign lava flows or domes can suddenly fragment into deadly pyroclastic flows(1-3). Several criteria have been proposed to define the point of magma fragmentation or foam stability(4-7). The criterion of Papale(7) is based on melt relaxation theory and equates magma strain rate with the rate of increase of flow velocity with distance. It ignores, however, the role of bubble pressure in causing fragmentation. Two empirical approaches(4,5) consider the role of high bubble pressure in causing fragmentation but do not address the underlying physics of magma fragmentation. Here I develop a fragmentation criterion for bubbly magma based on brittle failure theory and apply it to the fragmentation of lava domes and flows. On the basis of this theory, a bubbly magma will fragment when the tensile stress at the inner walls of bubbles exceeds the tensile strength of the magma. The fragmentation conditions depend strongly on initial water content, with calculated vesicularity and final water levels coinciding reasonably well with those in observed pumices. This suggests that the proposed criterion captures the essence of the fragmentation process in bubbly magma.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/62724/1/402648a0.pd

Submarine Fernandina : magmatism at the leading edge of the Galapagos hot spot

Author Posting. © American Geophysical Union, 2006. 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 7 (2006): Q12007, doi:10.1029/2006GC001290.New multibeam and side-scan sonar surveys of Fernandina volcano and the geochemistry of lavas provide clues to the structural and magmatic development of Galápagos volcanoes. Submarine Fernandina has three well-developed rift zones, whereas the subaerial edifice has circumferential fissures associated with a large summit caldera and diffuse radial fissures on the lower slopes. Rift zone development is controlled by changes in deviatoric stresses with increasing distance from the caldera. Large lava flows are present on the gently sloping and deep seafloor west of Fernandina. Fernandina's submarine lavas are petrographically more diverse than the subaerial suite and include picrites. Most submarine glasses are similar in composition to aphyric subaerially erupted lavas, however. These rocks are termed the “normal” series and are believed to result from cooling and crystallization in the subcaldera magma system, which buffers the magmas both thermally and chemically. These normal-series magmas are extruded laterally through the flanks of the volcano, where they scavenge and disaggregate olivine-gabbro mush to produce picritic lavas. A suite of lavas recovered from the terminus of the SW submarine rift and terraces to the south comprises evolved basalts and icelandites with MgO = 3.1 to 5.0 wt.%. This “evolved series” is believed to form by fractional crystallization at 3 to 5 kb, involving extensive crystallization of clinopyroxene and titanomagnetite in addition to plagioclase. “High-K” lavas were recovered from the southwest rift and are attributed to hybridization between normal-series basalt and evolved-series magma. The geochemical and structural findings are used to develop an evolutionary model for the construction of the Galápagos Platform and better understand the petrogenesis of the erupted lavas. The earliest stage is represented by the deep-water lava flows, which over time construct a broad submarine platform. The deep-water lavas originate from the subcaldera plumbing system of the adjacent volcano. After construction of the platform, eruptions focus to a point source, building an island with rift zones extending away from the adjacent, buttressing volcanoes. Most rift zone magmas intrude laterally from the subcaldera magma chamber, although a few evolve by crystallization in the upper mantle and deep crust.This work was supported by the National Science Foundation grants OCE0002818 and EAR0207605 (D.G.), OCE0002461 (D.J.F. and M.K.), OCE9811504 (D.J.F. and M.R.P.), and EAR0207425 (K.H.) and WHOI postdoctoral support for Soule

Primitive layered gabbros from fast-spreading lower oceanic crust

Three-quarters of the oceanic crust formed at fast-spreading ridges is composed of plutonic rocks whose mineral assemblages, textures and compositions record the history of melt transport and crystallization between the mantle and the sea floor. Despite the importance of these rocks, sampling them in situ is extremely challenging owing to the overlying dykes and lavas. This means that models for understanding the formation of the lower crust are based largely on geophysical studies and ancient analogues (ophiolites) that did not form at typical mid-ocean ridges. Here we describe cored intervals of primitive, modally layered gabbroic rocks from the lower plutonic crust formed at a fast-spreading ridge, sampled by the Integrated Ocean Drilling Program at the Hess Deep rift. Centimetre-scale, modally layered rocks, some of which have a strong layering-parallel foliation, confirm a long-held belief that such rocks are a key constituent of the lower oceanic crust formed at fast-spreading ridges. Geochemical analysis of these primitive lower plutonic rocks-in combination with previous geochemical data for shallow-level plutonic rocks, sheeted dykes and lavas-provides the most completely constrained estimate of the bulk composition of fast-spreading oceanic crust so far. Simple crystallization models using this bulk crustal composition as the parental melt accurately predict the bulk composition of both the lavas and the plutonic rocks. However, the recovered plutonic rocks show early crystallization of orthopyroxene, which is not predicted by current models of melt extraction from the mantle and mid-ocean-ridge basalt differentiation. The simplest explanation of this observation is that compositionally diverse melts are extracted from the mantle and partly crystallize before mixing to produce the more homogeneous magmas that erupt

Late Cenozoic tephrostratigraphy offshore the southern Central American Volcanic Arc: 1. Tephra ages and provenance

We studied the tephra inventory of 18 deep sea drill sites from six DSDP/ODP legs (Legs 84, 138, 170, 202, 205, 206) and two IODP legs (Legs 334 and 344) offshore the southern Central American Volcanic Arc (CAVA). Eight drill sites are located on the incoming Cocos plate and ten drill sites on the continental slope of the Caribbean plate. In total we examined ∼840 ash-bearing horizons and identified ∼650 of these as primary ash beds of which 430 originated from the CAVA. Correlations of ash beds were established between marine cores and with terrestrial tephra deposits, using major and trace element glass compositions with respect to relative stratigraphic order. As a prerequisite for marine-terrestrial correlations we present a new geochemical data set for significant Neogene and Quaternary Costa Rican tephras. Moreover, new Ar/Ar ages for marine tephras have been determined and marine ash beds are also dated using the pelagic sedimentation rates. The resulting correlations and provenance analyses build a tephrochronostratigraphic framework for Costa Rica and Nicaragua that covers the last >8 Myr. We define 39 correlations of marine ash beds to specific tephra formations in Costa Rica and Nicaragua; from the 4.15 Ma Lower Sandillal Ignimbrite to the 3.5 ka Rincón de la Vieja Tephra from Costa Rica, as well as another 32 widely distributed tephra layers for which their specific region of origin along Costa Rica and Nicaragua can be constrained