Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Travel 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.

137 records · Page 8Linked to original sources
USGS2012

Seismic guided waves trapped in the fault zone of the Landers, California, earthquake of 1992

A mobile seismic array of seven stations was deployed at 11 sites along the fault trace of the M 7.4 Landers earthquake of June 28, 1992, with a maximum offset of 1 km from the trace. We found a distinct wave train with a relatively long period following the S waves that shows up only when both the stations and the events are close to the fault trace. This wave train is interpreted as a seismic guided wave trapped in a low‐velocity fault zone. To study the distribution of amplitude of the guided waves with distance from the fault trace and also their attenuation with travel distance along the fault zone, we eliminated source and recording site effects by the coda normalization method. The normalized amplitudes of guided waves show a spectral peak at 3–4 Hz, which decays sharply with distance from the fault trace. Spectral amplitudes at high frequencies (8–15 Hz) show an opposite trend, increasing with distance from the fault trace. The normalized amplitudes of guided waves at 3–4 Hz also show a systematic decrease with hypocentral distance along the fault zone, from which we infer an apparent Q of 50. In order to confirm the existence of the guided waves, a dense array of 31 stations was deployed at one of the 11 sites. The resultant records revealed unequivocal evidence for the existence of guided waves associated with the fault zone. By modeling the waveforms as S waves trapped in a low‐velocity waveguide sandwiched between two homogeneous half‐spaces with velocity V s = 3.0 km/s, we infer a waveguide width of about 180 m, a shear velocity of 2.0–2.2 km/s, and a Q of ∼50. Hypocenters of aftershocks with clear guided waves show a systematic distribution both laterally and with depth delineating the extent of the low‐velocity fault zone in three dimensions. We find that the zone extends to a depth of at least 10 km. This zone apparently continues to the south across the Pinto Mountain fault because guided waves are observed at stations north of the Pinto Mountain fault for earthquakes with epicenters south of it. On the other hand, the zone appears to be discontinuous at the fault bend located about 20 km north of the mainshock epicenter; guided waves were observed for stations and epicenters which are located on the same sides of the fault bend but not for those on the opposite sides.

California↗
USGS2012

The giant submarine alika debris slide, Mauna Loa, Hawaii

A 4000‐km 2 area of submarine slump and slide deposits along the west flank of Mauna Loa volcano has been mapped with GLORIA side‐scan sonar images, seismic reflection profiles, and new bathymetry. The youngest deposits are two debris avalanche lobes that travelled from their breakaway area near the present shoreline as much as 100 km into the Hawaiian Deep at water depths of 4800 m. The two lobes partly overlap and together are designated the Alika slide. They were derived from the same source area and probably formed in rapid succession. Distinction hummocky topography, marginal levees, and other features on lower slopes (0.3°–0.6°) of these deposits resemble subaerial volcanic debris avalanche deposits such as 1980 Mount St. Helens and suggest high emplacement velocities. The breakaway area for the Alika slide (10°–15° slopes) is characterized by large block slumps, bounded by normal faults, that probably represent multiple subsidence events before, during, and after the debris avalanches. Lower slopes of the slide contain distinctive lobate‐terraced deposits that are interpreted as having been emplaced more slowly, prior to the debris avalanches. Estimated thicknesses of 50–200 m suggest volumes of 200–600 km 3 for the two lobes. The combined volume of the entire slide and slump terrane is probably 1500–2000 km 3 . The slide deposits predate a 13‐ka coral reef and probably postdate the block‐faulted Ninole Basalt, roughly dated as a few hundred thousand years old. The Alika slide, or a similar deposit recognized on GLORIA images further north along the Hawaiian Ridge, probably triggered a giant wave that washed 325 m high on Lanai at about 100 ka. Slumping on Mauna Loa has been most intense adjacent to the large arcuate bend in its southwest rift zone, as the rift zone migrated westward away from the growing Kilauea volcano. Slumping events were probably triggered by seismic activity accompanying dike injection along the rift zone. Such massive slumps, landslides, and distal submarine turbidity flows appear to be widespread on the flanks of Hawaiian volcanoes.

