USGS Science⌕ Search

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.

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

At least 649 records · Page 36Linked to original sources
USGS2005

Debris-bed friction of hard-bedded glaciers

[1] Field measurements of debris-bed friction on a smooth rock tablet at the bed of Engabreen, a hard-bedded, temperate glacier in northern Norway, indicated that basal ice containing 10% debris by volume exerted local shear traction of up to 500 kPa. The corresponding bulk friction coefficient between the dirty basal ice and the tablet was between 0.05 and 0.08. A model of friction in which nonrotating spherical rock particles are held in frictional contact with the bed by bed-normal ice flow can account for these measurements if the power law exponent for ice flowing past large clasts is 1. A small exponent (n < 2) is likely because stresses in ice are small and flow is transient. Numerical calculations of the bed-normal drag force on a sphere in contact with a flat bed using n = 1 show that this force can reach values several hundred times that on a sphere isolated from the bed, thus drastically increasing frictional resistance. Various estimates of basal friction are obtained from this model. For example, the shear traction at the bed of a glacier sliding at 20 m a-1 with a geothermally induced melt rate of 0.006 m a-1 and an effective pressure of 300 kPa can exceed 100 kPa. Debris-bed friction can therefore be a major component of sliding resistance, contradicting the common assumption that debris-bed friction is negligible. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2005

Generation and verification of theoretical rating curves in the Whitewater River basin, Kansas

[1] A new method for generating stage-discharge relations (rating curves) for geomorphically stable channels is presented and applied to two streams in the Whitewater River basin, Kansas. The approach converts measurements of stage into discharge using a fluid mechanically based model. The model does not use empirical roughness coefficients, such as Manning coefficients, but rather determines channel roughness from field measurements of the (1) channel geometry, (2) the physical roughness of the bed, banks, and floodplain, and (3) the vegetation density on the banks and floodplain. These measurements are used to calculate explicitly the drag on the small-scale topographic features on the boundary, the drag on the stems and branches of woody vegetation, and the friction on the bed, banks, and floodplain. The theoretical rating curves produced by the model for two study reaches, which are near U.S. Geological Survey (USGS) streamflow-gauging stations, are in good agreement with direct measurements of discharge made by the USGS. Our method has the potential of providing accurate estimates of stream flows less expensively than conventional gauging methods. In addition, the method can be used to obtain more accurate discharge estimates than conventional indirect methods for determining discharge, which are based on estimates of Manning's coefficient.

Journal of Geophysical Research F: Earth Surface↗
USGS2006

Effects of tributary debris on the longitudinal profile of the Colorado River in Grand Canyon

The Colorado River in Grand Canyon has long been known as a "rapids-and-pools" river, with the rapids owing their existence primarily to tributary debris flows. The debris flows deposit subaerial debris fans that constrict the channel laterally and, when they enter the river, raise the bed elevation. The rapids are short-wavelength (???0.1 to ???1 km), small-amplitude (??????5 m) convexities in the river's longitudinal profile, arising from the shallow gradient in the upstream pool and the steep gradient through the rapid itself. Analysis of the entire longitudinal profile through Grand Canyon reveals two long-wavelength (???100 km), large-amplitude (15-30 m) river profile convexities: the eastern canyon convexity between river mile (RM) 30 and RM 80 and the western canyon convexity between RM 150 and RM 250. Convexities of intermediate scale are also identified in the longitudinal profile. These longer-wavelength, larger-amplitude convexities have strong spatial correlations with high rates of debris flow occurrence, high densities of Holocene debris fans, the largest debris fans along the river, and alluvial thicknesses of 10 m or more. River profile convexities are unstable and require an active and powerful geologic process to maintain them, in this case the abundant, frequent, and voluminous Holocene debris flow activity in Grand Canyon. At all wavelengths the most likely cause for these river profile convexities is Holocene aggradation of the riverbed beneath them, driven by the coarse particles of tributary debris flows. Large enough debris flows will slow river flow for kilometers upstream, causing it to drop much of its suspended load. Integrated over time and all of the tributary point source contributions, this process will build short-wavelength convexities into long-wavelength convexities. For most if not all of the Holocene the Colorado River has been dissipating most of its energy in the rapids and expending the remainder in transporting fine sediment through Grand Canyon, with little or no regional incision of bedrock.

