Background Slip Research Articles

Frictional Strengthening Explored During Non‐Steady State Shearing: Implications for Fault Stability and Slip Event Recurrence Time

AbstractOn natural faults that host repeating slip events, the inter‐event loading time is quite large compared to the slip event duration. Since most friction studies focus on steady‐state frictional behavior, the fault loading phase is not typically examined. Here, we employ a method specifically designed to evaluate fault strength evolution during active loading, under shear driving rates as low as 10−10 m/s, on natural fault gouge samples from the Waikukupa Thrust in southern New Zealand. These tests reveal that in the early stages of loading following a slip event, there is a period of increased stability, which fades with accumulated slip. In the framework of rate‐ and state‐dependent friction laws, this temporary stable phase exists as long as slip is less than the critical slip distance and the elapsed time is less than the value of the state variable at steady state. These observations indicate a minimum earthquake recurrence time, which depends on the field value of the critical slip distance and the background slip rate. We compare estimates of minimum earthquake recurrence times with the recurrence times of repeating large earthquakes on the Alpine Fault in southern New Zealand and repeating small‐magnitude earthquakes on the San Andreas Fault system in California. We find that the observed recurrence times are mostly longer than the predicted minimum values, and exceptions in the San Andreas system may be explained by elevated slip rates due to larger earthquakes in this region.

Generation and Validation of Broadband Synthetic P Waves in Semistochastic Models of Large Earthquakes

ABSTRACTWe present an approach for generating stochastic scenario rupture models and semistochastic broadband seismic waveforms that include validated P waves, an important feature for application to early warning systems testing. There are few observations of large magnitude earthquakes available for development and refinement of early warning procedures; thus, simulated data are a valuable supplement. We demonstrate the advantage of using the Karhunen–Loève expansion method for generating stochastic scenario rupture models, as it allows the user to build in desired spatial qualities, such as a slip inversion, as a mean background slip model. For waveform computation, we employ a deterministic approach at low frequencies (<1 Hz) and a semistochastic approach at high frequencies (>1 Hz). Our approach follows Graves and Pitarka (2010) and extends to model P waves. We present the first validation of semistochastic broadband P waves, comparing our waveforms against observations of the 2014 Mw 8.1 Iquique, Chile, earthquake in the time domain and across frequencies of interest. We then consider the P waves in greater detail, using a set of synthetic waveforms generated for scenario ruptures in the Cascadia subduction zone. We confirm that the time-dependent synthetic P-wave amplitude growth is consistent with previous analyses and demonstrate how the data could be used to simulate earthquake early warning procedures.

Estimation of Strong Motion Generation Area for the 2004 Parkfield Earthquake Using Empirical Green’s Function Method

The empirical Green’s function (EGF) method is used to simulate the 2004 Parkfield earthquake (Mw = 6.1). The strong motion generation areas (SMGA) are estimated to reproduce near-source ground motions in a broadband frequency range of 0.25–10 Hz. In this simulation, the EGF uses the small events that occurred in 2005 (Mw = 4.5), 2012 (Mw = 4.4), 2010 (Mw = 4), and 2007 (Mw = 4). The source spectral ratio estimates the seismic moment ratio between large and small events whose corner frequencies obey the omega-squared scaling law. The estimated SMGA is 12 km in length by 8 km in width and is located near the large slip area of the waveform inversion. The result shows that the rupture propagates radially from the southeastern bottom part of the SMGA toward the northwestern shallow direction. In addition, a characterized source model comprising two asperities and a background slip area is proposed. The rupture propagates from asperity1 toward asperity2 with 2 s time delay. The characterized source model simulates more accurately than the single SMGA. The comparison between the 5% damped pseudo acceleration response spectra of our simulations verifies a better compatibility for more than 2 Hz.

Numerical Simulation of M9 Megathrust Earthquakes in the Cascadia Subduction Zone