Journal of Geophysical Research↗
USGS2009

Movements of juvenile Gyrfalcons from western and interior Alaska following departure from their natal areas

J uvenile raptors often travel thousands of kilometers from the time they leave their natal areas to the time they enter a breeding population. Documenting movements and identifying areas used by raptors before they enter a breeding population is important for understanding the factors that influence their survival. In North America, juvenile Gyrfalcons ( Falco rusticolus ) are routinely observed outside the species' breeding range during the nonbreeding season, but the natal origins of these birds are rarely known. We used satellite telemetry to track the movements of juvenile Gyrfalcons during their first months of independence. We instrumented nestlings with lightweight satellite transmitters within 10 d of estimated fledging dates on the Seward Peninsula in western Alaska and in Denali National Park (Denali) in interior Alaska. Gyrfalcons spent an average of 41.4 ± 6.1 d (range = 30–50 d) in their natal areas after fledging. The mean departure date from natal areas was 27 August ± 6.4 d. We tracked 15 individuals for an average of 70.5 ± 28.1 d post-departure; Gyrfalcons moved from 105 to 4299 km during this period and tended to move greater distances earlier in the tracking period than later in the tracking period. Gyrfalcons did not establish temporary winter ranges within the tracking period. We identified several movement patterns among Gyrfalcons, including unidirectional long-distance movements, multidirectional long- and short-distance movements, and shorter movements within a local region. Gyrfalcons from the Seward Peninsula remained in western Alaska or flew to eastern Russia with no movements into interior Alaska. In contrast, Gyrfalcons from Denali remained in interior Alaska, flew to northern and western Alaska, or flew to northern Alberta. Gyrfalcons from both study areas tended to move to coastal, riparian, and wetland areas during autumn and early winter. Because juvenile Gyrfalcons dispersed over a large geographic area and across three international boundaries, conservation efforts should focus on both regional and international scales.

Journal of Raptor Research↗
USGS2000

Nuées ardentes of 22 November 1994 at Merapi volcano, Java, Indonesia

Nuées ardentes associated with dome collapse on 22 November 1994, at Merapi volcano traveled to the south–southwest as far as 6.5 km, and collectively accumulated roughly 2.5–3 million cubic meters of deposits. The damaged area comprises 9.5 km 2 and is covered by two nuée ardente facies, a conventional “Merapi-type”, valley-fill block-and-ash flow facies and a pyroclastic surge facies. The proximal deposits reflect the accumulation of dozens of nuées ardentes, with many subsidiary flow units. The distal deposits are more simply organized, as only a few individual events reached to distances >3.5 km. The stratigraphic relationships north of Turgo hill indicate that the surge deposits are a facies of particularly mobile nuées ardentes that also deposited channeled block-and-ash flow facies. They further suggest that the surge facies beyond the channel margins correlate laterally with a finer-grained sublayer locally developed at the base of the block-and-ash flow facies. Eyewitness reports suggest that the emplacement of the block-and-ash flow facies in the distal part of the Boyong river may have followed, by a short time interval, the destruction and deposition of the surge facies at Turgo village. The stratigraphy is in accord with the eyewitness reports. The surge facies was emplaced by a dilute surge current, detached from the same dome-collapse nuée ardente that, as a separate flow unit, subsequently emplaced the distal block-and-ash deposit in the Boyong valley. The detachment occurred at higher elevations, likely at or above the slope break at about 2000 m elevation. This flow separation enabled the surge current to shortcut over the landscape and to emplace its deposit even as the block-and-ash flow continued its tortuous southward movement in the Boyong channel. Dome-collapse nuée ardente activity formed the bulk of the eruption, which was accompanied by virtually no significant vertical summit explosive activity.

Journal of Volcanology and Geothermal Research↗
USGS2020

The 1951 eruption of Mount Lamington, Papua New Guinea: Devastating directed blast triggered by small-scale edifice failure

