SEARCH · USGS Science

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,387 records · Page 77Linked to original sources

USGS2012

Basin and range crustal and upper mantle structure, northwest to central Nevada

We present an interpretation of the crustal and uppermost mantle structure of the Basin and Range of northwestern Nevada based on seismic refraction/wide‐angle reflection, near‐vertical reflection, and gravity data. In comparison to most previous estimates, we find that the crust is somewhat thicker (32–36 km versus 22–30 km), and the uppermost mantle velocity is somewhat higher (8.0 km/s versus 7.3–7.9 km/s). Along our transects, the crust is thinnest (32 km) in the Carson Sink‐Buena Vista Valley region and increases by 2–4 km to the west and east, respectively. There is considerable complexity throughout the crust where velocities range from of 2.5 km/s at the surface to 7.4 km/s in the lowermost crust. Variations in velocity and structure of the upper crustal layers reveal apparent basement velocity depressions (areas of lower velocities extending up to 10 km in depth) that underlie some surface ranges as well as the basins. The middle crust rises from about 20 km beneath central Nevada to within 12 km of the surface beneath the area of thinnest crust and is characterized by a modest (∼0.1 km/s) change in velocity and low‐velocity gradients. These midcrustal layers mark the onset of high crustal reflectivity and the apparent limiting depth to which Basin and Range faults can be traced in near‐vertical reflection profiles, suggesting that these midcrustal layers represent the transition between the brittle and ductile zones of the crust. The lower crust is more structurally complex, with layers thickening and thinning in a systematic manner with the upper crustal layers; generally, where there are velocity depressions in the upper crust, the lower crust is thickest and shallowest. The geometry of these lower crustal layers (derived from refraction modeling) coincides with changes in the crustal reflectivity, determined from the Consortium of Continental Reflection Profiling reflection data. The lower crustal layer is unusually high in velocity (7.4 km/s) and is likely the layer identified as mantle in some previous studies. We do not identify the 7.4 km/s layer as mantle because (1) there is an underlying layer with a velocity (8.0 km/s) that is more consistent with the worldwide average velocity for the upper mantle, and (2) the 7.4 km/s layer does not correspond to the “reflection” Moho. Gravity modeling and comparison to existing seismic models show a general consensus in many aspects with respect to crustal structure. This new model forms the basis for speculation on some of the processes associated with rifting of the Basin and Range Province. One such process, lithospheric magmatism, is inferred from the strong attenuation of transmitted seismic waves, which occurs at the same interface at which high‐amplitude, bright spot reflections originate. Unlike previous models, the overall structure and velocity of the crust and uppermost mantle of our new model are similar to other regions worldwide which have undergone high degrees of extension.

Nevada↗

USGS2012

The velocity field along the San Andreas Fault in central and southern California

The velocity field within a 100‐km‐broad zone centered on the San Andreas fault between the Mexican border and San Francisco Bay has been inferred from repeated surveys of trilateration networks in the 1973–1989 interval. The velocity field has the appearance of a shear flow that remains parallel to the local strike of the fault even through such major deflections as the big bend of the San Andreas fault in the Transverse Ranges of southern California. Across‐strike profiles of the fault‐parallel component of velocity exhibit the expected sigmoidal shape, whereas across‐strike profiles of the fault‐normal component of velocity are flat and featureless. No significant convergence upon the fault is observed even along the big bend sector of the fault. Simple dislocation models can explain most of the features of the observed velocity field, but those explanations are not unique. About 35 mm/yr of relative plate motion is accounted for within the span of the trilateration networks. Geologic studies indicate that the secular slip rate on the San Andreas fault is about 35 mm/yr. The agreement between these two estimates implies that most of the strain accumulation is elastic and will be recovered in subsequent earthquakes. The relative motion observed across the San Andreas fault (35 mm/yr) plus that observed across the Eastern California shear zone (8 mm/yr) accounts for most (43 mm/yr) of the observed North America‐Pacific relative plate motion (47 mm/yr).

California↗

USGS2012

