INDUCED SEISMICITY. Seismicity triggered by fluid injection-induced aseismic slip.

PubMed

Guglielmi, Yves; Cappa, Frédéric; Avouac, Jean-Philippe; Henry, Pierre; Elsworth, Derek

2015-06-12

Anthropogenic fluid injections are known to induce earthquakes. The mechanisms involved are poorly understood, and our ability to assess the seismic hazard associated with geothermal energy or unconventional hydrocarbon production remains limited. We directly measure fault slip and seismicity induced by fluid injection into a natural fault. We observe highly dilatant and slow [~4 micrometers per second (μm/s)] aseismic slip associated with a 20-fold increase of permeability, which transitions to faster slip (~10 μm/s) associated with reduced dilatancy and micro-earthquakes. Most aseismic slip occurs within the fluid-pressurized zone and obeys a rate-strengthening friction law μ = 0.67 + 0.045ln(v/v₀) with v₀ = 0.1 μm/s. Fluid injection primarily triggers aseismic slip in this experiment, with micro-earthquakes being an indirect effect mediated by aseismic creep. Copyright © 2015, American Association for the Advancement of Science.

Seafloor seismological/geodetic observations in the rupture area of the 2011 Tohoku-oki Earthquake

NASA Astrophysics Data System (ADS)

Hino, Ryota; Shinohara, Masanao; Ito, Yoshihiro

2016-04-01

A number of important aspects of the 2011 Tohoku-oki earthquake (Mw 9.0) were clarified by the seafloor seismological and geodetic observation above the rupture area of the earthquake. Besides the extraordinarily large coseismic displacements, various kinds of slow slip phenomena associated with intensive micro-seismicity on the plate boundary fault were identified by near field ocean bottom seismographs and seafloor geodetic observation networks. The Tohoku-oki earthquake was preceded by evident foreshock activity with a spatial expansion of this seismicity. The activity became significantly intense after the occurrence of the largest foreshock two days before the mainshock rupture. During the period, clear continuous seafloor deformation was identified caused by the aseismic slip following the largest foreshock. Another different type of aseismic slip event had occurred before this pre-imminent activity had started about a month before the largest foreshock happened. The observed increased seismicity associated with aseismic slip suggests that there must have been some chain reaction like interplay of seismic and interseismic slips before the large earthquake broke out. However, no evident deformation signals were observed indicating acceleration of fault slip immediately before the mainshock. Seafloor geodetic measurements reveals that the postseismic deformation around the rupture area of the Tohoku-oki earthquake shows complex spatial pattern and the complexity is mostly due to significant viscoelastic relaxation induced by the huge coseismic slip. The effects of viscoelastic deformation makes it difficult to identify the deformation associated with the after slip or regaining of interplate coupling and requires us to enhance the abilities of seafloor monitoring to detect the slip activities on the fault. We started an array of seismometer arrays observation including broad-band seismographs to detect and locate slow-slip events and low-frequency tremors. Another observation we started is direct-path acoustic ranging across the trench axis. Slip rate of the shallow fault can be measured by monitoring the change in distance between the benchmarks on the incoming and overrding plates.

Earthquake Clustering on Normal Faults: Insight from Rate-and-State Friction Models

NASA Astrophysics Data System (ADS)

Biemiller, J.; Lavier, L. L.; Wallace, L.

2016-12-01

Temporal variations in slip rate on normal faults have been recognized in Hawaii and the Basin and Range. The recurrence intervals of these slip transients range from 2 years on the flanks of Kilauea, Hawaii to 10 kyr timescale earthquake clustering on the Wasatch Fault in the eastern Basin and Range. In addition to these longer recurrence transients in the Basin and Range, recent GPS results there also suggest elevated deformation rate events with recurrence intervals of 2-4 years. These observations suggest that some active normal fault systems are dominated by slip behaviors that fall between the end-members of steady aseismic creep and periodic, purely elastic, seismic-cycle deformation. Recent studies propose that 200 year to 50 kyr timescale supercycles may control the magnitude, timing, and frequency of seismic-cycle earthquakes in subduction zones, where aseismic slip transients are known to play an important role in total deformation. Seismic cycle deformation of normal faults may be similarly influenced by its timing within long-period supercycles. We present numerical models (based on rate-and-state friction) of normal faults such as the Wasatch Fault showing that realistic rate-and-state parameter distributions along an extensional fault zone can give rise to earthquake clusters separated by 500 yr - 5 kyr periods of aseismic slip transients on some portions of the fault. The recurrence intervals of events within each earthquake cluster range from 200 to 400 years. Our results support the importance of stress and strain history as controls on a normal fault's present and future slip behavior and on the characteristics of its current seismic cycle. These models suggest that long- to medium-term fault slip history may influence the temporal distribution, recurrence interval, and earthquake magnitudes for a given normal fault segment.

Preslip and cascade processes initiating laboratory stick slip

USGS Publications Warehouse

McLaskey, Gregory C.; Lockner, David A.

2014-01-01

Recent modeling studies have explored whether earthquakes begin with a large aseismic nucleation process or initiate dynamically from the rapid growth of a smaller instability in a “cascade-up” process. To explore such a case in the laboratory, we study the initiation of dynamic rupture (stick slip) of a smooth saw-cut fault in a 76mm diameter cylindrical granite laboratory sample at 40–120MPa confining pressure. We use a high dynamic range recording system to directly compare the seismic waves radiated during the stick-slip event to those radiated from tiny (M _6) discrete seismic events, commonly known as acoustic emissions (AEs), that occur in the seconds prior to each large stick slip. The seismic moments, focal mechanisms, locations, and timing of the AEs all contribute to our understanding of their mechanics and provide us with information about the stick-slip nucleation process. In a sequence of 10 stick slips, the first few microseconds of the signals recorded from stick-slip instabilities are nearly indistinguishable from those of premonitory AEs. In this sense, it appears that each stick slip begins as an AE event that rapidly (~20 μs) grows about 2 orders of magnitude in linear dimension and ruptures the entire 150mm length of the simulated fault. We also measure accelerating fault slip in the final seconds before stick slip. We estimate that this slip is at least 98% aseismic and that it both weakens the fault and produces AEs that will eventually cascade-up to initiate the larger dynamic rupture.

Detecting aseismic strain transients from seismicity data

USGS Publications Warehouse

Llenos, A.L.; McGuire, J.J.

2011-01-01

Aseismic deformation transients such as fluid flow, magma migration, and slow slip can trigger changes in seismicity rate. We present a method that can detect these seismicity rate variations and utilize these anomalies to constrain the underlying variations in stressing rate. Because ordinary aftershock sequences often obscure changes in the background seismicity caused by aseismic processes, we combine the stochastic Epidemic Type Aftershock Sequence model that describes aftershock sequences well and the physically based rate- and state-dependent friction seismicity model into a single seismicity rate model that models both aftershock activity and changes in background seismicity rate. We implement this model into a data assimilation algorithm that inverts seismicity catalogs to estimate space-time variations in stressing rate. We evaluate the method using a synthetic catalog, and then apply it to a catalog of M???1.5 events that occurred in the Salton Trough from 1990 to 2009. We validate our stressing rate estimates by comparing them to estimates from a geodetically derived slip model for a large creep event on the Obsidian Buttes fault. The results demonstrate that our approach can identify large aseismic deformation transients in a multidecade long earthquake catalog and roughly constrain the absolute magnitude of the stressing rate transients. Our method can therefore provide a way to detect aseismic transients in regions where geodetic resolution in space or time is poor. Copyright 2011 by the American Geophysical Union.

Dual megathrust slip behaviors of the 2014 Iquique earthquake sequence

NASA Astrophysics Data System (ADS)

Meng, Lingsen; Huang, Hui; Bürgmann, Roland; Ampuero, Jean Paul; Strader, Anne

2015-02-01

The transition between seismic rupture and aseismic creep is of central interest to better understand the mechanics of subduction processes. A Mw 8.2 earthquake occurred on April 1st, 2014 in the Iquique seismic gap of northern Chile. This event was preceded by a long foreshock sequence including a 2-week-long migration of seismicity initiated by a Mw 6.7 earthquake. Repeating earthquakes were found among the foreshock sequence that migrated towards the mainshock hypocenter, suggesting a large-scale slow-slip event on the megathrust preceding the mainshock. The variations of the recurrence times of the repeating earthquakes highlight the diverse seismic and aseismic slip behaviors on different megathrust segments. The repeaters that were active only before the mainshock recurred more often and were distributed in areas of substantial coseismic slip, while repeaters that occurred both before and after the mainshock were in the area complementary to the mainshock rupture. The spatiotemporal distribution of the repeating earthquakes illustrates the essential role of propagating aseismic slip leading up to the mainshock and illuminates the distribution of postseismic afterslip. Various finite fault models indicate that the largest coseismic slip generally occurred down-dip from the foreshock activity and the mainshock hypocenter. Source imaging by teleseismic back-projection indicates an initial down-dip propagation stage followed by a rupture-expansion stage. In the first stage, the finite fault models show an emergent onset of moment rate at low frequency (< 0.1 Hz), while back-projection shows a steady increase of high frequency power (> 0.5 Hz). This indicates frequency-dependent manifestations of seismic radiation in the low-stress foreshock region. In the second stage, the rupture expands in rich bursts along the rim of a semi-elliptical region with episodes of re-ruptures, suggesting delayed failure of asperities. The high-frequency rupture remains within an area of local high trench-parallel gravity anomaly (TPGA), suggesting the presence of subducting seamounts that promote high-frequency generation. Our results highlight the complexity of the interactions between large-scale aseismic slow-slip and dynamic ruptures of megathrust earthquakes.

Spectrum of Slip Processes on the Subduction Interface in a Continuum Framework Resolved by Rate-and State Dependent Friction and Adaptive Time Stepping

NASA Astrophysics Data System (ADS)

Herrendoerfer, R.; van Dinther, Y.; Gerya, T.

2015-12-01

To explore the relationships between subduction dynamics and the megathrust earthquake potential, we have recently developed a numerical model that bridges the gap between processes on geodynamic and earthquake cycle time scales. In a self-consistent, continuum-based framework including a visco-elasto-plastic constitutive relationship, cycles of megathrust earthquake-like ruptures were simulated through a purely slip rate-dependent friction, albeit with very low slip rates (van Dinther et al., JGR, 2013). In addition to much faster earthquakes, a range of aseismic slip processes operate at different time scales in nature. These aseismic processes likely accommodate a considerable amount of the plate convergence and are thus relevant in order to estimate the long-term seismic coupling and related hazard in subduction zones. To simulate and resolve this wide spectrum of slip processes, we innovatively implemented rate-and state dependent friction (RSF) and an adaptive time-stepping into our continuum framework. The RSF formulation, in contrast to our previous friction formulation, takes the dependency of frictional strength on a state variable into account. It thereby allows for continuous plastic yielding inside rate-weakening regions, which leads to aseismic slip. In contrast to the conventional RSF formulation, we relate slip velocities to strain rates and use an invariant formulation. Thus we do not require the a priori definition of infinitely thin, planar faults in a homogeneous elastic medium. With this new implementation of RSF, we succeed to produce consistent cycles of frictional instabilities. By changing the frictional parameter a, b, and the characteristic slip distance, we observe a transition from stable sliding to stick-slip behaviour. This transition is in general agreement with predictions from theoretical estimates of the nucleation size, thereby to first order validating our implementation. By incorporating adaptive time-stepping based on a fraction of characteristic slip distance over maximum slip velocity, we are able to resolve stick-slip events and increase computational speed. In this better resolved framework, we examine the role of aseismic slip on the megathrust cycle and its dependence on subduction velocity.

Slow slip pulses driven by thermal pressurization of pore fluid: theory and observational constraints

NASA Astrophysics Data System (ADS)

Garagash, D.

2012-12-01

We discuss recently developed solutions for steadily propagating self-healing slip pulses driven by thermal pressurization (TP) of pore fluid [Garagash, 2012] on a fault with a constant sliding friction. These pulses are characterized by initial stage of undrained weakening of the fault (when fluid/heat can not yet escape the frictionally heated shear zone), which gives way to partial restrengthening due to increasing hydrothermal diffusion under conditions of diminished rate of heating, leading to eventual locking of the slip. The rupture speed of these pulses is decreasing function of the thickness (h) of the principal shear zone. We find that "thick" shear zones, h >> hdyna, where hdyna = (μ/τ0) (ρc/fΛ)(4α/cs), can support aseismic TP pulses propagating at a fraction hdyna/h of the shear wave speed cs, while "thin" shear zones, h˜hdyna or thinner, can only harbor seismic slip. (Here μ - shear modulus, τ0 - the nominal fault strength, f - sliding friction, ρc - the heat capacity of the fault gouge, Λ - the fluid thermal pressurization factor, α - hydrothermal diffusivity parameter of the gouge). For plausible range of fault parameters, hdyna is between 10s to 100s of micrometers, suggesting that slow slip transients propagating at 1 to 10 km/day may occur in the form of a TP slip pulse accommodated by a meter-thick shear zone. We verify that this is, indeed, a possibility by contrasting the predictions for aseismic, small-slip TP pulses operating at seismologically-constrained, near-lithostatic pore pressure (effective normal stress ≈ 3 to 10 MPa) with the observations (slip duration at a given fault location ≈ week, propagation speed ≈ 15 km/day, and the inferred total slip ≈ 2 to 3 cm) for along-strike propagation of the North Cascadia slow slip events of '98-99 [Dragert et al., 2001, 2004]. Furthermore, we show that the effect of thermal pressurization on the strength of the subduction interface is comparable to or exceeds that of the rate-dependence of friction, previously suggested as a mechanism for aseismic transients [e.g., Liu and Rice, 2009; Segall et al., 2010], if the frictional properties of gabbro [He et al., 2007] under the hydrothermal conditions for the North Cascadia slab [Hacker et al., 2003] are used. It therefore appears that while some friction weakening with the slip rate may be required to nucleate a slow slip event in the first place, thermal pressurization mechanism has to be included in realistic models of dynamics of aseismic slip transients, as long as the source of the transients is linked to the conditionally-stable part of the interface (with near velocity-neutral friction). The results of this study point to the importance of the principal shear zone thickness during a slip event and its possible change with the slip rate [e.g., Platt et al., AGU FM 2010]. The insight into how stable creep or a slow slip event may transition into a seismic rupture and how an earthquake rupture "selects" its principal shear zone, which is shown to largely define the TP slip dynamics, may require addressing the slip localization as a phenomena concurrent to the development of transient slip, and therefore coupled to other relevant source mechanisms.

Is Slow Slip a Cause or a Result of Tremor?

NASA Astrophysics Data System (ADS)

Luo, Y.; Ampuero, J. P.

2017-12-01