The catastrophic explosion of Mount Lamington volcano, Papua New Guinea on January 21, 1951 produced a devastating pyroclastic density current (PDC) that knocked down dense tropical rainforest over an area of 230 km 2 and killed approximately 3000 people. We present results of a field reinvestigation of the 1951 PDC deposit combined with an analysis of the available photographs and eyewitness accounts of the eruption first published in the fundamental work of G. A. M. Taylor (1958). We have concluded that the six-days-long pre-climactic activity before the 1951 eruption (which included felt local seismicity, frequent ash-laden explosions of vulcanian type, bulging of the volcano slope accompanied with landslides) was associated with shallow-level intrusion of a highly viscous magma body (cryptodome/dome) of andesitic composition with a volume of approximately 0.01 km 3 . This intrusion destabilized Mount Lamington's prehistoric intra-crater lava dome. On January 21 the destabilized dome gravitationally collapsed and produced a relatively small-volume debris avalanche, the deposit of which was not recognized during Taylor's original investigation. The debris avalanche had a volume of approximately 0.02–0.04 km 3 , travelled a distance (L) of 8.5 km and had the ratio of vertical drop (H) to runout (L) of 0.14. The edifice collapse decompressed the intruding cryptodome and triggered its explosive fragmentation. Photographs of the climactic explosion show that the eruptive cloud initially rose vertically but subsequently collapsed upon the terrain around the vent, and formed a PDC which flowed radially outward. The enhanced northward propagation of the PDC to a maximum distance of 13 km reveals that the northern breach in the ancient crater's high walls influenced the distribution of the deposit. In the studied NE-N-NW sector of the devastated area, in the zone proximal to the volcano, the PDC emplaced a normally graded layer of coarse ash and lapilli mixed in the base with picked-up soil and plant fragments. The layer gradually becomes thinner and finer-grained with distance from the volcano. The PDC deposit has a volume of approximately 0.025 km 3 and consists of approximately 80% juvenile rock fragments derived from the explosively fragmented cryptodome. The remaining 20% consists of accidental clasts derived from the old volcanic edifice. The juvenile material is crystal-rich andesite with a unimodal vesicularity distribution (4 to 36%). The reconstructed eruption sequence, the PDC tree blowdown pattern and characteristics of the PDC deposit are similar to those of catastrophic laterally-directed blasts of volcanoes Bezymianny in 1956, Mount St.Helens in 1980, and Soufriere Hills, Montserrat in 1997. In contrast to the cases of these “classic” lateral blasts, the blast cloud of Lamington was initially vertically-directed before collapsing to produce a PDC. We speculate that the climactic explosion of Mount Lamington was initially vertical because the rupture surface of the triggering sector collapse intersected the apex of the intruding cryptodome (it exposed a subhorizontal surface of the cryptodome apex), while at Bezymianny, Mount St.Helens, and Soufriere Hills the rupture intersected the main body of the cryptodome/dome, and exposed their steeply inclined surfaces.

Journal of Volcanology and Geothermal Research↗
USGS1984

The eruption of Mount Pagan volcano, Mariana Islands, 15 May 1981