We estimate ground motions in the Pacific Northwest urban areas during M9 subduction scenario earthquakes on the Cascadia megathrust by simulating wave propagation from an ensemble of kinematic source descriptions. Velocities and densities in our computational mesh are defined by integrating the regional Cascadia Community Velocity Model (CVM) v1.6 (Stephenson et al. P-and S-wave velocity models incorporating the Cascadia subduction zone for 3D earthquake ground motion simulations—update for open-file report 2007–1348, US Geological Survey, 2017) including the ocean water layer with a local velocity model of the Georgia basin (Molnar, Predicting earthquake ground shaking due to 1D soil layering and 3D basin structure in SW British Columbia, Canada, 2011), including additional near-surface velocity information. We generate six source realizations, each consisting of a background slip distribution with correlation lengths, rise times and rupture velocities consistent with data from previous megathrust earthquakes (e.g., 2011 M 9 Tohoku or 2010 M 8.8 Maule). We then superimpose M ~ 8 subevents, characterized by short rise times and high stress drops on the background slip model to mimic high-frequency strong ground motion generation areas in the deeper portion of the rupture (Frankel, Bull Seismol Soc Am 107(1):372–386, 2017). The wave propagation is simulated using the discontinuous mesh (DM) version of the AWP finite difference code. We simulate frequencies up to 1.25 Hz, using a spatial discretization of 100 m in the fine grid, resulting in surface grid dimensions of 6540 × 10,728 mesh points. At depths below 8 km, the grid step increases to 300 m. We obtain stable and accurate results for the DM method throughout the simulation time of 7.5 min as verified against a solution obtained with a uniform 100 m grid spacing. Peak ground velocities (PGVs) range between 0.57 and 1.0 m/s in downtown Seattle and between 0.25 and 0.54 m/s in downtown Vancouver, while spectral accelerations at 2 s range between 1.7 and 3.6 m/s2 and 1.0 and 1.3 m/s2, respectively. These long-period ground motions are not significantly reduced if plastic Drucker-Prager yielding in shallow cohesionless sediments is taken into account. Effects of rupture directivity are significant at periods of ~ 10 s, but almost absent at shorter periods. We find that increasing the depth extent of the subducting slab from the truncation at 60 km in the Cascadia CVM version 1.6 to ~ 100 km increases the PGVs by 15% in Seattle and by 40% in Vancouver.

Variations in precursory slip behavior resulting from frictional heterogeneity

Precursory seismicity is often observed before a large earthquake. Small foreshocks occur within the mainshock rupture area, which cannot be explained by simple models that assume homogeneous friction on the entire fault. In this study, we consider a frictionally heterogeneous fault model, motivated by recent observations of geologic faults and slow earthquakes. This study investigates slip behavior on faults governed by a rate- and state-dependent friction law. We consider a finite linear fault consisting of alternating velocity-weakening zones (VWZs) and velocity-strengthening zones (VSZs). Our model generates precursory slip before the mainshock that ruptures the entire fault, though the activity level of the precursory slip depends on the frictional parameters. We investigate variations in precursory slip behavior, which we characterize quantitatively by the background slip acceleration and seismic radiation, using parameter studies of the a value of the VWZs and VSZs. The results reveal that precursory slip is very small when VWZs are strongly locked and when VSZs consume only a small amount of energy during seismic slip. Precursory slip is significant around the stability boundary of the fault. Furthermore, the type of precursory slip (seismic or aseismic) is controlled by the amplitude of the frictional heterogeneity. Active foreshocks obeying an inverse Omori law associated with background aseismic slip can be interpreted as the nucleation of the mainshock, though this is different from classical nucleation because the monotonic increase in slip velocity is significantly perturbed by the occurrence of foreshocks. Frictional heterogeneity also affects interseismic slip behavior. Modeled variations in precursory slip behavior and interseismic activity can qualitatively explain the along-dip and among-subduction-zone variations in real seismicity patterns. Because even simple frictional heterogeneity produces complex seismicity, it is necessary to further investigate the slip behavior of frictionally heterogeneous faults, which could be utilized for modeling various real seismicity patterns.

