USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research: Earth Surface”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11Linked to original sources
USGS2017

Rock friction under variable normal stress

This study is to determine the detailed response of shear strength and other fault properties to changes in normal stress at room temperature using dry initially bare rock surfaces of granite at normal stresses between 5 and 7 MPa. Rapid normal stress changes result in gradual, approximately exponential changes in shear resistance with fault slip. The characteristic length of the exponential change is similar for both increases and decreases in normal stress. In contrast, changes in fault normal displacement and the amplitude of small high-frequency elastic waves transmitted across the surface follow a two stage response consisting of a large immediate and a smaller gradual response with slip. The characteristic slip distance of the small gradual response is significantly smaller than that of shear resistance. The stability of sliding in response to large step decreases in normal stress is well predicted using the shear resistance slip length observed in step increases. Analysis of the shear resistance and slip-time histories suggest nearly immediate changes in strength occur in response to rapid changes in normal stress; these are manifested as an immediate change in slip speed. These changes in slip speed can be qualitatively accounted for using a rate-independent strength model. Collectively, the observations and model show that acceleration or deceleration in response to normal stress change depends on the size of the change, the frictional characteristics of the fault surface, and the elastic properties of the loading system.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

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

Crystallographic controls on the frictional behavior of dry and water-saturated sheet structure minerals

We compare the frictional strengths of 17 sheet structure mineral powders, measured under dry and water-saturated conditions, to identify the factors that cause many of them to be relatively weak. The dry coefficient of friction μ ranges upward from 0.2 for graphite, leveling off at 0.8 for margarite, clintonite, gibbsite, kaolinite, and lizardite. The values of μ (dry) correlate directly with calculated (001) interlayer bond strengths of the minerals. This correlation occurs because shear becomes localized along boundary and Riedel shears and the platy minerals in them rotate into alignment with the shear planes. For those gouges with μ (dry) < 0.8, shear occurs by breaking the interlayer bonds to form new cleavage surfaces. Where μ (dry) = 0.8, consistent with Byerlee's law, the interlayer bonds are sufficiently strong that other frictional processes dominate. The transition in dry friction mechanisms corresponds to calculated surface energies of 2–3 J/m 2 . Adding water causes μ to decrease for every mineral tested except graphite. If the minerals are separated into groups with similar crystal structures, μ (wet) increases with increasing interlayer bond strength within each group. This relationship also holds for the swelling clay montmorillonite, whose water-saturated strength is consistent with the strengths of nonswelling clays of similar crystal structure. Water in the saturated gouges forms thin, structured films between the plate surfaces. The polar water molecules are bonded to the plate surfaces in proportion to the mineral's surface energy, and μ (wet) reflects the stresses required to shear through the water films.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

The strength and rheology of methane clathrate hydrate

Methane clathrate hydrate (structure I) is found to be very strong, based on laboratory triaxial deformation experiments we have carried out on samples of synthetic, high-purity, polycrystalline material. Samples were deformed in compressional creep tests (i.e., constant applied stress, ??), at conditions of confining pressure P = 50 and 100 MPa, strain rate 4.5 ?? 10-8 ??? ?? ??? 4.3 ?? 10-4 s-1, temperature 260 ??? T ??? 287 K, and internal methane pressure 10 ??? PCH4 ??? 15 MPa. At steady state, typically reached in a few percent strain, methane hydrate exhibited strength that was far higher than expected on the basis of published work. In terms of the standard high-temperature creep law, ?? = A??ne-(E*+PV*)/RT the rheology is described by the constants A = 108.55 MPa-n s-1, n = 2.2, E* = 90,000 J mol-1, and V* = 19 cm3 mol-1. For comparison at temperatures just below the ice point, methane hydrate at a given strain rate is over 20 times stronger than ice, and the contrast increases at lower temperatures. The possible occurrence of syntectonic dissociation of methane hydrate to methane plus free water in these experiments suggests that the high strength measured here may be only a lower bound. On Earth, high strength in hydrate-bearing formations implies higher energy release upon decomposition and subsequent failure. In the outer solar system, if Titan has a 100-km-thick near-surface layer of high-strength, low-thermal conductivity methane hydrate as has been suggested, its interior is likely to be considerably warmer than previously expected.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Usefulness of weak bands in midinfrared remote sensing of particulate planetary surfaces

