Geology Reports⌕ Search

Geology topics

J.C. Savage

Publications and source records attributed to J.C. Savage.

At least 37 records · Page 2Linked to original sources
USGS2001

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↗
USGS2001

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↗
USGS2000

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.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

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↗
USGS1998

Strain accumulation rates in the San Francisco Bay area, 1972–1989

Maps of the strain accumulation rate in the San Francisco Bay area have been constructed from trilateration observations extending from about 1972 until the Loma Prieta earthquake in late 1989. The observations were corrected to remove offsets imposed by shallow fault creep and by four M ∼6 earthquakes that occurred in the Bay area during that time interval. The Bay area was divided into 32 contiguous polygons, and the uniform (in both space and time) strain rates that best explain the changes in the corrected (earthquake and shallow fault creep offsets removed) distances within each polygon were calculated. In a coordinate system with the 1 axis directed N58°E and the 2 axis N32°W (perpendicular and parallel to the local tangent to the small circle drawn about the Pacific-Sierra Nevada pole of rotation) the averages of these 32 strain rates (each weighted by the area of the polygon) are ε11=9.2±7.4, ε12=−160.7±4.6, and ε22=8.2±6.2 nanostrain/yr, where extension is reckoned positive and quoted uncertainties are standard deviations. As expected from the Pacific-Sierra Nevada relative plate motion, the overall strain rate is predominantly right-lateral shear across a vertical plane striking N32°W. The net increase in the 12,225 km2 area of the trilateration network is only 212±110 m2/yr, which arises from almost equal extensions in the N32°W and N58°E directions. Within the network the strain rates vary from polygon to polygon. Significant areal dilatation rates are observed in almost ⅓ of the individual polygons and the N32°W extension rates tend to be negative to the west of the Hayward-Rodgers Creek fault trend and positive east of it. The N58E extension rate is positive in 22 out of 32 polygons, a proportion that is significantly larger than would be expected by chance if the N58°E extension rate were zero or negative. The pre-1989 strain accumulation across the eventual site of the Loma Prieta rupture involves fault normal contraction as well as right-lateral shear, consistent with the rupture mechanism.

California↗
USGS1998

Viscoelastic coupling model of the San Andreas fault along the Big Bend, southern California

The big bend segment of the San Andreas fault is the 300-km-long segment in southern California that strikes about N65°W, roughly 25° counterclockwise from the local tangent to the small circle about the Pacific-North America pole of rotation. The broad distribution of deformation of trilateration networks along this segment implies a locking depth of at least 25 km as interpreted by the conventional model of strain accumulation (continuous slip on the fault below the locking depth at the rate of relative plate motion), whereas the observed seismicity and laboratory data on fault strength suggest that the locking depth should be no greater than 10 to 15 km. The discrepancy is explained by the viscoelastic coupling model which accounts for the viscoelastic response of the lower crust. Thus the broad distribution of deformation observed across the big bend segment can be largely associated with the San Andreas fault itself, not subsidiary faults distributed throughout the region. The Working Group on California Earthquake Probabilities [1995] in using geodetic data to estimate the seismic risk in southern California has assumed that strain accumulated off the San Andreas fault is released by earthquakes located off the San Andreas fault. Thus they count the San Andreas contribution to total seismic moment accumulation more than once, leading to an overestimate of the seismicity for magnitude 6 and greater earthquakes in their Type C zones.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Displacement field for an edge dislocation in a layered half-space

The displacement field for an edge dislocation in an Earth model consisting of a layer welded to a half-space of different material is found in the form of a Fourier integral following the method given by Weeks et al. [1968]. There are four elementary solutions to be considered: the dislocation is either in the half-space or the layer and the Burgers vector is either parallel or perpendicular to the layer. A general two-dimensional solution for a dip-slip faulting or dike injection (arbitrary dip) can be constructed from a superposition of these elementary solutions. Surface deformations have been calculated for an edge dislocation located at the interface with Burgers vector inclined 0°, 30°, 60°, and 90° to the interface for the case where the rigidity of the layer is half of that of the half-space and the Poisson ratios are the same. Those displacement fields have been compared to the displacement fields generated by similarly situated edge dislocations in a uniform half-space. The surface displacement field produced by the edge dislocation in the layered half-space is very similar to that produced by an edge dislocation at a different depth in a uniform half-space. In general, a low-modulus (high-modulus) layer causes the half-space equivalent dislocation to appear shallower (deeper) than the actual dislocation in the layered half-space.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

