J. L. Svarc

Postseismic deformation associated with the 1992 M ω=7.3 Landers earthquake, southern California

Following the 1992 Mω=7.3 Landers earthquake, a linear array of 10 geodetic monuments at roughly 5-km spacing was established across the Emerson fault segment of the Landers rupture. The array trends perpendicular to the local strike of the fault segment and extends about 30 km on either side of it. The array was surveyed by Global Positioning System 0.034, 0.048, 0.381, 1.27, 1.88, 2.60, and 3.42 years after the Landers earthquake to measure both the spatial and temporal character of the postearthquake relaxation. The temporal behavior is described roughly by a short-term (decay time 84±23 days) exponential relaxation superimposed upon an apparently linear trend. Because the linear trend represents motions much more rapid than the observed preseismic motions, we attribute that trend to a slower (decay time greater than 5 years) postseismic relaxation, the curvature of which cannot be resolved in the short run (3.4 years) of postseismic data. About 100 mm of right-lateral displacement and 50 mm of fault-normal displacement accumulated across the geodetic array in the 3.4-year interval covered by the postseismic surveys. Those displacements are attributed to postseismic, right-lateral slip in the depth interval 10 to 30 km on the downward extension of the rupture trace. The right-lateral slip amounted to about 1 m directly beneath the geodetic array, and the fault-normal displacement is apparently primarily a consequence of the curvature of the rupture. These conclusions are based upon dislocation models fit to the observed deformation. However, no dislocation model was found with rms residuals as small as the expected observational error.

Uplift, thermal unrest and magma intrusion at Yellowstone caldera

The Yellowstone caldera, in the western United States, formed ∼640,000 years ago when an explosive eruption ejected ∼1,000 km3 of material. It is the youngest of a series of large calderas that formed during sequential cataclysmic eruptions that began ∼16 million years ago in eastern Oregon and northern Nevada. The Yellowstone caldera was largely buried by rhyolite lava flows during eruptions that occurred from ∼150,000 to ∼70,000 years ago. Since the last eruption, Yellowstone has remained restless, with high seismicity, continuing uplift/subsidence episodes with movements of ∼70 cm historically to several metres since the Pleistocene epoch, and intense hydrothermal activity. Here we present observations of a new mode of surface deformation in Yellowstone, based on radar interferometry observations from the European Space Agency ERS-2 satellite. We infer that the observed pattern of uplift and subsidence results from variations in the movement of molten basalt into and out of the Yellowstone volcanic system.

Postseismic strain following the 1989 Loma Prieta earthquake from GPS and leveling measurements

Postseismic deformation in the 5 years following the 1989 Loma Prieta earthquake has been measured with the Global Positioning System and precise leveling. Postearthquake velocities at distances greater than ∼20 km from the coseismic rupture are not significantly different from those observed in the 20 years prior to the earthquake. However, velocities at stations within ∼20 km of the rupture exceed preearthquake rates and exhibit unanticipated contraction normal to the strike of the San Andreas fault system. A combination of forward modeling and nonlinear optimization suggests that the observed postseismic deformations were caused by aseismic oblique reverse slip averaging 2.9 cm/yr on the San Andreas fault and/or the Loma Prieta rupture zone and 2.4 cm/yr reverse slip along a buried fault within the Foothills thrust belt. The best fitting sources of postseismic deformation are all located at depths of less than 15 km. We find no evidence for accelerated flow or shear below the Loma Prieta rupture in the first 5 years following the earthquake. The inferred postseismic slip is likely to have been caused by the coseismic stress change updip of the 1989 rupture.

Strain accumulation across the Eastern California Shear Zone at latitude 36°30′N

