Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Geophysics Research”

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 757 records · Page 42Linked to original sources

Bathymetry and absorptivity of Titan's Ontario Lacus

Ontario Lacus is the largest and best characterized lake in Titan's south polar region. In June and July 2009, the Cassini RADAR acquired its first Synthetic Aperture Radar (SAR) images of the area. Together with closest approach altimetry acquired in December 2008, these observations provide a unique opportunity to study the lake's nearshore bathymetry and complex refractive properties. Average radar backscatter is observed to decrease exponentially with distance from the local shoreline. This behavior is consistent with attenuation through a deepening layer of liquid and, if local topography is known, can be used to derive absorptive dielectric properties. Accordingly, we estimate nearshore topography from a radar altimetry profile that intersects the shoreline on the East and West sides of the lake. We then analyze SAR backscatter in these regions to determine the imaginary component of the liquid's complex index of refraction (κ). The derived value, κ = (6.1 −1.3 +1.7 ) × 10 −4 , corresponds to a loss tangent of tan Δ = (9.2 −2.0 +2.5 ) × 10 −4 and is consistent with a composition dominated by liquid hydrocarbons. This value can be used to test compositional models once the microwave optical properties of candidate materials have been measured. In areas that do not intersect altimetry profiles, relative slopes can be calculated assuming the index of refraction is constant throughout the liquid. Accordingly, we construct a coarse bathymetry map for the nearshore region by measuring bathymetric slopes for eleven additional areas around the lake. These slopes vary by a factor of ∼5 and correlate well with observed shoreline morphologies.

Journal of Geophysical Research E: Planets

A rain splash transport equation assimilating field and laboratory measurements

Process-based models of hillslope evolution require transport equations relating sediment flux to its major controls. An equation for rain splash transport in the absence of overland flow was constructed by modifying an approach developed by Reeve (1982) and parameterizing it with measurements from single-drop laboratory experiments and simulated rainfall on a grassland in East Africa. The equation relates rain splash to hillslope gradient, the median raindrop diameter of a storm, and ground cover density; the effect of soil texture on detachability can be incorporated from other published results. The spatial and temporal applicability of such an equation for rain splash transport in the absence of overland flow on uncultivated hillslopes can be estimated from hydrological calculations. The predicted transport is lower than landscape-averaged geologic erosion rates from Kenya but is large enough to modify short, slowly eroding natural hillslopes as well as microtopographic interrill surfaces between which overland flow transports the mobilized sediment. Copyright 2010 by the American Geophysical Union. Copyright 2010 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface

Flow structures and sandbar dynamics in a canyon river during a controlled flood, Colorado River, Arizona

In canyon rivers, debris fan constrictions create rapids and downstream pools characterized by secondary flow structures that are closely linked to channel morphology. In this paper we describe detailed measurements of the three‐dimensional flow structure and sandbar dynamics of two pools along the Colorado River in the Grand Canyon during a controlled flood release from Glen Canyon Dam. Results indicate that the pools are characterized by large lateral recirculation zones (eddies) resulting from flow separation downstream from the channel constrictions, as well as helical flow structures in the main channel and eddy. The lateral recirculation zones are low‐velocity areas conducive to fine sediment deposition, particularly in the vicinity of the separation and reattachment points and are thus the dominant flow structures controlling sandbar dynamics. The helical flow structures also affect morphology but appear secondary in importance to the lateral eddies. During the controlled flood, sandbars in the separation and reattachment zones at both sites tended to build gradually during the rising limb and peak flow. Deposition in shallow water on the sandbars was accompanied by erosion in deeper water along the sandbar slope at the interface with the main channel. Erosion occurred via rapid mass failures as well as by gradual boundary shear stress driven processes. The flow structures and morphologic links at our study sites are similar to those identified in other river environments, in particular sharply curved meanders and channel confluences where the coexistence of lateral recirculation and helical flows has been documented.

Journal of Geophysical Research F: Earth Surface

The perfect debris flow? Aggregated results from 28 large-scale experiments