Midinfrared (2.5–25 μm) reflectance spectra of minerals are often used to predict emittance qualitatively. These spectra display weak overtone and combination tone bands, which may be as diagnostic of composition as the strong fundamental molecular vibration bands usually considered for remote sensing applications, but which have been widely ignored. However, unlike the strong bands, the contrast of weak bands relative to the continuum in the midinfrared usually does not decrease with decreasing particle size, but typically increases. To illustrate this behavior, transmittance and reflectance spectra of calcite and quartz are presented for the wavelength range from 4000 to 400 cm −1 (2.5–25 μm). It is the purpose of this paper to point out that these weak bands are potentially useful for compositional remote sensing of particulate planetary surfaces; this will require completion of supporting laboratory studies to document the occurrence of weak bands.

Journal of Geophysical Research Solid Earth↗
USGS1997

Timing of detachment faulting in the Bullfrog Hills and Bare Mountain area, southwest Nevada: Inferences from 40Ar/39Ar, K-Ar, U-Pb, and fission track thermochronology

Crustal extension in the Bullfrog Hills and Bare Mountain area of southwest Nevada is associated with movement along a regional detachment fault. Normal faulting in the upper plate and rapid cooling (denudation) of the lower plate were coeval with Miocene silicic volcanism and with west-northwest transport along the detachment fault. A west-northwest progression of tilting along upper plate normal faults is indicated by ages of the volcanic rocks in relation to angular unconformities. Near the breakaway, tilting in the upper plate occurred between 12.7 and 11.6 Ma, continued less strongly past 10.7 Ma, and was over by 8.2 Ma. Ten to 20 km west of the breakaway, tilting occurred between 10.7 and 10.33 Ma, continued less strongly after 10.33 Ma, and was over by 8.1 Ma. The cooling histories of the lower plate metamorphic rocks were determined by thermochronologic dating methods: K-Ar and 40 Ar/ 39 Ar on muscovite, biotite, and hornblende, 40 Ar/ 39 Ar on K-feldspar, U-Pb on apatite, zircon, and sphene, and fission track on apatite, zircon, and sphene. Lower plate rocks 10 km west of the breakaway cooled slowly from Early Cretaceous lower-amphibolite facies conditions through 350±50° to 300±50°C between 57 and 38 Ma, then cooled rapidly from 205±50° to 120±50°C between 12.6±1.6 and 11.1±1.9 Ma. Lower plate rocks 20 km west of the breakaway cooled slowly from Early Cretaceous upper-amphibolite facies conditions through 500±50°C at 78–67 Ma, passed through 350±50° to 300±50°C between 16.3±0.4 and 10.5±0.3 Ma, then cooled rapidly from 285±50° to 120±50°C between 10.2 and 8.6 Ma. Upper plate tilting and rapid cooling (denudation) of the lower plate occurred simultaneously in the respective areas. The early slow-cooling part of the lower plate thermal histories was probably related to erosion at the Earth's surface, which stripped off about 9 km of material in 50 to 100 m.y. The results indicate an initial fault dip ≥30° and a 12 mm yr −1 west-northwest migration of the locus of rapid tilting in the upper plate.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Columbia Glacier, Alaska: Changes in velocity 1977-1986

The Columbia Glacier, a grounded, iceberg-calving tidewater glacier near Valdez, Alaska, began to retreat about 1977. Drastic retreat occurred in 1984, and by early 1986, retreat amounted to 2 km. The glacier has thinned more than 100 m since 1974 at a point 4 km behind the 1974 terminus position. Between 1977 and 1985 the lower glacier ice velocity increased from 3–8 m/d to 10–15 m/d. Ice velocity in the region 0.5 km above the terminus was highest near the time the glacier was most receded (late fall), and lowest near the time of maximum length (early summer), for years 1977–1982. Velocity in the region 52–57 from the head of the glacier was highest in mid-spring, and lowest in early fall from 1977 to 1985. Through the years 1983–1985, the dates of maximum and minimum velocities within 0.5 km of the receding terminus tended toward the dates of the 52–57 km maximum and minimums. This occurred because as the terminus receded, it was no longer strongly influenced by the reverse slope of the terminal moraine shoal. Velocities near the terminus fluctuated by 2–3 m/d during summer and fall, when liquid water input was variable, and were relatively constant during winter. Hourly variations in ice velocities are controlled by liquid water input to the glacier hydraulic system and tide stage. Velocity increases near periods of high surface water input and decreases during periods of high tide as a result of hydrostatic back pressure.

