On the role of fluids in stick-slip dynamics of saturated granular fault gouge using a coupled computational fluid dynamics-discrete element approach

NASA Astrophysics Data System (ADS)

Dorostkar, Omid; Guyer, Robert A.; Johnson, Paul A.; Marone, Chris; Carmeliet, Jan

2017-05-01

The presence of fault gouge has considerable influence on slip properties of tectonic faults and the physics of earthquake rupture. The presence of fluids within faults also plays a significant role in faulting and earthquake processes. In this paper, we present 3-D discrete element simulations of dry and fluid-saturated granular fault gouge and analyze the effect of fluids on stick-slip behavior. Fluid flow is modeled using computational fluid dynamics based on the Navier-Stokes equations for an incompressible fluid and modified to take into account the presence of particles. Analysis of a long time train of slip events shows that the (1) drop in shear stress, (2) compaction of granular layer, and (3) the kinetic energy release during slip all increase in magnitude in the presence of an incompressible fluid, compared to dry conditions. We also observe that on average, the recurrence interval between slip events is longer for fluid-saturated granular fault gouge compared to the dry case. This observation is consistent with the occurrence of larger events in the presence of fluid. It is found that the increase in kinetic energy during slip events for saturated conditions can be attributed to the increased fluid flow during slip. Our observations emphasize the important role that fluid flow and fluid-particle interactions play in tectonic fault zones and show in particular how discrete element method (DEM) models can help understand the hydromechanical processes that dictate fault slip.

Secular Variation in the Storage and Dissipation of Elastic Strain Energy Along the Central Altyn Tagh Fault (86-88.5°E), NW China

NASA Astrophysics Data System (ADS)

Cowgill, E.; Gold, R. D.; Arrowsmith, R.; Friedrich, A. M.

2015-12-01

In elastic rebound theory, hazard increases as interseismic strain rebuilds after rupture. This model is challenged by the temporal variation in the pacing of major earthquakes that is both predicted by mechanical models and suggested by some long paleoseismic records (e.g., 1-3). However, the extent of such behavior remains unclear due to a lack of long (5-25 ky) records of fault slip. Using Monte Carlo analysis of 11 offset landforms, we determined a 16-ky record of fault slip for the active, left-lateral Altyn Tagh fault, which bounds the NW margin of the Tibetan Plateau. This history reveals a pulse of accelerated slip between 6.4 and 6.0 ka, during which the fault slipped 9 +14/-2 m at a rate of 23 +35/-5 mm/y, or ~3x the 16 ky average of 8.1 +1.2/-0.9mm/y. These two modes of earthquake behavior suggest temporal variation in the rates of stress storage and release. The simplest explanation for the pulse is a cluster of 2-8 Mw > 7.5 earthquakes. Such supercyclicity has been reported for the Sunda (4) and Cascadia (3) megathrusts, but contrasts with steady slip along the strike-slip Alpine fault (5), for example. A second possibility is that the pulse reflects a single, unusually large rupture. However, this Black Swan event is unlikely: empirical scaling relationships require a Mw 8.2 rupture of the entire 1200-km-long ATF to produce 7 m of average slip. Likewise, Coulomb stress change from rupture on the adjacent North Altyn fault is of modest magnitude and overlap with the ATF. Poor temporal correlation between precipitation and the slip pulse argues against climatically modulated changes in surface loading (lakes/ice) or pore-fluid pressure. "Paleoslip" studies such as this sacrifice the single-event resolution of paleoseismology in exchange for long records that quantify both the timing and magnitude of fault slip averaged over multiple ruptures, and are essential for documenting temporal variations in fault slip as we begin to use calibrated physical models of the earthquake cycle to forecast time-dependent earthquake hazard (e.g., 6,7). 1. Weldon et al., 2004 GSA Today 14, 4; 2. Rockwell et al., 2015, PAGEOPH, 172, 1143; 3. Goldfinger et al., 2013, SRL, 84, 24; 4. Sieh et al., 2008, Science, 322, 1674; 5. Berryman et l., 2012, Science, 336, 1690; 6. Barbot et al., 2012, Science, 336, 707; 7. Field, 2015, BSSA, 105, 544.

Observed source parameters for dynamic rupture with non-uniform initial stressand relatively high fracture energy

USGS Publications Warehouse

Beeler, Nicholas M.; Kilgore, Brian D.; McGarr, Arthur F.; Fletcher, Jon Peter B.; Evans, John R.; Steven R. Baker,

2012-01-01

We have conducted dynamic rupture propagation experiments to establish the relations between in-source stress drop, fracture energy and the resulting particle velocity during slip of an unconfined 2 m long laboratory fault at normal stresses between 4 and 8 MPa. To produce high fracture energy in the source we use a rough fault that has a large slip weakening distance. An artifact of the high fracture energy is that the nucleation zone is large such that precursory slip reduces fault strength over a large fraction of the total fault length prior to dynamic rupture, making the initial stress non-uniform. Shear stress, particle velocity, fault slip and acceleration were recorded coseismically at multiple locations along strike and at small fault-normal distances. Stress drop increases weakly with normal stress. Average slip rate depends linearly on the fault strength loss and on static stress drop, both with a nonzero intercept. A minimum fracture energy of 1.8 J/m2 and a linear slip weakening distance of 33 μm are inferred from the intercept. The large slip weakening distance also affects the average slip rate which is reduced by in-source energy dissipation from on-fault fracture energy.Because of the low normal stress and small per event slip (∼86 μm), no thermal weakening such as melting or pore fluid pressurization occurs in these experiments. Despite the relatively high fracture energy, and the very low heat production, energy partitioning during these laboratory earthquakes is very similar to typical earthquake source properties. The product of fracture energy and fault area is larger than the radiated energy. Seismic efficiency is low at ∼2%. The ratio of apparent stress to static stress drop is ∼27%, consistent with measured overshoot. The fracture efficiency is ∼33%. The static and dynamic stress drops when extrapolated to crustal stresses are 2–7.3 MPa and in the range of typical earthquake stress drops. As the relatively high fracture energy reduces the slip velocities in these experiments, the extrapolated average particle velocities for crustal stresses are 0.18–0.6 m/s. That these experiments are consistent with typical earthquake source properties suggests, albeit indirectly, that thermal weakening mechanisms such as thermal pressurization and melting which lead to near complete stress drops, dominate earthquake source properties only for exceptional events unless crustal stresses are low.

Late Quaternary history of the Owens Valley fault zone, eastern California, and surface rupture associated with the 1872 earthquake

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

Beanland, S.; Clark, M.M.

1993-04-01

The right-lateral Owens Valley fault zone (OVFZ) in eastern California extends north about 100 km from near the northwest shore of Owens Lake to beyond Big Pine. It passes through Lone Pine near the eastern base of the Alabama Hills and follows the floor of Owens Valley northward to the Poverty Hills, where it steps 3 km to the left and continues northwest across Crater Mountain and through Big Pine. Data from one site suggest an average net slip rate for the OVFZ of 1.5 [+-] 1 mm/yr for the past 300 ky. Several other sites yield an average Holocenemore » net slip rate of 2 [+-] 1 mm/yr. The OVFZ apparently has experienced three major Holocene earthquakes. The minimum average recurrence interval is 5,000 years at the subsidiary Lone Pine fault, whereas it is 3,300 to 5,000 years elsewhere along the OVFZ. The prehistoric earthquakes are not dated, so an average recurrence interval need not apply. However, roughly equal (characteristic) displacement apparently happened during each Holocene earthquake. The Owens Valley fault zone accommodates some of the relative motion (dextral shear) between the North American and Pacific plates along a discrete structure. This shear occurs in the Walker Lane belt of normal and strike-slip faults within the mainly extensional Basin and Range Province. In Owens Valley displacement is partitioned between the OVFZ and the nearby, subparallel, and purely normal range-front faults of the Sierra Nevada. Compared to the OVFZ, these range-front normal faults are very discontinuous and have smaller Holocene slip rates of 0.1 to 0.8 mm/yr, dip slip. Contemporary activity on adjacent faults of such contrasting styles suggests large temporal fluctuations in the relative magnitudes of the maximum and intermediate principal stresses while the extension direction remains consistently east-west.« less

Surface faulting along the Superstition Hills fault zone and nearby faults associated with the earthquakes of 24 November 1987

USGS Publications Warehouse

Sharp, R.V.

1989-01-01

The M6.2 Elmore Desert Ranch earthquake of 24 November 1987 was associated spatially and probably temporally with left-lateral surface rupture on many northeast-trending faults in and near the Superstition Hills in western Imperial Valley. Three curving discontinuous principal zones of rupture among these breaks extended northeastward from near the Superstition Hills fault zone as far as 9km; the maximum observed surface slip, 12.5cm, was on the northern of the three, the Elmore Ranch fault, at a point near the epicenter. Twelve hours after the Elmore Ranch earthquake, the M6.6 Superstition Hills earthquake occurred near the northwest end of the right-lateral Superstition Hills fault zone. We measured displacements over 339 days at as many as 296 sites along the Superstition Hills fault zone, and repeated measurements at 49 sites provided sufficient data to fit with a simple power law. The overall distributions of right-lateral displacement at 1 day and the estimated final slip are nearly symmetrical about the midpoint of the surface rupture. The average estimated final right-lateral slip for the Superstition Hills fault zone is ~54cm. The average left-lateral slip for the conjugate faults trending northeastward is ~23cm. The southernmost ruptured member of the Superstition Hills fault zone, newly named the Wienert fault, extends the known length of the zone by about 4km. -from Authors

Coseismic temporal changes of slip direction: the effect of absolute stress on dynamic rupture

USGS Publications Warehouse

Guatteri, Mariagiovanna; Spudich, P.

1998-01-01

We investigate the dynamics of rupture at low-stress level. We show that one main difference between the dynamics of high- and low-stress events is the amount of coseismic temporal rake rotation occurring at given points on the fault. Curved striations on exposed fault surfaces and earthquake dislocation models derived from ground-motion inversion indicate that the slip direction may change with time at a point on the fault during dynamic rupture. We use a 3D boundary integral method to model temporal rake variations during dynamic rupture propagation assuming a slip-weakening friction law and isotropic friction. The points at which the slip rotates most are characterized by an initial shear stress direction substantially different from the average stress direction over the fault plane. We show that for a given value of stress drop, the level of initial shear stress (i.e., the fractional stress drop) determines the amount of rotation in slip direction. We infer that seismic events that show evidence of temporal rake rotations are characterized by a low initial shear-stress level with spatially variable direction on the fault (possibly due to changes in fault surface geometry) and an almost complete stress drop.Our models motivate a new interpretation of curved and cross-cutting striations and put new constraints on their analysis. The initial rake is in general collinear with the initial stress at the hypocentral zone, supporting the assumptions made in stress-tensor inversion from first-motion analysis. At other points on the fault, especially away from the hypocenter, the initial slip rake may not be collinear with the initial shear stress, contradicting a common assumption of structural geology. On the other hand, the later part of slip in our models is systematically more aligned with the average stress direction than the early slip. Our modeling suggests that the length of the straight part of curved striations is usually an upper bound of the slip-weakening distance if this parameter is uniform over the fault plane, and the direction of the late part of slip of curved striations should have more weight in the estimate of initial stress direction.

A New Geological Slip Rate Estimate for the Calico Fault, Eastern California: Implications for Geodetic Versus Geologic Rate Estimates in the Eastern California Shear Zone

NASA Astrophysics Data System (ADS)

Wetmore, P. H.; Xie, S.; Gallant, E.; Owen, L. A.; Dixon, T. H.

2017-12-01