Aggregation of data collected in 28 controlled experiments reveals reproducible debris-flow behavior that provides a clear target for model tests. In each experiment ∼10 m 3 of unsorted, water-saturated sediment composed mostly of sand and gravel discharged from behind a gate, descended a steep, 95-m flume, and formed a deposit on a nearly horizontal runout surface. Experiment subsets were distinguished by differing basal boundary conditions (1 versus 16 mm roughness heights) and sediment mud contents (1 versus 7 percent dry weight). Sensor measurements of evolving flow thicknesses, basal normal stresses, and basal pore fluid pressures demonstrate that debris flows in all subsets developed dilated, coarse-grained, high-friction snouts, followed by bodies of nearly liquefied, finer-grained debris. Mud enhanced flow mobility by maintaining high pore pressures in flow bodies, and bed roughness reduced flow speeds but not distances of flow runout. Roughness had these effects because it promoted debris agitation and grain-size segregation, and thereby aided growth of lateral levees that channelized flow. Grain-size segregation also contributed to development of ubiquitous roll waves, which had diverse amplitudes exhibiting fractal number-size distributions. Despite the influence of these waves and other sources of dispersion, the aggregated data have well-defined patterns that help constrain individual terms in a depth-averaged debris-flow model. The patterns imply that local flow resistance evolved together with global flow dynamics, contradicting the hypothesis that any consistent rheology applied. We infer that new evolution equations, not new rheologies, are needed to explain how characteristic debris-flow behavior emerges from the interactions of debris constituents.

Journal of Geophysical Research F: Earth Surface

Iceberg calving as a primary source of regional‐scale glacier‐generated seismicity in the St. Elias Mountains, Alaska

Since the installation of the Alaska Regional Seismic Network in the 1970s, data analysts have noted nontectonic seismic events thought to be related to glacier dynamics. While loose associations with the glaciers of the St. Elias Mountains have been made, no detailed study of the source locations has been undertaken. We performed a two-step investigation surrounding these events, beginning with manual locations that guided an automated detection and event sifting routine. Results from the manual investigation highlight characteristics of the seismic waveforms including single-peaked (narrowband) spectra, emergent onsets, lack of distinct phase arrivals, and a predominant cluster of locations near the calving termini of several neighboring tidewater glaciers. Through these locations, comparison with previous work, analyses of waveform characteristics, frequency-magnitude statistics and temporal patterns in seismicity, we suggest calving as a source for the seismicity. Statistical properties and time series analysis of the event catalog suggest a scale-invariant process that has no single or simple forcing. These results support the idea that calving is often a response to short-lived or localized stress perturbations. Our results demonstrate the utility of passive seismic instrumentation to monitor relative changes in the rate and magnitude of iceberg calving at tidewater glaciers that may be volatile or susceptible to ensuing rapid retreat, especially when existing seismic infrastructure can be used.

Alaska

Testing statistical self-similarity in the topology of river networks

Recent work has demonstrated that the topological properties of real river networks deviate significantly from predictions of Shreve's random model. At the same time the property of mean self-similarity postulated by Tokunaga's model is well supported by data. Recently, a new class of network model called random self-similar networks (RSN) that combines self-similarity and randomness has been introduced to replicate important topological features observed in real river networks. We investigate if the hypothesis of statistical self-similarity in the RSN model is supported by data on a set of 30 basins located across the continental United States that encompass a wide range of hydroclimatic variability. We demonstrate that the generators of the RSN model obey a geometric distribution, and self-similarity holds in a statistical sense in 26 of these 30 basins. The parameters describing the distribution of interior and exterior generators are tested to be statistically different and the difference is shown to produce the well-known Hack's law. The inter-basin variability of RSN parameters is found to be statistically significant. We also test generator dependence on two climatic indices, mean annual precipitation and radiative index of dryness. Some indication of climatic influence on the generators is detected, but this influence is not statistically significant with the sample size available. Finally, two key applications of the RSN model to hydrology and geomorphology are briefly discussed.

Journal of Geophysical Research F: Earth Surface

Fluorescence-based proxies for lignin in freshwater dissolved organic matter

Lignin phenols have proven to be powerful biomarkers in environmental studies; however, the complexity of lignin analysis limits the number of samples and thus spatial and temporal resolution in any given study. In contrast, spectrophotometric characterization of dissolved organic matter (DOM) is rapid, noninvasive, relatively inexpensive, requires small sample volumes, and can even be measured in situ to capture fine-scale temporal and spatial detail of DOM cycling. Here we present a series of cross-validated Partial Least Squares models that use fluorescence properties of DOM to explain up to 91% of lignin compositional and concentration variability in samples collected seasonally over 2 years in the Sacramento River/San Joaquin River Delta in California, United States. These models were subsequently used to predict lignin composition and concentration from fluorescence measurements collected during a diurnal study in the San Joaquin River. While modeled lignin composition remained largely unchanged over the diurnal cycle, changes in modeled lignin concentrations were much greater than expected and indicate that the sensitivity of fluorescence-based proxies for lignin may prove invaluable as a tool for selecting the most informative samples for detailed lignin characterization. With adequate calibration, similar models could be used to significantly expand our ability to study sources and processing of DOM in complex surface water systems.