Journal of Geophysical Research Solid Earth↗
USGS2012

On the state of lithospheric stress in the absence of applied tectonic forces

Numerous published analyses of the nontectonic state of stress are based on Hooke's law and the boundary condition of zero horizontal deformation. This approach has been used to determine the gravitational stress state as well as the effects of processes such as erosion and temperature changes on the state of lithospheric stress. The major disadvantage of these analyses involves the assumption of lateral constraint which seems unrealistic in view of the observational fact that the crust can deform horizontally in response to applied loads. If the same problems are addressed by assuming that the remote stress state is constant, instead of the condition of zero horizontal deformation, then the resulting stress states are entirely different and in good accord with observations. In the absence of applied tectonic forces the only likely gravitational stress states are those for which all three principal stresses are nearly equal. To the contrary, the gravitational stress states developed on the basis of the lateral constraint assumption can be ruled out. The processes of erosion and sedimentation have slight tendencies to increase and decrease, respectively, the state of deviatoric stress. In particular, for initial stress states in the range of slightly extensional to compressional, erosion has the effect of enhancing the ratio of average horizontal to vertical stress, which may explain, at least in part, the common observation of high near-surface horizontal stresses. Temperature changes have only minor effects on the stress state, as averaged over the thickness of the lithosphere.

Journal of Geophysical Research Solid Earth↗
USGS2005

Applying petrophysical models to radar travel time and electrical resistivity tomograms: Resolution-dependent limitations

[1] Geophysical imaging has traditionally provided qualitative information about geologic structure; however, there is increasing interest in using petrophysical models to convert tomograms to quantitative estimates of hydrogeologic, mechanical, or geochemical parameters of interest (e.g., permeability, porosity, water content, and salinity). Unfortunately, petrophysical estimation based on tomograms is complicated by limited and variable image resolution, which depends on (1) measurement physics (e.g., electrical conduction or electromagnetic wave propagation), (2) parameterization and regularization, (3) measurement error, and (4) spatial variability. We present a framework to predict how core‐scale relations between geophysical properties and hydrologic parameters are altered by the inversion, which produces smoothly varying pixel‐scale estimates. We refer to this loss of information as “correlation loss.” Our approach upscales the core‐scale relation to the pixel scale using the model resolution matrix from the inversion, random field averaging, and spatial statistics of the geophysical property. Synthetic examples evaluate the utility of radar travel time tomography (RTT) and electrical‐resistivity tomography (ERT) for estimating water content. This work provides (1) a framework to assess tomograms for geologic parameter estimation and (2) insights into the different patterns of correlation loss for ERT and RTT. Whereas ERT generally performs better near boreholes, RTT performs better in the interwell region. Application of petrophysical models to the tomograms in our examples would yield misleading estimates of water content. Although the examples presented illustrate the problem of correlation loss in the context of near‐surface geophysical imaging, our results have clear implications for quantitative analysis of tomograms for diverse geoscience applications.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Historic surface slip along the San Andreas Fault near Parkfield, California

The Parkfield Earthquake Prediction Experiment is focusing close attention on the 44-km-long section of the San Andreas fault that last ruptured seismically in 1966 ( M s 6.0). The 20-km-long central segment of the 1966 Parkfield rupture, extending from the mainshock epicenter at Middle Mountain southeastward to Gold Hill, forms a 1- to 2-km salient northeastward away from the dominant N40°W strike. Following the 1966 earthquake afterslip, aseismic slip has been nearly constant. Moderate Parkfield earthquakes have recurred on average every 21 years since 1857, when a great earthquake ( M ≈8) ruptured at least as far north as the southern Parkfield segment. Many measurements of slip have been made near Parkfield since 1966. Nevertheless, much of the history of surface slip remained uncertain, especially the total amount associated with the 1966 event. In 1985 we measured accumulated slip on the four oldest cultural features offset by the fault along the 1966 Parkfield rupture segment. We interpret net slip on each feature as a sum of event slip (sum of coseismic and rapid preseismic and postseismic slip) from Parkfield earthquakes and steady interseismic slip as measured over the last 20 years on nearby alinement arrays, creep meters, and trilateration lines. We assumed for each site that event slip was identical for the 1922, 1934, and 1966 Parkfield events and that long-term average rate of interseismic slip was constant between all events. Two fences on the southern segment, southeast of Gold Hill, indicate event slip of 13 and 15 cm and interseismic slip rate of 0.36 and 0.30 cm/yr since 1959 and 1908, respectively. At these sites, redundant independent data support our assumption that both event and interseismic slip occur uniformly. On the central segment, near Parkfield, both the 1934 and the 1966 ruptures offset a bridge built in 1932. Interseismic slip rate near the bridge has been about 1.1 cm/yr since 1966; thus we deduce an average event slip of 31 cm for the 1934 and 1966 earthquakes. On a parallel fault trace, 1 km to the southwest, slip was about 8 cm in 1966; thus total event slip summed across the entire fault zone near Parkfield was nearly 40 cm. On Middle Mountain, 4 km north of the 1966 mainshock epicenter, an offset fence indicates 17 cm of slip in 1966 and a 2.26-cm/yr interseismic slip rate since circa 1946. Thus the central segment of the 1966 rupture is characterized by much larger event slip (∼40 cm) than both distal segments (∼15 cm). This amount of surface slip per event is about twice what had been previously assumed. Larger 1966 surface slip in the central part of the rupture is geodetically compatible with a coseismic slip of 65±10 cm slip on a narrow, buried asperity between Middle Mountain and Gold Hill that has been inferred from the depth distribution of early aftershocks. Assuming our characteristic surface slip model, one can further deduce a deficit in slip since the great 1857 earthquake. Taking the long-term slip rate as 3.3 cm/yr, the surface slip deficit is 3±0.2 m south of Gold Hill but only 0.3±0.3 m northward from Parkfield.