High-frequency spectral falloff of earthquakes, fractal dimension of complex rupture, b value, and the scaling of strength on faults

The high-frequency falloff ω −γ of earthquake displacement spectra and the b value of aftershock sequences are attributed to the character of spatially varying strength along fault zones. I assume that the high frequency energy of a main shock is produced by a self-similar distribution of subevents, where the number of subevents with radii greater than R is proportional to R −D , D being the fractal dimension. In this model, an earthquake is composed of a hierarchical set of smaller earthquakes. The static stress drop is parameterized to be proportional to R η , and strength is assumed to be proportional to static stress drop. I find that a distribution of subevents with D = 2 and stress drop independent of seismic moment (η = 0) produces a main shock with an ω −2 falloff, if the subevent areas fill the rupture area of the main shock. By equating subevents to “islands” of high stress of a random, self-similar stress field on a fault, I relate D to the scaling of strength on a fault, such that D = 2 − η. Thus D = 2 corresponds to constant stress drop scaling (η = 0) and scale-invariant fault strength. A self-similar model of aftershock rupture zones on a fault is used to determine the relationship between the b value, the size distribution of aftershock rupture zones, and the scaling of strength on a fault. The b value for aftershock sequences on a fault is found to equal (3 − 1.5η)/(3 + η). Therefore this model indicates that the typically observed spectral falloffs of ω −2 and b values of 1 can be entirely caused by scale-invariant strength (η = 0) along fault zones.

Journal of Geophysical Research↗

USGS2012

Crustal resistivity structure from magnetotelluric soundings in the Colorado Plateau-Basin and Range provinces, central and western Arizona

Resistivity structure to about 25 km depth is defined from two-dimensional modeling of 29 magnetotelluric (MT) soundings (0.002–5 Hz) that traverse 280 km of the southwestern Colorado Plateau, transition zone, and Basin and Range provinces in Arizona. From the surface to 5 km depth, the MT model suggests structural relationships between low-resistivity sedimentary and volcanic rocks (50–300 ohm m) and high-resistivity granitic and gneissic basement (500–9000 ohm m). In the Basin and Range province, the MT model and a seismic reflection section show a generally consistent distribution of supracrustal rocks that have relatively low to moderate resistivity (MT) and relatively strong, locally coherent reflectivity. The supracrustal zone defined by these physical properties is inferred to be composed of upper plate rocks above a middle Tertiary detachment fault system which has been mapped in surrounding ranges. Some low-angle fault zones inferred from seismic reflections to extend into high resistivity basement below the supracrustal rocks are not resolved by the MT model. A low-resistivity zone with a conductance of 500 S or more is modeled in the crust at an average depth of about 15 km in the Basin and Range province and transition zone and may deepen below the southwestern part of the Colorado Plateau. In the Basin and Range province, the top of the low resistivity may correspond to a reflective layer with a 6-s two-way-travel time. This deep low-resistivity zone might be caused by a small fraction of connected hydrous solutions or silicic melts.

Journal of Geophysical Research↗

USGS2012

On plate tectonics and the geologic evolution of southwestern North America