Journal of Geophysical Research G: Biogeosciences

Projected changes in atmospheric heating due to changes in fire disturbance and the snow season in the western Arctic, 2003–2100

In high latitudes, changes in climate impact fire regimes and snow cover duration, altering the surface albedo and the heating of the regional atmosphere. In the western Arctic, under four scenarios of future climate change and future fire regimes (2003–2100), we examined changes in surface albedo and the related changes in regional atmospheric heating due to: (1) vegetation changes following a changing fire regime, and (2) changes in snow cover duration. We used a spatially explicit dynamic vegetation model (Alaskan Frame-based Ecosystem Code) to simulate changes in successional dynamics associated with fire under the future climate scenarios, and the Terrestrial Ecosystem Model to simulate changes in snow cover. Changes in summer heating due to the changes in the forest stand age distributions under future fire regimes showed a slight cooling effect due to increases in summer albedo (mean across climates of −0.9 W m −2 decade −1 ). Over this same time period, decreases in snow cover (mean reduction in the snow season of 4.5 d decade −1 ) caused a reduction in albedo, and a heating effect (mean across climates of 4.3 W m −2 decade −1 ). Adding both the summer negative change in atmospheric heating due to changes in fire regimes to the positive changes in atmospheric heating due to changes in the length of the snow season resulted in a 3.4 W m −2 decade −1 increase in atmospheric heating. These findings highlight the importance of gaining a better understanding of the influences of changes in surface albedo on atmospheric heating due to both changes in the fire regime and changes in snow cover duration.

Journal of Geophysical Research: Biogeosciences

Source water controls on the character and origin of dissolved organic matter in streams of the Yukon River basin, Alaska

Climate warming and permafrost degradation at high latitudes will likely impact watershed hydrology, and consequently, alter the concentration and character of dissolved organic carbon (DOC) in northern rivers. We examined seasonal variation of DOC chemistry in 16 streams of the Yukon River basin, Alaska. Our primary objective was to evaluate the relationship between source water (shallow versus deep groundwater flow paths) and DOC chemical composition. Using base cation chemistry and principal component analysis, we observed high contributions of deep groundwater to glacial and clearwater streams, whereas blackwater streams received larger contributions from shallow groundwater sources. DOC concentration and specific ultraviolet absorbance peaked during spring snowmelt in all streams, and were consistently higher in blackwater streams than in glacial and clearwater streams. The hydrophobic acid fraction of DOC dominated across all streams and seasons, comprising between 35% and 56% of total DOC. The hydrophilic acid fraction of DOC was more prominent in glacial (23% ± 3%) and clearwater streams (19% ± 1%) than in blackwater streams (16% ± 1%), and was enriched during winter base flow (29% ± 1%) relative to snowmelt and summer base flow. We observed that an increase in the contribution of deep groundwater to streamflow resulted in decreased DOC concentration, aromaticity, and DOC-to-dissolved organic nitrogen ratio, and an increase in the proportion of hydrophilic acids relative to hydrophobic acids. Our findings suggest that future permafrost degradation and higher contributions of groundwater to streamflow may result in a higher fraction of labile DOM in streams of the Yukon basin.

Journal of Geophysical Research G: Biogeosciences

Controls on mangrove forest‐atmosphere carbon dioxide exchanges in western Everglades National Park