Fault slip rate is fundamental to accurate seismic hazard assessment. In the Mojave Desert section of the Eastern California Shear Zone previous studies have suggested a discrepancy between short-term geodetic and long-term geologic slip rate estimates. Understanding the origin of this discrepancy could lead to better understanding of stress evolution, and improve earthquake hazard estimates in general. We measured offsets in alluvial fans along the Calico fault near Newberry Springs, California, and used exposure age dating based on the cosmogenic nuclide 10Be to date the offset landforms. We derive a mean slip rate of 3.6 mm/yr, representing an average over the last few hundred thousand years, significantly faster than previous estimates. Considering numerous faults in the Mojave Desert and limited geologic slip rate estimates, it is premature to claim a geologic versus geodetic "discrepancy" for the ECSZ. More slip rate data, from all faults with the ECSZ, are needed to provide a statistically meaningful assessment of the geologic rates for each of the faults comprising the ECSZ.

Geomorphology, kinematic history, and earthquake behavior of the active Kuwana wedge thrust anticline, central Japan

NASA Astrophysics Data System (ADS)

Ishiyama, Tatsuya; Mueller, Karl; Togo, Masami; Okada, Atsumasa; Takemura, Keiji

2004-12-01

We combine surface mapping of fault and fold scarps that deform late Quaternary alluvial strata with interpretation of a high-resolution seismic reflection profile to develop a kinematic model and determine fault slip rates for an active blind wedge thrust system that underlies Kuwana anticline in central Japan. Surface fold scarps on Kuwana anticline are closely correlated with narrow fold limbs and angular hinges on the seismic profile that suggest at least ˜1.3 km of fault slip completely consumed by folding in the upper 4 km of the crust. The close coincidence and kinematic link between folded terraces and the underlying thrust geometry indicate that Kuwana anticline has accommodated slip at an average rate of 2.2 ± 0.5 mm/yr on a 27°, west dipping thrust fault since early-middle Pleistocene time. In contrast to classical fault bend folds the fault slip budget in the stacked wedge thrusts also indicates that (1) the fault tip propagated upward at a low rate relative to the accrual of fault slip and (2) fault slip is partly absorbed by numerous bedding plane flexural-slip faults above the tips of wedge thrusts. An historic earthquake that occurred on the Kuwana blind thrust system possibly in A.D. 1586 is shown to have produced coseismic surface deformation above the doubly vergent wedge tip. Structural analyses of Kuwana anticline coupled with tectonic geomorphology at 103-105 years timescales illustrate the significance of active folds as indicators of slip on underlying blind thrust faults and thus their otherwise inaccessible seismic hazards.

Earthquake scaling laws for rupture geometry and slip heterogeneity

NASA Astrophysics Data System (ADS)

Thingbaijam, Kiran K. S.; Mai, P. Martin; Goda, Katsuichiro

2016-04-01

We analyze an extensive compilation of finite-fault rupture models to investigate earthquake scaling of source geometry and slip heterogeneity to derive new relationships for seismic and tsunami hazard assessment. Our dataset comprises 158 earthquakes with a total of 316 rupture models selected from the SRCMOD database (http://equake-rc.info/srcmod). We find that fault-length does not saturate with earthquake magnitude, while fault-width reveals inhibited growth due to the finite seismogenic thickness. For strike-slip earthquakes, fault-length grows more rapidly with increasing magnitude compared to events of other faulting types. Interestingly, our derived relationship falls between the L-model and W-model end-members. In contrast, both reverse and normal dip-slip events are more consistent with self-similar scaling of fault-length. However, fault-width scaling relationships for large strike-slip and normal dip-slip events, occurring on steeply dipping faults (δ~90° for strike-slip faults, and δ~60° for normal faults), deviate from self-similarity. Although reverse dip-slip events in general show self-similar scaling, the restricted growth of down-dip fault extent (with upper limit of ~200 km) can be seen for mega-thrust subduction events (M~9.0). Despite this fact, for a given earthquake magnitude, subduction reverse dip-slip events occupy relatively larger rupture area, compared to shallow crustal events. In addition, we characterize slip heterogeneity in terms of its probability distribution and spatial correlation structure to develop a complete stochastic random-field characterization of earthquake slip. We find that truncated exponential law best describes the probability distribution of slip, with observable scale parameters determined by the average and maximum slip. Applying Box-Cox transformation to slip distributions (to create quasi-normal distributed data) supports cube-root transformation, which also implies distinctive non-Gaussian slip distributions. To further characterize the spatial correlations of slip heterogeneity, we analyze the power spectral decay of slip applying the 2-D von Karman auto-correlation function (parameterized by the Hurst exponent, H, and correlation lengths along strike and down-slip). The Hurst exponent is scale invariant, H = 0.83 (± 0.12), while the correlation lengths scale with source dimensions (seismic moment), thus implying characteristic physical scales of earthquake ruptures. Our self-consistent scaling relationships allow constraining the generation of slip-heterogeneity scenarios for physics-based ground-motion and tsunami simulations.

Coseismic slip variation assessed from terrestrial lidar scans of the El Mayor-Cucapah surface rupture

NASA Astrophysics Data System (ADS)

Gold, Peter O.; Oskin, Michael E.; Elliott, Austin J.; Hinojosa-Corona, Alejandro; Taylor, Michael H.; Kreylos, Oliver; Cowgill, Eric

2013-03-01

We analyze high-resolution (>103 points/m2) terrestrial lidar surveys of the 4 April 2010 El Mayor-Cucapah earthquake rupture (Baja California, Mexico), collected at three sites 12-18 days after the event. Using point cloud-based tools in an immersive visualization environment, we quantify coseismic fault slip for hundreds of meters along strike and construct densely constrained along-strike slip distributions from measurements of offset landforms. Uncertainty bounds for each offset, determined empirically by repeatedly measuring offsets at each site sequentially, illuminate measurement uncertainties that are difficult to quantify in the field. These uncertainties are used to define length scales over which variability in slip distributions may be assumed to reflect either recognizable earthquake mechanisms or measurement noise. At two sites characterized by 2-3 m of concentrated right-oblique slip, repeat measurements yield 2σ uncertainties of ±11-12%. Each site encompasses ∼200 m along strike, and a smoothed linear slip gradient satisfies all measurement distributions, implying along-fault strains of ∼10-3. Conversely, the common practice of defining the slip curve by the local slip maxima distorts the curve, overestimates along-fault strain, and may overestimate actual fault slip by favoring measurements with large, positive, uncertainties. At a third site characterized by 1-2.5 m of diffuse normal slip, repeat measurements of fault throw summed along fault-perpendicular profiles yield 2σ uncertainties of ±17%. Here, a low order polynomial fit through the measurement averages best approximates surface slip. However independent measurements of off-fault strain accommodated by hanging wall flexure suggest that over the ∼200 m length of this site, a linear interpolation through the average values for the slip maxima at either end of this site most accurately represents subsurface displacement. In aggregate, these datasets show that given uncertainties of greater than ±11% (2σ), slip distributions over shorter scales are likely to be less uneven than those derived from a single set of field- or lidar-based measurements. This suggests that the relatively smooth slip curves we obtain over ∼102 m distances reflect real physical phenomena, whereas short wavelength variability over ∼100-101 m distances can be attributed to measurement uncertainty.

Variable rates of late Quaternary strike slip on the San Jacinto fault zone, southern California.

USGS Publications Warehouse

Sharp, R.V.

1981-01-01

3 strike slip displacements of strata with known approximate ages have been measured at 2 locations on the San Jacinto fault zone. Minimum horizontal offset between 5.7 and 8.6km in no more than 0.73Myr NE of Anza indicates 8-12 mm/yr average slip rate since late Pleistocene time. Horizontal slip of 1.7m has been calculated for the youngest sediment of Lake Cahuilla since its deposition 271- 510 yr BP. The corresponding slip rate is 2.8-5.0 mm/yr. Right lateral offset of 10.9m measured on a buried stream channel older than 5060 yr BP but younger than 6820 yr BP yields average slip rates for the intermediate time periods, 400 to 6000 yr BP of 1-2 mm/yr. The rates of slip suggest a relatively quiescent period from about 4000 BC to about 1600 AD.-from Author

Surface slip and off-fault deformation patterns in the 2013 MW 7.7 Balochistan, Pakistan earthquake: Implications for controls on the distribution of near-surface coseismic slip

NASA Astrophysics Data System (ADS)

Zinke, Robert; Hollingsworth, James; Dolan, James F.

2014-12-01

Comparison of 398 fault offsets measured by visual analysis of WorldView high-resolution satellite imagery with deformation maps produced by COSI-Corr subpixel image correlation of Landsat-8 and SPOT5 imagery reveals significant complexity and distributed deformation along the 2013 Mw 7.7 Balochistan, Pakistan earthquake. Average slip along the main trace of the fault was 4.2 m, with local maximum offsets up to 11.4 m. Comparison of slip measured from offset geomorphic features, which record localized slip along the main strand of the fault, to the total displacement across the entire width of the surface deformation zone from COSI-Corr reveals ˜45% off-fault deformation. While previous studies have shown that the structural maturity of the fault exerts a primary control on the total percentage of off-fault surface deformation, large along-strike variations in the percentage of strain localization observed in the 2013 rupture imply the influence of important secondary controls. One such possible secondary control is the type of near-surface material through which the rupture propagated. We therefore compared the percentage off-fault deformation to the type of material (bedrock, old alluvium, and young alluvium) at the surface and the distance of the fault to the nearest bedrock outcrop (a proxy for sediment thickness along this hybrid strike slip/reverse slip fault). We find significantly more off-fault deformation in younger and/or thicker sediments. Accounting for and predicting such off-fault deformation patterns has important implications for the interpretation of geologic slip rates, especially for their use in probabilistic seismic hazard assessments, the behavior of near-surface materials during coseismic deformation, and the future development of microzonation protocols for the built environment.

Synthetic velocity gradient map of the San Francisco Bay region, California, supports use of average block velocities to estimate fault slip rate where effective locking depth is small relative to inter-fault distance

NASA Astrophysics Data System (ADS)

Graymer, R. W.; Simpson, R. W.

2014-12-01

Graymer and Simpson (2013, AGU Fall Meeting) showed that in a simple 2D multi-fault system (vertical, parallel, strike-slip faults bounding blocks without strong material property contrasts) slip rate on block-bounding faults can be reasonably estimated by the difference between the mean velocity of adjacent blocks if the ratio of the effective locking depth to the distance between the faults is 1/3 or less ("effective" locking depth is a synthetic parameter taking into account actual locking depth, fault creep, and material properties of the fault zone). To check the validity of that observation for a more complex 3D fault system and a realistic distribution of observation stations, we developed a synthetic suite of GPS velocities from a dislocation model, with station location and fault parameters based on the San Francisco Bay region. Initial results show that if the effective locking depth is set at the base of the seismogenic zone (about 12-15 km), about 1/2 the interfault distance, the resulting synthetic velocity observations, when clustered, do a poor job of returning the input fault slip rates. However, if the apparent locking depth is set at 1/2 the distance to the base of the seismogenic zone, or about 1/4 the interfault distance, the synthetic velocity field does a good job of returning the input slip rates except where the fault is in a strong restraining orientation relative to block motion or where block velocity is not well defined (for example west of the northern San Andreas Fault where there are no observations to the west in the ocean). The question remains as to where in the real world a low effective locking depth could usefully model fault behavior. Further tests are planned to define the conditions where average cluster-defined block velocities can be used to reliably estimate slip rates on block-bounding faults. These rates are an important ingredient in earthquake hazard estimation, and another tool to provide them should be useful.

Variable slip-rate and slip-per-event on a plate boundary fault: The Dead Sea fault in northern Israel

NASA Astrophysics Data System (ADS)

Wechsler, Neta; Rockwell, Thomas K.; Klinger, Yann

2018-01-01

We resolved displacement on buried stream channels that record the past 3400 years of slip history for the Jordan Gorge (JGF) section of the Dead Sea fault in Israel. Based on three-dimensional (3D) trenching, slip in the past millennium amounts to only 2.7 m, similar to that determined in previous studies, whereas the previous millennium experienced two to three times this amount of displacement with nearly 8 m of cumulative slip, indicating substantial short term variations in slip rate. The slip rate averaged over the past 3400 years, as determined from 3D trenching, is 4.1 mm/yr, which agrees well with geodetic estimates of strain accumulation, as well as with longer-term geologic slip rate estimates. Our results indicate that: 1) the past 1200 years appear to significantly lack slip, which may portend a significant increase in future seismic activity; 2) short-term slip rates for the past two millennia have varied by more than a factor of two and suggest that past behavior is best characterized by clustering of earthquakes. From these observations, the earthquake behavior of the Jordan Gorge fault best fits is a "weak segment model" where the relatively short fault section (20 km), bounded by releasing steps, fails on its own in moderate earthquakes, or ruptures with adjacent segments.