Very rapid subduction of the Farallon plate under southwestern North America between 60 and 40 Ma was accompanied by a relatively low volume of magmatism throughout the southwestern United States and northern Mexico. Between 40 and 20 Ma, when subduction slowed significantly and in one area may have even stopped, magmatism became widespread and voluminous from Nevada and Utah to central Mexico. This correlation of rapid subduction with a relatively low volume of magmatism can be explained by the observation that subduction-related andesitic arc volcanism, often formed in a Laramide-style compressional regime, is relatively low volume compared to continental volcanism. The shallow roots of arc volcanic systems are clearly exposed in the porphyry copper deposits found in currently active arcs and common throughout southwestern North America between 60 and 50 Ma. By 43 Ma, worldwide plate motions changed, the Pacific plate began moving away from North America, and subduction of the Farallon plate slowed. By around 36 Ma, the easternmost part of the East Pacific Rise, which was located between the Pioneer and Murray fracture zones, approached the trench and the young, hot, buoyant lithosphere appears to have clogged part of the subduction zone. Uplift on land became widespread. Voluminous continental magmatism formed the Sierra Madre Occidental (SMO) of Mexico, one of the largest batholiths in the world, as well as volcanic centers now exposed in the San Juan Mountains of Colorado and the Rio Grande Rift of New Mexico. Vectors of motion of the Pacific plate relative to the North American plate determined by Stock and Molnar (1988) are consistent with formation of a transtensional environment along the plate boundary sufficient to create a 100- to 200-km-wide void just landward of the old volcanic arc. While the SMO batholith was forming within this void, the Monterey and Arguello microplates just offshore to the west were broken off from the Farallon plate and rotated so that the East Pacific Rise in this immediate area became nearly perpendicular to the trench and perpendicular to the vector of motion of the Pacific plate relative to North America. Formation of the SMO batholith was followed between 24 and 20 Ma by a major increase in the rate of subduction of the Guadalupe plate, a fragment of the former Farallon plate, and by increasing mylonitization, extension, and uplift in the metamorphic core complexes that extend northwestward through southern Arizona from the northern end of the SMO batholith. The plate margin underwent another major change between 12.5 and 10 Ma when subduction again stopped, strike slip faulting became dominant along the coast, the Basin and Range Province opened, and numerous tectonostratigraphic terranes in southern California underwent large rotations. By 3 Ma a large, new terrane had been severed from North America immediately west of the SMO batholith as the Gulf of California opened. These observations can be explained by a model for the weakening and ultimate falling apart of the uppermost part of the subducted oceanic plate in the 20–30 m.y. after the end of rapid subduction. As the plate falls apart, not only is compressional stress relieved, but significant backslip along the old subduction zone is also possible, perhaps bringing blueschists rapidly upward from 20- to 30-km depths.

Journal of Geophysical Research↗

USGS2012

Fractal properties of tremor and gas piston events observed at Kilauea Volcano, Hawaii

Studies the fractal properties of shallow volcanic tremor and gas piston events associated with magma degassing at Kilauea Volcano, Hawaii, using data from two dense short-baseline arrays of seismographs deployed near the active crater of Puu Oo on the east rift of the volcano. The existence of a categorically stable attractor characterizing both types of activities strongly suggests that the excitation mechanism of tremor is similar to that of gas piston events, which in turn are correlated with visual observations at the volcanic vent. Low values derived for the dimension of the attractor in phase space points to significant self-organization in the process of generation of tremor and offers general constraints on the dimensionality of attractors derived from models of acoustic emission associated with magma flow, vesiculation, and degassing. -from Authors

Hawaii↗

USGS2012

Low intensity of the geomagnetic field in early Jurassic time

From a large collection of Jurassic continental tholeiites cropping out in Europe and Africa, we selected 90 samples for paleointensity determinations. The samples were carefully selected to avoid any secondary magnetizations, especially viscous magnetization. Use of the Thellier method reveals that magnetic property changes due to heating begin often at quite low temperatures but fortunately without modifying noticeably their natural remanent magnetization-thermoremanent magnetization ratio. Twenty-eight well-clustered paleointensity estimates were obtained from two European dikes that were emplaced during Early Jurassic time: the Kerforne dike at Brenterc'h in Brittany (northwestern France) and the Messejana dike on the Iberian Peninsula (Spain and Portugal). Virtual dipole moments calculated from both magmatic units are similar and only about one-third of present-day values. These new data lend support to the recently postulated low dipole moment of the Mesozoic geomagnetic field.

Journal of Geophysical Research↗

USGS2012

Seismic reflection/refraction mapping of faulting and regional dips in the eastern Alaska Range