We report on net ecosystem production (NEP) and key environmental controls on net ecosystem exchange (NEE) of carbon dioxide (CO 2 ) between a mangrove forest and the atmosphere in the coastal Florida Everglades. An eddy covariance system deployed above the canopy was used to determine NEE during January 2004 through August 2005. Maximum daytime NEE ranged from −20 to −25 μ mol (CO 2 ) m −2 s −1 between March and May. Respiration (R d ) was highly variable (2.81 ± 2.41 μ mol (CO 2 ) m −2 s −1 ), reaching peak values during the summer wet season. During the winter dry season, forest CO 2 assimilation increased with the proportion of diffuse solar irradiance in response to greater radiative transfer in the forest canopy. Surface water salinity and tidal activity were also important controls on NEE. Daily light use efficiency was reduced at high (>34 parts per thousand (ppt)) compared to low (<17 ppt) salinity by 46%. Tidal inundation lowered daytime R d by ∼0.9 μ mol (CO 2 ) m −2 s −1 and nighttime R d by ∼0.5 μ mol (CO 2 ) m −2 s −1 . The forest was a sink for atmospheric CO 2 , with an annual NEP of 1170 ± 127 g C m −2 during 2004. This unusually high NEP was attributed to year‐round productivity and low ecosystem respiration which reached a maximum of only 3 g C m −2 d −1 . Tidal export of dissolved inorganic carbon derived from belowground respiration likely lowered the estimates of mangrove forest respiration. These results suggest that carbon balance in mangrove coastal systems will change in response to variable salinity and inundation patterns, possibly resulting from secular sea level rise and climate change.

Florida

A model‐data intercomparison of CO2 exchange across North America: Results from the North American Carbon Program site synthesis

Our current understanding of terrestrial carbon processes is represented in various models used to integrate and scale measurements of CO 2 exchange from remote sensing and other spatiotemporal data. Yet assessments are rarely conducted to determine how well models simulate carbon processes across vegetation types and environmental conditions. Using standardized data from the North American Carbon Program we compare observed and simulated monthly CO 2 exchange from 44 eddy covariance flux towers in North America and 22 terrestrial biosphere models. The analysis period spans ∼220 site‐years, 10 biomes, and includes two large‐scale drought events, providing a natural experiment to evaluate model skill as a function of drought and seasonality. We evaluate models' ability to simulate the seasonal cycle of CO 2 exchange using multiple model skill metrics and analyze links between model characteristics, site history, and model skill. Overall model performance was poor; the difference between observations and simulations was ∼10 times observational uncertainty, with forested ecosystems better predicted than nonforested. Model‐data agreement was highest in summer and in temperate evergreen forests. In contrast, model performance declined in spring and fall, especially in ecosystems with large deciduous components, and in dry periods during the growing season. Models used across multiple biomes and sites, the mean model ensemble, and a model using assimilated parameter values showed high consistency with observations. Models with the highest skill across all biomes all used prescribed canopy phenology, calculated NEE as the difference between GPP and ecosystem respiration, and did not use a daily time step.

Journal of Geophysical Research: Biogeosciences

On the resolution of shallow mantle viscosity structure using post-earthquake relaxation data: Application to the 1999 Hector Mine, California, earthquake

Most models of lower crust/mantle viscosity inferred from postearthquake relaxation assume one or two uniform-viscosity layers. A few existing models possess apparently significant radially variable viscosity structure in the shallow mantle (e.g., the upper 200 km), but the resolution of such variations is not clear. We use a geophysical inverse procedure to address the resolving power of inferred shallow mantle viscosity structure using postearthquake relaxation data. We apply this methodology to 9 years of GPS-constrained crustal motions after the 16 October 1999 M = 7.1 Hector Mine earthquake. After application of a differencing method to isolate the postearthquake signal from the “background” crustal velocity field, we find that surface velocities diminish from ∼20 mm/yr in the first few months to ≲2 mm/yr after 2 years. Viscoelastic relaxation of the mantle, with a time-dependent effective viscosity prescribed by a Burgers body, provides a good explanation for the postseismic crustal deformation, capturing both the spatial and temporal pattern. In the context of the Burgers body model (which involves a transient viscosity and steady state viscosity), a resolution analysis based on the singular value decomposition reveals that at most, two constraints on depth-dependent steady state mantle viscosity are provided by the present data set. Uppermost mantle viscosity (depth ≲ 60 km) is moderately resolved, but deeper viscosity structure is poorly resolved. The simplest model that explains the data better than that of uniform steady state mantle viscosity involves a linear gradient in logarithmic viscosity with depth, with a small increase from the Moho to 220 km depth. However, the viscosity increase is not statistically significant. This suggests that the depth-dependent steady state viscosity is not resolvably different from uniformity in the uppermost mantle.

California

Seismic swarm associated with the 2008 eruption of Kasatochi Volcano, Alaska: earthquake locations and source parameters