Journal of Geophysical Research F: Earth Surface↗
USGS2005

Regulation of landslide motion by dilatancy and pore pressure feedback

A new mathematical model clarifies how diverse styles and rates of landslide motion can result from regulation of Coulomb friction by dilation or contraction of water‐saturated basal shear zones. Normalization of the model equations shows that feedback due to coupling between landslide motion, shear zone volume change, and pore pressure change depends on a single dimensionless parameter α, which, in turn, depends on the dilatancy angle ψ and the intrinsic timescales for pore pressure generation and dissipation. If shear zone soil contracts during slope failure, then α < 0, and positive pore pressure feedback and runaway acceleration are inevitable. If the shear zone dilates, then α > 0, and negative feedback permits slow, steady landslide motion to occur while positive pore pressure is supplied by rain infiltration. Steady state slip velocities v 0 obey v 0 = −( K /ψ) p * e , where K is the hydraulic conductivity and p * e is the normalized (dimensionless) negative pore pressure generated by dilation. If rain infiltration and attendant pore pressure growth continue unabated, however, their influence ultimately overwhelms the stabilizing influence of negative p * e . Then, unbounded landslide acceleration occurs, accentuated by an instability that develops if ψ diminishes as landslide motion proceeds. Nonetheless, numerical solutions of the model equations show that slow, nearly steady motion of a clay‐rich landslide may persist for many months as a result of negative pore pressure feedback that regulates basal Coulomb friction. Similarly stabilized motion is less likely to occur in sand‐rich landslides that are characterized by weaker negative feedback.

Journal of Geophysical Research F: Earth Surface↗
USGS2006

The 1923 Kanto earthquake reevaluated using a newly augmented geodetic data set

This study revisits the mechanism of the 1923 Ms = 7.9 Kanto earthquake in Japan. We derive a new source model and use it to assess quantitative and qualitative aspects of the accommodation of plate motion in the Kanto region. We use a new geodetic data set that consists of displacements from leveling and angle changes from triangulation measurements obtained in surveys between 1883 and 1927. Two unique aspects of our analysis are the inclusion of a large number of second-order triangulation measurements and the application of a correction to remove interseismic deformation. The geometry of the fault planes is adopted from a recent seismic reflection study of the Kanto region. We evaluate the minimum complexity necessary in the model to fit the data optimally. Our final uniform-slip elastic dislocation model consists of two adjacent ???20?? dipping low-angle planes accommodating reverse dextral slip of 6.0 in on the larger, eastern plane and 9.5 m on the smaller, western plane with azimuths of 163?? and 121??, respectively. The earthquake was located in the Sagami trough, where the Philippine Sea plate subducts under Honshu. Compared to the highly oblique angle of plate convergence, the coseismic slip on the large fault plane has a more orthogonal orientation to the strike of the plate boundary, suggesting that slip partitioning plays a role in accommodation of plate motion. What other structure is involved in the partitioning is unclear. Uplift records of marine coastal terraces in Sagami Bay document 7500 years of earthquake activity and predict average recurrence intervals of 400 years for events with vertical displacement profiles similar to those of the 1923 earthquake. This means that the average slip deficit per recurrence interval is ???50% of the relative plate convergence. These findings of plate motion partitioning and slip deficit lead us to suggest that instead of a simple recurrence model with characteristic earthquakes, additional mechanisms are necessary to describe the accommodation of deformation in the Kanto region. So far, obvious candidates for these alternative mechanisms have not been discovered. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Coseismic slip distribution of the 1923 Kanto earthquake, Japan

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

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Earthquake fracture energy inferred from kinematic rupture models on extended faults

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

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Spatiotemporal evolution of a transient slip event on the San Andreas fault near Parkfield, California

