Search USGSSearch

SEARCH · Search USGS

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 379 records · Page 21Linked to original sources

A domain decomposition approach to implementing fault slip in finite-element models of quasi-static and dynamic crustal deformation

We employ a domain decomposition approach with Lagrange multipliers to implement fault slip in a finite-element code, PyLith, for use in both quasi-static and dynamic crustal deformation applications. This integrated approach to solving both quasi-static and dynamic simulations leverages common finite-element data structures and implementations of various boundary conditions, discretization schemes, and bulk and fault rheologies. We have developed a custom preconditioner for the Lagrange multiplier portion of the system of equations that provides excellent scalability with problem size compared to conventional additive Schwarz methods. We demonstrate application of this approach using benchmarks for both quasi-static viscoelastic deformation and dynamic spontaneous rupture propagation that verify the numerical implementation in PyLith.

Journal of Geophysical Research B: Solid Earth

Seismic evidence for a slab tear at the Puerto Rico Trench

The fore-arc region of the northeast Caribbean plate north of Puerto Rico and the Virgin Islands has been the site of numerous seismic swarms since at least 1976. A 6 month deployment of five ocean bottom seismographs recorded two such tightly clustered swarms, along with additional events. Joint analyses of the ocean bottom seismographs and land-based seismic data reveal that the swarms are located at depths of 50–150 km. Focal mechanism solutions, found by jointly fitting P wave first-motion polarities and S/P amplitude ratios, indicate that the broadly distributed events outside the swarm generally have strike- and dip-slip mechanisms at depths of 50–100 km, while events at depths of 100–150 km have oblique mechanisms. A stress inversion reveals two distinct stress regimes: The slab segment east of 65°W longitude is dominated by trench-normal tensile stresses at shallower depths (50–100 km) and by trench-parallel tensile stresses at deeper depths (100–150 km), whereas the slab segment west of 65°W longitude has tensile stresses that are consistently trench normal throughout the depth range at which events were observed (50–100 km). The simple stress pattern in the western segment implies relatively straightforward subduction of an unimpeded slab, while the stress pattern observed in the eastern segment, shallow trench-normal tension and deeper trench-normal compression, is consistent with flexure of the slab due to rollback. These results support the hypothesis that the subducting North American plate is tearing at or near these swarms. The 35 year record of seismic swarms at this location and the recent increase in seismicity suggest that the tear is still propagating.

Puerto Rico

Foreshocks during the nucleation of stick-slip instability

We report on laboratory experiments which investigate interactions between aseismic slip, stress changes, and seismicity on a critically stressed fault during the nucleation of stick-slip instability. We monitor quasi-static and dynamic changes in local shear stress and fault slip with arrays of gages deployed along a simulated strike-slip fault (2 m long and 0.4 m deep) in a saw cut sample of Sierra White granite. With 14 piezoelectric sensors, we simultaneously monitor seismic signals produced during the nucleation phase and subsequent dynamic rupture. We observe localized aseismic fault slip in an approximately meter-sized zone in the center of the fault, while the ends of the fault remain locked. Clusters of high-frequency foreshocks (M w ~ −6.5 to −5.0) can occur in this slowly slipping zone 5–50 ms prior to the initiation of dynamic rupture; their occurrence appears to be dependent on the rate at which local shear stress is applied to the fault. The meter-sized nucleation zone is generally consistent with theoretical estimates, but source radii of the foreshocks (2 to 70 mm) are 1 to 2 orders of magnitude smaller than the theoretical minimum length scale over which earthquake nucleation can occur. We propose that frictional stability and the transition between seismic and aseismic slip are modulated by local stressing rate and that fault sections, which would typically slip aseismically, may radiate seismic waves if they are rapidly stressed. Fault behavior of this type may provide physical insight into the mechanics of foreshocks, tremor, repeating earthquake sequences, and a minimum earthquake source dimension.

Journal of Geophysical Research B: Solid Earth

Eruptions at Lone Star Geyser, Yellowstone National Park, USA, part 1: energetics and eruption dynamics

