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 577 records · Page 32Linked to original sources
USGS2005

Magmatic effects of the Cobb hot spot on the Juan de Fuca Ridge

The interaction of the Juan de Fuca Ridge with the Cobb hot spot has had a considerable influence on the magmatism of the Axial Segment of the ridge, the second-order segment that overlies the hot spot. In addition to the construction of the large volcanic edifice of Axial Seamount, the Axial Segment has shallow bathymetry and a prevalence of constructional volcanic features along its 100-km length, suggesting that hot spot-derived magmas supplement and oversupply the ridge. Lavas are generally more primitive at Axial Seamount and more evolved in the Axial Segment rift zones, suggesting that fractional crystallization is enhanced with increasing distance from the hot spot because of a reduced magma supply and more rapid cooling. Although the Cobb hot spot is not an isotopically enriched plume, it produces lavas with some distinct geochemical characteristics relative to normal mid-ocean ridge basalt, such as enrichments in alkalis and highly incompatible trace elements, that can be used as tracers to identify the presence and prevalence of the hot spot influence along the ridge. These characteristics are most prominent at Axial Seamount and decline in gradients along the Axial Segment. The physical model that can best explain the geochemical observations is a scenario in which hot spot and mid-ocean ridge basalt (MORB) magmas mix to varying degrees, with the proportions controlled by the depth to the MORB source. Modeling of two-component mixing suggests that MORB is the dominant component in most Axial Segment basalts. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Noise in two-color electronic distance meter measurements revisited

Frequent, high-precision geodetic data have temporally correlated errors. Temporal correlations directly affect both the estimate of rate and its standard error; the rate of deformation is a key product from geodetic measurements made in tectonically active areas. Various models of temporally correlated errors are developed and these provide relations between the power spectral density and the data covariance matrix. These relations are applied to two-color electronic distance meter (EDM) measurements made frequently in California over the past 15-20 years. Previous analysis indicated that these data have significant random walk error. Analysis using the noise models developed here indicates that the random walk model is valid for about 30% of the data. A second 30% of the data can be better modeled with power law noise with a spectral index between 1 and 2, while another 30% of the data can be modeled with a combination of band-pass-filtered plus random walk noise. The remaining 10% of the data can be best modeled as a combination of band-pass-filtered plus power law noise. This band-pass-filtered noise is a product of an annual cycle that leaks into adjacent frequency bands. For time spans of more than 1 year these more complex noise models indicate that the precision in rate estimates is better than that inferred by just the simpler, random walk model of noise.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Upper crustal structure of southwestern British Columbia from the 1998 Seismic Hazards Investigation in Puget Sound

This paper applies nonlinear three-dimensional travel time tomography to refraction data recorded during the 1998 Seismic Hazards Investigation in Puget Sound (SHIPS) to derive the first large-scale, high-resolution upper crustal velocity model for southwestern British Columbia. A minimum structure P wave velocity model is constructed using 175,000 first arrival travel times picked from data recorded by 58 temporary onshore stations. The model details forearc crustal structures related to terrane accretion and subsequent basin formation to a depth of about 10 km. The Metchosin igneous complex (correlative with the Eocene Crescent-Siletz Terrane in Washington) is imaged as a laterally extensive WNW trending high-velocity anomaly underlying southernmost Vancouver Island and much of the Strait of Juan du Fuca. Northeast of the Strait of Georgia, the southwesterly dip of the contact between the Wrangellia terrane rocks of Vancouver Island and the Coast Plutonic Complex suggests Wrangellia rocks are down-faulted against the plutonic complex. At the southwestern end of the Strait of Juan de Fuca, the 50 km long WNW trending Clallam basin has a maximum thickness of 5-6 km. Near the eastern end of the Strait of Juan de Fuca, Port Townsend basin has an inferred thickness of approximately 4-5 km. The southern end of the 9 km thick Georgia basin is bounded by a high-velocity basement ridge. Beneath the Strait of Georgia, clusters of well-located earthquakes have a prominent NW trend and coincide spatially with rapid lateral velocity changes. Clusters of microearthquakes there are associated with the intersection of several east trending structural highs within this NW trend. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

New constraints on the active tectonic deformation of the Aegean