An energetic seismic swarm accompanied an eruption of Kasatochi Volcano in the central Aleutian volcanic arc in August of 2008. In retrospect, the first earthquakes in the swarm were detected about 1 month prior to the eruption onset. Activity in the swarm quickly intensified less than 48 h prior to the first large explosion and subsequently subsided with decline of eruptive activity. The largest earthquake measured as moment magnitude 5.8, and a dozen additional earthquakes were larger than magnitude 4. The swarm exhibited both tectonic and volcanic characteristics. Its shear failure earthquake features were b value = 0.9, most earthquakes with impulsive P and S arrivals and higher-frequency content, and earthquake faulting parameters consistent with regional tectonic stresses. Its volcanic or fluid-influenced seismicity features were volcanic tremor, large CLVD components in moment tensor solutions, and increasing magnitudes with time. Earthquake location tests suggest that the earthquakes occurred in a distributed volume elongated in the NS direction either directly under the volcano or within 5-10 km south of it. Following the M W 5.8 event, earthquakes occurred in a new crustal volume slightly east and north of the previous earthquakes. The central Aleutian Arc is a tectonically active region with seismicity occurring in the crusts of the Pacific and North American plates in addition to interplate events. We postulate that the Kasatochi seismic swarm was a manifestation of the complex interaction of tectonic and magmatic processes in the Earth's crust. Although magmatic intrusion triggered the earthquakes in the swarm, the earthquakes failed in context of the regional stress field.

Alaska

Analysis of nonvolcanic tremor on the San Andreas Fault near Parkfield, CA using U.S. Geological Survey Parkfield Seismic Array

Reports by Nadeau and Dolenc (2005) that tremor had been detected near Cholame Valley spawned an effort to use UPSAR (U. S. Geological Survey Parkfield Seismic Array) to study characteristics of tremor. UPSAR was modified to record three channels of velocity at 40–50 sps continuously in January 2005 and ran for about 1 month, during which time we recorded numerous episodes of tremor. One tremor, on 21 January at 0728, was recorded with particularly high signal levels as well as another episode 3 days later. Both events were very emergent, had a frequency content between 2 and 8 Hz, and had numerous high-amplitude, short-duration arrivals within the tremor signal. Here using the first episode as an example, we discuss an analysis procedure, which yields azimuth and apparent velocity of the tremor at UPSAR. We then provide locations for both tremor episodes. The emphasis here is how the tremor episode evolves. Twelve stations were operating at the time of recording. Slowness of arrivals was determined using cross correlation of pairs of stations; the same method used in analyzing the main shock data from 28 September 2004. A feature of this analysis is that 20 s of the time series were used at a time to calculate correlation; the longer windows resulted in more consistent estimates of slowness, but lower peak correlations. These values of correlation (peaks of about 0.25), however, are similar to that obtained for the S wave of a microearthquake. Observed peaks in slowness were traced back to source locations assumed to lie on the San Andreas fault. Our inferred locations for the two tremor events cluster near the locations of previously observed tremor, south of the Cholame Valley. Tremor source depths are in the 14–24 km range, which is below the seismogenic brittle zone, but above the Moho. Estimates of error do not preclude locations below the Moho, however. The tremor signal is very emergent but contains packets that are several times larger than the background tremor signal and lasts about 5 s. These impulsive wavelets are similar to low-frequency earthquakes signals seen in Japan but appear to be broader band rather than just higher in low-frequency energy. They may be more appropriately called high-energy tremor (HET). HET signals at UPSAR correlate well with the record of this event from station GHIB of the HRSN borehole array at Parkfield and HETs typically have a higher cross-correlation coefficient than the rest of the tremor event. The amplitudes of a large HET are consistent with a magnitude of 0.1 when compared with a M2.3 event that had about the same epicenter. Polarizations of the tremor episode at UPSAR are mostly just north of east. Both linearity and azimuth evolve over time suggesting a change in tremor source location over time and linearity is typically higher at the HETs.

California

Detecting aseismic strain transients from seismicity data

Aseismic deformation transients such as fluid flow, magma migration, and slow slip can trigger changes in seismicity rate. We present a method that can detect these seismicity rate variations and utilize these anomalies to constrain the underlying variations in stressing rate. Because ordinary aftershock sequences often obscure changes in the background seismicity caused by aseismic processes, we combine the stochastic Epidemic Type Aftershock Sequence model that describes aftershock sequences well and the physically based rate‐ and state‐dependent friction seismicity model into a single seismicity rate model that models both aftershock activity and changes in background seismicity rate. We implement this model into a data assimilation algorithm that inverts seismicity catalogs to estimate space‐time variations in stressing rate. We evaluate the method using a synthetic catalog, and then apply it to a catalog of M ≥ 1.5 events that occurred in the Salton Trough from 1990 to 2009. We validate our stressing rate estimates by comparing them to estimates from a geodetically derived slip model for a large creep event on the Obsidian Buttes fault. The results demonstrate that our approach can identify large aseismic deformation transients in a multidecade long earthquake catalog and roughly constrain the absolute magnitude of the stressing rate transients. Our method can therefore provide a way to detect aseismic transients in regions where geodetic resolution in space or time is poor.

