Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Atmospheres”

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 415 records · Page 23Linked to original sources
USGS2003

Importance of a stochastic distribution of floods and erosion thresholds in the bedrock river incision problem

Fluvial erosion of bedrock occurs during occasional flood events when boundary shear stress exceeds a critical threshold to initiate incision. Therefore efforts to model the evolution of topography over long timescales should include an erosion threshold and should be driven by a stochastic distribution of erosive events. However, most bedrock incision models ignore the threshold as a second‐order detail. In addition, climate is poorly represented in most landscape evolution models, so the quantitative relationship between erosion rate and measurable climatic variables has been elusive. Here we show that the presence of an erosion threshold, when combined with a well‐constrained, probabilistic model of storm and flood occurrence, has first‐order implications for the dynamics of river incision in tectonically active areas. First, we make a direct calculation of the critical shear stress required to pluck bedrock blocks for a field site in New York. Second, we apply a recently proposed stochastic, threshold, bedrock incision model to a series of streams in California, with known tectonic and climatic forcing. Previous work in the area has identified a weak relationship between channel gradient or relief and rock uplift rate that is not easily explained by simpler detachment‐limited models. The results with the stochastic threshold model show that even low erosion thresholds, which are exceeded in steep channels during high‐frequency flood events, fundamentally affect the predicted relationship between gradient and uplift rate in steady state rivers, in a manner consistent with the observed topography. This correspondence between theory and data is, however, nonunique; models in which a thin alluvial cover may act to inhibit channel incision in the low uplift rate zone also provide plausible explanations for the observed topography. Third, we explore the broader implications of the stochastic threshold model to the development of fluvial topography in active tectonic settings. We suggest that continued field applications of geomorphic models, including physically meaningful thresholds and stochastic climate distributions, are required to advance our knowledge of interactions among surficial, climatic, and crustal processes.

California↗
USGS2002

Waveform inversion of oscillatory signatures in long-period events beneath volcanoes

The source mechanism of long-period (LP) events is examined using synthetic waveforms generated by the acoustic resonance of a fluid-filled crack. We perform a series of numerical tests in which the oscillatory signatures of synthetic LP waveforms are used to determine the source time functions of the six moment tensor components from waveform inversions assuming a point source. The results indicate that the moment tensor representation is valid for the odd modes of crack resonance with wavelengths 2 L / n , 2 W / n , n = 3, 5, 7, …, where L and W are the crack length and width, respectively. For the even modes with wavelengths 2 L / n , 2 W / n , n = 2, 4, 6, …, a generalized source representation using higher-order tensors is required, although the efficiency of seismic waves radiated by the even modes is expected to be small. We apply the moment tensor inversion to the oscillatory signatures of an LP event observed at Kusatsu-Shirane Volcano, central Japan. Our results point to the resonance of a subhorizontal crack located a few hundred meters beneath the summit crater lakes. The present approach may be useful to quantify the source location, geometry, and force system of LP events, and opens the way for moment tensor inversions of tremor.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Thermal infrared spectroscopy of experimentally shocked anorthosite and pyroxenite: Implications for remote sensing of Mars

The feldspar and pyroxene mineralogies on Mars revealed by the Thermal Emission Spectrometer (TES) on Mars Global Surveyor likely record a variety of shock effects, as suggested by petrologic analyses of the Martian meteorites and the abundance of impact craters on the planet's surface. To study the effects of shock pressures on thermal infrared spectra of these minerals, we performed shock recovery experiments on orthopyroxenite and anorthosite samples from the Stillwater Complex (Montana) over peak pressures from 17 to 63 GPa. We acquired emissivity and hemispherical reflectance spectra (350–1400 cm −1 ; ∼7–29 μm) of both coherent chips and fine-grained powders of shocked and unshocked samples. These spectra are more directly comparable to remotely sensed data of Mars (e.g., TES) than previously acquired absorption or transmission spectra of shocked minerals. The spectra of experimentally shocked feldspar show systematic changes with increasing pressure due to depolymerization of the silica tetrahedra. For the spectra of chips, this includes the disappearance of small bands in the 500–650 cm −1 region and a strong band at 1115 cm −1 , and changes in positions of a strong band near 940 cm −1 and the Christiansen feature near 1250 cm −1 . Spectra of the shocked powders show the gradual disappearance of a transparency feature near 830 cm −1 . Fewer changes are observed in the pyroxene spectra at pressures as high as 63 GPa. Spectra of experimentally shocked minerals will help identify more precisely the mineralogy of rocks and soils not only from TES but also from Mars instruments such as miniTES and THEMIS.