Geysers provide a natural laboratory to study multiphase eruptive processes. We present results from a four&ndash;day experiment at Lone Star Geyser in Yellowstone National Park, USA. We simultaneously measured water discharge, acoustic emissions, infraredintensity, and visible and infrared video to quantify the energetics and dynamics of eruptions, occurring approximately every three hours. We define four phases in the eruption cycle: 1) a 28&thinsp;&plusmn;&thinsp;3 minute phase with liquid and steam fountaining, with maximum jet velocities of 16&ndash;28&thinsp;m s &minus;&thinsp;1 , steam mass fraction of less than &sim;&thinsp;0.01. Intermittently choked flow and flow oscillations with periods increasing from 20 to 40&thinsp;s are coincident with a decrease in jet velocity and an increase of steam fraction; 2) a 26&thinsp;&plusmn;&thinsp;8 minute post&ndash;eruption relaxation phase with no discharge from the vent, infrared (IR) and acoustic power oscillations gliding between 30 and 40&thinsp;s; 3) a 59&thinsp;&plusmn;&thinsp;13 minute recharge period during which the geyser is quiescent and progressively refills, and 4) a 69&thinsp;&plusmn;&thinsp;14 minute pre&ndash;play period characterized by a series of 5&ndash;10&thinsp;minute&ndash;long pulses of steam, small volumes of liquid water discharge and 50&ndash;70&thinsp;s flow oscillations. The erupted waters ascend froma 160&thinsp;&minus;&thinsp;170&deg; C reservoir and the volume discharged during the entire eruptive cycle is 20.8&thinsp;&plusmn;&thinsp;4.1 m 3 . Assuming isentropic expansion, we calculate a heat output from the geyser of 1.4&ndash;1.5&thinsp;MW, which is <&thinsp;0.1% of the total heat output from Yellowstone Caldera.

Wyoming

Seismicity around Parkfield correlates with static shear stress changes following the 2003 M w 6.5 San Simeon earthquake

Earthquakes trigger other earthquakes, but the physical mechanism of the triggering is currently debated. Most studies of earthquake triggering rely on earthquakes listed in catalogs, which are known to be incomplete around the origin times of large earthquakes and therefore missing potentially triggered events. Here we apply a waveform matched-filter technique to systematically detect earthquakes along the Parkfield section of the San Andreas Fault from 46 days before to 31 days after the nearby 2003 M w 6.5 San Simeon earthquake. After removing all possible false detections, we identify ~8 times more earthquakes than in the Northern California Seismic Network catalog. The newly identified events along the creeping section of the San Andreas Fault show a statistically significant decrease following the San Simeon main shock, which correlates well with the negative static stress changes (i.e., stress shadow) cast by the main shock. In comparison, the seismicity rate around Parkfield increased moderately where the static stress changes are positive. The seismicity rate changes correlate well with the static shear stress changes induced by the San Simeon main shock, suggesting a low friction in the seismogenic zone along the Parkfield section of the San Andreas Fault.

California

Clustering of velocities in a GPS network spanning the Sierra Nevada Block, the northern Walker Lane Belt, and the Central Nevada Seismic Belt, California-Nevada

The deformation across the Sierra Nevada Block, the Walker Lane Belt, and the Central Nevada Seismic Belt (CNSB) between 38.5°N and 40.5°N has been analyzed by clustering GPS velocities to identify coherent blocks. Cluster analysis determines the number of clusters required and assigns the GPS stations to the proper clusters. The clusters are shown on a fault map by symbols located at the positions of the GPS stations, each symbol representing the cluster to which the velocity of that GPS station belongs. Fault systems that separate the clusters are readily identified on such a map. Four significant clusters are identified. Those clusters are strips separated by (from west to east) the Mohawk Valley-Genoa fault system, the Pyramid Lake-Wassuk fault system, and the Central Nevada Seismic Belt. The strain rates within the westernmost three clusters approximate simple right-lateral shear (~13 nstrain/a) across vertical planes roughly parallel to the cluster boundaries. Clustering does not recognize the longitudinal segmentation of the Walker Lane Belt into domains dominated by either northwesterly trending, right-lateral faults or northeasterly trending, left-lateral faults.

California;Nevada

Semiautomated tremor detection using a combined cross-correlation and neural network approach

Despite observations of tectonic tremor in many locations around the globe, the emergent phase arrivals, low‒amplitude waveforms, and variable event durations make automatic detection a nontrivial task. In this study, we employ a new method to identify tremor in large data sets using a semiautomated technique. The method first reduces the data volume with an envelope cross‒correlation technique, followed by a Self‒Organizing Map (SOM) algorithm to identify and classify event types. The method detects tremor in an automated fashion after calibrating for a specific data set, hence we refer to it as being “semiautomated”. We apply the semiautomated detection algorithm to a newly acquired data set of waveforms from a temporary deployment of 13 seismometers near Cholame, California, from May 2010 to July 2011. We manually identify tremor events in a 3 week long test data set and compare to the SOM output and find a detection accuracy of 79.5%. Detection accuracy improves with increasing signal‒to‒noise ratios and number of available stations. We find detection completeness of 96% for tremor events with signal‒to‒noise ratios above 3 and optimal results when data from at least 10 stations are available. We compare the SOM algorithm to the envelope correlation method of Wech and Creager and find the SOM performs significantly better, at least for the data set examined here. Using the SOM algorithm, we detect 2606 tremor events with a cumulative signal duration of nearly 55 h during the 13 month deployment. Overall, the SOM algorithm is shown to be a flexible new method that utilizes characteristics of the waveforms to identify tremor from noise or other seismic signals.