Journal of Geophysical Research Solid Earth↗
USGS2003

New insights into Kilauea's volcano dynamics brought by large-scale relative relocation of microearthquakes

We investigated the microseismicity recorded in an active volcano to infer information concerning the volcano structure and long-term dynamics, by using relative relocations and focal mechanisms of microearthquakes. There were 32,000 earthquakes of the Mauna Loa and Kilauea volcanoes recorded by more than eight stations of the Hawaiian Volcano Observatory seismic network between 1988 and 1999. We studied 17,000 of these events and relocated more than 70%, with an accuracy ranging from 10 to 500 m. About 75% of these relocated events are located in the vicinity of subhorizontal decollement planes, at a depth of 8-11 km. However, the striking features revealed by these relocation results are steep southeast dipping fault planes working as reverse faults, clearly located below the decollement plane and which intersect it. If this decollement plane coincides with the pre-Mauna Loa seafloor, as hypothesized by numerous authors, such reverse faults rupture the pre-Mauna Loa oceanic crust. The weight of the volcano and pressure in the magma storage system are possible causes of these ruptures, fully compatible with the local stress tensor computed by Gillard et al. [1996]. Reverse faults are suspected of producing scarps revealed by kilometer-long horizontal slip-perpendicular lineations along the decollement surface and therefore large-scale roughness, asperities, and normal stress variations. These are capable of generating stick-slip, large-magnitude earthquakes, the spatial microseismic pattern observed in the south flank of Kilauea volcano, and Hilina-type instabilities. Rupture intersecting the decollement surface, causing its large-scale roughness, may be an important parameter controlling the growth of Hawaiian volcanoes.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Coseismic deformation of the Wrights tunnel during the 1906 San Francisco earthquake: A key to understanding 1906 fault slip and 1989 surface ruptures in the southern Santa Cruz Mountains, California

The Wrights tunnel is an abandoned railroad tunnel that crosses the San Andreas fault in the southern Santa Cruz Mountains in the vicinity of the 1989 Loma Prieta earthquake. The tunnel was damaged and deformed during the 1906 San Francisco earthquake and a plot showing postearthquake measurements made in the tunnel is given by Lawson [1908]. The amount of offset shown on this plot (1.5 m) has been used in several studies as being representative of the amount of fault offset along this segment of the San Andreas fault in 1906. Our historical research shows that different observers reported different amounts of fault offset in the tunnel and that the 1.5 m given on the plot is not a surveyed measurement. In addition, the plot of the tunnel has been interpreted in several previous studies as evidence of a broad (1.5 km) zone of faulting beneath Summit Ridge. Our analysis shows that this plot need not indicate a broad zone of deformation. Our historical research and modeling of the tunnel measurements indicate that faulting was confined to a zone less than 400 m wide and that 60–85% of the coseismic slip occurred across a single fault plane. There is no evidence for offset across a second shear zone beneath Summit Ridge in 1906. This implies that surface fractures reported on Summit Ridge in 1906 were not associated with significant deformation of the tunnel, implying that they were shallow, surficial features. By analogy, the very similar fractures that occurred on Summit Ridge in 1989 were also probably the result of shallow gravitational, rather than deep-seated tectonic, processes. Our modeling also indicates that total coseismic, near-surface slip across the San Andreas fault zone in the Wrights tunnel in 1906 was at least 1.7–1.8 m.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Coseismic slip distribution of the 1923 Kanto earthquake, Japan