The motion of a linear array of monuments extending across the Eastern California Shear Zone (ECSZ) has been measured from 1994 to 1999 with the Global Positioning System. The linear array is oriented N54°E, perpendicular to the tangent to the local small circle drawn about the Pacific-North America pole of rotation, and the observed motion across the ECSZ is approximated by differential rotation about that pole. The observations suggest uniform deformation within the ECSZ (strike N23°W) (26 nstrain yr−1 extension normal to the zone and 39 nstrain yr−1 simple right-lateral shear across it) with no significant deformation in the two blocks (the Sierra Nevada mountains and southern Nevada) on either side. The deformation may be imposed by right-lateral slip at depth on the individual major fault systems within the zone if the slip rates are: Death Valley-Furnace Creek fault 3.2±0.9 mm yr−1, Hunter Mountain-Panamint Valley fault 3.3±1.6 mm yr−1, and Owens Valley fault 6.9±1.6 mm yr−1. However, this estimate of the slip rate on the Owens Valley fault is 3 times greater than the geologic estimate.

Inversion of GPS data for spatially variable slip‐rate on the San Andreas Fault near Parkfield, CA

We analyze GPS data collected from 1991-1998 at 35 sites near the Parkeld segment of the San Andreas Fault. Inverting the resultant site velocities for the distri- bution of interseismic slip-rate on the San Andreas reveals an area of low slip-rate on thefault extendingfrom between Middle Mountain and Carr Hill to southeast of Gold Hill. This slip-rate patternis similar to thatfoundbyHarris and Segall(1987) using trilateration data collected between 1966 and 1984. We infer a deep slip-rate (33 mm/yr) and depth of the transition between seismogenic and non-seismogenic slip(14km)thatagree betterwithindependentgeologic ev- idence than those found in the 1987 study. In contrast to Harris and Segall(1987), wendnoevidenceoffault-normal contraction.

Geodetic estimates of fault slip rates in the San Francisco Bay area

Bourne et al. [1998] have suggested that the interseismic velocity profile at the surface across a transform plate boundary is a replica of the secular velocity profile at depth in the plastosphere. On the other hand, in the viscoelastic coupling model the shape of the interseismic surface velocity profile is a consequence of plastosphere relaxation following the previous rupture of the faults that make up the plate boundary and is not directly related to the secular flow in the plastosphere. The two models appear to be incompatible. If the plate boundary is composed of several subparallel faults and the interseismic surface velocity profile across the boundary known, each model predicts the secular slip rates on the faults which make up the boundary. As suggested by Bourne et al., the models can then be tested by comparing the predicted secular slip rates to those estimated from long-term offsets inferred from geology. Here we apply that test to the secular slip rates predicted for the principal faults (San Andreas, San Gregorio, Hayward, Calaveras, Rodgers Creek, Green Valley and Greenville faults) in the San Andreas fault system in the San Francisco Bay area. The estimates from the two models generally agree with one another and to a lesser extent with the geologic estimate. Because the viscoelastic coupling model has been equally successful in estimating secular slip rates on the various fault strands at a diffuse plate boundary, the success of the model of Bourne et al. [1998] in doing the same thing should not be taken as proof that the interseismic velocity profile across the plate boundary at the surface is a replica of the velocity profile at depth in the plastosphere.

Deformation across the Pacific-North America plate boundary near San Francisco, California

We have detected a narrow zone of compression between the Coast Ranges and the Great Valley, and we have estimated slip rates for the San Andreas, Rodgers Creek, and Green Valley faults just north of San Francisco. These results are based on an analysis of campaign and continuous Global Positioning System (GPS) data collected between 1992 and 2000 in central California. The zone of compression between the Coast Ranges and the Great Valley is 25 km wide. The observations clearly show 3.8±1.5 mm yr−1 of shortening over this narrow zone. The strike slip components are best fit by a model with 20.8±1.9 mm yr−1 slip on the San Andreas fault, 10.3±2.6 mm yr−1 on the Rodgers Creek fault, and 8.1±2.1 mm yr−1 on the Green Valley fault. The Pacific-Sierra Nevada-Great Valley motion totals 39.2±3.8 mm yr−1 across a zone that is 120 km wide (at the latitude of San Francisco). Standard deviations are one σ. The geodetic results suggest a higher than geologic rate for the Green Valley fault. The geodetic results also suggest an inconsistency between geologic estimates of the San Andreas rate and seismologic estimates of the depth of locking on the San Andreas fault. The only convergence observed is in the narrow zone along the border between the Great Valley and the Coast Ranges.