California

A fluid-driven earthquake swarm on the margin of the Yellowstone caldera

Over the past several decades, the Yellowstone caldera has experienced frequent earthquake swarms and repeated cycles of uplift and subsidence, reflecting dynamic volcanic and tectonic processes. Here, we examine the detailed spatial-temporal evolution of the 2010 Madison Plateau swarm, which occurred near the northwest boundary of the Yellowstone caldera. To fully explore the evolution of the swarm, we integrated procedures for seismic waveform-based earthquake detection with precise double-difference relative relocation. Using cross-correlation of continuous seismic data and waveform templates constructed from cataloged events, we detected and precisely located 8710 earthquakes during the three-week swarm, nearly four times the number of events included in the standard catalog. This high-resolution analysis reveals distinct migration of earthquake activity over the course of the swarm. The swarm initiated abruptly on January 17, 2010 at about 10 km depth and expanded dramatically outward (both shallower and deeper) over time, primarily along a NNW-striking, ~55º ENE-dipping structure. To explain these characteristics, we hypothesize that the swarm was triggered by the rupture of a zone of confined high-pressure aqueous fluids into a pre-existing crustal fault system, prompting release of accumulated stress. The high-pressure fluid injection may have been accommodated by hybrid shear and dilatational failure, as is commonly observed in exhumed hydrothermally affected fault zones. This process has likely occurred repeatedly in Yellowstone as aqueous fluids exsolved from magma migrate into the brittle crust, and it may be a key element in the observed cycles of caldera uplift and subsidence.

Wyoming

Very long period conduit oscillations induced by rockfalls at Kilauea Volcano, Hawaii

Eruptive activity at the summit of Kilauea Volcano, Hawaii, beginning in 2010 and continuing to the present time is characterized by transient outgassing bursts accompanied by very long period (VLP) seismic signals triggered by rockfalls from the vent walls impacting a lava lake in a pit within the Halemaumau pit crater. We use raw data recorded with an 11-station broadband network to model the source mechanism of signals accompanying two large rockfalls on 29 August 2012 and two smaller average rockfalls obtained by stacking over all events with similar waveforms to improve the signal-to-noise ratio. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. We apply a new waveform inversion method that accounts for the contributions from both translation and tilt in horizontal seismograms through the use of Green's functions representing the seismometer response to translation and tilt ground motions. This method enables a robust description of the source mechanism over the period range 1–1000 s. The VLP signals associated with the rockfalls originate in a source region ∼1 km below the eastern perimeter of the Halemaumau pit crater. The observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks including an east striking crack (dike) dipping 80° to the north, intersecting a north striking crack (another dike) dipping 65° to the east. Each rockfall is marked by a similar step-like inflation trailed by decaying oscillations of the volumetric source, attributed to the efficient coupling at the source centroid location of the pressure and momentum changes induced by the rock mass impacting the top of the lava column. Assuming a simple lumped parameter representation of the shallow magmatic system, the observed pressure and volume variations can be modeled with the following attributes: rockfall volume (200–4500 m 3 ), length of magma column (120–210 m), diameter of pipe connecting the Halemaumau pit crater to the subjacent dike system (6 m), average thickness of the two underlying dikes (3–6 m), and effective magma viscosity (30–210 Pa s). Most rockfalls occur during episodes of sustained deflation of the Kilauea summit. The mass loss rate in the shallow magmatic system is estimated to be 1400–15,000 kg s −1 based on measurements of the temporal variation of VLP period in the two large rockfalls that occurred on 29 August 2012.

Hawaii

Abrupt along-strike change in tectonic style: San Andreas fault zone, San Francisco Peninsula