In 1993 several baselines of the two-color electronic distance meter (EDM) network at Parkfield, California, deviated from their long-term rates, coincident with anomalous observations from nearby strain meters and a creep meter, as well as an increase in microseismicity. Between October 1992 and December 1994, three M ??? 4.5 earthquakes occurred beneath Middle Mountain, near the hypocenter of the 1934 and 1966 Parkfield M6 events. We analyzed the two-color EDM data using a Kalman-filtering based technique to image the spatiotemporal evolution of slip on the fault at Parkfield between the mid-1980s and 2003. This method accounts for localized random walk motion of the geodetic monuments and a prominent seasonal signal that affects many baselines. We find that a slip rate increase occurred between January 1993 and July 1996 on the upper 8 km of the fault near Middle Mountain. The peak estimated slip rate during this time was 49 mm/yr, which exceeds the long-term geologic rate of ???35 mm/yr. The slip rate evolution appears episodic, with an initial modest increase after the M4.3 earthquake and a much larger jump following the shallower M4.7 event in December 1994. This temporal correlation between inferred slip and seismicity suggests that the moderate earthquakes triggered the aseismic fault slip. The EDM data cannot resolve whether transient slip propagated across the nucleation zone of the 1934 and 1966 M6 Parkfield earthquakes. However, transient slip and its associated stress release in the hypocentral area of previous Parkfield events is consistent with the nucleation of the 2004 M6 Parkfield earthquake elsewhere on the fault. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Transient rheology of the upper mantle beneath central Alaska inferred from the crustal velocity field following the 2002 Denali earthquake

The M7.9 2002 Denali earthquake, Alaska, is one of the largest strike-slip earthquakes ever recorded. The postseismic GPS velocity field around the 300-km-long rupture is characterized by very rapid horizontal velocity up to ???300 mm/yr for the first 0.1 years and slower but still elevated horizontal velocity up to ???100 mm/yr for the succeeding 1.5 years. I find that the spatial and temporal pattern of the displacement field may be explained by a transient mantle rheology. Representing the regional upper mantle as a Burghers body, I infer steady state and transient viscosities of ??1 = 2.8 ?? 1018 Pa s and ??2 = 1.0 ?? 1017 Pa s, respectively, corresponding to material relaxation times of 1.3 and 0.05 years. The lower crustal viscosity is poorly constrained by the considered horizontal velocity field, and the quoted mantle viscosities assume a steady state lower crust viscosity that is 7??1. Systematic bias in predicted versus observed velocity vectors with respect to a fixed North America during the first 3-6 months following the earthquake is reduced when all velocity vectors are referred to a fixed site. This suggests that the post-Denali GPS time series for the first 1.63 years are shaped by a combination of a common mode noise source during the first 3-6 months plus viscoelastic relaxation controlled by a transient mantle rheology.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Northwest Basin and Range tectonic deformation observed with the Global Positioning System, 1999-2003

