Search

The Kisseynew domain (KD) is the largest component of the Paleoproterozoic Trans-Hudson orogen (THO) in Saskatchewan and Manitoba. It is bounded to the north by the Lynn Lake – Leaf Rapids (LL-LR) domain and to the south by the Glennie – Flin Flon (G-FF) complex. The THO was the focus of one of the study areas of Lithoprobe, Canada’s national Earth science research project (1984–2005). To further investigate the crustal structure of the KD and its bounding domains, this study reprocesses reflection line S3a across its northern boundary, analyses four 2.5-D gravity profiles, carries out 3-D gravity inversions for two areas, and replicates results from reflection lines 7 and 10 across the southern boundary of the KD. The reprocessed seismic section enhances the continuity of reflections within the crust. The reflectivity is representative of the complex tectonic development of the boundary zone and clearly identifies a subsurface deformation zone consistent with the boundary. The reflection section also shows that a lower plate (at about 30 to 50 km depth), interpreted as remnant lower crust of the G-FF complex, extends 30 km further northward than in the original section. Lines 7 and 10 illustrate the complex nature of the transition from the KD to the G-FF complex. The gravity analyses show that the variability and complexities of the boundary region between the LL-LR domain and KD, and the G-FF complex and KD, as indicated by the geological and Bouguer gravity maps, extend at depth throughout the crust.

EON-ROSE (Earth-System Observing Network - Réseau d’Observation du Système terrestrE) is a new initiative for a pan-Canadian research collaboration to holistically examine Earth systems from the ionosphere into the core. The Canadian Cordillera Array (CC Array) is the pilot phase, and will extend across the Cordillera from the Beaufort Sea to the U.S. border. The vision for EON-ROSE is to install a network of telemetered observatories to monitor solid Earth, environmental and atmospheric processes. EON-ROSE is an inclusive, combined effort of Canadian universities, federal, provincial and territorial government agencies, industry, and international collaborators. Brainstorming sessions and several workshops have been held since May 2016. The first station will be installed at Kluane Lake Research Station in southwestern Yukon during the summer of 2018. The purpose of this report is to provide a framework for continued discussion and development.RÉSUMÉEON-ROSE (Earth-System Observing Network - Réseau d’Observation du Système terrestrE) est une nouvelle initiative de collaboration de recherche pancanadienne visant à étudier de manière holistique les systèmes terrestres, depuis l’ionosphère jusqu’au noyau. Le Réseau canadien de la cordillère (CC Array) en est la phase pilote, laquelle couvrira toute la Cordillère, de la mer de Beaufort jusqu’à la frontière étasunienne. L’objectif d’EON-ROSE est d’installer un réseau d’observatoires télémétriques pour suivre en continu les processusterrestres, environnementaux et atmosphériques. EON-ROSE est un effort combiné et inclusif des universités canadiennes, des organismes gouvernementaux fédéraux, provinciaux et territoriaux, de l’industrie et de collaborateurs internationaux. Des séances de remue-méninges et plusieurs ateliers ont été tenus depuis mai 2016. La première station sera installée à la station de recherche du lac Kluane, dans le sud-ouest du Yukon, au cours de l’été 2018. Le but du présent rapport est de fournir un cadre de discussion et de développement continu.

Lithoprobe (1984–2005), Canada’s national, collaborative, multidisciplinary, Earth Science research project, investigated the structure and evolution of the Canadian landmass and its margins. It was a highly successful project that redefined the nature of Earth science research in Canada. One of many contributions deriving from the project was the demonstration by example that Earth scientists from geophysics and geology, including all applicable sub-disciplines within these general study areas, must work together to achieve thorough and comprehensive interpretations of all available data sets. In Part 1, this statement was exemplified through studies involving lithospheric structures. In Part 2, it is exemplified by summarizing interpretations from six exploration-related studies derived from journal publications. In the first example, subsurface structures associated with the Guichon Creek batholith in south-central British Columbia, which hosts porphyry copper and molybdenum deposits, are better defined and related to different geological phases of the batholith. Reprocessed seismic reflection data and 2.5-D and 3-D inversions of magnetic and gravity data are combined with detailed geological mapping and drillhole information to generate the revised and improved subsurface interpretation. Research around the Bell Allard volcanogenic massive sulphide deposit in the Matagami region of northern Quebec provides the second example. A seismic reflection line over the deposit shortly after it was discovered by drilling, aided by core and geophysical logs, was acquired to test whether the deposit could be imaged. Direct detection of the ore body from the seismic section would be difficult if its location were not already known; however, structural characteristics that can be tied to lithologies from boreholes and logs were well identified. Nickel deposits and associated structures in the Thompson belt at the western limit of the Superior Province in northern Manitoba were the focus of seismic and electromagnetic (EM) studies combined with geology and physical property measurements. The combined seismic/EM image indicates that the rocks of the prospective Ospwagan Group, which have low resistivity, extend southeastward beneath the Archean gneiss and that structural culminations control the subsurface geometry of the Ospwagan Group. The Sudbury structure in Ontario is famous for its nickel deposits, the largest in the world, which formed as the result of a catastrophic meteorite impact. To help reconcile some of the enigmas and apparent contradictions surrounding studies of the structure and to develop more effective geophysical techniques to locate new deposits, Lithoprobe partnered with industry to carry out geophysical surveys combined with the extensive geological information available. A revised structural model for the Sudbury structure was generated and a 3-D seismic reflection survey identified a nickel deposit, known from drilling results, prior to any mine development. The Athabasca Basin of northwestern Saskatchewan and northeastern Alberta is one of the world’s most prolific producers of uranium from its characteristically high-grade unconformity-type deposits and is the only current uranium producer in Canada. An extensive database of geology, drillhole data and physical properties exists. Working with industry collaborators, Lithoprobe demonstrated the value of high-resolution seismic for imaging the unconformity and faults associated with the deposits. The final example involves a unique seismic reflection experiment to image the diamondiferous Snap Lake kimberlite dyke in the Slave Province of the Northwest Territories. The opportunity to study geological samples of the kimberlite dyke and surrounding rocks and to ground-truth the seismic results with drillhole data made available by the two industry collaborators enabled a case history study that was highly successful.RÉSUMÉLithoprobe (1984-2005), ce projet de recherche pancanadien, multidisciplinaire et concerté en sciences de la Terre, a étudié la structure et l'évolution de la croûte continentale canadienne et de ses marges. Ça a été un projet très réussi et qui a redéfini la nature de la recherche en sciences de la Terre au Canada. L'une des nombreuses retombées de ce projet a démontré par l'exemple que les spécialistes des sciences de la Terre en géophysique et en géologie, y compris toutes les sous-disciplines applicables dans ces domaines d'étude généraux, doivent travailler de concert afin de parvenir à une interprétation exhaustive de tous les ensembles de données disponibles. Dans la partie 1, cette approche s'est concrétisée par des études portant sur les structures lithosphériques. Dans la partie 2, elle a produit un résumé des interprétations tirées de six études liées à l'exploration à partir de publications dans des revues scientifiques. Dans le premier exemple, les structures souterraines associées au batholite du ruisseau Guichon, dans le centre-sud de la Colombie-Britannique, et qui renferme des gisements porphyriques de cuivre et de molybdène, sont maintenant mieux définies et mieux reliées aux différentes phases géologiques du batholite. Un retraitement des données de sismique réflexion, et d’inversion magnétique et gravimétrique 2,5-D et 3-D combiné à une cartographie géologique détaillée et à des données de forage ont permis une interprétation révisée et améliorée du de subsurface. La recherche autour du gisement de sulfures massifs volcanogéniques de Bell Allard de la région de Matagami, dans le nord du Québec, est un deuxième exemple. Un levé de sismique réflexion réalisé au-dessus du gisement, peu après sa découverte par forage, couplé avec des diagraphies géophysiques et de carottes, a été réalisé pour vérifier si l'ensemble pouvait donner une image du gisement. La détection directe du gisement de minerai à partir de la coupe sismique serait difficile si son emplacement n'était pas déjà connu; cependant, les caractéristiques structurales qui peuvent être liées aux lithologies déduites des forages et des diagraphies ont été bien définies. Les gisements de nickel et les structures qui y sont reliées dans la bande de Thompson, à la limite ouest de la province du Supérieur, dans le nord du Manitoba, ont fait l'objet d'études sismiques et électromagnétiques (EM), combinés à des mesures de caractéristiques géologiques et physiques. L'image sismique/EM combinée indique que les roches du groupe d’intérêt d’Ospwagan, lesquelles ont une résistivité faible, s'étendent vers le sud-est sous le gneiss archéen et, les culminations structurales contrôlent la géométrie souterraine du groupe d’Ospwagan. La structure de Sudbury, en Ontario, est réputée pour ses gisements de nickel, les plus importants au monde, lesquels se sont formés à la suite d'un impact météoritique catastrophique. Pour aider à comprendre certaines des énigmes et résoudre d’apparentes contradictions entourant les études de la structure, et pour développer des techniques géophysiques plus efficaces afin de localiser de nouveaux gisements, Lithoprobe s'est associé à l'entreprise privée pour réaliser des levés géophysiques, et les comparer aux très nombreuses informations géologiques disponibles. Une révision du modèle structural du gisement de Sudbury, ajouté à un levé sismique réflexion tridimensionnelle, ont permis de circonscrire un gisement de nickel, avant tout autre travail de développement minier. Le bassin de l'Athabasca, dans le nord-ouest de la Saskatchewan et le nord-est de l'Alberta, est l'un des producteurs d'uranium les plus prolifiques au monde provenant de gisements à haute teneur de type discordant, et est le seul producteur d'uranium au Canada. Une volumineuse base de données sur la géologie, les forages et les propriétés physiques est disponible. En collaboration avec des entreprises privées, Lithoprobe a démontré la valeur de la sismique à haute résolution pour l'imagerie de la discordance et des failles associées aux gisements. Le dernier exemple est celui d'une expérience de sismique réflexion unique visant à représenter le dyke de kimberlite diamantifère du lac Snap dans la province des Esclaves, dans les Territoires du Nord-Ouest. L'occasion d'étudier des échantillons géologiques du dyke de kimberlite, et des roches environnantes, et de valider les résultats sismiques à l'aide des données de forage mises à disposition par les deux partenaires privés, a permis une étude de cas très fructueuse.

