Fault-zone fluids control effective normal stress and fault strength. While most earthquake models assume a fixed pore fluid pressure distribution, geologists have documented fault valving behavior, that is, cyclic changes in pressure and unsteady fluid migration along faults. Here we quantify fault valving through 2-D antiplane shear simulations of earthquake sequences on a strike-slip fault with rate-and-state friction, upward Darcy flow along a permeable fault zone, and permeability evolution. Fluid overpressure develops during the interseismic period, when healing/sealing reduces fault permeability, and is released after earthquakes enhance permeability. Coupling between fluid flow, permeability and pressure evolution, and slip produces fluid-driven aseismic slip near the base of the seismogenic zone and earthquake swarms within the seismogenic zone, as ascending fluids pressurize and weaken the fault. This model might explain observations of late interseismic fault unlocking, slow slip and creep transients, swarm seismicity, and rapid pressure/stress transmission in induced seismicity sequences.
The vibrant evolutionary patterns made by earthquake swarms are incompatible with standard, effectively two-dimensional (2D) models for general fault architecture. We leverage advances in earthquake monitoring with a deep-learning algorithm to image a fault zone hosting a 4-year-long swarm in southern California. We infer that fluids are naturally injected into the fault zone from below and diffuse through strike-parallel channels while triggering earthquakes. A permeability barrier initially limits up-dip swarm migration but ultimately is circumvented. This enables fluid migration within a shallower section of the fault with fundamentally different mechanical properties. Our observations provide high-resolution constraints on the processes by which swarms initiate, grow, and arrest. These findings illustrate how swarm evolution is strongly controlled by 3D variations in fault architecture.
more » « less- PAR ID:
- 10163086
- Publisher / Repository:
- American Association for the Advancement of Science (AAAS)
- Date Published:
- Journal Name:
- Science
- Volume:
- 368
- Issue:
- 6497
- ISSN:
- 0036-8075
- Page Range / eLocation ID:
- p. 1357-1361
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
SUMMARY Earthquakes come in clusters formed of mostly aftershock sequences, swarms and occasional foreshock sequences. This clustering is thought to result either from stress transfer among faults, a process referred to as cascading, or from transient loading by aseismic slip (pre-slip, afterslip or slow slip events). The ETAS statistical model is often used to quantify the fraction of clustering due to stress transfer and to assess the eventual need for aseismic slip to explain foreshocks or swarms. Another popular model of clustering relies on the earthquake nucleation model derived from experimental rate-and-state friction. According to this model, earthquakes cluster because they are time-advanced by the stress change imparted by the mainshock. This model ignores stress interactions among aftershocks and cannot explain foreshocks or swarms in the absence of transient loading. Here, we analyse foreshock, swarm and aftershock sequences resulting from cascades in a Discrete Fault Network model governed by rate-and-state friction. We show that the model produces realistic swarms, foreshocks and aftershocks. The Omori law, characterizing the temporal decay of aftershocks, emerges in all simulations independently of the assumed initial condition. In our simulations, the Omori law results from the earthquake nucleation process due to rate and state friction and from the heterogeneous stress changes due to the coseismic stress transfers. By contrast, the inverse Omori law, which characterizes the accelerating rate of foreshocks, emerges only in the simulations with a dense enough fault system. A high-density complex fault zone favours fault interactions and the emergence of an accelerating sequence of foreshocks. Seismicity catalogues generated with our discrete fault network model can generally be fitted with the ETAS model but with some material differences. In the discrete fault network simulations, fault interactions are weaker in aftershock sequences because they occur in a broader zone of lower fault density and because of the depletion of critically stressed faults. The productivity of the cascading process is, therefore, significantly higher in foreshocks than in aftershocks if fault zone complexity is high. This effect is not captured by the ETAS model of fault interactions. It follows that a foreshock acceleration stronger than expected from ETAS statistics does not necessarily require aseismic slip preceding the mainshock (pre-slip). It can be a manifestation of a cascading process enhanced by the topological properties of the fault network. Similarly, earthquake swarms might not always imply transient loading by aseismic slip, as they can emerge from stress interactions.
-
more » « less
Abstract The Kolumbo submarine volcano in the southern Aegean (Greece) is associated with repeated seismic unrest since at least two decades and the causes of this unrest are poorly understood. We present a ten‐month long microseismicity data set for the period 2006–2007. The majority of earthquakes cluster in a cone‐shaped portion of the crust below Kolumbo. The tip of this cone coincides with a low Vp‐anomaly at 2–4 km depth, which is interpreted as a crustal melt reservoir. Our data set includes several earthquake swarms, of which we analyze the four with the highest events numbers in detail. Together the swarms form a zone of fracturing elongated in the SW‐NE direction, parallel to major regional faults. All four swarms show a general upward migration of hypocenters and the cracking front propagates unusually fast, compared to swarms in other volcanic areas. We conclude that the swarm seismicity is most likely triggered by a combination of pore‐pressure perturbations and the re‐distribution of elastic stresses. Fluid pressure perturbations are induced likely by obstructions in the melt conduits in a rheologically strong layer between 6 and 9 km depth. We conclude that the zone of fractures below Kolumbo is exploited by melts ascending from the mantle and filling the crustal melt reservoir. Together with the recurring seismic unrest, our study suggests that a future eruption is probable and monitoring of the Kolumbo volcanic system is highly advisable.
-
Earthquake swarms are ubiquitous in volcanic systems, being manifestations of underlying nontectonic processes such as magma intrusions or volatile fluid transport. The Long Valley caldera, California, is one such setting where episodic earthquake swarms and persistent uplift suggest the presence of active magmatism. We quantify the long-term spatial and temporal characteristics of seismicity in the region using cluster analysis on a 25-year high-resolution earthquake catalog derived using leading-edge deep-learning algorithms. Our results show that earthquake swarms beneath the caldera exhibit enlarged families with statistically significant tendency for upward migration patterns. The ascending swarms tend to nucleate at the base of the seismogenic zone with a spatial footprint that is laterally constrained by the southern rim of the caldera. We suggest that these swarms are driven by the transport of volatile-rich fluids released from deep volcanic processes. The observations highlight the potential for extreme spatial segmentation of earthquake triggering processes in magmatic systems.more » « less
-
more » « less
Abstract In 2013 and 2018, earthquake swarms with a maximum moment magnitude of 4.5 occurred ~5 km from the northern section of the Dead Sea Transform Fault. Here we show that aquifer pressure data, interferometric synthetic aperture radar surface deformation time series, and seismic monitoring suggest that groundwater withdrawal triggered these earthquakes. Continuous groundwater extraction from several wells located ~10 km west of the swarms has accelerated since 2010 and resulted in a total decrease of ~50 m of the groundwater level at the time of the 2018 earthquake swarm. The withdrawal also corresponds to surface subsidence of ~10 mm/year based on repeat interferometric synthetic aperture radar measurements. The temporal correlation, extensive subsidence, anomalous swarm characteristics, and normal faulting orientation suggest a connection between the groundwater withdrawal and recent earthquakes. Poroelastic modeling demonstrates that pumping‐induced pore pressure decrease west of the earthquake could have caused significant dilatational stresses that led to normal faulting events outside the aquifer.