We use geodetic velocities obtained with the Global Positioning System (GPS) to quantify tectonic deformation of the northwest Basin and Range province of the western United States. The results are based on GPS data collected in 1999 and 2003 across five new quasi-linear networks in northern Nevada, northeast California, and southeast Oregon. The velocities show ???3 mm/yr westward movement of northern Nevada with respect to stable North America. West of longitude 119??W the velocities increase and turn northwest, parallel to Sierra Nevada/Great Valley microplate motion, and similar to velocities previously obtained to the south. The observations are explained by a kinematic model with three domains that rotate around Euler poles in eastern Oregon and western Idaho. Northeast California experiences internal dextral shear deformation (11.2 ?? 3.6 nstrain/yr) subparallel to Pacific/North America motion. Relative motions of the domains imply 2-5 mm/yr approximately east-west extension in northwest Nevada and 1-4 mm/yr approximately north-south contraction near the California/Oregon border. The northward decreasing approximately east-west extension in northwest Nevada is consistent with the northern termination of Basin and Range deformation, faulting and characteristic topography. No significant extension is detected in the Oregon Basin and Range. The Oregon Cascade arc moves north at ???3.5 mm/yr and is possibly influenced by the approximately eastward motion of the Juan de Fuca plate. These results disagree with secular northwest trenchward motion of the Oregon forearc inferred from paleomagnetic rotations. South of latitude 43??, however, trenchward motion exists and is consistent with block rotations, approximately east-west Basin and Range extension, and northwest Sierra Nevada translation. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Postseismic relaxation and transient creep

Postseismic deformation has been observed in the epicentral area following the 1992 Landers (M = 7.3), 1999 Chi-Chi (M = 7.6), 1999 Hector Mine (M = 7.1), 2002 Denali (M = 7.9), 2003 San Simeon (M = 6.5), and 2004 Parkfield (M = 6.0) earthquakes. The observations consist of repeated GPS measurements of the position of one monument relative to another (separation ???100 km). The early observations (t < 0.1 year) are well fit by the function a' + c'log(t), where t is the time after the earthquake and a' and c' are constants chosen to fit the data. Because a log(t) time dependence is characteristic of transient (primary) creep, the early postseismic response may be governed by transient creep as Benioff proposed in 1951. That inference is provisional as the stress conditions prevailing in postseismic relaxation are not identical to the constant stress condition in creep experiments. The observed logarithmic time dependence includes no characteristic time that might aid in identifying the micromechanical cause.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Postseismic deformation following the June 2000 earthquake sequence in the south Iceland seismic zone

We observe postseismic deformation on two spatiotemporal scales following Mw = 6.5 earthquakes in the south Iceland seismic zone on 17 and 21 June 2000. We see a rapidly decaying deformation transient lasting no more than 2 months and extending about 5 km away from the two main shock ruptures. This local, month-scale transient is captured by several radar interferograms and is also observed at a few campaign GPS sites located near the faults. A slower transient with a characteristic timescale of about a year is detected only by GPS measurements. The month-scale deformation pattern has been explained by poroelastic rebound due to postearthquake pore pressure changes. In contrast, the year-scale deformation can be explained by either afterslip at 8-14 km depth or viscoelastic relaxation of the lower crust and upper mantle in response to the coseismic stress changes. The optimal viscoelastic models have lower crustal viscosities of 0.5-1 ?? 1019 Pa s and upper mantle viscosity of ???3 ?? 1018 Pa s. Because of the limitations of our GPS campaign data, we consider both afterslip and viscoelastic relaxation as plausible mechanisms explaining the deformation field. Both types of postseismic deformation models suggest that the areas of large coseismic stress increase east of the 17 June and west of the 21 June ruptures continue to be loaded by the postseismic deformation. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Gas slug ascent through changes in conduit diameter: Laboratory insights into a volcano-seismic source process in low-viscosity magmas

Seismic signals generated during the flow and degassing of low-viscosity magmas include long-period (LP) and very-long-period (VLP) events, whose sources are often attributed to dynamic fluid processes within the conduit. We present the results of laboratory experiments designed to investigate whether the passage of a gas slug through regions of changing conduit diameter could act as a suitable source mechanism. A vertical, liquid-filled glass tube featuring a concentric diameter change was used to provide canonical insights into potentially deep or shallow seismic sources. As gas slugs ascend the tube, we observe systematic pressure changes varying with slug size, liquid depth, tube diameter, and liquid viscosity. Gas slugs undergoing an abrupt flow pattern change upon entering a section of significantly increased tube diameter induce a transient pressure decrease in and above the flare and an associated pressure increase below it, which stimulates acoustic and inertial resonant oscillations. When the liquid flow is not dominantly controlled by viscosity, net vertical forces on the apparatus are also detected. The net force is a function of the magnitude of the pressure transients generated and the tube geometry, which dictates where, and hence when, the traveling pressure pulses can couple into the tube. In contrast to interpretations of related volcano-seismic data, where a single downward force is assumed to result from an upward acceleration of the center of mass in the conduit, our experiments suggest that significant downward forces can result from the rapid deceleration of relatively small volumes of downward-moving liquid. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