Lithoprobe (1984–2005), Canada’s national, collaborative, multidisciplinary, Earth Science research project, investigated the structure and evolution of the Canadian landmass and its margins. It was a highly successful project that redefined the nature of Earth science research in Canada. One of many contributions deriving from the project was the demonstration by example that Earth scientists from geophysics and geology, including all applicable sub-disciplines within these general study areas, must work together to achieve thorough and comprehensive interpretations of all available data sets. In Part 1, this statement was exemplified through studies involving lithospheric structures. In Part 2, it is exemplified by summarizing interpretations from six exploration-related studies derived from journal publications. In the first example, subsurface structures associated with the Guichon Creek batholith in south-central British Columbia, which hosts porphyry copper and molybdenum deposits, are better defined and related to different geological phases of the batholith. Reprocessed seismic reflection data and 2.5-D and 3-D inversions of magnetic and gravity data are combined with detailed geological mapping and drillhole information to generate the revised and improved subsurface interpretation. Research around the Bell Allard volcanogenic massive sulphide deposit in the Matagami region of northern Quebec provides the second example. A seismic reflection line over the deposit shortly after it was discovered by drilling, aided by core and geophysical logs, was acquired to test whether the deposit could be imaged. Direct detection of the ore body from the seismic section would be difficult if its location were not already known; however, structural characteristics that can be tied to lithologies from boreholes and logs were well identified. Nickel deposits and associated structures in the Thompson belt at the western limit of the Superior Province in northern Manitoba were the focus of seismic and electromagnetic (EM) studies combined with geology and physical property measurements. The combined seismic/EM image indicates that the rocks of the prospective Ospwagan Group, which have low resistivity, extend southeastward beneath the Archean gneiss and that structural culminations control the subsurface geometry of the Ospwagan Group. The Sudbury structure in Ontario is famous for its nickel deposits, the largest in the world, which formed as the result of a catastrophic meteorite impact. To help reconcile some of the enigmas and apparent contradictions surrounding studies of the structure and to develop more effective geophysical techniques to locate new deposits, Lithoprobe partnered with industry to carry out geophysical surveys combined with the extensive geological information available. A revised structural model for the Sudbury structure was generated and a 3-D seismic reflection survey identified a nickel deposit, known from drilling results, prior to any mine development. The Athabasca Basin of northwestern Saskatchewan and northeastern Alberta is one of the world’s most prolific producers of uranium from its characteristically high-grade unconformity-type deposits and is the only current uranium producer in Canada. An extensive database of geology, drillhole data and physical properties exists. Working with industry collaborators, Lithoprobe demonstrated the value of high-resolution seismic for imaging the unconformity and faults associated with the deposits. The final example involves a unique seismic reflection experiment to image the diamondiferous Snap Lake kimberlite dyke in the Slave Province of the Northwest Territories. The opportunity to study geological samples of the kimberlite dyke and surrounding rocks and to ground-truth the seismic results with drillhole data made available by the two industry collaborators enabled a case history study that was highly successful. RESUME Lithoprobe (1984-2005), ce projet de recherche pancanadien, multidisciplinaire et concerte en sciences de la Terre, a etudie la structure et l'evolution de la croute continentale canadienne et de ses marges. Ca a ete un projet tres reussi et qui a redefini la nature de la recherche en sciences de la Terre au Canada. L'une des nombreuses retombees de ce projet a demontre par l'exemple que les specialistes des sciences de la Terre en geophysique et en geologie, y compris toutes les sous-disciplines applicables dans ces domaines d'etude generaux, doivent travailler de concert afin de parvenir a une interpretation exhaustive de tous les ensembles de donnees disponibles. Dans la partie 1, cette approche s'est concretisee par des etudes portant sur les structures lithospheriques. Dans la partie 2, elle a produit un resume des interpretations tirees de six etudes liees a l'exploration a partir de publications dans des revues scientifiques. Dans le premier exemple, les structures souterraines associees au batholite du ruisseau Guichon, dans le centre-sud de la Colombie-Britannique, et qui renferme des gisements porphyriques de cuivre et de molybdene, sont maintenant mieux definies et mieux reliees aux differentes phases geologiques du batholite. Un retraitement des donnees de sismique reflexion, et d’inversion magnetique et gravimetrique 2,5-D et 3-D combine a une cartographie geologique detaillee et a des donnees de forage ont permis une interpretation revisee et amelioree du de subsurface. La recherche autour du gisement de sulfures massifs volcanogeniques de Bell Allard de la region de Matagami, dans le nord du Quebec, est un deuxieme exemple. Un leve de sismique reflexion realise au-dessus du gisement, peu apres sa decouverte par forage, couple avec des diagraphies geophysiques et de carottes, a ete realise pour verifier si l'ensemble pouvait donner une image du gisement. La detection directe du gisement de minerai a partir de la coupe sismique serait difficile si son emplacement n'etait pas deja connu; cependant, les caracteristiques structurales qui peuvent etre liees aux lithologies deduites des forages et des diagraphies ont ete bien definies. Les gisements de nickel et les structures qui y sont reliees dans la bande de Thompson, a la limite ouest de la province du Superieur, dans le nord du Manitoba, ont fait l'objet d'etudes sismiques et electromagnetiques (EM), combines a des mesures de caracteristiques geologiques et physiques. L'image sismique/EM combinee indique que les roches du groupe d’interet d’Ospwagan, lesquelles ont une resistivite faible, s'etendent vers le sud-est sous le gneiss archeen et, les culminations structurales controlent la geometrie souterraine du groupe d’Ospwagan. La structure de Sudbury, en Ontario, est reputee pour ses gisements de nickel, les plus importants au monde, lesquels se sont formes a la suite d'un impact meteoritique catastrophique. Pour aider a comprendre certaines des enigmes et resoudre d’apparentes contradictions entourant les etudes de la structure, et pour developper des techniques geophysiques plus efficaces afin de localiser de nouveaux gisements, Lithoprobe s'est associe a l'entreprise privee pour realiser des leves geophysiques, et les comparer aux tres nombreuses informations geologiques disponibles. Une revision du modele structural du gisement de Sudbury, ajoute a un leve sismique reflexion tridimensionnelle, ont permis de circonscrire un gisement de nickel, avant tout autre travail de developpement minier. Le bassin de l'Athabasca, dans le nord-ouest de la Saskatchewan et le nord-est de l'Alberta, est l'un des producteurs d'uranium les plus prolifiques au monde provenant de gisements a haute teneur de type discordant, et est le seul producteur d'uranium au Canada. Une volumineuse base de donnees sur la geologie, les forages et les proprietes physiques est disponible. En collaboration avec des entreprises privees, Lithoprobe a demontre la valeur de la sismique a haute resolution pour l'imagerie de la discordance et des failles associees aux gisements. Le dernier exemple est celui d'une experience de sismique reflexion unique visant a representer le dyke de kimberlite diamantifere du lac Snap dans la province des Esclaves, dans les Territoires du Nord-Ouest. L'occasion d'etudier des echantillons geologiques du dyke de kimberlite, et des roches environnantes, et de valider les resultats sismiques a l'aide des donnees de forage mises a disposition par les deux partenaires prives, a permis une etude de cas tres fructueuse.

Seismic inversion is applied to generate physical property models (P-wave velocity and numerical attenuation) for four profiles in the Nechako–Chilcotin plateau region of south-central British Columbia, Canada. A newly developed method that combines three-dimensional (3-D) travel-time inversion and 2.5-dimensional (2.5-D) viscoacoustic full-waveform inversion was applied to generate the geophysical models from vibroseis data acquired along the preexisting crooked roads. These models are useful for the characterization of rock types in terms of their positions and thicknesses, which may be used in conjunction with geological ground truth to infer the extent of lithostratigraphic units in the subsurface. The velocity structures also may be used for future reprocessing of the seismic reflection data to derive improved images based on the better near-surface velocity models. The subsurface geology of the Nechako–Chilcotin plateau region is complex, resulting from multiple stages of tectonic compression and extension, contemporaneous with the deposition of sediments and volcanic material. Several basin structures are identified from the joint interpretation of the waveform tomography velocity models and post-stack time migration images. The combination of these results enables the extrapolation and characterization of geological structures to ∼3 km depth, particularly within the Cenozoic volcanic units that dominate near-surface stratigraphy. Based on the seismic profiles, a fence-diagram geological interpretation that extends to ∼3 km depth illustrates the complex structure of the Jurassic to Neogene stratigraphic sequence.

[1] A method for 2.5D viscoacoustic waveform tomography that can be applied to generate 2D models of velocity and attenuation from inversion of refraction waveforms on land seismic reflection data acquired along crooked roads is developed. It is particularly useful for typical crustal reflection surveys. First-arrival travel time tomography is applied using a 3D method, but with constraints on the intermediate 3D velocity model; the result is the starting model for the next step. A 2.5D frequency-domain full-waveform inversion stage parameterizes 3D geometry in the seismic source and receiver arrays, with the assumption that the velocity and attenuation models are homogeneous in the out-of-plane direction. This approach results in superior results compared to a strictly 2D approach when the acquisition line is crooked, with a moderate increase in computational cost. A case study using data acquired in the Nechako Basin in south-central British Columbia, Canada, exemplifies and validates the procedure. The velocity model derived from 2.5D waveform tomography is compared with that from a previous study in which 2D waveform tomography was applied to the same data set and with results from 3D travel time tomography. The resolution and accuracy of the velocity model from 2.5D waveform tomography are demonstrated to be greater than those from travel time or 2D waveform tomography. A model of viscoacoustic attenuation, which was not possible in the 2D case, is also generated. These models are interpreted jointly to highlight features of geological interest, such as a sedimentary basin, basement rocks, and faults, from surface to about 3 km depth.

In the BATHOLITHSonland seismic project, a refraction – wide-angle reflection survey was shot in 2009 across the Coast Mountains and Interior Plateau of central British Columbia. Part of the seismic profile crossed the Nechako Basin, a Jurassic–Cretaceous basin with potential for hydrocarbons within sedimentary strata that underlies widespread volcanic rocks. Along this 205 km-long line segment, eight large explosive shots were fired into 980 seismometers. Forward and inverse modelling of the traveltime data were conducted with two independent methods: ray-tracing based modelling of first and secondary arrivals, and a higher resolution wavefront-based first-arrival seismic tomography. Material with velocities less than 5.0 km/s is interpreted as sedimentary rocks of the Nechako Basin, while velocities from 5.0–6.0 km/s may correspond to interlayered sedimentary and volcanic rocks. The greatest thickness of sedimentary rocks in the basin is found in the central 110 km of the profile. Two sub-basins were identified in this region, with widths of 20–50 km and maximum sedimentary depths of 2.5 and 3.3 km. Such features are well-defined in the velocity model, since resolution tests indicate that features with widths greater than ∼13 km are reliable. Beneath the sedimentary rocks, seismic velocities increase more slowly with depth — from 6.0 km/s just below the basin to 6.3 km/s at ∼17 km in depth, and then to 6.8–7.0 km/s at the base of the crust. The Moho is found at a depth of 33.5–35 km beneath the profile, and mantle velocities are high at 8.05–8.10 km/s.

Three lithospheric cross sections provide a continental-scale synthesis of more than two decades of coordinated multidisciplinary research during the Canadian Lithoprobe project. The sections are based on seismic reflection and refraction data combined with a broad range of geological, geochemical, geochronological, and geophysical data. The dataset is derived from remnants of nearly every kind of tectonic regime, and the geologic history of the entrained rocks spans the Present to the Mesoarchean. The longest of the three cross sections is located within a 6000 km long Trans-Canada corridor traversing the North American continent at 45°N–55°N. From west to east, the profile crosses the Juan de Fuca ridge and active Cascadia subduction zone, the Cordilleran, Albertan, and Trans-Hudson orogens, the Superior Province, the Midcontinent rift, the Grenville and Appalachian orogens, and the Atlantic passive margin. The two northern cross sections include (i) a 2000 km long corridor in northwestern Canada (54°N–63°N) crossing the Cordilleran, Wopmay, and Slave orogens; and (ii) a 1600 km long corridor in northeastern Canada (52°N–61°N) crossing the New Quebec and Torngat orogens, the Nain craton, and the Makkovik and Grenville orogens. The unprecedented scale of the cross sections illuminates the assembly of the North American continent. Relationships between orogens are emphasized; plate collisions and accretions have sequentially stacked orogen upon orogen such that the older crust forms basement to the next younger. The large-scale perspective of these regional sections highlights the subhorizontal Moho that is indicative of either structural or thermal re-equilibration (or both), as few crustal roots beneath orogens are preserved. In contrast, heterogeneities in the lithospheric mantle suggest that, in certain situations, relict subducted or delaminated lithosphere can remain intact beneath and eventually within cratonic lithospheric mantle.

Application of regional geophysical and geological methods throughout two decades of Canada’s Lithoprobe project provides new opportunities to analyze the Mohorovičić discontinuity (Moho) and crust–mantle transition. The transect format employed during Lithoprobe, in which 10 specified regions of Canada were targeted for approximately a decade each, between 1984 and 2003, permitted teams of scientists to focus on geological, geophysical, and tectonic issues for each transect. As a primary objective was to enhance knowledge of the structure of the crust and lithosphere, an obvious target in each transect was the nature and origin of the Moho and crust–mantle transition. Accordingly, the combined results provide new perspectives on the Moho and the relationship of the Moho to the crust–mantle transition. Perhaps the most important result is that the continental geophysical Moho is a deceptively simple feature; it has a variety of signatures at different scales that preclude a single, universally applicable interpretation. In methods that provide large-scale information, such as regional seismic studies, it is a relatively abrupt refraction velocity contrast that often displays a dramatic downward decrease in seismic reflectivity. However, its origin in a geological or tectonic sense is perhaps best determined by careful analyses of structural details near the geophysical Moho, which are complex and varied. In some areas within Canada, it appears that the geophysical Moho may be old and perhaps remains from the time the crust formed; in other areas, it appears to be a relatively young feature that was superimposed onto older crustal fabrics.

Lithoprobe is Canada’s national, collaborative, multidisciplinary, Earth science research project established to develop a comprehensive understanding of the evolution of northern North America. It is regarded internationally as one of the most successful national geoscience projects ever undertaken. Part of Lithoprobe’s success derives from its history and organizational structure. A one-year Phase I program (1984–1985) was highly successful. Phases II to V carried the project forward until 2005, when funding terminated. Lithoprobe was organized as a decentralized research network of multidisciplinary scientific studies and collaborating scientists. About 1500 scientific publications were generated. Lithoprobe’s legacy also includes economic and social benefits. New and improved understanding of Earth history provides petroleum and mining companies with an enhanced knowledge base. Lithoprobe demonstrated the applicability of high-resolution seismic reflection studies to mineral exploration. Development of a portable seismic refraction recorder and a long-period magnetotelluric system led to technology transfer to Canadian companies. Very high-resolution seismic reflection studies, applied first by Lithoprobe, are proving effective for exploration for uranium deposits in Saskatchewan. In the cratonic areas of Canada, Lithoprobe seismic and magnetotelluric studies have provided significant new information relevant to exploration for diamonds. On the west coast of Canada, Lithoprobe studies provided data and a framework for better understanding the mega-thrust earthquake hazard in the region. Lithoprobe established a highly effective public outreach strategy that involved the print and electronic media as well as material for educational purposes. Unequivocally, Lithoprobe has enhanced the already strong reputation of the Earth sciences in Canada.

Within Lithoprobe’s 10 transects, data from more than 20 000 km of multichannel seismic (MCS) reflection profiling and 12 refraction – wide-angle reflection (R/WAR) surveys were acquired. While the main results related to crustal structure, the data also indicated substantial heterogeneity in the lithospheric mantle. Images of fossilized subduction zones from the Eocene to the Neoarchean demonstrate that current plate tectonic processes have been active for more than 2.6 Ga. The Canadian Cordillera has a thin (50–60 km) lithosphere that is likely receiving some dynamic support from the asthenosphere below. Vestiges of the last stage of accretionary tectonic processes that formed the Archean Superior craton are indicated by an unusual anisotropic high velocity layer that may represent relic oceanic lithosphere. Within the Paleoproterozoic Trans-Hudson Orogen, a restricted region of upper mantle P-wave velocity anisotropy is identified with the continental collision between the bounding Hearne and Superior cratons. In the Archean Hearne and Wyoming provinces, two dipping structures within the sub-crustal lithosphere are interpreted as subduction features related to the assembly of the two cratons. Finite-difference modeling of long-offset data (over 1300 km) reveals fine-scale heterogeneities within a layer between 90 and 150 km in the continental lithosphere, perhaps formed through lateral flow or deformation within the upper mantle. Based on Lithoprobe data, heterogeneities within the lithospheric mantle are reasonably common. They have a wide range of seismic signatures, include many different types and show differing scales. Nevertheless, their extent in the lithospheric mantle is considerably less than in the crust.