The slip distribution associated with the 1923 M = 7.9 Kanto, Japan, earthquake is reexamined in light of new data and modeling. We utilize a combination of first-order triangulation, second-order triangulation, and leveling data in order to constrain the coseismic deformation. The second-order triangulation data, which have not been utilized in previous studies of 1923 coseismic deformation, are associated with only slightly smaller errors than the first-order triangulation data and expand the available triangulation data set by about a factor of 10. Interpretation of these data in terms of uniform-slip models in a companion study by Nyst et al. shows that a model involving uniform coseismic slip on two distinct rupture planes explains the data very well and matches or exceeds the fit obtained by previous studies, even one which involved distributed slip. Using the geometry of the Nyst et al. two-plane slip model, we perform inversions of the same geodetic data set for distributed slip. Our preferred model of distributed slip on the Philippine Sea plate interface has a moment magnitude of 7.86. We find slip maxima of ???8-9 m beneath Odawara and ???7-8 m beneath the Miura peninsula, with a roughly 2:1 ratio of strike-slip to dip-slip motion, in agreement with a previous study. However, the Miura slip maximum is imaged as a more broadly extended feature in our study, with the high-slip region continuing from the Miura peninsula to the southern Boso peninsula region. The second-order triangulation data provide good evidence for ???3 m right-lateral strike slip on a 35-km-long splay structure occupying the volume between the upper surface of the descending Philippine Sea plate and the southern Boso peninsula. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Mechanical properties of sand, silt, and clay containing tetrahydrofuran hydrate

The mechanical behavior of hydrate-bearing sediments subjected to large strains has relevance for the stability of the seafloor and submarine slopes, drilling and coring operations, and the analysis of certain small-strain properties of these sediments (for example, seismic velocities). This study reports on the results of comprehensive axial compression triaxial tests conducted at up to 1 MPa confining pressure on sand, crushed silt, precipitated silt, and clay specimens with closely controlled concentrations of synthetic hydrate. The results show that the stress-strain behavior of hydrate-bearing sediments is a complex function of particle size, confining pressure, and hydrate concentration. The mechanical properties of hydrate-bearing sediments at low hydrate concentration (probably < 40% of pore space) appear to be determined by stress-dependent soil stiffness and strength. At high hydrate concentrations (>50% of pore space), the behavior becomes more independent of stress because the hydrates control both stiffness and strength and possibly the dilative tendency of sediments by effectively increasing interparticle coordination, cementing particles together, and filling the pore space. The cementation contribution to the shear strength of hydrate-bearing sediments decreases with increasing specific surface of soil minerals. The lower the effective confining stress, the greater the impact of hydrate formation on normalized strength.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Earthquake fracture energy inferred from kinematic rupture models on extended faults

We estimate fracture energy on extended faults for several recent earthquakes by retrieving dynamic traction evolution at each point on the fault plane from slip history imaged by inverting ground motion waveforms. We define the breakdown work (Wb) as the excess of work over some minimum traction level achieved during slip. Wb is equivalent to "seismological" fracture energy (G) in previous investigations. Our numerical approach uses slip velocity as a boundary condition on the fault. We employ a three-dimensional finite difference algorithm to compute the dynamic traction evolution in the time domain during the earthquake rupture. We estimate Wb by calculating the scalar product between dynamic traction and slip velocity vectors. This approach does not require specifying a constitutive law and assuming dynamic traction to be collinear with slip velocity. If these vectors are not collinear, the inferred breakdown work depends on the initial traction level. We show that breakdown work depends on the square of slip. The spatial distribution of breakdown work in a single earthquake is strongly correlated with the slip distribution. Breakdown work density and its integral over the fault, breakdown energy, scale with seismic moment according to a power law (with exponent 0.59 and 1.18, respectively). Our estimates of breakdown work range between 4 ?? 105 and 2 ?? 107 J/m2 for earthquakes having moment magnitudes between 5.6 and 7.2. We also compare our inferred values with geologic surface energies. This comparison might suggest that breakdown work for large earthquakes goes primarily into heat production. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2013

