Search USGS⌕ Search

Geology topics

James C. Savage

Publications and source records attributed to James C. Savage.

At least 55 records · Page 3Linked to original sources
USGS2012

Deformation across the Salton Trough, California, 1973-1977

A trilateration network extending across the San Andreas, San Jacinto, and Elsinore faults in the vicinity of the Salton Sea, California, has been surveyed to very high precision several times in the 5‐year interval 1973–1977. The average strain across the entire network is essentially a uniaxial north‐south contraction at the rate of about 0.3 μstrain/a. There is no substantial strain perpendicular to the Salton trough, indicating no tendency to either open or close that rift. The observed uniaxial north‐south contraction differs from a pure shear parallel to the major faults by a uniform dilatation of about −0.25 μstrain/a that is only partly explained. The shear strain across the network reaches a maximum near the San Jacinto fault and could be explained by right lateral slip at depth on that fault. The displacement pattern also suggests right lateral slip at depth on the San Jacinto fault with total right lateral relative movement near 50 mm/a across the 120‐km breadth of the network. The deformation appears to be uniform in time in the 1973–1977 interval. A dislocation model of the deformation suggests about 50±15 mm/a of relative right lateral slip at depth distributed between the San Andreas, San Jacinto, and Elsinore faults and demonstrates that a significant average dilatation can be generated by pure strike slip on several subparallel faults.

California↗
USGS2012

Strain accumulation rates in the western United States between 1970 and 1978

The rate of dilatation and the rate and direction of shear have been determined from trilateration data for 23 Geodolite networks in the western United States. Sixteen nets are located along the San Andreas fault system between Point Reyes, California, and the United States‐Mexico border. Other locations are across the Garlock fault in California; across Puget Sound near Seattle, Washington; near Hanford in eastern Washington; near Hebgen Lake in Montana; across the Wasatch fault at Ogden, Utah; across the Rio Grande rift at Socorro, New Mexico; and Dixie Valley in Nevada; and at the northern end of Owens Valley on the California‐Nevada border. Implicit in the treatment are the assumptions that the strain was accumulating at a constant rate over the time period (within the interval 1970–1978) and over the local area (usually about 50‐km diameter) covered by the surveys. Of the nets located away from the San Andreas fault, only Ogden and Hebgen show significant strain accumulation. At Ogden the deformation is principally an east‐west compression of 0.23±0.05 μstrain/yr and at Hebgen Lake a northeast‐southwest extension of 0.17±0.03 μstrain/yr. Along the San Andreas fault system the rate of shear is 0.2 to 0.4μ/yr. The direction of shear agrees very well with the surface strike of nearby faults. This agreement is maintained even in regions like the ‘big bend,’ where both the fault strike and the observed shear direction are more westerly than they are elsewhere. Shear strain in northern California appears to be concentrated more closely on the faults, whereas in southern California the strain is a broader, smoother feature. In the San Francisco Bay area the strain data indicate slip at depth on both the San Andreas and the Calaveras faults. In addition to the observed shear the nets in California indicate a negative dilatation (areal decrease) of about 0.2 μstrain/yr. This dilatation is unexplained, but the following sources appear unlikely: (1) systematic survey error, (2) an association with the southern California uplift, (3) an association with the big bend in the San Andreas fault in Southern California, or (4) the result of the superposition of a uniaxial strain on the Pacific‐North American plate boundary shear.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Geodetic measurement of crustal deformation on the San Andreas, Hayward, and Calaveras faults near San Francisco, California

Analysis of a geodetic network of 115 lines crossing the San Andreas, Hayward, and Calaveras faults in the vicinity of San Francisco Bay and measured repeatedly between 1970 and 1980 has revealed details about the accommodation of relative plate motion in this area. The most striking result is that the deformation is not uniformly distributed across the area. In the east bay, along the Hayward and Calaveras faults, all motion appears to take place as slip directly on the fault, with no accumulation of strain in the adjacent crust. On both the Calaveras and the Hayward faults the rate obtained for the 1970 to 1980 period agrees with geologic rates spanning a few million years and with creep rates spanning a few decades. The Hayward fault slip rate is 7 ± 1 mm/yr. The Calaveras fault slip rate is 7 ± 1 mm/yr, with perhaps half of this slip distributed across a zone a few kilometers wide, probably as inelastic deformation of weak near‐surface material. The absence of strain accumulation in the east bay is surprising since the Hayward and Calaveras faults have been the site of large earthquakes in the past. A block located east of the Calaveras fault and south of the Las Positas fault has been rotating clockwise at a rate of 0.3 ± 0.1 μrad/yr with very little internal deformation. Along the San Francisco peninsula no detectable slip occurs (less than 1.5 mm/yr) at the surface, but appreciable strain is accumulating. Near fault shear strain rates are 0.6 ± 0.1 μstrain/yr (engineering) with direction N47°W ± 9. The slip rate near the San Andreas fault is 12.2 ± 3.9 mm/yr distributed across a broad zone. The relative motion across the whole region during the period 1970–1980 is 32.1 ± 7.4 mm/yr.