Near-vertical incidence and refraction – wide-angle reflection seismic data, recorded as part of Lithoprobe studies in the Paleoproterozoic–Archean domains of Canada’s Northwest Territories, show remarkable reflections from within the upper mantle. A parallel pair of reflectors imaged by the near-vertical data can be traced from Moho levels (∼33 km) down to ∼70 km depth. In a previous study, the reflectors were interpreted as the top and bottom of an ∼1.8 Ga subducted oceanic crust beneath the Hottah terrane. Further inboard, where the seismic line changes its direction from east–west to nearly north–south, another pair of reflectors extends subhorizontally for about 100 km at ∼70 km depth before dipping downward. The subhorizontal reflectors were not correlated with the dipping slab; instead they were interpreted as a separate feature. However, they roughly coincide with a horizontal interface modeled from wide-angle data by an earlier study. Considering the crooked line acquisition geometry, we re-examined both near-vertical incidence and wide-angle reflection data using 2-dimensional (2-D) and 3-D forward and inverse modeling algorithms. Our results demonstrate that the subhorizontal reflectors are the continuation of the relict subducted slab, which now extends laterally for 300 km. Its base is the source of the wide-angle data. The apparent flattening for the near-vertical data is most likely an artifact of projecting a 3-D geometry onto a 2-D cross section. The shallowly subducted slab probably contributed to the thickening and stabilization of the subcrustal lithosphere below the Wopmay orogen.

SUMMARY Using a 3-D seismic reflection data set recorded by the Canadian petroleum exploration industry to 5.1 s two-way traveltime (∼15 km depth) in northwestern Alberta, the Winagami reflection sequence (WRS), an interpreted sill complex previously identified on Lithoprobe 2-D multichannel seismic reflection lines, is investigated to determine its 3-D geometry and reflective characteristics. Clear evidence of the top of the WRS emerges from the data set. Data sections outline a 3-D reflective sheet dipping to the southeast. From polarity comparisons with reflections from the sedimentary sequence, the reflective signature of the deep body is inferred to result from higher impedance material than the surrounding host rocks, thus reinforcing earlier interpretations that the deep reflections result from doleritic sills intruded into gneissic crystalline basement. Simple 1-D forward modelling reveals that the thickness of the sheet is on the order of 100 m throughout the survey region. From 3-D Kirchhoff forward modelling, a model of the reflective sheet is constrained between 11 and 16 km depth along the northwestern edge of the survey area. The absence of the Winagami reflections, and basement reflectivity in general, in the southeast of the survey region coincides with a positive aeromagnetic anomaly inferred to be caused by magmatic rocks. The presence of the magmatic rocks may have either influenced the geometry of the intrusion of the sills, overprinted their reflective signature or been co-magmatic, dependent upon the relative timing of emplacement of the two features. The Winagami sill complex appears to have been intruded horizontally into the crust during a period of tectonic compression. The emplacement of such deep Precambrian sills represents an under-appreciated mode of magmatic addition into the crystalline crust.

Shear-wave seismic data recorded along four profiles during the SNoRE 97 (1997 Slave – Northern Cordillera Refraction Experiment) refraction – wide-angle reflection experiment in northwestern Canada are analyzed to provide S-wave velocity (Vs) models. These are combined with previous P-wave velocity (Vp) models to produce cross sections of the ratio Vp/Vs for the crust and upper mantle. The Vp/Vs values are related to rock types through comparisons with published laboratory data. The Slave craton has low Vp/Vs values of 1.68–1.72, indicating a predominantly silicic crustal composition. Higher values (1.78) for the Great Bear and eastern Hottah domains of the Wopmay orogen imply a more mafic than average crustal composition. In the western Hottah and Fort Simpson arc, values of Vp/Vs drop to ∼1.69. These low values continue westward for 700 km into the Foreland and Omineca belts of the Cordillera, providing support for the interpretation from coincident seismic reflection studies that much of the crust from east of the Cordilleran deformation front to the Stikinia terrane of the Intermontane Belt consists of quartzose metasedimentary rocks. Stikinia shows values of 1.78–1.73, consistent with its derivation as a volcanic arc terrane. Upper mantle velocity and ratio values beneath the Slave craton indicate an ultramafic peridotitic composition. In the Wopmay orogen, the presence of low Vp/Vs ratios beneath the Hottah – Fort Simpson transition indicates the presence of pyroxenite in the upper mantle. Across the northern Cordillera, low Vp values and a moderate-to-high ratio in the uppermost mantle are consistent with the region's high heat flow and the possible presence of partial melt.

The SNORCLE refraction – wide-angle reflection (R/WAR) experiment, SNORE'97, included four individual lines along the three transect corridors. A combination of SNORE'97 results with those from earlier studies permits generation of a 2000 km long lithospheric velocity model that extends from the Archean Slave craton to the present Pacific basin. Using this model and coincident near-vertical incidence (NVI) reflection data and geological information, an interpreted cross section that exemplifies 4 Ga of lithospheric development is generated. The velocity structural models correlate well with the reflection sections and provide additional structural, compositional, and thermal constraints. Geological structures and some faults are defined in the upper crust. At a larger scale, the seismic data identify a variety of orogenic styles ranging from thin- to thick-skinned accretion in the Cordillera and crustal-scale tectonic wedging associated with both Paleoproterozoic and Mesozoic collisions. Models of Poisson's ratio support the NVI interpretation that a thick wedge of cratonic metasediments underlies the eastern accreted Cordilleran terranes. Despite the variety of ages, orogenic styles, and tectono-magmatic deformations that are spanned by the seismic corridors, the Moho remains remarkably flat and shallow (33–36 km) across the majority of the transect. Significant variations only occur at major tectonic boundaries. Laterally variable crustal velocities are consistently slower beneath the Cordillera than beneath the cratonic crust. This is consistent with the high temperatures (800–900 °C) required by the slow upper mantle velocities (7.8–7.9 km/s) observed beneath much of the Cordillera. Heterogeneity of the lithospheric mantle is indicated by wide-angle reflections below the Precambrian domains and the western Cordillera.

A three-dimensional model of the regional crustal architecture of the western Trans-Hudson Orogen, based on the interpretation of 590 km of deep-sounding seismic reflection data and a comparable length of existing seismic reflection information, is presented. The seismic images identify the regional geometry of the basal detachment zone (Pelican thrust) that separates juvenile allochthonous terranes from the underlying Archean microcontinent (Sask craton). The Sask Craton is inferred to have a minimum spatial extent of over 100 000 km2 with an associated crustal root that extends for 200 km along strike. During terminal collision, complete convergence of the Rae–Hearne and Superior continental blocks was precluded by the presence of the Sask Craton, resulting in the preservation of anomalous amounts of oceanic and associated sedimentary juvenile material. Along regional tectonic strike, consistency of crustal structure across the Rae–Hearne margin – Reindeer zone boundary is established. Several phases of tectonic development, including multistage subduction and continent–continent collision, are inferred for the western margin of the orogen. A bright, shallow (2–3.5 s two-way traveltime) band of reflectivity (Wollaston Lake reflector) imaged over ~150 000 km2 area is inferred to be a large post-orogenic mafic intrusion. A highly reflective, well-defined and structurally disturbed Moho discontinuity is mapped throughout the western Trans-Hudson Orogen. The present-day crustal architecture of the western Trans-Hudson Orogen is described in terms of the tectonic evolution within the region.

The Trans-Hudson Orogen (THO) is the world's largest Paleoproterozoic orogenic belt. Data from three refraction profiles are used to investigate its lithospheric structure in Saskatchewan and Manitoba. R1 crosses the orogen from the Hearne craton on the west to the Superior craton on the east; R2 and R3 are along the orogen. P-wave velocity structural models are generated using a ray-based technique. On line R1, higher crustal velocities in its eastern part coincide with rocks of the Flin Flon – Namew gneiss complex. Depth to Moho is in the 40–45 km range and equates to that from the reflection data, including a small crustal root below the Sask minicontinent. Along lines R2 and R3, depth to Moho varies from about 40 km up to 55 km at the north end of R2 and south end of R3. In general, variations in crustal velocity and depth to Moho do not correlate with the location and extent of geological domains; they appear to reflect the complex deformation and metamorphic history of the crustal rocks. Mantle velocities are high, ~8.2 km/s. However a limited area shows prominent velocity anisotropy, with values of 8.6 km/s along R2 and R3 and 8.1 km/s along R1. We speculate that the observed anisotropy represents an ~100-km-wide mantle suture zone resulting from the collision of Archean plates. The suture zone accommodated limited extensional deformation, associated with a counterclockwise rotation of the Superior plate, to generate the anisotropy. In this model, the lithospheric mantle of the THO internal domains and Sask craton are detached.