We present the results of a Trans‐Alaska Crustal Transect (TACT) investigation of the upper 2–5 km of the eastern Alaska Range in the vicinity of the Denali fault based on seismic reflection/refraction data, laboratory measurements of rock velocities, and structural mapping. The Denali fault is a major dextral slip structure mappable for more than 2000 km separating the Wrangellia and adjacent terranes to the south from the Yukon‐Tanana and adjacent terranes to the north. Geologic mapping suggests over 400 km of dextral slip has occurred on the fault, yet within the upper 1.5 km of the crust along the TACT corridor, basement rocks juxtaposed along the Denali fault reveal no significant seismic velocity differences, although the fault zone itself is associated with a minor lowering of velocity. The lack of seismic velocity contrast adjacent to the fault is in agreement with laboratory measurements of elastic wave velocities of samples from terranes bordering the fault. Laboratory measurements of elastic wave velocities of the metasedimentary mica‐quartz schists comprising the Yukon‐Tanana basement are highly anisotropic because of preferred orientation of mica and predict significant variations in velocity accompanying variations in foliation dip. Although other interpretations are possible, the northward shallowing of foliation dip of basement rocks in the Yukon‐Tanana terrane combined with the strong anisotropy associated with these highly foliated rocks can explain an observed northward increase in seismic velocity within this terrane. Seismic reflections from basement rocks within the Yukon‐Tanana terrane may originate from variations in anisotropy with depth and/or changes in composition reflecting different proportions of sandstone and shale in the protolith.

Alaska↗

USGS2012

Lg and Rg waves on the California regional networks from the December 23, 1985 Nahanni earthquake

We investigate Lg and Rg propagation in California using the central and southern California regional networks. Approximately 550 stations constitute these two short-period networks providing a dense coverage of almost the entire state. The waveforms recorded from the December 23, 1985, Nahanni, Canada, earthquake are used to construct three profiles along the propagation path (almost N-S) and three perpendicular to the propagation path (almost E-W) to look at the nature of propagation of these two types of surface waves. Groups of records from stations in various geological and tectonic provinces in California are also examined in order to establish regional characteristics of the surface waves. We find that the propagation characteristics of Lg differ from those of Rg across California; Lg waves are apparently more sensitive to crustal heterogeneities. The most striking observations are the similarity of coda for both the Lg and the Rg waves within geologic provinces and the marked difference in coda between regions. These differences are seen in the amplitudes, coda duration, shape of the energy envelope, frequency content, and sharpness of the phase initiation. In general, a decrease in the Moho depth near the Pacific Coast is correlated with a decrease in the surface wave amplitude, especially at higher frequencies (0.15–0.2 Hz). Most interesting is the association of the San Andreas fault with abrupt changes in the wave train amplitudes. The surface waves are amplified in the vicinity of the fault zone and then decrease in amplitude after the zone is crossed. In the Coast Ranges, amplitudes are low and waveform coherence is poor. The Rg phase dominates the record in the Sierra Nevada, and both surface waves are amplified by the thick sedimentary sequence of the Great Valley.

Journal of Geophysical Research↗

USGS2012

An axial view of a metamorphic core complex: Crustal structure of the Whipple and Chemehuevi Mountains, southeastern California

A 135‐km‐long, NW‐SE trending, seismic refraction/wide‐angle reflection profile provides a unique along‐strike view of the crustal structure of a belt of metamorphic core complexes in southeastern California: the Whipple, Chemehuevi, and Sacramento mountains metamorphic core complexes. Interpretation of the seismic data was done by two‐dimensional forward modeling of travel times and amplitudes. The final model consists of (1) a thin (< 1.5 km) veneer of upper plate and fractured lower plate rocks (velocities of 1.5–5.3 km s −1 ) overlying a fairly homogeneous basement with velocities of 6.0 km s −1 ; (2) a localized, high‐velocity (6.4 km s −1 ) body, situated directly beneath the Whipple Mountains; (3) a 6.3–6.4 km s −1 middle crust that is thickest beneath the core complexes; (4) a 6.65±0.15 km s −1 lower crust; (5) crustal thickness of 27 km with a deeper crustal root (3 km) beneath the Whipple Mountains metamorphic core complex; and (6) a P n velocity of 8.0±0.10 km s −1 . The crustal structure that underlies the belt of metamorphic core complexes provides new insights into the processes that control extension in the deep crust. Upper crustal velocities are higher beneath the Whipple Mountains (where velocities increase to 6.4 km s −1 at ∼5 km depth) than beneath the Chemehuevi and Sacramento mountains. In addition, midcrustal discontinuities rise 2–5 km beneath the Whipple complex compared to the other complexes. These observations support greater uplift and a slightly deeper midcrustal origin for the rocks now exposed in the core of the Whipple Mountains compared to rocks in the Chemehuevi and Sacramento mountains. Despite the enhanced uplift and extension in the Whipple Mountains, the crust is thicker here (30 km) than anywhere else along the Colorado River extensional corridor. This may be in part a relic of compressional and magmatic thickening during the Mesozoic. However, we suggest that inflation of the crust during Tertiary extension was the dominant mechanism. Both mantle‐derived magmatism and lateral ductile inflow in the crust are proposed.