Site velocities from six separate Global Positioning System (GPS) networks comprising 374 stations have been referred to a single common Eurasia-fixed reference frame to map the velocity distribution over the entire Aegean. We use the GPS velocity field to identify deforming regions, rigid elements, and potential microplate boundaries, and build upon previous work by others to initially specify rigid elements in central Greece, the South Aegean, Anatolia, and the Sea of Marmara. We apply an iterative approach, tentatively defining microplate boundaries, determining best fit rigid rotations, examining misfit patterns, and revising the boundaries to achieve a better match between model and data. Short-term seismic cycle effects are minor contaminants of the data that we remove when necessary to isolate the long-term kinematics. We find that present day Aegean deformation is due to the relative motions of four microplates and straining in several isolated zones internal to them. The RMS misfit of model to data is about 2-sigma, very good when compared to the typical match between coseismic fault models and GPS data. The simplicity of the microplate description of the deformation and its good fit to the GPS data are surprising and were not anticipated by previous work, which had suggested either many rigid elements or broad deforming zones that comprise much of the Aegean region. The isolated deforming zones are also unexpected and cannot be explained by the kinematics of the microplate motions. Strain rates within internally deforming zones are extensional and range from 30 to 50 nanostrain/year (nstrain/year, 10-9/year), 1 to 2 orders of magnitude lower than rates observed across the major microplate boundaries. Lower strain rates may exist elsewhere withi the microplates but are only resolved in Anatolia, where extension of 13 ?? 4 nstrain/ year is required by the data. Our results suggest that despite the detailed complexity of active continental deformation revealed by seismicity, active faulting, fault geomorphology, and earthquake fault plane solutions, continental tectonics, at least in the Aegean, is to first order very similar to global plate tectonics and obeys the same simple kinematic rules. Although the widespread distribution of Aegean seismicity and active faulting might suggest a rather spatially homogeneous seismic hazard, the focusing of deformation near microplate boundaries implies the highest hazard is comparably localized.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Crustal deformation measurements in Guerrero, Mexico

GPS measurements of crustal deformation in Guerrero, southern Mexico, include surveys collected between 1992 and 2001 as well as continuous GPS measurements at a few sites. These geodetic observations are used to calculate interseismic deformation rates and assess the presence and possible location of transient deformation during the period encompassing 1992.25 to 2001.75. The data are used to examine transient deformation in 1998 previously described from data at a single site by Lowry et al. [2001]. Survey measurements and continuous data from a site near Popocate??petl volcano confirm the 1998 transient, and survey data also suggest another transient occurred following the 14 September 1995 (Mw = 7.3) Copala earthquake. All of the available GPS position estimates have been inverted for a combined model of slip during each event plus the steady state slip on the plate interface. Modeling of the steady state deformation rates confirms that the Guerrero seismic gap is partially frictionally locked at depths shallower than about 25 km and accumulating strain that may eventually be released in a great earthquake. The data also suggest that there is frictional coupling to much greater (>40 km) depths, which releases more frequently in aseismic slip events. The locations and sizes of the transient events are only partially constrained by the available data. However, the transient models which best fit the GPS coordinate time series suggest that aseismic slip was centered downdip of the seismogenic portion of the plate-bounding thrust in both events, and the moment release had equivalent magnitudes Mw = 7.1 + 1.3/-1.0 in 1995-1996 and 7.1 + 0.4/-0.1 in 1998. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

New insights on stress rotations from a forward regional model of the San Andreas fault system near its Big Bend in southern California