Journal of Geophysical Research E: Planets↗
USGS2003

Slip history and dynamic implications of the 1999 Chi-Chi, Taiwan, earthquake

[1] We investigate the rupture process of the 1999 Chi‐Chi, Taiwan, earthquake using extensive near‐source observations, including three‐component velocity waveforms at 36 strong motion stations and 119 GPS measurements. A three‐plane fault geometry derived from our previous inversion using only static data [ Ji et al. , 2001 ] is applied. The slip amplitude, rake angle, rupture initiation time, and risetime function are inverted simultaneously with a recently developed finite fault inverse method that combines a wavelet transform approach with a simulated annealing algorithm [ Ji et al. , 2002b ]. The inversion results are validated by the forward prediction of an independent data set, the teleseismic P and SH ground velocities, with notable agreement. The results show that the total seismic moment release of this earthquake is 2.7 × 10 20 N m and that most of the slip occurred in a triangular‐shaped asperity involving two fault segments, which is consistent with our previous static inversion. The rupture front propagates with an average rupture velocity of ∼2.0 km s −1 , and the average slip duration (risetime) is 7.2 s. Several interesting observations related to the temporal evolution of the Chi‐Chi earthquake are also investigated, including (1) the strong effect of the sinuous fault plane of the Chelungpu fault on spatial and temporal variations in slip history, (2) the intersection of fault 1 and fault 2 not being a strong impediment to the rupture propagation, and (3) the observation that the peak slip velocity near the surface is, in general, higher than on the deeper portion of the fault plane, as predicted by dynamic modeling.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Stress-induced, time-dependent fracture closure at hydrothermal conditions

Time-dependent closure of fractures in quartz was measured in situ at 22–530°C temperature and 0.1-150 MPa water pressure. Unlike previous crack healing and rock permeability studies, in this study, fracture aperture is monitored directly and continuously using a windowed pressure vessel, a long-working-distance microscope, and reflected-light interferometry. Thus the fracture volume and geometry can be measured as a function of time, temperature, and water pressure. Relatively uniform closure occurs rapidly at temperatures and pressures where quartz becomes significantly soluble in water. During closure the aperture is reduced by as much as 80% in a few hours. We infer that this closure results from the dissolution of small particles or asperities that prop the fracture open. The driving force for closure via dissolution of the prop is the sum of three chemical potential terms: (1) the dissolution potential, proportional to the logarithm of the degree of undersaturation of the solution; (2) the coarsening potential, proportional to the radius of curvature of the prop; and (3) the pressure solution potential, proportional to the effective normal stress at the contact between propping particles and the fracture wall. Our observations suggest that closure is controlled by a pressure solution-like process. The aperture of dilatant fractures and microcracks in the Earth that are similar to those in our experiments, such as ones generated from thermal stressing or brittle failure during earthquake rupture and slip, will decrease rapidly with time, especially if the macroscopic stress is nonhydrostatic.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Three-dimensional velocity structure of crust and upper mantle in southwestern China and its tectonic implications

Using P and S arrival times from 4625 local and regional earthquakes recorded at 174 seismic stations and associated geophysical investigations, this paper presents a three‐dimensional crustal and upper mantle velocity structure of southwestern China (21°–34°N, 97°–105°E). Southwestern China lies in the transition zone between the uplifted Tibetan plateau to the west and the Yangtze continental platform to the east. In the upper crust a positive velocity anomaly exists in the Sichuan Basin, whereas a large‐scale negative velocity anomaly exists in the western Sichuan Plateau, consistent with the upper crustal structure under the southern Tibetan plateau. The boundary between these two anomaly zones is the Longmen Shan Fault. The negative velocity anomalies at 50‐km depth in the Tengchong volcanic area and the Panxi tectonic zone appear to be associated with temperature and composition variations in the upper mantle. The Red River Fault is the boundary between the positive and negative velocity anomalies at 50‐km depth. The overall features of the crustal and the upper mantle structures in southwestern China are a low average velocity, large crustal thickness variations, the existence of a high‐conductivity layer in the crust or/and upper mantle, and a high heat flow value. All these features are closely related to the collision between the Indian and the Asian plates.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Basin-centered asperities in great subduction zone earthquakes: A link between slip, subsidence, and subduction erosion?