Arizona, California↗

USGS2012

Anatomy of a metamorphic core complex: seismic refraction/wide-angle reflection profiling in southeastern California and western Arizona

The metamorphic core complex belt in southeastern California and western Arizona is a NW-SE trending zone of unusually large Tertiary extension and uplift. Midcrustal rocks exposed in this belt raise questions about the crustal thickness, crustal structure, and the tectonic evolution of the region. Three seismic refraction/wide-angle reflection profiles, acquired and analyzed as a part of the U.S. Geological Survey's Pacific to Arizona Crustal Experiment, were collected to address these issues. The results presented here, which focus on the Whipple and Buckskin-Rawhide mountains, yield a consistent three-dimensional image of this part of the metamorphic core complex belt. The seismic refraction/wide-angle reflection data are of excellent quality and are characterized by six principal phases that can be observed on all three profiles. These phases include refractions from the near-surface and crystalline basement, reflections from boundaries in the middle and lower crust, and reflections and refractions from the upper mantle. The final model consists of a thin veneer (<2 km) of upper plate and fractured lower plate rocks (1.5–5.5 km s−1 ) overlying a fairly homogeneous basement (∼6.0 km s −1 ) and a localized high-velocity (6.4 km s −1 ) body situated beneath the western Whipple Mountains. A prominent midcrustal reflection is identified beneath the Whipple and Buckskin-Rawhide mountains between 10 and 20 km depth. This reflector has an arch-like shape and is centered beneath, or just west of, the metamorphic core complex belt. This event is underlain by a weaker, approximately subhorizontal reflection at 24 km depth. Together, these two discontinuities define a lens-shaped midcrustal layer with a velocity of 6.35–6.5 km s −1 . The apex of this midcrustal layer corresponds roughly to a region of major tectonic denudation and uplift (∼10 km) defined by surface geologic mapping and petrologic barometry studies. The layer thins to the northeast and is absent in the Transition Zone. The 6.35–6.5 km s −1 velocities are compatible with a diorite composition or a mixture of mafic and silicic rocks. This midcrustal layer is underlain by a higher-velocity lower crustal layer that is modeled as only 3–6 km thick beneath the metamorphic core complex belt and regions to the southwest. To the northeast, however, this layer thickens to 8–10 km as the midcrustal layer pinches out above it. The velocity of the lower crust is constrained by traveltime modeling and is 6.6±0.15 kms −1 beneath the western Transition Zone and the metamorphic core complex belt; higher velocities may be present farther to the southwest where the layer is thin. The velocity of the lower crust is too low to accommodate significant amount of mafic underplating at the base of the crust. Instead, we interpret the velocities to indicate that the lower crust is passively thinned beneath these regions without significant addition of mafic mantle-derived intrusions. The crust-mantle boundary does not dome up beneath the core complexes but remains approximately subhorizontal at a depth of 26–28 km or, in the case of the Whipple Mountains, actually deepens; a 3-km crustal root is modeled. This lack of upward doming of the Moho, together with the vertical alignment of the metamorphic core complex belt over what are believed to be extension-related structures in the middle and lower crust, suggest that there is no lateral offset of upper crustal deformation from deeper zones of extension, as one would expect if extension occurred along crust-penetrating shear zones (Wernicke, 1981; Wernicke et al., 1985). Instead, domed and inflated middle crust and thinned lower crust directly underlie the region of greatest thinning of the upper crust.