California↗
USGS2012

A dislocation model of strain accumulation and release at a subduction zone

Strain accumulation and release at a subduction zone are attributed to stick slip on the main thrust zone and steady aseismic slip on the remainder of the plate interface. This process can be described as a superposition of steady state subduction and a repetitive cycle of slip on the main thrust zone, consisting of steady normal slip at the plate convergence rate plus occasional thrust events that recover the accumulated normal slip. Because steady state subduction does not contribute to the deformation at the free surface, deformation observed there is completely equivalent to that produced by the slip cycle alone. The response to that slip is simply the response of a particular earth model to embedded dislocations. For a purely elastic earth model, the deformation cycle consists of a coseismic offset followed by a linear‐in‐time recovery to the initial value during the interval between earthquakes. For an elastic‐viscoelastic earth model (elastic lithosphere over a viscoelastic asthenosphere), the postearthquake recovery is not linear in time. Records of local uplift as a function of time indicate that the long‐term postseismic recovery is approximately linear, suggesting that elastic earth models are adequate to describe the deformation cycle. However, the deformation predicted for a simple elastic half‐space earth model does not reproduce the deformation observed along the subduction zones in Japan at all well if stick slip is restricted to the main thrust zone. As recognized earlier by Shimazaki, Seno, and Kato, the uplift profiles could be explained if stick slip were postulated to extend along the plate interface beyond the main thrust zone to a depth of perhaps 100 km, but independent evidence suggests that stick slip at such depths is unlikely.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Precision of geodolite surveys: A reply to Jackson and Cheng

Jackson and Cheng have suggested that the changes in areal dilatation measured in Geodolite surveys by the U.S. Geological Survey may be simply an artifact of the measuring system. Although systematic error could conceivably account for the observed excursions in dilatation, we maintain that the specific criticisms by Jackson and Cheng are incorrect: the excursions in dilatation cannot be attributed to the offset correction nor to proportional error associated with temperature. The absence of both errors is demonstrated by using data that are particularly sensitive to the two effects.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Regional deformation near Palmdale, California, 1973-1983 (USA)

The Tehachapi trilateration network spans the intersection of the San Andreas and Garlock faults in southern California in the “Big Bend” region of the San Andreas fault. Analysis of data from 1973–1983 shows strain differences between the northwest and southeast regions of the network and slip at depth on both faults. The Palmdale network, spanning the San Andreas fault entirely within the Tehachapi network, showed increases of about 1 μstrain in both east‐west and north‐south extension in late 1979. The Tehachapi strains also jumped at this time, but the magnitude of the increase was only about one third that of Palmdale. The principal strain rates for Tehachapi over the time interval 1973–1983 are μstrain/yr and μstrain/yr, with the 1 axis directed N76°E. Strains were also computed for two subregions. The principal strain rates for the southeast Tehachapi region are μstrain/yr and μStrain/yr, with the 1 axis directed N73°E. This result differs significantly from the principal strain rates at Palmdale ( μstrain/yr, μstrain/yr, with the 1 axis directed N71°E). For the northwest Tehachapi region, μstrain/yr and μstrain/yr, with the 1 axis directed N93°E.

California↗
USGS2012

Deformation in the White Mountain seismic gap, California-Nevada, 1972-1982

A 100×40 km trilateration network extending from Bishop, California, to near Hawthorne, Nevada, crosses the east end of the Long Valley caldera, site of renewed magma inflation in the 1979–1980 interval, and spans most of the White Mountain seismic gap. The network was surveyed in 1972, 1973, 1976, 1979, 1980, and 1982. The 1980 survey may be contaminated by a scale error. In addition, leveling surveys across the caldera have been run in 1932, 1957, 1975, 1980, 1982, and 1983. Interpretation of the deformation is complicated by the occurrence of the May 1980 Mammoth Lakes earthquake sequence (four earthquakes M L ≥6) at the south edge of the caldera as well as other moderate earthquakes within the White Mountain seismic gap. The vertical deformation is largely accounted for by 0.10‐ to 0.15‐km 3 expansion of a spherical magma chamber 8–10 km beneath the resurgent dome within the Long Valley caldera sometime between July 1979 and September 1980 with an additional expansion of perhaps 0.05 km 3 between September 1980 and July 1982. Some additional sources of deformation within the aftershock zone of the Mammoth Lakes earthquakes seem to be required to explain the horizontal deformation. We show that right‐lateral slip on vertical faults extending WNW from each of the three largest earthquakes in the Mammoth Lakes sequence provides the required additional deformation, but this solution is by no means unique. There are simply too few data to define the rather complex deformation that apparently occurred within the aftershock zone. There is little doubt, however, that inflation of a magma chamber beneath the resurgent dome within the Long Valley caldera was involved in the deformation.