Published areas of high coseismic slip, or asperities, for 29 of the largest Circum-Pacific megathrust earthquakes are compared to forearc structure revealed by satellite free-air gravity, bathymetry, and seismic profiling. On average, 71% of an earthquake's seismic moment and 79% of its asperity area occur beneath the prominent gravity low outlining the deep-sea terrace; 57% of an earthquake's asperity area, on average, occurs beneath the forearc basins that lie within the deep-sea terrace. In SW Japan, slip in the 1923, 1944, 1946, and 1968 earthquakes was largely centered beneath five forearc basins whose landward edge overlies the 350??C isotherm on the plate boundary, the inferred downdip limit of the locked zone. Basin-centered coseismic slip also occurred along the Aleutian, Mexico, Peru, and Chile subduction zones but was ambiguous for the great 1964 Alaska earthquake. Beneath intrabasin structural highs, seismic slip tends to be lower, possibly due to higher temperatures and fluid pressures. Kilometers of late Cenozoic subsidence and crustal thinning above some of the source zones are indicated by seismic profiling and drilling and are thought to be caused by basal subduction erosion. The deep-sea terraces and basins may evolve not just by growth of the outer arc high but also by interseismic subsidence not recovered during earthquakes. Basin-centered asperities could indicate a link between subsidence, subduction erosion, and seismogenesis. Whatever the cause, forearc basins may be useful indicators of long-term seismic moment release. The source zone for Cascadia's 1700 A.D. earthquake contains five large, basin-centered gravity lows that may indicate potential asperities at depth. The gravity gradient marking the inferred downdip limit to large coseismic slip lies offshore, except in northwestern Washington, where the low extends landward beneath the coast. Transverse gravity highs between the basins suggest that the margin is seismically segmented and could produce a variety of large earthquakes. Published in 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Finite element model predictions of static deformation from dislocation sources in a subduction zone: Sensitivities to homogeneous, isotropic, Poisson-solid, and half-space assumptions

Dislocation models can simulate static deformation caused by slip along a fault. These models usually take the form of a dislocation embedded in a homogeneous, isotropic, Poisson-solid half-space (HIPSHS). However, the widely accepted HIPSHS assumptions poorly approximate subduction zone systems of converging oceanic and continental crust. This study uses three-dimensional finite element models (FEMs) that allow for any combination (including none) of the HIPSHS assumptions to compute synthetic Green's functions for displacement. Using the 1995 Mw = 8.0 Jalisco-Colima, Mexico, subduction zone earthquake and associated measurements from a nearby GPS array as an example, FEM-generated synthetic Green's functions are combined with standard linear inverse methods to estimate dislocation distributions along the subduction interface. Loading a forward HIPSHS model with dislocation distributions, estimated from FEMs that sequentially relax the HIPSHS assumptions, yields the sensitivity of predicted displacements to each of the HIPSHS assumptions. For the subduction zone models tested and the specific field situation considered, sensitivities to the individual Poisson-solid, isotropy, and homogeneity assumptions can be substantially greater than GPS. measurement uncertainties. Forward modeling quantifies stress coupling between the Mw = 8.0 earthquake and a nearby Mw = 6.3 earthquake that occurred 63 days later. Coulomb stress changes predicted from static HIPSHS models cannot account for the 63-day lag time between events. Alternatively, an FEM that includes a poroelastic oceanic crust, which allows for postseismic pore fluid pressure recovery, can account for the lag time. The pore fluid pressure recovery rate puts an upper limit of 10-17 m2 on the bulk permeability of the oceanic crust. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Magma supply dynamics at Westdahl volcano, Alaska, modeled from satellite radar interferometry

A group of satellite radar interferograms that span the time period from 1991 to 2000 shows that Westdahl volcano, Alaska, deflated during its 1991-1992 eruption and is reinflating at a rate that could produce another eruption within the next several years. The rates of inflation and deflation are approximated by exponential decay functions having time constants of about 6 years and a few days, respectively. This behavior is consistent with a deep, constant-pressure magma source connected to a shallow reservoir by a magma-filled conduit. An elastic deformation model indicates that the reservoir is located about 6 km below sea level and beneath Westdahl Peak. We propose that the magma flow rate through the conduit is governed by the pressure gradient between the deep source and the reservoir. The pressure gradient, and hence the flow rate, are greatest immediately after eruptions. Pressurization of the reservoir decreases both the pressure gradient and the flow rate, but eventually the reservoir ruptures and an eruption or intrusion ensues. The eruption rate is controlled partly by the pressure gradient between the reservoir and surface, and therefore it, too, decreases with time. When the supply of eruptible magma is exhausted, the eruption stops, the reservoir begins to repressurize at a high rate, and the cycle repeats. This model might also be appropriate for other frequently active volcanoes with stable magma sources and relatively simple magma storage systems.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Interaction of lateral baroclinic forcing and turbulence in an estuary