Journal of Geophysical Research↗

USGS2012

The high-pressure electronic structure of magnesiowustite (Mg, Fe)O: applications to the physics and chemistry of the lower mantle

The electronic structure of magnesiowustite is investigated using self-consistent field X α scattered wave (SCF- X α-SW) molecular orbital calculations on (FeO 6 ) 10− and (FeMg 12 O 14 ) 2− clusters. Calculated one-electron transition energies are used to interpret the optical spectrum of (Mg, Fe)O. The results are applied to the electrical and thermal conductivity of the lower mantle. The spin pairing of Fe 2+ and the effect of pressure on bonding in magnesiowustite, with some inferences regarding the incorporation of oxygen in the outer core, is also addressed. The approach used here appears to give a reliable description of the energy and pressure dependence of the spin-allowed 5 T 2 g → 5 E g ligand field transition and the spin-pairing transition of Fe 2+ in (Mg, Fe)O. However, the oxygen to metal charge transfer transitions in (Mg, Fe)O are not as reliably determined insofar as the p - d band gap varies with cluster size and the energies of the charge transfer states cannot be found without including configurational interaction. Nevertheless, it is argued that the charge transfer transitions that are intrinsic to (Fe, Mg)O are of a sufficiently high energy to be irrelevant to the electrical and thermal conductivity of the lower mantle. This is especially true if Fe 2+ adopts the low-spin configuration. The geophysically significant properties of (Fe, Mg)O probably result from defect Fe 3+ .

Journal of Geophysical Research↗

USGS2012

Strain accumulation in western Washington

The Juan de Fuca plate is subducted beneath the North American plate off the coast of Washington at a rate of about 40 mm/yr N68°E. The average principal strain rates (extension reckoned positive) measured in northwestern Washington are as follows: Olympic peninsula 25 km south of Port Angeles from 1982 through 1990, and and near Seattle from 1972 through 1985, and . Both strain measurements are consistent with uniaxial contraction in the direction of plate convergence. Uplift rates inferred from tide gage recordings are about 4 mm/yr on the Pacific coast and near 0 mm/yr farther inland near Seattle. These deformation rates are consistent with a model of the Cascadia subduction zone in which the plate interface beneath the continental slope and outer continental shelf is locked but free to slip farther landward. The limited downdip extent of the locked segment of the plate interface is consistent with a shallow depth (∼20 km) of the isotherm (∼450°C) that defines the brittle‐ductile transition. Small thrust events diagnostic of seismic subduction should then occur only offshore and at shallow depths. The principal strain rates measured from 1972 through 1983 in the back arc region near Richland, Washington, are and .

Washington↗

USGS2012

Strain accumulation along the Denali Fault at the Nenana River and Delta River Crossings, Alaska

Surveys of trilateration networks across the Denali fault at the Nenana River in 1982, 1984, and 1988 and at the Delta River in 1975, 1979, 1982, and 1984 indicate a minor (0.10±0.04 μstrain/yr) northeastward uniaxial extension. The component of right‐lateral shear‐strain accumulation across the fault is not significant at the two‐standard‐deviation level. At the Delta River network the strain accumulation rate decreases rapidly with distance from the fault, but evidence for a similar decrease with distance from the fault is lacking at the Nenana River network. The strain accumulation rates inferred from trilateration are consistent with the very long baseline interferometry (VLBI) measurement reported by Ma et al. (1990) and support their contention that significant right‐lateral shear is not accumulating along the Denali fault at the present time. Savage et al. (1981) had earlier concluded erroneously that preliminary geodetic measurements at the Delta River network demonstrated right‐lateral shear strain accumulation. The absence of significant right‐lateral deformation across the Denali fault in the 1975–1988 interval is in marked contrast with the abundant geomorphic evidence for Holocene right‐lateral secular slip at the rate of 10–20 mm/yr on the Denali fault in this sector.

Alaska↗

USGS2012

Multicycle slip distribution along a laboratory fault

