Crustal velocity structure of the southern Nechako basin, British Columbia, from wide-angle seismic traveltime inversion1This article is one of a series of papers published in this Special Issue on the theme of New insights in Cordilleran Intermontane geoscience: reducing exploration risk in the mountain pine beetle-affected area, British Columbia.

2011 ◽  
Vol 48 (6) ◽  
pp. 1050-1063 ◽  
Author(s):  
A.L. Stephenson ◽  
G.D. Spence ◽  
K. Wang ◽  
J.A. Hole ◽  
K.C. Miller ◽  
...  

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.

1998 ◽  
Vol 35 (5) ◽  
pp. 556-561 ◽  
Author(s):  
P J Patchett ◽  
G E Gehrels ◽  
C E Isachsen

Nd isotopic data are presented for a suite of metamorphic and plutonic rocks from a traverse across the Coast Mountains between Terrace and Prince Rupert, British Columbia, and for three contrasting batholiths in the Omineca Belt of southern Yukon. A presumed metamorphic equivalent of Jurassic volcanic rocks of the Stikine terrane gives epsilon Nd = +6, and a number of other metaigneous and metasedimentary rocks in the core of the Coast Mountains give epsilon Nd values from +3 to +7. A single metasedimentary rock approximately 3 km east of the Work Channel shear zone gives a epsilon Nd value of -9. Coast Belt plutons in the traverse yield epsilon Nd from -1 to +2. The Omineca Belt plutons give epsilon Nd from -10 to -17. All results are consistent with published data in demonstrating that (i) juvenile origins for both igneous and metamorphic rocks are common in the Coast Belt; (ii) representatives of a continental-margin sedimentary sequence with Precambrian crustal Nd are tectonically interleaved in the Coast Mountains; (iii) Coast Mountains plutons can be interpreted as derived from a blend of metamorphic rocks like those seen at the surface, or as arc-type melts contaminated with the older crustal component; and (iv) Omineca Belt plutons are dominated by remelted Precambrian crustal rocks.


2019 ◽  
Vol 218 (1) ◽  
pp. 45-56 ◽  
Author(s):  
C Nur Schuba ◽  
Jonathan P Schuba ◽  
Gary G Gray ◽  
Richard G Davy

SUMMARY We present a new approach to estimate 3-D seismic velocities along a target interface. This approach uses an artificial neural network trained with user-supplied geological and geophysical input features derived from both a 3-D seismic reflection volume and a 2-D wide-angle seismic profile that were acquired from the Galicia margin, offshore Spain. The S-reflector detachment fault was selected as the interface of interest. The neural network in the form of a multilayer perceptron was employed with an autoencoder and a regression layer. The autoencoder was trained using a set of input features from the 3-D reflection volume. This set of features included the reflection amplitude and instantaneous frequency at the interface of interest, time-thicknesses of overlying major layers and ratios of major layer time-thicknesses to the total time-depth of the interface. The regression model was trained to estimate the seismic velocities of the crystalline basement and mantle from these features. The ‘true’ velocities were obtained from an independent full-waveform inversion along a 2-D wide-angle seismic profile, contained within the 3-D data set. The autoencoder compressed the vector of inputs into a lower dimensional space, then the regression layer was trained in the lower dimensional space to estimate velocities above and below the targeted interface. This model was trained on 50 networks with different initializations. A total of 37 networks reached minimum achievable error of 2 per cent. The low standard deviation (<300  m s−1) between different networks and low errors on velocity estimations demonstrate that the input features were sufficient to capture variations in the velocity above and below the targeted S-reflector. This regression model was then applied to the 3-D reflection volume where velocities were predicted over an area of ∼400 km2. This approach provides an alternative way to obtain velocities across a 3-D seismic survey from a deep non-reflective lithology (e.g. upper mantle) , where conventional reflection velocity estimations can be unreliable.


Geophysics ◽  
2008 ◽  
Vol 73 (3) ◽  
pp. S63-S71 ◽  
Author(s):  
Rongrong Lu ◽  
Mark Willis ◽  
Xander Campman ◽  
Jonathan Ajo-Franklin ◽  
M. Nafi Toksöz