Observations of density and velocity in a channel in northern San Francisco Bay show that the onset of vertical density stratification during flood tides is controlled by the balance between the cross-channel baroclinic pressure gradient and vertical mixing due to turbulence. Profiles of velocity, salinity, temperature, and suspended sediment concentration were measured in transects across Suisun Cutoff, in northern San Francisco Bay, on two days over the 12.5-hour tidal cycle. During flood tides an axial density front developed between fresher water flowing from the shallows of Grizzly Bay into the northern side of Suisun Cutoff and saltier water flowing up the channel. North of the front, transverse currents were driven by the lateral salinity gradient, with a top-to-bottom velocity difference greater than 30 cm/s. South of the front, the secondary circulation was weak, and along-channel velocities were greater than to the north. The gradient Richardson number shows that stratification was stable north of the front, while the water column was turbulently mixed south of the front. Time-series measurements of velocity and salinity demonstrate that the front develops during each tidal cycle. In estuaries, longitudinal dynamics predict less stratification during flood than ebb tides. These data show that stratification can develop during flood tides due to a lateral baroclinic pressure gradient in estuaries with complex bathymetry.

Journal of Geophysical Research C: Oceans↗
USGS2003

Resurfacing history of the northern plains of Mars based on geologic mapping of Mars Global Surveyor data

Geologic mapping of the northern plains of Mars, based on Mars Orbiter Laser Altimeter topography and Viking and Mars Orbiter Camera images, reveals new insights into geologic processes and events in this region during the Hesperian and Amazonian Periods. We propose four successive stages of lowland resurfacing likely related to the activity of near-surface volatiles commencing at the highland-lowland boundary (HLB) and progressing to lower topographic levels as follows (highest elevations indicated): Stage 1, upper boundary plains, Early Hesperian, <-2.0 to -2.9 km; Stage 2, lower boundary plains and outflow channel dissection, Late Hesperian, <-2.7 to -4.0 km; Stage 3, Vastitas Borealis Formation (VBF) surface, Late Hesperian to Early Amazonian, <-3.1 to -4.1 km; and Stage 4, local chaos zones, Early Amazonian, <-3.8 to -5.0 km. At Acidalia Mensa, Stage 2 and 3 levels may be lower (<-4.4 and -4.8 km, respectively). Contractional ridges form the dominant structure in the plains and developed from near the end of the Early Hesperian to the Early Amazonian. Geomorphic evidence for a northern-plains-filling ocean during Stage 2 is absent because one did not form or its evidence was destroyed by Stage 3 resurfacing. Remnants of possible Amazonian dust mantles occur on top of the VBF. The north polar layered deposits appear to be made up of an up to kilometer-thick lower sequence of sandy layers Early to Middle Amazonian in age overlain by Late Amazonian ice-rich dust layers; both units appear to have outliers, suggesting that they once were more extensive.

Journal of Geophysical Research E: Planets↗
USGS2003

Orogenic tectonism on Io

We catalog 143 Ionian mountains (montes) and mountain‐like features (mensae, tholi, plana, and small peaks) in order to investigate orogenic tectonism on Io. From this comprehensive list, we select 96 mountains for which there are sufficient coverage and resolution to discern spatial relationships with surrounding geologic features. Three of the 96 mountains are probably volcanoes, 92 appear to be tectonic massifs, and 1 is ambiguous. Of the 92 tectonic mountains, 38 abut paterae (volcanic or volcano‐tectonic craters with irregular or scalloped margins). This juxtaposition is unlikely to be a coincidence as the probability of it occurring by chance is ∼0.1%. We propose instead that orogenic faults may act as conduits for magma ascent, thus fueling patera formation near mountains. As resurfacing buries a shell of material from Io's surface to the base of the lithosphere, its effective radius is reduced and it heats up. We calculate the lithospheric volume change due to subsidence and thermal expansion as a function of lithospheric thickness. Conservation of volume dictates that this material must be uplifted at Io's surface. By estimating the total volume of the mountains, we are able to place a lower limit of 12 km on Io's lithospheric thickness. We hypothesize that, in some cases, mountain formation may be facilitated by asthenospheric diapirs impinging on the base of the lithosphere. The resulting lithospheric swell could focus the compressive stresses that drive orogenic tectonism. This model is one of several possible mechanisms for uplifting isolated mountains such as are observed on Io.