California↗
USGS2012

Earthquake swarm in Long Valley caldera, California, January 1983: Evidence for dike inflation

The 1982–1983 deformation observed by trilateration and leveling surveys across the Long Valley caldera is apparently related to the 8.5‐km‐long by 8‐km‐deep vertical rupture surface defined by the January 1983 earthquake swarm that occurred in the south moat of the caldera. The observed deformation can be explained as follows. In late 1982, 0.03 km 3 of magma was injected into a dike that dips 30° northward from the bottom of the rupture surface. The downdip dimension of this dike is 8 km. The dike inflation accounts for the uplift observed across the caldera as well as some of the horizontal deformation. Inflation of the dike generated a tension of about 3 bars across the vertical plane that was to become the rupture surface of the January swarm. This reduced the frictional stress on the rupture plane and perhaps triggered the slip that caused the January swarm. Right‐lateral slip (0.22 m) on the uppermost 2 km of the rupture plane during and after the January swarm accounts for the additional horizontal deformation observed. The model can be improved marginally if strike slip is admitted over the entire rupture surface and 0.006 km 3 of magma is injected along that surface in the depth interval 3–8 km. The improvement in the model fit, however, is not sufficient to require shallow injection of magma. Thus we conclude that inflation of a dike at depth (8–12 km) dipping northward beneath the resurgent dome plus shallow right‐lateral slip on the rupture surface is a simple, but not unique, explanation of the observed deformation and seismicity.

California↗
USGS2012

A plate flexure approximation to postseismic and interseismic deformation

The rather large postseismic deformation that is associated with two‐dimensional dip‐slip faulting in the lithosphere is related to the bending of a free plate generated by dip‐slip faulting. In the absence of gravity, asthenosphere relaxation eventually permits the faulted lithosphere to assume the dihedral configuration of a faulted free plate. For thrust faulting, the faulted area is depressed into the asthenosphere, and the flanks of the plate slope uniformly upward. In the presence of gravity, buoyancy forces act upon the plate, and the ultimate ( t →∞) postseismic uplift is approximated ( y /α > 0.2) by w 1 exp (− y /α) (cos y /α ‐ sin y /α), where w 1 is an explicit function of the fault parameters and lithosphere properties, y is the horizontal distance from the downdip end of the fault, and α is the flexural parameter for the lithosphere plate. The relaxed‐asthenosphere response is the sum of the coseismic deformation (unrelaxed‐astheriosphere response) and this ultimate postseismic deformation. The annual, steady deformation associated with strain accumulation at a subduction zone is simply the relaxed‐asthenosphere response to virtual, normal slip on the main thrust zone of an amount equal to the annual plate convergence. From these relations we have estimated the deformation (strain and uplift) expected along the Pacific coast of northern Honshu. The estimated deformation is about twice the observed deformation, but the predicted geographical distribution of the deformation is similar to the observed deformation.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Effect of crustal layering upon dislocation modeling

Slip distribution at depth on a fault may be inferred from the deformation observed on the surface. In inverting the surface deformation data to obtain the slip distribution, the Earth is generally approximated by an elastic half‐space. Slip distributions inferred from a half‐space model may contain artifacts, including zones of reversed slip, due solely to effects of layering in the real Earth. This effect is demonstrated for a vertical strike‐slip fault in an Earth consisting of an elastic layer overlying an elastic half‐space. Slip on the fault is taken to be independent of the along‐strike coordinate (i.e., antiplane strain is assumed). For a given slip distribution in this model the slip distribution on a similar fault in an elastic half‐space is found that produces the identical surface deformation. Comparison of the two slip distributions reveals structure introduced into the half‐space equivalent slip profile by crustal layering. The comparisons suggest that low‐resolution inversion schemes (e.g., single screw dislocation models) are not drastically affected by Earth structure, but attempts at detailed inversion are likely to produce profiles contaminated by artifacts of Earth structure.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