Slip distribution along a laboratory fault, which consists of eight spring-connected blocks that are elastically driven to slide on a frictional surface, has been examined for a “long” sequence of slip events to test the applicability of some conceptual models proposed recently in the literature. The distributions of large slip events are found to be quite variable and do not fit the uniform slip or characteristic earthquake models. The rupture initiation points are usually not near the corresponding maximum slip points, in contrast to observations by Thatcher (1990) and by Fukao and Kikuchi (1987) that earthquake hypocenters are commonly near corresponding regions of maximum slip in the fault planes. This contrast may suggest that either the present observations or theirs are not representative or the teleseismically determined hypocenters may not always be true rupture initiation points as usually assumed. Large slip events are also found to be a stress-roughening process. They are triggered by some small events after the stresses have been adjusted by some earlier small-to-moderate events to be near the critical levels at most locations along the fault. This suggests that earthquake prediction monitoring efforts should not be limited to a small region near an asperity but should be spread out to cover the entire fault segment in a seismic gap in order to detect the condition of simultaneous strain buildup.

Journal of Geophysical Research↗

USGS2012

Observations constraining near-source ground motion estimated from locally recorded seismograms

To estimate the seismic hazard to underground facilities or operations in the environs of a mining-induced tremor or a natural earthquake, it is useful to be able to relate locally recorded seismic waveforms to peak ground velocity and slip at the causative fault. For this purpose, far-field S wave pulses are analyzed to define the faulting slip D and near-fault peak ground velocity D /2 that give rise to the most significant ground motion. This most intense region of faulting, an assumed circular asperity, has radius r within a broader source zone of radius r 0 , which is traditionally calculated from the corner frequency of the S wave spectrum. In developing relationships between peak far-field velocity v and peak acceleration a, and the source processes of the asperity, D and D , as well as its radius r , the key model assumption is that r = k β/ω, where ω is the angular frequency of the sinusoidal velocity pulse of maximum amplitude, β is the sheaf wave speed, and k is a constant. Observations in deep-level gold mines of fault slip and slip velocity as well as laboratory observations of slip rate as a function of stress drop for stick-slip failure support a choice of about k = 2.34, the value commonly used for estimating r 0 using the Brune model. In particular, observations of fault slip up to 410 mm for mining-induced tremors in the moment magnitude range 4–5 are consistent with D = 8.1 R v/β, where R is hypocentral distance. Moreover, estimates based on underground damage of near-fault ground velocities ranging up to 3.5 m/s are in accord with D /2 = 1.28(β/μ) ρ R a, where μ is the modulus of rigidity and ρ is the density. Alternatively, the average slip velocity 〈 D 〉 can be expressed in terms of the stress drop Δσ a of the asperity as 〈 D 〉 = 0.51 β Δσ a /μ, and the agreement of this relationship with measurements made during stick-slip failure in the laboratory is good. To the extent that seismic slip exterior to the asperity is a consequence of preevent suppression of slip due to the asperity, the broader-scale( r 0 ) slip can be related to that of the asperity. Just as the asperity radius r can be estimated from r = 2.34 βv/a, an alternative estimate for r 0 is given by r 0 = ρ R a M 0 /[75.8ρμ( R v) 2 ], the results of which are generally in good agreement with estimates based on the spectral corner frequency method.

Journal of Geophysical Research↗

USGS2012

Seismicity and shear strain in the southern Great Basin of Nevada and California