Deterministic estimate of hypocentral pore fluid pressure of the M5.8 Pawnee, Oklahoma earthquake: Lower pre-injection pressure requires lower resultant pressure for slip

NASA Astrophysics Data System (ADS)

Levandowski, W. B.; Walsh, F. R. R.; Yeck, W.

2016-12-01

Quantifying the increase in pore-fluid pressure necessary to cause slip on specific fault planes can provide actionable information for stakeholders to potentially mitigate hazard. Although the M5.8 Pawnee earthquake occurred on a previously unmapped fault, we can retrospectively estimate the pore-pressure perturbation responsible for this event. We first estimate the normalized local stress tensor by inverting focal mechanisms surrounding the Pawnee Fault. Faults are generally well oriented for slip, with instabilities averaging 96% of maximum. Next, with an estimate of the weight of local overburden we solve for the pore pressure needed at the hypocenters. Specific to the Pawnee fault, we find that hypocentral pressure 43-104% of hydrostatic (accounting for uncertainties in all relevant parameters) would have been sufficient to cause slip. The dominant source of uncertainty is the pressure on the fault prior to fluid injection. Importantly, we find that lower pre-injection pressure requires lower resultant pressure to cause slip, decreasing from a regional average of 30% above hydrostatic pressure if the hypocenters begin at hydrostatic pressure to 6% above hydrostatic pressure with no pre-injection fluid. This finding suggests that underpressured regions such as northern Oklahoma are predisposed to injection-induced earthquakes. Although retrospective and forensic, similar analyses of other potentially induced events and comparisons to natural earthquakes will provide insight into the relative importance of fault orientation, the magnitude of the local stress field, and fluid-pressure migration in intraplate seismicity.

Tsunami Source Inversion Using Tide Gauge and DART Tsunami Waveforms of the 2017 Mw8.2 Mexico Earthquake

NASA Astrophysics Data System (ADS)

Adriano, Bruno; Fujii, Yushiro; Koshimura, Shunichi; Mas, Erick; Ruiz-Angulo, Angel; Estrada, Miguel

2018-01-01

On September 8, 2017 (UTC), a normal-fault earthquake occurred 87 km off the southeast coast of Mexico. This earthquake generated a tsunami that was recorded at coastal tide gauge and offshore buoy stations. First, we conducted a numerical tsunami simulation using a single-fault model to understand the tsunami characteristics near the rupture area, focusing on the nearby tide gauge stations. Second, the tsunami source of this event was estimated from inversion of tsunami waveforms recorded at six coastal stations and three buoys located in the deep ocean. Using the aftershock distribution within 1 day following the main shock, the fault plane orientation had a northeast dip direction (strike = 320°, dip = 77°, and rake =-92°). The results of the tsunami waveform inversion revealed that the fault area was 240 km × 90 km in size with most of the largest slip occurring on the middle and deepest segments of the fault. The maximum slip was 6.03 m from a 30 × 30 km2 segment that was 64.82 km deep at the center of the fault area. The estimated slip distribution showed that the main asperity was at the center of the fault area. The second asperity with an average slip of 5.5 m was found on the northwest-most segments. The estimated slip distribution yielded a seismic moment of 2.9 × 10^{21} Nm (Mw = 8.24), which was calculated assuming an average rigidity of 7× 10^{10} N/m2.

Late quaternary slip-rate variations along the Warm Springs Valley fault system, northern Walker Lane, California-Nevada border

USGS Publications Warehouse

Gold, Ryan; dePolo, Craig; Briggs, Richard W.; Crone, Anthony

2013-01-01

The extent to which faults exhibit temporally varying slip rates has important consequences for models of fault mechanics and probabilistic seismic hazard. Here, we explore the temporal behavior of the dextral‐slip Warm Springs Valley fault system, which is part of a network of closely spaced (10–20 km) faults in the northern Walker Lane (California–Nevada border). We develop a late Quaternary slip record for the fault using Quaternary mapping and high‐resolution topographic data from airborne Light Distance and Ranging (LiDAR). The faulted Fort Sage alluvial fan (40.06° N, 119.99° W) is dextrally displaced 98+42/-43 m, and we estimate the age of the alluvial fan to be 41.4+10.0/-4.8 to 55.7±9.2  ka, based on a terrestrial cosmogenic 10Be depth profile and 36Cl analyses on basalt boulders, respectively. The displacement and age constraints for the fan yield a slip rate of 1.8 +0.8/-0.8 mm/yr to 2.4 +1.2/-1.1 mm/yr (2σ) along the northern Warm Springs Valley fault system for the past 41.4–55.7 ka. In contrast to this longer‐term slip rate, shorelines associated with the Sehoo highstand of Lake Lahontan (~15.8  ka) adjacent to the Fort Sage fan are dextrally faulted at most 3 m, which limits a maximum post‐15.8 ka slip rate to 0.2  mm/yr. These relations indicate that the post‐Lahontan slip rate on the fault is only about one‐tenth the longer‐term (41–56 ka) average slip rate. This apparent slip‐rate variation may be related to co‐dependent interaction with the nearby Honey Lake fault system, which shows evidence of an accelerated period of mid‐Holocene earthquakes.

Horizontal surface-slip distribution through several seismic cycles: The Eastern Bogd fault, Gobi-Altai, Mongolia

NASA Astrophysics Data System (ADS)

Kurtz, R.; Klinger, Y.; Ferry, M.; Ritz, J.-F.

2018-06-01

The 1957, MW 8.1, Gobi-Altai earthquake, Southern Mongolia, produced a 360-km-long surface rupture along the Eastern Bogd fault. Cumulative offsets of geomorphic features suggest that the Eastern Bogd fault might produce characteristic slip over the last seismic cycles. Using orthophotographs derived from a dataset of historical aerial photographs acquired in 1958, we measured horizontal offsets along two thirds ( 170 km) of the 1957 left-lateral strike-slip surface rupture. We propose a new empirical methodology to extract the average slip for each past earthquake that could be recognized along successive fault segments, to determine the slip distribution associated with successive past earthquakes. Our results suggest that the horizontal slip distribution of the 1957 Gobi-Altai earthquake is fairly flat, with an average offset of 3.5 m ± 1.3 m. A combination of our lateral measurements with vertical displacements derived from the literature, allows us to re-assess the magnitude of the Gobi-Altai earthquake to be between MW 7.8 and MW 8.2, depending on the depth of the rupture, and related value of the shear modulus. When comparing this magnitude to magnitudes derived from seismic data, it suggests that the rupture may have extended deeper than the 15 km to 20 km usually considered for the seismogenic crust. We observe that some fault segments are more likely than others to record seismic deformation through several seismic cycles, depending on the local rupture complexity and geomorphology. Additionally, our results allow us to model the horizontal slip function for the 1957 Gobi-Altai earthquake and for three previous paleoseismic events along 70% of the studied area. Along about 50% of the fault sections where we could recognize three past earthquakes, our results suggest that the slip per event was similar for each earthquake.

Constraints on the source parameters of low-frequency earthquakes on the San Andreas Fault

USGS Publications Warehouse

Thomas, Amanda M.; Beroza, Gregory C.; Shelly, David R.

2016-01-01

Low-frequency earthquakes (LFEs) are small repeating earthquakes that occur in conjunction with deep slow slip. Like typical earthquakes, LFEs are thought to represent shear slip on crustal faults, but when compared to earthquakes of the same magnitude, LFEs are depleted in high-frequency content and have lower corner frequencies, implying longer duration. Here we exploit this difference to estimate the duration of LFEs on the deep San Andreas Fault (SAF). We find that the M ~ 1 LFEs have typical durations of ~0.2 s. Using the annual slip rate of the deep SAF and the average number of LFEs per year, we estimate average LFE slip rates of ~0.24 mm/s. When combined with the LFE magnitude, this number implies a stress drop of ~104 Pa, 2 to 3 orders of magnitude lower than ordinary earthquakes, and a rupture velocity of 0.7 km/s, 20% of the shear wave speed. Typical earthquakes are thought to have rupture velocities of ~80–90% of the shear wave speed. Together, the slow rupture velocity, low stress drops, and slow slip velocity explain why LFEs are depleted in high-frequency content relative to ordinary earthquakes and suggest that LFE sources represent areas capable of relatively higher slip speed in deep fault zones. Additionally, changes in rheology may not be required to explain both LFEs and slow slip; the same process that governs the slip speed during slow earthquakes may also limit the rupture velocity of LFEs.

Effects of fault dip and slip rake angles on near-source ground motions: Why rupture directivity was minimal in the 1999 Chi-Chi, Taiwan, earthquake

USGS Publications Warehouse

Aagaard, Brad T.; Hall, J.F.; Heaton, T.H.

2004-01-01

We study how the fault dip and slip rake angles affect near-source ground velocities and displacements as faulting transitions from strike-slip motion on a vertical fault to thrust motion on a shallow-dipping fault. Ground motions are computed for five fault geometries with different combinations of fault dip and rake angles and common values for the fault area and the average slip. The nature of the shear-wave directivity is the key factor in determining the size and distribution of the peak velocities and displacements. Strong shear-wave directivity requires that (1) the observer is located in the direction of rupture propagation and (2) the rupture propagates parallel to the direction of the fault slip vector. We show that predominantly along-strike rupture of a thrust fault (geometry similar in the Chi-Chi earthquake) minimizes the area subjected to large-amplitude velocity pulses associated with rupture directivity, because the rupture propagates perpendicular to the slip vector; that is, the rupture propagates in the direction of a node in the shear-wave radiation pattern. In our simulations with a shallow hypocenter, the maximum peak-to-peak horizontal velocities exceed 1.5 m/sec over an area of only 200 km2 for the 30??-dipping fault (geometry similar to the Chi-Chi earthquake), whereas for the 60??- and 75??-dipping faults this velocity is exceeded over an area of 2700 km2 . These simulations indicate that the area subjected to large-amplitude long-period ground motions would be larger for events of the same size as Chi-Chi that have different styles of faulting or a deeper hypocenter.

Near-fault peak ground velocity from earthquake and laboratory data

USGS Publications Warehouse

McGarr, A.; Fletcher, Joe B.

2007-01-01

We test the hypothesis that peak ground velocity (PGV) has an upper bound independent of earthquake magnitude and that this bound is controlled primarily by the strength of the seismogenic crust. The highest PGVs, ranging up to several meters per second, have been measured at sites within a few kilometers of the causative faults. Because the database for near-fault PGV is small, we use earthquake slip models, laboratory experiments, and evidence from a mining-induced earthquake to investigate the factors influencing near-fault PGV and the nature of its scaling. For each earthquake slip model we have calculated the peak slip rates for all subfaults and then chosen the maximum of these rates as an estimate of twice the largest near-fault PGV. Nine slip models for eight earthquakes, with magnitudes ranging from 6.5 to 7.6, yielded maximum peak slip rates ranging from 2.3 to 12 m/sec with a median of 5.9 m/sec. By making several adjustments, PGVs for small earthquakes can be simulated from peak slip rates measured during laboratory stick-slip experiments. First, we adjust the PGV for differences in the state of stress (i.e., the difference between the laboratory loading stresses and those appropriate for faults at seismogenic depths). To do this, we multiply both the slip and the peak slip rate by the ratio of the effective normal stresses acting on fault planes measured at 6.8 km depth at the KTB site, Germany (deepest available in situ stress measurements), to those acting on the laboratory faults. We also adjust the seismic moment by replacing the laboratory fault with a buried circular shear crack whose radius is chosen to match the experimental unloading stiffness. An additional, less important adjustment is needed for experiments run in triaxial loading conditions. With these adjustments, peak slip rates for 10 stick-slip events, with scaled moment magnitudes from -2.9 to 1.0, range from 3.3 to 10.3 m/sec, with a median of 5.4 m/sec. Both the earthquake and laboratory results are consistent with typical maximum peak slip rates averaging between 5 and 6 m/sec or corresponding maximum near-fault PGVs between 2.5 and 3 m/sec at seismogenic depths, independent of magnitude. Our ability to replicate maximum slip rates in the fault zones of earthquakes by adjusting the corresponding laboratory rates using the ratio of effective normal stresses acting on the fault planes suggests that the strength of the seismogenic crust is the important factor limiting the near-fault PGV.

Temporal slip-rate stability and variations on the Hope Fault, New Zealand, during the late Quaternary

NASA Astrophysics Data System (ADS)

Khajavi, Narges; Nicol, Andrew; Quigley, Mark C.; Langridge, Robert M.

2018-07-01

The Hope Fault transfers slip from Hikurangi subduction to the Alpine Fault in the northern South Island of New Zealand. It accommodates mainly dextral strike slip and currently carries the highest slip rate in the Marlborough Fault System. Displacements, displacement rates and earthquake recurrence intervals have been determined using a combination of high resolution LiDAR for 59 dextral displacements ( 2.5-200 m) together with calibrated radiocarbon ages ( 130 yr to 13,000 yr) for abandoned stream channels, terrace risers and alluvial fans. Mean single-event displacement (SED) of 3 ± 0.6 m (2.2 to 4.6 m for 21 measurements) and mean recurrence interval of 266 ± 100 yr (range 128 to 560 yr) have been determined for the five most recent surface-rupturing earthquakes. On time scales ≥2300 yr the dextral slip rate is uniform at 12.2 ± 2.4 mm/yr, however, when averaged over time intervals of 230 to 1700 yr slip rates range from 4 to 46.4 mm/yr. This order-of-magnitude variability in slip rate over shorter timescales cannot be fully attributed to errors in displacement and age data, and is at least partly due to variations in earthquake recurrence interval and inferred SED. Short-term non-characteristic earthquake behaviour may be due to changes in fault loading arising from stress interactions between different segments of the Hope Fault and nearby faults.

Fault geometric complexity and how it may cause temporal slip-rate variation within an interacting fault system

NASA Astrophysics Data System (ADS)

Zielke, Olaf; Arrowsmith, Ramon

2010-05-01

Slip-rates along individual faults may differ as a function of measurement time scale. Short-term slip-rates may be higher than the long term rate and vice versa. For example, vertical slip-rates along the Wasatch Fault, Utah are 1.7+/-0.5 mm/yr since 6ka, <0.6 mm/yr since 130ka, and 0.5-0.7 mm/yr since 10Ma (Friedrich et al., 2003). Following conventional earthquake recurrence models like the characteristic earthquake model, this observation implies that the driving strain accumulation rates may have changed over the respective time scales as well. While potential explanations for such slip-rate variations may be found for example in the reorganization of plate tectonic motion or mantle flow dynamics, causing changes in the crustal velocity field over long spatial wavelengths, no single geophysical explanation exists. Temporal changes in earthquake rate (i.e., event clustering) due to elastic interactions within a complex fault system may present an alternative explanation that requires neither variations in strain accumulation rate or nor changes in fault constitutive behavior for frictional sliding. In the presented study, we explore this scenario and investigate how fault geometric complexity, fault segmentation and fault (segment) interaction affect the seismic behavior and slip-rate along individual faults while keeping tectonic stressing-rate and frictional behavior constant in time. For that, we used FIMozFric--a physics-based numerical earthquake simulator, based on Okada's (1992) formulations for internal displacements and strains due to shear and tensile faults in a half-space. Faults are divided into a large number of equal-sized fault patches which communicate via elastic interaction, allowing implementation of geometrically complex, non-planar faults. Each patch has assigned a static and dynamic friction coefficient. The difference between those values is a function of depth--corresponding to the temperature-dependence of velocity-weakening that is observed in laboratory friction experiments and expressed in an [a-b] term in Rate-State-Friction (RSF) theory. Patches in the seismic zone are incrementally loaded during the interseismic phase. An earthquake initiates if shear stress along at least one (seismic) patch exceeds its static frictional strength and may grow in size due to elastic interaction with other fault patches (static stress transfer). Aside from investigating slip-rate variations due to the elastic interactions within a fault system with this tool, we want to show how such modeling results can be very useful in exploring the physics underlying the patterns that the paleoseismology sees and that those methods (simulation and observations) can be merged, with both making important contributions. Using FIMozFric, we generated synthetic seismic records for a large number of fault geometries and structural scenarios to investigate along-fault slip accumulation patterns and the variability of slip at a point. Our simulations show that fault geometric complexity and the accompanied fault interactions and multi-fault ruptures may cause temporal deviations from the average fault slip-rate, in other words phases of earthquake clustering or relative quiescence. Slip-rates along faults within an interacting fault system may change even when the loading function (stressing rate) remains constant and the magnitude of slip rate change is suggested to be proportional to the magnitude of fault interaction. Thus, spatially isolated and structurally mature faults are expected to experience less slip-rate changes than strongly interacting and less mature faults. The magnitude of slip-rate change may serve as a proxy for the magnitude of fault interaction and vice versa.

Finding the buried memory of past earthquakes with geophysical, GPR-based paleoseismology

NASA Astrophysics Data System (ADS)

Manighetti, I.; Beaupretre, S.; Garambois, S.; Malavieille, J.; Chatton, M.; Sénéchal, G.

2011-12-01

We hypothesized that, in places where sedimentation and erosion compete at fast rates, part of the memory of past earthquakes on faults may be buried, hence hidden, in the first tens meters of the ground. We test that hypothesis on a fast slipping, large, strike-slip fault (Hope, New Zealand), at a site where marked alluvial conditions prevail (Terako). We first use LiDAR data to analyze the ground surface morphology of the 2 km2 site at the greatest resolution. About twenty clear, distinct, morphological markers are observed -mainly alluvial terrace risers and small stream channels, all are laterally offset by the fault. The measured offsets range between 3 and 200 m, yet are discrete and showing several large slip gaps. The measurements are well-constrained and allow estimating the mean slip per event amplitude to 3.9 ± 1.4 m, and the last earthquake slip to 3 ± 0.5 m. About 10 past earthquakes are well documented in the surface data, while about 50 are requested to account for the 200 m largest cumulative slip. We then investigate the zone (on smaller area, 400 x 600 m2) with dense, pseudo-3D Ground Penetrating Radar (GPR) data. We measured 56, ~400 m-long, 5-10 m spaced GPR profiles (250 MHz), parallel to the fault and evenly distributed on either sides. Their analysis reveals the existence of several tens morphological markers buried in the first 3 m of the ground, most of them do not imprint the ground surface as they are blanketed with a 0.1-3 m-thick poorly reflective layer. A few buried markers exhibit however surface expressions. All buried markers are laterally offset by the fault. Based on a number of evidence, we interpret these buried markers as stream channels, most were decapitated by the repeated fault slips and abandoned. We measured ~50 lateral offsets in the buried channel network, almost three times more than at the surface. These offsets range between 2.5 and 106 m, as observed at the surface, yet provide a more continuous record of the fault slip. The similarity of the successive slip increments suggests a slip per event averaging 3.9 ± 1.9 m, similar to that estimated from surface data. From the total 'surface and buried' offset collection, we infer that a minimum of 30 large earthquakes have broken the Hope fault at the Terako site in the last 5 kyrs, with an average slip per event of 3.8 ± 1.3 m, an average recurrence time of 100-250 yrs, and a likely magnitude of at least Mw 7.2-7.7. The last major earthquake likely occurred at 1875 ± 15 AD, in agreement with previous suggestions. Our study therefore confirms that part of the memory of past earthquakes may indeed reside in the first tens meters of the ground, where it may be explored with a novel type, geophysical and GPR-based, paleoseismology. We emphasize that developing such a new paleoseismology will provide a rich information complementary to surface observation, and help documenting the past earthquakes on faults.

Slip rates and spatially variable creep on faults of the northern San Andreas system inferred through Bayesian inversion of Global Positioning System data

USGS Publications Warehouse

Murray, Jessica R.; Minson, Sarah E.; Svarc, Jerry L.

2014-01-01

Fault creep, depending on its rate and spatial extent, is thought to reduce earthquake hazard by releasing tectonic strain aseismically. We use Bayesian inversion and a newly expanded GPS data set to infer the deep slip rates below assigned locking depths on the San Andreas, Maacama, and Bartlett Springs Faults of Northern California and, for the latter two, the spatially variable interseismic creep rate above the locking depth. We estimate deep slip rates of 21.5 ± 0.5, 13.1 ± 0.8, and 7.5 ± 0.7 mm/yr below 16 km, 9 km, and 13 km on the San Andreas, Maacama, and Bartlett Springs Faults, respectively. We infer that on average the Bartlett Springs fault creeps from the Earth's surface to 13 km depth, and below 5 km the creep rate approaches the deep slip rate. This implies that microseismicity may extend below the locking depth; however, we cannot rule out the presence of locked patches in the seismogenic zone that could generate moderate earthquakes. Our estimated Maacama creep rate, while comparable to the inferred deep slip rate at the Earth's surface, decreases with depth, implying a slip deficit exists. The Maacama deep slip rate estimate, 13.1 mm/yr, exceeds long-term geologic slip rate estimates, perhaps due to distributed off-fault strain or the presence of multiple active fault strands. While our creep rate estimates are relatively insensitive to choice of model locking depth, insufficient independent information regarding locking depths is a source of epistemic uncertainty that impacts deep slip rate estimates.

The rigid-plate and shrinking-plate hypotheses: Implications for the azimuths of transform faults

NASA Astrophysics Data System (ADS)

Mishra, Jay Kumar; Gordon, Richard G.

2016-08-01

The rigid-plate hypothesis implies that oceanic lithosphere does not contract horizontally as it cools (hereinafter "rigid plate"). An alternative hypothesis, that vertically averaged tensional thermal stress in the competent lithosphere is fully relieved by horizontal thermal contraction (hereinafter "shrinking plate"), predicts subtly different azimuths for transform faults. The size of the predicted difference is as large as 2.44° with a mean and median of 0.46° and 0.31°, respectively, and changes sign between right-lateral (RL)-slipping and left-lateral (LL)-slipping faults. For the MORVEL transform-fault data set, all six plate pairs with both RL- and LL-slipping faults differ in the predicted sense, with the observed difference averaging 1.4° ± 0.9° (95% confidence limits), which is consistent with the predicted difference of 0.9°. The sum-squared normalized misfit, r, to global transform-fault azimuths is minimized for γ = 0.8 ± 0.4 (95% confidence limits), where γ is the fractional multiple of the predicted difference in azimuth between the shrinking-plate (γ = 1) and rigid-plate (γ = 0) hypotheses. Thus, observed transform azimuths differ significantly between RL-slipping and LL-slipping faults, which is inconsistent with the rigid-plate hypothesis but consistent with the shrinking-plate hypothesis, which indicates horizontal shrinking rates of 2% Ma-1 for newly created lithosphere, 1% Ma-1 for 0.1 Ma old lithosphere, 0.2% Ma-1 for 1 Ma old lithosphere, and 0.02% Ma-1 for 10 Ma old lithosphere, which are orders of magnitude higher than the mean intraplate seismic strain rate of 10-6 Ma-1 (5 × 10-19 s-1).

Detailed fault structure of the Tarutung Pull-Apart Basin in Sumatra, Indonesia, derived from local earthquake data

NASA Astrophysics Data System (ADS)

Muksin, Umar; Haberland, Christian; Nukman, Mochamad; Bauer, Klaus; Weber, Michael

2014-12-01

The Tarutung Basin is located at a right step-over in the northern central segment of the dextral strike-slip Sumatran Fault System (SFS). Details of the fault structure along the Tarutung Basin are derived from the relocations of seismicity as well as from focal mechanism and structural geology. The seismicity distribution derived by a 3D inversion for hypocenter relocation is clustered according to a fault-like seismicity distribution. The seismicity is relocated with a double-difference technique (HYPODD) involving the waveform cross-correlations. We used 46,904 and 3191 arrival differences obtained from catalogue data and cross-correlation analysis, respectively. Focal mechanisms of events were analyzed by applying a grid search method (HASH code). Although there is no significant shift of the hypocenters (10.8 m in average) and centroids (167 m in average), the application of the double difference relocation sharpens the earthquake distribution. The earthquake lineation reflects the fault system, the extensional duplex fault system, and the negative flower structure within the Tarutung Basin. The focal mechanisms of events at the edge of the basin are dominantly of strike-slip type representing the dextral strike-slip Sumatran Fault System. The almost north-south striking normal fault events along extensional zones beneath the basin correlate with the maximum principal stress direction which is the direction of the Indo-Australian plate motion. The extensional zones form an en-echelon pattern indicated by the presence of strike-slip faults striking NE-SW to NW-SE events. The detailed characteristics of the fault system derived from the seismological study are also corroborated by structural geology at the surface.

Deformation of the 2002 Denali Fault earthquakes, mapped by Radarsat-1 interferometry

USGS Publications Warehouse

Lu, Zhong; Wright, Tim; Wicks, Chuck

2003-01-01

The magnitude 7.9 earthquake that struck central Alaska on 3 November 2002 was the largest strike-slip earthquake in North America for more than 150 years. The earthquake ruptured about 340 km of the Denali Fault system with observed right-lateral offsets of up to 9 m [Eberhart-Phillips et al., 2003] (Figure l). The rupture initiated with slip on a previously unknown thrust fault, the 40-km-long Susitna Glacier Fault. The rupture propagated eastward for about 220 km along the right-lateral Denali Fault where right-lateral slip averaged ˜5 m, before stepping southeastward onto the Totschunda Fault for about 70 km, with offsets as large as 2 m. The 3 November earthquake was preceded by a magnitude 6.7 shock on 23 October—the Nenana Mountain Earthquake—which was located about 25 km to the west of the 3 November earthquake.

Accelerating slip rates on the puente hills blind thrust fault system beneath metropolitan Los Angeles, California, USA

USGS Publications Warehouse

Bergen, Kristian J.; Shaw, John H.; Leon, Lorraine A.; Dolan, James F.; Pratt, Thomas L.; Ponti, Daniel J.; Morrow, Eric; Barrera, Wendy; Rhodes, Edward J.; Murari, Madhav K.; Owen, Lewis A.

2017-01-01

Slip rates represent the average displacement across a fault over time and are essential to estimating earthquake recurrence for proba-bilistic seismic hazard assessments. We demonstrate that the slip rate on the western segment of the Puente Hills blind thrust fault system, which is beneath downtown Los Angeles, California (USA), has accel-erated from ~0.22 mm/yr in the late Pleistocene to ~1.33 mm/yr in the Holocene. Our analysis is based on syntectonic strata derived from the Los Angeles River, which has continuously buried a fold scarp above the blind thrust. Slip on the fault beneath our field site began during the late-middle Pleistocene and progressively increased into the Holocene. This increase in rate implies that the magnitudes and/or the frequency of earthquakes on this fault segment have increased over time. This challenges the characteristic earthquake model and presents an evolving and potentially increasing seismic hazard to metropolitan Los Angeles.

Microplate model for the present-day deformation of Tibet

USGS Publications Warehouse

Thatcher, W.

2007-01-01

Site velocities from 349 Global Positioning System (GPS) stations are used to construct an 11-element quasi-rigid block model of the Tibetan Plateau and its surroundings. Rigid rotations of five major blocks are well determined, and average translation velocities of six smaller blocks can be constrained. Where data are well distributed the velocity field can be explained well by rigid block motion and fault slip across block boundaries. Residual misfits average 1.6 mm/yr compared to typical one standard deviation velocity uncertainties of 1.3 mm/yr. Any residual internal straining of the blocks is small and heterogeneous. However, residual substructure might well represent currently unresolved motions of smaller blocks. Although any smaller blocks must move at nearly the same rate as the larger blocks within which they lie, undetected relative motions between them could be significant, particularly where there are gaps in GPS coverage. Predicted relative motions between major blocks agree with the observed sense of slip and along-strike partitioning of motion across major faults. However, predicted slip rates across Tibet's major strike-slip faults are low, only 5-12 mm/yr, a factor of 2-3 smaller than most rates estimated from fault offset features dated by radiometric methods as ???2000 to ???100,000 year old. Previous work has suggested that both GPS data and low fault slip rates are incompatible with rigid block motions of Tibet. The results reported here overcome these objections.

Crustal deformation associated with east Mediterranean strike-slip earthquakes: The 8 June 2008 Movri (NW Peloponnese), Greece, earthquake (M w6.4)

NASA Astrophysics Data System (ADS)

Papadopoulos, Gerassimos A.; Karastathis, Vassilis; Kontoes, Charalambos; Charalampakis, Marinos; Fokaefs, Anna; Papoutsis, Ioannis

2010-09-01

The 2008 mainshock ( Mw = 6.4) was the first modern, strong strike-slip earthquake in the Greek mainland. The fault strikes NE-SW, dips ˜ 85°NW while the motion was right-lateral with small reverse component. Historical seismicity showed no evidence that the fault ruptured in the last 300 years. For rectangular planar fault we estimated fault dimensions from aftershock locations. Dimensions are consistent with that a buried fault was activated, lateral expansion occurred only along length and the rupture stopped at depth ˜ 20 km implying that more rupture along length was favoured. We concluded that no major asperities remained unbroken and that the aftershock activity was dominated rather by creeping mechanism than by the presence of locked patches. For Mo = 4.56 × 10 25 dyn cm we calculated average slip of 76 cm and stress drop Δσ ˜ 13 bars. This Δσ is high for Greek strike-slip earthquakes, due rather to increased rigidity because of the relatively long recurrence ( Τ > 300 years) of strong earthquakes in the fault, than to high slip. Values of Δσ and Τ indicated that the fault is neither a typical strong nor a typical weak fault. Dislocation modeling of a buried fault showed uplift of ˜ 8.0 cm in Kato Achaia ( Δ ˜ 20 km) at the hanging wall of the reverse fault component. DInSAR analysis detected co-seismic motion only in Kato Achaia where interferogram fringes pattern showed vertical displacement from 3.0 to 6.0 cm. From field-surveys we estimated maximum intensity of VIII in Kato Achaia. The most important liquefaction spots were also observed there. These observations are attributable neither to surface fault-breaks nor to site effects but possibly to high ground acceleration due to the co-seismic uplift. The causal association between displacement and earthquake damage in the hanging wall described for dip-slip faults in Taiwan, Greece and elsewhere, becomes possible also for strike-slip faults with dip-slip component, as the 2008 earthquake.

Viscoelastic shear zone model of a strike-slip earthquake cycle

USGS Publications Warehouse

Pollitz, F.F.

2001-01-01

I examine the behavior of a two-dimensional (2-D) strike-slip fault system embedded in a 1-D elastic layer (schizosphere) overlying a uniform viscoelastic half-space (plastosphere) and within the boundaries of a finite width shear zone. The viscoelastic coupling model of Savage and Prescott [1978] considers the viscoelastic response of this system, in the absence of the shear zone boundaries, to an earthquake occurring within the upper elastic layer, steady slip beneath a prescribed depth, and the superposition of the responses of multiple earthquakes with characteristic slip occurring at regular intervals. So formulated, the viscoelastic coupling model predicts that sufficiently long after initiation of the system, (1) average fault-parallel velocity at any point is the average slip rate of that side of the fault and (2) far-field velocities equal the same constant rate. Because of the sensitivity to the mechanical properties of the schizosphere-plastosphere system (i.e., elastic layer thickness, plastosphere viscosity), this model has been used to infer such properties from measurements of interseismic velocity. Such inferences exploit the predicted behavior at a known time within the earthquake cycle. By modifying the viscoelastic coupling model to satisfy the additional constraint that the absolute velocity at prescribed shear zone boundaries is constant, I find that even though the time-averaged behavior remains the same, the spatiotemporal pattern of surface deformation (particularly its temporal variation within an earthquake cycle) is markedly different from that predicted by the conventional viscoelastic coupling model. These differences are magnified as plastosphere viscosity is reduced or as the recurrence interval of periodic earthquakes is lengthened. Application to the interseismic velocity field along the Mojave section of the San Andreas fault suggests that the region behaves mechanically like a ???600-km-wide shear zone accommodating 50 mm/yr fault-parallel motion distributed between the San Andreas fault system and Eastern California Shear Zone. Copyright 2001 by the American Geophysical Union.

Rapid acceleration leads to rapid weakening in earthquake-like laboratory experiments

NASA Astrophysics Data System (ADS)

Chang, J. C.; Lockner, D. A.; Reches, Z.

2012-12-01

We simulated the slip of a fault-patch during a large earthquake by rapidly loading an experimental, ring-shaped fault with energy stored in a spinning flywheel. The flywheel abruptly delivers a finite amount of energy by spinning the fault-patch that spontaneously dissipates the energy without operator intervention. We conducted 42 experiments on Sierra White granite (SWG) samples, and 24 experiments on Kasota dolomite (KD) samples. Each experiment starts by spinning a 225 kg disk-shaped flywheel to a prescribed angular velocity. We refer to this experiment as an "earthquake-like slip-event" (ELSE). The strength-evolution in ELSE experiments is similar to the strength-evolution proposed for earthquake models and observed in stick-slip experiments. Further, we found that ELSE experiments are similar to earthquakes in at least three ways: (1) slip driven by the release of a finite amount of stored energy; (2) pattern of fault strength evolution; and (3) seismically observed values, such as average slip, peak-velocity and rise-time. By assuming that the measured slip, D, in ELSE experiments is equivalent to the average slip during an earthquake, we found that ELSE experiments (D = 0.003-4.6 m) correspond to earthquakes in moment-magnitude range of Mw = 4-8. In ELSE experiments, the critical-slip-distance, dc, has mean values of 2.7 cm and 1.2 cm for SWG and KD, that are much shorter than the 1-10 m in steady-state classical experiments in rotary shear systems. We attribute these dc values, to ELSE loading in which the fault-patch is abruptly loaded by impact with a spinning flywheel. Under this loading, the friction-velocity relations are strikingly different from those under steady-state loading on the same rock samples with the same shear system (Reches and Lockner, Nature, 2010). We further note that the slip acceleration in ELSE evolves systematically with fault strength and wear-rate, and that the dynamic weakening is restricted to the period of intense acceleration (up to 25 m/s2 during ~0.1 s). Thus, the weakening distance, dc, is reached within the initial acceleration spike. These observations are not unique, and similar weakening-acceleration associations were reported in stick-slip, rotary shear, and impact shear experiments. These studies greatly differ from each other in slip distance, normal stress, acceleration, and slip-velocities with the outstanding commonality of abrupt loading and intense acceleration. We propose that impact loading induces extremely high strain-rates that significantly increase rock brittleness, fracture tendency, and fragmentation. We envision that these processes intensify fault wear as manifested in ELSE experiments by extremely high initial wear-rates. This intense, early wear generates a layer of fine-grain gouge that reduces the fault strength by powder-lubrication. Our analysis indicates that rapid acceleration associated with earthquake rupture accelerates fault weakening and shortens the weakening-distance.

Apparent late Quaternary fault slip rate increase in the southwestern Lower Rhine Graben, central Europe

USGS Publications Warehouse

Gold, Ryan D.; Friedrich, Anke M.; Kubler, Simon; Salamon, Martin

2017-01-01

In regions of low strain, long earthquake recurrence intervals (104–106  yrs) and erosive processes limit preservation of Quaternary markers suitable for distinguishing whether faults slip at uniform or secularly varying rates. The Lower Rhine graben in the border region of Germany, The Netherlands, and Belgium provides a unique opportunity to explore Quaternary slip‐rate variations in a region of low strain using the basal (2.29±0.29  Ma) and surface (700±80  ka) contacts of the regionally extensive main terrace (“Hauptterrasse”), deposited by the Rhine and Maas Rivers. These surfaces are vertically offset 3–140 m and 0–68 m, respectively, across individual fault strands within a distributed network of northwest‐trending, slow‐slipping (<0.1  mm/yr) normal faults. In this investigation, we construct Quaternary slip histories for the southern Lower Rhine graben faults using new main terrace surface vertical offset measurements made from light detection and ranging (lidar)‐derived bare‐earth digital terrain models, which we synthesize with existing constraints on the offset basal contact of this fluvial deposit (n=91 collocated sites with displacement constraints). We find that >80% of the sites record an apparent increase in slip rate for the more recent interval from 700 ka to present, which corresponds to a period of increased uplift of the nearby Rhenish Massif and regional volcanism. However, the apparent increase in slip rate could result, in part, from erosion of the footwall surface below the main terrace, leading to an apparent displacement that is smaller than the total vertical offset since the start of the Quaternary. Prior work focused on characterization of these faults as seismic sources in the Lower Rhine graben has preferentially relied on the average fault‐slip rate constrained using the base of the main terrace. We suggest that average fault‐slip rates calculated using the ∼700  ka main terrace surface are subjected to fewer uncertainties and sample a time interval that is more relevant for seismic‐hazard analysis.

Secular Variation in Slip (Invited)

NASA Astrophysics Data System (ADS)

Cowgill, E.; Gold, R. D.

2010-12-01

Faults show temporal variations in slip rate at time scales ranging from the hours following a major rupture to the millions of years over which plate boundaries reorganize. One such behavior is secular variation in slip (SVS), which we define as a pulse of accelerated strain release along a single fault that occurs at a frequency that is > 1 order of magnitude longer than the recurrence interval of earthquakes within the pulse. Although numerous mechanical models have been proposed to explain SVS, it has proven much harder to measure long (5-500 kyr) records of fault displacement as a function of time. Such fault-slip histories may be obtained from morphochronologic data, which are measurements of offset and age obtained from faulted landforms. Here we describe slip-history modeling of morphochronologic data and show how this method holds promise for obtaining long records of fault slip. In detail we place SVS in the context of other types of time-varying fault-slip phenomena, explain the importance of measuring fault-slip histories, summarize models proposed to explain SVS, review current approaches for measuring SVS in the geologic record, and illustrate the slip-history modeling approach we advocate here using data from the active, left-slip Altyn Tagh fault in NW Tibet. In addition to SVS, other types of temporal variation in fault slip include post-seismic transients, discrepancies between geologic slip rates and those derived from geodetic and/or paleoseismic data, and single changes in slip rate resulting from plate reorganization. Investigating secular variation in slip is important for advancing understanding of long-term continental deformation, fault mechanics, and seismic risk. Mechanical models producing such behavior include self-driven mode switching, changes in pore-fluid pressure, viscoelasticity, postseismic reloading, and changes in local surface loads (e.g., ice sheets, large lakes, etc.) among others. However, a key problem in testing these models is the paucity of long records of fault slip. Paleoseismic data are unlikely to yield such histories because measurements of the slip associated with each event are generally unavailable and long records require large accumulated offsets, which can result in structural duplication or omission of the stratigraphic records of events. In contrast, morphochronologic data capture both the age and offset of individual piercing points, although this approach generally does not resolve individual earthquake events. Because the uncertainties in both age and offset are generally large (5-15%) for individual markers, SVS is best resolved by obtaining suites of such measurements, in which case the errors can be used to reduce the range of slip histories common to all such data points. A suite of such data from the central Altyn Tagh fault reveals a pulse of accelerated strain release in the mid Holocene, with ~20 m of slip being released from ~6.7 to ~5.9 ka at a short-term rate (~28 mm/yr) that is 3 times greater than the average rate (~9 mm/yr). We interpret this pulse to represent a cluster of two to six, Mw > 7.2 earthquakes. To our knowledge, this is the first possible earthquake cluster detected using morphochronologic techniques.

Seven big strike-slip earthquakes

NASA Astrophysics Data System (ADS)

Lohman, R. B.; Simons, M.; Pritchard, M. E.

2003-12-01

We examine seven large (Mw > 7) strike-slip earthquakes that occurred since the beginning of ERS 1 and 2 missions. We invert GPS observations and InSAR interferograms and azimuth offsets for coseismic slip distributions. We explore two refinements to the traditional least-squares inversion technique with roughness constraints. First, we diverge from the usual definition of ``roughness'' as the average roughness over the entire fault plane, and allow ``variable smoothing'' constraints. Variable smoothing allows our inversion to select models that are more complex in regions that are well-resolved by the data, while still damping regions that are poorly resolved. Second, we choose our smoothing parameters using the jR_i criterion. The jR_i criterion draws on the theory behind cross-validation and the bootstrap method. We examine the theoretical basis behind such methods and use an analytical approximation technique for linear problems. We provide maps of model variance and spatial averaging scale over the fault plane, to explicitly show which features in our slip models are robust. We examine the 1992 Landers (CA), 1995 Sakhalin (Russia), 1995 Kobe (Japan), 1997 Ardekul (Iran), 1997 Manyi (Tibet), 1999 Hector Mine (CA), and 2001 Kunlun (Tibet) earthquakes. We compare features of the slip distributions such as the depth distribution of slip, the inferred magnitude and the degree of heterogeneity of slip over the fault plane, as resolved by the available InSAR and GPS data. We end with a brief description of the data coverage required for future earthquakes of similar size if we want to infer some of the above quantities to within a given confidence interval. We describe both the number of InSAR scenes and the distribution of GPS points that would be required, based on theoretical treatments of the fault plane/data point geometry using the jR_i method.

Source parameters of the 2016 Menyuan earthquake in the northeastern Tibetan Plateau determined from regional seismic waveforms and InSAR measurements

NASA Astrophysics Data System (ADS)

Liu, Yunhua; Zhang, Guohong; Zhang, Yingfeng; Shan, Xinjian

2018-06-01

On January 21st, 2016, a Ms 6.4 earthquake hit Menyuan County, Qinghai province, China. The nearest known fault is the Leng Long Ling (LLL) fault which is located approximately 7 km north of the epicenter. This fault has mainly shown sinistral strike-slip movement since the late Quaternary Period. However, the focal mechanism indicates that it is a thrust earthquake, which is different from the well-known strike-slip feature of the LLL fault. In this study, we determined the focal mechanism and primary nodal plane through multi-step inversions in the frequency and time domain by using the broadband regional seismic waveforms recorded by the China Digital Seismic Network (CDSN). Our results show that the rupture duration was short, within 0-2 s after the earthquake, and the rupture expanded upwards along the fault plane. Based on these fault parameters, we then solve for variable slip distribution on the fault plane using the InSAR data. We applied a three-segment fault model to simulate the arc-shaped structure of the northern LLL fault, and obtained a detailed slip distribution on the fault plane. The inversion results show that the maximum slip is 0.43 m, and the average slip angle is 78.8°, with a magnitude of Mw 6.0 and a focal depth of 9.38 km. With the geological structure and the inversion results taken into consideration, it can be suggested that this earthquake was caused by the arc-shaped secondary fault located at the north side of the LLL fault. The secondary fault, together with the LLL fault, forms a normal flower structure. The main LLL fault extends almost vertically into the base rock and the rocks between the two faults form a bulging fault block. Therefore, we infer that this earthquake is the manifestation of a neotectonics movement, in which the bulging fault block is lifted further up under the compresso-shear action caused by the Tibetan Plateau pushing towards the northwest direction.

Strength variation along the Altyn Tagh and the Kunlun fault, northern Tibetan plateau, inferred from 3D mechanical modeling

NASA Astrophysics Data System (ADS)

Zhu, X.; He, J.; Xiao, J.

2017-12-01

The Altyn Tagh (ATF) and the Kunlun (KLF) fault distribute around the northern Tibetan plateau from west to east about 2000 km and 1200 km in length, and deform predominately with left-lateral strike-slip motion. Previous geological and geodetic observations suggested that over at least 800-km length of the faults, the slip rate averaged on active deformation period is quite uniform, for the ATF being about 9-10 mm/yr and the KLF about 10-12mm/yr. Strike-slip deformation of these faults is undoubtedly result from regional loading by ongoing collision between the India and the Eurasia continent. Whereas, dense GPS measurements show that along the central Tibetan plateau from west to east, the GPS velocity field changes greatly both on magnitude and on direction, suggesting that tectonic loading to the ATF and the KLF could be changed along their strike directions. To investigate how a non-uniform tectonic loading condition as documented by the GPS velocity field could cause a relatively uniform slip rate of the two active faults, we built a three-dimensional viscoelastic finite element model, in which motion of the strike-slip fault is governed by frictional strength. Given a reasonable bound of model parameters, we at first test the numerical calculation with uniform frictional coefficient of the faults. At this condition, the predicted slip rate is inevitably largest near center of the faults and gradually decreasing to the fault ends. To better fitting the observed uniform slip rate along the faults over 1000km length, variation of fault strength along the ATF and the KLF must be invoked. By testing numerous models, an optimum result was obtained, among which the frictional coefficient for the ATF is varied from 0.02 to 0.12 between 820E and 1000E with its maximum at 840E, and for the KLF from 0.02 to 0.10 with its maximum between 950E and 970E. This means that the strength of the two large-scale strike-slip faults exists significant difference along their strikes. We believe that the predicted fault pattern could play an important role on partitioning strain aside the fault, together on determination of potential rupture during an earthquake.

Surface Rupture and Slip Distribution Resulting from the 2013 M7.7 Balochistan, Pakistan Earthquake

NASA Astrophysics Data System (ADS)

Reitman, N. G.; Gold, R. D.; Briggs, R. W.; Barnhart, W. D.; Hayes, G. P.

2014-12-01

The 24 September 2013 M7.7 earthquake in Balochistan, Pakistan, produced a ~200 km long left-lateral strike-slip surface rupture along a portion of the Hoshab fault, a moderately dipping (45-75º) structure in the Makran accretionary prism. The rupture is remarkably continuous and crosses only two (0.7 and 1.5 km wide) step-overs along its arcuate path through southern Pakistan. Displacements are dominantly strike-slip, with a minor component of reverse motion. We remotely mapped the surface rupture at 1:5,000 scale and measured displacements using high resolution (0.5 m) pre- and post-event satellite imagery. We mapped 295 laterally faulted stream channels, terrace margins, and roads to quantify near-field displacement proximal (±10 m) to the rupture trace. The maximum near-field left-lateral offset is 15±2 m (average of ~7 m). Additionally, we used pre-event imagery to digitize 254 unique landforms in the "medium-field" (~100-200 m from the rupture) and then measured their displacements compared to the post-event imagery. At this scale, maximum left-lateral offset approaches 17 m (average of ~8.5 m). The width (extent of observed surface faulting) of the rupture zone varies from ~1 m to 3.7 km. Near- and medium-field offsets show similar slip distributions that are inversely correlated with the width of the fault zone at the surface (larger offsets correspond to narrow fault zones). The medium-field offset is usually greater than the near-field offset. The along-strike surface slip distribution is highly variable, similar to the slip distributions documented for the 2002 Denali M7.9 earthquake and 2001 Kunlun M7.8 earthquake, although the Pakistan offsets are larger in magnitude. The 2013 Pakistan earthquake ranks among the largest documented continental strike-slip displacements, possibly second only to the 18+ m surface displacements attributed to the 1855 Wairarapa M~8.1 earthquake.