Postseismic deformation following the 1989 (M = 7.1) Loma Prieta, California, earthquake

Postseismic deformation along a 90-km profile bisecting the projected surface trace of the coseismic rupture of the 1989 Loma Prieta earthquake has been monitored by frequent GPS surveys for 3.3 years following the earthquake. In addition to the expected deformation associated with secular strain accumulation on the San Andreas and Calaveras faults, deformation associated with postseismic readjustment has been detected. Most of that deformation can be attributed to 1.5 m right-lateral and 0.9 m reverse postseismic slip on a 5-km-wide downdip extension of the Loma Prieta rupture. In addition, there seems to be a 0.1 m postseismic collapse of the Loma Prieta rupture zone in the direction perpendicular to the plane of the rupture. The fault-normal (N48°E) surface displacements plotted as a function of time exhibit a curvature suggesting a relaxation time of about 1.4 years. Similar plots of the fault-parallel (N42°W) displacement components do not exhibit significant curvature. Presumably, the deformation shown in those plots is dominated by secular strain accumulation along the San Andreas and Calaveras faults rather than postseismic relaxation.

Deformation across the forearc of the Cascadia subduction zone at Cape Blanco, Oregon

Over the interval 1992–1999 the U.S. Geological Survey measured the deformation of a geodetic array extending N80°E (approximate direction of plate convergence) from Cape Blanco on the Oregon coast to the volcanic arc near Newberry Crater (55 and 350 km, respectively, from the deformation front). Within about 150 km from the deformation front, the forearc is being compressed arcward (N80°E) by coupling to the subducting Juan de Fuca plate. Dislocation modeling of the observed N80°E compression suggests that the main thrust zone (the locked portion of the Juan de Fuca-forearc interface) is about 40 km wide in the downdip direction. The transverse (N10°W) velocity component of the forearc measured with respect to the fixed interior of North America decreases with distance from the deformation front at a rate of about 0.03 mm yr−1 km−1. That gradient appears to be a consequence of rigid rotation of the forearc block relative to fixed interior North America (Euler vector of 43.4°±0.1°N, 120.0°±0.4°W, and −1.67±0.17 ° (m.y.)−1; quoted uncertainties are standard deviations). The rotation rate is similar to the paleomagnetically measured rotation rate (−1.0±0.2 °(m.y.)−1) of the 15 Ma lava flows along the Columbia River 250 km farther north. The back arc does not appear to participate in this rotation but rather is migrating at a rate of about 3.6 mm yr−1 northward with respect to fixed North America. That migration could be partly an artifact of an imperfect tie of our reference coordinate system to the interior of North America.

Geodetic Constraints on the 2014 M 6.0 South Napa Earthquake

On 24 August 2014, the M 6.0 South Napa earthquake shook much of the San Francisco Bay area, leading to significant damage in the Napa Valley. The earthquake occurred in the vicinity of the West Napa fault (122.313° W, 38.22° N, 11.3 km), a mapped structure located between the Rodger’s Creek and Green Valley faults, with nearly pure right‐lateral strike‐slip motion (strike 157°, dip 77°, rake –169°; http://comcat.cr.usgs.gov/earthquakes/eventpage/nc72282711#summary, last accessed December 2014) (Fig. 1). The West Napa fault previously experienced an M 5 strike‐slip event in 2000 but otherwise exhibited no previous definitive evidence of historic earthquake rupture (Rodgers et al., 2008; Wesling and Hanson, 2008). Evans et al. (2012) found slip rates of ∼9.5  mm/yr along the West Napa fault, with most slip rate models for the Bay area placing higher slip rates and greater earthquake potential on the Rodger’s Creek and Green Valley faults, respectively (e.g., Savage et al., 1999; d’Alessio et al., 2005; Funning et al., 2007).