The U-tube: A novel system for acquiring borehole fluid samples from a deep geologic CO2 sequestration experiment

A novel system has been deployed to obtain geochemical samples of water and gas, at in situ pressure, during a geologic CO 2 sequestration experiment conducted in the Frio brine aquifer in Liberty County, Texas. Project goals required high-frequency recovery of representative and uncontaminated aliquots of a rapidly changing two-phase fluid (supercritical CO 2 and brine) fluid from 1.5 km depth. The data sets collected, using both the liquid and gas portions of the downhole samples, provide insights into the coupled hydrogeochemical issues affecting CO 2 sequestration in brine-filled formations. While the basic premise underlying the U-tube sampler is not new, the system is unique because careful consideration was given to the processing of the recovered two-phase fluids. In particular, strain gauges mounted beneath the high-pressure surface sample cylinders measured the ratio of recovered brine to supercritical CO 2 . A quadrupole mass spectrometer provided real-time gas analysis for perfluorocarbon and noble gas tracers that were injected along with the CO 2 . The U-tube successfully acquired frequent samples, facilitating accurate delineation of the arrival of the CO 2 plume, and on-site analysis revealed rapid changes in geochemical conditions.

Journal of Geophysical Research Solid Earth↗
USGS2006

Seismic characteristics of central Brazil crust and upper mantle: A deep seismic refraction study

[1] A two‐dimensional model of the Brazilian central crust and upper mantle was obtained from the traveltime interpretation of deep seismic refraction data from the Porangatu and Cavalcante lines, each approximately 300 km long. When the lines were deployed, they overlapped by 50 km, forming an E‐W transect approximately 530 km long across the Tocantins Province and western São Francisco Craton. The Tocantins Province formed during the Neoproterozoic when the São Francisco, the Paranapanema, and the Amazon cratons collided, following the subduction of the former Goiás ocean basin. Average crustal V P and V P /V S ratios, Moho topography, and lateral discontinuities within crustal layers suggest that the crust beneath central Brazil can be associated with major geological domains recognized at the surface. The Moho is an irregular interface, between 36 and 44 km deep, that shows evidences of first‐order tectonic structures. The 8.05 and 8.23 km s −1 P wave velocities identify the upper mantle beneath the Porangatu and Cavalcante lines, respectively. The observed seismic features allow for the identification of (1) the crust has largely felsic composition in the studied region, (2) the absence of the mafic‐ultramafic root beneath the Goiás magmatic arc, and (3) block tectonics in the foreland fold‐and‐thrust belt of the northern Brasília Belt during the Neoproterozoic. Seismic data also suggested that the Bouguer gravimetric discontinuities are mainly compensated by differences in mass distribution within the lithospheric mantle. Finally, the Goiás‐Tocantins seismic belt can be interpreted as a natural seismic alignment related to the Neoproterozoic mantle domain.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

The instantaneous rate dependence in low temperature laboratory rock friction and rock deformation experiments