Slip rate and slip magnitudes of past earthquakes along the Bogd left-lateral strike-slip fault (Mongolia)

USGS Publications Warehouse

Rizza, M.; Ritz, J.-F.; Braucher, R.; Vassallo, R.; Prentice, C.; Mahan, S.; McGill, S.; Chauvet, A.; Marco, S.; Todbileg, M.; Demberel, S.; Bourles, D.

2011-01-01

We carried out morphotectonic studies along the left-lateral strike-slip Bogd Fault, the principal structure involved in the Gobi-Altay earthquake of 1957 December 4 (published magnitudes range from 7.8 to 8.3). The Bogd Fault is 260 km long and can be subdivided into five main geometric segments, based on variation in strike direction. West to East these segments are, respectively: the West Ih Bogd (WIB), The North Ih Bogd (NIB), the West Ih Bogd (WIB), the West Baga Bogd (WBB) and the East Baga Bogd (EBB) segments. Morphological analysis of offset streams, ridges and alluvial fans-particularly well preserved in the arid environment of the Gobi region-allows evaluation of late Quaternary slip rates along the different faults segments. In this paper, we measure slip rates over the past 200 ka at four sites distributed across the three western segments of the Bogd Fault. Our results show that the left-lateral slip rate is ~1 mm yr-1 along the WIB and EIB segments and ~0.5 mm yr-1 along the NIB segment. These variations are consistent with the restraining bend geometry of the Bogd Fault. Our study also provides additional estimates of the horizontal offset associated with the 1957 earthquake along the western part of the Bogd rupture, complementing previously published studies. We show that the mean horizontal offset associated with the 1957 earthquake decreases progressively from 5.2 m in the west to 2.0 m in the east, reflecting the progressive change of kinematic style from pure left-lateral strike-slip faulting to left-lateral-reverse faulting. Along the three western segments, we measure cumulative displacements that are multiples of the 1957 coseismic offset, which may be consistent with a characteristic slip. Moreover, using these data, we re-estimate the moment magnitude of the Gobi-Altay earthquake at Mw 7.78-7.95. Combining our slip rate estimates and the slip distribution per event we also determined a mean recurrence interval of ~2500-5200 yr for past earthquakes along the different segments of the western Bogd Fault. This suggests that the three western segments of the Bogd Fault and the Gurvan Bulag thrust fault (a reverse fault bounding the southern side of the Ih Bogd range that ruptured during the 1957 earthquake) have similar average recurrence times, and therefore may have ruptured together in previous earthquakes as they did in 1957. These results suggest that the western part of the Bogd Fault system, including the Gurvan Bulag thrust fault, usually behaves in a 'characteristic earthquake' mode. ?? 2011 The Authors Geophysical Journal International ?? 2011 RAS.

Slip rate and slip magnitudes of past earthquakes along the Bogd left-lateral strike-slip fault (Mongolia)

USGS Publications Warehouse

Prentice, Carol S.; Rizza, M.; Ritz, J.F.; Baucher, R.; Vassallo, R.; Mahan, S.

2011-01-01

We carried out morphotectonic studies along the left-lateral strike-slip Bogd Fault, the principal structure involved in the Gobi-Altay earthquake of 1957 December 4 (published magnitudes range from 7.8 to 8.3). The Bogd Fault is 260 km long and can be subdivided into five main geometric segments, based on variation in strike direction. West to East these segments are, respectively: the West Ih Bogd (WIB), The North Ih Bogd (NIB), the West Ih Bogd (WIB), the West Baga Bogd (WBB) and the East Baga Bogd (EBB) segments. Morphological analysis of offset streams, ridges and alluvial fans—particularly well preserved in the arid environment of the Gobi region—allows evaluation of late Quaternary slip rates along the different faults segments. In this paper, we measure slip rates over the past 200 ka at four sites distributed across the three western segments of the Bogd Fault. Our results show that the left-lateral slip rate is∼1 mm yr–1 along the WIB and EIB segments and∼0.5 mm yr–1 along the NIB segment. These variations are consistent with the restraining bend geometry of the Bogd Fault. Our study also provides additional estimates of the horizontal offset associated with the 1957 earthquake along the western part of the Bogd rupture, complementing previously published studies. We show that the mean horizontal offset associated with the 1957 earthquake decreases progressively from 5.2 m in the west to 2.0 m in the east, reflecting the progressive change of kinematic style from pure left-lateral strike-slip faulting to left-lateral-reverse faulting. Along the three western segments, we measure cumulative displacements that are multiples of the 1957 coseismic offset, which may be consistent with a characteristic slip. Moreover, using these data, we re-estimate the moment magnitude of the Gobi-Altay earthquake at Mw 7.78–7.95. Combining our slip rate estimates and the slip distribution per event we also determined a mean recurrence interval of∼2500–5200 yr for past earthquakes along the different segments of the western Bogd Fault. This suggests that the three western segments of the Bogd Fault and the Gurvan Bulag thrust fault (a reverse fault bounding the southern side of the Ih Bogd range that ruptured during the 1957 earthquake) have similar average recurrence times, and therefore may have ruptured together in previous earthquakes as they did in 1957. These results suggest that the western part of the Bogd Fault system, including the Gurvan Bulag thrust fault, usually behaves in a ‘characteristic earthquake’ mode.

New slip rate estimates for the Mission Creek strand of the San Andreas fault zone

NASA Astrophysics Data System (ADS)

Blisniuk, K.; Scharer, K. M.; Sharp, W. D.; Burgmann, R.; Rymer, M. J.; Williams, P. L.

2013-12-01

The potential for a large-magnitude earthquake (Mw ≥ 6.7) on the southern San Andreas fault zone (SAFZ) is generally considered high (Working Group on California Earthquake Probabilities, 2007). However, the proportion of slip accommodated by each of its three major fault strands (Mission Creek, Banning, and Garnet Hill, from north to south) in the Indio Hills is poorly constrained. Each of these strands cut through San Gorgonio Pass west to the Los Angeles metropolitan region. To better assess the relative importance of these faults and their potential for a major earthquake, we dated offsets at two sites on the Mission Creek fault in the central Indio Hills, an offset channel at Pushawalla Canyon and an offset debris cone at a small unnamed canyon located ~1.5 km farther southeast. Previous work on this strand at Biskra Palms, in the southern Indio Hills, demonstrated a slip rate between 12 and 22 mm/yr, with a preferred rate of 14-17 mm/yr (Behr et al., GSAB, 2010). It is generally assumed that the slip rate on the Mission Creek fault decreases northwestwards from Biskra Palms (e.g. Fumal et al., BSSA, 2002) towards these two sites in the central Indio Hills. However, our initial results from uranium-series dating of pedogenic carbonate and 10Be cosmogenic exposure dating of surface clasts from deposits offset 1.3-1.6 km since ~70 ka and 44-50 m since ~2.5 ka indicate that during the late Pleistocene and Holocene slip on the Mission Creek fault in the central Indio Hills has occurred at a relatively constant and unexpectedly high rate of ~20 mm/yr. Combined with published paleoseismic studies for the Mission Creek fault, which show an average earthquake recurrence interval of 225 years for the past 5 events since 900 AD (Fumal et al., 2002), these data imply an average slip-per-event of ~4.5 m. The last earthquake to rupture this section of the Mission Creek fault occurred over 300 years ago (ca. 1690), which indicates that ca. 5.0 to 7.5 m of strain may have accumulated since the last surface-rupturing event. While additional work is needed to better understand how slip along the SAFZ is partitioned in the northwestern Indio Hills, the new data underscore the seismic hazard posed by the Mission Creek fault in this region.

Constraining the slip distribution and fault geometry of the Mw 7.9, 3 November 2002, Denali fault earthquake with Interferometric Synthetic Aperture Radar and Global Positioning System data

USGS Publications Warehouse

Wright, Tim J.; Lu, Z.; Wicks, Charles

2004-01-01

The Mw 7.9, Denali fault earthquake (DFE) is the largest continental strike-slip earthquake to occur since the development of Interferometric Synthetic Aperture Radar (InSAR). We use five interferograms, constructed using radar images from the Canadian Radarsat-1 satellite, to map the surface deformation at the western end of the fault rupture. Additional geodetic data are provided by displacements observed at 40 campaign and continuous Global Positioning System (GPS) sites. We use the data to determine the geometry of the Susitna Glacier fault, thrusting on which initiated the DFE, and to determine a slip model for the entire event that is consistent with both the InSAR and GPS data. We find there was an average of 7.3 ± 0.4 m slip on the Susitna Glacier fault, between 1 and 9.5 km depth on a 29 km long fault that dips north at 41 ± 0.7° and has a surface projection close to the mapped rupture. On the Denali fault, a simple model with large slip patches finds a maximum of 8.7 ± 0.7 m of slip between the surface and 14.3 ± 0.2 km depth. A more complex distributed slip model finds a peak of 12.5 ± 0.8 m in the upper 4 km, significantly higher than the observed surface slip. We estimate a geodetic moment of 670 ± 10 × 1018 N m (Mw 7.9), consistent with seismic estimates. Lack of preseismic data resulted in an absence of InSAR coverage for the eastern half of the DFE rupture. A dedicated geodetic InSAR mission could obviate coverage problems in the future.

Southeastward increase of the late Quaternary slip-rate of the Xianshuihe fault, eastern Tibet. Geodynamic and seismic hazard implications

NASA Astrophysics Data System (ADS)

Bai, Mingkun; Chevalier, Marie-Luce; Pan, Jiawei; Replumaz, Anne; Leloup, Philippe Hervé; Métois, Marianne; Li, Haibing

2018-03-01

The left-lateral strike-slip Xianshuihe fault system located in the eastern Tibetan Plateau is considered as one of the most tectonically active intra-continental fault system in China, along which more than 20 M > 6.5 and more than 10 M > 7 earthquakes occurred since 1700. Therefore, studying its activity, especially its slip rate at different time scales, is essential to evaluate the regional earthquake hazard. Here, we focus on the central segment of the Xianshuihe fault system, where the Xianshuihe fault near Kangding city splays into three branches: the Selaha, Yalahe and Zheduotang faults. In this paper we use precise dating together with precise field measurements of offsets to re-estimate the slip rate of the fault that was suggested without precise age constraints. We studied three sites where the active Selaha fault cuts and left-laterally offsets moraine crests and levees. We measured horizontal offsets of 96 ± 20 m at Tagong levees (TG), 240 ± 15 m at Selaha moraine (SLH) and 80 ± 5 m at Yangjiagou moraine (YJG). Using 10Be cosmogenic dating, we determined abandonment ages at Tagong, Selaha and Yangjiagou of 12.5 (+ 2.5 / - 2.2) ka, 22 ± 2 ka, and 18 ± 2 ka, respectively. By matching the emplacement age of the moraines or levees with their offsets, we obtain late Quaternary horizontal average slip-rates of 7.6 (+ 2.3 / - 1.9) mm/yr at TG and 10.7 (+ 1.3 / - 1.1) mm/yr at SLH, i.e., 5.7-12 mm/yr or between 9.6 and 9.9 mm/yr assuming that the slip rate should be constant between the nearby TG and SLH sites. At YJG, we obtain a lower slip rate of 4.4 ± 0.5 mm/yr, most likely because the parallel Zheduotang fault shares the slip rate at this longitude, therefore suggesting a ∼5 mm/yr slip rate along the Zheduotang fault. The ∼10 mm/yr late Quaternary rate along the Xianshuihe fault is higher than that along the Ganzi fault to the NW (6-8 mm/yr). This appears to be linked to the existence of the Longriba fault system that separates the Longmenshan and Bayan Har blocks north of the Xianshuihe fault system. A higher slip rate along the short (∼60 km) and discontinuous Selaha fault compared to that along the long (∼300 km) and linear Ganzi fault suggests a high hazard for a M > 6 earthquake in the Kangding area in the near future, which could devastate that densely populated city.

