
Possible control of subduction zone slow-earthquake periodicity by silica enrichment
- Select a language for the TTS:
- UK English Female
- UK English Male
- US English Female
- US English Male
- Australian Female
- Australian Male
- Language selected: (auto detect) - EN
Play all audios:

ABSTRACT Seismic and geodetic observations in subduction zone forearcs indicate that slow earthquakes, including episodic tremor and slip, recur at intervals of less than six months to more
than two years1,2. In Cascadia, slow slip is segmented along strike3 and tremor data show a gradation from large, infrequent slip episodes to small, frequent slip events with increasing
depth of the plate interface4. Observations5,6,7 and models8,9 of slow slip and tremor require the presence of near-lithostatic pore-fluid pressures in slow-earthquake source regions;
however, direct evidence of factors controlling the variability in recurrence times is elusive. Here we compile seismic data from subduction zone forearcs exhibiting recurring slow
earthquakes and show that the average ratio of compressional (P)-wave velocity to shear (S)-wave velocity (_v_P/_v_S) of the overlying forearc crust ranges between 1.6 and 2.0 and is
linearly related to the average recurrence time of slow earthquakes. In northern Cascadia, forearc _v_P/_v_S values decrease with increasing depth of the plate interface and with decreasing
tremor-episode recurrence intervals. Low _v_P/_v_S values require a large addition of quartz in a mostly mafic forearc environment10,11. We propose that silica enrichment varying from 5 per
cent to 15 per cent by volume from slab-derived fluids and upward mineralization in quartz veins12 can explain the range of observed _v_P/_v_S values as well as the downdip decrease in
_v_P/_v_S. The solubility of silica depends on temperature13, and deposition prevails near the base of the forearc crust11. We further propose that the strong temperature dependence of
healing and permeability reduction in silica-rich fault gouge via dissolution–precipitation creep14 can explain the reduction in tremor recurrence time with progressive silica enrichment.
Lower gouge permeability at higher temperatures leads to faster fluid overpressure development and low effective fault-normal stress, and therefore shorter recurrence times. Our results also
agree with numerical models of slip stabilization under fault zone dilatancy strengthening15 caused by decreasing fluid pressure as pore space increases. This implies that
temperature-dependent silica deposition, permeability reduction and fluid overpressure development control dilatancy and slow-earthquake behaviour. Access through your institution Buy or
subscribe This is a preview of subscription content, access via your institution ACCESS OPTIONS Access through your institution Subscribe to this journal Receive 51 print issues and online
access $199.00 per year only $3.90 per issue Learn more Buy this article * Purchase on SpringerLink * Instant access to full article PDF Buy now Prices may be subject to local taxes which
are calculated during checkout ADDITIONAL ACCESS OPTIONS: * Log in * Learn about institutional subscriptions * Read our FAQs * Contact customer support SIMILAR CONTENT BEING VIEWED BY OTHERS
PHYSICAL CONDITIONS AND FRICTIONAL PROPERTIES IN THE SOURCE REGION OF A SLOW-SLIP EVENT Article 06 May 2021 OSCILLATIONS IN FLUID PRESSURE CAUSED BY SILICA PRECIPITATION IN A FRACTURE
Article Open access 20 February 2025 A LIKELY GEOLOGICAL RECORD OF DEEP TREMOR AND SLOW SLIP EVENTS FROM A SUBDUCTED CONTINENTAL BROKEN FORMATION Article Open access 16 March 2022 REFERENCES
* Schwartz, S. Y. & Rokosky, J. M. Slow slip events and seismic tremor at circum-Pacific subduction zones. _Rev. Geophys._ 45, RG3004 (2007) ADS Google Scholar * Beroza, G. C. &
Ide, S. Slow earthquakes and nonvolcanic tremor. _Annu. Rev. Earth Planet. Sci._ 39, 271–296 (2011) ADS CAS Google Scholar * Brudzinski, M. R. & Allen, R. M. Segmentation in episodic
tremor and slip all along Cascadia. _Geology_ 35, 907–910 (2007) ADS Google Scholar * Wech, A. G. & Creager, K. C. A continuum of stress, strength and slip in the Cascadia subduction
zone. _Nature Geosci._ 4, 624–628 (2011) ADS CAS Google Scholar * Shelly, D. R., Beroza, G. C., Ide, S. & Nakamula, S. Low-frequency earthquakes in Shikoku, Japan, and their
relationship to episodic tremor and slip. _Nature_ 442, 188–191 (2006) ADS CAS PubMed Google Scholar * Audet, P., Bostock, M. G., Christensen, N. I. & Peacock, S. M. Seismic evidence
for overpressured subducting oceanic crust and megathrust fault sealing. _Nature_ 457, 76–78 (2009) ADS CAS PubMed Google Scholar * Song, T.-R. A. et al. Subducting slab ultra-slow
velocity layer coincident with silent earthquakes in southern Mexico. _Science_ 324, 502–506 (2009) ADS CAS PubMed Google Scholar * Liu, Y. & Rice, J. R. Spontaneous and triggered
aseismic deformation transients in a subduction fault model. _J. Geophys. Res._ 112, B09404 (2007) ADS Google Scholar * Segall, P., Rubin, A. M., Bradley, A. M. & Rice, J. R. Dilatant
strengthening as a mechanism for slow slip events. _J. Geophys. Res._ 115, B12305 (2010) ADS Google Scholar * Christensen, N. I. Poisson’s ratio and crustal seismology. _J. Geophys. Res._
101, 3139–3156 (1996) ADS CAS Google Scholar * Ramachandran, K. & Hyndman, R. D. The fate of fluids released from subducting slab in northern Cascadia. _Solid Earth_ 3, 121–129 (2012)
ADS Google Scholar * Breeding, C. M. & Ague, J. J. Slab-derived fluids and quartz-vein formation in an accretionary prism, Otago Schist, New Zealand. _Geology_ 30, 499–502 (2002) ADS
CAS Google Scholar * Manning, C. E. The solubility of quartz in H2O in the lower crust and upper mantle. _Geochim. Cosmochim. Acta_ 58, 4831–4839 (1994) ADS CAS Google Scholar *
Giger, S. B., Tenthorey, E., Cox, S. F. & Fitz Gerald, J. D. Permeability evolution in quartz gouges under hydrothermal conditions. _J. Geophys. Res._ 112, B07202 (2007) ADS Google
Scholar * Liu, Y. Numerical simulations on megathrust rupture stabilization under strong dilatancy strengthening in slow slip region. _Geophys. Res. Lett._ 40, 1311–1316 (2013) ADS Google
Scholar * Shen, Z.-K. et al. Pole-tide modulation of slow slip events at circum-Pacific subduction zones. _Bull. Seismol. Soc. Am._ 95, 2009–2015 (2005) Google Scholar * Lowry, A. R.
Resonant slow fault slip in subduction zones forced by climatic load stress. _Nature_ 442, 802–805 (2006) ADS CAS PubMed Google Scholar * Pollitz, F. F., Wech, A. G., Kao, H. &
Bürgmann, R. Annual modulation of non-volcanic tremor in northern Cascadia. _J. Geophys. Res._ 118, 2445–2459 (2013) ADS Google Scholar * Hansen, R. T. J., Bostock, M. G. &
Christensen, N. I. Nature of the low velocity zone in Cascadia from receiver function waveform inversion. _Earth Planet. Sci. Lett._ 337–388, 25–38 (2012) ADS Google Scholar * Audet, P.
& Schwartz, S. Y. Hydrologic control of forearc strength and seismicity in the Costa Rican subduction zone. _Nature Geosci._ 6, 852–855 (2013) ADS CAS Google Scholar * Audet, P. et
al. Slab morphology in the Cascadia fore arc and its relation to episodic tremor and slip. _J. Geophys. Res._ 115, B00A16 (2010) Google Scholar * Peacock, S. M., Christensen, N. I.,
Bostock, M. G. & Audet, P. High pore pressures and porosity at 35 km depth in the Cascadia subduction zone. _Geology_ 39, 471–474 (2011) ADS Google Scholar * Lowry, A. R. &
Pérez-Gussinyé, M. The role of crustal quartz in controlling Cordilleran deformation. _Nature_ 471, 353–357 (2011) ADS CAS PubMed Google Scholar * Kern, H. Elastic-wave velocity in
crustal and mantle rocks at high pressure and temperature: the role of the high-low quartz transition and of dehydration reactions. _Phys. Earth Planet. Inter._ 29, 12–23 (1982) ADS CAS
Google Scholar * Sillitoe, R. H. Major gold deposits and belts of the North and South American Cordillera: distribution, tectonomagmatic settings, and metallogenic considerations. _Econ.
Geol._ 103, 663–687 (2008) CAS Google Scholar * Wada, I. Wang. K., He, J. & Hyndman, R. D. Weakening of the subduction interface and its effects on surface heat flow, slab dehydration,
and mantle wedge serpentinization. _J. Geophys. Res._ 113, B04402 (2008) ADS Google Scholar * Dempsey, D. E., Rowland, J. V., Zyvoloski, G. A. & Archer, R. A. Modeling the effects of
silica deposition and fault rupture on natural geothermal systems. _J. Geophys. Res._ 117, B05207 (2012) ADS Google Scholar * Fisher, D. M. & Brantley, S. L. The role of silica
redistribution in the evolution of slip instabilities along subduction interfaces: constraints from the Kodiak accretionary complex, Alaska. _J. Struct. Geol._ doi:10.1016/j.jsg.2014.03.010.
(in the press) ADS Google Scholar * Frederiksen, A. W. & Bostock, M. G. Modelling teleseismic waves in dipping, anisotropic structures. _Geophys. J. Int._ 141, 401–412 (2000) ADS
Google Scholar * Hayes, G. P., Wald, D. J. & Johnson, R. L. Slab1.0: a three-dimensional model of global subduction zone geometries. _J. Geophys. Res._ 117, B01302 (2012) ADS Google
Scholar * Frederiksen, A. W., Folsom, H. & Zandt, G. Neighbourhood inversion of teleseismic _Ps_ conversions for anisotropy and layer dip. _Geophys. J. Int._ 155, 200–212 (2003) ADS
Google Scholar * Wallace, L. M. & Beavan, J. Diverse slow slip behaviour at the Hukirangi subduction margin, New Zealand. _J. Geophys. Res._ 115, B12402 (2010) ADS Google Scholar *
Heki, K. & Kataoka, T. On the biannually repeating slow-slip events at the Ryukyu Trench, southwestern Japan. _J. Geophys. Res._ 113, B11402 (2008) ADS Google Scholar * Jiang, Y. et
al. Slow slip events in Costa Rica detected by continuous GPS observations, 2002-2011. _Geochem. Geophys. Geosyst._ 13, Q04006 (2012) ADS Google Scholar * Obara, K. Characteristics and
interactions between non-volcanic tremor and related slow earthquakes in the Nankai subduction zone, southwest Japan. _J. Geodyn._ 52, 229–248 (2011) Google Scholar * Vergnolle, M. et al.
Slow slip events in Mexico revised from the processing of 11 year GPS observations. _J. Geophys. Res._ 115, B08403 (2010) ADS Google Scholar * DeShon, H. R. et al. Seismogenic zone
structure beneath the Nicoya Peninsula, Costa Rica, from 3D local earthquake P- and S- wave tomography. _Geophys. J. Int._ 164, 109–124 (2006) ADS Google Scholar * Reyners, M.,
Eberhart-Phillips, D., Stuart, G. & Nishimura, Y. Imaging subduction from the trench to 300 km depth beneath the central North Island, New Zealand, with Vp and Vp/Vs. _Geophys. J. Int._
165, 565–583 (2006) ADS Google Scholar * Matsubara, M., Obara, K. & Kasahara, K. High Vp/Vs zone accompanying non-volcanic tremors and slow slip events beneath southwestern Japan.
_Tectonophysics_ 472, 6–17 (2009) ADS Google Scholar * Huesca-Perez, E. & Husker, A. Shallow travel-time tomography below southern Mexico. _Geofis. Int._ 51, 281–291 (2012) Google
Scholar * Ansell, J. H. & Bannister, S. C. Shallow morphology of the subducted Pacific plate along the Hikurangi margin. _Phys. Earth Planet. Inter._ 93, 3–20 (1996) ADS Google Scholar
Download references ACKNOWLEDGEMENTS Data used in this study come from the Japan Meteorological Agency, New Zealand National Seismograph Network, and the Tectonic Observatory (Caltech).
Funding for this work comes from the Natural Science and Engineering Research Council (Canada) and the Miller Institute for Basic Research in Science (University of California, Berkeley). We
thank N. Beeler, A. Rubin and P. Ampuero for discussions and comments. AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * Department of Earth Sciences, University of Ottawa, Ottawa, Ontario K1N
6N5, Canada, Pascal Audet * Department of Earth and Planetary Science and Berkeley Seismological Laboratory, University of California, Berkeley, 94720-4767, California, USA Roland Bürgmann
Authors * Pascal Audet View author publications You can also search for this author inPubMed Google Scholar * Roland Bürgmann View author publications You can also search for this author
inPubMed Google Scholar CONTRIBUTIONS P.A. performed data processing and inversion. Both P.A. and R.B. contributed to the interpretations and preparation of the final manuscript.
CORRESPONDING AUTHOR Correspondence to Pascal Audet. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing financial interests. EXTENDED DATA FIGURES AND TABLES EXTENDED
DATA FIGURE 1 EXAMPLES OF RECEIVER FUNCTIONS AND INVERSION RESULTS FOR EACH SUBDUCTION ZONE. In A–C the slab contours from refs 30 and 41 are in yellow, contours of slow-slip patches from
ref. 2 are in light green, contours and epicentres of tremors from ref. 2 are in purple, and station locations used in this study are shown as inverted red triangles. For a subset of
stations (PLAY, PXZ and IGK, identified by the blue squares on the maps) we show the observed (top, A) and modelled (bottom, B) radial receiver functions ordered by back-azimuth and, for
each back-azimuth, by slowness of the incoming P wave. C, A slice through model misfits, with warm colours indicating low values, showing the _v_P/_v_S of forearc crust versus the _v_P/_v_S
of the low-velocity zone. The star shows the minimum value of the misfit plot (best-fitting value). POWERPOINT SLIDES POWERPOINT SLIDE FOR FIG. 1 POWERPOINT SLIDE FOR FIG. 2 POWERPOINT SLIDE
FOR FIG. 3 RIGHTS AND PERMISSIONS Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Audet, P., Bürgmann, R. Possible control of subduction zone slow-earthquake periodicity by
silica enrichment. _Nature_ 510, 389–392 (2014). https://doi.org/10.1038/nature13391 Download citation * Received: 03 December 2013 * Accepted: 08 April 2014 * Published: 18 June 2014 *
Issue Date: 19 June 2014 * DOI: https://doi.org/10.1038/nature13391 SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a
shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative