The LITHOPROBE seismic reflection project on Vancouver Island was designed to study the large-scale structure of several accreted terranes exposed on the island and to determine the geometry and structural characteristics of the subducting Juan de Fuca plate. In this paper, we interpret two LITHOPROBE profiles from southernmost Vancouver Island that were shot across three important terrane-bounding faults—Leech River, San Juan, and Survey Mountain—to determine their subsurface geometry and relationship to deeper structures associated with modem subduction.The structure beneath the island can be divided into an upper crustal region, consisting of several accreted terranes, and a deeper region that represents a landward extension of the modern offshore subduction complex. In the upper region, the Survey Mountain and Leech River faults are imaged as northeast-dipping thrusts that separate Wrangellia, a large Mesozoic–Paleozoic terrane, from two smaller accreted terranes: the Leech River schist, Mesozoic rocks that were metamorphosed in the Late Eocene; and the Metchosin Formation, a Lower Eocene basalt and gabbro unit. The Leech River fault, which was clearly imaged on both profiles, dips 35–45 °northeast and extends to about 10 km depth. The Survey Mountain fault lies parallel to and above the Leech River fault and extends to similar depths. The San Juan fault, the western continuation of the Survey Mountain fault, was not imaged, although indirect evidence suggests that it also is a thrust fault. These faults accommodated the Late Eocene amalgamation of the Leech River and Metchosin terranes along the southern perimeter of Wrangellia. Thereafter, these terranes acted as a relatively coherent lid for a younger subduction complex that has formed during the modem (40 Ma to present) convergent regime.Within this subduction complex, the LITHOPROBE profiles show three prominent bands of differing reflectivity that dip gently northeast. These bands represent regionally extensive layers lying beneath the lid of older accreted terranes. We interpret them as having formed by underplating of oceanic materials beneath the leading edge of an overriding continental place. The upper reflective layer can be projected updip to the south, where it is exposed in the Olympic Mountains as the Core rocks, an uplifted Cenozoic subduction complex composed dominantly of accreted marine sedimentary rocks. A middle zone of low reflectivity is not exposed at the surface, but results from an adjacent refraction survey indicate it is probably composed of relatively high velocity materials (~ 7.7 km/s). We consider two possibilities for the origin of this zone: (1) a detached slab of oceanic lithosphere accreted during an episodic tectonic event or (2) an imbricated package of mafic rocks derived by continuous accretion from the top of the subducting oceanic crust. The lower reflective layer is similar in reflection character to the upper layer and, therefore, is also interpreted as consisting dominantly of accreted marine sedimentary rocks. It represents the active zone of decoupling between the overriding and underthrusting plates and, thus, delimits present accretionary processes occurring directly above the descending Juan de Fuca plate. These results provide the first direct evidence for the process of subduction underplating or subcretion and illustrate a process that is probably important in the evolution and growth of continents.

Analysis of the Lithoprobe Deep Probe and Southern Alberta Refraction Experiment data sets, focusing on the region between Deep Probe shots 43 and 55, has resulted in a continental-scale velocity structural model of the lithosphere of platformal western Laurentia reaching depths of ~150 km. Three major lithospheric blocks were investigated: (i) the Hearne Province, a typical continental Archean cratonic province lying beneath the Western Canada Sedimentary Basin; (ii) the Wyoming Province, an even older block of Phanerozoic-modified Archean crust with an enigmatic lower lithosphere; and (iii) the Yavapai–Mazatzal Province, Proterozoic terranes underlying the Colorado Plateau and Southern Rocky Mountains. In this study, the northern two of these regions are investigated with a modified ray-theoretical traveltime inversion routine that respects the spherical geometry of the Earth. The resulting crustal velocity structure, combined with supporting geological and geophysical data, reveals that the Medicine Hat block (MHB), lying between the Hearne and Wyoming provinces, is a third independent Archean crustal block. The subcrustal lithosphere along the profile is homogeneous in velocity structure, but two significant northward-dipping reflectors are apparent and interpreted as relic subduction zones associated with sutures between the three Archean blocks. The Hearne crust is typical of an Archean shield or platform both in its thickness of 34–50 km and its seismic velocity structure. The crust of the Archean MHB and Wyoming Province, which ranges in thickness from 49 to 60 km, includes a 10–30 km thick high-velocity layer, interpreted to be Proterozoic in age. Such a feature is unexpected beneath Archean crustal provinces, but if the region is considered to be the remanent marginal portion of a larger Archean continent, then the interpreted Proterozoic underplating and lack of an Archean lithospheric root can be explained. The variable topography along the reflective upper and lower boundaries of this layer, especially within the MHB, suggests considerable variability in its emplacement and subsequent tectonic history.

A time domain polarization filter, originally proposed by Flinn, has been modified and is used to increase the signal to noise ratio of teleseismic body phases. Both the rectilinearity and direction of particle motion is obtained from the co-variance matrix for three components of ground motion over a small time interval. The estimate of rectilinearity is found by diagonalizing the matrix and forming a function that involves the ratio of the largest and intermediate principle axes of the matrix. The direction of polarization is found from the eigenvector of the largest principle axis. A set of time-varying operators are then obtained which act as a point by point gain control to modulate the digital seismic records.This non-linear filter is useful for enhancing P, p P , sP, PP, S and other compressional or shear phases. When applied to an array of three-component stations it appears to be possible to identify multiple events in the source function of some earthquakes.