The 12 June 2017 Mw 6.3 Lesvos Island (Aegean Sea) earthquake: Slip model and directivity estimated with finite-fault inversion

NASA Astrophysics Data System (ADS)

Kiratzi, Anastasia

2018-01-01

On 12 June 2017 (UTC 12:28:38.26) a magnitude Mw 6.3 earthquake occurred offshore Lesvos Island in SE Aegean Sea, which was widely felt, caused 1 fatality, and partially ruined the village of Vrisa on the south-eastern coast of the island. I invert broad band and strong motion waveforms from regional stations to obtain the source model and the distribution of slip onto the fault plane. The hypocentre is located at a depth of 7 km in the upper crust. The mainshock ruptured a WNW-ESE striking, SW dipping, normal fault, projecting offshore and bounding the Lesvos Basin. The strongest and most aftershocks clustered away from the hypocentre, at the eastern edge of the activated area. This cluster indicates the activation of a different fault segment, exhibiting sinistral strike-slip motions, along a plane striking WNW-ESE. The slip of the mainshock is confined in a single large asperity, WNW from the hypocentre, with dimensions 20 km × 10 km along fault strike and dip, respectively. The average slip of the asperity is 50 cm and the peak slip is 1 m. The rupture propagated unilaterally towards WNW to the coastline of Lesvos island at a relatively high speed ( 3.1 km/s). The imaged slip model and forward modelling was used to calculate peak ground velocities (PGVs) in the near-field. The damage pattern produced by this earthquake, especially in the village of Vrisa is compatible with the combined effect of rupture directivity, proximity to the slip patch and the fault edge, spectral content of motions, and local site conditions.

A low-angle detachment fault revealed: Three-dimensional images of the S-reflector fault zone along the Galicia passive margin

NASA Astrophysics Data System (ADS)

Schuba, C. Nur; Gray, Gary G.; Morgan, Julia K.; Sawyer, Dale S.; Shillington, Donna J.; Reston, Tim J.; Bull, Jonathan M.; Jordan, Brian E.

2018-06-01

A new 3-D seismic reflection volume over the Galicia margin continent-ocean transition zone provides an unprecedented view of the prominent S-reflector detachment fault that underlies the outer part of the margin. This volume images the fault's structure from breakaway to termination. The filtered time-structure map of the S-reflector shows coherent corrugations parallel to the expected paleo-extension directions with an average azimuth of 107°. These corrugations maintain their orientations, wavelengths and amplitudes where overlying faults sole into the S-reflector, suggesting that the parts of the detachment fault containing multiple crustal blocks may have slipped as discrete units during its late stages. Another interface above the S-reflector, here named S‧, is identified and interpreted as the upper boundary of the fault zone associated with the detachment fault. This layer, named the S-interval, thickens by tens of meters from SE to NW in the direction of transport. Localized thick accumulations also occur near overlying fault intersections, suggesting either non-uniform fault rock production, or redistribution of fault rock during slip. These observations have important implications for understanding how detachment faults form and evolve over time. 3-D seismic reflection imaging has enabled unique insights into fault slip history, fault rock production and redistribution.

Earthquake geology and paleoseismology of major strands of the San Andreas fault system: Chapter 38

USGS Publications Warehouse

Rockwell, Thomas; Scharer, Katherine M.; Dawson, Timothy E.

2016-01-01

The San Andreas fault system in California is one of the best-studied faults in the world, both in terms of the long-term geologic history and paleoseismic study of past surface ruptures. In this paper, we focus on the Quaternary to historic data that have been collected from the major strands of the San Andreas fault system, both on the San Andreas Fault itself, and the major subparallel strands that comprise the plate boundary, including the Calaveras-Hayward- Rogers Creek-Maacama fault zone and the Concord-Green Valley-Bartlett Springs fault zone in northern California, and the San Jacinto and Elsinore faults in southern California. The majority of the relative motion between the Pacific and North American lithospheric plates is accommodated by these faults, with the San Andreas slipping at about 34 mm/yr in central California, decreasing to about 20 mm/yr in northern California north of its juncture with the Calaveras and Concord faults. The Calaveras-Hayward-Rogers Creek-Maacama fault zone exhibits a slip rate of 10-15 mm/yr, whereas the rate along the Concord-Green Valley-Bartlett Springs fault zone is lower at about 5 mm/yr. In southern California, the San Andreas exhibits a slip rate of about 35 mm/yr along the Mojave section, decreasing to as low as 10-15 mm/yr along its juncture with the San Jacinto fault, and about 20 mm/yr in the Coachella Valley. The San Jacinto and Elsinore fault zones exhibit rates of about 15 and 5 mm/yr, respectively. The average recurrence interval for surface-rupturing earthquakes along individual elements of the San Andreas fault system range from 100-500 years and is consistent with slip rate at those sites: higher slip rates produce more frequent or larger earthquakes. There is also evidence of short-term variations in strain release (slip rate) along various fault sections, as expressed as “flurries” or clusters of earthquakes as well as periods of relatively fewer surface ruptures in these relatively short records. This is reflected by non-periodic coefficients of variation in earthquake recurrence of 0.4 to 0.7 for the various paleoseismic sites.

Spatial and Temporal Variations in Slip Partitioning During Oblique Convergence Experiments

NASA Astrophysics Data System (ADS)

Beyer, J. L.; Cooke, M. L.; Toeneboehn, K.

2017-12-01

Physical experiments of oblique convergence in wet kaolin demonstrate the development of slip partitioning, where two faults accommodate strain via different slip vectors. In these experiments, the second fault forms after the development of the first fault. As one strain component is relieved by one fault, the local stress field then favors the development of a second fault with different slip sense. A suite of physical experiments reveals three styles of slip partitioning development controlled by the convergence angle and presence of a pre-existing fault. In experiments with low convergence angles, strike-slip faults grow prior to reverse faults (Type 1) regardless of whether the fault is precut or not. In experiments with moderate convergence angles, slip partitioning is dominantly controlled by the presence of a pre-existing fault. In all experiments, the primarily reverse fault forms first. Slip partitioning then develops with the initiation of strike-slip along the precut fault (Type 2) or growth of a secondary reverse fault where the first fault is steepest. Subsequently, the slip on the first fault transitions to primarily strike-slip (Type 3). Slip rates and rakes along the slip partitioned faults for both precut and uncut experiments vary temporally, suggesting that faults in these slip-partitioned systems are constantly adapting to the conditions produced by slip along nearby faults in the system. While physical experiments show the evolution of slip partitioning, numerical simulations of the experiments provide information about both the stress and strain fields, which can be used to compute the full work budget, providing insight into the mechanisms that drive slip partitioning. Preliminary simulations of precut experiments show that strain energy density (internal work) can be used to predict fault growth, highlighting where fault growth can reduce off-fault deformation in the physical experiments. In numerical simulations of uncut experiments with a first non-planar oblique slip fault, strain energy density is greatest where the first fault is steepest, as less convergence is accommodated along this portion of the fault. The addition of a second slip-partitioning fault to the system decreases external work indicating that these faults increase the mechanical efficiency of the system.

Trench investigation along the Merida section of the Bocono fault (central Venezuelan Andes), Venezuela

USGS Publications Warehouse

Audemard, F.; Pantosti, D.; Machette, M.; Costa, C.; Okumura, K.; Cowan, H.; Diederix, H.; Ferrer, C.

1999-01-01

The Bocono fault is a major NE-SW-trending, dextral fault that extends for about 500 km along the backbone of the Venezuelan Andes. Several large historical earthquakes in this region have been attributed to the Bocono fault, and some of these have been recently associated with specific parts through paleoseismologic investigations. A new trench study has been performed, 60 km to the northeast of Merida in the central Venezuelan Andes, where the fault forms a releasing bend, comprising two conspicuous late Holocene fault strands that are about 1 km apart. The southern and northern strands carry about 70% and 30% (respectively) of the 7-10 mm/yr net slip rate measured in this sector, which is based on a 40 vs. 85-100 m right-lateral offset of the Late Pleistocene Los Zerpa moraines. A trench excavated on the northern strand of the fault (near Morros de los Hoyos, slightly northeast of Apartaderos) across a twin shutter ridge and related sag pond exposed two main fault zones cutting Late Pleistocene alluvial and Holocene peat deposits. Each zone forms a shutter ridge with peat deposits ponded against the uplifted block. The paleoearthquake reconstruction derived from this trench allow us to propose the occurrence of at least 6-8 earthquakes in the past 9000 yr, yielding a maximum average recurrence interval of about 1100-1500 yr. Based on the northern strands average slip rate (2.6 mm/yr), such as earthquake sequence should have accommodated about 23 m of slip since 9 ka, suggesting that the maximum slip per event ranges between 3 and 4 m. No direct evidence for the large 1812 earthquake has been found in the trench, although this earthquake may have ruptured this section of the fault. Further paleoseismic studies will investigate the possibility that this event occurred in the Bocono fault, but ruptured mainly its southern strand in this region.

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

Slip on the San Andreas fault at Parkfield, California, over two earthquake cycles, and the implications for seismic hazard

USGS Publications Warehouse

Murray, J.; Langbein, J.

2006-01-01

Parkfield, California, which experienced M 6.0 earthquakes in 1934, 1966, and 2004, is one of the few locales for which geodetic observations span multiple earthquake cycles. We undertake a comprehensive study of deformation over the most recent earthquake cycle and explore the results in the context of geodetic data collected prior to the 1966 event. Through joint inversion of the variety of Parkfield geodetic measurements (trilateration, two-color laser, and Global Positioning System), including previously unpublished two-color data, we estimate the spatial distribution of slip and slip rate along the San Andreas using a fault geometry based on precisely relocated seismicity. Although the three most recent Parkfield earthquakes appear complementary in their along-strike distributions of slip, they do not produce uniform strain release along strike over multiple seismic cycles. Since the 1934 earthquake, more than 1 m of slip deficit has accumulated on portions of the fault that slipped in the 1966 and 2004 earthquakes, and an average of 2 m of slip deficit exists on the 33 km of the fault southeast of Gold Hill to be released in a future, perhaps larger, earthquake. It appears that the fault is capable of partially releasing stored strain in moderate earthquakes, maintaining a disequilibrium through multiple earthquake cycles. This complicates the application of simple earthquake recurrence models that assume only the strain accumulated since the most recent event is relevant to the size or timing of an upcoming earthquake. Our findings further emphasize that accumulated slip deficit is not sufficient for earthquake nucleation.

Holocene slip rates along the San Andreas Fault System in the San Gorgonio Pass and implications for large earthquakes in southern California

NASA Astrophysics Data System (ADS)

Heermance, Richard V.; Yule, Doug

2017-06-01

The San Gorgonio Pass (SGP) in southern California contains a 40 km long region of structural complexity where the San Andreas Fault (SAF) bifurcates into a series of oblique-slip faults with unknown slip history. We combine new 10Be exposure ages (Qt4: 8600 (+2100, -2200) and Qt3: 5700 (+1400, -1900) years B.P.) and a radiocarbon age (1260 ± 60 years B.P.) from late Holocene terraces with scarp displacement of these surfaces to document a Holocene slip rate of 5.7 (+2.7, -1.5) mm/yr combined across two faults. Our preferred slip rate is 37-49% of the average slip rates along the SAF outside the SGP (i.e., Coachella Valley and San Bernardino sections) and implies that strain is transferred off the SAF in this area. Earthquakes here most likely occur in very large, throughgoing SAF events at a lower recurrence than elsewhere on the SAF, so that only approximately one third of SAF ruptures penetrate or originate in the pass.