Distant earthquakes can trigger the eruption of mud volcanoes. We document the response of the Davis‐Schrimpf mud volcanoes, California, to two earthquakes and non‐response to four additional events. ...We show that the Davis‐Schrimpf mud volcanoes are more sensitive to long period seismic waves than to shorter period waves of the same amplitude. Our observations are consistent with models for dislodging bubbles and particles by time varying flows produced by seismic waves. Mobilizing trapped bubbles or particles increases permeability or fluid mobility, increasing discharge.
Key Points
Mud volcanoes sometimes respond to earthquakes
In addition to dissipated energy, the frequency of shaking affects response
Triggering occurs by mechanisms that depend on frequency
Young flood channels emanate from Cerberus Fossae, Mars. In order to determine whether subsurface aquifers can discharge water sufficiently fast to explain these flood features, I develop a model ...that couples groundwater flow in a sub‐cryosphere aquifer with flow through a fissure that penetrates from the surface to the aquifer. The model is constrained by estimates of peak discharge, and the distance water flowed uphill from the fissure. Aquifer permeabilities similar to those of young basalt aquifers on the Earth, ∼10−9 m2, can produce inferred discharges of 106 m3/s Burr et al., 2002a but probably require multiple floods to create the channels.
Magma mixing can occur in a fluid manner to produce banded pumice or in a brittle manner to form enclaves. We propose that the critical control on mixing style is a competition between developing ...networks of crystals in the intruding magma that impart a strength to the magma and melting and disrupting those networks in the host. X-ray computed tomography analysis demonstrates that banded pumice from the 1915 Mt. Lassen eruption lacks crystal networks. In contrast, rhyodacite hosts with mafic enclaves from Chaos Crags contain well-developed networks of large crystals. We present a one-dimensional conductive cooling model that predicts mixing style, either ductile or brittle, as a function of magma compositions, temperatures, and the size of the intruding dike. Our model relies on three assumptions: (1) Mixing is initiated by the injection of a hot dike into a cooler magma body with a yield strength; (2) when magma crystallinity exceeds a critical value, 13 vol% plagioclase, the magma develops a yield strength; and (3) when total crystallinity exceeds 40 vol%, the magma has a penetrative crystal network and is effectively solid. Importantly, because the two magmas are of different compositions, their crystallinities and viscosities do not have the same variations with temperature. As the intruding magma cools, it crystallizes from the outside in, while simultaneously, host magma temperature near the intruder rises. Mixing of the two magmas begins when the host magma is heated sufficiently to (1) disrupt the crystal network and (2) initiate convection. If the shear stress exerted by the convecting host magma on the dike is greater than the yield strength of the dike margin (and dike crystallinity does not exceed 40 %), then fluid mixing occurs, otherwise enclaves form by brittle deformation of the dike. Application of the model to magma compositions representative of Lassen and Chaos Crags shows that emplacement of dikes <1 m thick should produce enclaves, whereas thicker dikes should generate fluid mixing and form banded pumice within days to weeks of emplacement. Similar relationships apply to other modeled magmatic systems, including Pinatubo, Unzen, and Ksudach/Shtuybel’ volcanoes. For all studied systems, the absolute size of the intruding dike, not just its proportion relative to the host, influences mixing style.
Large, sustained well water level changes (>10 cm) in response to distant (more than hundreds of kilometers) earthquakes have proven enigmatic for over 30 years. Here we use high sampling rates at a ...well near Grants Pass, Oregon, to perform the first simultaneous analysis of both the dynamic response of water level and sustained changes, or steps. We observe a factor of 40 increase in the ratio of water level amplitude to seismic wave ground velocity during a sudden coseismic step. On the basis of this observation we propose a new model for coseismic pore pressure steps in which a temporary barrier deposited by groundwater flow is entrained and removed by the more rapid flow induced by the seismic waves. In hydrothermal areas, this mechanism could lead to 4 × 10−2 MPa pressure changes and triggered seismicity.
A long-standing conceptual model for deep submarine eruptions is that high hydrostatic pressure hinders degassing and acceleration, and suppresses magma fragmentation. The 2012 submarine rhyolite ...eruption of Havre volcano in the Kermadec arc provided constraints on critical parameters to quantitatively test these concepts. This eruption produced a >1 km3 raft of floating pumice and a 0.1 km3 field of giant (>1 m) pumice clasts distributed down-current from the vent. We address the mechanism of creating these clasts using a model for magma ascent in a conduit. We use water ingestion experiments to address why some clasts float and others sink. We show that at the eruption depth of 900 m, the melt retained enough dissolved water, and hence had a low enough viscosity, that strain-rates were too low to cause brittle fragmentation in the conduit, despite mass discharge rates similar to Plinian eruptions on land. There was still, however, enough exsolved vapor at the vent depth to make the magma buoyant relative to seawater. Buoyant magma was thus extruded into the ocean where it rose, quenched, and fragmented to produce clasts up to several meters in diameter. We show that these large clasts would have floated to the sea surface within minutes, where air could enter pore space, and the fate of clasts is then controlled by the ability to trap gas within their pore space. We show that clasts from the raft retain enough gas to remain afloat whereas fragments from giant pumice collected from the seafloor ingest more water and sink. The pumice raft and the giant pumice seafloor deposit were thus produced during a clast-generating effusive submarine eruption, where fragmentation occurred above the vent, and the subsequent fate of clasts was controlled by their ability to ingest water.
•Havre magma entered the ocean before fragmenting.•Clasts were produced by quenching buoyant magma in the ocean.•Buoyant >1 m diameter pumice blocks floated to the ocean surface.•Clasts with enough isolated porosity and trapped gas floated in a raft while the rest sank.
The El Tatio geothermal field in the Chilean Altiplano contains hydrothermal silica sinter deposits overlaying glacial and volcanic units, providing an opportunity to constrain the timing of ...deglaciation and volcanic activity in an area with sparse absolute chronologies. We obtained 51 new radiocarbon ages and δ13C values on the organic material trapped in these sinter deposits. Based on the δ13C values, we exclude 29 samples for possible contamination with bacterial mats that incorporate old carbon. We infer that hydrothermal activity initiated ~27 ka ago and has been nearly continuous ever since. The ages of the oldest sinter deposits coincide with ages of moraines that stabilized after the most recent deglaciation. Whereas late Pleistocene sinters are broadly distributed in the field, Holocene deposits are found around active hydrothermal features. Although recent volcanism is absent in the vicinity of El Tatio, persistent hydrothermal discharge implies a long‐lived magmatic heat source.
Plain Language Summary
The El Tatio geothermal field in the Chilean Altiplano contains large deposits of silica precipitated by hot springs, called sinter or geyserite. Silica sinter traps organic material that is suitable for radiocarbon dating; however, previous attempts using this methodology have shown inconsistencies that still need to be critically assessed. Here, we found that bacteria can contaminate the radiocarbon measurements leading to older radiocarbon ages. Stable carbon isotopes (13C and 12C) in dated material can be used to evaluate when the organic matter is contaminated and thus provides an unreliable age. Sinter deposits in El Tatio overlay glacial deposits and volcanic units; thus, sinter ages provide a bound for volcanism and deglaciation in the area. Uncontaminated radiocarbon ages indicate that El Tatio is a long‐lived geothermal system that started at least ~27 ka ago and has been active since. However, the spatial distribution of thermal water discharge in the field has changed over time. The onset of hydrothermal activity is consistent with the onset of deglaciation in the Chilean Altiplano. Even though recent volcanic deposits are absent in the vicinity of El Tatio, the persistent discharge of hot water implies a significant long‐lived magmatic source providing heat to the geothermal system.
Key Points
El Tatio is a long‐lived hydrothermal system, active for at least 27 ka
Dated sinter deposits provide bounds on recent volcanic activity and deglaciation in the vicinity of El Tatio
Bacterial mats trapped in sinter lead to old apparent radiocarbon ages
Spherulites are spherical clusters of radiating crystals that occur naturally in rhyolitic obsidian. The growth of spherulites requires diffusion and uptake of crystal forming components from the ...host rhyolite melt or glass, and rejection of non-crystal forming components from the crystallizing region. Water concentration profiles measured by synchrotron-source Fourier transform spectroscopy reveal that water is expelled into the surrounding matrix during spherulite growth, and that it diffuses outward ahead of the advancing crystalline front. We compare these profiles to models of water diffusion in rhyolite to estimate timescales for spherulite growth. Using a diffusion-controlled growth law, we find that spherulites can grow on the order of days to months at temperatures above the glass transition. The diffusion-controlled growth law also accounts for spherulite size distribution, spherulite growth below the glass transition, and why spherulitic glasses are not completely devitrified.
Earthquake-induced crustal deformation and ground shaking can alter stream flow and water levels in wells through consolidation of surficial deposits, fracturing of solid rocks, aquifer deformation, ...and the clearing of fracture-filling material. Although local conditions affect the type and amplitude of response, a compilation of reported observations of hydrological response to earthquakes indicates that the maximum distance to which changes in stream flow and water levels in wells have been reported is related to earthquake magnitude. Detectable streamflow changes occur in areas within tens to hundreds of kilometers of the epicenter, whereas changes in groundwater levels in wells can occur hundreds to thousands of kilometers from earthquake epicenters.
Continuous monitoring of in situ reservoir responses to stress transients provides insights into the evolution of geothermal reservoirs. By exploiting the stress dependence of seismic velocity ...changes, we investigate the temporal evolution of the reservoir stress state of the Salton Sea geothermal field (SSGF), California. We find that the SSGF experienced a number of sudden velocity reductions (~0.035 to 0.25%) that are most likely caused by openings of fractures due to dynamic stress transients (as small as 0.08 MPa and up to 0.45 MPa) from local and regional earthquakes. Depths of velocity changes are estimated to be about 0.5 to 1.5 km, similar to the depths of the injection and production wells. We derive an empirical in situ stress sensitivity of seismic velocity changes by relating velocity changes to dynamic stresses. We also observe systematic velocity reductions (0.04 to 0.05%) during earthquake swarms in mid-November 2009 and late-December 2010. On the basis of volumetric static and dynamic stress changes, the expected velocity reductions from the largest earthquakes with magnitude ranging from 3 to 4 in these swarms are less than 0.02%, which suggests that these earthquakes are likely not responsible for the velocity changes observed during the swarms. Instead, we argue that velocity reductions may have been induced by poroelastic opening of fractures due to aseismic deformation. We also observe a long-term velocity increase (~0.04%/year) that is most likely due to poroelastic contraction caused by the geothermal production. Our observations demonstrate that seismic interferometry provides insights into in situ reservoir response to stress changes.