This study examines the relationship between the distribution of small earthquakes (M L ≤4.3) and mechanisms of strain accumulation and relaxation in an area with long repeat times between large events, the Southern Great Basin Seismic Network (SGBSN) region. The Great Basin is a unique continental extensional province characterized by normal and strike-slip faulting, high heat flow, crust of thin to normal thickness, and high elevations. The SGBSN is operated to provide data to address suitability issues pertaining to Yucca Mountain, Nevada which is being evaluated as a potential site for a national mined geologic nuclear waste repository. Suitability issues include estimation of the probability of occurrence of future damaging earthquakes, the characterization of the mechanisms that drive hydrologic flow, and the identification of fractures (faults) that might act as flow conduits or barriers. This study attempts to explain the distribution of small earthquakes in terms of spatial variations in the shear strain field; where strain concentrates there should be a greater number of small earthquakes. Strain field models are constructed under the assumption that long term fault behavior perturbs an otherwise uniform strain field. These strain field models are then interpreted with regard to the regional tectonics and site suitability issues. Modeling results provide one possible explanation of why earthquake clusters cover regions much larger than the surface projections of any of mapped major faults; clusters in a wide band along and extending northeast of the northern half of the Furnace Creek fault may correspond to elevated shear strains along the fault and a broad cluster in the Pahranagat Shear Zone may be associated with shear strain arising from a distribution of smaller localized faults. The relatively large number of small earthquakes in the southern and eastern portions of the Nevada Test Site is consistent with the strain field models. A minimum in shear strain at Yucca Mountain is predicted by all models consistent with an almost total lack of earthquakes observed there. The region to the west of the Death Valley/Furnace Creek fault system, the portion of the study area with the most active deformation but few small earthquakes, is an area of low shear strain. A possible reason for this is that the fault configuration in the area is optimal for accommodating regional deformation via large earthquakes or creep. While there is also a relative lack of earthquakes at Yucca Mountain, this may be indicative of a lack of accumulating strain energy and thus, a lower potential for a large earthquake.

Journal of Geophysical Research↗

USGS2012

Historic creep rate and potential for seismic slip along the Hayward Fault, California

The Hayward fault is considered the most likely source of one or more major earthquakes in the San Francisco Bay area in the next few decades. Historically, at least one, and probably two, major earthquakes (about M 6.8) occurred along the Hayward fault, one in 1836 and another in 1868. Little is known about the 1836 event, but the 1868 earthquake was accompanied by a surface rupture that extended as much as 41 km along the southern part of the fault. Although the amount of surface slip in 1868 is uncertain, right slip (including afterslip) reached at least several centimeters, and possibly several decimeters in places. This paper documents the spatial variation of creep rate along the Hayward fault since the 1868 earthquake. Creep (aseismic fault slip) occurs over at least 66 km and may extend over the fault's entire 82-km length, of which about 13 km lies underwater. Creep rate seems nearly constant over decades, but short-term variations occur. We derive creep rate mainly from our own systematic surveying of offset cultural features (curbs, fences, and buildings). On each feature we solve directly for accumulated creep by using multiple linear regression. Creep rate mostly falls in the range of 3.5–6.5 mm/yr; but systematic variation occurs along strike. Fault segments with distinctly higher and lower rates generally correspond to parts of the fault most salient from the overall average alinement of the fault. Most distinctive is a 4-km-long section near the south end of the fault that creeps at about 9 mm/yr. Such a high rate has occurred there at least since the 1920s and probably since the 1868 earthquake, as indicated by an offset railroad track built in 1869. We suggest that this 9 mm/yr slip rate may approach the long-term or deep slip rate that controls average recurrence interval between major earthquakes. If so, assuming an elastic rebound model, the potential for slip in large earthquakes below the surficial creeping zone is now &sim;1.1 m in the southern (1868) segment of the fault and ≥ 1.4 m in the northern (1836?) segment. Subtracting surface creep rates from a long-term slip rate of 9 mm/yr gives present potential for surface slip in large earthquakes of up to 0.8 m, with an average of 0.6 m in the northern segment and 0.4 m in the southern segment. We present a simple hypothesis for rupture potential that is compatible with historic creep rate, microseismicity distribution, and geodetic data. If seismic rupture occurs on segments 41 km long by 10 km deep (7 km fully locked, 3 km creeping), today's potential for seismic moment release is 1.4 × 10 19 and 1.1 × 10 19 N m for both 1836? and 1868 segments, respectively, and 2.5 × 10 19 N m for both segments jointly. Converting moment to magnitude gives M L 6.8 in the northern segment, M L 6.7 in the southern segment, and M L 7.0 for simultaneous rupture of both.

California↗