No: The L.A. array is not ready for prime time

Although much interest will focus upon the temporal behavior of observed deformation, the principal justification for the SCIGN array is that within a 5‐year interval it will provide an accurate and detailed determination of the velocity field in the Los Angeles basin that can be used to identify the active faults and estimate their secular slip rates. Obviously, the accuracy of the measurements will determine the success of the SCIGN array in reaching its objective. Over the past several years, Duncan Agnew, Hadley Johnson, and Frank Wyatt have developed arguments that quantify the accuracy likely to be obtained in the measurements and the accuracy that will be required to resolve slip rates on individual faults. In view of those arguments I do not believe that the SCIGN strategy can accomplish its objective. Instead, a mix of annual and continuous GPS surveys may be a more cost‐effective way to accomplish what actually can be achieved.

California↗
USGS2011

Earthquake mechanism and displacement fields close to fault zones: Report on the Sixth GEOP Research Conference

The Sixth Geodesy/Solid Earth and Ocean Physics (GEOP) Research Conwas held on February 4–5, 1974, at the Institute of Geophysics and Planetary Physics, University of California, San Diego, in La Jolla, California. It was attended by about 100 persons. James N. Brune, program chairman, opened the conference and delivered the introductory address, a somewhat extended version of which is printed elsewhere in this issue. Brune's paper and the following summaries of the sessions constitute a report of the conference.

Eos, Transactions, American Geophysical Union↗
USGS2010

Postseismic relaxation following the 1989 MS7.1 Loma Prieta earthquake, central California

The postseismic relaxation (postseismic displacement less displacement that would have occurred at the preseismic rate) measured by GPS and leveling following the 1989 M S 7.1 Loma Prieta earthquake is reexamined. The temporal dependence of the relaxation over the first 1200 days postseismic is well described by 1 − e − t / τ , where τ = 414 ± 92 days. (Quoted uncertainties are standard deviations.) That temporal dependence appears to be a linear function of the cumulative number of M > 2.5 aftershocks that have occurred. The relaxation is attributed to afterslip (1.56 ± 0.20 m dextral strike slip and 0.60 ± 0.04 m reverse slip) on the downdip extension (depth, 16–21 km) of the coseismic rupture plus a collapse (0.11 ± 0.02 m fault-normal displacement) of the rupture zone (depth, 5–16 km). Because the postseismic uplift was determined by leveling over a route with relief in excess of 1000 m, an allowance (18 ± 2 ppm of height above the base elevation) for excess (beyond corrections already applied) height-dependent error in the measured uplift was estimated simultaneously with the afterslip and collapse parameters. This new solution for afterslip and collapse on the plane of the rupture provides an alternative explanation to the suggestion by Bürgmann et al. (1997) that the postseismic deformation was due to afterslip on the coseismic rupture and on a nearby, shallow, 30° dipping thrust.

California↗
USGS2010

Calculation of aftershock accumulation from observed postseismic deformation: M6 2004 Parkfield, California, earthquake

[1] The postseismic stress accumulation τ ( t ) over the interval 0.004 to 880 days following the 2004 Parkfield earthquake (M6) can be inferred from GPS measurements of postseismic deformation. The stress relaxation τ ( t ) − τ ′ l t , where τ ′ l is the interseismic loading rate and t is the time after the earthquake, plotted as a function of the number of M > 1.5 aftershocks N a ( t ) that have occurred by time t is bilinear with the slope of the fit to the first half of the aftershock sequence less than the slope of the fit to the second half. Thus, the aftershock seismicity rate is not proportional to the stress relaxation rate τ ′( t ) − τ ′ l over the entire sequence, but rather exhibits two distinct proportionalities. The observed postearthquake accumulation of M > 1.5 earthquakes in the aftershock zone as a function of time can be explained with the rate‐and‐state friction relation proposed by Dieterich (1994) between the cumulative number of earthquakes and τ ( t ).

California↗
USGS2010

Crustal movement investigations

