Search USGS⌕ Search

Geology topics

J.C. Savage

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

At least 55 records · Page 3Linked to original sources

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↗

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↗

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↗

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↗

Is internal friction friction?

Mogi [1974] proposed a simple model of the incipient rupture surface to explain the Coulomb failure criterion. We show here that this model can plausibly be extended to explain the Mohr failure criterion. In Mogi's model the incipient rupture surface immediately before fracture consists of areas across which material integrity is maintained (intact areas) and areas across which it is not (cracks). The strength of the incipient rupture surface is made up of the inherent strength of the intact areas plus the frictional resistance to sliding offered by the cracked areas. Although the coefficient of internal friction (slope of the strength versus normal stress curve) depends upon both the frictional and inherent strengths, the phenomenon of internal friction can be identified with the frictional part. The curvature of the Mohr failure envelope is interpreted as a consequence of differences in damage (cracking) accumulated in prefailure loading at different confining pressures.

Geophysical Research Letters↗

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↗

The paradigm that failed

Advances in science are memorialized; diversion are soon forgotten. Yet it is a diversion that I commemorate here: the v p /v s anomaly, the so-called seismic velocity ratio, as a method of predicting the time, place, and magnitude of an impending earthquake.

Earthquakes & Volcanoes (USGS)↗

Recent crustal subsidence at Yellowstone Caldera, Wyoming

Following a period of net uplift at an average rate of 15??1 mm/year from 1923 to 1984, the east-central floor of Yellowstone Caldera stopped rising during 1984-1985 and then subsided 25??7 mm during 1985-1986 and an additional 35??7 mm during 1986-1987. The average horizontal strain rates in the northeast part of the caldera for the period from 1984 to 1987 were: {Mathematical expression}1 = 0.10 ?? 0.09 ??strain/year oriented N33?? E??9?? and {Mathematical expression}2 = 0.20 ?? 0.09 ??strain/year oriented N57?? W??9?? (extension reckoned positive). A best-fit elastic model of the 1985-1987 vertical and horizontal displacements in the eastern part of the caldera suggests deflation of a horizontal tabular body located 10??5 km beneath Le Hardys Rapids, i.e., within a deep hydrothermal system or within an underlying body of partly molten rhyolite. Two end-member models each explain most aspects of historical unrest at Yellowstone, including the recent reversal from uplift to subsidence. Both involve crystallization of an amount of rhyolitic magma that is compatible with the thermal energy requirements of Yellowstone's vigorous hydrothermal system. In the first model, injection of basalt near the base of the rhyolitic system is the primary cause of uplift. Higher in the magmatic system, rhyolite crystallizes and releases all of its magmatic volatiles into the shallow hydrothermal system. Uplift stops and subsidence starts whenever the supply rate of basalt is less than the subsidence rate produced by crystallization of rhyolite and associated fluid loss. In the second model, uplift is caused primarily by pressurization of the deep hydrothermal system by magmatic gas and brine that are released during crystallization of rhyolite and them trapped at lithostatic pressure beneath an impermeable self-sealed zone. Subsidence occurs during episodic hydrofracturing and injection of pore fluid from the deep lithostatic-pressure zone into a shallow hydrostatic-pressure zone. Heat input from basaltic intrusions is required to maintain Yellowstone's silicic magmatic system and shallow hydrothermal system over time scales longer than about 105 years, but for the historical time period crystallization of rhyolite can account for most aspects of unrest at Yellowstone, including seismicity, uplift, subsidence, and hydrothermal activity. ?? 1990 Springer-Verlag.

Bulletin of Volcanology↗

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↗

Deformation near the epicenter of the 1984 Round Valley, California, earthquake

A trilateration network extending from near Mammoth Lakes to Bishop, California, was resurveyed following the 23 November 1984, Round Valley earthquake (ML = 5.8). The network had previously been surveyed in 1982. Deformation apparently associated with the Round Valley earthquake was detected as well as deformation due to the expansion of a magma chamber 8 km beneath the resurgent dome in the Long Valley caldera and right-lateral slip on the uppermost 2 km of the 1983 rupture surface in the south moat of the caldera. The deformation associated with Round Valley earthquake suggests left-lateral slip on the north-northeasterly striking vertical plane defined by the aftershock hypocenters located by A. S. Ryall. The earthquake moment implied by the deformation is about 3.8·1017 N-m, a value equivalent to an earthquake magnitude ML = 5.7 in good agreement with the observed magnitude of 5.8. A 0.053 km3 expansion of the magma chamber and 0.32 m slip on the 1983 rupture surface in the 1982-1985 interval was also required to account for the observed deformation.

California↗

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↗

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↗

Tectonic state: its significance and characterization in the assessment of seismic effects associated with reservoir impounding

Any analysis of seismicity associated with the filling of large reservoirs requires an evaluation of the natural tectonic state in order to determine whether impoundment is the basic source, a mechanically unrelated companion feature, or a triggering stimulus of the observed seismicity. Several arguments indicate that the associated seismicity is usually a triggered effect. Among the elements of tectonic state considered here (existing fractures, accumulated elastic strain, and deformational style), deformational style is especially critical in forecasting the occurrence of impoundment-induced seismicity. The observational evidence indicates that seismicity associated with impounding generally occurs in areas that combine steeply dipping faults, relatively high strain rates, and either extensional or horizontal-shear strain. Simple physical arguments suggest: (1) that increased fluid pressures resulting from increased reservoir head should enhance the likelihood of seismic activity, whatever the tectonic environment; (2) that stress changes resulting from surface loading may increase the likelihood of crustal failure in areas of normal and transcurrent faulting, whereas they generally inhibit failure in areas of thrust faulting. Comparisons with other earthquake-producing artificial and natural processes (underground explosions, fluid injection, underground mining, fluid extraction, volcanic emissions) indicate that reservoir loading may similarly modify the natural tectonic state. Subsurface loading resulting from fluid extraction may be a particularly close analogue of reservoir loading; "seismotectonic" events associated with fluid extraction have been recognized in both seismically active and otherwise aseismic regions. Because the historic record of seismicity and surface faulting commonly is short in comparison with recurrence intervals of earthquake and fault-slip events, tectonic state is most reliably appraised through combined studies of historic seismicity and faulting, instrumentally measured strain, and the geological record, especially that of the Quaternary. Experience in California and elsewhere demonstrates that the character and activity of recognized faults can be assessed by means of: instrumental earthquake investigations, repeated geodetic measurements, written history, archeological studies, fault topography, and local stratigraphic relations. Where faults are less easily distinguished, appraisals of tectonic state may be based on both the regional seismicity and the regional history of vertical movement as shown by: repeated levelling and sea-level measurements, written history, archeologic investigations, terrace and shoreline deformation, and denudation and sedimentation studies.

Engineering Geology↗

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↗

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↗