A major explosive eruption occurred 15 May 1981 at Mount Pagan Volcano, the larger of two historic eruptive centers on Pagan Island, Mariana Islands. The eruption was preceded by increased numbers of locally felt earthquakes beginning in late March or early April and by new ground cracks, new sublimates, and increased gas emissions. A swarm of felt earthquakes began at 0745h (local time = UCT+10 hours) 15 May, and at 0915 h, closely following a loud sonic boom, a strong plinian column issued from the volcano. The high-altitude ash cloud (at least 13.5 km) travelled south-southeast, but ash and scoria deposits were thickest (> 2 m) in the NW sector of the island because of the prevailing low-altitude southeasterly winds. The early activity of 15 May probably involved magmatic eruption along a fissure system oriented about N10??E. However, the eruption became hydromagmatic, possibly within minutes, and was largely restricted to three long-lived vents. The northernmost of these built a substantial new scoria-ash cinder cone. Flows and air-fall deposits, consisting almost entirely of juvenile material, exceeded 105 ?? 106 m3 in volume (75 ?? 106 m3 of magma) on land and at least 70-100 ?? 606 m3 at sea. An unknown volume was carried away by stratospheric winds. Lithic blocks and juvenile bombs as large as 1 m in diameter were thrown more than 2 km from the summit, and evidence for base-surge was observed in restricted corridors as low as 200 m elevation on the north and south slopes of the volcano. Neither of these events resulted in serious injuries to the 54 residents of the island, nor did the eruption produce serious chemical hazards in their water supply. Weak eruptions occurred during the ensuing month, and some of these were monitored by ground observations, seismic monitoring, and deformation studies. Precursory seismicity and possibly deformation occurred with some of the observed eruptions. More vigorous eruptions were reported by visiting residents in late 1981 and early 1982, but these were of lesser magnitude than the 15 May 1981 event. The 15 May lava is predominantly aa and ranges from 3 to > 30 m in thickness. In composition, it is a high-alumina basalt with small (< 1 mm long) phenocrysts of plagioclase and clinopyroxene (7%) that is more or less typical of basalt of the northern Marianas volcanoes. It contains slightly more SiO2 (52%), K2O, TiO2, and less Al2O3 and CaO than does the basalt of the last eruptive event of Mount Pagan Volcano in 1925. Gas analyses indicate that a large portion of air was introduced into the vent system through the porous volcanic edifice and that the carbon gases were not in equilibrium with the magma or each other. ?? 1984.

Journal of Volcanology and Geothermal Research↗
USGS2007

Three-dimensional P-wave velocity structure derived from local earthquakes at the Katmai group of volcanoes, Alaska

The three-dimensional P-wave velocity structure beneath the Katmai group of volcanoes is determined by inversion of more than 10,000 rays from over 1000 earthquakes recorded on a local 18 station short-period network between September 1996 and May 2001. The inversion is well constrained from sea level to about 6 km below sea level and encompasses all of the Katmai volcanoes; Martin, Mageik, Trident, Griggs, Novarupta, Snowy, and Katmai caldera. The inversion reduced the average RMS travel-time error from 0.22 s for locations from the standard one-dimensional model to 0.13 s for the best three-dimensional model. The final model, from the 6th inversion step, reveals a prominent low velocity zone (3.6–5.0 km/s) centered at Katmai Pass and extending from Mageik to Trident volcanoes. The anomaly has values about 20–25% slower than velocities outboard of the region (5.0–6.5 km/s). Moderately low velocities (4.5–6.0 km/s) are observed along the volcanic axis between Martin and Katmai Caldera. Griggs volcano, located about 10 km behind (northwest of) the volcanic axis, has unremarkable velocities (5.0–5.7 km/s) compared to non-volcanic regions. The highest velocities are observed between Snowy and Griggs volcanoes (5.5–6.5 km/s). Relocated hypocenters for the best 3-D model are shifted significantly relative to the standard model with clusters of seismicity at Martin volcano shifting systematically deeper by about 1 km to depths of 0 to 4 km below sea level. Hypocenters for the Katmai Caldera are more tightly clustered, relocating beneath the 1912 scarp walls. The relocated hypocenters allow us to compare spatial frequency-size distributions ( b -values) using one-dimensional and three-dimensional models. We find that the distribution of b is significantly changed for Martin volcano, which was characterized by variable values (0.8 < b < 2.0) with standard locations and more uniform values (0.8 < b < 1.2) after relocation. Other seismic clusters at Mageik (1.2 < b < 2.2), Trident (0.5 < b < 1.5) and Katmai Caldera (0.8 < b < 1.8) had stable b -values indicating the robustness of the observations. The strong high b -value region at Mageik volcano is mainly associated with an earthquake swarm in October, 1996 that possibly indicates a shallow intrusion or influx of gas. The new velocity and spatial b -value results, in conjunction with prior gravity (Bouguer anomalies up to − 40 mgal) and interferometry (several cm uplift) data, provide strong evidence in favor of partially molten rock at shallow depths beneath the Mageik–Katmai–Novarupta region. Moderately low velocities beneath Martin and Katmai suggest that old, mostly solidified intrusions exist beneath these volcanoes. Higher relative velocities beneath the Griggs and Snowy vents suggest that no magma is resident in the shallow crust beneath these volcanoes.

Journal of Volcanology and Geothermal Research↗
USGS2002

Fault structure and mechanics of the Hayward Fault, California from double-difference earthquake locations

The relationship between small-magnitude seismicity and large-scale crustal faulting along the Hayward Fault, California, is investigated using a double-difference (DD) earthquake location algorithm. We used the DD method to determine high-resolution hypocenter locations of the seismicity that occurred between 1967 and 1998. The DD technique incorporates catalog travel time data and relative P and S wave arrival time measurements from waveform cross correlation to solve for the hypocentral separation between events. The relocated seismicity reveals a narrow, near-vertical fault zone at most locations. This zone follows the Hayward Fault along its northern half and then diverges from it to the east near San Leandro, forming the Mission trend. The relocated seismicity is consistent with the idea that slip from the Calaveras Fault is transferred over the Mission trend onto the northern Hayward Fault. The Mission trend is not clearly associated with any mapped active fault as it continues to the south and joins the Calaveras Fault at Calaveras Reservoir. In some locations, discrete structures adjacent to the main trace are seen, features that were previously hidden in the uncertainty of the network locations. The fine structure of the seismicity suggests that the fault surface on the northern Hayward Fault is curved or that the events occur on several substructures. Near San Leandro, where the more westerly striking trend of the Mission seismicity intersects with the surface trace of the (aseismic) southern Hayward Fault, the seismicity remains diffuse after relocation, with strong variation in focal mechanisms between adjacent events indicating a highly fractured zone of deformation. The seismicity is highly organized in space, especially on the northern Hayward Fault, where it forms horizontal, slip-parallel streaks of hypocenters of only a few tens of meters width, bounded by areas almost absent of seismic activity. During the interval from 1984 to 1998, when digital waveforms are available, we find that fewer than 6.5% of the earthquakes can be classified as repeating earthquakes, events that rupture the same fault patch more than one time. These most commonly are located in the shallow creeping part of the fault, or within the streaks at greater depth. The slow repeat rate of 2–3 times within the 15-year observation period for events with magnitudes around M = 1.5 is indicative of a low slip rate or a high stress drop. The absence of microearthquakes over large, contiguous areas of the northern Hayward Fault plane in the depth interval from ∼5 to 10 km and the concentrations of seismicity at these depths suggest that the aseismic regions are either locked or retarded and are storing strain energy for release in future large-magnitude earthquakes.

California↗
USGS2021

Deepwater debrites and linked megaturbidites in confined basins: An example from the Onnuri Basin, East Sea of Korea