We describe a new shortcut strategy for imaging the sediments and salt edge around a salt flank through an overburden salt canopy. We tested its performance and capabilities on 2D synthetic acoustic seismic data from a Gulf of Mexico style model. We first redatumed surface shots, using seismic interferometry, from a walkaway vertical seismic profile survey as if the source and receiver pairs had been located in the borehole at the positions of the receivers. This process creates effective downhole shot gathers by completely moving surface shots through the salt canopy, without any knowledge of overburden velocity structure. After redatuming, we can apply multiple passes of prestack migration from the reference datum of the bore-hole. In our example, first-pass migration, using only a simple vertical velocity gradient model, reveals the outline of the salt edge. A second pass of reverse-time, prestack depth migration using full two-way wave equation was performed with an updated velocity model that consisted of the velocity gradient and salt dome. The second-pass migration brings out dipping sediments abutting the salt flank because these reflectors were illuminated by energy that bounced off the salt flank, forming prismatic reflections. In this target-oriented strategy, the computationally fast redatuming process eliminates the need for the traditional complex process of velocity estimation, model building, and iterative depth migration to remove effects of the salt canopy and surrounding overburden. This might allow this strategy to be used in the field in near real time.


Geophysics ◽  
2001 ◽  
Vol 66 (6) ◽  
pp. 1925-1936 ◽  
Author(s):  
Moritz M. Fliedner ◽  
Robert S. White

We use the wide‐angle wavefield to constrain estimates of the seismic velocity and thickness of basalt flows overlying sediments. Wide angle means the seismic wavefield recorded at offsets beyond the emergence of the direct wave. This wide‐angle wavefield contains arrivals that are returned from within and below the basalt flows, including the diving wave through the basalts as the first arrival and P‐wave reflections from the base of the basalts and from subbasalt structures. The velocity structure of basalt flows can be determined to first order from traveltime information by ray tracing the basalt turning rays and the wide‐angle base‐basalt reflection. This can be refined by using the amplitude variation with offset (AVO) of the basalt diving wave. Synthetic seismogram models with varying flow thicknesses and velocity gradients demonstrate the sensitivity to the velocity structure of the basalt diving wave and of reflections from the base of the basalt layer and below. The diving‐wave amplitudes of the models containing velocity gradients show a local amplitude minimum followed by a maximum at a greater range if the basalt thickness exceeds one wavelength and beyond that an exponential amplitude decay. The offset at which the maximum occurs can be used to determine the basalt thickness. The velocity gradient within the basalt can be determined from the slope of the exponential amplitude decay. The amplitudes of subbasalt reflections can be used to determine seismic velocities of the overburden and the impedance contrast at the reflector. Combining wide‐angle traveltimes and amplitudes of the basalt diving wave and subbasalt reflections enables us to obtain a more detailed velocity profile than is possible with the NMO velocities of small‐offset reflections. This paper concentrates on the subbasalt problem, but the results are more generally applicable to situations where high‐velocity bodies overlie a low‐velocity target, such as subsalt structures.


1998 ◽  
Vol 35 (11) ◽  
pp. 1238-1251 ◽  
Author(s):  
Deping Chian ◽  
François Marillier ◽  
Jeremy Hall ◽  
Garry Quinlan