An image segmentation based algorithm for imaging of slow slip earthquakes

Laboratory experiments next to a variety of observations, especially in subduction zones, have explored the existence of a premonitory stable slow slip growth phase preceding large earthquakes. These phenomena play an important role in the earthquake cycle and thus precise imaging and monitoring of these events are of great significance. In the literature, ENIF (extended network inversion filter) has been proposed as a rigorous algorithm capable of isolating signal from different types of noise and thereby provides us with deep insight into spatio-temporal evolution of slow slip events. Despite its considerable advantages, ENIF still suffers from some limitations. ENIF applies Tikhonov method of regularization with a quadratic form of cost function. While anomalous slip regions have clear contrast with the background slip in reality, Tikhonov regularization tends to over smooth (globally smooth) the slipping portion on the estimated images. In order to avoid over smoothing phenomenon, we have incorporated into ENIF an image segmentation step which tries to preserve edges of slow-slip event. As a second limitation, due to the nonlinearity imposed by such constraint as non-negativity of slip rate, uncertainty propagation through model is not simple. As the core of ENIF, EKF (extended Kalman filter), performs uncertainty propagation by linearization of nonlinear model using Jacobian and Hessian matrices. As an alternative for EKF, we have also investigated the application of UKF (unscented Kalman filter) which uses UT (unscented transform) for uncertainty propagation. Finally, we tested our proposed algorithm using a low signal to noise ratio synthetic data set. The results show a significant improvement in the performance of ENIF when the segmentation step is incorporated into the algorithm.

Evaluating a Kinematic Method for Generating Broadband Ground Motions for Great Subduction Zone Earthquakes: Application to the 2003 M w 8.3 Tokachi‐Oki Earthquake