Seismicity and high-resolution aeromagnetic data are used to define an abrupt change from compressional to extensional tectonism within a 10- to 15-km-wide zone along the San Andreas fault on the San Francisco Peninsula and offshore from the Golden Gate. This 100-km-long section of the San Andreas fault includes the hypocenter of the M w = 7.8 1906 San Francisco earthquake as well as the highest level of persistent microseismicity along that ∼470-km-long rupture. We define two distinct zones of deformation along this stretch of the fault using well-constrained relocations of all post-1969 earthquakes based a joint one-dimensional velocity/hypocenter inversion and a redetermination of focal mechanisms. The southern zone is characterized by thrust- and reverse-faulting focal mechanisms with NE trending P axes that indicate “fault-normal” compression in 7- to 10-km-wide zones of deformation on both sides of the San Andreas fault. A 1- to 2-km-wide vertical zone beneath the surface trace of the San Andreas is characterized by its almost complete lack of seismicity. The compressional deformation is consistent with the young, high topography of the Santa Cruz Mountains/Coast Ranges as the San Andreas fault makes a broad restraining left bend (∼10°) through the southernmost peninsula. A zone of seismic quiescence ∼15 km long separates this compressional zone to the south from a zone of combined normal-faulting and strike-slip-faulting focal mechanisms (including a M L = 5.3 earthquake in 1957) on the northernmost peninsula and offshore on the Golden Gate platform. Both linear pseudogravity gradients, calculated from the aeromagnetic data, and seismic reflection data indicate that the San Andreas fault makes an abrupt ∼3-km right step less than 5 km offshore in this northern zone. A similar right-stepping (dilatational) geometry is also observed for the subparallel San Gregorio fault offshore. Persistent seismicity and extensional tectonism occur within the San Andreas right stepover region and at least 15 km along-strike both to the SE and NW. The 1906 San Francisco earthquake may have nucleated within the San Andreas right stepover, which may help explain the bilateral nature of rupture of this event. Our analysis suggests two seismic hazards for the San Francisco Peninsula in addition to the hazard associated with a M = 7 to 8 strike-slip earthquake along the San Andreas fault: the potential for a M ≈ 6 normal-faulting earthquake just 5–8 km west of San Francisco and a M = 6+ thrust faulting event in the southern peninsula.

Journal of Geophysical Research B: Solid Earth

Analysis of stress-induced oval fractures in a borehole at Deep Sea Drilling Project Site 504, eastern equatorial Pacific

Deep Sea Drilling Project (DSDP) Hole 504B is located in the eastern equatorial Pacific Ocean and extends to a total depth of 2111 m beneath the seafloor (mbsf). Several acoustic televiewer logs have been obtained in this well during successive stages of drilling, and the resulting digital images have revealed numerous oval-shaped fractures seemingly etched into the borehole wall. A theoretical examination of these stress-induced features identifies a unique and ephemeral set of stress distributions and magnitudes that are necessary for their production. Consequently, the ovals provide a basis for quantifying the magnitudes and orientations of the maximum and minimum horizontal principal stresses, S H and S h at this site. Vertical, truncated breakouts and horizontal tensile fractures define the spatial boundaries of the ovals. Explicit criteria for their occurrence are combined with estimates for various physical properties of the rock to yield a range of possible values for the horizontal principal stresses. The conspicuous oval geometry is completed by a curved fracture that joins the vertical and horizontal components. Its degree of curvature is delineated by the modified Griffith failure criterion and is directly related to the principal stress difference ( S H − S h ). Matching a series of type curves corresponding to specific values for ( S H − S h ) with the actual undistorted well bore images allows the magnitude of the stress difference to be further constrained. With a value for ( S H − S h ) of 45 ± 5 MPa the individual magnitudes of S H and S h are determined more precisely. Final estimates for the horizontal principal stresses in DSDP Hole 504B at a depth of 1200 mbsf are 141 MPa ≤ S H ≤ 149 MPa and 91 MPa ≤ S h ≤ 109 MPa. Stress magnitudes derived from this approach rely heavily upon the values of a variety of physical properties, and complementary laboratory measurements performed on relevant rock samples provide critical information. Uncertainties in estimating these properties translate into less precise predictions of principal stresses.

Journal of Geophysical Research B: Solid Earth

P wave crustal velocity structure in the greater Mount Rainier area from local earthquake tomography