New modelling of wide-angle reflection-refraction data of the Canadian Lithoprobe East profile 91-1 along the central mobile belt of the Newfoundland Appalachian orogen reveals new features of the upper mantle, and establishes links in the crust and upper mantle between existing land and marine wide-angle data sets by combining onshore-offshore recordings. The revised model provides detailed velocity structure in the 30-34 km thick crust and the top 30 km of upper mantle. The lower crust is characterized by a velocity of 6.6-6.8 km/s onshore, increasing by 0.2 km/s to the northeast offshore beneath the sedimentary basins. This seaward increase in velocity may be caused by intrusion of about 4 km of basic rocks into the lower crust during the extension that formed the overlying Carboniferous basins. The Moho is found at 34 km depth onshore, rising to 30 km offshore to the northeast with a local minimum of 27 km. The data confirm the absence of deep crustal roots under the central mobile belt of Newfoundland. Our long-range (up to 450 km offset) wide-angle data define a bulk velocity of 8.1-8.3 km/s within the upper 20 km of mantle. The data also contain strong reflective phases that can be correlated with two prominent mantle reflectors. The upper reflector is found at 50 km depth under central Newfoundland, rising abruptly towards the northeast where it reaches a minimum depth of 36 km. This reflector is associated with a thin layer (1-2 km) unlikely to coincide with a discontinuity with a large cross-boundary change in velocity. The lower reflector at 55-65 km depths is much stronger, and may have similar origins to reflections observed below the Appalachians in the Canadian Maritimes which are associated with a velocity increase to 8.5 km/s. Our data are insufficient for discriminating among various interpretations for the origins of these mantle reflectors.


2011 ◽  
Vol 48 (6) ◽  
pp. 1021-1037 ◽  
Author(s):  
A.J. Calvert ◽  
N.E. Hayward ◽  
J.E. Spratt ◽  
J.A. Craven

In 2008, a Vibroseis seismic reflection survey was acquired by Geoscience BC across the eastern part of the volcanic-covered Nechako basin in central British Columbia, where Cretaceous sedimentary rocks have been exhumed along a NNW trend. Good signal penetration through the volcanic cover is indicated by lower crustal reflections at 8–12 s, which were recorded by the entire seismic survey. Comparison of the 2008 seismic survey with data from a previous survey indicates that the lack of reflectivity in the earlier surveys is generally representative of the subsurface geology. The seismic data show that ∼1700 and ∼2900 m thick sub-basins are present at the northern and southern ends of this trend, but the intervening Cretaceous rocks are discontinuous and relatively thin. The creation of a passive-roof duplex by Campanian or later low-angle thrusting is inferred within the thickest Cretaceous strata, but elsewhere faulting is likely related to Eocene extension or transtension. Seismic reflections are also recorded from folded volcanic stratigraphy, the base of the surface volcanic rocks, an underlying volcaniclastic stratigraphy, and intrusions projecting into a Quaternary volcanic cone. Seismic interpretation is complemented by coincident audiofrequency magnetotelluric surveys, from which faulting is inferred at offsets in a regional conductor. No regionally extensive stratigraphy can be identified within the seismic data, and the central Nechako basin appears to be a complex network of small, deformed sub-basins, rather than a single large basin.


1997 ◽  
Vol 34 (12) ◽  
pp. 1603-1618 ◽  
Author(s):  
Carmala N. Garzione ◽  
P. Jonathan Patchett ◽  
Gerald M. Ross ◽  
JoAnne Nelson

Nd isotopes and trace elements in sedimentary rocks of the Yukon, the Northwest Territories, and northern British Columbia are used to examine the source of sediments in the Canadian Cordilleran miogeocline. Previous Nd isotope studies in southern Alberta demonstrated that strata of Neoproterozoic to Late Ordovician age were derived from Archean and Proterozoic Canadian Shield sources, whereas by the Late Devonian, a shift of 6 εNd units to younger crustal sources (εNd (T) = −6 to −9) had occurred. In this study, we found that the shift to younger crustal Nd isotopic signatures in the Yukon and Northwest Territories occurred much earlier than in southern Alberta. Cambrian and older strata have εNd(T) values of −10.0 to −21.1, consistent with derivation from Canadian Shield sources. Lower Ordovician through Permian strata in the Yukon and Northwest Territories, including the Innuitian-derived Imperial Assemblage, have εNd(T) values of −5 to −11.4. In northern British Columbia, the shift to a younger source reflects a wider range of εNd(T) values, from -−8.7 to −14.6 in Middle Ordovician through Middle Devonian strata, suggesting continued input from Canadian Shield sources. By the Middle Devonian, a complete shift to younger crustal signatures (εNd(T) = −5.9 to −10.5) had occurred in northern British Columbia. Several sources for the more juvenile sediments include (1) a mixture of locally erupted volcanic rocks with Canadian Shield sources, (2) a Grenville source, and (3) an Innuitian source. We propose that Ordovician to Lower Devonian strata were derived from a mixture of locally erupted, juvenile volcanics and pre-Cambrian Canadian Shield sources, and post-Middle Devonian strata were sourced from the Innuitian orogen in the Canadian Arctic.


