USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research: Solid Earth”

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

Why earthquakes correlate weakly with the solid Earth tides: Effects of periodic stress on the rate and probability of earthquake occurrence

We provide an explanation why earthquake occurrence does not correlate well with the daily solid Earth tides. The explanation is derived from analysis of laboratory experiments in which faults are loaded to quasiperiodic failure by the combined action of a constant stressing rate, intended to simulate tectonic loading, and a small sinusoidal stress, analogous to the Earth tides. Event populations whose failure times correlate with the oscillating stress show two modes of response; the response mode depends on the stressing frequency. Correlation that is consistent with stress threshold failure models, e.g., Coulomb failure, results when the period of stress oscillation exceeds a characteristic time tn; the degree of correlation between failure time and the phase of the driving stress depends on the amplitude and frequency of the stress oscillation and on the stressing rate. When the period of the oscillating stress is less than tn, the correlation is not consistent with threshold failure models, and much higher stress amplitudes are required to induce detectable correlation with the oscillating stress. The physical interpretation of tn is the duration of failure nucleation. Behavior at the higher frequencies is consistent with a second-order dependence of the fault strength on sliding rate which determines the duration of nucleation and damps the response to stress change at frequencies greater than 1/tn. Simple extrapolation of these results to the Earth suggests a very weak correlation of earthquakes with the daily Earth tides, one that would require >13,000 earthquakes to detect. On the basis of our experiments and analysis, the absence of definitive daily triggering of earthquakes by the Earth tides requires that for earthquakes, tn exceeds the daily tidal period. The experiments suggest that the minimum typical duration of earthquake nucleation on the San Andreas fault system is ???1 year.

Journal of Geophysical Research B: Solid Earth↗
USGS2000

Geochemical evidence for a magmatic CO 2 degassing event at Mammoth Mountain, California, September-December 1997

Recent time series soil CO 2 concentration data from monitoring stations in the vicinity of Mammoth Mountain, California, reveal strong evidence for a magmatic degassing event during the fall of 1997 lasting more than 2 months. Two sensors at Horseshoe Lake first recorded the episode on September 23, 1997, followed 10 days later by a sensor on the north flank of Mammoth Mountain. Direct degassing from shallow intruding magma seems an implausible cause of the degassing event, since the gas released at Horseshoe Lake continued to be cold and barren of other magmatic gases, except for He. We suggest that an increase in compressional strain on the area south of Mammoth Mountain driven by movement of major fault blocks in Long Valley caldera may have triggered an episode of increased degassing by squeezing additional accumulated CO 2 from a shallow gas reservoir to the surface along faults and other structures where it could be detected by the CO 2 monitoring network. Recharge of the gas reservoir by CO 2 emanating from the deep intrusions that probably triggered deep long-period earthquakes may also have contributed to the degassing event. The nature of CO 2 discharge at the soil-air interface is influenced by the porous character of High Sierra soils and by meteorological processes. Solar insolation is the primary source of energy for the Earth atmosphere and plays a significant role in most diurnal processes at the Earth surface. Data from this study suggest that external forcing due largely to local orographic winds influences the fine structure of the recorded CO 2 signals.

California↗
USGS2012

An aeromagnetic study of Mount St. Helens

Aeromagnetic data from surveys flown by the U.S. Geological Survey over Mount St. Helens, Washington, before and after the climactic May 18, 1980, eruption were used to determine the bulk magnetic properties of the volcano and to delineate a buried source. We assumed that most of the edifice of preeruption and posteruption Mount St. Helens has a magnetization direction near the present earth's field of 69° inclination and 20° declination and calculated its intensity as 4.1 A/m. After subtraction of magnetic anomalies due to topogaphy magnetized with this direction and intensity, the preeruption and posteruption surveys revealed nearly identical residual magnetic highs and lows, indicative that their sources were not altered or removed by the May 1980 eruption. The residual highs were explained by a 200-m-deep source lying mostly within the edifice of Mount St. Helens. The source could be terrain that predates Mount St. Helens, such as a buried ridge or a cone or a valley filled with lava. We calculated the magnetization of the material removed by the May 1980 eruption and found its intensity to be 4.2 A/m in a direction near the present earth's field, similar to that assumed for the volcano as a whole. This similar result confirmed the validity of the magnetization vector assumed for the entire edifice.

Washington↗
USGS2004

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

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

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Paleomagnetism and K-Ar ages of volcanic rocks from Long Valley caldera, California

Paleomagnetic measurements and K‐Ar age determinations on volcanic rocks from Long Valley caldera, California, have enabled further refinement of eruptive activity within this large silicic volcanic center. K‐Ar age determinations show that postcaldera volcanic eruptions began 0.73 m.y. ago and continued periodically until about 50,000 years ago. The eruptions were not temporally random but tended to occur in distinct episodes separated by periods of quiescence. Volcanism in the western half of the caldera was particularly intense between 0.15 and 0.50 m.y. ago, when many units ranging in composition from basalt to rhyolite were erupted. An average of the virtual geomagnetic poles for 33 units from the caldera yields a paleomagnetic pole at 89.7°N, 138.4°E, (α 95 = 5.1°), which is indistinguishable from the earth's rotational axis. The ancient geomagnetic field dispersion about this mean pole is 16.0°, with upper and lower limits of 19.3° and 13.6°, respectively.

California↗
USGS1997

Timing of detachment faulting in the Bullfrog Hills and Bare Mountain area, southwest Nevada: Inferences from 40Ar/39Ar, K-Ar, U-Pb, and fission track thermochronology

Crustal extension in the Bullfrog Hills and Bare Mountain area of southwest Nevada is associated with movement along a regional detachment fault. Normal faulting in the upper plate and rapid cooling (denudation) of the lower plate were coeval with Miocene silicic volcanism and with west-northwest transport along the detachment fault. A west-northwest progression of tilting along upper plate normal faults is indicated by ages of the volcanic rocks in relation to angular unconformities. Near the breakaway, tilting in the upper plate occurred between 12.7 and 11.6 Ma, continued less strongly past 10.7 Ma, and was over by 8.2 Ma. Ten to 20 km west of the breakaway, tilting occurred between 10.7 and 10.33 Ma, continued less strongly after 10.33 Ma, and was over by 8.1 Ma. The cooling histories of the lower plate metamorphic rocks were determined by thermochronologic dating methods: K-Ar and 40 Ar/ 39 Ar on muscovite, biotite, and hornblende, 40 Ar/ 39 Ar on K-feldspar, U-Pb on apatite, zircon, and sphene, and fission track on apatite, zircon, and sphene. Lower plate rocks 10 km west of the breakaway cooled slowly from Early Cretaceous lower-amphibolite facies conditions through 350±50° to 300±50°C between 57 and 38 Ma, then cooled rapidly from 205±50° to 120±50°C between 12.6±1.6 and 11.1±1.9 Ma. Lower plate rocks 20 km west of the breakaway cooled slowly from Early Cretaceous upper-amphibolite facies conditions through 500±50°C at 78–67 Ma, passed through 350±50° to 300±50°C between 16.3±0.4 and 10.5±0.3 Ma, then cooled rapidly from 285±50° to 120±50°C between 10.2 and 8.6 Ma. Upper plate tilting and rapid cooling (denudation) of the lower plate occurred simultaneously in the respective areas. The early slow-cooling part of the lower plate thermal histories was probably related to erosion at the Earth's surface, which stripped off about 9 km of material in 50 to 100 m.y. The results indicate an initial fault dip ≥30° and a 12 mm yr −1 west-northwest migration of the locus of rapid tilting in the upper plate.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Age and character of basaltic rocks of the Yucca Mountain region, southern Nevada

Volcanism in the Yucca Mountain region of southern Nevada in the last 5 m.y. is restricted to moderate-to-small volumes of subalkaline basaltic magmas, produced during at least 6 intervals, and spanning an age range from 4.6 Ma to about 125 ka. Where paleomagnetic evidence is available, the period of volcanism at individual eruptive centers apparently was geologically short-lived, even where multiple eruptions involved different magma types. K-Ar studies are consistent with most other geochronologic information, such as the minimum ages of exposure-dating techniques, and show no evidence of renewed volcanism after a significant quiescence at any of the centers in the Yucca Mountain region. A volcanic recurrence interval of 860 ± 350 kyr is computed from a large K-Ar data set and an evaluation of their uncertainties. Monte Carlo error propagations demonstrate the validity of uncertainties obtained for weighted-mean ages when modified using the goodness of fit parameter, MSWD. Elevated 87 Sr/ 86 Sr initial ratios (Sr i ) in the basalts, nearly constant at 0.707, combined with low SiO 2 and Rb/Sr ratios indicate a subcontinental, lithospheric mantle source, previously enriched in radiogenic Sr and depleted in Rb. Beginning with eruptions of the most voluminous eruptive center, the newly dated Pliocene Thirsty Mountain volcano, basaltic magmas have decreased in eruptive volume, plagioclase-phenocryst content, various trace element ratios, and TiO 2 , while increasing in light rare earth elements, U, Th, P 2 O 5 , and light REE/heavy REE ratios. These time-correlated changes are consistent with either increasing depths of melting or a decreasing thermal gradient in the Yucca Mountain region during the last 5 m.y.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Imaging the mantle beneath Iceland using integrated seismological techniques

Using a combination of body wave and surface wave data sets to reveal the mantle plume and plume head, this study presents a tomographic image of the mantle structure beneath Iceland to 400 km depth. Data comes primarily from the PASSCAL-HOTSPOT deployment of 30 broadband instruments over a period of 2 years, and is supplemented by data from the SIL and ICEMELT networks. Three sets of relative teleseismic body wave arrival times are generated through cross correlation: S and SKS arrivals at 0.03–0.1 Hz, and P and PKIKP arrivals at 0.03–0.1 and 0.8–2.0 Hz. Prior to inversion the crustal portion of the travel time anomalies is removed using the crustal model ICECRTb. This step has a significant effect on the mantle velocity variations imaged down to a depth of ∼250 km. Inversion of relative arrival times only provides information on lateral velocity variations. Surface waves are therefore used to provide absolute velocity information for the uppermost mantle beneath Iceland. The average wave number for the Love wave fundamental mode at 0.020 and 0.024 Hz is measured and used to invert for the average S velocity. Combination of the body wave and surface wave information reveals a predominantly horizontal low-velocity anomaly extending from the Moho down to ∼250 km depth, interpreted as a plume head. Below the plume head a near-cylindrical low-velocity anomaly with a radius of ∼100 km and peak V P and V S anomalies of −2% and −4%, respectively, extends down to the maximum depth of resolution at 400 km. Within the plume head, in the uppermost mantle above the core of the plume, there is a relatively high velocity with a maximum V P and V S anomaly of +2%. This high-velocity anomaly may be the result of the extreme degree of melt extraction necessary to generate the thick (46 km) crust in central Iceland. Comparison of the plume volumetric flux implied by our images, the crustal generation rate, and the degree of melting suggested by rare earth element inversions, suggests that (1) mantle material must be flowing horizontally away from the plume core faster than the overlying lithosphere and (2) the bulk of the plume material does not participate in melting beneath Iceland.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Seismically observed seiching in the Panama Canal

A large portion of the seismic noise spectrum is dominated by water wave energy coupled into the solid Earth. Distinct mechanisms of water wave induced ground motions are distinguished by their spectral content. For example, cultural noise is generally <1 s period, microseisms dominate the seismic spectrum from periods of 2 to 20 s, and the Earth's "hum" is in the range of 50 to 600 s. We show that in a large lake in the Panama Canal there is an additional source of long-period noise generated by standing water waves, seiches, induced by disturbances such as passing ships and wind pressure. We compare seismic waveforms to water level records and relate these observations to changes in local tilt and gravity due to an oscillating seiche. The methods and observations discussed in this paper provide a first step toward quantifying the impact of water inundation as recorded by seismometers. This type of quantified understanding of water inundation will help in future estimates of similar phenomena such as the seismic observations of tsunami impact. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

The strength and rheology of methane clathrate hydrate