We present results from a local earthquake tomographic imaging experiment in the greater Mount Rainier area. We inverted P wave arrival times from local earthquakes recorded at permanent and temporary Pacific Northwest Seismograph Network seismographs between 1980 and 1996. We used a method similar to that described by Lees and Crosson [1989], modified to incorporate the parameter separation method for decoupling the hypocenter and velocity problems. In the upper 7 km of the resulting model there is good correlation between velocity anomalies and surface geology. Many focal mechanisms within the St. Helens seismic zone have nodal planes parallel to the epicentral trend as well as to a north-south trending low-velocity trough, leading us to speculate that the trough represents a zone of structural weakness in which a moderate ( M 6.5–7.0) earthquake could occur. In contrast, the western Rainier seismic zone does not correlate in any simple way with anomaly patterns or focal mechanism fault planes, leading us to infer that it is less likely to experience a moderate earthquake. A ∼10 km-wide low-velocity anomaly occurs 5 to 18 km beneath the summit of Mount Rainier, which we interpret to be a signal of a region composed of hot, fractured rock with possible small amounts of melt or fluid. No systematic velocity pattern is observed in association with the southern Washington Cascades conductor. A midcrustal anomaly parallels the Olympic-Wallowa lineament as well as several other geophysical trends, indicating that it may play an important role in regional tectonics.

Journal of Geophysical Research B: Solid Earth

Stress sensitivity of fault seismicity: A comparison between limited-offset oblique and major strike-slip faults

We present a new three-dimensional inventory of the southern San Francisco Bay area faults and use it to calculate stress applied principally by the 1989 M = 7.1 Loma Prieta earthquake and to compare fault seismicity rates before and after 1989. The major high-angle right-lateral faults exhibit a different response to the stress change than do minor oblique (right-lateral/thrust) faults. Seismicity on oblique-slip faults in the southern Santa Clara Valley thrust belt increased where the faults were undamped. The strong dependence of seismicity change on normal stress change implies a high coefficient of static friction. In contrast, we observe that faults with significant offset (>50–100 km) behave differently; microseismicity on the Hayward fault diminished where right-lateral shear stress was reduced and where it was undamped by the Loma Prieta earthquake. We observe a similar response on the San Andreas fault zone in southern California after the Landers earthquake sequence. Additionally, the offshore San Gregorio fault shows a seismicity rate increase where right-lateral/oblique shear stress was increased by the Loma Prieta earthquake despite also being clamped by it. These responses are consistent with either a low coefficient of static friction or high pore fluid pressures within the fault zones. We can explain the different behavior of the two styles of faults if those with large cumulative offset become impermeable through gouge buildup; coseismically pressurized pore fluids could be trapped and negate imposed normal stress changes, whereas in more limited offset faults, fluids could rapidly escape. The difference in behavior between minor and major faults may explain why frictional failure criteria that apply intermediate coefficients of static friction can be effective in describing the broad distributions of aftershocks that follow large earthquakes, since many of these events occur both inside and outside major fault zones.

Journal of Geophysical Research B: Solid Earth

Paleomagnetism of an east-west transect across the Cascade arc in southern Washington: Implications for regional tectonism

Paleomagnetic data from a transect across the Cascade arc in southern Washington were collected to test models of vertical axis rotations for this part of the Pacific Northwest and to provide information on the tectonic history of the St. Helens seismic zone. The 75 site-mean directions are divided into three groups based on isotopic, fission track, and stratigraphic ages. Group 1 consists of samples from rocks deposited between 39 and 30 Ma and shows a mean clockwise vertical axis rotation of 34° ± 13°, group 2 consists of samples from rocks deposited between 30 and 24 Ma and shows a mean clockwise rotation of 17° ± 11°, and group 3 consists of samples from rocks between 24 and 16 Ma and shows a mean clockwise rotation of 20° ± 12°. Although these three values of rotation are statistically indistinguishable at the 95% confidence level, we interpret them in combination with other data to indicate differential rotation across the St. Helens seismic zone (SHZ). The available paleomagnetic data for Eocene-Oligocene rocks west of the seismic zone show clockwise vertical axis rotations of 30° ± 8° and 35° ± 9° consistent with the value for group 1. The rotational values for groups 2 and 3 and the intrusive suite of Kidd Creek (13 Ma), despite their different ages, have similar values east of the SHZ. Comparing these groups of data east and west of the seismic zone indicates a differential rotation of 10° ± 3° across it. In addition, the rates of rotation in southern Washington are similar to those for rocks in the western Cascades of Oregon and indicate that rotational deformation of Miocene age extends northward into southern Washington and eastward into the axis of the arc.

Journal of Geophysical Research B: Solid Earth

Stress transfer by the 1988-1989 M=5.3 and 5.4 Lake Elsman foreshocks to the Loma Prieta fault: Unclamping at the site of peak mainshock slip