[1] Earthquake occurrence probabilities that account for stress transfer and time-dependent failure depend on the product of the effective normal stress and a lab-derived dimensionless coefficient a . This coefficient describes the instantaneous dependence of fault strength on deformation rate, and determines the duration of precursory slip. Although an instantaneous rate dependence is observed for fracture, friction, crack growth, and low temperature plasticity in laboratory experiments, the physical origin of this effect during earthquake faulting is obscure. We examine this rate dependence in laboratory experiments on different rock types using a normalization scheme modified from one proposed by Tullis and Weeks [1987]. We compare the instantaneous rate dependence in rock friction with rate dependence measurements from higher temperature dislocation glide experiments. The same normalization scheme is used to compare rate dependence in friction to rock fracture and to low-temperature crack growth tests. For particular weak phyllosilicate minerals, the instantaneous friction rate dependence is consistent with dislocation glide. In intact rock failure tests, for each rock type considered, the instantaneous rate dependence is the same size as for friction, suggesting a common physical origin. During subcritical crack growth in strong quartzofeldspathic and carbonate rock where glide is not possible, the instantaneous rate dependence measured during failure or creep tests at high stress has long been thought to be due to crack growth; however, direct comparison between crack growth and friction tests shows poor agreement. The crack growth rate dependence appears to be higher than the rate dependence of friction and fracture by a factor of two to three for all rock types considered.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Frictional strength heterogeneity and surface heat flow: Implications for the strength of the creeping San Andreas fault

Heat flow measurements along much of the San Andreas fault (SAF) constrain the apparent coefficient of friction (??app) of the fault to <0.2, much lower than laboratory-derived friction values for most geologic materials. However, heat flow data are sparse near the creeping section of the SAF, a frictional "asperity" where the fault slips almost exclusively by aseismic creep. We test the hypothesis that the creeping section has a substantially higher or lower ?? app than adjacent sections of the SAF. We use numerical models to explore the effects of faults with spatially and temporally heterogeneous frictional strength on the spatial distribution of surface heat flow. Heat flow from finite length asperities is uniformly lower than predicted by assuming an infinitely long fault. Over geologic time, lateral offset from strike-slip faulting produces heat flow patterns that are asymmetric across the fault and along strike. We explore a range of asperity sizes, slip rates, and displacement histories for comparing predicted spatial patterns of heat flow with existing measurements. Models with ??app ??? 0.1 fit the data best. For most scenarios, heat flow anomalies from a frictional asperity with ??app > 0.2 should be detectable even with the sparse existing observations, implying that ??app for the creeping section is as low as the surrounding SAF. Because the creeping section does not slip in large earthquakes, the mechanism controlling its weakness is not related to dynamic processes resulting from high slip rate earthquake ruptures. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

V P and V S structure of the Yellowstone hot spot from teleseismic tomography: Evidence for an upper mantle plume

The movement of the lithosphere over a stationary mantle magmatic source, often thought to be a mantle plume, explains key features of the 16 Ma Yellowstone–Snake River Plain volcanic system. However, the seismic signature of a Yellowstone plume has remained elusive because of the lack of adequate data. We employ new teleseismic P and S wave traveltime data to develop tomographic images of the Yellowstone hot spot upper mantle. The teleseismic data were recorded with two temporary seismograph arrays deployed in a 500 km by 600 km area centered on Yellowstone. Additional data from nearby regional seismic networks were incorporated into the data set. The V P and V S models reveal a strong low-velocity anomaly from ∼50 to 200 km directly beneath the Yellowstone caldera and eastern Snake River Plain, as has been imaged in previous studies. Peak anomalies are −2.3% for V P and −5.5% for V S . A weaker, anomaly with a velocity perturbation of up to −1.0% V P and −2.5% V S continues to at least 400 km depth. This anomaly dips 30° from vertical, west-northwest to a location beneath the northern Rocky Mountains. We interpret the low-velocity body as a plume of upwelling hot, and possibly wet rock, from the mantle transition zone that promotes small-scale convection in the upper ∼200 km of the mantle and long-lived volcanism. A high-velocity anomaly, 1.2% V P and 1.9% V S , is located at ∼100 to 250 km depth southeast of Yellowstone and may represent a downwelling of colder, denser mantle material.

Idaho, Montana, Nevada, Oregon, Utah, Wyoming↗