Seismic reflection data from the south central Canadian Cordillera covering the interval from the easternmost metamorphic core complexes near Arrow Lakes to the Fraser River fault system along the Fraser River reveal a highly reflective and complex crust. The base of the crustal reflectivity, interpreted as the reflection Moho, is clearly delineated by a continuous sharp boundary that is essentially planar and slopes uniformly over a distance of 250 km from about 12.0 s in the east to about 10.5 s in the west. This virtual lack of relief at the base of the crust contrasts sharply with surface structures that involve 25 km or more of structural relief. Some of these surface structures can be readily correlated to structures that are outlined by the reflection data and that can be followed into the middle and lower crust. Even though part of this area was subjected to large amounts of Eocene extension, the crust is not divisible into transparent upper and reflective lower layers as it is in parts of the U.S. Cordillera. Three structural culminations, the Monashee complex, the Vernon antiform, and the Central Nicola horst, are interpreted on the basis of the reflection configuration and the surface geological relationships to have formed initially during Jurassic to Eocene compression and then to have been modified and exposed during early and middle Eocene extension. An example of a compressional structure observed on the profiles is the Monashee decollement, which can be traced from the surface westward into the lower crust. Extension is manifested along a variety of normal faults, including the regionally extensive low angle Okanagan Valley‐Eagle River fault system, moderately dipping faults such as the Columbia River and Slocan Lake faults, and high‐angle faults such as the Quilchena Creek and Coldwater faults. Both Jurassic to Eocene compressional shear zones and early to middle Eocene extensional shear zones are listric into the lower crust or Moho under the Intermontane belt.

A particle velocity‐stress, finite‐difference method is developed for the simulation of wave propagation in 2-D heterogeneous poroelastic media. Instead of the prevailing second‐order differential equations, we consider a first‐order hyperbolic system that is equivalent to Biot’s equations. The vector of unknowns in this system consists of the solid and fluid particle velocity components, the solid stress components, and the fluid pressure. A MacCormack finite‐difference scheme that is fourth‐order accurate in space and second‐order accurate in time forms the basis of the numerical solutions for Biot’s hyperbolic system. An original analytic solution for a P‐wave line source in a uniform poroelastic medium is derived for the purposes of source implementation and algorithm testing. In simulations with a two‐layer model, additional “slow” compressional incident, transmitted, and reflected phases are recorded when the damping coefficient is small. This “slow” compressional wave is highly attenuated in porous media saturated by a viscous fluid. From the simulation we also verified that the attenuation mechanism introduced in Biot’s theory is of secondary importance for “fast” compressional and rotational waves. The existence of seismically observable differences caused by the presence of pores has been examined through synthetic experiments that indicate that amplitude variation with offset may be observed on receivers and could be diagnostic of the matrix and fluid parameters. This method was applied in simulating seismic wave propagation over an expanded steam‐heated zone in Cold Lake, Alberta in an area of enhanced oil recovery (EOR) processing. The results indicate that a seismic surface survey can be used to monitor thermal fronts.

Lithoprobe's Southern Alberta Refraction Experiment, SAREX, extends 800 km from east-central Alberta to central Montana. It was designed to investigate crustal velocity structure of the Archean domains underlying the Western Canada Sedimentary Basin. From north to south, SAREX crosses the Loverna domain of the Hearne Province, the Vulcan structure, the Medicine Hat block (previously considered part of the Hearne Province), the Great Falls tectonic zone, and the northern Wyoming Province. Ten shot points along the profile in Canada were recorded on 521 seismographs deployed at 1 km intervals. To extend the line, an additional 140 seismographs were deployed at intervals of 1.25–2.50 km in Montana. Data interpretation used an iterative application of damped least-squares inversion of traveltime picks and forward modeling. Results show different velocity structures for the major blocks (Loverna, Medicine Hat, and Wyoming), indicating that each is distinct. Wavy undulations in the velocity structure of the Loverna block may be associated with internal crustal deformation. The most prominent feature of the model is a thick (10–25 km) lower crustal layer with high velocities (7.5–7.9 km/s) underlying the Medicine Hat and Wyoming blocks. Based on data from lower crustal xenoliths in the region, this layer is interpreted to be the result of Paleoproterozoic magmatic underplating. Crustal thickness varies from 40 km in the north to almost 60 km in the south, where the high-velocity layer is thickest. Uppermost mantle velocities range from 8.05 to 8.2 km/s, with the higher values below the thicker crust. Results from SAREX and other recent studies are synthesized to develop a schematic representation of Archean to Paleoproterozoic tectonic development for the region encompassing the profile. Tectonic processes associated with this development include collisions of continental blocks, subduction, crustal thickening, and magmatic underplating.

Variations in crustal thickness provide important clues as to the formation of the crust, present-day isostatic equilibrium, and crustal stress. The map of the topography of the Mohorovicic discontinuity in Alberta is revised using 1900 km of reanalysed seismic reflection profiles acquired as part of the Lithoprobe Alberta Basement Transect. Time sections were depth migrated using a parallelized algorithm that accounts for steeply dipping structures. The migration process employed a geologically consistent velocity model of the metamorphic crust derived from earlier refraction experiments and constrained by compilations of high-pressure rock velocity measurements. We found that knowledge of the velocities in the sedimentary column strongly influenced the quality of the migration calculations. The Mohorovičić discontinuity is generally distinguished in these profiles, on the basis of sharp changes in reflectivity, at depths of 35–48 km. Sharp reflections from this boundary are rare. A number of geologic features are of note in these lines. A localized (~50 km extent) crustal thinning is observed in the Peace River region; this thinning is consistent with an adjacent oxygen isotopic anomaly indicative of crustal extension. In central Alberta, the Mohorovičić discontinuity topography is suggestive of a sharp jump of 10 km indicative of mantle faulting associated with the Snowbird tectonic zone. In southern Alberta, the crust thickens substantially across the Vulcan structure, with the greatest thickness correlating with the Vulcan structure itself indicating a collisional origin as noted by other authors.

No abstract