Methane clathrate hydrate (structure I) is found to be very strong, based on laboratory triaxial deformation experiments we have carried out on samples of synthetic, high-purity, polycrystalline material. Samples were deformed in compressional creep tests (i.e., constant applied stress, ??), at conditions of confining pressure P = 50 and 100 MPa, strain rate 4.5 ?? 10-8 ??? ?? ??? 4.3 ?? 10-4 s-1, temperature 260 ??? T ??? 287 K, and internal methane pressure 10 ??? PCH4 ??? 15 MPa. At steady state, typically reached in a few percent strain, methane hydrate exhibited strength that was far higher than expected on the basis of published work. In terms of the standard high-temperature creep law, ?? = A??ne-(E*+PV*)/RT the rheology is described by the constants A = 108.55 MPa-n s-1, n = 2.2, E* = 90,000 J mol-1, and V* = 19 cm3 mol-1. For comparison at temperatures just below the ice point, methane hydrate at a given strain rate is over 20 times stronger than ice, and the contrast increases at lower temperatures. The possible occurrence of syntectonic dissociation of methane hydrate to methane plus free water in these experiments suggests that the high strength measured here may be only a lower bound. On Earth, high strength in hydrate-bearing formations implies higher energy release upon decomposition and subsequent failure. In the outer solar system, if Titan has a 100-km-thick near-surface layer of high-strength, low-thermal conductivity methane hydrate as has been suggested, its interior is likely to be considerably warmer than previously expected.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

A dislocation model of strain accumulation and release at a subduction zone

Strain accumulation and release at a subduction zone are attributed to stick slip on the main thrust zone and steady aseismic slip on the remainder of the plate interface. This process can be described as a superposition of steady state subduction and a repetitive cycle of slip on the main thrust zone, consisting of steady normal slip at the plate convergence rate plus occasional thrust events that recover the accumulated normal slip. Because steady state subduction does not contribute to the deformation at the free surface, deformation observed there is completely equivalent to that produced by the slip cycle alone. The response to that slip is simply the response of a particular earth model to embedded dislocations. For a purely elastic earth model, the deformation cycle consists of a coseismic offset followed by a linear‐in‐time recovery to the initial value during the interval between earthquakes. For an elastic‐viscoelastic earth model (elastic lithosphere over a viscoelastic asthenosphere), the postearthquake recovery is not linear in time. Records of local uplift as a function of time indicate that the long‐term postseismic recovery is approximately linear, suggesting that elastic earth models are adequate to describe the deformation cycle. However, the deformation predicted for a simple elastic half‐space earth model does not reproduce the deformation observed along the subduction zones in Japan at all well if stick slip is restricted to the main thrust zone. As recognized earlier by Shimazaki, Seno, and Kato, the uplift profiles could be explained if stick slip were postulated to extend along the plate interface beyond the main thrust zone to a depth of perhaps 100 km, but independent evidence suggests that stick slip at such depths is unlikely.

Journal of Geophysical Research B: Solid Earth↗
USGS2019

Evaluation of temporally correlated noise in global navigation satellite system time series: Geodetic monument performance

Estimates of background noise of Global Positioning System‐derived time series of positions for 740 sites in the western United States are examined. These data consist of daily epochs of three components of displacements that are at least 9.75 years long within the interval between 2000 and 2018. We find that these time series have significant temporal correlations that could be represented as a combination of white, flicker, random‐walk, and band‐pass filtered noise. From this noise model, two other metrics are computed: the root‐mean‐square of seasonal noise, that is, the integrated power spectrum between 0.5 and 2 cycles per year, and the standard error in position rate for a 10‐year‐long time series. These two metrics are used to evaluate potential correlations with different geographic regions and with different methods of construction of monuments used to attach the Global Positioning System antenna to the Earth's surface. The sites with the lowest noise, both in terms of rate error and seasonal root‐mean‐square, are located in semiarid regions east of the rain shadow provided by the Cascade and Sierra Nevada mountain ranges. In addition, according to statistical rank tests, monuments known as drilled‐braced monuments perform 30% to 50% better than other monument types (buildings, boreholes, piers, etc.) in terms of having smaller rate errors and lower seasonal noise.

Journal of Geophysical Research B: Solid Earth↗
USGS2014

Slip rates and spatially variable creep on faults of the northern San Andreas system inferred through Bayesian inversion of Global Positioning System data