We study the stress transferred by the June 27, 1988, M = 5.3 and August 8, 1989, M = 5.4 Lake Elsman earthquakes, the largest events to strike within 15 km of the future Loma Prieta rupture zone during 74 years before the 1989 M = 6.9 Loma Prieta earthquake. We find that the first Lake Elsman event brought the rupture plane of the second event 0.3–1.6 bars (0.03–0.16 MPa) closer to Coulomb failure but that the Lake Elsman events did not bring the future Loma Prieta hypocentral zone closer to failure. Instead, the Lake Elsman earthquakes are calculated to have reduced the normal stress on (or “undamped”) the Loma Prieta rupture surface by 0.5–1.0 bar (0.05–0.10 MPa) at the site where the greatest slip subsequently occurred in the Loma Prieta earthquake. This association between the sites of peak unclamping and slip suggests that the Lake Elsman events did indeed influence the Loma Prieta rupture process. Unclamping the fault would have locally lowered the resistance to sliding. Such an effect could have been enhanced if the lowered normal stress permitted fluid infusion into the undamped part of the fault. Although less well recorded, the M L = 5.0 1964 and M L = 5.3 1967 Corralitos events struck within 10 km of the southwest end of the future Loma Prieta rupture. No similar relationship between the normal stress change and subsequent Loma Prieta slip is observed, although the high-slip patch southwest of the Loma Prieta epicenter corresponds roughly to the site of calculated Coulomb stress increase for a low coefficient of friction. The Lake Elsman-Loma Prieta result is similar to that for the 1987 M = 6.2 Elmore Ranch and M = 6.7 Superstition Hills earthquakes, suggesting that foreshocks might influence the distribution of mainshock slip rather than the site of mainshock nucleation.

Journal of Geophysical Research B: Solid Earth

Three-dimensional velocity structure of Siletzia and other accreted terranes in the Cascadia forearc of Washington

Eocene mafic crust with high seismic velocities underlies much of the Oregon and Washington forearc and acts as a backstop for accretion of marine sedimentary rocks from the obliquely subducting Juan de Fuca slab. Arc-parallel migration of relatively strong blocks of this terrane, known as Siletzia, focuses upper crustal deformation along block boundaries, which are potential sources of earthquakes. In a three-dimensional velocity model of coastal Washington, we have combined surface geology, well data, and travel times from earthquakes and controlled source seismic experiments to resolve the major boundaries of the Siletz terrane with the adjacent accreted sedimentary prism and volcanic arc. In southern Washington and northern Oregon the Siletz terrane appears to be a thick block (???20 km) that extends west of the coastline and makes a high-angle contact with the offshore accreted sedimentary prism. On its east flank the high-velocity Siletz terrane boundary coincides with an en echelon zone of seismicity in the arc. In northern Washington the western edge of Siletzia makes a lower-angled, fault-bound contact with the accretionary prism. In addition, alternating, east-west trending uplifts and downwarps of the Siletz terrane centered on the antiformal Olympic Mountains may reflect focusing of north-south compression in the northern part of the Siletz terrane. This compressional strain may result from northward transport and clockwise rotation of the Siletz terrane into the relatively fixed Canadian Coast Mountains restraining bend along the coast.

Oregon, Washington

Small scatterers in the lower mantle observed at German broadband arrays

Seismograms of earthquakes from the South Pacific recorded at a German broadband array and network show precursors to PKPdf . These precursors mainly originate from off-path scattering of PKPab or a nearby PKPbc to P (for receiver-side scattering) or from scattering of P to PKPab or PKPbc on the PKPdf path (for source-side scattering). Standard array processing techniques based on plane wave approximations (such as vespagram or frequency-wavenumber analysis) are inadequate for investigating these precursors since scattered waves cannot be approximated as plane waves for arrays and networks larger than 300 × 300 km for short-period waves. We therefore develop a migration method to estimate the location of scatterers in the mantle, at the core-mantle boundary and at the top of the outer core. With our method we are able to find isolated scatterers at the source side and the receiver side, although the depth of the scatterer is not well constrained. However, from looking at the first possible arrival time of precursors at different depth and the region where scattering can take place (scattering volume), we believe that the location of the scatterers is in the lowermost mantle. Since we have detected scatterers in regions where ultralow-velocity zones have been discovered recently, we think that the precursor energy possibly originates from scattering at partial melt at the base of the mantle. Comparing results from broadband and band-pass-filtered data the detection of small-scale structure of the ultralow-velocity zones becomes possible.

Journal of Geophysical Research B: Solid Earth