Seismic imaging of Santorini: subsurface constraints on caldera collapse and present-day magma recharge

Volcanic calderas are surface depressions formed by roof collapse following evacuation of magma from an underlying reservoir. The mechanisms of caldera formation are debated and predict differences in the evolution of the caldera floor and distinct styles of magma recharge. Here we use a dense, active source, seismic tomography study to reveal the sub-surface physical properties of the Santorini caldera in order to understand caldera formation. We find a ∼3-km-wide, cylindrical low-velocity anomaly in the upper 3 km beneath the north-central portion of the caldera, that lies directly above the pressure source of the 2011-2012 inflation. We interpret this anomaly as a low-density volume caused by excess porosities of between 4% and 28%, with pore spaces filled with hot seawater. Vents that were formed during the first three phases of the 3.6 ka Late Bronze Age (LBA) eruption are located close to the edge of the imaged structure. The correlation between older volcanic vents and the low-velocity anomaly suggests that this feature may be long-lived. We infer that collapse of a limited area of the caldera floor resulted in a high-porosity, low-density cylindrical volume, which formed by either chaotic collapse along reverse faults, wholesale subsidence and infilling with tuffs and ignimbrites, phreatomagmatic fracturing, or a combination of these processes. Phase 4 eruptive vents are located along the margins of the topographic caldera and the velocity structure indicates that coherent down-drop of the wider topographic caldera followed the more limited collapse in the northern caldera. This progressive collapse sequence is consistent with models for multi-stage formation of nested calderas along conjugate reverse and normal faults. The upper crustal density differences inferred from the seismic velocity model predict differences in subsurface gravitational loading that correlate with the location of 2011-2012 edifice inflation. This result supports the hypothesis that sub-surface density anomalies may influence present-day magma recharge events. We postulate that past collapses and the resulting topographical and density variations at Santorini influence magma focusing between eruptive cycles, a feedback process that may be important in other volcanoes

Protracted timescales of lower crustal growth at the fast-spreading East Pacific Rise

Author Posting. © The Author(s), 2011. This is the author's version of the work. It is posted here by permission of Nature Publishing Group for personal use, not for redistribution. The definitive version was published in Nature Geoscience 5 (2012): 275-278, doi:10.1038/ngeo1378.Formation of the oceanic crust at mid-ocean ridges is a fundamental component of plate tectonics. A majority of the crust at many ridges is composed of plutonic rocks that form by crystallization of mantle-derived magmas within the crust. Recent application of U/Pb dating to samples from in-situ oceanic crust has begun to provide exciting new insight into the timing, duration and distribution of magmatism during formation of the plutonic crust1-4. Previous studies have focused on samples from slow-spreading ridges, however, the time scales and processes of crustal growth are expected to vary with plate spreading rate. Here we present the first high-precision dates from plutonic crust formed at the fast-spreading East Pacific Rise (EPR). Individual zircon minerals yielded dates from 1.420–1.271 million years ago, with uncertainties of ± 0.006–0.081 million years. Within individual samples, zircons record a range of dates of up to ~0.124 million years, consistent with protracted crystallization or assimilation of older zircons from adjacent rocks. The variability in dates is comparable to data from the Vema lithospheric section on the Mid-Atlantic Ridge (MAR)3, suggesting that time scales of magmatic processes in the lower crust may be similar at slow- and fast-spreading ridges.This research was partially funded by NSF grant OCE-0727914 (SAB), a Cardiff University International Collaboration Award (CJL) and NERC grant NE/C509023/1 (CJM).2012-07-2

Mantle thermal pulses below the Mid-Atlantic Ridge and temporal variations in the formation of oceanic lithosphere

A 20-Myr record of creation of oceanic lithosphere is exposed along a segment of the central Mid-Atlantic Ridge on an uplifted sliver of lithosphere. The degree of melting of the mantle that is upwelling below the ridge, estimated from the chemistry of the exposed mantle rocks, as well as crustal thickness inferred from gravity measurements, show oscillations of ∼3–4 Myr superimposed on a longer-term steady increase with time. The time lag between oscillations of mantle melting and crustal thickness indicates that the mantle is upwelling at an average rate of ∼25 mm yr-1, but this appears to vary through time. Slow-spreading lithosphere seems to form through dynamic pulses of mantle upwelling and melting, leading not only to along-axis segmentation but also to across-axis structural variability. Also, the central Mid-Atlantic Ridge appears to have become steadily hotter over the past 20 Myr, possibly owing to north–south mantle flow