While various modeling efforts have been conducted to reproduce subsets of observations of tremor and slow-slip events (SSE), a fundamental but yet unanswered question is whether slow slip is a cause or a result of tremor. Tremor is commonly regarded as driven by SSE. This view is mainly based on observations of SSE without detected tremors and on (frequency-limited) estimates of total tremor seismic moment being lower than 1% of their concomitant SSE moment. In previous studies we showed that models of heterogeneous faults, composed of seismic asperities embedded in an aseismic fault zone matrix, reproduce quantitatively the hierarchical patterns of tremor migration observed in Cascadia and Shikoku. To address the title question, we design two end-member models of a heterogeneous fault. In the SSE-driven-tremor model, slow slip events are spontaneously generated by the matrix (even in the absence of seismic asperities) and drive tremor. In the Tremor-driven-SSE model the matrix is stable (it slips steadily in the absence of asperities) and slow slip events result from the collective behavior of tremor asperities interacting via transient creep (local afterslip fronts). We study these two end-member models through 2D quasi-dynamic multi-cycle simulations of faults governed by rate-and-state friction with heterogeneous frictional properties and effective normal stress, using the earthquake simulation software QDYN (https://zenodo.org/record/322459). We find that both models reproduce first-order observations of SSE and tremor and have very low seismic to aseismic moment ratio. However, the Tremor-driven-SSE model assumes a simpler rheology than the SSE-driven-tremor model and matches key observations better and without fine tuning, including the ratio of propagation speeds of forward SSE and rapid tremor reversals and the decay of inter-event times of Low Frequency Earthquakes. These modeling results indicate that, in contrast to a common view, SSE could be a result of tremor activity. We also find that, despite important interactions between asperities, tremor activity rates are proportional to the underlying aseismic slip rate, supporting an approach to estimate SSE properties with high spatial-temporal resolutions via tremor activity.

A kinematic model of patchy slip at depth explains observed tremor waveforms on the San Andreas fault near Parkfield, California

NASA Astrophysics Data System (ADS)

Gottschaemmer, E.; Harrington, R. M.; Cochran, E. S.; Bohlen, T.

2011-12-01

Recent observations of both triggered and ambient tremor suggest that tremor results from simple shear-failure. Tremor episodes on the San Andreas fault near Parkfield are thought to be comprised of clusters of individual events with frequencies between 2-8 Hz. Such low frequency earthquakes (LFEs) occur at depths where the frictional properties of the fault surface are primarily slip-strengthening with imbedded patches of slip weakening material that slip seismically when the surrounding fault creeps in a slow-slip event. Here we show new tremor waveforms from a temporary deployment of 13 broadband seismometers spaced at a maximum on the order of 30 km near Cholame, California are consistent with a series of small seismically slipping patches surrounded by an aseismic region along a fault surface. We model individual seismic events kinematically as small shear failures (M ~ 1) at depths exceeding 15 km. We use stress drop values of 1 MPa, based on a slip to fault area ratio. We simulate tremor recorded at the surface by our temporary array centered near Cholame, for frequencies up to 8 Hz using a staggered-grid finite-difference scheme to solve the elastic equations of motion, and the 3D velocity and density model from Thurber et al. (2006). Our simulations indicate that multiple seismically slipping patches in an aseismic region successfully recreate tremor characteristics observed in multiple studies, including individual tremor bursts, individual events, and episodic behavior. The kinematic model presented here will help to constrain the distribution and amplitude of the seismically slipping patches at depth, which will then be used in a dynamic model with variable frictional properties.

Constraints on Slow Slip from Landsliding and Faulting

NASA Astrophysics Data System (ADS)

Delbridge, Brent Gregory

The discovery of slow-slip has radically changed the way we understand the relative movement of Earth's tectonic plates and the accumulation of stress in fault zones that fail in large earthquakes. Prior to the discovery of slow-slip, faults were thought to relieve stress either through continuous aseismic sliding, as is the case for continental creeping faults, or in near instantaneous failure. Aseismic deformation reflects fault slip that is slow enough that both inertial forces and seismic radiation are negligible. The durations of observed aseismic slip events range from days to years, with displacements of up to tens of centimeters. These events are not unique to a specific depth range and occur on faults in a variety of tectonic settings. This aseismic slip can sometimes also trigger more rapid slip somewhere else on the fault, such as small embedded asperities. This is thought to be the mechanism generating observed Low Frequency Earthquakes (LFEs) and small repeating earthquakes. I have preformed a series of studies to better understanding the nature of tectonic faulting which are compiled here. The first is entitled "3D surface deformation derived from airborne interferometric UAVSAR: Application to the Slumgullion Landslide", and was originally published in the Journal of Geophysical Research in 2016. In order to understand how landslides respond to environmental forcing, we quantify how the hydro-mechanical forces controlling the Slumgullion Landslide express themselves kinematically in response to the infiltration of seasonal snowmelt. The well-studied Slumgullion Landslide, which is 3.9 km long and moves persistently at rates up to 2 cm/day is an ideal natural laboratory due to its large spatial extent and rapid deformation rates. The lateral boundaries of the landslide consist of strike-slip fault features, which over time have built up large flank ridges. The second study compiled here is entitled "Temporal variation of intermediate-depth earthquakes around the time of the M9.0 Tohoku-oki earthquake" and was originally published in Geophysical Research Letters in 2017. The temporal evolution of intermediate depth seismicity before and after the 2011 M 9.0 Tohoku-oki earthquake reveals interactions between plate interface slip and deformation in the subducting slab. I investigate seismicity rate changes in the upper and lower planes of the double seismic zone beneath northeast Japan. The average ratio of upper plane to lower plane activity and the mean deep aseismic slip rate both increased by factor of two. An increase of down-dip compression in the slab resulting from coseismic and postseismic deformation enhanced seismicity in the upper plane, which is dominated by events accommodating down-dip shortening from plate unbending. In the third and final study included here I use geodetic measurements to place a quantitative upper bound on the size of the slow slip accompanying large bursts of quasi-periodic tremors and LFEs on the Parkfield section of the SAF. We use a host of analysis methods to try to isolate the small signal due to the slow slip and characterize noise properties. We find that in addition to subduction zones, transform faults are also capable of producing ETSs. However, given the upper-bounds from our analysis, surface geodetic measurements of this slow slip is likely to remain highly challenging.

Transient Aseismic Slip in the Cascadia Subduction Zone: From Monitoring to Useful Real-time Hazards Information

NASA Astrophysics Data System (ADS)

Roeloffs, E. A.; Beeler, N. M.

2010-12-01

The Cascadia subduction zone, extending from northern California to Vancouver Island, has a 10,000 year record of earthquakes > M8.5 at intervals of several hundred years, with the last major event (~M9) in 1700. Agencies in CA, OR, WA, and BC are raising public awareness of the hazards posed by a repeat Cascadia earthquake and its ensuing tsunami. Because most of the subduction interface is now seismically quiet, an interface event M6 or larger would generate intense public concern that it could be a potential foreshock of a great earthquake. Cascadia residents are also interested in the episodic tremor and slip (ETS) events that recur months to years apart: strong evidence implies these aseismic events represent accelerated interface slip downdip of the seismogenic zone. Simple mechanics implies ETS events temporarily increase the stressing rate on the locked zone. ETS events in northern Cascadia recur at fairly regular intervals and produced roughly similar patterns of deformation. However, an unusually large ETS event or increased interface seismicity would certainly prompt public officials and local residents to expect scientists to quickly determine the implications for a major Cascadia earthquake. Earthquake scientists generally agree that such “situations of concern” warrant close monitoring, but attempts to quantify potential probability changes are in very early stages. With >30 borehole strainmeters and >100 GPS stations of the NSF-funded Plate Boundary Observatory (PBO) in Cascadia, geodesists must develop a well-organized real-time monitoring scheme for interpreting aseismic deformation, with an accompanying public communication strategy. Two previously-exercised monitoring and communication protocols could be adapted for Cascadia. During the Parkfield, California, Earthquake Experiment, geodetic signals were assigned alert levels based on their rareness in the past record, on confirmation by more than one instrument, and on consistency with aseismic slip near the expected Parkfield hypocenter. In Japan, the Tokai Earthquake Prediction Experiment assesses strainmeter and tiltmeter anomalies based on their consistency with slip near the anticipated nucleation point of the next Tokai earthquake. Earthquake scientists unfamiliar with these two projects often presume that releasing uncertain “pre-event” information will have negative consequences, such as dangerous, unnecessary evacuations. Emergency managers are better qualified to plan effective communication, but many have experience only with post-earthquake information, and multi-state and international discussions are stymied by lack of funds for non-federal officials to travel outside their states or countries. Both the Parkfield and Tokai efforts have included pre-planning with emergency management officials. The critical public message remains that communities must plan for major Cascadia earthquakes to occur without warning. But every effort should still be made to recognize a foreshock or aseismic precursor, which could save lives in Cascadia coastal communities facing tsunami impact 10-20 minutes after a seismic rupture. Even if no pre-earthquake signals are observed, geodetic data will track post-seismic deformation, which may contain clues to the timing of large aftershocks.

Low-frequency earthquakes in Shikoku, Japan, and their relationship to episodic tremor and slip.

PubMed

Shelly, David R; Beroza, Gregory C; Ide, Satoshi; Nakamula, Sho

2006-07-13

Non-volcanic seismic tremor was discovered in the Nankai trough subduction zone in southwest Japan and subsequently identified in the Cascadia subduction zone. In both locations, tremor is observed to coincide temporally with large, slow slip events on the plate interface downdip of the seismogenic zone. The relationship between tremor and aseismic slip remains uncertain, however, largely owing to difficulty in constraining the source depth of tremor. In southwest Japan, a high quality borehole seismic network allows identification of coherent S-wave (and sometimes P-wave) arrivals within the tremor, whose sources are classified as low-frequency earthquakes. As low-frequency earthquakes comprise at least a portion of tremor, understanding their mechanism is critical to understanding tremor as a whole. Here, we provide strong evidence that these earthquakes occur on the plate interface, coincident with the inferred zone of slow slip. The locations and characteristics of these events suggest that they are generated by shear slip during otherwise aseismic transients, rather than by fluid flow. High pore-fluid pressure in the immediate vicinity, as implied by our estimates of seismic P- and S-wave speeds, may act to promote this transient mode of failure. Low-frequency earthquakes could potentially contribute to seismic hazard forecasting by providing a new means to monitor slow slip at depth.

The Effect of Semi-Brittle Rheology on the Seismicity at the Subduction Interface: Coseismic and Aseismic Events

NASA Astrophysics Data System (ADS)

Tong, X.; Lavier, L.

2017-12-01

Cold and warm subduction zones usually have different seismicity and tectonic structure. Aseismic events like episodic tremor and slip (ETS) and slow slip event (SSE) are often observed in warm and young slabs which typically have less megathrust seismicity and smaller seismogenic area (e.g. southwest Japan). On the other hand, cold and old slabs (e.g. Northeast Japan) have more megathrust events and larger seismogenic area and few aseismic events. Recent studies have try to model the differences in seismic behaviors with different approaches, includes rheological heterogeneity (e.g. frictional vs. viscous), petrological heterogeneity (e.g. hydration-dehydration process and mineral phase changes), and the frictional heterogeneity (e.g. rate-and-state dependent friction). Following previous works, we proposed a new model in which the subduction channel has a temperature dependent material assembly which composed of an explicit mixture of basalt/eclogite and mantle peridotite. Our model also take into account rate and state dependent friction and pore fluid pressure. Depending on the temperature, the basalt and peridotite mixture can behave either as an elastoplastic frictional or a Maxwell viscoelastic material. To model the mixture numerically, we use DynEarthSol3D (DES3D). DES3D is a robust, adaptive, multi-dimensional, finite element method solver which has a composite Elasto-Visco-Plastic rheology. We vary the temperature profile, the ratio of basalt vs. peridotite, the rheology of the mantle peridotites and the loading rate of the subduction interface. Over multiple earthquake cycles, our two end member experiments show that megathrust earthquakes are dominate the seismicity for cold condition (e.g. Japan trench) while both coseismic and aseismic events account for the seismicity for warm condition (e.g. Nankai trench).

Major and micro seismo-volcanic crises in the Asal Rift, Djibouti

NASA Astrophysics Data System (ADS)

Peltzer, G.; Doubre, C.; Tomic, J.

2009-05-01

The Asal-Ghoubbet Rift is located on the eastern branch of the Afar triple junction between the Arabia, Somalia, and Nubia tectonic plates. The last major seismo-volcanic crisis on this segment occurred in November 1978, involving two earthquakes of mb=5+, a basaltic fissure eruption, the development of many open fissures across the rift and up to 80 cm of vertical slip on the bordering faults. Geodetic leveling revealed ~2 m of horizontal opening of the rift accompanied by ~70 cm of subsidence of the inner-floor, consistent with models of the elastic deformation produced by the injection of magma in a system of two dykes. InSAR data acquired at 24-day intervals during the last 12 years by the Canadian Radarsat satellite over the Asal Rift show that the two main faults activated in 1978 continue to slip with periods of steady creep at rates of 0.3-1.3 mm/yr, interrupted by sudden slip events of a few millimeters, in 2000 and 2003. Slip events are coincident with bursts of micro earthquakes distributed around and over the Fieale volcanic center in the eastern part of the Asal Rift. In both cases (the 1978 crisis and micro-slip events), the observed geodetic moment released by fault slip exceeds by a few orders of magnitude the total seismic moment released by earthquakes over the same period. Aseismic fault slip is likely to be the faults response to a changing stress field associated with a volcanic process and not due to dry friction on faults. Sustained injection of magma (1978 crisis) and/or crustal fluids (micro-slip events) in dykes and fissures is a plausible mechanism to control fluid pressure in the basal parts of faults and trigger aseismic slip. In this respect, the micro-events observed by InSAR during a 12-year period of low activity in the rift and the 1978 seismo-volcanic episode are of same nature.

Pressure Wave Propagation along the Décollement of the Nankai Accretionary Wedge: Implications for Aseismic Slip Events

NASA Astrophysics Data System (ADS)

Joshi, A.; Appold, M. S.

2015-12-01

Seismic and hydrologic observations of the Nankai subduction zone made by the Ocean Drilling Program suggest that pore fluid pressures within the accretionary wedge décollement are highly overpressured to near lithostatic values below depths of 2 km beneath the sea floor as a result of sediment diagenesis and dehydration of the subducting oceanic plate. This overpressured zone is also observed to discharge pulses of high fluid pressure that migrate up-dip along the décollement at rates of 1's of km/day. These high pressure pulses along the décollement may cause large enough reductions in the local effective stress to account for aseismic slip events that have been found to propagate also at rates of 1's of km/day. Because elevated fluid pressure and correspondingly decreased effective stress can lead to a dilation of porosity, the pressure waves may become effective agents of fluid transport that can travel more quickly than fluids flowing in the background Darcian flow regime. The purpose of the present study was to seek theoretical confirmation that pressure waves are able to travel quickly enough to account for the seismic and hydrological observations documented. This confirmation was sought through a transient one-dimensional numerical solution to the differential fluid mass conservation equation for an elastic porous medium. Results of the numerical simulations show that when overpressures at depths greater than 2 km in the décollement exceed lithostatic pressure by at least 3%, pressure waves are formed that migrate up-dip at rates fast enough to account for aseismic slip over a broad range of geologic conditions. Pressure waves spawned from these depths in the décollement may travel fast enough to account for aseismic slip when overpressures there are as low as 99% of lithostatic pressure, but require low specific storage of 3×10-6 m-1, high sensitivity of permeability to effective stress, low permeability no higher than about 10-21 m2 at depths below 2 km in the décollement, and an accurate accounting of the decrease in fluid viscosity with increasing depth. Thus, pressure waves could account for aseismic slip in the Nankai accretionary wedge if conditions were near the limits of geologically reasonable ranges.

From Geodetic Imaging of Seismic and Aseismic Fault Slip to Dynamic Modeling of the Seismic Cycle

NASA Astrophysics Data System (ADS)

Avouac, Jean-Philippe

2015-05-01

Understanding the partitioning of seismic and aseismic fault slip is central to seismotectonics as it ultimately determines the seismic potential of faults. Thanks to advances in tectonic geodesy, it is now possible to develop kinematic models of the spatiotemporal evolution of slip over the seismic cycle and to determine the budget of seismic and aseismic slip. Studies of subduction zones and continental faults have shown that aseismic creep is common and sometimes prevalent within the seismogenic depth range. Interseismic coupling is generally observed to be spatially heterogeneous, defining locked patches of stress accumulation, to be released in future earthquakes or aseismic transients, surrounded by creeping areas. Clay-rich tectonites, high temperature, and elevated pore-fluid pressure seem to be key factors promoting aseismic creep. The generally logarithmic time evolution of afterslip is a distinctive feature of creeping faults that suggests a logarithmic dependency of fault friction on slip rate, as observed in laboratory friction experiments. Most faults can be considered to be paved with interlaced patches where the friction law is either rate-strengthening, inhibiting seismic rupture propagation, or rate-weakening, allowing for earthquake nucleation. The rate-weakening patches act as asperities on which stress builds up in the interseismic period; they might rupture collectively in a variety of ways. The pattern of interseismic coupling can help constrain the return period of the maximum- magnitude earthquake based on the requirement that seismic and aseismic slip sum to match long-term slip. Dynamic models of the seismic cycle based on this conceptual model can be tuned to reproduce geodetic and seismological observations. The promise and pitfalls of using such models to assess seismic hazard are discussed.

Introduction to special section on phenomenology, underlying processes, and hazard implications of aseismic slip and nonvolcanic tremor

USGS Publications Warehouse

Gomberg, Joan

2010-01-01

This paper introduces the special section on the "phenomenology, underlying processes, and hazard implications of aseismic slip and nonvolcanic tremor" by highlighting key results of the studies published in it. Many of the results indicate that seismic and aseismic manifestations of slow slip reflect transient shear displacements on the plate interface, with the outstanding exception of northern Cascadia where tremor sources have been located on and above the plate interface (differing models of the plate interface there also need to be reconciled). Slow slip phenomena appear to result from propagating deformation that may develop with persistent gaps and segment boundaries. Results add to evidence that when tectonic deformation is relaxed via slow slip, most relaxation occurs aseismically but with seismic signals providing higher-resolution proxies for the aseismic slip. Instead of two distinct slip modes as suggested previously, lines between "fast" and "slow" slip more appropriately may be described as blurry zones. Results reported also show that slow slip sources do not coincide with a specific temperature or metamorphic reaction. Their associations with zones of high conductivity and low shear to compressional wave velocity ratios corroborate source models involving pore fluid pressure buildup and release. These models and spatial anticorrelations between earthquake and tremor activity also corroborate a linkage between slow slip and frictional properties transitional between steady state and stick-slip. Finally, this special section highlights the benefits of global and multidisciplinary studies, which demonstrate that slow phenomena are not confined to beneath the locked zone but exist in many settings.

Reflection signature of seismic and aseismic slip on the northern Cascadia subduction interface.

PubMed

Nedimović, Mladen R; Hyndman, Roy D; Ramachandran, Kumar; Spence, George D

2003-07-24

At the northern Cascadia margin, the Juan de Fuca plate is underthrusting North America at about 45 mm x yr(-1) (ref. 1), resulting in the potential for destructive great earthquakes. The downdip extent of coupling between the two plates is difficult to determine because the most recent such earthquake (thought to have been in 1700) occurred before instrumental recording. Thermal and deformation studies indicate that, off southern Vancouver Island, the interplate interface is presently fully locked for a distance of approximately 60 km downdip from the deformation front. Great thrust earthquakes on this section of the interface (with magnitudes of up to 9) have been estimated to occur at an average interval of about 590 yr (ref. 3). Further downdip there is a transition from fully locked behaviour to aseismic sliding (where high temperatures allow ductile deformation), with the deep aseismic zone exhibiting slow-slip thrust events. Here we show that there is a change in the reflection character on seismic images from a thin sharp reflection where the subduction thrust is inferred to be locked, to a broad reflection band at greater depth where aseismic slip is thought to be occurring. This change in reflection character may provide a new technique to map the landward extent of rupture in great earthquakes and improve the characterization of seismic hazards in subduction zones.

Seismic wave triggering of nonvolcanic tremor, episodic tremor and slip, and earthquakes on Vancouver Island

NASA Astrophysics Data System (ADS)

Rubinstein, Justin L.; Gomberg, Joan; Vidale, John E.; Wech, Aaron G.; Kao, Honn; Creager, Kenneth C.; Rogers, Garry

2009-02-01

We explore the physical conditions that enable triggering of nonvolcanic tremor and earthquakes by considering local seismic activity on Vancouver Island, British Columbia during and immediately after the arrival of large-amplitude seismic waves from 30 teleseismic and 17 regional or local earthquakes. We identify tremor triggered by four of the teleseismic earthquakes. The close temporal and spatial proximity of triggered tremor to ambient tremor and aseismic slip indicates that when a fault is close to or undergoing failure, it is particularly susceptible to triggering of further events. The amplitude of the triggering waves also influences the likelihood of triggering both tremor and earthquakes such that large amplitude waves triggered tremor in the absence of detectable aseismic slip or ambient tremor. Tremor and energy radiated from regional/local earthquakes share the same frequency passband so that tremor cannot be identified during these smaller, more frequent events. We confidently identify triggered local earthquakes following only one teleseism, that with the largest amplitude, and four regional or local events that generated vigorous aftershock sequences in their immediate vicinity. Earthquakes tend to be triggered in regions different from tremor and with high ambient seismicity rates. We also note an interesting possible correlation between large teleseismic events and episodic tremor and slip (ETS) episodes, whereby ETS events that are "late" and have built up more stress than normal are susceptible to triggering by the slight nudge of the shaking from a large, distant event, while ETS events that are "early" or "on time" are not.

Migrating slow slip detected by slow and repeating earthquakes along the Nankai trough, Japan

NASA Astrophysics Data System (ADS)

Uchida, N.; Obara, K.; Takagi, R.; Asano, Y.

2017-12-01

In the western part of the Nankai trough region, successive occurrences of deep non-volcanic tremors and shallow very low frequency earthquakes (VLFEs) associated with long-term slow slip events (SSEs) are reported in 2003 and 2010. To understand the link between the two seismic slow earthquakes, we identify small repeating earthquake in and around the region from the waveform similarity of earthquakes observed by NIED Hi-net. The result shows the repeaters are located in 15-30 km depth that is in between the depth range of the shallow VLFEs (depth <=15 km) and deep SSEs (depth>= 25km). They are also located outside of the source area of the 1946 Mw8.3 Nankai earthquake, consistent with the hypothesis that repeaters occur due to stress concentration to a locked patch by aseismic slip (creep) in the surrounding area. The long-term trend of aseismic slip estimated from the repeaters shows that the slip rate were faster during 2-3 years period before the 2003 and 2010 episodes. We also found short-term (days to month) accelerations of aseismic slip during the episode of 2010 that migrated toward north. The migration detected from repeaters follows shallow migration of VLFEs and precedes the deep migration of tremors. Therefore we consider that during the period of the long-term SSE of 3 years period, short-term slow slip migrated about 300 km length in 1 month from shallower and south part to deeper and north part of the plate boundary near the edge of the slip area of the Nankai earthquake. Such long-distance migration probably related to large-scale locking of plate boundary that is responsible to the Nankai earthquake and the interseismic stress concentration to the locked area.

Aftershock Comparisons of the Tehuantepec and Puebla Earthquakes: Implications for the Transition between Aseismic and Seismic Behavior?

NASA Astrophysics Data System (ADS)

Richardson, E.

2017-12-01

Reduced aftershock productivity has been observed in subduction zones where slow slip events and aseismic transients have also been observed. A comparison of the aftershock productivity of the recent Tehuantepec and Puebla earthquakes corroborates such observations. The Tehuantepec earthquake of 8 September 2017 produced hundreds of aftershocks and arguably still continues to produce them as of late October 2017, whereas the Puebla earthquake of 19 September 2017 notably lacks aftershocks. This difference in productivity cannot simply be ascribed to differences in mainshock magnitude or detection thresholds. The Puebla earthquake occurred downdip from and just adjacent to the eastern edges of previously observed slow slip events in the Guerrero Gap, whereas the Tehuantepec event is quite removed along strike from the Guerrero Gap and ruptured a patch of fault adjacent to other previous ruptures that also produced standard aftershock sequences. In order to compare aftershock productivity of earthquakes near the Guerrero Gap slow slip region with adjacent regions I used the Advanced National Seismic System catalog and counted aftershocks within a 14-day 100-km window of 42 M>=6.0 slab earthquakes that occurred since 2001 in a box bounded by 13°N and 20°N, and between 91°W and 103°W. This box includes the Guerrero Gap and significant portions of the plate boundary on either side. Preliminary results indicate that ordinary fast-rupturing earthquake productivity in general is much reduced near the location of known SSEs and aftershock productivity of those events that do occur is low compared to earthquakes outside that zone. Earthquakes with low aftershock productivity may represent transitional behavior from aseismic to seismic and in terms of frictional rheology, may represent the transition from velocity weakening to velocity strengthening.

An investigation of deformation and fluid flow at subduction zones using newly developed instrumentation and finite element modeling

NASA Astrophysics Data System (ADS)

Labonte, Alison Louise

Detecting seafloor deformation events in the offshore convergent margin environment is of particular importance considering the significant seismic hazard at subduction zones. Efforts to gain insight into the earthquake cycle have been made at the Cascadia and Costa Rica subduction margins through recent expansions of onshore GPS and seismic networks. While these studies have given scientists the ability to quantify and locate slip events in the seismogenic zone, there is little technology available for adequately measuring offshore aseismic slip. This dissertation introduces an improved flow meter for detecting seismic and aseismic deformation in submarine environments. The value of such hydrologic measurements for quantifying the geodetics at offshore margins is verified through a finite element modeling (FEM) study in which the character of deformation in the shallow subduction zone is determined from previously recorded hydrologic events at the Costa Rica Pacific margin. Accurately sensing aseismic events is one key to determining the stress state in subduction zones as these slow-slip events act to load or unload the seismogenic zone during the interseismic period. One method for detecting seismic and aseismic strain events is to monitor the hydrogeologic response to strain events using fluid flow meters. Previous instrumentation, the Chemical Aqueous Transport (CAT) meter which measures flow rates through the sediment-water interface, can detect transient events at very low flowrates, down to 0.0001 m/yr. The CAT meter performs well in low flow rate environments and can capture gradual changes in flow rate, as might be expected during ultra slow slip events. However, it cannot accurately quantify high flow rates through fractures and conduits, nor does it have the temporal resolution and accuracy required for detecting transient flow events associated with rapid deformation. The Optical Tracer Injection System (OTIS) developed for this purpose is an electronic flow meter that can measure flow rates of 0.1 to >500 m/yr at a temporal resolution of 30 minutes to 0.5 minutes, respectively. Test deployments of the OTIS at cold seeps in the transpressional Monterey Bay demonstrated the OTIS functionality over this range of flow environments. Although no deformation events were detected during these test deployments, the OTIS's temporally accurate measurements at the vigorously flowing Monterey Bay cold seep rendered valuable insight into the plumbing of the seep system. In addition to the capability to detect transient flow events, a primary functional requirement of the OTIS was the ability to communicate and transfer data for long-term real-time monitoring deployments. Real-time data transfer from the OTIS to the desktop was successful during a test deployment of the Nootka Observatory, an acoustically-linked moored-buoy system. A small array of CAT meters was also deployed at the Nootka transform-Cascadia subduction zone triple junction. Four anomalous flow rate events were observed across all four meters during the yearlong deployment. Although the records have low temporal accuracy, a preliminary explanation for the regional changes in flow rate is made through comparison between flow rate records and seismic records. The flow events are thought to be a result of a tectonic deformation event, possibly with an aseismic component. Further constraints are not feasible given the unknown structure of faulting near the triple junction. In a final proof of concept study, I find that use these hydrologic instruments, which capture unique aseismic flow rate patterns, is a valuable method for extracting information about deformation events on the decollement in the offshore subduction zone margin. Transient flow events observed in the frontal prism during a 1999--2000 deployment of CAT meters on the Costa Rica Pacific margin suggest episodic slow-slip deformation events may be occurring in the shallow subduction zone. The FEM study to infer the character of the hypothetical deformation event driving flow transients verify that indeed, a shallow slow-slip event can reproduce the unique flow rate patterns observed. Along (trench) strike variability in the rupture initiation location, and bidirectional propagation, is one way to explain the opposite sign of flow rate transients observed at different along-strike distances. The larger question stimulated by this dissertation project, is: What are the controls on fault mechanics in offshore subduction zone environments? It appears the shallow subduction zone plate interface doesn't behave solely in response to frictional properties of the sediment lining the decollement. Shallow episodic slip at the Costa Rica Pacific margin and further north off Nicaragua, where a slow earthquake broke through the shallow 'stable-sliding' zone and resulted in a tsunami, are potentially conceived through the normally faulted incoming basement topography. Scientists should seek to map out the controls of faulting mechanics, whatever they may be, at all temporal and spatial scales in order to understand these dynamic subduction zone systems. The quest to understanding these controls, in part, requires the characterization of aseismic and seismic strain occurring over time and space. The techniques presented in this dissertation advance scientists' capability for quantifying such strains. With the new instrumentation presented here, long-term real-time observatory networks on the seafloor, and modeling for characterization of deformation events, the pieces of the subduction zone earthquake cycle puzzle may start to come together.

Foreshock triggering of the 1 April 2014 Mw 8.2 Iquique, Chile, earthquake

NASA Astrophysics Data System (ADS)

Herman, Matthew W.; Furlong, Kevin P.; Hayes, Gavin P.; Benz, Harley M.

2016-08-01

On April 1st, 2014, a Mw 8.2 (U.S. Geological Survey moment magnitude) earthquake occurred in the subduction zone offshore northern Chile. In the two weeks leading up to the earthquake, a sequence of foreshocks, starting with a Mw 6.7 earthquake on March 16th and including three more Mw 6.0+ events, occurred predominantly south of the April 1st mainshock epicenter and up-dip of the area of significant slip during the mainshock. Using earthquake locations and source parameters derived in a previous study (Hayes et al., 2014) and a Coulomb failure stress change analysis of these events, we assess in detail the hypothesis that the earthquakes occurred as a cascading sequence, each event successively triggering the next, ultimately triggering the rupture of the mainshock. Following the initial Mw 6.7 event, each of the three largest foreshocks (Mw 6.4, 6.2 and 6.3), as well as the hypocenter of the mainshock, occurred in a region of positive Coulomb stress change produced by the preceding events, indicating these events were brought closer to failure by the prior seismicity. In addition, we reexamine the possibility that aseismic slip occurred and what role it may have played in loading the plate boundary. Using horizontal GPS displacements from along the northern Chile coast prior to the mainshock, we find that the foreshock seismicity alone likely does not account for the observed signals. We perform a grid search for the location and magnitude of an aseismic slip patch that can account for the difference between observed signals and foreshock-related displacement, and find that a slow slip region with slip corresponding to a Mw ∼ 6.8 earthquake located coincident with or up-dip of the foreshock seismicity can best explain this discrepancy. Additionally, such a slow slip region positively loads the mainshock hypocentral area, enhancing the positive loading produced by the foreshock seismicity.

Automated identification and modeling aseismic slip events on Kilauea Volcano, Hawaii

NASA Astrophysics Data System (ADS)

Desmarais, E. K.; Segall, P.; Miklius, A.

2006-12-01

Several aseismic slip events have been observed on the south flank of Kilauea volcano, Hawaii (Cervelli et al., Nature, 2002; Brooks et al., EPSL, 2006; Segall et al., Nature, 2006). These events are identified as spatially coherent offsets in GPS time series. We have interpreted the events as slip on a sub-horizontal surface at depths consistent with a decollement under Kilauea's south flank. In order to determine whether smaller slow slip events are present in the time series, we developed an algorithm that searches for coherent displacement patterns similar to the known slow slip events. We compute candidate displacements by taking a running difference of the mean position 6 days before and after a window of 6 days centered on the candidate time step. The candidate displacements are placed in a 3N dimensional data vector, where N is the number of stations. We then compute the angle, in the 3N dimensional data space, between the candidate displacement and a reference vector at each time step. The reference vector is a stack of displacements due to the four largest known slow slip events. Small angles indicate similar displacement patterns, regardless of amplitude. The algorithm strongly identifies four events (September 20, 1998, November 9, 2000, December 16, 2002, and January 26, 2005), each separated by approximately 2.11 years. The algorithm also identified one smaller event (March 3, 1998) that preceeded the September 1998 event by ~ 200 days, and another event (July 4, 2003) that followed the December 2002 event by ~ 200 days. These smaller, 'paired' events appear to alternate rupturing of the eastern and western parts of the south flank. Each of the slow slip events is correlated with an increase, sometimes slight, in microseismicity on the south flank of Kilauea. The temporal evolution of the microseismicity for the 2005 event is well explained by increased stress due to the slow slip (Segall et al., Nature, 2006). The microearthquakes, at depths of 6.5 - 8.5 km, thus constrain the depth of the slow earthquakes to comparable depths. The triggering of microearthquakes implies that there is finite probability that a larger earthquake could be triggered, given appropriate stress conditions. In order to better constrain the locations of the slow slip events based solely on geodetic observations, we expand on the simple uniform slip models by adding the effects of distributed slip, layered elastic structure, and topography. There are many difficulties in observing slow slip events on Kilauea volcano. The GPS network only provides displacements on land, which is primarily to the north of the largest slip. The vertical displacement field is essential to understanding the northward extent of the slip, however, the GPS observations of slow slip events are primarily observed in the horizontal component, which have smaller noise levels (~ 3 mm). The maximum vertical deformation from the largest event (2005) was very small (± 9 mm), about the same size as the typical vertical noise. We are exploring the possibility that tiltmeters will allow sufficiently accurate measurements to help identify the northern extent of the slip surface.

Creep and slip: Seismic precursors to the Nuugaatsiaq landslide (Greenland)

NASA Astrophysics Data System (ADS)

Poli, Piero

2017-09-01

Precursory signals to material's failure are predicted by numerical models and observed in laboratory experiments or using field data. These precursory signals are a marker of slip acceleration on weak regions, such as crustal faults. Observation of these precursory signals of catastrophic natural events, such as earthquakes and landslides, is necessary for improving our knowledge about the physics of the nucleation process. Furthermore, observing such precursory signals may help to forecast these catastrophic events or reduce their hazard. I report here the observation of seismic precursors to the Nuugaatsiaq landslide in Greenland. Time evolution of the detected precursors implies that an aseismic slip event is taking place for hours before the landslide, with an exponential increase of slip velocity. Furthermore, time evolution of the precursory signals' amplitude sheds light on the evolution of the fault physics during the nucleation process.

Dynamic Triggering of Seismic Events and Their Relation to Slow Slip in Interior Alaska

NASA Astrophysics Data System (ADS)

Sims, N. E.; Holtkamp, S. G.

2017-12-01

We conduct a search for dynamically triggered events in the Minto Flats Fault Zone (MFFZ), a left-lateral strike-slip zone expressed as multiple, partially overlapping faults, in central Alaska. We focus on the MFFZ because we have observed slow slip processes (earthquake swarms and Very Low Frequency Earthquakes) and interaction between earthquake swarms and larger main-shock (MS) events in this area before. We utilize the Alaska Earthquake Center catalog to identify potential earthquake swarms and dynamically triggered foreshock and mainshock events along the fault zone. We find 30 swarms occurring in the last two decades, five of which we classify as foreshock (FS) swarms due to their close proximity in both time and space to MS events. Many of the earthquake swarms cluster around 15-20 km depth, which is near the seismic-aseismic transition along this fault zone. Additionally, we observe instances of large teleseismic events such as the M8.6 2012 Sumatra earthquake and M7.4 2012 Guatemala earthquake triggering seismic events within the MFFZ, with the Sumatra earthquake triggering a mainshock event that was preceded by an ongoing earthquake swarm and the Guatemala event triggering earthquake swarms that subsequently transition into a larger mainshock event. In both cases an earthquake swarm transitioned into a mainshock-aftershock event and activity continued for several days after the teleseismic waves had passed, lending some evidence to delayed dynamic triggering of seismic events. We hypothesize that large dynamic transient strain associated with the passage of teleseismic surface waves is triggering slow slip processes near the base of the seismogenic zone. These triggered aseismic transient events result in earthquake swarms, which sometimes lead to the nucleation of larger earthquakes. We utilize network matched filtering to build more robust catalogs of swarm earthquake families in this region to search for additional swarm-like or triggered activity in response to teleseismic surface waves, and to test dynamic triggering hypotheses.

Revealing the cluster of slow transients behind a large slow slip event.

PubMed

Frank, William B; Rousset, Baptiste; Lasserre, Cécile; Campillo, Michel

2018-05-01

Capable of reaching similar magnitudes to large megathrust earthquakes [ M w (moment magnitude) > 7], slow slip events play a major role in accommodating tectonic motion on plate boundaries through predominantly aseismic rupture. We demonstrate here that large slow slip events are a cluster of short-duration slow transients. Using a dense catalog of low-frequency earthquakes as a guide, we investigate the M w 7.5 slow slip event that occurred in 2006 along the subduction interface 40 km beneath Guerrero, Mexico. We show that while the long-period surface displacement, as recorded by Global Positioning System, suggests a 6-month duration, the motion in the direction of tectonic release only sporadically occurs over 55 days, and its surface signature is attenuated by rapid relocking of the plate interface. Our proposed description of slow slip as a cluster of slow transients forces us to re-evaluate our understanding of the physics and scaling of slow earthquakes.

Paleogeodetic records of seismic and aseismic subduction from central Sumatran microatolls, Indonesia

USGS Publications Warehouse

Natawidjaja, D.H.; Sieh, K.; Ward, S.N.; Cheng, H.; Edwards, R. Lawrence; Galetzka, J.; Suwargadi, B.W.

2004-01-01

We utilize coral microatolls in western Sumatra to document vertical deformation associated with subduction. Microatolls are very sensitive to fluctuations in sea level and thus act as natural tide gauges. They record not only the magnitude of vertical deformation associated with earthquakes (paleoseismic data), but also continuously track the long-term aseismic deformation that occurs during the intervals between earthquakes (paleogeodetic data). This paper focuses on the twentieth century paleogeodetic history of the equatorial region. Our coral paleogeodetic record of the 1935 event reveals a classical example of deformations produced by seismic rupture of a shallow subduction interface. The site closest to the trench rose 90 cm, whereas sites further east sank by as much as 35 cm. Our model reproduces these paleogeodetic data with a 2.3 m slip event on the interface 88 to 125 km from the trench axis. Our coral paleogeodetic data reveal slow submergence during the decades before and after the event in the areas of coseismic emergence. Likewise, interseismic emergence occurred before and after the 1935 event in areas of coseismic submergence. Among the interesting phenomenon we have discovered in the coral record is evidence of a large aseismic slip or "silent even" in 1962, 27 years after the 1935 event. Paleogeodetic deformation rates in the decades before, after, and between the 1935 and 1962 events have varied both temporally and spatially. During the 25 years following the 1935 event, submergence rates were dramatically greater than in prior decades. During the past four decades, however, rates have been lower than in the preceding decades, but are still higher than they were prior to 1935. These paleogeodetic records enable us to model the kinematics of the subduction interface throughout the twentieth century. Copyright 2004 by the American Geophysical Union.

Elastic stress transfer as a diffusive process due to aseismic fault slip in response to fluid injection

NASA Astrophysics Data System (ADS)

Viesca, R. C.

2015-12-01

Subsurface fluid injection is often followed by observations of an enlarging cloud of microseismicity. The cloud's diffusive growth is thought to be a direct response to the diffusion of elevated pore fluid pressure reaching pre-stressed faults, triggering small instabilities; the observed high rates of this growth are interpreted to reflect a relatively high permeability of a fractured subsurface [e.g., Shapiro, GJI 1997]. We investigate an alternative mechanism for growing a microseismic cloud: the elastic transfer of stress due to slow, aseismic slip on a subset of the pre-existing faults in this damaged subsurface. We show that the growth of the slipping region of the fault may be self-similar in a diffusive manner. While this slip is driven by fluid injection, we show that, for critically stressed faults, the apparent diffusion of this slow slip may quickly exceed the poroelastically driven diffusion of the elevated pore fluid pressure. Under these conditions, microseismicity can be first triggered by the off-fault stress perturbation due to the expanding region of slip on principal faults. This provides an alternative interpretation of diffusive growth rates in terms of the subsurface stress state rather than an enhanced hydraulic diffusivity. That such aseismic slip may occur, outpace fluid diffusion, and in turn trigger microseismic events, is also suggested by on- and near-fault observations in past and recently reported fluid injection experiments [e.g., Cornet et al., PAGEOPH 1997; Guglielmi et al., Science 2015]. The model of injection-induced slip assumes elastic off-fault behavior and a fault strength determined by the product of a constant friction coefficient and the local effective normal stress. The sliding region is enlarged by the pore pressure increase resolved on the fault plane. Remarkably, the rate of self-similar expansion may be determined by a single parameter reflecting both the initial stress state and the magnitude of the pore pressure increase.

Seismic and geodetic signatures of fault slip at the Slumgullion Landslide Natural Laboratory

USGS Publications Warehouse

Gomberg, J.; Schulz, W.; Bodin, P.; Kean, J.

2011-01-01

We tested the hypothesis that the Slumgullion landslide is a useful natural laboratory for observing fault slip, specifically that slip along its basal surface and side-bounding strike-slip faults occurs with comparable richness of aseismic and seismic modes as along crustal- and plate-scale boundaries. Our study provides new constraints on models governing landslide motion. We monitored landslide deformation with temporary deployments of a 29-element prism array surveyed by a robotic theodolite and an 88-station seismic network that complemented permanent extensometers and environmental instrumentation. Aseismic deformation observations show that large blocks of the landslide move steadily at approximately centimeters per day, possibly punctuated by variations of a few millimeters, while localized transient slip episodes of blocks less than a few tens of meters across occur frequently. We recorded a rich variety of seismic signals, nearly all of which originated outside the monitoring network boundaries or from the side-bounding strike-slip faults. The landslide basal surface beneath our seismic network likely slipped almost completely aseismically. Our results provide independent corroboration of previous inferences that dilatant strengthening along sections of the side-bounding strike-slip faults controls the overall landslide motion, acting as seismically radiating brakes that limit acceleration of the aseismically slipping basal surface. Dilatant strengthening has also been invoked in recent models of transient slip and tremor sources along crustal- and plate-scale faults suggesting that the landslide may indeed be a useful natural laboratory for testing predictions of specific mechanisms that control fault slip at all scales.

Revealing the cascade of slow transients behind a large slow slip event

NASA Astrophysics Data System (ADS)

Frank, W.; Rousset, B.; Lasserre, C.; Campillo, M.

2017-12-01

Capable of reaching similar magnitudes to large megathrust earthquakes (Mw > 7), slow slip events play a major role in accommodating tectonic motion on plate boundaries. These slip transients are the slow release of built-up tectonic stress that are geodetically imaged as a predominantly aseismic rupture, which is smooth in both time and space. We demonstrate here that large slow slip events are in fact a complex cascade of short-duration slow transients. Using a dense catalog of low-frequency earthquakes as a guide, we investigate the Mw 7.5 slow slip event that occurred in 2006 along the subduction interface 40 km beneath Guerrero, Mexico. We show that while the long-period surface displacement as recorded by GPS suggests a six month duration, motion in the direction of tectonic release only sporadically occurs over 55 days and its surface signature is attenuated by rapid relocking of the plate interface. These results demonstrate that our current conceptual model of slow and continuous rupture is an artifact of low-resolution geodetic observations of a superposition of small, clustered slip events. Our proposed description of slow slip as a cascade of slow transients implies that we systematically overestimate the duration T and underestimate the moment magnitude M of large slow slip events.

Inter-plate aseismic slip on the subducting plate boundaries estimated from repeating earthquakes

NASA Astrophysics Data System (ADS)

Igarashi, T.

2015-12-01

Sequences of repeating earthquakes are caused by repeating slips of small patches surrounded by aseismic slip areas at plate boundary zones. Recently, they have been detected in many regions. In this study, I detected repeating earthquakes which occurred in Japan and the world by using seismograms observed in the Japanese seismic network, and investigated the space-time characteristics of inter-plate aseismic slip on the subducting plate boundaries. To extract repeating earthquakes, I calculate cross-correlation coefficients of band-pass filtering seismograms at each station following Igarashi [2010]. I used two data-set based on USGS catalog for about 25 years from May 1990 and JMA catalog for about 13 years from January 2002. As a result, I found many sequences of repeating earthquakes in the subducting plate boundaries of the Andaman-Sumatra-Java and Japan-Kuril-Kamchatka-Aleutian subduction zones. By applying the scaling relations among a seismic moment, recurrence interval and slip proposed by Nadeau and Johnson [1998], they indicate the space-time changes of inter-plate aseismic slips. Pairs of repeating earthquakes with the longest time interval occurred in the Solomon Islands area and the recurrence interval was about 18.5 years. The estimated slip-rate is about 46 mm/year, which correspond to about half of the relative plate motion in this area. Several sequences with fast slip-rates correspond to the post-seismic slips after the 2004 Sumatra-Andaman earthquake (M9.0), the 2006 Kuril earthquake (M8.3), the 2007 southern Sumatra earthquake (M8.5), and the 2011 Tohoku-oki earthquake (M9.0). The database of global repeating earthquakes enables the comparison of the inter-plate aseismic slips of various plate boundary zones of the world. I believe that I am likely to detect more sequences by extending analysis periods in the area where they were not found in this analysis.

Dynamics of Rifting in two Active Rift Segments in Afar - Geodetic and Structural Studies - DoRA Project

NASA Astrophysics Data System (ADS)

Doubre, C.; Socquet, A.; Masson, F.; Jacques, E.; Grandin, R.; Nercessian, A.; Kassim, M.; Vergne, J.; Diament, M.; Hinderer, J.; Ayele, A.; Lewi, E.; Calais, E.; Peltzer, G.; Toussaint, R.; de Chaballier, J.; Ballu, V. S.; Luck, B.; King, G. C.; Vigny, C.; Cattin, R.; Tiberi, C.; Kidane, T.; Jalludin, M.; Maggi, A.; Dorbath, C.; Manatschal, G.; Schmittbuhl, J.; Le Moigne, N.; Deroussi, S.

2009-12-01

The DoRA project aims to conduct complementary studies in two volcano-tectonic rifts in the Afar Depression. In Northern Afar, the Wal’is Dabbahu Rift (WD, Ethiopia) is currently undergoing a major rifting episode. This event started in September 2005 with a significant seismic activity. InSAR data revealed the injection of a 65 km-long mega-dyke that opened by up to 8 m, the slip of numerous normal faults and opening of fissures, and a rhyolitic eruption. Similarly, the Asal-Ghoubbet Rift (AG, Djibouti) was affected in 1978 by a smaller episode of rifting associated with the intrusion of a 2 m wide dyke into the crust. Since then, a large catalog of geodetic data that includes recent InSAR time series reveals the importance of non-steady deformation controlling the rift dynamics. Our goal is to gain an understanding of such volcano-tectonic segments on several time scales, including the dyking period itself and the post-event period. The study of the behavior of the AG Rift during its whole post-rifting period offers an image at t+30 years of the WD segment, while keeping in mind important structural and scale differences. First, we propose to build a complete and accurate set of geodetic data (InSAR, cGPS, GPS), covering the period under study. With a narrow temporal sample window, we will precisely describe the aseismic slip affecting the normal faults of these rifts, the periods of sudden slip and/or slip acceleration but also measure the deformation associated with probable future dyke intrusion. Second, we aim to constrain the origin of these displacements and their relation with mass transfers within the crust. Series of gravity measurements will be pursue or initiated in both rifts. Third, the recording of seismic activity is essential to constrain the relative importance of seismic and aseismic deformation. This will also help to evaluate the thickness of the seismogenic layer. Together with structural data collected during a seismic survey in the AG Rift, these results will offer crucial constraints on modeling the rifting dynamics in order to test the relative influences of the rheology, the fault/dyke geometry and fluids on the rupture mechanics, the viscous relaxation, dyke intrusion/inflation and aseismic slip and their interactions. Our multidisciplinary approach should provide important new constraints on the dynamics of rifting along divergent plate boundaries, and ultimately, in other geodynamical contexts affected by aseismic fault slip transients.

Aseismic Slip Events along the Southern San Andreas Fault System Captured by Radar Interferometry

DOE Office of Scientific and Technical Information (OSTI.GOV)

Vincent, P

2001-10-01

A seismic slip is observed along several faults in the Salton Sea and southernmost Landers rupture zone regions using interferometric synthetic aperture radar (InSAR) data spanning different time periods between 1992 and 1997. In the southernmost Landers rupture zone, projecting south from the Pinto Mountain Fault, sharp discontinuities in the interferometric phase are observed along the sub-parallel Burnt Mountain and Eureka Peak Faults beginning three months after the Landers earthquake and is interpreted to be post-Landers after-slip. Abrupt phase offsets are also seen along the two southernmost contiguous 11 km Durmid Hill and North Shore segments of the San Andreasmore » Fault with an abrupt termination of slip near the northern end of the North Shore Segment. A sharp phase offset is seen across 20 km of the 30 km-long Superstition Hills Fault before phase decorrelation in the Imperial Valley along the southern 10 km of the fault prevents coherent imaging by InSAR. A time series of deformation interferograms suggest most of this slip occurred between 1993 and 1995 and none of it occurred between 1992 and 1993. A phase offset is also seen along a 5 km central segment of the Coyote Creek fault that forms a wedge with an adjoining northeast-southwest trending conjugate fault. Most of the slip observed on the southern San Andreas and Superstition Hills Faults occurred between 1993 and 1995--no slip is observed in the 92-93 interferograms. These slip events, especially the Burnt Mountain and Eureka Peak events, are inferred to be related to stress redistribution from the June, 1992 M{sub w} = 7.3 Landers earthquake. Best-fit elastic models of the San Andreas and Superstition Hills slip events suggest source mechanisms with seismic moments over three orders of magnitude larger than a maximum possible summation of seismic moments from all seismicity along each fault segment during the entire 4.8-year time interval spanned by the InSAR data. Aseismic moment releases of this magnitude (equivalent to M{sub w} = 5.3 and 5.6 events on the Superstition Hills and San Andreas Faults respectively) are hitherto unknown and have not been captured previously by any geodetic technique.« less

Foreshocks during the nucleation of stick-slip instability

USGS Publications Warehouse

McLaskey, Gregory C.; Kilgore, Brian D.

2013-01-01

We report on laboratory experiments which investigate interactions between aseismic slip, stress changes, and seismicity on a critically stressed fault during the nucleation of stick-slip instability. We monitor quasi-static and dynamic changes in local shear stress and fault slip with arrays of gages deployed along a simulated strike-slip fault (2 m long and 0.4 m deep) in a saw cut sample of Sierra White granite. With 14 piezoelectric sensors, we simultaneously monitor seismic signals produced during the nucleation phase and subsequent dynamic rupture. We observe localized aseismic fault slip in an approximately meter-sized zone in the center of the fault, while the ends of the fault remain locked. Clusters of high-frequency foreshocks (Mw ~ −6.5 to −5.0) can occur in this slowly slipping zone 5–50 ms prior to the initiation of dynamic rupture; their occurrence appears to be dependent on the rate at which local shear stress is applied to the fault. The meter-sized nucleation zone is generally consistent with theoretical estimates, but source radii of the foreshocks (2 to 70 mm) are 1 to 2 orders of magnitude smaller than the theoretical minimum length scale over which earthquake nucleation can occur. We propose that frictional stability and the transition between seismic and aseismic slip are modulated by local stressing rate and that fault sections, which would typically slip aseismically, may radiate seismic waves if they are rapidly stressed. Fault behavior of this type may provide physical insight into the mechanics of foreshocks, tremor, repeating earthquake sequences, and a minimum earthquake source dimension.

A silent slip event on the deeper Cascadia subduction interface.

PubMed

Dragert, G; Wang, K; James, T S

2001-05-25

Continuous Global Positioning System sites in southwestern British Columbia, Canada, and northwestern Washington state, USA, have been moving landward as a result of the locked state of the Cascadia subduction fault offshore. In the summer of 1999, a cluster of seven sites briefly reversed their direction of motion. No seismicity was associated with this event. The sudden displacements are best explained by approximately 2 centimeters of aseismic slip over a 50-kilometer-by-300-kilometer area on the subduction interface downdip from the seismogenic zone, a rupture equivalent to an earthquake of moment magnitude 6.7. This provides evidence that slip of the hotter, plastic part of the subduction interface, and hence stress loading of the megathrust earthquake zone, can occur in discrete pulses.

Sudden aseismic fault slip on the south flank of Kilauea volcano.

PubMed

Cervelli, Peter; Segall, Paul; Johnson, Kaj; Lisowski, Michael; Miklius, Asta

2002-02-28

One of the greatest hazards associated with oceanic volcanoes is not volcanic in nature, but lies with the potential for catastrophic flank failure. Such flank failure can result in devastating tsunamis and threaten not only the immediate vicinity, but coastal cities along the entire rim of an ocean basin. Kilauea volcano on the island of Hawaii, USA, is a potential source of such flank failures and has therefore been monitored by a network of continuously recording geodetic instruments, including global positioning system (GPS) receivers, tilt meters and strain meters. Here we report that, in early November 2000, this network recorded transient southeastward displacements, which we interpret as an episode of aseismic fault slip. The duration of the event was about 36 hours, it had an equivalent moment magnitude of 5.7 and a maximum slip velocity of about 6[?]cm per day. Inversion of the GPS data reveals a shallow-dipping thrust fault at a depth of 4.5[?]km that we interpret as the down-dip extension of the Hilina Pali--Holei Pali normal fault system. This demonstrates that continuously recording geodetic networks can detect accelerating slip, potentially leading to warnings of volcanic flank collapse.

Coseismic and aseismic deformations of the rock mass around deep level mining in South Africa - Joint South African and Japanese study

NASA Astrophysics Data System (ADS)

Milev, A. M.; Yabe, Y.; Naoi, M. M.; Nakatani, M.; Durrheim, R. J.; Ogasawara, H.; Scholz, C. H.

2010-12-01

Two underground sites in a deep level gold mine in South Africa were instrumented by the Council for Scientific and Industrial Research (CSIR) with tilt meters and seismic monitors. One of the sites was also instrumented by JApanese-German Underground Acoustic emission Research in South Africa (JAGUARS) with a small network, approx. 40 m span, of eight Acoustic Emission (AE) sensors. The rate of tilt, defined as quasi-static deformations, and the seismic ground motion, defined as dynamic deformations, were analysed in order to understand the rock mass behavior around deep level mining. In addition the high frequency AE events recorded at hypocentral distances of about 50m were analysed. This was the first implementation of high frequency AE events at such a great depth (3300m below the surface). A good correspondence between the dynamic and quasi-static deformations was found. The rate of coseismic and aseismic tilt, as well as seismicity recorded by the mine seismic network, are approximately constant until the daily blasting time, which takes place from about 19:30 until shortly before 21:00. During the blasting time and the subsequent seismic events the coseismic and aseismic tilt shows a rapid increase indicated by a rapid change of the tilt during the seismic event. Much of the quasi-static deformation, however, occurs independently of the seismic events and was described as ‘slow’ or aseismic events. During the monitoring period a seismic event with MW 1.9 (2.1) occurred in the vicinity of the instrumented site. This event was recorded by both the CSIR integrated monitoring system and JAGUARS acoustic emotion network. The tilt changes associated with this event showed a well pronounced after-tilt. More than 21,000 AE aftershocks were located in the first 150 hours after the main event. Using the distribution of the AE events the position of the fault in the source area was successfully delineated. The distribution of the AE events following the main shock was related to after tilt in order to quantify post slip behavior of the source. There was no evidence found for coseismic expansion of the source after the main slip. An attempt to associate the different type of deformations with the various fracture regions and geological structures around the stopes was carried out. A model, was introduced in which the coseismic deformations are associated with the stress regime outside the stope fracture envelope and very often located on existing geological structures, while the aseismic deformations are associated with mobilization of fractures and stress relaxation within the fracture envelope.

Recurrent slow slip events as a barrier to the northward rupture propagation of the 2016 Pedernales earthquake (Central Ecuador)

NASA Astrophysics Data System (ADS)

Vaca, Sandro; Vallée, Martin; Nocquet, Jean-Mathieu; Battaglia, Jean; Régnier, Marc

2018-01-01

The northern Ecuador segment of the Nazca/South America subduction zone shows spatially heterogeneous interseismic coupling. Two highly coupled zones (0.4° S-0.35° N and 0.8° N-4.0° N) are separated by a low coupled area, hereafter referred to as the Punta Galera-Mompiche Zone (PGMZ). Large interplate earthquakes repeatedly occurred within the coupled zones in 1958 (Mw 7.7) and 1979 (Mw 8.1) for the northern patch and in 1942 (Mw 7.8) and 2016 (Mw 7.8) for the southern patch, while the whole segment is thought to have rupture during the 1906 Mw 8.4-8.8 great earthquake. We find that during the last decade, the PGMZ has experienced regular and frequent seismic swarms. For the best documented sequence (December 2013-January 2014), a joint seismological and geodetic analysis reveals a six-week-long Slow Slip Event (SSE) associated with a seismic swarm. During this period, the microseismicity is organized into families of similar earthquakes spatially and temporally correlated with the evolution of the aseismic slip. The moment release (3.4 × 1018 Nm, Mw 6.3), over a 60 × 40 km area, is considerably larger than the moment released by earthquakes (5.8 × 1015 Nm, Mw 4.4) during the same time period. In 2007-2008, a similar seismic-aseismic episode occurred, with higher magnitudes both for the seismic and aseismic processes. Cross-correlation analyses of the seismic waveforms over a 15 years-long period further suggest a 2-year repeat time for seismic swarms, which also implies that SSEs recurrently affect this area. Such SSEs contribute to release the accumulated stress, likely explaining why the 2016 Pedernales earthquake did not propagate northward into the PGMZ.

Slip rate and tremor genesis in Cascadia

USGS Publications Warehouse

Wech, Aaron G.; Bartlow, Noel M.

2014-01-01

At many plate boundaries, conditions in the transition zone between seismogenic and stable slip produce slow earthquakes. In the Cascadia subduction zone, these events are consistently observed as slow, aseismic slip on the plate interface accompanied by persistent tectonic tremor. However, not all slow slip at other plate boundaries coincides spatially and temporally with tremor, leaving the physics of tremor genesis poorly understood. Here we analyze seismic, geodetic, and strainmeter data in Cascadia to observe for the first time a large, tremor-generating slow earthquake change from tremor-genic to silent and back again. The tremor falls silent at reduced slip speeds when the migrating slip front pauses as it loads the stronger adjacent fault segment to failure. The finding suggests that rheology and slip-speed-regulated stressing rate control tremor genesis, and the same section of fault can slip both with and without detectable tremor, limiting tremor's use as a proxy for slip.

Explaining postseismic and aseismic transient deformation in subduction zones with rate and state friction modeling constrained by lab and geodetic observations

NASA Astrophysics Data System (ADS)

Liu, Y.; Dedontney, N. L.; Rice, J. R.

2007-12-01

Rate and state friction, as applied to modeling subduction earthquake sequences, routinely predicts postseismic slip. It also predicts spontaneous aseismic slip transients, at least when pore pressure p is highly elevated near and downdip from the stability transition [Liu and Rice, 2007]. Here we address how to make such postseismic and transient predictions more fully compatible with geophysical observations. For example, lab observations can determine the a, b parameters and state evolution slip L of rate and state friction as functions of lithology and temperature and, with aid of a structural and thermal model of the subduction zone, as functions of downdip distance. Geodetic observations constrain interseismic, postseismic and aseismic transient deformations, which are controlled in the modeling by the distributions of a \\barσ and b \\barσ (parameters which also partly control the seismic rupture phase), where \\barσ = σ - p. Elevated p, controlled by tectonic compression and dehydration, may be constrained by petrologic and seismic observations. The amount of deformation and downdip extent of the slipping zone associated with the spontaneous quasi- periodic transients, as thus far modeled [Liu and Rice, 2007], is generally smaller than that observed during episodes of slow slip events in northern Cascadia and SW Japan subduction zones. However, the modeling was based on lab data for granite gouge under hydrothermal conditions because data is most complete for that case. We here report modeling based on lab data on dry granite gouge [Stesky, 1975; Lockner et al., 1986], involving no or lessened chemical interaction with water and hence being a possibly closer analog to dehydrated oceanic crust, and limited data on gabbro gouge [He et al., 2007], an expected lithology. Both data sets show a much less rapid increase of a-b with temperature above the stability transition (~ 350 °C) than does wet granite gouge; a-b increases to ~ 0.08 for wet granite at 600 °C, but to only ~ 0.01 in the dry granite and gabbro cases. We find that the lessened high-T a - b does, for the same \\barσ, modestly extend the transient slip episodes further downdip, although a majority of slip is still contributed near and in the updip rate-weakening region. However, postseismic slip, for the same \\barσ, propagates much further downdip into the rate-strengthening region. To better constrain the downdip distribution of (a - b) \\barσ, and possibly a \\barσ and L, we focus on the geodetically constrained [Hutton et al., 2001] space-time distribution of postseismic slip for the 1995 Mw = 8.0 Colima-Jalisco earthquake. This is a similarly shallow dipping subduction zone with a thermal profile [Currie et al., 2001] comparable to those that have thus far been shown to exhibit aseismic transients and non-volcanic tremor [Peacock et al., 2002]. We extrapolate the modeled 2-D postseismic slip, following a thrust earthquake with a coseismic slip similar to the 1995 event, to a spatial-temporal 3-D distribution. Surface deformation due to such slips on the thrust fault in an elastic half space is calculated and compared to that observed at western Mexico GPS stations, to constrain the above depth-variable model parameters.

Distribution of aseismic slip rate on the Hayward fault inferred from seismic and geodetic data

USGS Publications Warehouse

Schmidt, D.A.; Burgmann, R.; Nadeau, R.M.; d'Alessio, M.

2005-01-01

We solve for the slip rate distribution on the Hayward fault by performing a least squares inversion,of geodetic and seismic data sets. Our analysis focuses on the northern 60 km of the fault. Interferometric synthetic aperture radar (InSAR) data from 13 independent ERS interferograms are stacked to obtain range change rates from 1992 to 2000. Horizontal surface displacement rates at 141 bench marks are measured using GPS from 1994 to 2003. Surface creep observations and estimates of deep slip rates determined from characteristic repeating earthquake sequences are also incorporated in the inversion. The fault is discretized into 283 triangular dislocation elements that approximate the nonplanar attributes of the fault surface. South of the city of Hayward, a steeply, east dipping fault geometry accommodates the divergence of the surface trace and the microseismicity at depth. The inferred slip rate distribution is consistent with a fault that creeps aseismically at a rate of ???5 mm/yr to a depth of 4-6 km. The interferometric synthetic aperture radar (InSAR) data require an aseismic slip rate that approaches the geologic slip rate on the northernmost fault segment beneath Point Pinole, although the InSAR data might be complicated by a small dip-slip component at this location. A low slip rate patch of <1 mm/yr is inferred beneath San Leandro consistent with the source location of the 1868 earthquake. We calculate that the entire fault is accumulating a slip rate deficit equivalent to a Mw = 6.77 ?? 0.05 per century. However, this estimate of potential coseismic moment represents an upper bound because we do not know how much of the accumulated strain will be released through aseismic processes such as afterslip. Copyright 2005 by the American Geophysical Union.

Pore Pressure Evolution in Shallow Subduction Earthquake Sequences and Effects on Aseismic Slip Transients -- Numerical Modeling With Rate and State Friction

NASA Astrophysics Data System (ADS)

Liu, Y.; Rice, J. R.

2005-12-01

In 3D modeling of long tectonic loading and earthquake sequences on a shallow subduction fault [Liu and Rice, 2005], with depth-variable rate and state friction properties, we found that aseismic transient slip episodes emerge spontaneously with only a simplified representation of effects of metamorphic fluid release. That involved assumption of a constant in time but uniformly low effective normal stress in the downdip region. As suggested by observations in several major subduction zones [Obara, 2002; Rogers and Dragert, 2003; Kodaira et al, 2004], the presence of fluids, possibly released from dehydration reactions beneath the seismogenic zone, and their pressurization within the fault zone may play an important role in causing aseismic transients and associated non-volcanic tremors. To investigate the effects of fluids in the subduction zone, particularly on the generation of aseismic transients and their various features, we develop a more complete physical description of the pore pressure evolution (specifically, pore pressure increase due to supply from dehydration reactions and shear heating, decrease due to transport and dilatancy during slip), and incorporate that into the rate and state based 3D modeling. We first incorporated two important factors, dilatancy and shear heating, following Segall and Rice [1995, 2004] and Taylor [1998]. In the 2D simulations (slip varies with depth only), a dilatancy-stabilizing effect is seen which slows down the seismic rupture front and can prevent rapid slip from extending all the way to the trench, similarly to Taylor [1998]. Shear heating increases the pore pressure, and results in faster coseismic rupture propagation and larger final slips. In the 3D simulations, dilatancy also stabilizes the along-strike rupture propagation of both seismic and aseismic slips. That is, aseismic slip transients migrate along the strike faster with a shorter Tp (the characteristic time for pore pressure in the fault core to re-equilibrate with that of its surroundings). This is consistent with our previous simulations, which show that the aseismic transients migrate along the strike at a higher speed under a lower, constant in time, effective normal stress. As a combination of the two factors, we show the pore pressure evolution with drops (due to dilatancy during slip) and then rises (due to shear heating) on the fault over multiple time scales. We next plan to formulate, and merge with the slip-rupture analysis, fuller fluid release models based on phase equilibria and models of transport in which the average fault-parallel permeability is a decreasing function of the effective normal stress. The thrust fault zone, at seismogenic depths and slightly downdip, is represented in a conceptually similar manner to the well-studied major continental faults, assuming the fault core materials have a lower permeability than the neighboring damaged zone. Heat diffusion in the fault core and damaged zone will also be considered in the modeling. The simulation results may help to improve our understanding of the processes of the aseismic transients observed within a transform plate boundary along the SAF near Cholame, California [Nadeau and Dolenc, 2005].

Episodic Tremor and Slip Explained by Fluid-Enhanced Microfracturing and Sealing

NASA Astrophysics Data System (ADS)

Bernaudin, M.; Gueydan, F.

2018-04-01

Episodic tremor and slow-slip events at the deep extension of plate boundary faults illuminate seismic to aseismic processes around the brittle-ductile transition. These events occur in volumes characterized by overpressurized fluids and by near failure shear stress conditions. We present a new modeling approach based on a ductile grain size-sensitive rheology with microfracturing and sealing, which provides a mechanical and field-based explanation of such phenomena. We also model pore fluid pressure variation as a function of changes in porosity/permeability and strain rate-dependent fluid pumping. The fluid-enhanced dynamic evolution of microstructures defines cycles of ductile strain localization and implies increase in pore fluid pressure. We propose that slow-slip events are ductile processes related to transient strain localization, while nonvolcanic tremor corresponds to fracturing of the whole rock at the peak of pore fluid pressure. Our model shows that the availability of fluids and the efficiency of fluid pumping control the occurrence and the P-T conditions of episodic tremor and slip.

Interpretation of Microseismicity Observed From Surface and Borehole Seismic Arrays During Hydraulic Fracturing in Shale - Bedding Plane Slip Model

NASA Astrophysics Data System (ADS)

Stanek, F.; Jechumtalova, Z.; Eisner, L.

2017-12-01

We present a geomechanical model explaining microseismicity induced by hydraulic fracturing in shales developed from many datasets acquired with two most common types of seismic monitoring arrays, surface and dual-borehole arrays. The geomechanical model explains the observed source mechanisms and locations of induced events from two stimulated shale reservoirs. We observe shear dip-slip source mechanisms with nodal planes aligned with location trends. We show that such seismicity can be explained as a shearing along bedding planes caused by aseismic opening of vertical hydraulic fractures. The source mechanism inversion was applied only to selected high-quality events with sufficient signal-to-noise ratio. We inverted P- and P- and S-wave arrival amplitudes to full-moment tensor and decomposed it to shear, volumetric and compensated linear vector dipole components. We also tested an effect of noise presented in the data to evaluate reliability of non-shear components. The observed seismicity from both surface and downhole monitoring of shale stimulations is very similar. The locations of induced microseismic events are limited to narrow depth intervals and propagate along distinct trend(s) showing fracture propagation in direction of maximum horizontal stress from injection well(s). The source mechanisms have a small non-shear component which can be partly explained as an effect of noise in the data, i.e. events represent shearing on faults. We observe predominantly dip-slip events with a strike of the steeper (almost vertical) nodal plane parallel to the fracture propagation. Therefore the other possible nodal plane is almost horizontal. The rake angles of the observed mechanisms divide these dip-slips into two groups with opposite polarities. It means that we observe opposite movements on the nearly identically oriented faults. Realizing a typical structural weakness of shale in horizontal planes, we interpret observed microseismicity as a result of shearing along bedding planes caused by seismically silent (aseismic) vertical fracture opening.

Simulating spontaneous aseismic and seismic slip events on evolving faults

NASA Astrophysics Data System (ADS)

Herrendörfer, Robert; van Dinther, Ylona; Pranger, Casper; Gerya, Taras

2017-04-01

Plate motion along tectonic boundaries is accommodated by different slip modes: steady creep, seismic slip and slow slip transients. Due to mainly indirect observations and difficulties to scale results from laboratory experiments to nature, it remains enigmatic which fault conditions favour certain slip modes. Therefore, we are developing a numerical modelling approach that is capable of simulating different slip modes together with the long-term fault evolution in a large-scale tectonic setting. We extend the 2D, continuum mechanics-based, visco-elasto-plastic thermo-mechanical model that was designed to simulate slip transients in large-scale geodynamic simulations (van Dinther et al., JGR, 2013). We improve the numerical approach to accurately treat the non-linear problem of plasticity (see also EGU 2017 abstract by Pranger et al.). To resolve a wide slip rate spectrum on evolving faults, we develop an invariant reformulation of the conventional rate-and-state dependent friction (RSF) and adapt the time step (Lapusta et al., JGR, 2000). A crucial part of this development is a conceptual ductile fault zone model that relates slip rates along discrete planes to the effective macroscopic plastic strain rates in the continuum. We test our implementation first in a simple 2D setup with a single fault zone that has a predefined initial thickness. Results show that deformation localizes in case of steady creep and for very slow slip transients to a bell-shaped strain rate profile across the fault zone, which suggests that a length scale across the fault zone may exist. This continuum length scale would overcome the common mesh-dependency in plasticity simulations and question the conventional treatment of aseismic slip on infinitely thin fault zones. We test the introduction of a diffusion term (similar to the damage description in Lyakhovsky et al., JMPS, 2011) into the state evolution equation and its effect on (de-)localization during faster slip events. We compare the slip spectrum in our simulations to conventional RSF simulations (Liu and Rice, JGR, 2007). We further demonstrate the capability of simulating the evolution of a fault zone and simultaneous occurrence of slip transients. From small random initial distributions of the state variable in an otherwise homogeneous medium, deformation localizes and forms curved zones of reduced states. These spontaneously formed fault zones host slip transients, which in turn contribute to the growth of the fault zone.

Reconsidering earthquake scaling

USGS Publications Warehouse

Gomberg, Joan S.; Wech, Aaron G.; Creager, Kenneth; Obara, K.; Agnew, Duncan

2016-01-01

The relationship (scaling) between scalar moment, M0, and duration, T, potentially provides key constraints on the physics governing fault slip. The prevailing interpretation of M0-T observations proposes different scaling for fast (earthquakes) and slow (mostly aseismic) slip populations and thus fundamentally different driving mechanisms. We show that a single model of slip events within bounded slip zones may explain nearly all fast and slow slip M0-T observations, and both slip populations have a change in scaling, where the slip area growth changes from 2-D when too small to sense the boundaries to 1-D when large enough to be bounded. We present new fast and slow slip M0-T observations that sample the change in scaling in each population, which are consistent with our interpretation. We suggest that a continuous but bimodal distribution of slip modes exists and M0-T observations alone may not imply a fundamental difference between fast and slow slip.

Aseismic Transform Fault Slip at the Mendocino Triple Junction From Characteristically Repeating Earthquakes

NASA Astrophysics Data System (ADS)

Materna, Kathryn; Taira, Taka'aki; Bürgmann, Roland

2018-01-01

The Mendocino Triple Junction (MTJ), at the northern terminus of the San Andreas Fault system, is an actively deforming plate boundary region with poorly constrained estimates of seismic coupling on most offshore fault surfaces. Characteristically repeating earthquakes provide spatial and temporal descriptions of aseismic creep at the MTJ, including on the oceanic transform Mendocino Fault Zone (MFZ) as it subducts beneath North America. Using a dataset of earthquakes from 2008 to 2017, we find that the easternmost segment of the MFZ displays creep during this period at about 65% of the long-term slip rate. We also find creep at slower rates on the shallower strike-slip interface between the Pacific plate and the North American accretionary wedge, as well as on a fault that accommodates Gorda subplate internal deformation. After a nearby M5.7 earthquake in 2015, we observe a possible decrease in aseismic slip on the near-shore MFZ that lasts from 2015 to at least early 2017.

Multiscale Dynamics of Aseismic Slip on Central San Andreas Fault

NASA Astrophysics Data System (ADS)

Khoshmanesh, M.; Shirzaei, M.

2018-03-01

Understanding the evolution of aseismic slip enables constraining the fault's seismic budget and provides insight into dynamics of creep. Inverting the time series of surface deformation measured along the Central San Andreas Fault obtained from interferometric synthetic aperture radar in combination with measurements of repeating earthquakes, we constrain the spatiotemporal distribution of creep during 1992-2010. We identify a new class of intermediate-term creep rate variations that evolve over decadal scale, releasing stress on the accelerating zone and loading adjacent decelerating patches. We further show that in short-term (<2 year period), creep avalanches, that is, isolated clusters of accelerated aseismic slip with velocities exceeding the long-term rate, govern the dynamics of creep. The statistical properties of these avalanches suggest existence of elevated pore pressure in the fault zone, consistent with laboratory experiments.

Change-point detection of induced and natural seismicity

NASA Astrophysics Data System (ADS)

Fiedler, B.; Holschneider, M.; Zoeller, G.; Hainzl, S.

2016-12-01

Earthquake rates are influenced by tectonic stress buildup, earthquake-induced stress changes, and transient aseismic sources. While the first two sources can be well modeled due to the fact that the source is known, transient aseismic processes are more difficult to detect. However, the detection of the associated changes of the earthquake activity is of great interest, because it might help to identify natural aseismic deformation patterns (such as slow slip events) and the occurrence of induced seismicity related to human activities. We develop a Bayesian approach to detect change-points in seismicity data which are modeled by Poisson processes. By means of a Likelihood-Ratio-Test, we proof the significance of the change of the intensity. The model is also extended to spatiotemporal data to detect the area of the transient changes. The method is firstly tested for synthetic data and then applied to observational data from central US and the Bardarbunga volcano in Iceland.

Slow slip events and the 2016 Te Araroa Mw 7.1 earthquake interaction: Northern Hikurangi subduction, New Zealand

NASA Astrophysics Data System (ADS)

Koulali, A.; McClusky, S.; Wallace, L.; Allgeyer, S.; Tregoning, P.; D'Anastasio, E.; Benavente, R.

2017-08-01

Following a sequence of three Slow Slip Events (SSEs) on the northern Hikurangi Margin, between June 2015 and August 2016, a Mw 7.1 earthquake struck 30 km offshore of the East Cape region in the North Island of New Zealand on the 2 September 2016 (NZ local time). The earthquake was also followed by a transient deformation event (SSE or afterslip) northeast of the North Island, closer to the earthquake source area. We use data from New Zealand's continuous Global Positioning System networks to invert for the SSE slip distribution and evolution on the Hikurangi subduction interface. Our slip inversion results show an increasing amplitude of the slow slip toward the Te Araroa earthquake foreshock and main shock area, suggesting a possible triggering of the Mw 7.1 earthquake by the later stage of the slow slip sequence. We also show that the transient deformation following the Te Araroa earthquake ruptured a portion of the Hikurangi Trench northeast of the North Island, farther north than any previously observed Hikurangi margin SSEs. Our slip inversion and the coulomb stress calculation suggest that this transient may have been induced as a response to the increase in the static coulomb stress change downdip of the rupture plane on the megathrust. These observations show the importance of considering the interaction between slow slip events, seismic, and aseismic events, not only on the same megathrust interface but also on faults within the surrounding crust.

Aseismic Deformation Associated with an Earthquake Swarm in the Northern Apennines (Italy)

NASA Astrophysics Data System (ADS)

Gualandi, A.; Nichele, C.; Serpelloni, E.; Chiaraluce, L.; Anderlini, L.; Latorre, D.; Belardinelli, M. E.; Avouac, J. P.

2017-12-01

Analyzing the displacement time series from continuous GPS (cGPS) with an Independent Component Analysis (ICA) we detect a transient deformation signal that correlates both in space and time with a seismic swarm activity (maximum Mw = 3.69 ± 0.09) occurred in the hanging wall of the Altotiberina normal fault (Northern Apennines, Italy) in 2013-2014. The geodetic transient lasted ˜6 months and produced a NW-SE trending extension of ˜ 5.3 mm, consistent with the regional tectonic regime. The seismicity and the geodetic signal are consistent with slip on two splay faults in the ATF hanging wall. Comparing the seismic moment associated with the geodetic transient and the seismic events, we observe that seismicity accounts for only a fraction of the measured geodetic deformation. The combined seismic and aseismic slip decreased the Coulomb stress on the locked shallow portion of the ATF, while the transition region to the creeping section has been loaded.

The Slip Behavior and Source Parameters for Spontaneous Slip Events on Rough Faults Subjected to Slow Tectonic Loading

NASA Astrophysics Data System (ADS)

Tal, Yuval; Hager, Bradford H.

2018-02-01

We study the response to slow tectonic loading of rough faults governed by velocity weakening rate and state friction, using a 2-D plane strain model. Our numerical approach accounts for all stages in the seismic cycle, and in each simulation we model a sequence of two earthquakes or more. We focus on the global behavior of the faults and find that as the roughness amplitude, br, increases and the minimum wavelength of roughness decreases, there is a transition from seismic slip to aseismic slip, in which the load on the fault is released by more slip events but with lower slip rate, lower seismic moment per unit length, M0,1d, and lower average static stress drop on the fault, Δτt. Even larger decreases with roughness are observed when these source parameters are estimated only for the dynamic stage of the rupture. For br ≤ 0.002, the source parameters M0,1d and Δτt decrease mutually and the relationship between Δτt and the average fault strain is similar to that of a smooth fault. For faults with larger values of br that are completely ruptured during the slip events, the average fault strain generally decreases more rapidly with roughness than Δτt.

Source parameters of a M4.8 and its accompanying repeating earthquakes off Kamaishi, NE Japan: Implications for the hierarchical structure of asperities and earthquake cycle

USGS Publications Warehouse

Uchida, N.; Matsuzawa, T.; Ellsworth, W.L.; Imanishi, K.; Okada, T.; Hasegawa, A.

2007-01-01

We determine the source parameters of a M4.9 ?? 0.1 'characteristic earthquake' sequence and its accompanying microearthquakes at ???50 km depth on the subduction plate boundary offshore of Kamaishi, NE Japan. The microearthquakes tend to occur more frequently in the latter half of the recurrence intervals of the M4.9 ?? 0.1 events. Our results show that the microearthquakes are repeating events and they are located not only around but also within the slip area for the 2001 M4.8 event. From the hierarchical structure of slip areas and smaller stress drops for the microearthquakes compared to the M4.8 event, we infer the small repeating earthquakes rupture relatively weak patches in and around the slip area for the M4.8 event and their activity reflects a stress concentration process and/or change in frictional property (healing) at the area. We also infer the patches for the M4.9 ?? 0.1 and other repeating earthquakes undergo aseismic slip during their interseismic period. Copyright 2007 by the American Geophysical Union.

Alternative source models of very low frequency events

NASA Astrophysics Data System (ADS)

Gomberg, J.; Agnew, D. C.; Schwartz, S. Y.

2016-09-01

We present alternative source models for very low frequency (VLF) events, previously inferred to be radiation from individual slow earthquakes that partly fill the period range between slow slip events lasting thousands of seconds and low-frequency earthquakes (LFE) with durations of tenths of a second. We show that VLF events may emerge from bandpass filtering a sum of clustered, shorter duration, LFE signals, believed to be the components of tectonic tremor. Most published studies show VLF events occurring concurrently with tremor bursts and LFE signals. Our analysis of continuous data from Costa Rica detected VLF events only when tremor was also occurring, which was only 7% of the total time examined. Using analytic and synthetic models, we show that a cluster of LFE signals produces the distinguishing characteristics of VLF events, which may be determined by the cluster envelope. The envelope may be diagnostic of a single, dynamic, slowly slipping event that propagates coherently over kilometers or represents a narrowly band-passed version of nearly simultaneous arrivals of radiation from slip on multiple higher stress drop and/or faster propagating slip patches with dimensions of tens of meters (i.e., LFE sources). Temporally clustered LFE sources may be triggered by single or multiple distinct aseismic slip events or represent the nearly simultaneous chance occurrence of background LFEs. Given the nonuniqueness in possible source durations, we suggest it is premature to draw conclusions about VLF event sources or how they scale.

Alternative source models of very low frequency events

USGS Publications Warehouse

Gomberg, Joan S.; Agnew, D.C.; Schwartz, S.Y.

2016-01-01

We present alternative source models for very low frequency (VLF) events, previously inferred to be radiation from individual slow earthquakes that partly fill the period range between slow slip events lasting thousands of seconds and low-frequency earthquakes (LFE) with durations of tenths of a second. We show that VLF events may emerge from bandpass filtering a sum of clustered, shorter duration, LFE signals, believed to be the components of tectonic tremor. Most published studies show VLF events occurring concurrently with tremor bursts and LFE signals. Our analysis of continuous data from Costa Rica detected VLF events only when tremor was also occurring, which was only 7% of the total time examined. Using analytic and synthetic models, we show that a cluster of LFE signals produces the distinguishing characteristics of VLF events, which may be determined by the cluster envelope. The envelope may be diagnostic of a single, dynamic, slowly slipping event that propagates coherently over kilometers or represents a narrowly band-passed version of nearly simultaneous arrivals of radiation from slip on multiple higher stress drop and/or faster propagating slip patches with dimensions of tens of meters (i.e., LFE sources). Temporally clustered LFE sources may be triggered by single or multiple distinct aseismic slip events or represent the nearly simultaneous chance occurrence of background LFEs. Given the nonuniqueness in possible source durations, we suggest it is premature to draw conclusions about VLF event sources or how they scale.

Recognising Paleoseismic Events and Slip Styles in Vein Microstructures - is Incrementality Enough?

NASA Astrophysics Data System (ADS)

Fagereng, A.; Sibson, R. H.

2008-12-01

'Subduction channels', containing highly sheared, fluid-saturated, trench-fill sediments, are commonly present along subduction thrust interfaces. These shear zones accommodate fast plate boundary slip rates (1~-~10~cm/yr) and exhibit high levels of seismicity, accomplishing slip in a broad range of styles including standard earthquakes, slow slip, non-volcanic tremor and aseismic creep. Exhumed subduction channel fault rocks provide a time-integrated record of these varied slip modes though the degree of overprinting may be considerable. The Chrystalls Beach accretionary mélange, within the Otago Schist accretion-collision assemblage, New Zealand, is analogous to an active subduction channel assemblage. It contains asymmetric lenses of sandstone, chert and minor basalt enclosed within a relatively incompetent, cleaved pelitic matrix. This assemblage has been intensely sheared in a mixed continuous/discontinuous style within a flat-lying, <~4~km thick, shear zone. Ductile structures such as folds, S/C-like structures, and asymmetric boudins and clasts formed by soft sediment deformation and pressure solution creep. An extensive anastomosing vein system can be divided into mutually cross-cutting extension fractures (V1) and slickenfibre shear veins (V2). V1 commonly cut competent lenses within the mélange, while V2 mostly follow lithological contacts. Both vein sets are predominantly elongate-blocky with 'crack-seal' extension and shear increments of 10~- ~100~μm. Little sign of wall rock alteration or heating is present adjacent to V1 veins, which likely formed by incremental hydrofracture with episodic fluid influx. Post-fracture drop in Pf promoted solute precipitation from advecting fluids. This process may reflect fracture and fluid flow in a distributed fault-fracture mesh, an often inferred mechanism of non-volcanic tremor. In contrast, wall rock alteration and pressure solution seams are common adjacent to V2 veins. Slickenfibres on these shear surfaces likely formed by relatively slow dissolution and precipitation of wall rock material, which may translate to a slip mode of rise-time intermediate between earthquakes (seconds - minutes) and aseismic creep (years - infinite). Fibres are typically ≤ 10 cm long, similar to slip observed in slow slip events (rise-time weeks - months). We propose that these veins are possible records of slow slip along weak, fluid-saturated and highly overpressured planes. No definite record of large, fast earthquakes is observed in the complex, either because the rocks never experienced such events, or because significant shear heating was inhibited by thermal pressurisation. The only record of fast events would be discrete planes of cataclasite, easily overprinted by slow interseismic material diffusion. The mélange is a record of episodic, distributed deformation over a range of time- and length-scales, which may reflect distributed seismic activity accommodated by a range of slip modes including episodic tremor and slow slip.

Slow slip events and seismic tremor at circum-Pacific subduction zones

NASA Astrophysics Data System (ADS)

Schwartz, Susan Y.; Rokosky, Juliana M.

2007-09-01

It has been known for a long time that slip accompanying earthquakes accounts for only a fraction of plate tectonic displacements. However, only recently has a fuller spectrum of strain release processes, including normal, slow, and silent earthquakes (or slow slip events) and continuous and episodic slip, been observed and generated by numerical simulations of the earthquake cycle. Despite a profusion of observations and modeling studies the physical mechanism of slow slip events remains elusive. The concurrence of seismic tremor with slow slip episodes in Cascadia and southwestern Japan provides insight into the process of slow slip. A perceived similarity between subduction zone and volcanic tremor has led to suggestions that slow slip involves fluid migration on or near the plate interface. Alternatively, evidence is accumulating to support the notion that tremor results from shear failure during slow slip. Global observations of the location, spatial extent, magnitude, duration, slip rate, and periodicity of these aseismic slip transients indicate significant variation that may be exploited to better understand their generation. Most slow slip events occur just downdip of the seismogenic zone, consistent with rate- and state-dependent frictional modeling that requires unstable to stable transitional properties for slow slip generation. At a few convergent margins the occurrence of slow slip events within the seismogenic zone makes it highly likely that transitions in frictional properties exist there and are the loci of slow slip nucleation. Slow slip events perturb the surrounding stress field and may either increase or relieve stress on a fault, bringing it closer to or farther from earthquake failure, respectively. This paper presents a review of slow slip events and related seismic tremor observed at plate boundaries worldwide, with a focus on circum-Pacific subduction zones. Trends in global observations of slow slip events suggest that (1) slow slip is a common phenomena observed at almost all subduction zones with instrumentation capable of recording it, (2) different frictional properties likely control fast versus slow slip, (3) the depth range of slow slip may be related to the thermal properties of the plate interface, and (4) the equivalent seismic moment of slow slip events is proportional to their duration (Moατ), different from the Moατ3 scaling observed for earthquakes.

Slow Slip and Earthquake Nucleation in Meter-Scale Laboratory Experiments

NASA Astrophysics Data System (ADS)

Mclaskey, G.

2017-12-01

The initiation of dynamic rupture is thought to be preceded by a quasistatic nucleation phase. Observations of recent earthquakes sometimes support this by illuminating slow slip and foreshocks in the vicinity of the eventual hypocenter. I describe laboratory earthquake experiments conducted on two large-scale loading machines at Cornell University that provide insight into the way earthquake nucleation varies with normal stress, healing time, and loading rate. The larger of the two machines accommodates a 3 m long granite sample, and when loaded to 7 MPa stress levels, we observe dynamic rupture events that are preceded by a measureable nucleation zone with dimensions on the order of 1 m. The smaller machine accommodates a 0.76 m sample that is roughly the same size as the nucleation zone. On this machine, small variations in nucleation properties result in measurable differences in slip events, and we generate both dynamic rupture events (> 0.1 m/s slip rates) and slow slip events ( 0.001 to 30 mm/s slip rates). Slow events occur when instability cannot fully nucleate before reaching the sample ends. Dynamic events occur after long healing times or abrupt increases in loading rate which suggests that these factors shrink the spatial and temporal extents of the nucleation zone. Arrays of slip, strain, and ground motion sensors installed on the sample allow us to quantify seismic coupling and study details of premonitory slip and afterslip. The slow slip events we observe are primarily aseismic (less than 1% of the seismic coupling of faster events) and produce swarms of very small M -6 to M -8 events. These mechanical and seismic interactions suggest that faults with transitional behavior—where creep, small earthquakes, and tremor are often observed—could become seismically coupled if loaded rapidly, either by a slow slip front or dynamic rupture of an earthquake that nucleated elsewhere.

Measuring Aseismic Slip through Characteristically Repeating Earthquakes at the Mendocino Triple Junction, Northern California

NASA Astrophysics Data System (ADS)

Materna, K.; Taira, T.; Burgmann, R.

2016-12-01

The Mendocino Triple Junction (MTJ), at the transition point between the San Andreas fault system, the Mendocino Transform Fault, and the Cascadia Subduction Zone, undergoes rapid tectonic deformation and produces more large (M>6.0) earthquakes than any region in California. Most of the active faults of the triple junction are located offshore, making it difficult to characterize both seismic slip and aseismic creep. In this work, we study aseismic creep rates near the MTJ using characteristically repeating earthquakes (CREs) as indicators of creep rate. CREs are generally interpreted as repeated failures of the same seismic patch within an otherwise creeping fault zone; as a consequence, the magnitude and recurrence time of the CREs can be used to determine a fault's creep rate through empirically calibrated scaling relations. Using seismic data from 2010-2016, we identify CREs as recorded by an array of eight 100-Hz PBO borehole seismometers deployed in the Cape Mendocino area. For each event pair with epicenters less than 30 km apart, we compute the cross-spectral coherence of 20 seconds of data starting one second before the P-wave arrival. We then select pairs with high coherence in an appropriate frequency band, which is determined uniquely for each event pair based on event magnitude, station distance, and signal-to-noise ratio. The most similar events (with median coherence above 0.95 at two or more stations) are selected as CREs and then grouped into CRE families, and each family is used to infer a local creep rate. On the Mendocino Transform Fault, we find relatively high creep rates of >5 cm/year that increase closer to the Gorda Ridge. Closer to shore and to the MTJ itself, we find many families of repeaters on and off the transform fault with highly variable creep rates, indicative of the complex deformation that takes place there.

Seismic and aseismic fault slip in response to fluid injection observed during field experiments at meter scale

NASA Astrophysics Data System (ADS)

Cappa, F.; Guglielmi, Y.; De Barros, L.; Wynants-Morel, N.; Duboeuf, L.

2017-12-01

During fluid injection, the observations of an enlarging cloud of seismicity are generally explained by a direct response to the pore pressure diffusion in a permeable fractured rock. However, fluid injection can also induce large aseismic deformations which provide an alternative mechanism for triggering and driving seismicity. Despite the importance of these two mechanisms during fluid injection, there are few studies on the effects of fluid pressure on the partitioning between seismic and aseismic motions under controlled field experiments. Here, we describe in-situ meter-scale experiments measuring synchronously the fluid pressure, the fault motions and the seismicity directly in a fault zone stimulated by controlled fluid injection at 280 m depth in carbonate rocks. The experiments were conducted in a gallery of an underground laboratory in south of France (LSBB, http://lsbb.eu). Thanks to the proximal monitoring at high-frequency, our data show that the fluid overpressure mainly induces a dilatant aseismic slip (several tens of microns up to a millimeter) at the injection. A sparse seismicity (-4 < Mw < -3) is observed several meters away from the injection, in a part of the fault zone where the fluid overpressure is null or very low. Using hydromechanical modeling with friction laws, we simulated an experiment and investigated the relative contribution of the fluid pressure diffusion and stress transfer on the seismic and aseismic fault behavior. The model reproduces the hydromechanical data measured at injection, and show that the aseismic slip induced by fluid injection propagates outside the pressurized zone where accumulated shear stress develops, and potentially triggers seismicity. Our models also show that the permeability enhancement and friction evolution are essential to explain the fault slip behavior. Our experimental results are consistent with large-scale observations of fault motions at geothermal sites (Wei et al., 2015; Cornet, 2016), and suggest that controlled field experiments at meter-scale are important for better assessing the role of fluid pressure in natural and human-induced earthquakes.

Retardations in fault creep rates before local moderate earthquakes along the San Andreas fault system, central California

USGS Publications Warehouse

Burford, R.O.

1988-01-01

Records of shallow aseismic slip (fault creep) obtained along parts of the San Andreas and Calaveras faults in central California demonstrate that significant changes in creep rates often have been associated with local moderate earthquakes. An immediate postearthquake increase followed by gradual, long-term decay back to a previous background rate is generally the most obvious earthquake effect on fault creep. This phenomenon, identified as aseismic afterslip, usually is characterized by above-average creep rates for several months to a few years. In several cases, minor step-like movements, called coseismic slip events, have occurred at or near the times of mainshocks. One extreme case of coseismic slip, recorded at Cienega Winery on the San Andreas fault 17.5 km southeast of San Juan Bautista, consisted of 11 mm of sudden displacement coincident with earthquakes of ML=5.3 and ML=5.2 that occurred 2.5 minutes apart on 9 April 1961. At least one of these shocks originated on the main fault beneath the winery. Creep activity subsequently stopped at the winery for 19 months, then gradually returned to a nearly steady rate slightly below the previous long-term average. The phenomena mentioned above can be explained in terms of simple models consisting of relatively weak material along shallow reaches of the fault responding to changes in load imposed by sudden slip within the underlying seismogenic zone. In addition to coseismic slip and afterslip phenomena, however, pre-earthquake retardations in creep rates also have been observed. Onsets of significant, persistent decreases in creep rates have occurred at several sites 12 months or more before the times of moderate earthquakes. A 44-month retardation before the 1979 ML=5.9 Coyote Lake earthquake on the Calaveras fault was recorded at the Shore Road creepmeter site 10 km northwest of Hollister. Creep retardation on the San Andreas fault near San Juan Bautista has been evident in records from one creepmeter site for the past 5 years. Retardations with durations of 21 and 19 months also occurred at Shore Road before the 1974 and 1984 earthquakes of ML=5.2 and ML=6.2, respectively. Although creep retardation remains poorly understood, several possible explanations have been discussed previously. (1) Certain onsets of apparent creep retardation may be explained as abrupt terminations of afterslip generated from previous moderate-mainshock sequences. (2) Retardations may be related to significant decreases in the rate of seismic and/or aseismic slip occurring within or beneath the underlying seismogenic zone. Such decreases may be caused by changes in local conditions related to growth of asperities, strain hardening, or dilatancy, or perhaps by passage of stress-waves or other fluctuations in driving stresses. (3) Finally, creep rates may be lowered (or increased) by stresses imposed on the fault by seismic or aseismic slip on neighboring faults. In addition to causing creep-rate increases or retardations, such fault interactions occasionally may trigger earthquakes. Regardless of the actual mechanisms involved and the current lack of understanding of creep retardation, it appears that shallow fault creep is sensitive to local and regional effects that promote or accompany intermediate-term preparation stages leading to moderate earthquakes. A strategy for more complete monitoring of fault creep, wherever it is known to occur, therefore should be assigned a higher priority in our continuing efforts to test various hypotheses concerning the mechanical relations between seismic and aseismic slip. ?? 1988 Birkha??user Verlag.

Retardations in fault creep rates before local moderate earthquakes along the San Andreas fault system, central California

NASA Astrophysics Data System (ADS)

Burford, Robert O.

1988-06-01

Records of shallow aseismic slip (fault creep) obtained along parts of the San Andreas and Calaveras faults in central California demonstrate that significant changes in creep rates often have been associated with local moderate earthquakes. An immediate postearthquake increase followed by gradual, long-term decay back to a previous background rate is generally the most obvious earthquake effect on fault creep. This phenomenon, identified as aseismic afterslip, usually is characterized by above-average creep rates for several months to a few years. In several cases, minor step-like movements, called coseismic slip events, have occurred at or near the times of mainshocks. One extreme case of coseismic slip, recorded at Cienega Winery on the San Andreas fault 17.5 km southeast of San Juan Bautista, consisted of 11 mm of sudden displacement coincident with earthquakes of M L =5.3 and M L =5.2 that occurred 2.5 minutes apart on 9 April 1961. At least one of these shocks originated on the main fault beneath the winery. Creep activity subsequently stopped at the winery for 19 months, then gradually returned to a nearly steady rate slightly below the previous long-term average. The phenomena mentioned above can be explained in terms of simple models consisting of relatively weak material along shallow reaches of the fault responding to changes in load imposed by sudden slip within the underlying seismogenic zone. In addition to coseismic slip and afterslip phenomena, however, pre-earthquake retardations in creep rates also have been observed. Onsets of significant, persistent decreases in creep rates have occurred at several sites 12 months or more before the times of moderate earthquakes. A 44-month retardation before the 1979 M L =5.9 Coyote Lake earthquake on the Calaveras fault was recorded at the Shore Road creepmeter site 10 km northwest of Hollister. Creep retardation on the San Andreas fault near San Juan Bautista has been evident in records from one creepmeter site for the past 5 years. Retardations with durations of 21 and 19 months also occurred at Shore Road before the 1974 and 1984 earthquakes of M L =5.2 and M L =6.2, respectively. Although creep retardation remains poorly understood, several possible explanations have been discussed previously. (1) Certain onsets of apparent creep retardation may be explained as abrupt terminations of afterslip generated from previous moderate-mainshock sequences. (2) Retardations may be related to significant decreases in the rate of seismic and/or aseismic slip occurring within or beneath the underlying seismogenic zone. Such decreases may be caused by changes in local conditions related to growth of asperities, strain hardening, or dilatancy, or perhaps by passage of stress-waves or other fluctuations in driving stresses. (3) Finally, creep rates may be lowered (or increased) by stresses imposed on the fault by seismic or aseismic slip on neighboring faults. In addition to causing creep-rate increases or retardations, such fault interactions occasionally may trigger earthquakes. Regardless of the actual mechanisms involved and the current lack of understanding of creep retardation, it appears that shallow fault creep is sensitive to local and regional effects that promote or accompany intermediate-term preparation stages leading to moderate earthquakes. A strategy for more complete monitoring of fault creep, wherever it is known to occur, therefore should be assigned a higher priority in our continuing efforts to test various hypotheses concerning the mechanical relations between seismic and aseismic slip.

A comparison of long-baseline strain data and fault creep records obtained near Hollister, California

USGS Publications Warehouse

Slater, L.E.; Burford, R.O.

1979-01-01

A comparison of creepmeter records from nine sites along a 12-km segment of the Calaveras fault near Hollister, California and long-baseline strain changes for nine lines in the Hollister multiwavelength distance-measuring (MWDM) array has established that episodes of large-scale deformation both preceded and accompanied periods of creep activity monitored along the fault trace during 1976. A concept of episodic, deep-seated aseismic slip that contributes to loading and subsequent aseismic failure of shallow parts of the fault plane seems attractive, implying that the character of aseismic slip sensed along the surface trace may be restricted to a relatively shallow (~ 1-km) region on the fault plane. Preliminary results from simple dislocation models designed to test the concept demonstrate that extending the time-histories and amplitudes of creep events sensed along the fault trace to depths of up to 10 km on the fault plane cannot simulate adequately the character and amplitudes of large-scale episodic movements observed at points more than 1 km from the fault. Properties of a 2-3-km-thick layer of unconsolidated sediments present in Hollister Valley, combined with an essentially rigid-block behavior in buried basement blocks, might be employed in the formulation of more appropriate models that could predict patterns of shallow fault creep and large-scale displacements much more like those actually observed. ?? 1979.

Using Low-Frequency Earthquake Families on the San Andreas Fault as Deep Creepmeters

NASA Astrophysics Data System (ADS)

Thomas, A. M.; Beeler, N. M.; Bletery, Q.; Burgmann, R.; Shelly, D. R.

2018-01-01

The central section of the San Andreas Fault hosts tectonic tremor and low-frequency earthquakes (LFEs) similar to subduction zone environments. LFEs are often interpreted as persistent regions that repeatedly fail during the aseismic shear of the surrounding fault allowing them to be used as creepmeters. We test this idea by using the recurrence intervals of individual LFEs within LFE families to estimate the timing, duration, recurrence interval, slip, and slip rate associated with inferred slow slip events. We formalize the definition of a creepmeter and determine whether this definition is consistent with our observations. We find that episodic families reflect surrounding creep over the interevent time, while the continuous families and the short time scale bursts that occur as part of the episodic families do not. However, when these families are evaluated on time scales longer than the interevent time these events can also be used to meter slip. A straightforward interpretation of episodic families is that they define sections of the fault where slip is distinctly episodic in well-defined slow slip events that slip 16 times the long-term rate. In contrast, the frequent short-term bursts of the continuous and short time scale episodic families likely do not represent individual creep events but rather are persistent asperities that are driven to failure by quasi-continuous creep on the surrounding fault. Finally, we find that the moment-duration scaling of our inferred creep events are inconsistent with the proposed linear moment-duration scaling. However, caution must be exercised when attempting to determine scaling with incomplete knowledge of scale.

Spatial and temporal variations in creep rate along the El Pilar fault at the Caribbean-South American plate boundary (Venezuela), from InSAR

NASA Astrophysics Data System (ADS)

Pousse Beltran, Léa.; Pathier, Erwan; Jouanne, François; Vassallo, Riccardo; Reinoza, Carlos; Audemard, Franck; Doin, Marie Pierre; Volat, Matthieu

2016-11-01

In eastern Venezuela, the Caribbean-South American plate boundary follows the El Pilar fault system. Previous studies based on three GPS campaigns (2003-2005-2013) demonstrated that the El Pilar fault accommodates the whole relative displacement between the two tectonic plates (20 mm/yr) and proposed that 50-60% of the slip is aseismic. In order to quantify the possible variations of the aseismic creep in time and space, we conducted an interferometric synthetic aperture radar (InSAR) time series analysis, using the (NSBAS) New Small BAseline Subset method, on 18 images from the Advanced Land Observing Satellite (ALOS-1) satellite spanning the 2007-2011 period. During this 3.5 year period, InSAR observations show that aseismic slip decreases eastward along the fault: the creep rate of the western segment reaches 25.3 ± 9.4 mm/yr on average, compared to 13.4 ± 6.9 mm/yr on average for the eastern segment. This is interpreted, through slip distribution models, as being related to coupled and uncoupled areas between the surface and 20 km in depth. InSAR observations also show significant temporal creep rate variations (accelerations) during the considered time span along the western segment. The transient behavior of the creep is not consistent with typical postseismic afterslip following the 1997 Ms 6.8 earthquake. The creep is thus interpreted as persistent aseismic slip during an interseismic period, which has a pulse- or transient-like behavior.

A 20-year catalog comparing smooth and sharp estimates of slow slip events in Cascadia

NASA Astrophysics Data System (ADS)

Molitors Bergman, E. G.; Evans, E. L.; Loveless, J. P.

2017-12-01

Slow slip events (SSEs) are a form of aseismic strain release at subduction zones resulting in a temporary reversal in interseismic upper plate motion over a period of weeks, frequently accompanied in time and space by seismic tremor at the Cascadia subduction zone. Locating SSEs spatially along the subduction zone interface is essential to understanding the relationship between SSEs, earthquakes, and tremor and assessing megathrust earthquake hazard. We apply an automated slope comparison-based detection algorithm to single continuously recording GPS stations to determine dates and surface displacement vectors of SSEs, then apply network-based filters to eliminate false detections. The main benefits of this algorithm are its ability to detect SSEs while they are occurring and track the spatial migration of each event. We invert geodetic displacement fields for slip distributions on the subduction zone interface for SSEs between 1997 and 2017 using two estimation techniques: spatial smoothing and total variation regularization (TVR). Smoothing has been frequently used in determining the location of interseismic coupling, earthquake rupture, and SSE slip and yields spatially coherent but inherently blurred solutions. TVR yields compact, sharply bordered slip estimates of similar magnitude and along-strike extent to previously presented studied events, while fitting the constraining geodetic data as well as corresponding smoothing-based solutions. Slip distributions estimated using TVR have up-dip limits that align well with down-dip limits of interseismic coupling on the plate interface and spatial extents that approximately correspond to the distribution of tremor concurrent with each event. TVR gives a unique view of slow slip distributions that can contribute to understanding of the physical properties that govern megathrust slip processes.

DOE Office of Scientific and Technical Information (OSTI.GOV)

McNally, K.C.; Minster, J.B.

Revised estimates of seismic slip rates along the Middle America Trench are lower on the average than plate convergence rates but match them locally (for example, Oaxaca). Along the Cocos-North American plate boundary this can be explained by nonuniformities in slip at points of aseismic ridge or fracture zone subduction. For at least 81 yr (and possibly several hundred years), no major (M/sub s/> or =7.5) shallow earthquake is known to have occurred near the Orozco Fracture Zone and Tehuantepec Ridge areas. Compared with the average recurrence periods for large earthquakes (33 +- 8 yr since 1898 and 35 +-more » 24 yr between 1542 and 1979), this suggests that either a large (M> or =8.4) event may be anticipated at such locations, or that these are points of aseismic subduction. Large coastal terraces and evidence suggesting tectonic uplift are found onshore near the Orozco Fracture zone. The larger discrepancy between plate convergence and seismic slip rates along the Cocos-Carribbean plate boundary is more likely due to decoupling and downbending of the subducted plate. We used the limited statistical evidence available to characterize both spatial and temporal deficiencies in recent seismic slip. The observations appear consistent with a possible forthcoming episode of more intense seismic activity. Based on a series of comparisons with carefully delineated aftershock zones, we conclude that the zones of anomalous seismic activity can be indentified by a systematic, automated analysis of the worldwide earthquake catalog (m/sub b/> or =4).« less

Active accommodation of plate convergence in Southern Iran: Earthquake locations, triggered aseismic slip, and regional strain rates

NASA Astrophysics Data System (ADS)

Barnhart, William D.; Lohman, Rowena B.; Mellors, Robert J.

2013-10-01

We present a catalog of interferometric synthetic aperture radar (InSAR) constraints on deformation that occurred during earthquake sequences in southern Iran between 1992 and 2011, and explore the implications on the accommodation of large-scale continental convergence between Saudi Arabia and Eurasia within the Zagros Mountains. The Zagros Mountains, a salt-laden fold-and-thrust belt involving ~10 km of sedimentary rocks overlying Precambrian basement rocks, have formed as a result of ongoing continental collision since 10-20 Ma that is currently occurring at a rate of ~3 cm/yr. We first demonstrate that there is a biased misfit in earthquake locations in global catalogs that likely results from neglect of 3-D velocity structure. Previous work involving two M ~ 6 earthquakes with well-recorded aftershocks has shown that the deformation observed with InSAR may represent triggered slip on faults much shallower than the primary earthquake, which likely occurred within the basement rocks (>10 km depth). We explore the hypothesis that most of the deformation observed with InSAR spanning earthquake sequences is also due to shallow, triggered slip above a deeper earthquake, effectively doubling the moment release for each event. We quantify the effects that this extra moment release would have on the discrepancy between seismically and geodetically constrained moment rates in the region, finding that even with the extra triggered fault slip, significant aseismic deformation during the interseismic period is necessary to fully explain the convergence between Eurasia and Saudi Arabia.

Salton Trough Post-seismic Afterslip, Viscoelastic Response, and Contribution to Regional Hazard

NASA Astrophysics Data System (ADS)

Parker, J. W.; Donnellan, A.; Lyzenga, G. A.

2012-12-01

The El Mayor-Cucapah M7.2 April 4 2010 earthquake in Baja California may have affected accumulated hazard to Southern California cities due to loading of regional faults including the Elsinore, San Jacinto and southern San Andreas, faults which already have over a century of tectonic loading. We examine changes observed via multiple seismic and geodetic techniques, including micro seismicity and proposed seismicity-based indicators of hazard, high-quality fault models, the Plate Boundary Observatory GNSS array (with 174 stations showing post-seismic transients with greater than 1 mm amplitude), and interferometric radar maps from UAVSAR (aircraft) flights, showing a network of aseismic fault slip events at distances up to 60 km from the end of the surface rupture. Finite element modeling is used to compute the expected coseismic motions at GPS stations with general agreement, including coseismic uplift at sites ~200 km north of the rupture. Postseismic response is also compared, with GNSS and also with the CIG software "RELAX." An initial examination of hazard is made comparing micro seismicity-based metrics, fault models, and changes to coulomb stress on nearby faults using the finite element model. Comparison of seismicity with interferograms and historic earthquakes show aseismic slip occurs on fault segments that have had earthquakes in the last 70 years, while other segments show no slip at the surface but do show high triggered seismicity. UAVSAR-based estimates of fault slip can be incorporated into the finite element model to correct Coloumb stress change.

Spatiotemporal evolution of a transient slip event on the San Andreas fault near Parkfield, California

USGS Publications Warehouse

Murray, J.R.; Segall, P.

2005-01-01

In 1993 several baselines of the two-color electronic distance meter (EDM) network at Parkfield, California, deviated from their long-term rates, coincident with anomalous observations from nearby strain meters and a creep meter, as well as an increase in microseismicity. Between October 1992 and December 1994, three M ??? 4.5 earthquakes occurred beneath Middle Mountain, near the hypocenter of the 1934 and 1966 Parkfield M6 events. We analyzed the two-color EDM data using a Kalman-filtering based technique to image the spatiotemporal evolution of slip on the fault at Parkfield between the mid-1980s and 2003. This method accounts for localized random walk motion of the geodetic monuments and a prominent seasonal signal that affects many baselines. We find that a slip rate increase occurred between January 1993 and July 1996 on the upper 8 km of the fault near Middle Mountain. The peak estimated slip rate during this time was 49 mm/yr, which exceeds the long-term geologic rate of ???35 mm/yr. The slip rate evolution appears episodic, with an initial modest increase after the M4.3 earthquake and a much larger jump following the shallower M4.7 event in December 1994. This temporal correlation between inferred slip and seismicity suggests that the moderate earthquakes triggered the aseismic fault slip. The EDM data cannot resolve whether transient slip propagated across the nucleation zone of the 1934 and 1966 M6 Parkfield earthquakes. However, transient slip and its associated stress release in the hypocentral area of previous Parkfield events is consistent with the nucleation of the 2004 M6 Parkfield earthquake elsewhere on the fault. Copyright 2005 by the American Geophysical Union.

Coseismic and postseismic deformation due to the 2007 M5.5 Ghazaband fault earthquake, Balochistan, Pakistan

NASA Astrophysics Data System (ADS)

Fattahi, H.; Amelung, F.; Chaussard, E.; Wdowinski, S.

2015-05-01

Time series analysis of interferometric synthetic aperture radar data reveals coseismic and postseismic surface displacements associated with the 2007 M5.5 earthquake along the southern Ghazaband fault, a major but little studied fault in Pakistan. Modeling indicates that the coseismic surface deformation was caused by ~9 cm of strike-slip displacement along a shallow subvertical fault. The earthquake was followed by at least 1 year of afterslip, releasing ~70% of the moment of the main event, equivalent to a M5.4 earthquake. This high aseismic relative to the seismic moment release is consistent with previous observations for moderate earthquakes (M < 6) and suggests that smaller earthquakes are associated with a higher aseismic relative to seismic moment release than larger earthquakes.

Bayesian explorations of fault slip evolution over the earthquake cycle

NASA Astrophysics Data System (ADS)

Duputel, Z.; Jolivet, R.; Benoit, A.; Gombert, B.

2017-12-01

The ever-increasing amount of geophysical data continuously opens new perspectives on fundamental aspects of the seismogenic behavior of active faults. In this context, the recent fleet of SAR satellites including Sentinel-1 and COSMO-SkyMED permits the use of InSAR for time-dependent slip modeling with unprecedented resolution in time and space. However, existing time-dependent slip models rely on spatial smoothing regularization schemes, which can produce unrealistically smooth slip distributions. In addition, these models usually do not include uncertainty estimates thereby reducing the utility of such estimates. Here, we develop an entirely new approach to derive probabilistic time-dependent slip models. This Markov-Chain Monte Carlo method involves a series of transitional steps to predict and update posterior Probability Density Functions (PDFs) of slip as a function of time. We assess the viability of our approach using various slow-slip event scenarios. Using a dense set of SAR images, we also use this method to quantify the spatial distribution and temporal evolution of slip along a creeping segment of the North Anatolian Fault. This allows us to track a shallow aseismic slip transient lasting for about a month with a maximum slip of about 2 cm.

Mixed-Mode Slip Behavior of the Altotiberina Low-Angle Normal Fault System (Northern Apennines, Italy) through High-Resolution Earthquake Locations and Repeating Events

NASA Astrophysics Data System (ADS)

Valoroso, Luisa; Chiaraluce, Lauro; Di Stefano, Raffaele; Monachesi, Giancarlo

2017-12-01

We generated a 4.5-year-long (2010-2014) high-resolution earthquake catalogue, composed of 37,000 events with ML < 3.9 and MC = 0.5 completeness magnitude, to report on the seismic activity of the Altotiberina (ATF) low-angle normal fault system and to shed light on the mechanical behavior and seismic potential of this fault, which is capable of generating a M7 event. Seismicity defines the geometry of the fault system composed of the low-angle (15°-20°) ATF, extending for 50 km along strike and between 4 and 16 km at depth showing an 1.5 km thick fault zone made of multiple subparallel slipping planes, and a complex network of synthetic/antithetic higher-angle segments located in the ATF hanging wall (HW) that can be traced along strike for up to 35 km. Ninety percent of the recorded seismicity occurs along the high-angle HW faults during a series of minor, sometimes long-lasting (months) seismic sequences with multiple MW3+ mainshocks. Remaining earthquakes (ML < 2.4) are released instead along the low-angle ATF at a constant rate of 2.2 events per day. Within the ATF-related seismicity, we found 97 clusters of repeating earthquakes (RE), mostly consisting of doublets occurring during short interevent time (hours). RE are located within the geodetically recognized creeping portions of the ATF, around the main locked asperity. The rate of occurrence of RE seems quite synchronous with the ATF-HW seismic release, suggesting that creeping may guide the strain partitioning in the ATF system. The seismic moment released by the ATF seismicity accounts for 30% of the geodetic one, implying aseismic deformation. The ATF-seismicity pattern is thus consistent with a mixed-mode (seismic and aseismic) slip behavior.

Possible deep fault slip preceding the 2004 Parkfield earthquake, inferred from detailed observations of tectonic tremor

USGS Publications Warehouse

Shelly, David R.

2009-01-01

Earthquake predictability depends, in part, on the degree to which sudden slip is preceded by slow aseismic slip. Recently, observations of deep tremor have enabled inferences of deep slow slip even when detection by other means is not possible, but these data are limited to certain areas and mostly the last decade. The region near Parkfield, California, provides a unique convergence of several years of high-quality tremor data bracketing a moderate earthquake, the 2004 magnitude 6.0 event. Here, I present detailed observations of tectonic tremor from mid-2001 through 2008 that indicate deep fault slip both before and after the Parkfield earthquake that cannot be detected with surface geodetic instruments. While there is no obvious short-term precursor, I find unidirectional tremor migration accompanied by elevated tremor rates in the 3 months prior to the earthquake, which suggests accelerated creep on the fault ∼16 km beneath the eventual earthquake hypocenter.

Historic surface slip along the San Andreas Fault near Parkfield, California

USGS Publications Warehouse

Lienkaemper, J.J.; Prescott, W.H.

1989-01-01

The Parkfield Earthquake Prediction Experiment is focusing close attention on the 44-km-long section of the San Andreas fault that last ruptured seismically in 1966 (Ms 6.0). The 20-km-long central segment of the 1966 Parkfield rupture, extending from the mainshock epicenter at Middle Mountain southeastward to Gold Hill, forms a 1- to 2-km salient northeastward away from the dominant N40??W strike. Following the 1966 earthquake afterslip, aseismic slip has been nearly constant. Moderate Parkfield earthquakes have recurred on average every 21 years since 1857, when a great earthquake (M ~ 8) ruptured at least as far north as the southern Parkfield segment. Many measurements of slip have been made near Parkfield since 1966. Nevertheless, much of the history of surface slip remained uncertain, especially the total amount associated with the 1966 event. In 1985 we measured accumulated slip on the four oldest cultural features offset by the fault along the 1966 Parkfield rupture segment. -from Authors

Numerical modeling of porosity waves in the Nankai accretionary wedge décollement, Japan: implications for aseismic slip

NASA Astrophysics Data System (ADS)

Joshi, Ajit; Appold, Martin S.

2017-01-01

Seismic and hydrologic observations of the Nankai accretionary wedge décollement, Japan, show that overpressures at depths greater than ˜2 km beneath the seafloor could have increased to near lithostatic values due to sediment compaction and diagenesis, clay dehydration, and shearing. The resultant high overpressures are hypothesized then to have migrated in rapid surges or pulses called `porosity waves' up the dip of the décollement. Such high velocities—much higher than expected Darcy fluxes—are possible for porosity waves if the porous media through which the waves travel are deformable enough for porosity and permeability to increase strongly with increasing fluid pressure. The present study aimed to test the hypothesis that porosity waves can travel at rates (kilometers per day) fast enough to cause aseismic slip in the Nankai décollement. The hypothesis was tested using a one-dimensional numerical solution to the fluid mass conservation equation for elastic porous media. Results show that porosity waves generated at depths of ˜2 km from overpressures in excess of lithostatic pressure can propagate at rates sufficient to account for aseismic slip along the décollement over a wide range of hydrogeological conditions. Sensitivity analysis showed porosity wave velocity to be strongly dependent on specific storage, fluid viscosity, and the permeability-depth gradient. Overpressure slightly less than lithostatic pressure could also produce porosity waves capable of traveling at velocities sufficient to cause aseismic slip, provided that hydrogeologic properties of the décollement are near the limits of their geologically reasonable ranges.

Slow slip events in the early part of the earthquake cycle

NASA Astrophysics Data System (ADS)

Voss, Nicholas K.; Malservisi, Rocco; Dixon, Timothy H.; Protti, Marino

2017-08-01

In February 2014 a Mw = 7.0 slow slip event (SSE) took place beneath the Nicoya Peninsula, Costa Rica. This event occurred 17 months after the 5 September 2012, Mw = 7.6, earthquake and along the same subduction zone segment, during a period when significant postseismic deformation was ongoing. A second SSE occurred in the middle of 2015, 21 months after the 2014 SSE and 38 months after the earthquake. The recurrence interval for Nicoya SSEs was unchanged by the earthquake. However, the spatial distribution of slip for the 2014 event differed significantly from previous events, having only deep ( 40 km) slip, compared to previous events, which had both deep and shallow slip. The 2015 SSE marked a return to the combination of deep plus shallow slip of preearthquake SSEs. However, slip magnitude in 2015 was nearly twice as large (Mw = 7.2) as preearthquake SSEs. We employ Coulomb Failure Stress change modeling in order to explain these changes. Stress changes associated with the earthquake and afterslip were highest near the shallow portion of the megathrust, where preearthquake SSEs had significant slip. Lower stress change occurred on the deeper parts of the plate interface, perhaps explaining why the deep ( 40 km) region for SSEs remained unchanged. The large amount of shallow slip in the 2015 SSE may reflect lack of shallow slip in the prior SSE. These observations highlight the variability of aseismic strain release rates throughout the earthquake cycle.