Journal of Geophysical Research B: Solid Earth

High‐resolution locations of triggered earthquakes and tomographic imaging of Kilauea Volcano's south flank

The spatiotemporal patterns of seismicity beneath Kilauea's south flank give insight to the structure and geometry of the decollement on which large, tsunamigenic earthquakes have occurred, and its relation to slow slip events (SSEs), which have been observed every 1 to 2 years since 1997. In order to record earthquakes triggered by a SSE that was predicted to occur in March 2007, a temporary network of 20 seismometers was deployed on Kilauea's south flank, termed the SEQ network. While the SSE did not occur until 17 June 2007, theSEQ network recorded over 3000 earthquakes, including those triggered by the SSE. We relocate hypocenters of volcano‐tectonic earthquakes and invert for P and S wave velocity structure using waveform cross‐correlation and double‐difference tomography using data from the SEQ network and the permanent Hawaii Volcano Observatory network (HVO) data, with additional data from other previous temporary arrays. The best‐constrained hypocenters, recorded by both the SEQ and HVO networks, indicate the decollement as a subhorizontal layer of seismicity at 8 km depth less than 1 km thick in most areas, with the western portion of the decollement dipping to the southeast. The seismicity triggered by the June 2007 SSE includes over 400 earthquakes overlapping with the southern edge of the decollement seismicity. A shallower swarm of earthquakes also occurred between 2 and 7 km depth in April 2007 near Apua Point, and may have been indirectly triggered by the M w 8.1 Solomon Islands earthquake at ∼6000 km distance, which occurred 48 h prior to the beginning of the swarm.

Hawaii

Magma flux at Okmok Volcano, Alaska, from a joint inversion of continuous GPS, campaign GPS, and interferometric synthetic aperture radar

Volcano deformation is usually measured using satellite geodetic techniques including interferometric synthetic aperture radar (InSAR), campaign GPS, and continuous GPS. Differences in the spatial and temporal sampling of each system mean that most appropriate inversion scheme to determine the source parameters from each data set is different. Most studies either compare results from independent inversions or subsample the data sets to the lowest common factor. It is unclear whether differences in the solution reflect differences in source behavior, differences in measurement bias, or differences in inversion technique. Here we develop a single inversion procedure that captures the benefits of each system, especially the daily sampling of continuous GPS and the high spatial resolution of InSAR. Okmok Volcano, Alaska, is an ideal target for such a test because a long series (<15 years) of InSAR and continuous GPS measurement exists and the source is almost continuously active and in a stable location.

Alaska

Aftershocks are well aligned with the background stress field, contradicting the hypothesis of highly-heterogeneous crustal stress

It has been proposed that the crustal stress field contains small-length-scale heterogeneity of much larger amplitude than the uniform background stress. This model predicts that earthquake focal mechanisms should reflect the loading stress rather than the uniform background stress. So, if the heterogeneous stress hypothesis is correct, focal mechanisms before and after a large earthquake should align with the tectonic loading and the earthquake-induced static stress perturbation, respectively. However, I show that the off-fault triggered aftershocks of the 1992 M7.3 Landers, California, earthquake align with the same stress field as the pre-Landers mechanisms. The aftershocks occurred on faults that were well oriented for failure in the pre-Landers stress field and then loaded by the Landers-induced static stress change. Aftershocks in regions experiencing a 0.05 to 5 MPa coseismic differential stress change align with the modeled Landers-induced static stress change, implying that they were triggered by the stress perturbation. Contrary to the heterogeneous stress hypothesis, these triggered aftershocks are also well aligned with the pre-Landers stress field obtained from inverting the pre-Landers focal mechanisms. Therefore, the inverted pre-Landers stress must represent the persistent background stress field. Earthquake focal mechanisms provide an unbiased sample of the spatially coherent background stress field, which is large relative to any small-scale stress heterogeneity. The counterexample provided by the Landers earthquake is strong evidence that the heterogeneous stress model is not widely applicable.

California