A test of the double-shearing model of flow for granular materials

The double-shearing model of flow attributes plastic deformation in granular materials to cooperative slip on conjugate Coulomb shears (surfaces upon which the Coulomb yield condition is satisfied). The strict formulation of the double-shearing model then requires that the slip lines in the material coincide with the Coulomb shears. Three different experiments that approximate simple shear deformation in granular media appear to be inconsistent with this strict formulation. For example, the orientation of the principal stress axes in a layer of sand driven in steady, simple shear was measured subject to the assumption that the Coulomb failure criterion was satisfied on some surfaces (orientation unspecified) within the sand layer. The orientation of the inferred principal compressive axis was then compared with the orientations predicted by the double-shearing model. The strict formulation of the model [ Spencer , 1982] predicts that the principal stress axes should rotate in a sense opposite to that inferred from the experiments. A less restrictive formulation of the double-shearing model by de Josselin de Jong [1971] does not completely specify the solution but does prescribe limits on the possible orientations of the principal stress axes. The orientations of the principal compression axis inferred from the experiments are probably within those limits. An elastoplastic formulation of the double-shearing model [ de Josselin de Jong , 1988] is reasonably consistent with the experiments, although quantitative agreement was not attained. Thus we conclude that the double-shearing model may be a viable law to describe deformation of granular materials, but the macroscopic slip surfaces will not in general coincide with the Coulomb shears.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Observed discrepancy between geodolite and GPS distance measurements

Comparison of contemporaneous measurements of 84 distances in the range of 10 to 50 km by both Global Positioning System (GPS) and Geodolite (an electro-optical distance-measuring instrument) indicates that the Geodolite measurements are systematically longer by 0.283 ± 0.100 parts per million of the measured distance. Quoted uncertainty is 1 standard deviation. This amounts to 11.3 ± 4.0 mm at 40 km, which is near the maximum Geodolite range. The systematic difference is within the random uncertainty of an individual GPS-Geodolite comparison and was detected only from an analysis of a large number (84) of such comparisons. The source of the systematic difference has not been identified.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Failure in laboratory fault models in triaxial tests

A model of a fault in the Earth is a sand-filled saw cut in a granite cylinder subjected to a triaxial test. The saw cut is inclined at an angle α to the cylinder axis, and the sand filling is intended to represent gouge. The triaxial test subjects the granite cylinder to a constant confining pressure and increasing axial stress to maintain a constant rate of shortening of the cylinder. The required axial stress increases at a decreasing rate to a maximum, beyond which a roughly constant axial stress is sufficient to maintain the constant rate of shortening. Such triaxial tests were run for saw cuts inclined at angles α of 20°, 25°, 30°, 35°, 40°, 45°, and 50° to the cylinder axis, and the apparent coefficient of friction μ a (ratio of the shear stress to the normal stress, both stresses resolved onto the saw cut) at failure was determined. Subject to the assumption that the observed failure involves slip on Coulomb shears (orientation unspecified), the orientation of the principal compression axis within the gouge can be calculated as a function of μ a for a given value of the coefficient of internal friction μ i . The rotation of the principal stress axes within the gouge in a triaxial test can then be followed as the shear strain across the gouge layer increases. For μ i ∼ 0.8, an appropriate value for highly sheared sand, the observed values μ a imply that the principal axis of compression within the gouge rotates so as to approach being parallel to the cylinder axis for all saw cut angles (20° < α < 50°). In the limiting state (principal compression axis parallel to cylinder axis) the stress state in the gouge layer would be the same as that in the granite cylinder, and the failure criterion would be independent of the saw cut angle.

Journal of Geophysical Research B: Solid Earth↗
USGS1991

Strain measurements and the potential for a great subduction earthquake off the coast of Washington

Geodetic measurements of deformation in northwestern Washington indicate that strain is accumulating at a rate close to that predicted by a model of the Cascadia subduction zone in which the plate interface underlying the continental slope and outer continental shelf is currently locked but the remainder of the interface slips continuously. Presumably this locked segment will eventually rupture in a great thrust earthquake with a down-dip extent greater than 100 kilometers.

Washington↗
USGS1986

Strain accumulation in the Shumagin and Yakataga seismic gaps, Alaska

