Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical 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 1,171 records · Page 65Linked to original sources
USGS2012

The detection of coastal-trapped waves

We outline a simple method for estimating the cross-spectral matrix of coastal-trapped wave amplitudes,A, from a set of oceanographic observations. Specifically, we propose that A may be estimated by (M'M) -1 M'UM(M'M) -1 where a prime denote conjugate transpose, U is the sample cross-spectral matrix of observations and M is a matrix which has the spatial form of the waves for columns. In general, M will be complex and frequency-dependent. We discuss the bias of this estimator and show how to estimate the variance of the power and cross spectra of wave amplitudes. We also outline an ad hoc scheme for assessing the predictive skill of the coastal trapped wave representation and finally give some advice on how to interpret A. Although the method is presented in the context of shelf circulation and coastal trapped waves, it may be applied to any linear system where the spatial forms of the waves are known and the cross-spectral matrix of their amplitudes is required.

Journal of Geophysical Research C: Oceans↗
USGS2012

Atmospheric organochlorine pollutants and air-sea exchange of hexachlorocyclohexane in the Bering and Chukchi Seas

Organochlorine pesticides have been found in Arctic fish, marine mammals, birds, and plankton for some time. The lack of local sources and remoteness of the region imply long-range transport and deposition of contaminants into the Arctic from sources to the south. While on the third Soviet-American Joint Ecological Expedition to the Bering and Chukchi Seas (August 1988), high-volume air samples were taken and analyzed for organochlorine pesticides. Hexachlorocyclohexane (HCH), hexachlorobenzene, polychlorinated camphenes, and chlordane (listed in order of abundance, highest to lowest) were quantified. The air-sea gas exchange of HCH was estimated at 18 stations during the cruise. Average alpha-HCH concentrations in concurrent atmosphere and surface water samples were 250 pg m-3 and 2.4 ng L-1, respectively, and average gamma-HCH concentrations were 68 pg m-3 in the atmosphere and 0.6 ng L-1 in surface water. Calculations based on experimentally derived Henry's law constants showed that the surface water was undersaturated with respect to the atmosphere at most stations (alpha-HCH, average 79% saturation; gamma-HCH, average 28% saturation). The flux for alpha-HCH ranged from -47 ng m-2 day-1 (sea to air) to 122 ng m-2 d-1 (air to sea) and averaged 25 ng m-2 d-1 air to sea. All fluxes of gamma-HCH were from air to sea, ranged from 17 to 54 ng m-2 d-1, and averaged 31 ng m-2 d-1.

Journal of Geophysical Research↗
USGS2012

Oligocene basaltic volcanism of the northern Rio Grande Rift: San Luis Hills, Colorado

The inception of the Rio Grande rift in northern New Mexico and southern Colorado was accompanied by voluminous mafic volcanism preserved in part as erosional remnants on an intrarift horst within the current axial rift graben of the San Luis Valley. Oligocene (∼26 Ma) volcanic rocks of the Hinsdale Formation at San Luis Hills range from 49 to 57 wt % SiO 2 and include nepheline and hypersthene normative lavas. A mildly alkalic series consisting of trachybasalt, basaltic trachyandesite, and trachyandesite is volumetrically dominant, olivine tholeiites are subordinate, and xenocrystic trachyandesites containing abundant quartz and plagioclase xenocrysts occur only locally. Relative to the San Luis Hills olivine tholeiites which have La/Sm n ∼ 2, the more alkaline series are enriched in light rare earth elements (LREE) and have La/Sm ratios that increase in the trachybasalt-basaltic trachyandesite suite (La/Sm n ∼ 3) to xenocrystic trachyandesites that are the most LREE enriched (La/Sm n ∼ 4). Chondrite-normalized, trace element patterns for the lavas in the San Luis Hills are similar in shape within the mildly alkaline to transitional series; they have characteristic Nb and Ta depletions and high K and Th relative to Ta, Nb, and LREE. Major and trace element constraints support a petrogenetic model of fractionation plus lower crustal assimilation for petrologic suites within the San Luis Hills rocks, although the model cannot relate lavas for the entire series to a common parent. Most mafic lavas of the San Luis Hills were evolved (Mg # <60) and contaminated by LREE-enriched silicic partial melts of granulitic lower crust depleted in Rb, Th, and U. Pb isotopes are the most sensitive indicators of crustal contamination, whereas shifts in Nd and Sr isotope ratios are associated with large amounts of assimilation. However, relatively noncontaminated lavas can be identified and indicate at least two mantle source regions were involved.

Journal of Geophysical Research↗
USGS2012

Intrusion of horizontal dikes: tectonic significance of Middle Proterozoic diabase sheets widespread in the upper crust of the southwestern United States

Initially horizontal sheet intrusions of Middle Proterozoic diabase are abundant in a region 650 by 300 km across in Arizona and California. The diabase forms discordant sheets in basement granite and gneiss and sills in overlying shelf sedimentary sequences. Massive granite is the most common basement host for the sheets, probably because it fractured more easily than foliated hosts during sheet emplacement. Steep feeder dikes are rare compared to the sheets. The diabase in many places is exposed in fault blocks that were tilted during Tertiary tectonic extension. Structure sections restored from the map patterns of upended blocks show that the sheets were intruded at levels throughout the upper crust, to depths of at least 13 km. Sheet intrusion implies a vertical orientation of the least compressive stress, so I conclude that the crust was under tectonic compression or in an isotropic state of stress at the time of diabase intrusion about 1.1 Ga. Magma overpressures, water encountered by rising magma, vertical changes in the crustal stress regime, and flotation of low-density granite all may be important factors for sheet intrusion. The stress conditions suggested by the presence of the sheets argue against an extensional tectonic regime earlier proposed for the diabase event. The diabase province contrasts structurally with similar-age provinces of basaltic magmatism elsewhere in North America that show evidence of tectonic extension, such as the midcontinent rift. The consistent orientations of the sheets commonly allow them to be used as structural markers for postdiabase deformation of basement blocks. Sheet intrusions in the geologic record may be seriously underreported. Their recognition is important for the interpretation of seismic reflection profiles of continental crust. Sheets in Arizona may be responsible for the Bagdad reflection sequence, which extends to depths of at least 15 km.

Journal of Geophysical Research↗
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

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

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

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

Crustal subsidence and extension and Medicine Lake volcano, northern California

The pattern of historical ground deformation, seismicity, and crustal structure near Medicine Lake volcano illustrates a close relation between magmatism and tectonism near the margin of the Cascade volcanic chain and the Basin and Range tectonic province. Between leveling surveys in 1954 and 1989 the summit of Medicine Lake volcano subsided 389±43 mm with respect to a reference bench mark 40 km to the southwest (average rate = 11.1±1.2 mm/yr). A smaller survey across the summit caldera in 1988 suggests that the subsidence rate was 15–28 mm/yr during 1988–1989. Swarms of shallow earthquakes ( M ≤ 4.6) occurred in the region during August 1978, January–February 1981, and September 1988. Except for the 1988 swarm, which occurred beneath Medicine Lake caldera, most historical earthquakes were located at least 25 km from the summit. The spatial relation between subsidence and seismicity indicates (1) radially symmetric downwarping of the volcano's summit and flanks centered near the caldera and (2) downfaulting of the entire edifice along regional faults located 25–30 km from the summit. We propose that contemporary subsidence, seismicity, and faulting are caused by (1) loading of the crust by more than 600 km 3 of erupted products plus a large volume of mafic intrusives; (2) east‐west extension in the western Basin and Range province; and, to a lesser extent, (3) crystallization or withdrawal of magma beneath the volcano. Thermal weakening of the subvolcanic crust by mafic intrusions facilitates subsidence and influences the distribution of earthquakes. Subsidence occurs mainly by aseismic creep within 25 km of the summit, where the crust has been heated and weakened by intrusions, and by normal faulting during episodic earthquake swarms in surrounding, cooler terrain.

California↗