Search USGS⌕ Search

Geology topics

J. L. Svarc

Publications and source records attributed to J. L. Svarc.

16 recordsLinked to original sources

Postseismic relaxation following the 1994 M w 6.7 Northridge earthquake, southern California

We have reexamined the postearthquake deformation of a 65 km long linear array of 11 geodetic monuments extending north–south across the rupture (reverse slip on a blind thrust dipping 40°S–20°W) associated with the 1994 M w 6.7 Northridge earthquake. That array was surveyed frequently in the interval from 4 to 2650 days after the earthquake. The velocity of each of the monuments over the interval 100–2650 days postearthquake appears to be constant. Moreover, the profile of those velocities along the length of the array is very similar to a preearthquake velocity profile for a nearby, similarly oriented array. We take this to indicate that significant postseismic relaxation is evident only in the first 100 days postseismic and that the subsequent linear trend is typical of the interseismic interval. The postseismic relaxation (postseismic displacement less displacement that would have occurred at the preseismic velocity) is found to be almost wholly parallel (N70°W) to the nearby (40 km) San Andreas Fault with only negligible relaxation in the direction of coseismic slip (N20°E) on the Northridge rupture. We suggest that the N70°W relaxation is caused by aseismic, right-lateral slip at depth on the San Andreas Fault, excess slip presumably triggered by the Northridge rupture. Finally, using the Dieterich (1994) stress-seismicity relation, we show that return to the preseismic deformation rate within 100 days following the earthquake could be consistent with the cumulative number of M > 2.5 earthquakes observed following the main shock.

California↗

Postseismic relaxation following the 1992 M7.3 Landers and 1999 M7.1 Hector Mine earthquakes, southern California

Postseismic relaxation (measured postseismic deformation less the deformation that would have occurred at the preseismic rate) has been monitored at the same 10 monuments over ???6 years following both the 1992 Landers and the 1999 Hector Mine earthquakes. For both earthquakes the displacement components of the observed relaxation are well described by ?? i + ?? i f 1 (t), where ?? i and ?? i are constants peculiar to each component at each monument, t is the time after the earthquake, and f 1 (t) is a temporal function common to all components and all monuments for that earthquake. The temporal fanction f 1 (t) can be approximated by bt + c log e (1 + t /??), where ?? = 38.7 ?? 15.2 days and 25.6 ?? 7.7 days for the Landers and Hector Mine relaxations, respectively. Because the estimated values of ?? do not differ significantly, the transient term log e (1 + t/??) in the temporal function may be the same for both earthquakes. The asymptotic (t ??? ???) relaxation rates ?? i b are only a few mm/a and do not appear to be significantly different following the two earthquakes. The asymptotic deformation rates appear to be slightly greater than the preseismic deformation rates, but the preseismic rates were not measured directly. Thus, the deformations of the Landers array measured over the first 5.6 years following the Landers earthquake and over the first 6.4 years following the Hector Mine earthquake are generally consistent with a simple model in which the transient postearthquake relaxation depends upon time as log e (1 + t/??) with ?? = 28 ?? 5 days and the asymptotic postseismic speeds of the monuments exceed the preseismic speeds by at most only a few millimeters per annum.

Journal of Geophysical Research B: Solid Earth↗

Postseismic relaxation and aftershocks

[1] Perfettini et al. (2005) suggested that the temporal dependence of surface displacements u ( t ) measured in the epicentral area following an earthquake is related to N ( t ), the cumulative number of aftershocks, by the equation u ( t ) = a + bt + cN ( t ) + d (1 − e − αt ), where a , b , c , d , and α are constants chosen to fit the data and t is the postearthquake time. N ( t ) appears in the expression for u ( t ) because both the aftershocks and a portion of u ( t ) are thought to be driven by the same source, postseismic fault creep at subseismogenic depths on the downdip extension of the coseismic rupture. We show that this equation with the actually observed N ( t ) fits the postseismic displacements recorded on several baselines following each of five earthquakes: 1999 M7.6 Chi-Chi (Taiwan), 1999 M7.1 Hector Mine (southern California), 2002 M7.9 Denali (central Alaska), 2003 M6.5 San Simeon (central California), and 2004 M6.0 Parkfield (central California) earthquakes. Although there are plausible physical interpretations for each of the terms in the expression for u ( t ), the large number of adjustable constants ( a , b , c , d , and α ) involved in fitting the rather simple postseismic displacements diminishes the significance of the fit. Because the observed N ( t ) is well fit by the modified Omori's law, fault creep at depth presumably exhibits the same temporal dependence. That dependence could be explained if the rheology of the fault downdip from the coseismic rupture is consistent with ordinary transient creep. Montesi (2004) demonstrated that power law creep across a shear zone at depth would also produce that temporal signal.