The Slave geological province, a relatively small craton in NW Canada that includes the oldest rocks in the world, is one of five cratons that form the Archaean continental core of North America. The Palaeoproterozoic Wopmay orogen, lying between the craton and the Phanerozoic northern Canadian Cordillera, represents crustal growth through addition of a series of north-south oriented magmatic arcs and accreted terranes. As part of LITHOPROBE's multidisciplinary studies, Line 11 of the SNoRE 97 seismic refraction/wide-angle reflection survey investigated the structure of the lithosphere in this region. Interpretation of crustal and mantle arrivals from 13 record sections followed an iterative procedure of traveltime inversion and forward modelling. Results show that crustal thickness at 32-35 km remains relatively constant from the Archaean regions in the east to Proterozoic domains in the west, contrary to average values of crustal thicknesses from global compilations that suggest thicker Proterozoic crust. At the westernmost end of the line, where crustal extension and basin formation are inferred, the crust thins to 30 km. Upper crustal velocities to about 12 km depth lie in the range of 5.8-6.2 km s - 1 . The middle crust of thickness 8-10 km has velocities between 6.2 and 6.6 km s - 1 . Average velocities for the lower crust, differentiated as a 10 km thick layer, range between 6.6 and 6.9 km s - 1 with a higher velocity region (7.1 km s - 1 ) in the transition zone between the Fort Simpson magmatic arc and Hottah accreted terrane in the Wopmay orogen, where a Palaeoproterozoic subduction zone has been inferred from coincident multichannel reflection data. In the same area, anomalously low upper-mantle velocities (7.5 km s - 1 ) have been found, indicating that the effects of collision, subduction and accompanying physical changes in the rocks comprising the mantle lithospheric wedge above the underthrusting structure, may have been retained for 1800 Ma. Data from most record sections with offset distances greater than 250 km show continuous wide-angle reflections from upper-mantle depths of about 75 km and almost 100 km. Although poorly constrained, velocities in this depth range vary from 8.4 to 8.8 km s - 1 . Interpretation of these phases and mantle reflections from equivalent depths in coincident, near-vertical incidence multichannel data relates them to subduction and underthrusting of oceanic crust and subsequent phase changes that occurred during the assembly of the craton and the evolution of the Palaeoproterozoic orogen.

The lithospheric structure of the transition from the craton to the Cordillera in northeastern British Columbia is interpreted from inversion of seismic refraction – wide-angle reflection data along a 460-km profile, and from 3-d (3-dimensional) inversion and 2.5-d forward modelling of Bouguer gravity data. The seismic profile extends westward from the sediment-covered edge of cratonic North America across the Foreland and Omineca morphogeological belts to the eastern boundary of accreted terranes, beyond the Tintina Fault. Across the ancient cratonic margin, the resultant models reveal a westward-thickening package of low upper crustal velocities (6.2 km/s and less) and low densities to almost 20 km depth below the Western Canada Sedimentary Basin, overlying a west-facing ramp of higher velocities and densities in the middle and lower crust. These features are inferred to represent passive-margin sediments deposited on the ancient rifted margin during the mid-to-late Proterozoic and early Paleozoic. A wedge-shaped high-velocity (7.3 km/s) crustal layer at the base of the crust beneath the edge of cratonic North America is interpreted to be the result of magmatic underplating during rifting. In the Cordilleran Foreland Belt, high velocities (6.4 km/s) in the upper 5 km of the crust indicate rocks upthrust from the middle crust. A narrow trench of low velocities in the near-surface, which is imaged ~20 km to the west of the inferred location of the Tintina Fault, is interpreted to represent the actual location of the fault or a major splay. From east to west, the Moho decreases in depth from ~40 km to ~34 km below the rifted margin of ancestral North America, then defines a small root at ~38 km depth below the high topography and upper crustal velocities of the eastern Foreland Belt, and gradually shallows to ~34 km beneath the Omineca belt. An enigmatic laterally heterogeneous upper mantle has anomalously high velocities (up to 8.3 km/s) beneath the Foreland Belt, flanked by regions of low velocities (7.7–7.8 km/s). Results indicate that the location of the Cordilleran deformation front west of the ramped cratonic margin directly affected the tectonic evolution of the region.

Three seismic refraction/wide-angle reflection profiles were recorded over the Trans-Hudson Orogen. The profiles display high quality coherent seismic signals associated with the Moho and lithospheric mantle to depths in excess of 160 km. As outlined by wide-angle reflections, depth to Moho varies between 40–54 km with a number of well-defined structural reliefs. Immediately beneath the Moho, in the central part of the orogen, an anomalous high-velocity zone (∼8.45 km/s) with lateral dimensions of 100 km (north-south) by 130 km (east-west) and maximum thickness of 50 km was delineated by mantle refraction arrivals. This zone may represent a locally preserved remnant of older mantle from micro-continental collisions. Mantle discontinuities with some structural relief are interpreted at average depths of 75 and 158 km. Based on wide-angle reflections at offsets beyond 400 km, they bound a layered zone with variable velocities, probably representing imbricated sequences of peridotite and eclogite. We postulate that these unusual features of the upper mantle are the result of lithospheric convergence of bounding Archean cratons and closure of the intervening ocean basin.

The northern Cascadia subduction zone is a region of convergence between the oceanic Explorer and northern Juan de Fuca plates and the continental North American plate. Potential field and new seismic reflection data coupled with previous seismic results and geology enable derivation of a series of density – magnetic susceptibility cross sections and a block density model from the ocean basin to the volcanic arc from 2.5- and 3-dimensional interpretations. The lateral extent and thickness of the accreted wedge vary significantly along the zone. The narrow, metasedimentary Pacific Rim terrane lies immediately west of Wrangellia and extends the length of Vancouver Island, consistent with its emplacement by strike-slip faulting following the accretion of Wrangellia. The ophiolitic Crescent terrane is a narrow slice lying seaward of the Pacific Rim terrane but not extending northward beyond the Juan de Fuca plate. In this region, the Crescent terrane was emplaced in a strike-slip or obliquely convergent style during the latter stages of emplacement of Pacific Rim terrane. Below the accreted terranes, the Explorer plate is shallower than Juan de Fuca plate, resulting in a thinner crust. High-density lower crustal material lies beneath the western edge of Vancouver Island, supporting interpretations of wide-scale underplating of Wrangellia. The shape of the boundary region between Wrangellia and the Coast belt to the east varies along strike and may be controlled by properties of preexisting plutonic rocks. The low-density Coast belt plutons extend throughout most of the crust and are underlain by a lowermost crustal high-density layer, which may be indicative of fractionation accompanying magma generation.

We applied an iterative combination of two-dimensional traveltime inversion and amplitude forward modelling to seismic refraction data along a 350 km along-strike profile in the Coast Belt of the southern Canadian Cordillera to determine crust and upper mantle P-wave velocity structure. The crustal model features a thin (0.5–3.0 km) near-surface layer with an average velocity of 4.4 km/s, and upper-, middle-, and lower-crustal strata which are each approximately 10 km thick and have velocities ranging from 6.2 to 6.7 km/s. The Moho appears as a 2 km thick transitional layer with an average depth of 35 km and overlies an upper mantle with a poorly constrained velocity of over 8 km/s. Other interpretations indicate that this profile lies within a collision zone between the Insular superterrane and the ancient North American margin and propose two collision-zone models: (i) crustal delamination, whereby the Insular superterrane was displaced along east-vergent faults over the terranes below; and (ii) crustal wedging, in which interfingering of Insular rocks occurs throughout the crust. The latter model involves thick layers of Insular material beneath the Coast Belt profile, but crustal velocities indicate predominantly non-Insular material, thereby favoring the crustal delamination model. Comparisons of the velocity model with data from the proximate reflection lines show that the top of the Moho transition zone corresponds with the reflection Moho. Comparisons with other studies suggest that likely sources for intracrustal wide-angle reflections observed in the refraction data are structural features, lithological contrasts, and transition zones surrounding a region of layered porosity in the crust.