Journal of Geophysical Research E: Planets↗
USGS2003

Effects of spectrometer band pass, sampling, and signal‐to‐noise ratio on spectral identification using the Tetracorder algorithm

[1] Estimates of spectrometer band pass, sampling interval, and signal‐to‐noise ratio required for identification of pure minerals and plants were derived using reflectance spectra convolved to AVIRIS, HYDICE, MIVIS, VIMS, and other imaging spectrometers. For each spectral simulation, various levels of random noise were added to the reflectance spectra after convolution, and then each was analyzed with the Tetracorder spectral identification algorithm [ Clark et al. , 2003]. The outcome of each identification attempt was tabulated to provide an estimate of the signal‐to‐noise ratio at which a given percentage of the noisy spectra were identified correctly. Results show that spectral identification is most sensitive to the signal‐to‐noise ratio at narrow sampling interval values but is more sensitive to the sampling interval itself at broad sampling interval values because of spectral aliasing, a condition when absorption features of different materials can resemble one another. The band pass is less critical to spectral identification than the sampling interval or signal‐to‐noise ratio because broadening the band pass does not induce spectral aliasing. These conclusions are empirically corroborated by analysis of mineral maps of AVIRIS data collected at Cuprite, Nevada, between 1990 and 1995, a period during which the sensor signal‐to‐noise ratio increased up to sixfold. There are values of spectrometer sampling and band pass beyond which spectral identification of materials will require an abrupt increase in sensor signal‐to‐noise ratio due to the effects of spectral aliasing. Factors that control this threshold are the uniqueness of a material's diagnostic absorptions in terms of shape and wavelength isolation, and the spectral diversity of the materials found in nature and in the spectral library used for comparison. Array spectrometers provide the best data for identification when they critically sample spectra. The sampling interval should not be broadened to increase the signal‐to‐noise ratio in a photon‐noise‐limited system when high levels of accuracy are desired. It is possible, using this simulation method, to select optimum combinations of band‐pass, sampling interval, and signal‐to‐noise ratio values for a particular application that maximize identification accuracy and minimize the volume of imaging data.

Journal of Geophysical Research: Planets↗
USGS2004

Morphometric properties of Martian volcanoes

Mars Orbiter Laser Altimeter (MOLA) data have been used to construct Digital Elevation Models (DEM) of the Martian volcanoes in order to determine height, flank slope, caldera depth, and volumes. Summit elevations range from 21.1 km to −0.5 km, and relief varies from 1.0 km to almost 22 km. Average flank slopes are in the range of <1° to ∼10°, consistent with basaltic shield volcanism. The very low slopes of highland patera are also consistent with pyroclastic volcanism. Minimum volumes range from <10 12 to 10 15 m 3 . Estimates of the time required to build these volcanoes, on the basis of long-term terrestrial eruption rates, range from hundreds of thousands to tens of millions of years.

Journal of Geophysical Research E: Planets↗
USGS2004

Posteruption suspended sediment transport at Mount St. Helens: Decadal‐scale relationships with landscape adjustments and river discharges

Widespread landscape disturbance by the cataclysmic 1980 eruption of Mount St. Helens abruptly increased sediment supply in surrounding watersheds. The magnitude and duration of the redistribution of sediment deposited by the eruption as well as decades‐ to centuries‐old sediment remobilized from storage have varied chiefly with the style of disturbance. Posteruption suspended sediment transport has been greater and more persistent from zones of channel disturbance than from zones of hillslope disturbance. Despite the severe landscape disturbances caused by the eruption, relationships between discharge magnitudes and frequencies and suspended sediment transport have been remarkably consistent. Discharges smaller than mean annual flows generally have transported <5%, but locally ∼15%, of the annual suspended sediment loads, and infrequent ( p < 0.01), large floods have transported as much as 50% of the annual suspended sediment loads in a single day. However, moderate‐magnitude discharges (those greater than mean annual flows but less than 2‐year floods) have transported the greatest amounts of sediment from all disturbance zones. Such discharges have transported, on average, 60% to ∼95% of the annual suspended sediment loads, usually within cumulative periods of 1–3 weeks each year. Although small‐magnitude and large‐magnitude discharges have locally and episodically transported considerable amounts of suspended sediment, there has not been any notable change in the overall nature of the effective discharges; moderate‐magnitude flows have been the predominant discharges responsible for transporting the majority of suspended sediment during 20 years of posteruption landscape adjustment.