InSAR Evidence for an active shallow thrust fault beneath the city of Spokane Washington, USA

In 2001, a nearly five month long sequence of shallow, mostly small magnitude earthquakes occurred beneath the city of Spokane, a city with a population of about 200,000, in the state of Washington. During most of the sequence, the earthquakes were not well located because seismic instrumentation was sparse. Despite poor-quality locations, the earthquake hypocenters were likely very shallow, because residents near the city center both heard and felt many of the earthquakes. The combination of poor earthquake locations and a lack of known surface faults with recent movement make assessing the seismic hazards related to the earthquake swarm difficult. However, the potential for destruction from a shallow moderate-sized earthquake is high, for example Christchurch New Zealand in 2011, so assessing the hazard potential of a seismic structure involved in the Spokane earthquake sequence is important. Using interferometric synthetic aperture radar (InSAR) data from the European Space Agency ERS2 and ENVISAT satellites and the Canadian Space Agency RADARSAT-1, satellite we are able to show that slip on a shallow previously unknown thrust fault, which we name the Spokane Fault, is the source of the earthquake sequence. The part of the Spokane Fault that slipped during the 2001 earthquake sequence underlies the north part of the city, and slip on the fault was concentrated between ~0.3 and 2 km depth. Projecting the buried fault plane to the surface gives a possible surface trace for the Spokane Fault that strikes northeast from the city center into north Spokane.

Washington↗
USGS2002

Prehistoric ruptures of the Gurvan Bulag fault, Gobi Altay, Mongolia

The 1957 Gobi Altay M8.3 earthquake in southern Mongolia was associated with the simultaneous rupture of several faults, including the Gurvan Bulag reverse fault, which is located about 25 km south of the main strike-slip Bogd fault. Our study of paleoseismic excavations across the Gurvan Bulag fault suggests that the penultimate surface rupture occurred after 6.0 ka, most likely between 2.6 and 4.4 ka, and a possible earlier rupture occurred after 7.3 ka. Our interpretation of the stratigraphic relations in one of the exposures suggests that at least five earthquakes have generated surface rupture of the Gurvan Bulag fault since the abandonment of an ancient alluvial fan surface. Luminescence dating of sediment associated with this surface indicates that it formed either 26.6 ± 2.1 ka or 16.1 ± 2.0 ka. These data imply that the recurrence intervals for surface faulting on the Gurvan Bulag and Bogd faults are similar, on the order of several thousands of years, but that the penultimate surface ruptures of the two faults did not occur during the same earthquake.

Gobi Altay↗
USGS1996

Frictional behavior of large displacement experimental faults

The coefficient of friction and velocity dependence of friction of initially bare surfaces and 1-mm-thick simulated fault gouges (<90 μm) of Westerly granite were determined as a function of displacement to >400 mm at 25°C and 25 MPa normal stress. Steady state negative friction velocity dependence and a steady state fault zone microstructure are achieved after ∼18 mm displacement, and an approximately constant strength is reached after a few tens of millimeters of sliding on initially bare surfaces. Simulated fault gouges show a large but systematic variation of friction, velocity dependence of friction, dilatancy, and degree of localization with displacement. At short displacement (<10 mm), simulated gouge is strong, velocity strengthening and changes in sliding velocity are accompanied by relatively large changes in dilatancy rate. With continued displacement, simulated gouges become progressively weaker and less velocity strengthening, the velocity dependence of dilatancy rate decreases, and deformation becomes localized into a narrow basal shear which at its most localized is observed to be velocity weakening. With subsequent displacement, the fault restrengthens, returns to velocity strengthening, or to velocity neutral, the velocity dependence of dilatancy rate becomes larger, and deformation becomes distributed. Correlation of friction, velocity dependence of friction and of dilatancy rate, and degree of localization at all displacements in simulated gouge suggest that all quantities are interrelated. The observations do not distinguish the independent variables but suggest that the degree of localization is controlled by the fault strength, not by the friction velocity dependence. The friction velocity dependence and velocity dependence of dilatancy rate can be used as qualitative measures of the degree of localization in simulated gouge, in agreement with previous studies. Theory equating the friction velocity dependence of simulated gouge to the sum of the friction velocity dependence of bare surfaces and the velocity dependence of dilatancy rate of simulated gouge fails to quantitatively account for the experimental observations.

Journal of Geophysical Research B: Solid Earth↗