Fault creep, depending on its rate and spatial extent, is thought to reduce earthquake hazard by releasing tectonic strain aseismically. We use Bayesian inversion and a newly expanded GPS data set to infer the deep slip rates below assigned locking depths on the San Andreas, Maacama, and Bartlett Springs Faults of Northern California and, for the latter two, the spatially variable interseismic creep rate above the locking depth. We estimate deep slip rates of 21.5 ± 0.5, 13.1 ± 0.8, and 7.5 ± 0.7 mm/yr below 16 km, 9 km, and 13 km on the San Andreas, Maacama, and Bartlett Springs Faults, respectively. We infer that on average the Bartlett Springs fault creeps from the Earth's surface to 13 km depth, and below 5 km the creep rate approaches the deep slip rate. This implies that microseismicity may extend below the locking depth; however, we cannot rule out the presence of locked patches in the seismogenic zone that could generate moderate earthquakes. Our estimated Maacama creep rate, while comparable to the inferred deep slip rate at the Earth's surface, decreases with depth, implying a slip deficit exists. The Maacama deep slip rate estimate, 13.1 mm/yr, exceeds long-term geologic slip rate estimates, perhaps due to distributed off-fault strain or the presence of multiple active fault strands. While our creep rate estimates are relatively insensitive to choice of model locking depth, insufficient independent information regarding locking depths is a source of epistemic uncertainty that impacts deep slip rate estimates.

California↗
USGS2002

Plume-driven plumbing and crustal formation in Iceland

Through combination of surface wave and body wave constraints we derive a three-dimensional (3-D) crustal S velocity model and Moho map for Iceland. It reveals a vast plumbing system feeding mantle plume melt into upper crustal magma chambers where crustal formation takes place. The method is based on the partitioned waveform inversion to which we add additional observations. Love waves from six local events recorded on the HOTSPOT- SIL networks are fitted, S n travel times from the same events measured, previous observations of crustal thickness are added, and all three sets of constraints simultaneously inverted for our 3-D model. In the upper crust (0–15 km) an elongated low-velocity region extends along the length of the Northern, Eastern and Western Neovolcanic Zones. The lowest velocities (−7%) are found at 5–10 km below the two most active volcanic complexes: Hekla and Bárdarbunga-Grímsvötn. In the lower crust (>15 km) the low-velocity region can be represented as a vertical cylinder beneath central Iceland. The low-velocity structure is interpreted as the thermal halo of pipe work which connects the region of melt generation in the uppermost mantle beneath central Iceland to active volcanoes along the neovolcanic zones. Crustal thickness in Iceland varies from 15–20 km beneath the Reykjanes Peninsula, Krafla and the extinct Snæfellsnes rift zone, to 46 km beneath central Iceland. The average crustal thickness is 29 km. The variations in thickness can be explained in terms of the temporal variation in plume productivity over the last ∼20 Myr, the Snæfellsnes rift zone being active during a minimum in plume productivity. Variations in crustal thickness do not depart significantly from an isostatically predicted crustal thickness. The best fit linear isostatic relation implies an average density jump of 4% across the Moho. Rare earth element inversions of basalt compositions on Iceland suggest a melt thickness (i.e., crustal thickness) of 15–20 km, given passive upwelling. The observed crustal thickness of up to 46 km implies active fluxing of source material through the melt zone by the mantle plume at up to 3 times the passive rate.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Volcanic plume height measured by seismic waves based on a mechanical model

In August 2008 an unmonitored, largely unstudied Aleutian volcano, Kasatochi, erupted catastrophically. Here we use seismic data to infer the height of large eruptive columns such as those of Kasatochi based on a combination of existing fluid and solid mechanical models. In so doing, we propose a connection between a common, observable, short-period seismic wave amplitude to the physics of an eruptive column. To construct a combined model, we estimate the mass ejection rate of material from the vent on the basis of the plume height, assuming that the height is controlled by thermal buoyancy for a continuous plume. Using the estimated mass ejection rate, we then derive the equivalent vertical force on the Earth through a momentum balance. Finally, we calculate the far-field surface waves resulting from the vertical force. The model performs well for recent eruptions of Kasatochi and Augustine volcanoes if v , the velocity of material exiting the vent, is 120-230 m s -1 . The consistency between the seismically inferred and measured plume heights indicates that in these cases the far-field ~1 s seismic energy radiated by fluctuating flow in the volcanic jet during the eruption is a useful indicator of overall mass ejection rates. Thus, use of the model holds promise for characterizing eruptions and evaluating ash hazards to aircraft in real time on the basis of far-field short-period seismic data. This study emphasizes the need for better measurements of eruptive plume heights and a more detailed understanding of the full spectrum of seismic energy radiated coeruptively.

Journal of Geophysical Research B: Solid Earth↗
USGS2014

Triggering and modulation of geyser eruptions in Yellowstone National Park by earthquakes, earth tides, and weather