Abstract We compare broadband synthetic seismograms with recordings of the 2003 M w 8.3 Tokachi‐Oki earthquake to evaluate a compound rupture model, in which slip on the fault consists of multiple high‐stress‐drop asperities superimposed on a background slip distribution with longer rise times. Low‐frequency synthetics ( 1 Hz) stochastic synthetics using a matched filter at 1 Hz. We show that this compound rupture model and overall approach accurately reproduces waveform envelopes and observed response spectral accelerations (SAs) from the Tokachi‐Oki event. We find that sufficiently short subfault rise times (i.e., Electronic Supplement: Additional waveforms and subevent slip distributions.

Seismicity rate increases associated with slow slip episodes prior to the 2012 Mw 7.4 Ometepec earthquake

The March 20, 2012 Mw 7.4 Ometepec earthquake in the Oaxaca region of Southern Mexico provides a unique opportunity to examine whether subtle changes in seismicity, tectonic tremor, or slow slip can be observed prior to a large earthquake that may illuminate changes in stress or background slip rate. Continuous Global Positioning System (cGPS) data reveal a 5-month-long slow slip event (SSE) between ∼20 and 35 km depth that migrated toward and reached the vicinity of the mainshock a few weeks prior to the earthquake. Seismicity in Oaxaca is examined using single station tectonic tremor detection and multi-station waveform template matching of earthquake families. An increase in seismic activity, detected with template matching using aftershock waveforms, is only observed in the weeks prior to the mainshock in the region between the SSE and mainshock. In contrast, a SSE ∼15 months earlier occurred at ∼25–40 km depth and was primarily associated with an increase in tectonic tremor. Together, these observations indicate that in the Oaxaca region of Mexico shallower slow slip promotes elevated seismicity rates, and deeper slow slip promotes tectonic tremor. Results from this study add to a growing number of published accounts that indicate slow slip may be a common pre-earthquake signature.

Modeling Strong‐Motion Recordings of the 2010Mw 8.8 Maule, Chile, Earthquake with High Stress‐Drop Subevents and Background Slip

Strong‐motion recordings of the M w 8.8 Maule earthquake were modeled using a compound rupture model consisting of (1) a background slip distribution with large correlation lengths, relatively low slip velocity, and long peak rise time of slip of about 10 s and (2) high stress‐drop subevents (asperities) on the deeper portion of the rupture with moment magnitudes 7.9–8.2, high slip velocity, and rise times of slip of about 2 s. In this model, the high‐frequency energy is not produced in the same location as the peak coseismic slip, but is generated in the deeper part of the rupture zone. Using synthetic seismograms generated for a plane‐layered velocity model, I find that the high stress‐drop subevents explain the observed Fourier spectral amplitude from about 0.1 to 1.0 Hz. Broadband synthetics (0–10 Hz) were calculated by combining deterministic synthetics derived from the background slip and asperities (≤1 Hz) with stochastic synthetics generated only at the asperities (≥1 Hz). The broadband synthetics produced response spectral accelerations with low bias compared to the data, for periods of 0.1–10 s. A subevent stress drop of 200–350 bars for the high‐frequency stochastic synthetics was found to bracket the observed spectral accelerations at frequencies greater than 1 Hz. For most of the stations, the synthetics had durations of the Arias intensity similar to the observed records.

A slower fault may produce a smaller preseismic moment rate: Non‐1/tf acceleration of moment rate during nucleation and dependency on the background slip rate

A recent global study has revealed that seismicity near the hypocenter prior to large earthquakes, which could be a proxy for preseismic moment rate, accelerates before interplate earthquakes, while it rarely does before intraplate earthquakes. Understanding the amplitude of preseismic deformation is important in assessing the possibility of its detection. For a class of rate‐state friction laws without a characteristic speed‐related parameter (e.g., aging law and slip law), a dimensional analysis has shown that if the moment rate increases more mildly than 1/tf where tf is the time‐to‐failure, then the amplitude of preseismic moment rate is smaller for a smaller quasistatic slip rate. Three‐dimensional numerical simulations have revealed that the aging law yields 1/tf acceleration, while the slip law causes milder acceleration. If the latter is the case, faults with very low long‐term slip rates (e.g., intraplate faults) may have very small preseismic moment rates.

Material Evidence for Multiple Firings of Ancient Athenian Red‐Figure Pottery

The production of Athenian fine ware pottery, produced between the 6th and 4th centuries B.C., required alternating the high‐temperature kiln between oxidative and reductive environments during a single firing to create the iconic red and black decorative scenes. Here, we show that the production of this pottery was even more complex, with vessels subjected to two, or possibly more, firings in the kiln, with applications of slip between each firing. On a representative sherd, we compared three painted black decorative features—relief line, contour line, and background slip. Scanning transmission electron microscopy (STEM) of the slips revealed that the relief line had a more melted microstructure than either the contour line or background slip. By characterizing the chemistry and micromorphology of the slips, we find that the relief line microstructure could only be produced through a separate firing, at a hotter temperature, than the other two decorative features.

The role of velocity‐neutral creep on the modulation of tectonic tremor activity by periodic loading

Slow slip events and associated non‐volcanic tremors are sensitive to oscillatory stress perturbations, such as those induced by tides or seismic surface waves. Slow slip events and tremors are thought to occur near the seismic‐aseismic transition regions of active faults, where the differencea − b = ∂μ/∂lnVbetween the sensitivity of friction to slip rate and fault state, which characterizes the stability of slip, can be arbitrarily small. We investigate the response of a velocity‐strengthening fault region to oscillatory loads through analytical approximations and spring‐slider simulations. We find that fault areas that are near velocity‐neutral at steady‐state, i.e., ∂μ / ∂lnV ≈ 0, are highly sensitive to cyclic loads: oscillatory stress perturbations in a certain range of periods induce large transient slip velocities. These aseismic transients can in turn trigger tremor activity with enhanced oscillatory modulation. In this sensitive regime, a harmonic Coulomb stress perturbation of amplitude ΔS causes a slip rate perturbation varying as eΔS/(a−b)σ, where σ is the effective normal stress. This result is in agreement with observations of the relationship between tremor rate and amplitude of stress perturbations for tremors triggered by passing seismic waves. Our model of tremor modulation mediated by transient creep does not require extremely high pore fluid pressure and provides a framework to interpret the sensitivity and phase of tidally modulated tremors observed in Parkfield and Shikoku in terms of spatial variations of friction parameters and background slip rate.

Source Characterization for Broadband Ground-Motion Simulation: Kinematic Heterogeneous Source Model and Strong Motion Generation Area

We estimate strong motion generation areas that reproduce near-source ground motions in a broadband frequency range (0.2-10 Hz) using the empirical Green's function technique. The strong motion generation areas are defined as extended areas with relatively large slip velocities within a total rupture area. Four M 6 class (the 1997 Kagoshima on March and May, 1997 Yamaguchi, and 1998 Iwate) and several moderate-size earthquakes in Japan were analyzed. We examine the relationship between the strong motion generation area and the area of asperities, which is characterized based on heterogeneous slip distributions estimated from low-frequency (<1 Hz) waveform inversions (Somerville et al. , 1999). We performed waveform fitting in acceleration, velocity, and displacement, then obtained the strong motion generation area occupying about a quarter of the total rupture area. The size and position of the strong motion generation area coincide with those of characterized asperities. We find self-similar scaling of seismic moment to both the size of the strong motion generation area and the rise time for a magnitude range of analyzed earthquakes. Based on our results, we propose a characterized source model for the broadband ground-motion simulation, which consists of strong motion generation areas with large slip velocities and a background slip area with a small slip velocity. Waveform modeling of the characterized source model suggests that the strong motion generation areas play important roles in simulating broadband ground motions and have the potential to be an about 10-MPa stress-release area on the fault. Manuscript received 22 August 2002.

Technological investigation of the decorative coatings on Yangshao pottery from Henan, China

Yangshao pottery sherds dated to c. 4200 bc from two archaeological sites in Henan province, China, were examined to determine the technological aspects of their decoration. The selected samples comprised plain and painted pottery, some of which have a white background slip and black or red decoration on the surface. Various analytical techniques including SEM/EDS, XRD and FTIR were applied in the study. Observations show that some pieces were wheel‐thrown. The refinement of the body fabric varies. Surface finishing includes simple burnishing, slipping and brushing. The black decoration on the surface is manganese black. The white background layer is kaolinitic clay slip. No organic binders or decorations were detected. The results of these analyses will help our understanding of the production steps involved in the craft production of Yangshao pottery.

Anteroposterior dynamic balance reactions to circular movement of the visual field without ocular pursuit.

It has been hypothesised that during natural vision, retinal image motion and eye motion information exert antagonistic effects on balance that cancel each other out. In a previous study, we found that pursuit movements of the eye unaccompanied by background slip on the retina induced strong, stereotyped anteroposterior sway in dynamic balance. In the present work, we investigated the influence of background retinal image slip alone on anteroposterior dynamic balance. We expected that retinal slip without ocular pursuit would induce large anteroposterior sway. Circular slip without rotation of the whole visual field image on the retina was obtained using a rotating prism. In this study, carried out on monocular vision, eye movement was prevented by fixating a stationary target seen through a small aperture in the prism. Anteroposterior dynamic balance conditions were elicited by a rocking platform. Strong, stereotyped anteroposterior sway accompanied by a visual illusion of target motion was observed. The induced sway was synchronised to the stimulus and led to a significant disturbance of balance. The similarity in magnitude of the effects of both retinal image motion studied here and eye motion information obtained previously, and their opposing phases, strongly support the initial hypothesis that the antagonistic effects of the two sources of information provided by normal ocular pursuit could interact in an additive algebraic mode.

Eureka Peak fault afterslip following the 28 June 1992 Landers earthquake

Abstract Although surface ruptures of the principal Landers faults apparently slipped co-seismically, surface slip on the 10-km-long Eureka Peak fault has continued at a decreasing rate during the year following the mainshock. Three creepmeters installed on the fault as early as 4 days after the mainshock have recorded afterslip amplitudes of more than 8 cm on the northern half of the mapped fault. The maximum observed slip 1 yr after the earthquake is 23 ± 2 cm, and an additional 1 to 3 cm of afterslip may occur on the fault in the next decade if the current rate of slip decrease remains steady. Should this occur, it will bring surface slip amplitudes close to the ∼25 cm of slip inferred from geodetic data to have occurred at depth on the northern 4 km of the fault. Creepmeters on the fault record both episodic creep events and slow background slip, but, although changes in creep rate occur, they do not correlate well with seismic moment release in volumes near the fault zone. As is the case with many afterslip processes and other creeping faults in California, the Eureka Peak fault cuts through sediments. It is possible that such sedimentary layers, and corresponding depths to which afterslip is occurring, are relatively shallow or that the fault lacks a well-developed gouge zone near sites where episodic creep is absent.

A two-layer model for aseismic slip on the Superstition Hills fault, California

Abstract A new differential creepmeter on the Superstition Hills fault reveals that afterslip at approximately 28 mm/year consists of episodic creep events super-imposed on a slow stable-slip component of 2.4 mm/year. The relative contribution from this background slip (8.5% of the total aseismic surface slip) has remained approximately constant for the past 2 years. To account for this observed behavior, we proposed a two-layer model for aseismic slip in which stable-sliding occurs from the surface to a transition depth, below which episodic creep events are initiated. These creep events propagate to the surface through the stable-sliding layer. From the ratio of background slip velocity to total afterslip velocity, we estimate the ratio of depths of stable-sliding to episodic-slip regions to be approximately 1:10. Thus, assuming that episodic slip extends to 3 km, we infer that background surface creep is moderated by processes in the uppermost 300 m. We speculate that the transition depth is sensitive to applied fault-normal stresses and suggest that the ratio of stable-sliding to episodic-slip velocities may provide an indication of secular variations in tectonic stress. Since the background creep is evidently continuous and little affected by episodic events, its velocity may provide a sensitive indication of applied tectonic strain in the period range days to years.

Responses of monkey nucleus of the optic tract neurons during pursuit and fixation

We recorded 101 neurons in the nucleus of the optic tract (NOT) of 3 rhesus monkeys. The neurons were tested in a variety of oculomotor paradigms. This report focusses on the modulation of NOT neuronal activity during smooth pursuit eye movements. A small horizontally moving spot (<1°) elicited a directionally specific response during fixation and revealed thereby the extent of the receptive fields. During pursuit NOT neurons are coding for target slip. If eye speed exceeds target speed the direction of retinal slip is reversed and in accordance with their directional sensitivity NOT neurons immediately change their activity. This result proves the slip transfer function as well as the independence from eye movement signals of NOT neurons. During pursuit across a structured background some neurons are still coding for target slip whereas other neurons are coding for background slip. These two groups of neurons can also be distinguished by their response during fixation. The response of a target slip neuron to a background movement is cancelled, whereas the response of a background neuron is not affected by fixation. There is no difference in size of receptive fields for these two groups of neurons. We conclude from our findings that directionally selective cells in the monkey NOT may provide input to the pursuit system as well as to the optokinetic system. This dichotomy may also be reflected in different efferent projections: to the nucleus reticularis tegmenti pontis and to the inferior olive, respectively. A similar notion was introduced by the late Maekawa for the rabbit's NOT.