Washington↗
USGS2004

Stress field variations in the Swiss Alps and the northern Alpine foreland derived from inversion of fault plane solutions

This study is devoted to a systematic analysis of the state of stress of the central European Alps and northern Alpine foreland in Switzerland based on focal mechanisms of 138 earthquakes with magnitudes between 1 and 5. The most robust feature of the results is that the azimuth of the minimum compressive stress, S 3 , is generally well constrained for all data subsets and always lies in the NE quadrant. However, within this quadrant, the orientation of S 3 changes systematically both along the structural strike of the Alpine chain and across it. The variation in stress along the mountain belt from NE to SW involves a progressive, counterclockwise rotation of S 3 and is most clear in the foreland, where it amounts to 45°–50°. This pattern of rotation is compatible with the disturbance to the stress field expected from the indentation of the Adriatic Block into the central European Plate, possibly together with buoyancy forces arising from the strongly arcuate structure of the Moho to the immediate west of our study area. Across the Alps, the variation in azimuth of S 3 is defined by a progressive, counterclockwise rotation of about 45° from the foreland in the north across the Helvetic domain to the Penninic nappes in the south and is accompanied by a change from a slight predominance of strike-slip mechanisms in the foreland to a strong predominance of normal faulting in the high parts of the Alps. The observed rotation can be explained by the perturbation of the large-scale regional stress by a local uniaxial deviatoric tension with a magnitude similar to that of the regional differential stress and with an orientation perpendicular to the strike of the Alpine belt. The tensile nature and orientation of this stress is consistent with the “spreading” stress expected from lateral density changes due to a crustal root beneath the Alps.

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

Stress triggering in thrust and subduction earthquakes and stress interaction between the southern San Andreas and nearby thrust and strike-slip faults

We argue that key features of thrust earthquake triggering, inhibition, and clustering can be explained by Coulomb stress changes, which we illustrate by a suite of representative models and by detailed examples. Whereas slip on surface-cutting thrust faults drops the stress in most of the adjacent crust, slip on blind thrust faults increases the stress on some nearby zones, particularly above the source fault. Blind thrusts can thus trigger slip on secondary faults at shallow depth and typically produce broadly distributed aftershocks. Short thrust ruptures are particularly efficient at triggering earthquakes of similar size on adjacent thrust faults. We calculate that during a progressive thrust sequence in central California the 1983 M w = 6.7 Coalinga earthquake brought the subsequent 1983 M w = 6.0 Nuñez and 1985 M w = 6.0 Kettleman Hills ruptures 10 bars and 1 bar closer to Coulomb failure. The idealized stress change calculations also reconcile the distribution of seismicity accompanying large subduction events, in agreement with findings of prior investigations. Subduction zone ruptures are calculated to promote normal faulting events in the outer rise and to promote thrust-faulting events on the periphery of the seismic rupture and its downdip extension. These features are evident in aftershocks of the 1957 M w = 9.1 Aleutian and other large subduction earthquakes. We further examine stress changes on the rupture surface imparted by the 1960 M w = 9.5 and 1995 M w = 8.1 Chile earthquakes, for which detailed slip models are available. Calculated Coulomb stress increases of 2–20 bars correspond closely to sites of aftershocks and postseismic slip, whereas aftershocks are absent where the stress drops by more than 10 bars. We also argue that slip on major strike-slip systems modulates the stress acting on nearby thrust and strike-slip faults. We calculate that the 1857 M w = 7.9 Fort Tejon earthquake on the San Andreas fault and subsequent interseismic slip brought the Coalinga fault ∼1 bar closer to failure but inhibited failure elsewhere on the Coast Ranges thrust faults. The 1857 earthquake also promoted failure on the White Wolf reverse fault by 8 bars, which ruptured in the 1952 M w = 7.3 Kern County shock but inhibited slip on the left-lateral Garlock fault, which has not ruptured since 1857. We thus contend that stress transfer exerts a control on the seismicity of thrust faults across a broad spectrum of spatial and temporal scales.

California↗