Journal of Geophysical Research B: Solid Earth↗

Postseismic relaxation and transient creep

Postseismic deformation has been observed in the epicentral area following the 1992 Landers (M = 7.3), 1999 Chi-Chi (M = 7.6), 1999 Hector Mine (M = 7.1), 2002 Denali (M = 7.9), 2003 San Simeon (M = 6.5), and 2004 Parkfield (M = 6.0) earthquakes. The observations consist of repeated GPS measurements of the position of one monument relative to another (separation ???100 km). The early observations (t < 0.1 year) are well fit by the function a' + c'log(t), where t is the time after the earthquake and a' and c' are constants chosen to fit the data. Because a log(t) time dependence is characteristic of transient (primary) creep, the early postseismic response may be governed by transient creep as Benioff proposed in 1951. That inference is provisional as the stress conditions prevailing in postseismic relaxation are not identical to the constant stress condition in creep experiments. The observed logarithmic time dependence includes no characteristic time that might aid in identifying the micromechanical cause.

Journal of Geophysical Research B: Solid Earth↗

Strain accumulation across the Coast Ranges at the latitude of San Francisco, 1994-2000

A 66-monument geodetic array spanning the Coast Ranges near San Francisco has been surveyed more than eight times by GIPS between late 1993 and early 2001. The measured horizontal velocities of the monuments are well represented by uniform, right-lateral, simple shear parallel to N29°W. (The local strike of the San Andreas Fault is ∼N34°W.) The observed areal dilatation rate of 6.9 ± 10.0 nstrain yr −1 (quoted uncertainty is one standard deviation and extension is reckoned positive) is not significantly different from zero, which implies that the observed strain accumulation could be released by strike-slip faulting alone. Our results are consistent with the slip rates assigned by the Working Group on California Earthquake Probabilities [2003] to the principal faults (San Gregorio, San Andreas, Hayward-Rodgers Creek, Calaveras-Concord-Green Valley, and Greenville Faults) cutting across the GPS array. The vector sum of those slip rates is is 39.8 ± 2.6 mm yr −1 N29.8°W ± 2.8°, whereas the motion across the GPS array (breadth 120 km) inferred from the uniform strain rate approximation is 38.7 ± 1.2 mm yr −1 N29.0°W ± 0.9° right-lateral shear and 0.4 ± 0.9 mm yr −1 N61°E ± 0.9° extension. We interpret the near coincidence of these rates and the absence of significant accumulation of areal dilatation to imply that right-lateral slip on the principal faults can release the accumulating strain; major strain release on reverse faults subparallel to the San Andreas Fault within the Coast Ranges is not required.

California↗

Interseismic strain and rotation rates in the northeast Mojave domain, eastern California

The northeast Mojave domain, a type locality for bookshelf faulting, is a region of east striking, left-lateral faults in the northeast corner of the Mojave block, a block otherwise dominated by ∼N40°W striking, right-lateral faults. Paleomagnetic evidence suggests that blocks within the domain have rotated clockwise about a vertical axis as much as 60° since 12.8 Ma [Schermer et al., 1996] . In 1994, and again in 2002, the U.S. Geological Survey surveyed an array of 14 geodetic monuments distributed across the northeast Mojave domain. The 2002 survey results were adjusted to remove the coseismic offsets imposed by the nearby Hector Mine earthquake (16 October 1999, M w = 7.1). The adjusted deformation across the array appears to be uniform and can be approximated by the principal strain rates ε 1 = 28.9 ± 9.1 N77.2°W ± 4.8° and ε 2 = −48.2 ± 8.9 N12.8°E ± 4.8° nstrain yr −1 ; extension reckoned positive, and quoted uncertainties are standard deviations. That strain accumulation could be released by slip on faults striking N32°W but not by bookshelf faulting on the east striking faults alone. The vertical axis rotation rate of the northeast Mojave domain as a whole relative to fixed North America is 71.0 ± 6.4 nrad yr −1 (4.07° ± 0.37° Myr −1 ) clockwise, about twice the maximum tensor shear strain rate. The observed rotation rate acting over 12.8 Myr would produce a clockwise rotation of 52.1° ± 4.7°, exclusive of possible coseismic rotations. That rotation is in rough agreement with the paleomagnetic rotation accumulated in the individual fault blocks within the northeast Mojave domain since 12.8 Ma.