Understanding the stress field surrounding and driving active fault systems is an important component of mechanistic seismic hazard assessment. We develop and present results from a time-forward three-dimensional (3-D) model of the San Andreas fault system near its Big Bend in southern California. The model boundary conditions are assessed by comparing model and observed tectonic regimes. The model of earthquake generation along two fault segments is used to target measurable properties (e.g., stress orientations, heat flow) that may allow inferences on the stress state on the faults. It is a quasi-static model, where GPS-constrained tectonic loading drives faults modeled as mostly sealed viscoelastic bodies embedded in an elastic half-space subjected to compaction and shear creep. A transpressive tectonic regime develops southwest of the model bend as a result of the tectonic loading and migrates toward the bend because of fault slip. The strength of the model faults is assessed on the basis of stress orientations, stress drop, and overpressures, showing a departure in the behavior of 3-D finite faults compared to models of 1-D or homogeneous infinite faults. At a smaller scale, stress transfers from fault slip transiently induce significant perturbations in the local stress tensors (where the slip profile is very heterogeneous). These stress rotations disappear when subsequent model earthquakes smooth the slip profile. Maps of maximum absolute shear stress emphasize both that (1) future models should include a more continuous representation of the faults and (2) that hydrostatically pressured intact rock is very difficult to break when no material weakness is considered. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Precise relocation of earthquakes following the 15 June 1991 eruption of Mount Pinatubo (Philippines)

The 15 June 1991 climactic eruption of Mount Pinatubo (Philippines) was followed by intense seismicity that remained at a high level for several months. We located 10,839 events recorded between 1 July and mid-December 1991. In contrast to the preeruptive seismicity which was focused in two groups below the summit area, posteruptive events were widely distributed below and around the volcano. The classification of the events indicates the presence of several large multiplets, and the application of relative relocation techniques to the similar events by calculating high-precision delays between traces outlines a number of clear seismogenic structures. We used different methods to confirm the validity of our results; these tests indicate that reliable features can be detected with a small monitoring network. While the main cluster of activity can be attributed to an intrusive process starting from below the 15 June crater, the volcanic origin of the seismic activity in the other areas is more difficult to establish. Away from the summit, relocations define streaks or planes which are oriented predominantly southwest-northeast, with in several cases the presence of northwest-southeast conjugate structures. Most of the composite focal mechanisms that we could determine indicate predominantly strike-slip, right-lateral faulting. Our results indicate that most of the seismicity that occurred after the 15 June eruption is related to the east-west regional compressional stress field related to the subduction. We suggest that the regional stress field induces seismicity along new or preexisting faults in the medium surrounding the volcano where the stress field was locally disturbed by the volcanic eruption.

Mount Pinatubo↗
USGS2004

A physical model for strain accumulation in the San Francisco Bay region: Stress evolution since 1838

Understanding of the behavior of plate boundary zones has progressed to the point where reasonably comprehensive physical models can predict their evolution. The San Andreas fault system in the San Francisco Bay region (SFBR) is dominated by a few major faults whose behavior over about one earthquake cycle is fairly well understood. By combining the past history of large ruptures on SFBR faults with a recently proposed physical model of strain accumulation in the SFBR, we derive the evolution of regional stress from 1838 until the present. This effort depends on (1) an existing compilation of the source properties of historic and contemporary SFBR earthquakes based on documented shaking, geodetic data, and seismic data (Bakun, 1999) and (2) a few key parameters of a simple regional viscoelastic coupling model constrained by recent GPS data (Pollitz and Nyst, 2004). Although uncertainties abound in the location, magnitude, and fault geometries of historic ruptures and the physical model relies on gross simplifications, the resulting stress evolution model is sufficiently detailed to provide a useful window into the past stress history. In the framework of Coulomb failure stress, we find that virtually all M ??? 5.8 earthquakes prior to 1906 and M ??? 5.5 earthquakes after 1906 are consistent with stress triggering from previous earthquakes. These events systematically lie in zones of predicted stress concentration elevated 5-10 bars above the regional average. The SFBR is predicted to have emerged from the 1906 "shadow" in about 1980, consistent with the acceleration in regional seismicity at that time. The stress evolution model may be a reliable indicator of the most likely areas to experience M ??? 5.5 shocks in the future.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Stress interaction between subduction earthquakes and forearc strike-slip faults: Modeling and application to the northern Caribbean plate boundary