We analyzed data from seven piston cores, multi-channel seismic-reflection (MCS) and chirp profiles, and multibeam echosounder (MBES) data to study the distribution, emplacement time, sedimentary facies, and depositional processes of sediment-gravity-flow deposits in the Onnuri Basin, a confined basin in the East Sea. These data reveal that debris flows have traveled ca. 30 km downslope, forming a seismic facies consisting of stacked, wedge-shaped, transparent units separated by high-amplitude continuous reflectors. Analysis of piston cores shows three distinct sedimentary units, throughout the basin. The lowest unit, I, is a debrite containing numerous mud clasts of varying size and color distributed in a mud-rich matrix; it is absent over elevated basinal highs or ridges, such as the Onnuri Ridge, suggesting that local topography controls its distribution. The debrite forms a recognizable acoustically transparent layer on subbottom chirp profiles (av. 7 m thick), covers approximately 500 km 2 , and has an estimated volume of ∼ 3.5 km 3 . The overlying unit, II, contains normally graded beds composed of massive sand, laminated and cross-laminated sand and silt, and a thick cap of structureless mud. This unit is interpreted to be a megaturbidite deposited from turbidity currents that originated from the flow transformation of debris flows on the upper continental slope. The megaturbidite covers the entire basin (at least 650 km 2 ), and has an average thickness of 2.8 m (maximum thickness of 4.35 m), and comprises a volume of 1.8 km 3 . Variations in grain size and sedimentary structures suggest that the megaturbidite was deposited by progressively waning flows that reflected off basin flanks and ridges. The thick (up to 3.65 m) structureless mud cap further indicates deposition in a confined basin. The sharp basal contact, together with the lack of hemipelagic sediments between debrite and overlying megaturbidite, suggest that both were deposited during the same flow event, likely to have originated from a single catastrophic slope failure. Collapsing slide material evolved into a debris flow, from which a turbidite formed by dilution of the debris flow. Radiocarbon dates suggest that the slope failure occurred about 13–11 ka, a time when sea level was ca. 50 m lower than at the present day. Hemipelagic sediments in the topmost unit, III-2, above the megaturbidite indicate that the basin has been stable since ca. 11 ka. We provide robust evidence that submarine slope failures evolve downslope into slides, debris flows, and finally, thick megaturbidites. This contribution highlights the importance of seafloor morphology on the distribution and stratigraphy of submarine flows in confined basins.

Onnuri Basin↗
USGS2001

Upper crustal structure in Puget Lowland, Washington: Results from the 1998 Seismic Hazards Investigation in Puget Sound

A new three-dimensional (3-D) model shows seismic velocities beneath the Puget Lowland to a depth of 11 km. The model is based on a tomographic inversion of nearly one million first-arrival travel times recorded during the 1998 Seismic Hazards Investigation in Puget Sound (SHIPS), allowing higher-resolution mapping of subsurface structures than previously possible. The model allows us to refine the subsurface geometry of previously proposed faults (e.g., Seattle, Hood Canal, southern Whidbey Island, and Devils Mountain fault zones) as well as to identify structures (Tacoma, Lofall, and Sequim fault zones) that warrant additional study. The largest and most important of these newly identified structures lies along the northern boundary of the Tacoma basin; we informally refer to this structure here as the Tacoma fault zone. Although tomography cannot provide information on the recency of motion on any structure, Holocene earthquake activity on the Tacoma fault zone is suggested by seismicity along it and paleoseismic evidence for abrupt uplift of tidal marsh deposits to its north. The tomography reveals four large, west to northwest trending low-velocity basins (Tacoma, Seattle, Everett, and Port Townsend) separated by regions of higher velocity ridges that are coincident with fault-bounded uplifts of Eocene Crescent Formation basalt and pre-Tertiary basement. The shapes of the basins and uplifts are similar to those observed in gravity data; gravity anomalies calculated from the 3-D tomography model are in close agreement with the observed anomalies. In velocity cross sections the Tacoma and Seattle basins are asymmetric: the basin floor dips gently toward a steep boundary with the adjacent high-velocity uplift, locally with a velocity "overhang" that suggests a basin vergent thrust fault boundary. Crustal fault zones grow from minor folds into much larger structures along strike. Inferred structural relief across the Tacoma fault zone increases by several kilometers westward along the fault zone to Lynch Cove, where we interpret it as a zone of south vergent faulting overthrusting Tacoma basin. In contrast, structural relief along the Seattle fault zone decreases west of Seattle, which we interpret as evidence that the N-S directed compression is being accommodated by slip transfer between the Seattle and Tacoma fault zones. Together, the Tacoma and Seattle fault zones raise the Seattle uplift, one of a series of east-west trending, pop-up structures underlying Puget Lowland from the Black Hills to the San Juan Islands.

Washington↗
USGS1998

Closing the gap between regional and global travel time tomography

Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al . [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [ Spakman and Bijwaard , 1998]. Second, we exploit the global data set of Engdahl et al . [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used ( P , pP , and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

Journal of Geophysical Research B: Solid Earth↗