1990 ◽  
Vol 27 (11) ◽  
pp. 1456-1461 ◽  
Author(s):  
R. M. Friedman ◽  
J. W. H. Monger ◽  
H. W. Tipper

A new U–Pb date of [Formula: see text] for foliated felsic metavolcanic rocks of the Bowen Island Group, from Mount Elphinstone in the southwesternmost Coast Mountains of British Columbia, indicates that there the age of this hitherto undated unit is early Middle Jurassic. These rocks grade along strike to the north-northwest into a more sedimentary facies, which north of Jervis Inlet contains a probable Sinemurian (Lower Jurassic) ammonite. The Bowen Island Group thus appears to include Lower and Middle Jurassic rocks and to be coeval in part with volcanic rocks of the Bonanza Formation on Vancouver Island to the west and the Harrison Lake Formation within the central Coast Mountains 75 km to the east.


2017 ◽  
Vol 54 (2) ◽  
pp. 163-172 ◽  
Author(s):  
Shutian Ma ◽  
Pascal Audet

Models of the seismic velocity structure of the crust in the seismically active northern Canadian Cordillera remain poorly constrained, despite their importance in the accurate location and characterization of regional earthquakes. On 29 August 2014, a moderate earthquake with magnitude 5.0, which generated high-quality Rayleigh wave data, occurred in the Northwest Territories, Canada, ∼100 km to the east of the Cordilleran Deformation Front. We carefully selected 23 seismic stations that recorded the Rayleigh waves and divided them into 13 groups according to the azimuth angle between the earthquake and the stations; these groups mostly sample the Cordillera. In each group, we measured Rayleigh wave group velocity dispersion, which we inverted for one-dimensional shear-wave velocity models of the crust. We thus obtained 13 models that consistently show low seismic velocities with respect to reference models, with a slow upper and lower crust surrounding a relatively fast mid crustal layer. The average of the 13 models is consistent with receiver function data in the central portion of the Cordillera. Finally, we compared earthquake locations determined by the Geological Survey of Canada using a simple homogenous crust over a mantle half space with those estimated using the new crustal velocity model, and show that estimates can differ by as much as 10 km.


2011 ◽  
Vol 48 (6) ◽  
pp. 1065-1089 ◽  
Author(s):  
M.D. Thomas ◽  
M. Pilkington ◽  
R.G. Anderson

The ability of airborne sensors to image the magnetic signatures of prospective Quesnel terrane rocks through ubiquitous Quaternary glacial sedimentary cover in central British Columbia helps target new areas for mineral exploration. Newly acquired high-resolution data provide new perspectives on the nature and probable areal distribution of many geological units, revealing detail and information unattainable by conventional geological mapping. In combination with gravity data, these magnetic data indicate the presence of a granitic intrusion and a development of Nicola Group volcanic rocks, both potential hosts for porphyry- and (or) vein-type mineralization, under younger Tertiary volcanic cover. At a finer scale, magnetic patterns and fabrics permit discrimination between volcanic rocks of the Tertiary Chilcotin and Kamloops groups, and detection of subtle compositional and (or) structural variations within the groups. Contacts between volcanic cover and basement rocks and between basement units are more accurately defined, significantly reducing locally the areal extent of volcanic cover and opening up more ground for exploration. The high resolution of features in images of magnetic vertical derivatives reveals the Naver pluton to be more complex than currently mapped, comprising several integral elements, one of which may be a large roof pendant. Internal subdivisions of the Thuya batholith are defined, and annular marginal phases are proposed within two large granodioritic intrusions. Several new intrusions are proposed within the extensive, mainly sedimentary Devonian–Triassic terrain northeast of Kamloops, internal composition variation is suggested for some larger mapped intrusions, and areas underlain by some intrusions are enlarged.


Sign in / Sign up

Export Citation Format

Share Document