Strike-slip faults in the forearc region of a subduction zone often present significant seismic hazard because of their proximity to population centers. We explore the interaction between thrust events on the subduction interface and strike-slip faults within the forearc region using three-dimensional models of static Coulomb stress change. Model results reveal that subduction earthquakes with slip vectors subparallel to the trench axis enhance the Coulomb stress on strike-slip faults adjacent to the trench but reduce the stress on faults farther back in the forearc region. In contrast, subduction events with slip vectors perpendicular to the trench axis enhance the Coulomb stress on strike-slip faults farther back in the forearc, while reducing the stress adjacent to the trench. A significant contribution to Coulomb stress increase on strike-slip faults in the back region of the forearc comes from "unclamping" of the fault, i.e., reduction in normal stress due to thrust motion on the subduction interface. We argue that although Coulomb stress changes from individual subduction earthquakes are ephemeral, their cumulative effects on the pattern of lithosphere deformation in the forearc region are significant. We use the Coulomb stress models to explain the contrasting deformation pattern between two adjacent segments of the Caribbean subduction zone. Subduction earthquakes with slip vectors nearly perpendicular to the Caribbean trench axis is dominant in the Hispaniola segment, where the strike-slip faults are more than 60 km inland from the trench. In contrast, subduction slip motion is nearly parallel to the Caribbean trench axis along the Puerto Rico segment, where the strike-slip fault is less than 15 km from the trench. This observed jump from a strike-slip fault close to the trench axis in the Puerto Rico segment to the inland faults in Hispaniola is explained by different distributions of Coulomb stress in the forearc region of the two segments, as a result of the change from the nearly trench parallel slip on the Puerto Rico subduction interface to the more perpendicular subduction slip beneath Hispaniola. The observations and modeling suggest that subduction-induced strike-slip seismic hazard to Puerto Rico may be smaller than previously assumed but the hazard to Hispaniola remains high. Copyright 2004 by the American Geophysical Union.

Puerto Rico↗
USGS2005

Interferometric synthetic aperture radar study of Okmok volcano, Alaska, 1992-2003: Magma supply dynamics and postemplacement lava flow deformation

Okmok volcano, located in the central Aleutian arc, Alaska, is a dominantly basaltic complex topped with a 10-km-wide caldera that formed circa 2.05 ka. Okmok erupted several times during the 20th century, most recently in 1997; eruptions in 1945, 1958, and 1997 produced lava flows within the caldera. We used 80 interferometric synthetic aperture radar (InSAR) images (interferograms) to study transient deformation of the volcano before, during, and after the 1997 eruption. Point source models suggest that a magma reservoir at a depth of 3.2 km below sea level, located beneath the center of the caldera and about 5 km northeast of the 1997 vent, is responsible for observed volcano-wide deformation. The preeruption uplift rate decreased from about 10 cm yr −1 during 1992–1993 to 2 ∼ 3 cm yr −1 during 1993–1995 and then to about −1 ∼ −2 cm yr −1 during 1995–1996. The posteruption inflation rate generally decreased with time during 1997–2001, but increased significantly during 2001–2003. By the summer of 2003, 30 ∼ 60% of the magma volume lost from the reservoir in the 1997 eruption had been replenished. Interferograms for periods before the 1997 eruption indicate consistent subsidence of the surface of the 1958 lava flows, most likely due to thermal contraction. Interferograms for periods after the eruption suggest at least four distinct deformation processes: (1) volcano-wide inflation due to replenishment of the shallow magma reservoir, (2) subsidence of the 1997 lava flows, most likely due to thermal contraction, (3) deformation of the 1958 lava flows due to loading by the 1997 flows, and (4) continuing subsidence of 1958 lava flows buried beneath 1997 flows. Our results provide insights into the postemplacement behavior of lava flows and have cautionary implications for the interpretation of inflation patterns at active volcanoes.

Alaska↗
USGS2004

Quantifying probabilities of volcanic events: The example of volcanic hazard at Mount Vesuvius

We describe an event tree scheme to quantitatively estimate both long- and short-term volcanic hazard. The procedure is based on a Bayesian approach that produces a probability estimation of any possible event in which we are interested and can make use of all available information including theoretical models, historical and geological data, and monitoring observations. The main steps in the procedure are (1) to estimate an a priori probability distribution based upon theoretical knowledge, (2) to modify that using past data, and (3) to modify it further using current monitoring data. The scheme allows epistemic and aleatoric uncertainties to be dealt with in a formal way, through estimation of probability distributions at each node of the event tree. We then describe an application of the method to the case of Mount Vesuvius. Although the primary intent of the example is to illustrate the methodology, one result of this application merits special mention. The present emergency response plan for Mount Vesuvius is referenced to a maximum expected event (MEE), the largest out of all the possible eruptions within the next few decades. Our calculation suggest that there is a nonnegligible (1-20%) chance that the next eruption could be larger than that stipulated in the present MEE. The methodology allows all assumptions and thresholds to be clearly identified and provides a rational means for their revision if new data or information are obtained. Copyright 2004 by the American Geophysical Union.