We analyze intervals between eruptions (IBEs) data acquired between 2001 and 2011 at Daisy and Old Faithful geysers in Yellowstone National Park. We focus our statistical analysis on the response of these geysers to stress perturbations from within the solid earth (earthquakes and earth tides) and from weather (air pressure and temperature, precipitation, and wind). We conclude that (1) the IBEs of these geysers are insensitive to periodic stresses induced by solid earth tides and barometric pressure variations; (2) Daisy (pool geyser) IBEs lengthen by evaporation and heat loss in response to large wind storms and cold air; and (3) Old Faithful (cone geyser) IBEs are not modulated by air temperature and pressure variations, wind, and precipitation, suggesting that the subsurface water column is decoupled from the atmosphere. Dynamic stress changes of 0.1−0.2 MPa resulting from the 2002 M-7.9 Denali, Alaska, earthquake surface waves caused a statistically significant shortening of Daisy geyser's IBEs. Stresses induced by other large global earthquakes during the study period were at least an order of magnitude smaller. In contrast, dynamic stresses of >0.5 MPa from three large regional earthquakes in 1959, 1975, and 1983 caused lengthening of Old Faithful's IBEs. We infer that most subannual geyser IBE variability is dominated by internal processes and interaction with other geysers. The results of this study provide quantitative bounds on the sensitivity of hydrothermal systems to external stress perturbations and have implications for studying the triggering and modulation of volcanic eruptions by external forces.

Wyoming↗
USGS1997

Seismic-geodynamic constraints on three-dimensional structure, vertical flow, and heat transfer in the mantle

Joint inversions of seismic and geodynamic data are carried out in which we simultaneously constrain global-scale seismic heterogeneity in the mantle as well as the amplitude of vertical mantle flow across the 670 km seismic discontinuity. These inversions reveal the existence of a family of three-dimensional (3-D) mantle models that satisfy the data while at the same time yielding predictions of layered mantle flow. The new 3-D mantle models we obtain demonstrate that the buoyancy forces due to the undulations of the 670 km phase-change boundary strongly inhibit the vertical flow between the upper and lower mantle. The strong stabilizing effect of the 670 km topography also has an important impact on the predicted dynamic topography of the Earth's solid surface and on the surface gravity anomalies. The new 3-D models that predict strongly or partially layered mantle flow provide essentially identical fits to the global seismic data as previous models that have, until now, predicted only whole-mantle flow. The convective vertical transport of heat across the mantle predicted on the basis of the new 3-D models shows that the heat flow is a minimum at 1000 km depth. This suggests the presence at this depth of a globally defined horizon across which the pattern of lateral heterogeneity changes rapidly.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Influence of welded boundaries in anelastic media on energy flow, and characteristics of P, S-I, and S-II waves: Observational evidence for inhomogeneous body waves in low-loss solids

A general computer code, developed to calculate anelastic reflection-refraction coefficients, energy flow, and the physical characteristics for general P , S -I, and S -II waves, quantitatively describes physical characteristics for wave fields in anelastic media that do not exist in elastic media. Consideration of wave fields incident on boundaries between anelastic media shows that scattered wave fields experience reductions in phase and energy speeds, increases in maximum attenuation and Q −1 , and directions of maximum energy flow distinct from phase propagation. Each of these changes in physical characteristics are shown to vary with angle of incidence. Finite relaxation times for anelastic media result in energy flow due to interaction of superimposed radiation fields and contribute to energy flow across anelastic boundaries for all angles of incidence. Agreement of theoretical and numerical results with laboratory measurements argues for the validity of the theoretical and numerical formulations incorporating inhomogeneous wave fields. The agreement attests to the applicability of the model and helps confirm the existence of inhomogeneous body waves and their associated set of distinct physical characteristics in the earth. The existence of such body waves in layered, low-loss anelastic solids implies the need to reformulate some seismological models of the earth. The exact anelastic formulation for a liquid-solid interface with no low-loss approximations predicts the existence of a range of angles of incidence or an anelastic Rayleigh window, through which significant amounts of energy are transmitted across the boundary. The window accounts for the discrepancy apparent between measured reflection data presented in early textbooks and predictions based on classical elasticity theory. Characteristics of the anelastic Rayleigh window are expected to be evident in certain sets of wide-angle, ocean-bottom reflection data and to be useful in estimating Q −1 for some ocean bottom reflectors.

Journal of Geophysical Research B: Solid Earth↗