The Queen Charlotte Basin consists of up to 6 km of Tertiary clastic sediments in a complex sequence of fault-bounded subbasins. The tectonic evolution of the basin in still being debated, with recent interpretations including distributed strike-slip extension, oblique or en echelon rifting, simple extension orthogonal to the plate margin, and block faulting and vertical tectonics. A combined seismic reflection and refraction survey was carried out in 1988 to investigate the structure and tectonic evolution of the basin and underlying crust. While the marine multichannel reflection data were being collected, refracted and wide-angle reflected energy from the large air-gun array was recorded at surrounding land sites in both two-dimensional (in-line) and three-dimensional (broadside) geometries. The broadside refraction data recorded on the Queen Charlotte Islands provide good three-dimensional coverage of western Hecate Strait. These data are modelled to determine the three-dimensional structure of the Queen Charlotte Basin. The reflection data indicate that the sedimentary Queen Charlotte Basin beneath the shotpoints varies rapidly in thickness and is highly three-dimensional. First-arrival traveltimes from the broadside refraction data are inverted to find the three-dimensional structure of the basement interface beneath the shots and out of the planes of the reflection sections. A map of basement depth is derived for a region several kilometres wide adjacent to the reflection lines. Basin thickness varies rapidly between ~ 200 m and ~ 6 km. The model is consistent with the seismic reflection and potential field data sets. Although most of the basin is modelled as sediments overlying rocks with crustal velocities, a thick sequence of interbedded sedimentary and volcanic rocks is interpreted to underlie the shot lines in one region that lies east of the central Queen Charlotte Islands. Four major faults are also interpreted. These are based on sharp vertical relief of over 2 km on the map of basement depth. The orientation and topography across the faults and the small lateral scale and large topographic changes of the subbasins support the distributed strike-slip extension evolutionary model for the basin.

As part of a multidisciplinary investigation of the structure and tectonics of the Queen Charlotte Basin and underlying crust, deep multichannel seismic reflection and coincident crustal refraction data were collected in 1988. Energy from the reflection air-gun array source was recorded at land sites at offsets appropriate to record crustal refraction and wide-angle reflection data. Refraction data recorded in a broadside geometry provide good three-dimensional coverage of western Hecate Strait. These data are modelled using tomographic inversion techniques to determine the three-dimensional velocity structure of the crust in this region. The one-dimensional average velocity increases rapidly with depth to 6.5 km/s at 7 km depth. Velocities from 7 to at least 12 km depth remain approximately constant and are associated with rocks of the Wrangellia terrane. Significant lateral velocity variations, including large differences in near-surface velocities attributable to surface features, relatively low velocities representing interbedded Tertiary sediments and volcanics, and a deep high-velocity anomaly that may represent the root of an igneous intrusion, are mapped. Wide-angle reflections from the Moho are used to determine the thickness of the crust. The Moho is at 29 km depth beneath the east coast of the Queen Charlotte Islands. This is deeper than the Moho observed below Queen Charlotte Sound and as deep as, or deeper than, that below Hecate Strait. Crustal thinning during Tertiary extension was thus greatest beneath the surface expression of the Queen Charlotte Basin, leaving the crust under the islands considerably thicker than under the basin. In an alternate or additional explanation, compression at the continental margin during the last 4 Ma may have been taken up by thickening or underplating of the continental crust beneath the islands. If the Pacific plate is subducting beneath the islands, the Moho observations constrain the slab to dip greater than 20–26°.

Seismic refraction data recorded along a 330 km cross-strike profile through the eastern Insular and southernmost Coast belts of the Canadian Cordillera are interpreted using an iterative combination of traveltime inversion and amplitude forward modelling. The resultant model is characterized by large lateral variations in velocity. The most significant of these variations is a decrease in upper and middle crustal velocities to the east of the surface trace of the Harrison fault, which likely represents the transition from crust of the Insular superterrane to that of the Intermontane superterrane. This interpretation is consistent with some present geological models that place the possible (probable) location of the suture between the two superterranes less than 20 km east of the Harrison fault. Velocities at the base of the upper crust average 6.4 and 6.2 km/s west and east of the fault, respectively. Mid-crustal velocities average 6.6–6.9 km/s to the west and 6.35–6.45 km/s to the east of the fault. Lower crustal velocities also decrease slightly to the east. Other features of the velocity model include (i) a thin near-surface layer with velocities between 2.5 and 6.1 km/s; (ii) upper crustal thickness of 12.5 km, thinning to 8 km at the eastern boundary of the Western Coast Belt (WCB); (iii) high velocity (6.6–6.9 km/s) mid-crustal layer west of the Harrison fault extending to 21 km depth; (iv) high-velocity (6.75–7.1 km/s) lower crustal layer; (v) low-velocity gradient upper mantle with depth to Moho at 34–37 km beneath most of the Coast Belt, decreasing to 30 km beneath the eastern Insular Belt, a depth much less than previous estimates. The inferred crustal velocity structure beneath the WCB is consistent with the three-layer electrical conductivity structure for this area derived from magnetotelluric surveys. The association of high resistivities with the upper crust suggests that the upper 8–12 km represents the massive cover of plutonic rocks which characterizes the WCB. Middle and lower crustal velocities beneath the WCB are consistent with Wrangellian velocities found beneath Vancouver Island, suggesting Wrangellia may extend at depth eastward as far as the Harrison fault.

COLLISIONS between continents generate the world's most striking mountains; the deep structure of these collision zones, inferred from seismic reflection profiles1–3, can be equally dramatic, if less well understood. Improved knowledge of the structure of ancient collisional belts, now known to have developed through plate tectonic processes4,5, can provide important new insights into the mechanisms of crustal growth. Although many such belts have been partly destroyed by severe post-collisional deformation6,7 or magmatism8, some have been preserved in the stable interiors of continents9: one example is the Early Proterozoic Trans-Hudson orogen9,10, in western Canada. Here we present seismic reflection images across the entire orogen, which reveal a broadly symmetric structure, with reflections linked to island-arc rocks dipping beneath both bounding Archaean cratons. The observed reflection geometries imply a doubling of crustal thickness during the collision, along crust-penetrating faults. The presence of a significant amount of Archaean crust at depth, the extent of which could not have been determined from surface mapping, suggests the need for caution in inferring crustal growth rate from the surface area of juvenile crust11,12.

Lithoprobe is Canada's national, collaborative, multidisciplinary earth science research program directed toward an enhanced understanding of how the North American continent evolved. Research in its eight transects or study areas, which span the country from Vancouver Island to Newfoundland and geological time from 4 Ga to the present, is spearheaded by seismic reflection surveys. These, combined with many other studies, are providing new insight into the varied tectonic processes that have been active in forming the continent. Results from the Southern Cordillera transect show that Mesozoic crustal growth occurred in the central and eastern Cordillera by the accretion and amalgamation of exotic terranes, the collision of which resulted in the generation of crustal-scale antiforms and duplexes. After the principal periods of compression, this area was affected by a major episode of extension that led to the unroofing of the metamorphic core complexes. Farther to the west, past and present subduction processes have eroded the lower lithosphere of accreted terranes and left underplated sediments and oceanic lithosphere. The Lithoprobe East transect, covering the Paleozoic Newfoundland Appalachians and Mesozoic rifted Atlantic margin, reveals three lower crustal blocks, each with distinctive reflection signatures on marine seismic data. Structures of the geologically established tectono-stratigraphic domains, imaged clearly by new onshore reflection data, sole at upper crustal to mid-crustal levels, suggesting that much of the surface stratigraphy is allochthonous to the lower crustal blocks. At the ocean–continent transition, interpretations suggest underplating of thinned continental crust by basaltic melt during the rifting process.In Lake Superior, data from the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE) transect reveal the complex structures of the late Middle Proterozoic Keweenawan rift, which is up to 35 km deep, that almost split North America. The GLIMPCE data in Lake Huron show a spectacular series of east-dipping crustal-scale reflections that coincide with the Grenville front tectonic zone. These and other data have led to a two-stage model involving collision of an exotic terrane with the southern Superior cratonic margin in the late Early Proterozoic followed by stacking–crustal penetrating imbrication and ramping associated with the Middle Proterozoic Grenvillian orogeny. The Archean Kapuskasing structural zone, a prominent northeast-trending feature that cuts obliquely across the dominant east-west structures of the Superior Province, is interpreted as a thin thrust sheet, soled by a variably reflective décollement, above which about 70 km of crustal shortening has occurred to bring mid-crustal to lower crustal rocks to the surface, and below which the Moho deepens. The shortening may have been accomplished by brittle faulting and erosion at levels above 20 km and ductile folding or faulting in the lower crust. Preliminary studies in the Archean Abitibi greenstone belt indicate that two major fault zones, the Larder Lake–Cadillac and Porcupine–Destor, which host significant mineralization, were generated by crustal-scale thrust and (or) strike-slip tectonics. Archean crustal sections are as structurally diverse and complex as their Proterozoic and Phanerozoic counterparts. The reflection Moho has highly variable characteristics as imaged within transects and among different transects. Crustal and Moho reflectivity observed in the various transects is caused by a wide range of features, including fault–shear zones, lithologic contacts, compositional layering, fluids in zones of high porosity, and anisotropy.