Mount Vesuvius↗
USGS2005

Significance of stress transfer in time-dependent earthquake probability calculations

A sudden change in stress is seen to modify earthquake rates, but should it also revise earthquake probability? Data used to derive input parameters permits an array of forecasts; so how large a static stress change is require to cause a statistically significant earthquake probability change? To answer that question, effects of parameter and philosophical choices are examined through all phases of sample calculations, Drawing at random from distributions of recurrence-aperiodicity pairs identifies many that recreate long paleoseismic and historic earthquake catalogs. Probability density funtions built from the recurrence-aperiodicity pairs give the range of possible earthquake forecasts under a point process renewal model. Consequences of choices made in stress transfer calculations, such as different slip models, fault rake, dip, and friction are, tracked. For interactions among large faults, calculated peak stress changes may be localized, with most of the receiving fault area changed less than the mean. Thus, to avoid overstating probability change on segments, stress change values should be drawn from a distribution reflecting the spatial pattern rather than using the segment mean. Disparity resulting from interaction probability methodology is also examined. For a fault with a well-understood earthquake history, a minimum stress change to stressing rate ratio of 10:1 to 20:1 is required to significantly skew probabilities with >80-85% confidence. That ratio must be closer to 50:1 to exceed 90-95% confidence levels. Thus revision to earthquake probability is achievable when a perturbing event is very close to the fault in question or the tectonic stressing rate is low.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Rupture dynamics with energy loss outside the slip zone

Energy loss in a fault damage zone, outside the slip zone, contributes to the fracture energy that determines rupture velocity of an earthquake. A nonelastic two-dimensional dynamic calculation is done in which the slip zone is modeled as a fault plane and material off the fault is subject to a Coulomb yield condition. In a mode 2 crack-like solution in which an abrupt uniform drop of shear traction on the fault spreads from a point, Coulomb yielding occurs on the extensional side of the fault. Plastic strain is distributed with uniform magnitude along the fault, and it has a thickness normal to the fault proportional to propagation distance. Energy loss off the fault is also proportional to propagation distance, and it can become much larger than energy loss on the fault specified by the fault constitutive relation. The slip velocity function could be produced in an equivalent elastic problem by a slip-weakening friction law with breakdown slip Dc increasing with distance. Fracture energy G and equivalent Dc will be different in ruptures with different initiation points and stress drops, so they are not constitutive properties; they are determined by the dynamic solution that arrives at a particular point. Peak slip velocity is, however, a property of a fault location. Nonelastic response can be mimicked by imposing a limit on slip velocity on a fault in an elastic medium.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

New constraints on mechanisms of remotely triggered seismicity at Long Valley Caldera

Regional-scale triggering of local earthquakes in the crust by seismic waves from distant main shocks has now been robustly documented for over a decade. Some of the most thoroughly recorded examples of repeated triggering of a single site from multiple, large earthquakes are measured in geothermal fields of the western United States like Long Valley Caldera. As one of the few natural cases where the causality of an earthquake sequence is apparent, triggering provides fundamental constraints on the failure processes in earthquakes. We show here that the observed triggering by seismic waves is inconsistent with any mechanism that depends on cumulative shaking as measured by integrated energy density. We also present evidence for a frequency-dependent triggering threshold. On the basis of the seismic records of 12 regional and teleseismic events recorded at Long Valley Caldera, long-period waves (>30 s) are more effective at generating local seismicity than short-period waves of comparable amplitude. If the properties of the system are stationary over time, the failure threshold for long-period waves is ~0.05 cm/s vertical shaking. Assuming a phase velocity of 3.5 km/s and an elastic modulus of 3.5 x 10 10 Pa, the threshold in terms of stress is 5 kPa. The frequency dependence is due in part to the attenuation of the surface waves with depth. Fluid flow through a porous medium can produce the rest of the observed frequency dependence of the threshold. If the threshold is not stationary with time, pore pressures that are >99.5% of lithostatic and vary over time by a factor of 4 could explain the observations with no frequency dependence of the triggering threshold.