Studies of horizontal crustal movement using conventional geodetic methods have been considerably expanded in the quadrennium 1971–1974. The basic fault monitor Geodimeter network now covers most of the major faults in California as well as the zone of faulting that extends into Nevada. Isolated Geodimeter networks in seismic areas of Montana, New Mexico, Utah, and Washington are also monitored. Part of the Geodimeter network along the San Andreas fault has been monitored for over 15 yr, and there appears to be a systematic deviation from a linear trend for most of the lines [ Greensfelder and Bennett , 1973; Savage et al. , 1973]. However, this deviation may be an artifact of a change in survey procedures in mid‐1969. The conclusions that can be drawn from the Geodimeter observations at present are the following: (1) The present apparent interplate motion across the San Andreas fault in central California is only 30–40 mm/yr [ Savage and Burford , 1973], in contrast to the average of 50–60 mm/yr as estimated from magnetic anomalies at the mouth of the Gulf of California. (2) Measurements of fault creep on the creeping segment of the San Andreas fault agree reasonably well with the plate movement indicated by geodetic measurements, this agreement suggesting that fault creep on this segment is the principal mode of accommodation [ Savage and Burford , 1971]. (3) The rate of strain accumulation along the San Andreas fault is not well measured anywhere, but the overall tensor strain rate appears to be about 0.3 microstrain/yr or less [ Savage et al. , 1973]. A considerable effort has also been expended reanalyzing old triangulation data for networks that cross the San Andreas fault. In such a study, Thatcher [1974] deduced evidence for several meters of afterslip on the San Andreas fault at depths greater than 10 km following the 1906 earthquake. This appears at the surface as a very rapid accumulation of strain in the years immediately following the earthquake. Meade [1974] has called attention to a remarkably uniform change of the astronomic azimuth with time for a line that crosses the San Andreas fault about 35 km south of the southern terminus of the 1906 rupture. This change of azimuth implies a uniform right lateral motion across the fault amounting to 32 mm/yr in the period 1885–1962. By comparing a 1942 triangulation survey and a 1970 Geodimeter survey, Page [1972] found neither lateral slip nor shear strain accumulation consistent with lateral slip across the Denali fault in Alaska. He did observe a north‐south extension that was possibly an effect of strain release at the time of the 1964 Alaska earthquake.

Reviews of Geophysics↗
USGS2007

Postseismic relaxation associated with transient creep rheology

[1] Perfettini and Avouac (2004) postulated that both the aftershock rate (assumed proportional to the local stressing rate) and the postseismic relaxation are driven by the loading imposed by postseismic slip on the brittle creep fault zone (BCFZ), the downdip extension of the fault zone below the coseismic rupture. I explore the consequences of that hypothesis for a long, strike-slip fault in the case where the BCFZ rheology is compatible with ordinary transient creep (creep strain proportional to log e (1 + t / τ 2 )). Because the important relaxation occurs near the bottom of the coseismic rupture, I calculate the postearthquake response with a model in which the BCFZ is represented by a viscoelastic half-space below the coseismic rupture. I find that both the predicted postseismic relaxation and the cumulative number of aftershocks can be approximated by the same temporal dependence N MO ( t ) = a MO (1−(1 + t / τ ) 1− p )/( p − 1), where t is the time after the earthquake and a MO , τ , and p are the constants chosen to fit either data set. Notice that d N MO ( t )/d t = ( a MO / τ )/(1 + t / τ ) p is the modified Omori law used to describe the rate of aftershock occurrence. Thus, the modified Omori law can be understood as a consequence of the Perfettini–Avouac hypothesis (aftershocks driven by slip on the BCFZ) and a BCFZ rheology compatible with ordinary transient creep. Moreover, the temporal dependence N MO ( t ) has been shown to fit postseismic surface deformation following at least 9 earthquakes. I also show that the conventional, one-dimensional, spring-block model of a BFCZ with a rheology compatible with ordinary transient creep leads to the same temporal dependence ( N MO ( t )).

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Consequences of viscous drag beneath a transform fault

A transform fault is modeled as a vertical cut through an elastic layer (schizosphere) of thickness overlying a viscous substrate (plastosphere). We consider a steady transform motion accommodated in the schizosphere wholly by slip on the fault and in the plastosphere, insofar as possible, by viscous flow. For the case where the viscosity in the plastosphere is strain rate dependent but independent of temperature, the velocity solution in the plastosphere is θ/π, where is the slip rate on the fault in the schizosphere and and θ are the cylindrical coordinates with the origin at the bottom of the fault. The viscous stress is singular at the bottom of the fault ( = 0) and exceeds the brittle (frictional) strength for . Equating the brittle strength to the viscous stress defines the brittle–ductile boundary in the plastosphere as a function of and viscosity. The additional condition that must be small allows the viscosity to be estimated from . For small , the temperature‐independent solution is a valid approximation to the temperature‐dependent solution, and the relation between viscosity and should remain valid. From the temperature‐independent model, we estimate that self‐heating due to dissipation in the plastosphere for reasonable Earth parameters is less than ∼20°C.

Journal of Geophysical Research↗