Strain accumulation during the 1980-85 interval has been measured by means of trilateration surveys in the Shumagin and Yakataga seismic gaps, which are the two regions identified as the most likely sites for the next great thrust earthquakes along the Alaska-Aleutian arc. No significant strain accumulation was detected in the Shumagin gap, but experience at similar subduction zones and simple models of the subduction process suggest that a measurable amount of strain should have accumulated. The most likely explanation of the observation is that subduction mere is either aseismic or episodic. The strain accumulation measured in the Yakataga gap is consistent with that expected for the plate convergence rate, although the direction of maximum compression may suggest a somewhat more oblique convergence than expected.

Alaska↗
USGS1982

Magmatic resurgence in Long Valley caldera, California: Possible cause of the 1980 Mammoth Lakes earthquakes

Changes in elevation between 1975 and October 1980 along a leveling line across the Long Valley caldera indicate a broad (half-width, 15 kilometers) uplift (maximum, 0.25 meter) centered on the old resurgent dome. This uplift is consistent with reinflation of a magma reservoir at a depth of about 10 kilometers. Stresses generated by this magmatic resurgence may have caused the sequence of four magnitude 6 earthquakes near Mammoth Lakes in May 1980.

California↗
USGS1981

Strain on the San Andreas fault near Palmdale, California: Rapid, aseismic change

Frequently repeated strain measurements near Palmdale, California, during the period from 1971 through 1980 indicate that, in addition to a uniform accumulation of right-lateral shear strain (engineering shear, 0.35 microradian per year) across the San Andreas fault, a 1-microstrain contraction perpendicular to the fault that accumulated gradually during the interval 1974 through 1978 was aseismically released between February and November 1979. Subsequently (November 1979 to March 1980), about half of the contraction was recovered. This sequence of strain changes can be explained in terms of south-southwestward migration of a slip event consisting of the south-southwestward movement of the upper crust on a horizontal detachment surface at a depth of 10 to 30 kilometers. The large strain change in 1979 corresponds to the passage of the slip event beneath the San Andreas fault.

California↗
USGS1978

Strain in southern California: Measured uniaxial north-south regional contraction

The plate tectonics model of the Pacific moving northwest relative to North America implies that the regional strain in California should be simple shear across a vertical plane striking N45°W or equivalently equal parts of north-south contraction and east-west extension. Measurements of the strain accumulation at seven separate sites in southern California in the interval 1972 through 1978 indicate a remarkably consistent uniaxial north-south contraction of about 0.3 part per million per year; the expected east-west extension is absent. It is not clear whether the period from 1972 through 1978 is anomalous or whether the secular strain in southern California is indeed a uniaxial north-south contraction.

California↗
USGS1978

Geodimeter measurements and the Southern California uplift

Modern surveying instruments, such as geodimeters, are capable of measuring distances in the range of 1 to 30 kilometers with remarkable precision. Indeed, the present limitation upon the precision of measurement is not the resolution of the instruments themselves but rather the uncertainty introduced by variations in the velocity of light in the atmosphere between the two endpoints of the measured distance. This capability in precise distance measurement can be applied to earthquake studies by using repeated distance measurements to determine changes in the distance between monuments located along the major faults. Such measurements have been made along the San Andreas fault since late 1959, and, as a result, a wealth of data on crustal deformation in that area is now available. In fact, geodimeter measurements of this type furnish a more stable measure of secular strain (change in the ratio of length to length) than any continuously recording strain meter. The superiority of the geodimeter measurement stems principally from the long base line measured which averages over local inhomogeneities (cracks and joints in the rock, inclusions, and so forth). This article describes the important features of the measurement technique as well as some results for the region of the southern California uplift.

California↗
USGS1974

Geodimeter measurements of slip and strain accumulation along the San Andreas fault

The U.S. Geological Survey conducts repeated geodimeter surveys of trilateration networks in central California in order to study the processes of slip and strain accumulation along the San Andreas fault. The precision of distance measurement is described by a standard deviation σ = (a 2 + b 2 L 2 ) 12 where a = (a 2 + b 2 L 2 ) 12 where a = 3mm, b = 2 · 10 −7 , and L is the line length. Within the precision of measurement, no anomalous strain episodes preceding earthquakes or even strain discontinuities at the time of earthquakes were detected from repeated measurements of lines near the epicenters of small (magnitude 4.5–5.1) earthquakes. Annual measurements of small (5-km aperture) strain polygons near the San Andreas fault have not proved strain accumulation in a 3-year period. Repeated measurements of longer lines over periods of 8 to 14 years indicate changes that cannot be attributed to fault slip and must represent strain accumulation at the level of a few parts in 10 7 per year.

California↗