California↗
USGS2005

Evolution of melt-vapor surface tension in silicic volcanic systems: Experiments with hydrous melts

We evaluate the melt‐vapor surface tension (σ) of natural, water‐saturated dacite melt at 200 MPa, 950–1055°C, and 4.8–5.7 wt % H 2 O. We experimentally determine the critical supersaturation pressure for bubble nucleation as a function of dissolved water and then solve for σ at those conditions using classical nucleation theory. The solutions obtained give dacite melt‐vapor surface tensions that vary inversely with dissolved water from 0.042 (±0.003) J m −2 at 5.7 wt % H 2 O to 0.060 (±0.007) J m −2 at 5.2 wt % H 2 O to 0.073 (±0.003) J m −2 at 4.8 wt % H 2 O. Combining our dacite results with data from published hydrous haplogranite and high‐silica rhyolite experiments reveals that melt‐vapor surface tension also varies inversely with the concentration of mafic melt components (e.g., CaO, FeO total , MgO). We develop a thermodynamic context for these observations in which melt‐vapor surface tension is represented by a balance of work terms controlled by melt structure. Overall, our results suggest that cooling, crystallization, and vapor exsolution cause systematic changes in σ that should be considered in dynamic modeling of magmatic processes.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Stress orientations at intermediate angles to the San Andreas Fault, California

There are currently two competing models for the frictional strength of the San Andreas Fault in California: the strong-fault model and the weak-fault model. The strong-fault model predicts the maximum horizontal compressive stress axis to be at low angles to the fault, while the relatively weak fault model predicts it to be at high angles. Previous studies have disagreed as to which model is supported by observed stress orientations. We review and compare these studies and present results from several new focal mechanism stress inversions. We find that the observed stress orientations of different studies are generally consistent, implying that the disagreement is one of interpretation. The majority of studies find compressive stress orientations at intermediate angles to the fault, not strictly consistent with either current model. The strong-fault model is acceptable if the San Andreas is assumed to be a nonoptimally orientated fault that fails because optimally oriented, preexisting planes are not present. The relatively weak fault model is not consistent with the stress orientations. We propose two alternative models to better explain the observed intermediate stress orientations: an intermediate-strength San Andreas model and a model in which all major active faults are weak.

California↗
USGS2005

Forearc structure beneath southwestern British Columbia: A three-dimensional tomographic velocity model

This paper presents a three-dimensional compressional wave velocity model of the forearc crust and upper mantle and the subducting Juan de Fuca plate beneath southwestern British Columbia and the adjoining straits of Georgia and Juan de Fuca. The velocity model was constructed through joint tomographic inversion of 50,000 first-arrival times from earthquakes and active seismic sources. Wrangellia rocks of the accreted Paleozoic and Mesozoic island arc assemblage underlying southern Vancouver Island in the Cascadia forearc are imaged at some locations with higher than average lower crustal velocities of 6.5-7.2 km/s, similar to observations at other island arc terranes. The mafic Eocene Crescent terrane, thrust landward beneath southern Vancouver Island, exhibits crustal velocities in the range of 6.0-6.7 km/s and is inferred to extend to a depth of more than 20 km. The Cenozoic Olympic Subduction Complex, an accretionary prism thrust beneath the Crescent terrane in the Olympic Peninsula, is imaged as a low-velocity wedge to depths of at least 20 km. Three zones with velocities of 7.0-7.5 km/s, inferred to be mafic and/or ultramafic units, lie above the subducting Juan de Fuca plate at depths of 25-35 km. The forearc upper mantle wedge beneath southeastern Vancouver Island and the Strait of Georgia exhibits low velocities of 7.2-7.5 km/s, inferred to correspond to ???20% serpentinization of mantle peridotites, and consistent with similar observations in other warm subduction zones. Estimated dip of the Juan de Fuca plate beneath southern Vancouver Island is ???11??, 16??, and 27?? at depths of 30, 40, and 50 km, respectively. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