California↗

Near-field postseismic deformation associated with the 1992 Landers and 1999 Hector Mine, California, earthquakes

After the Landers earthquake (Mw = 7.3, 1992.489) a linear array of 10 monuments extending about 30 km N50??E on either side of the earthquake rupture plus a nearby offtrend reference monument were surveyed frequently by GPS until 2003.2. The array also spans the rupture of the subsequent Hector Mine earthquake (Mw = 7.1, 1999.792 . The pre-Landers velocities of monuments in the array relative to interior North America were estimated from earlier trilateration and very long baseline interferometry measurements. Except at the reference monument, the post-Landers velocities of the individual monuments in the array relaxed to their preseismic values within 4 years. Following the Hector Mine earthquake the velocities of the monuments relaxed to steady rates within 1 year. Those steady rates for the east components are about equal to the pre-Landers rates as is the steady rate for the north component of the one monument east of the Hector Mine rupture. However, the steady rates for the north components of the 10 monuments west of the rupture are systematically ???10 mm yr1 larger than the pre-Landers rates. The relaxation to a steady rate is approximately exponential with decay times of 0.50 ?? 0.10 year following the Landers earthquake and 0.32 ?? 0.18 year following the Hector Mine earthquake. The postearthquake motions of the Landers array following the Landers earthquake are not well approximated by the viscoelastic-coupling model of Pollitz et al. [2000]. A similar viscoelastic-coupling model [Pollitz et al., 2001] is more successful in representing the deformation after the Hector Mine earthquake.

Journal of Geophysical Research B: Solid Earth↗

Initiation of deformation of the Eastern California Shear Zone: Constraints from Garlock fault geometry and GPS observations

We suggest a 2-stage deformation model for the Eastern California Shear Zone (ECSZ) to explain the geometry of the Garlock fault trace. We assume the Garlock fault was originally straight and then was gradually curved by right-lateral shear deformation across the ECSZ. In our 2-stage deformation model, the first stage involves uniform shear deformation across the eastern part of the shear zone, and the second stage involves uniform shear deformation across the entire shear zone. In addition to the current shape of the Garlock fault, our model incorporates constraints on contemporary deformation rates provided by GPS observations. We find that the best fitting age for initiation of shear in eastern part of the ECSZ is about 5.0 ?? 0.4 Ma, and that deformation of the western part started about 1.6 Myr later.

Geophysical Research Letters↗

Strain accumulation and rotation in western Nevada, 1993-2000

The positions of 44 GPS monuments in an array extending from the Sierra Nevada at the latitude of Reno to near Austin, Nevada, have been measured several times in the 1993–2000 interval. The western half of the array spans the Walker Lane belt, whereas the eastern half spans the central Nevada seismic zone (CNSZ). The principal strain rates in the Walker Lane belt are 29.6 ± 5.3 nstrain yr −1 N88.4°E ± 5.4° and −12.8 ± 6.0 nanostrain yr −1 N01.6°W ± 5.4°, extension reckoned positive, and the clockwise (as seen from above the Earth) rotation rate about a vertical axis is 13.6 ± 4.0 nrad yr −1 . The quoted uncertainties are standard deviations. The motion in the Walker Lane belt can then be represented by a zone striking N35°W subject to 16.8 ± 4.9 nstrain yr −1 extension perpendicular to it and 19.5 ± 4.0 nstrain yr −1 right-lateral, simple shear across it. The N35°W strike of the zone is the same as the direction of the local tangent to the small circle drawn about the Pacific-North America pole of rotation. The principal strain rates for the CNSZ are 46.2 ± 11.0 nstrain yr −1 N49.9°W ± 6.0° and −13.6 ± 6.1 nstrain yr −1 N40.1°E ± 6.0°, and the clockwise rotation rate about a vertical axis is 20.3 ± 6.3 nrad yr −1 . The motion across the CNSZ can then be represented by a zone striking N12°E subject to 32.6 ± 11.0 nstrain yr −1 extension perpendicular to it and 25.1 ± 6.3 nstrain yr −1 right-lateral, simple shear across it. The N12°E strike of the zone is similar to the strikes of the faults (Rainbow Mountain, Fairview Peak, and Dixie Valley) within it.

California, Nevada↗

Strain accumulation and rotation in western Oregon and southwestern Washington

Velocities of 75 geodetic monuments in western Oregon and southwestern Washington extending from the coast to more than 300 km inland have been determined from GPS surveys over the interval 1992–2000. The average standard deviation in each of the horizontal velocity components is ∼1 mm yr −1 . The observed velocity field is approximated by a combination of rigid rotation (Euler vector relative to interior North America: 43.40°N ± 0.14°, 119.33°W ± 0.28°, and 0.822 ± 0.057° Myr −1 clockwise; quoted uncertainties are standard deviations), uniform regional strain rate (ε EE = −7.4 ± 1.8, ε EN = −3.4 ± 1.0, and ε NN = −5.0 ± 0.8 nstrain yr −1 , extension reckoned positive), and a dislocation model representing subduction of the Juan de Fuca plate beneath North America. Subduction south of 44.5°N was represented by a 40‐km‐wide locked thrust and subduction north of 44.5°N by a 75‐km‐wide locked thrust.

Oregon, Washington↗

Strain accumulation and rotation in the Eastern California Shear Zone

Although the Eastern California Shear Zone (ECSZ) (strike ∼N25°W) does not quite coincide with a small circle drawn about the Pacific-North America pole of rotation, trilateration and GPS measurements demonstrate that the motion within the zone corresponds to right-lateral simple shear across a vertical plane (strike N33°W±5°) roughly parallel to the tangent to that local small circle (strike ∼N40°W). If the simple shear is released by slip on faults subparallel to the shear zone, the accumulated rotation is also released, leaving no secular rotation. South of the Garlock fault the principal faults (e.g., Calico-Blackwater fault) strike ∼N40°W, close enough to the strike of the vertical plane across which maximum right-lateral shear accumulates to almost wholly accommodate that accumulation of both strain and rotation by right-lateral slip. North of the Garlock fault dip slip as well as strike slip on the principal faults (strike ∼N20°W) is required to accommodate the simple shear accumulation. In both cases the accumulated rotation is released with the shear strain. The Garlock fault, which transects the ECSZ, is not offset by north-northwest striking faults nor, despite geological evidence for long-term left-lateral slip, does it appear at the present time to be accumulating left-lateral simple shear strain across the fault due to slip at depth. Rather the motion is explained by right-lateral simple shear across the orthogonal ECSZ. Left-lateral slip on the Garlock fault will release the shear strain accumulating there but would augment the accumulating rotation, resulting in a secular clockwise rotation rate ∼80 nrad yr −1 (4.6° Myr −1 ).

California↗

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.

California↗

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 N880°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.

Journal of Geophysical Research B: Solid Earth↗

Deformation across the Alaska-Aleutian Subduction Zone near Kodiak

The Kodiak-Katmai geodetic array, nine monuments distributed along a profile trending north-northwestward across Kodiak Island and the Alaska Peninsula, was surveyed in 1993, 1995 and 1997 to determine the deformation at the Alaska-Aleutian subduction zone. Velocities on Kodiak island measured relative to the stable North American plate decrease with distance from the Alaska-Aleutian trench (distance range 106 to 250 km), whereas no appreciable deformation was measured on the Alaska Peninsula (distances 250 to 370 km from the trench). The measured deformation is reasonably well predicted by the conventional dislocation representation of subduction with the model parameters determined independently (i.e., not simply by fitting the observations). The deformation of Kodiak Island is in striking contrast to the very minor deformation measured in the similarly situated Shumagin Islands, 450 km southwest of Kodiak along the Alaska-Aleutian trench.

Geophysical Research Letters↗

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 at. [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.

Journal of Geophysical Research B: Solid Earth↗