As part of the Lithoprobe Southern Cordillera transect, seismic refraction data were recorded along a 330 km long strike profile in the Intermontane belt. An iterative combination of two-dimensional traveltime inversion and amplitude forward modelling was used to interpret crust and upper mantle P-wave velocity structure. This region is characterized by (i) a thin near-surface layer with large variations in velocity between 2.8 and 5.4 km/s, and low-velocity regions that correlate well with surface expressions of Tertiary sedimentary and volcanic rocks; (ii) an upper and middle crust with low average velocity gradient, possibly a weak low-velocity zone, and lateral velocity variations between 6.0 and 6.4 km/s; (iii) a distinctive lower crust characterized by significantly higher average velocities relative to midcrustal values beginning at 23 km depth, approximately 8 km thick with average velocities of 6.5 and 6.7 km/s at top and base; (iv) a depth to Moho, as defined by wide-angle reflections, that averages 33 km with variations up to 2 km; and (v) a Moho transition zone of depth extent 1–3 km, below which lies the upper mantle with velocities decreasing from 7.9 km/s in the south to 7.7 km/s in the north. Where the refraction line obliquely crosses a Lithoprobe deep seismic-reflection profile, good agreement is obtained between the interpreted reflection section and the derived velocity structure model. In particular, depths to wide-angle reflectors in the upper crust agree with depths to prominent reflection events, and Moho depths agree within 1 km. From this comparison, the upper and middle crust probably comprise the upper part of the Quesnellia terrane. The lower crust from the refraction interpretation does not show the division into two components, parautochthonous and cratonic North America, that is inferred from the reflection data, indicating that their physical properties are not significantly different within the resolution of the refraction data. Based on these interpretations, the lower lithosphere of Quesnellia is absent and presumably was recycled in the mantle. At a depth of ~ 16 km below the Moho, an upper mantle reflector may represent the base of the present lithosphere.

Two multichannel seismic reflection lines recorded over the northern Juan de Fuca plate off western Canada provide images of the crustal structure of three pseudofaults, the traces of propagating rifts. Anomalous reflectivity, including both subhorizontal and dipping events at times below the normal crust-mantle transition, is prevalent in the region of the pseudofaults. Although the crustal thickness elsewhere on the lines varies by only about 10%, crust associated with the pseudofaults is as much as 30% thicker than average. This may be due to an increase in thickness of the interpreted oceanic crustal layer 3B or crustal underplating. The results indicate that magma supply at propagating rifts may at times be large, perhaps due to the effects of nearby hotspots or convective partial melting associated with the initiation of rifting.

Multichannel deep seismic reflection data from the subduction zone of western Canada delineate the wedge of accereted sediments and the principal terranes (Crescent, Pacific Rim, and Wrangellia) that form the convergent margin. The top of the igneous oceanic crust is defined by subhorizontal reflections extending at least 100 km landward of the deformation front. Upon incorporation into the accretionary wedge, the clearly defined stratigraphy of the incoming oceanic sedimentary section is destroyed over a distance of about 10 km. Initially, an unreflective zone, which correlates well with maximum fluid expulsion, is formed. Farther landward, a predominantly landward-dipping reflectivity exists. A number of reflections are thrust faults, which appear to merge at depth with the subhorizontal reflections, but most have another origin. These reflections may be related to the movement of fluids generated by the compaction of sediments or possibly by the dehydration of the subducting plate. They are strongest in a region of depressed seismic velocities beneath the continental slope, where an analysis of reflection amplitude with offset implies that a high Poisson's ratio exists; this is consistent with the presence of elevated pore pressures. Thus, pore pressure variations associated with the migration of fluids may be the cause of much of the reflectivity within the accreted wedge, although the precipitation of minerals from rising fluids could also be important. Evidence from the seismic data also indicates that fluids from the accretionary prism are being expelled into the sediments of the overlying Tofino basin. A number of anomalously strong reflections and disruption of the horizontally stratified sediments within the lower levels of the basin probably represent fluids that migrated upward from the accreted wedge and were trapped against impermeable barriers created through the deposition of sediments on the continental slope and in the basin.

[1] Factors controlling widespread cratonic arch formation across western North America from the Early Paleozoic to the Middle Mesozoic remain poorly understood. While a number of causes, such as mantle hot spots, lithospheric contraction, mantle flow, and lithospheric flexure due to surface loading, have been suggested to explain the localized uplift of basement blocks on otherwise tectonically quiescent margins, a variety of mechanisms may be at work. Here, we test the hypothesis that seismically detected packages of discrete midcrustal igneous sheet intrusions impart rheological contrasts within the crust, which act as seeds for cratonic arch formation in the presence of weak intraplate compression, using finite element modeling of intraplate deformation. We find that the seismically resolved sills can generate significant surface uplift in the presence of modest strain rates over reasonable orogenic timescales. Cratonic arch formation would be enhanced in the presence of additional thin, seismically unresolvable sills. While midcrustal sheet intrusions may be just one of many factors influencing cratonic arch formation, their presence in northern Alberta, Canada, can help explain the asymmetric configuration of the well-studied Peace River Arch.

[1] The Cordillera in northern Canada is underlain by westward tapering layers that can be followed from outcrops of Proterozoic strata in the Foreland belt to the lowermost crust of the orogenic interior, a distance of as much as 500 km across strike. They are interpreted as stratified Proterozoic rocks, including ∼1.8–0.7 Ga supracrustal rocks and their basement. The layering was discovered on two new deep seismic reflection profiles in the Yukon (Line 3; ∼650 km) and northern British Columbia (Line 2; ∼1245 km in two segments) that were acquired as part of the Lithoprobe Slave-Northern Cordillera Lithospheric Evolution (SNORCLE) transect. In the Mackenzie Mountains of the eastern Yukon, the layering in Line 3 is visible between 5.0 and 12.0 s (∼15 to 36 km depth). It is followed southwestward for nearly 650 km (∼500 km across strike) and thins to less than 1.0 s (∼3.0–3.5 km thickness) near the Moho at the Yukon-Alaska international boundary. In the northern Rocky Mountains of British Columbia, the upper part of the layering on Line 2 correlates with outcrops of Proterozoic (1.76–1.0 Ga) strata in the Muskwa anticlinorium. At this location, the layering is at least 15 km thick and is followed westward then southward into the middle and lower crust for ∼700 km (∼300 km across strike). It disappears as a thin taper at the base of the crust ∼150 km east of the coast of the Alaskan panhandle. The only significant disruption in the layering occurs at the Tintina fault zone, a late to postorogenic strike-slip fault with up to 800 km of displacement, which appears as a vertical zone of little reflectivity that disrupts the continuity of the deep layering on both profiles (∼300 km apart). The base of the layered reflection zone coincides with the Moho, which exhibits variable character and undulates in a series of broad arches with widths of ∼150 km. In general, the mantle appears to have few reflections. However, at the southwest end of Line 3 near the Alaska-British Columbia border, a reflection dips eastward from ∼14.0 s to ∼21.0 s (∼45 to 73 km depth) beneath exposed Eocene magmatic rocks. It is interpreted as a relict subduction surface of the Kula plate. Our interpretation of Proterozoic layered rocks beneath most of the northern Cordillera suggests a much different crustal structure than previously considered: (1) Ancient North American crust comprising up to 25 km of metamorphosed Proterozoic to Paleozoic sediments plus 5–10 km of pre-1.8 Ga crystalline basement projects westward beneath most of the northern Canadian Cordillera. (2) The lateral (500 km by at least 1000 km) and vertical (up to 25 km) extent of the Proterozoic layers and their internal deformation are consistent with a long-lived margin for northwestern North America with alternating episodes of extension and contraction. (3) The detachments that carry deformed rocks of the Mackenzie Mountains and northern Rocky Mountains are largely confined to the upper crustal region above the layering. (4) Accreted terranes include thin klippen that were thrust over North American pericratonic strata (e.g., Yukon-Tanana), and terranes such as Nisling and Stikinia that thicken westward as the underlying Proterozoic layers taper and disappear. (5) The ages of exposed rocks are not necessarily indicative of the ages of underlying crust, a frequent observation in Lithoprobe interpretations, so that estimates of crustal growth based on surface geology may not be representative.

Precambrian cratons—the relatively stable areas of the Earth's crust that form the core of continents—hold the key to understanding how continents evolve, but vast cratonic areas are concealed from scientific probing beneath Phanerozoic sediments. Working for the Alberta Basement Transect (ABT), a part of the Canadian Lithoprobe project, geologists are examining the evolution of the petroleum-rich Western Canada Sedimentary Basin (WCSB). Armed with a variety of techniques, the scientists are resolving features within the sedimentary basin and unravelling the genesis of the underlying Precambrian crust and upper mantle. This represents a new view of a significant and largely unknown part of cratonic North America.

Seismic reflection data show that a regionally extensive band of landward-dipping reflections exists above the subducting Juan de Fuca plate at the western Canadian continental margin. The reflections truncate at depth a major terrane boundary mapped near the surface, and must have a structural origin. Many reflections are very strong, having reflection coefficients locally as high as 10%-20%.These amplitudes are consistent with theoretical calculations of reflections from sharp poirosity contrasts, if the pores are very thin; however, seismic data from the accreted sedimentary melange farther seaward, from which much of the deep reflective zone is derived, suggest that significant porosity contrasts do not occur alone. Because higher reflection amplitudes are found in regions of inferred high strain, both geometric and amplitude constraints imply that the deep reflective layering represents intensely sheared rocks, possibly originating close to the subduction decollement.

The Strait of Georgia, a topographic depression between Vancouver Island and the mainland of British Columbia, is considered to be the boundary between two tectonic provinces: the Coast Plutonic Complex on the east and the Insular Belt to the west. The allochthonous nature of the Insular Belt has been established, mainly on the basis of paleomagnetic measurements. Various tectonic models to explain the geological differences between the two provinces have been proposed. One of these suggests that the boundary is an old transform fault zone and is represented currently by a thrust fault along the eastern side of the Strait of Georgia. Other models propose that the Coast Plutonic Complex is a feature superimposed by tectonic and metamorphic events after the accretion of the Insular Belt. Such models do not require a major crustal discontinuity along the Strait of Georgia.In May 1982, a seismic refraction survey using a 32 L air gun and a radio telemetering sonobuoy system was carried out in the Strait of Georgia with the objective of investigating the nature of this boundary and determining the upper crustal structure. Three reversed profiles across the strait were shot; these are supplemented by several high-resolution reflection profiles from previous experiments. Two-dimensional models of the crustal structure across the strait have been constructed using a forward modelling ray trace and synthetic seismogram algorithm to match the travel times and amplitude characteristics of the data.Three basic layers or strata form the models, for which the maximum depth of reliability is 3 km. The first layer consists of unconsolidated sediments and Pleistocene glacial deposits, and the second represents Late Cretaceous – early Tertiary basin fill sediments that form the Nanaimo Group, the Burrard–Kitsilano formations, and the Chuckanut Formation. The third layer is likely the extension of the Coast Plutonic Complex beneath the strait, but the westerly limit of this unit is undetermined because of seismic properties similar to those of the Insular Belt volcanics. A local fault is located ~15 km northeast of Galiano Island on the west side of the strait. However, our study shows no evidence for a major fault along the strait. Thus those aspects of tectonic models that require the existence of a major transform or transcurrent fault boundary along the Strait of Georgia. may have to be revised.

Summary. The Winona Basin is a young deep water basin along the base of the continental slope off central British Columbia which is filled primarily with turbidite sediments. Seismic reflection profiles in the basin indicate that in places the sediment thickness is greater than 4.0 s two-way travel time, and simple interpretations of seismic refraction data indicate that velocities up to 6.0 kni s-* exist within the sediment section. Because the turbidite sediments are acoustically layered with large velocity contrasts between layers, the sediments are probably highly anisotropic. A simple model to account for this has been used to estimate the errors involved in using refraction data to determine layer thicknesses and average velocities. If the effects of anisotropy are ignored, as they usually are, values for thickness and velocity may be overestimated by as much as 10-20 per cent. Even the corrected velocities are still high, however, and the inferred high densities and high degrees of lithification are remarkable when the youth of the sediments is considered. Nearly the entire depositional history is contained within the Pleistocene, and the age of the sediments whose velocity exceeds 4 km s-l is estimated to be only 1 .O Myr old. Corrected velocities in the deepest part of the basin, where the sediment section is over 6 kin thick, exceed 5 km s-'. The nature of samples dredged from the basin (mudstones to coarse-grained sandstones and conglomerates) suggests a mechanism for initial diagenesis: The mudstones are moderately porous, 40-60 per cent, and possess little carbonate, whereas the coarser grained rocks are all well cemented with calcium carbonate, with carbonate contents ranging up to 37 per cent. Measured porosities are as low as 8 per cent, but porosities calculated on a carbonate-free basis range from 35 to 55 per cent. It appears that the fine-grained sediments are first consolidated primarily by gravitational and/or tectonic compaction, and the coarsegrained sediments are cemented by carbonate which is precipitated from pore fluids escaping through the highly permeable coarse-grained layers. As lithification deeper in the basin continues, pore fluids must continue to escape; tectonic fracturing may facilitate this.

One of three explanations of a prominent time delay in the 6.5 km/s branch of a recent seismic refraction survey in the Rocky Mountain Trench suggested a high-angle crustal fault crossing the trench near Radium, British Columbia. If the density contrast between basement and cover rocks is 0.1 g/cm3, a gravity anomaly of approximately 18 mGal should be observed. To test the fault hypothesis, a gravity survey has been carried out in and adjacent to the trench in the Radium area. The resultant data are not consistent with the proposed fault model. The principal feature of the data is a pronounced low, which coincides with the trench throughout the survey area. The low is due to Cenozoic fill and interpretation by two-dimensional modelling indicates the thickness of fill is about 550 m to the north and 420 m to the south of Radium. As a result of this survey, the two alternative hypotheses to explain the seismic data must be reconsidered. These are (1) the existence of a crustal low velocity zone, and (2) a major deformation of the basement and overlying rocks due to the trench being an ancient zone of weakness, which coincides with the western limit of the continental Precambrian craton. As reflections from the top of a low velocity zone are not observed, the second alternative is preferred.

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.

Summary. Multichannel seismic reflection sections recorded across Vancouver Island have revealed two extensive zones of deep seismic reflections that dip gently to the northeast, and a number of moderate northeasterly dipping reflections that can be traced to the surface where major faults are exposed. Based on an integrated interpretation of these data with information from gravity, heat flow, seismicity, seismic refraction, magnetotelluric and geological studies it is concluded that the lower zone of gently dipping reflections is due to underplated oceanic sediments and igneous rocks associated with the current subduction of the Juan de Fuca plate, and that the upper zone represents a similar sequence of accreted rocks associated with an earlier episode of subduction. The high densitylhigh velocity material between the two reflection zones is either an underplated slab of oceanic lithosphere or an imbricated package of mafic rocks. Reprocessing of data from two of the seismic lines has produced a remarkable image of the terrane bounding Leech River fault, with its dip undulating from > 60" near the surface to 20" at 3 km depth and -38" at 6 km depth.

The Queen Charlotte Basin, which lies on the west coast of Canada, was the location of a petroleum exploration program in the 1960's. Recently published tectonic models indicate that estimates of hydrocarbon potential require re-evaluation, and a renewed exploration interest has been expressed. In 1981, a seismic refraction survey of the upper crust using radio-telemetering sonobuoys and a 27 L air-gun source was carried out. A particular objective of the study was to determine the existence and depth extent of any sedimentary layer, hypothesized on the basis of other studies, beneath or within the Tertiary Masset volcanics in which some of the exploration wells had terminated. Three reversed profiles and one unreversed profile up to 40 km long were recorded. Interpretation of the data made use of the travel-time and amplitude information of the seismic sections by comparison with theoretical sections computed by two-dimensional ray tracing and a new asymptotic ray theory synthetic seismogram algorithm.Consistent with the earlier industry results, sediment thicknesses vary considerably throughout the southern Queen Charlotte Basin. The Tertiary Masset volcanics appear to be pervasive throughout the study area, with thicknesses varying from less than 1 km to greater than 3 km. On three of the four profiles a low-velocity layer, interpreted as Mesozoic sediments or sediments interbedded with volcanics, was found to lie beneath the volcanics. Thicknesses ranged from about 1 km to zero at a pinchout. The lowermost layer of all models is considered to be crustal rocks and is identified with the top of Wrangellia, an allochthonous terrane proposed to underlie the southern Queen Charlotte Basin. Along the profile for which no low-velocity layer sediments were interpreted, the Wrangellia terrane forms a dome rising to within 2 km of the surface. Other recent studies suggest that hydrocarbon sources could be associated with the Mesozoic rocks of Wrangellia and with any sediments underlying the Tertiary lavas, as well as with Tertiary marine sediments above the volcanics. Thus further exploration is warranted.

Summary. A comprehensive interpretation of two shear wave refraction profiles recorded on ocean bottom seismometers provides velocity-depth models for the oceanic crust and upper mantle. The data were obtained as part of an investigation of lithospheric structure carried out in the Nootka fault zone area off the west coast of Canada. One-dimensional S-wave velocity-depth models are derived from the data by trial and error comparisons with WKBJ synthetic seismograms. The results of this study confirm crustal characteristics interpreted previously from the compressional wave data. Layers 3A and 2B in the upper crust are characterized by rapid increases in seismic velocities; shear-wave velocity increases from 2.2 to 3.6 km s-l over an interval of 3.5 km, although the exact nature of the shearwave velocity gradient is uncertain. In contrast, the shear velocity in layer 3A remains fairly constant with depth at values near 3.6km s-'. No velocity discontinuity is observed at the Moho. Instead, the transition from crust to upper mantle is marked by the presence of a velocity gradient zone, layer 3B. The shear-wave velocities in the upper mantle are 4.5 and 4.6 km s-l in directions approximately parallel to the spreading ridge and parallel to the direction of spreading respectively. This low degree of velocity anisotropy contrasts sharply with the 10 per cent anisotropic effect observed for the Pn velocities. Combining the P-wave interpretation with this S-wave study enables evaluation of Poisson's ratio as a function of depth. Except for low values of Poisson's ratio found in the upper crust, the seismic properties of the oceanic crust and upper mantle interpreted in this study agree remarkably well with those of laboratory measurements on ophiolite samples.