Research Article U-Pb and Hf Analyses of Detrital Zircons from Paleozoic and Cretaceous Strata on Vancouver Island, British Colu

GeoScienceWorld Lithosphere Volume 2021, Article ID 7866944, 20 pages https://doi.org/10.2113/2021/7866944

Research Article U-Pb and Hf Analyses of Detrital Zircons from Paleozoic and Cretaceous Strata on Vancouver Island, British Columbia: Constraints on the Paleozoic Tectonic Evolution of Southern Wrangellia

Daniel Alberts ,1 George E. Gehrels ,1 and Joanne Nelson 2

1Department of Geosciences, University of Arizona, Tucson, AZ 85721, USA 2British Columbia Geological Survey, Victoria, British Columbia V8W 9N3, Canada

Correspondence should be addressed to Daniel Alberts; [email protected]

Received 23 February 2020; Accepted 25 January 2021; Published 26 February 2021

Academic Editor: Sarah Roeske

Wrangellia is a late Paleozoic arc terrane that occupies two distinct coastal regions of western Canada and Alaska. The Skolai arc of northern Wrangellia in south-central Alaska and Yukon has been linked to the older, adjacent Alexander terrane by shared Late Devonian rift-related gabbros and also by Late Pennsylvanian postcollisional plutons. Late Devonian to Early Permian Sicker arc rocks of southern Wrangellia are exposed in uplifts on Vancouver Island, southwestern British Columbia, surrounded by younger strata and lacking physical connections to other terranes. Utilizing the detrital zircon record of Paleozoic and Cretaceous sedimentary rocks, we provide insight into the magmatic and depositional evolution of southern Wrangellia and its relationships to both northern Wrangellia and the Alexander terrane. 1422 U-Pb LA-ICPMS analyses from the Fourth Lake Formation (Mississippian–Permian) reveal syndepositional Carboniferous age peaks (344, 339, 336, 331, and 317 Ma), sourced from the Sicker arc of southern Wrangellia. These populations overlap in part known ages of volcanism, but the Middle Mississippian cumulative peak (337 Ma) documents a previously unrecognized magmatic episode. Paleozoic detrital zircons exhibit intermediate to juvenile ƐHf ðtÞ values between +15 and +5, indicating that southern Wrangellia was not strictly built on primitive oceanic crust, but instead on transitional crust with a small evolved component. The Fourth Lake samples yielded 49 grains (3.4% of the total grains analyzed) with ages between 2802 Ma and 442 Ma, and with corresponding ƐHf ðtÞ values ranging from +13 to -20. In age—ƐHf ðtÞ space, these grains fall within the Alexander terrane array. They were probably derived from sedimentary rocks in the basement of the Sicker arc. By analogy with northern Wrangellia, this basement incorporated rifted fragments of the Alexander terrane margin as the combined Sicker-Skolai arc system advanced ocean-ward due to slab rollback in Late Devonian to Early Mississippian time. Ultimately, data from detrital zircons preserved in the Fourth Lake Formation provides signiﬁcant information allowing for an updated tectonic model of Paleozoic Wrangellia.

1. Introduction Wrangellia and Karmutsen Formation of southern Wrangel- lia, that overlie Paleozoic arc-related sequences [1, 2]. Paleo- Wrangellia is one of the most outboard of the major northern zoic sequences of southern Wrangellia comprise the Upper Cordilleran terranes. It occupies two separate regions: south- Devonian to Lower Permian Sicker and Buttle Lake Groups central Alaska and southwestern Yukon (northern Wrangel- [3–7]. Those in northern Wrangellia comprise the Carbonif- lia) and Haida Gwaii and Vancouver Island (southern erous to Lower Permian Skolai Group, including the Station Wrangellia) (Figure 1 inset). It was originally deﬁned as a Creek and Hasen Creek formations [8–13]. coherent terrane based on characteristic thick piles of Until recently, the records of Paleozoic arc activity in Triassic ﬂood basalts, the Nikolai Greenstone of northern southern versus northern Wrangellia were considered to be

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 2 Lithosphere

Knight inlet

Bute inlet

129°W, 50°N

Coast Plutonic Complex

Tertiary Rocks Nanaimo Eocene Crescent terrane Dragon Property = Eocene Pacifc Rim terrane Bedingfeld Uplif = Buttle Lake Uplif = Cretaceous Nanaimo Group Cowichan Uplif =

Wrangellia Terrane Comox Samples Tis Study Jurassic Island Intrusions/ Fourth Lake Samples West Coast Crystalline Complex North Samples Saltspring Jurassic Bonanza Group Island South Samples Ruks, 2015 Triassic Karmutsen Formation Detrital Samples Victoria Devonian Sicker/ Samples CO1 and CO2 Matthews et al., [2017] Trust Fault Carboniferous Buttle Lake Groups 100 km 125°W, 48°N

Figure 1: Generalized geologic map of Vancouver Island (modiﬁed from Massey et al. [82, 83]). Inset map indicates relevant terranes of the northern Cordillera, location of study area, and general locations for three subdivisions of the Alexander terrane: SEM: Saint Elias Mountain region [17, 18]; SE: Alaska [65, 66, 71]; and BIM: Banks Island assemblage [65]. Abbreviations for displayed tectonic components on inset map: YA: Yakutat; CH: Chugach; CPC: Coast Plutonic Complex; PE: Peninsular; WR: Wrangellia; sWR: southern Wrangellia; nWR: northern Wrangellia; AX: Alexander. Sample locations for this study are denoted with squares. Colored circles are used to represent samples from Ruks [14]. Yellow hexagon represents sample location for samples from Matthews et al. [42]. Dashed ellipses mark the area of uplifts.

divergent, with Sicker volcanism confined to the Late Devo- to E-MORB basalts [7]. As such, the Sicker arc has been nian [7] and Station Creek volcanism poorly constrained as modeled as developing in an intraoceanic setting as a possible Carboniferous [11], with Pennsylvanian plutonism [12]. An extension of the Skolai arc, but with no direct connection to Early Mississippian U-Pb age from the base of the Station older terranes [15]. Creek Formation [13] and an extensive database of Devonian This study presents U-Pb ages and Hf isotope composi- to Permian U-Pb and microfossil ages from the Sicker Group tions of detrital zircons from upper Paleozoic strata of south- [14] now support consideration of the Sicker and Station ern Vancouver Island. These data are then used to refine Creek as parts of a single, evolving arc system, such as pro- interpretations on the tectonic evolution of southern Wran- posed by Beranek et al. [15], mainly based on data from gellia. First, the main grain populations will provide addi- northern Wrangellia. tional information on the nature and duration of Paleozoic The tectonic settings of northern and southern Wrangel- arc-related igneous activity. Second, minor grain populations lia are also distinct. Northern Wrangellia lies adjacent to the can be used to test for the possible presence of pre-Late Devo- older Alexander terrane and has been linked to it by Late nian basement to this part of Wrangellia. We analyzed nine Devonian gabbro complexes [13] and Pennsylvanian plu- samples from southern Vancouver Island, including five tonic suites that crosscut the terrane boundary [12, 15]. samples from the Fourth Lake Formation of the Buttle Lake Southern Wrangellia is isolated from other terranes by faults, Group and four samples from the Comox Formation at the by the Late Jurassic Coast Plutonic Complex, and by seaways, base of the Upper Cretaceous Nanaimo Group where it which render its Paleozoic tectonic context highly enigmatic. directly overlies Paleozoic strata (Figure 1). Nanaimo Group Moreover, the Sicker Group of Vancouver Island has been samples were analyzed primarily to gather additional infor- modeled as a Late Devonian nascent arc succession that mation about the Paleozoic history of Vancouver Island via developed on oceanic crust [4–7]. The lowest exposed rock primary and second-cycle zircons potentially preserved in unit, the lower Duck Lake Formation, is a pile of N-MORB these rocks.

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 Lithosphere 3

2. Geological Background Near the Duke River fault, the Station Creek Formation is absent, and the Hasen Creek nonconformably overlies the 2.1. Alexander Terrane. The late Paleozoic arc complex at the Late Devonian Steele Creek gabbro complex [13]. The Hasen base of northern Wrangellia has been linked to the adjacent Creek Formation represents postsubduction clastic sedimen- Alexander terrane in terms of probable basement and tec- tation following the Alexander-Wrangellia collision. tonic coevolution [13, 15]. This section begins with a descrip- tion of that older crustal fragment. The Alexander terrane 2.3. Paleozoic Stratigraphy and History of Southern extends over thousand kilometers from the St. Elias Moun- Wrangellia. The oldest rocks on Vancouver Island record tains on the Yukon-Alaska border, through southeast Alaska the early evolution of the southern Wrangellia arc. Included and along the north and central coast of British Columbia are Late Devonian–Carboniferous rocks of the Sicker Group ff (Figure 1 inset). It is a composite of two terranes with di er- and Mississippian–Permian strata of the Buttle Lake Group ent pre-Permian histories, the main Craig subterrane and the (Figures 1 and 3), which are exposed in the Buttle Lake and much smaller Admiralty subterrane of Admiralty Island, Cowichan uplifts of central and southern Vancouver Island southeast Alaska [15]. In northwestern Yukon/east-central (Figure 1) [4–6, 14, 20–22]. Alaska, Craig terrane volcanic, carbonate, and siliciclastic strata are subdivided into the Donjek (Upper Cambrian- 2.3.1. Sicker Group. Primarily, volcanogenic strata of the Middle Cambrian), Goatherd Mountain (Ordovician-Silu- Upper Devonian (to lowermost Mississippian?) Sicker Group fi rian), and Ice eld (Upper Silurian-Triassic) assemblages are exposed in the Cowichan and Buttle Lake uplifts – [16 18]. Detrital zircons from the Donjek assemblage show (Figure 1). Although stratigraphic nomenclature differs fl peaks at ca. 477 Ma, re ecting local arc/back-arc volcanism, between the two uplifts, overall similarities in the sections and ca. 565-760, 1000-1250, 1450, and 1650 Ma, from Balti- support correlations; therefore, they are combined in a single can cratonal sources and Timanian volcanic arcs of the east- column in Figure 2. fi ern Arctic region [17]. The lower Ice eld assemblage is of The oldest exposed unit, the Duck Lake Formation, is Late Silurian to Middle Devonian age and resembles Old restricted mainly to the southern part of the Cowichan uplift Red Sandstone successions of the Caledonides [18]. Consis- [4, 5, 7] It is absent in the Buttle Lake uplift. It comprises tent with this interpretation, main detrital zircon populations aphyric to plagioclase-phyric pillowed and massive basalts fl are 390-490 Ma, with Precambrian Hf model ages that re ect with minor tuff and chert and local felsic bodies near its progressive Silurian-Devonian orogenesis [15, 18]. The Craig top. Basalts of the lower Duck Lake Formation are tholeiites subterrane evolved into a passive margin environment from of E-MORB affinities, whereas the upper part of the Late Devonian to Early Pennsylvanian time [15]. formation comprises high-potassium calc-alkaline basalt and basaltic andesite with dacite dikes and felsic tuffs[4–7]. 2.2. Paleozoic Development of Northern Wrangellia and Its It represents the inception of arc volcanism in southern Interactions with Alexander Terrane. The oldest rocks in Wrangellia. A dacite flow near the top of the formation in northern Wrangellia belong to the Steele Creek gabbro com- its type area has yielded an LA-ICPMS age of 366:6±0:7 plex (Figure 2) along Duke River fault, which separates it Ma [14]. from the Alexander terrane [13]. Southwest of the fault, the The Nitinat Formation conformably overlies the Duck Mt. Constantine gabbro complex intrudes the Silurian- Lake Formation (Figure 2) [7]. It comprises augite- Devonian Bullion Creek limestone within the Alexander ter- plagioclase phyric calc-alkaline basalt to basaltic andesite rane [13] (Figure 2). Both complexes show N-MORB, nonarc breccias and flows, volcanic sandstone and siltstone, and geochemistry, and their U-Pb ages agree within error at ca. cherty tuffs[4–6]. The Price Formation in the Buttle Lake 363.5 Ma [13]. These data support the hypothesis that north- uplift, which comprises plagioclase-augite phyric andesite ern Wrangellia initiated as a Late Devonian arc that rifted flows and volcaniclastic deposits [23], is considered correla- away from the Alexander terrane due to slab rollback [13]. tive with the Nitinat Formation [7]. The Station Creek Formation records development of that Uppermost Sicker Group units are the predominantly arc-back-arc system, beginning in the Early Mississippian volcaniclastic McLaughlin Ridge and Myra formations, in according to a ca. 353 Ma U-Pb age of a felsic tuff in the lower the southern Cowichan and Buttle Lake uplifts, respectively part of the formation [13] (Figure 2). The Station Creek For- (Figure 2) [7, 23]. Igneous compositions range from mafic mation comprises a lower unit of mafic flows overlain by to felsic, with intermediate compositions most common [4– mafic to intermediate volcaniclastic strata [8] with arc type 7, 23]. Rhyolites are locally abundant, in some cases as well as BABB and N-MORB and E-MORB geochemical associated with volcanogenic massive sulphide deposits. signatures [13, 19]. A suite of Late Pennsylvanian plutons Myra Formation rhyolites have been dated as 366 ± 4 Ma intrudes both the Station Creek Formation and the Icefield [24] and 361:5±2:5 Ma (LA-ICPMS, Ruks [14]). McLaugh- Ranges assemblage of Alexander terrane, and one of them, lin Ridge rhyolites have yielded LA-ICPMS ages that span the Barnard Glacier pluton, cuts the terrane boundary [12, the Devonian-Mississippian boundary, between 363:0±6:7 15]. The Barnard Glacier suite is interpreted as the product and 353:1±3:4 Ma [14]. The coeval, cogenetic Saltspring of melting due to slab breakoff after Alexander terrane litho- Intrusive Suite yielded ages between 360:7±2:4 and sphere entered the subduction zone of the Skolai arc [15] The 355:0±1:5 Ma [7, 14, 25, 26]. Granitoids of the Saltspring Station Creek Formation passes upwards into the Lower Intrusive Suite are transitional to calc-alkaline, primarily Permian sedimentary Hasen Creek Formation (Figure 2). metaluminous, and yield trace element pattern characteristic

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 4 Lithosphere

250 252 Cowichan north (Alberni) Lopingian Cowichan south, Northern Wrangellia Northern Alexander (a) Buttle Lake (b) Vancouver Is. NW (c) (d) 260 259

Guadalupian LEGEND 270 St. Mary Lake St. Mary Lake 273 Limestone 280 PERMIAN F F F Hasen U F Mt. Mark Mt. Mark Sedimentary 290 U

Cisuralian Creek U Donjek suite Volcaniclastic 300 299 Upper F U U 307 F U U U F U Rhyolite, felsic tuf 310 Middle UpperBarnard Barnard U Glacier Glacier 315 U suite suite Flower Ridge Station Volcanic rocks Lower PENNSLVANIAN 320 Sicker GROUP SKOLAI (arc tholeiite, Creek 323 calc-alkalic)

Upper Group BUTTLE LAKE GROUP LAKE BUTTLE 330 331 (undivided) Lower Maﬁc volcanic rocks, Fourth Lk., U gabbro (N-MORB, Middle UU Station E-MORB, BABB) 340 U Telwood U

Creek assemblage Icefeld Felsic plutons 347 U 350 MISSISSIPPIAN U Mt. Lower U F U Constantine Diorite, granodiorite 359 complex 360 McLaughlin Ridge, U U-Pb age U Steele U Myra U U U U Creek Fossil Famennian u. Duck Lk.; Nitinat, PriceU UU F 370 Upper Duck Lake, Nitinat complex age

SICKER GP. Lower Duck Lake 372 F (ML Dev) Bullion DEVONIAN Frasnian Creek 380 limestone 383

Figure 2: Comparative stratigraphic columns from (a, b) southern Wrangellia on Vancouver Island, (c) northern Wrangellia, and (d) northern Alexander terrane. Vancouver Island data are from Juras [23], Massey [4–6], Yorath et al. [7], and Ruks [14]. Northern Wrangellia and Alexander data are from Israel et al. [13] and Beranek et al. [15]. Timescale is from Cohen et al. [84].

with rocks generated within an oceanic arc. However, there are enriched magmas during early phases of juvenile arc minor occurrences of peraluminous granitoids, which are construction [4–6, 14, 31]. more commonly associated with continental arcs [14]. McLaughlin Ridge and Myra formations represent a partially 2.3.2. Buttle Lake Group. The Sicker Group is overlain by the subaerial magmatic arc (Figure 2) [4–6]. mainly sedimentary lower Buttle Lake Group, including the The presence of volcanogenic massive sulphide occur- Fourth Lake Formation in the southern Cowichan uplift rences associated with small rhyolitic centers is consistent with and the Thelwood and Flower Ridge formations in the Buttle the extensional submarine arc to back-arc environments that Lake uplift (Figures 1–3) [7, 23]. host modern seabed massive sulphide deposits [27]. A unique The Thelwood Formation includes fine-grained, thin- occurrence of higher radiogenic lead within the lower accu- bedded siliceous and tuffaceous strata and penecontempora- mulations of VMS deposits in the Buttle Lake uplift indicates neous mafic sills. It is overlain by amygdaloidal plagioclase- a contribution of lead from a more evolved source [28–30]. pyroxene phyric basalt lapilli tuff, breccia, and flows of the Along with higher radiogenic lead levels in early VMS Flower Ridge Formation [23], which is correlated with deposits and granitoids with peraluminous occurrences, a basalts in the Fourth Lake Formation [6]. third piece of evidence from early Wrangellia indicative of The Fourth Lake Formation, the target unit for this study Ɛ fi an evolved component are Ndðt=360Þ values of ma c to felsic (Figures 2 and 3), consists of a 100-200-meter-thick sequence rocks of the Sicker Group yielding values from +5.9 to +4.5 of radiolarian ribbon chert that is overlain by cherty siltstone, [14]. These values are less radiogenic than the depleted thin siltstone-argillite beds, and thinly bedded fine- to mantle reservoir at 360 Ma of +9.2 and are interpreted to be medium-grained interlayered sandstone and mudstone isotopically juvenile to intermediate. This led to the interpre- (Figures 4(a), 4(b), and 4(e)–4(g)). The Fourth Lake Forma- tation that sediment with an evolved isotopic composition tion represents a marginal-basin assemblage that developed

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 Lithosphere 5

66 Nanaimo Group 17AVI06 17AVI07A 17AVI07B Cretaceous 17AVI09

145

Jurassic Bonanza Volcanics/Island Intrusions /West Coast Crystalline Complex 201

Quatsino and Parson Bay Fms Triassic Karmutsen Flood Basalts

252

Permian St. Mary’s Lake Fm Mount Mark Fm 299 Pennsylvanian Buttle Lake 17AVI08 323 Group 17AVI05 Fourth Lake Fm 17AVI04 Mississippian 17AVI03 17AVI02 359 McLaughlin Ridge Fm Sicker Nitinat Fm Group Duck Lake Fm Devonian

419

Intermediate to Sandstone Felsic Volcanic rock

Limestone Mafc Volcanic rock

Interbedded fne-grained Saltspring Intrusion clastic rock

Chert Basalt

Figure 3: Generalized stratigraphic column and interpreted sample positions (based on maximum peak ages, illustrated for stratigraphic context) for strata from Vancouver Island, B.C. (adapted from Massey and Friday, [20]). Stratigraphic nomenclature is from Massey [5]. Timescale is from Cohen et al. [84]. Stratigraphic column ages are displayed in millions of years.

with coeval VMS-type deposits in the back-arc of southern portion of the Cowichan uplift, within the Fourth Lake For- Wrangellia [4–6, 14]. mation, yields three distinct populations with age ranges of The Fourth Lake Formation has yielded Early to mid- 339-337 Ma, 322-309 Ma, and ca. 295 Ma, the latter of which Tournaisian and Middle to Late Pennsylvanian conodonts is interpreted to represent the depositional age of the tuﬀ [14]. [7] and one collection of Early Permian radiolaria [14]. Overlaying the Fourth Lake Formation is the Mount Previous detrital zircon U-Pb ages from the Fourth Lake For- Mark Formation, a massive limestone unit rich in marine mation yield grains that range in age from 355 Ma to 294 Ma fossils (Figures 2 and 3) [32]. The contact between the Mt. with dominant peak ages of 320, 312, and 304 Ma (Figures 1 Mark and underlying strata is diachronous: it contains faunal and 5) [14]. A heterolithic lapilli tuﬀ from the northern assemblages of Middle Pennsylvanian and Early Permian age

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 6 Lithosphere

(a)

(b)

(e) (f) (g)

Figure 4: Select outcrop and hand sample photos from the (a, b, e–f) Fourth Lake Formation and the (c, d) Comox Formation. Further sample descriptions are displayed in Table 1. (a) 17AVI02 hand sample displaying interbedding of chert and medium-grained sandstone. (b) Outcrop of sample 17AVI05 consisting of chert and interbedded layers of silt to fine-grained sandstone. (c) Large pebble clast of interbedded chert and sandstone, indicative of the Fourth Lake Formation, within the Comox Formation (sample 17AVI06). (d) Smaller pebble clast of interbedded chert and sandstone, indicative of the Fourth Lake Formation, within the Comox Formation (sample 17AVI07A). (e) Outcrop photo from sample 17AVI03 displaying layering of chert and interbedded fine-grained sandstone. (f) Close-up outcrop photo of sample 17AVI05, a layer of silt to fine sandstone at the centimeter scale. (g) Hand sample photo of sample 17AVI08, centimeter scale beds of fine- to medium-grained sandstone interbedded with chert, with a distinctive black color.

[7], overlapping those in the Fourth Lake Formation 2.3.3. Mississippian-Early Permian Sicker Group Strata, (Figures 2 and 3). Conformably overlying the Mount Mark Alberni (North Cowichan); Northwest Vancouver Island. Formation is a succession of clastic and volcaniclastic rocks Extensive U-Pb dating of igneous rocks at the northern end of the St. Mary’s Lake Formation (Figures 2 and 3), inter- of the Cowichan uplift near Port Alberni and on northwestern preted to be Early Permian in age [6, 21]. Vancouver Island (Bedingﬁeld uplift and Dragon property,

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 Lithosphere 7

Location on VI 317 17AVI08 (n = 307) 49.07, -124.38 (8xVE) 331 17AVI05 (n = 109) (8xVE) 48.87, -124.07

336 17AVI04 (n = 313) 48.90, -124.18 (8xVE) 339 17AVI02 (n = 313) 49.00, -124.42 344 17AVI03 (n = 380) 48.90, -124.19 (8xVE) Fourth Lake Formation 337 (n = 1422) (8xVE) 312 08TR019 (Ruks, 2015) n = 61 49.17, -124.09

320 08TR017 (Ruks, 2015) 304 n 353 = 26 49.17, -124.10 305 Ruks, (2015) North (N = 11, n = 180) 358 Ruks, (2015) South (N = 32, n = 565) 358 Ruks, (2015) 309 (N = 44, n = 832) 250 300 350 400 800 1200 1600 2000 2400 2800 Detrital zircon age (Ma)

Figure 5: Normalized age distribution diagram for detrital zircons from the Fourth Lake Formation and igneous zircon U-Pb ages from Ruks [14] separated into two curves based on geographic location as expressed in Figure 1. Two detrital samples (08TR017 and 08TR019) from Ruks [14] are separated from igneous samples in order to compare with detrital samples in this study. Main peaks are noted for each sample in millions of years. Number of samples within Ruks [14] distributions denoted with “N,” number of analysis in each curve denoted with “n.” Proportions of >400 Ma ages have been vertically exaggerated by a factor of eight relative to <400 Ma grains. Sample locations are generally placed on a map of Vancouver Island; latitude and longitude coordinates are displayed to the right for individual samples and in Table 1.

Figure 1; Ruks [14]) has shown that mainly volcaniclastic strata 293 Ma and 291 Ma (Figure 1) [14]. Nd isotope data for igne- assigned in previous mapping [4] to the Sicker Group are ous rocks of this age yield values ranging from +7.0 to +4.3 at largely time-equivalent with the sedimentary Fourth Lake For- an age of ca. 300 Ma. A single sample with an age of ca. mation in the classic South Cowichan and Buttle Lake Paleo- 317 Ma yields a ƐNdðtÞ value of +6.1 [14]. These values are zoic sections (Figure 2). Here, we assign these younger strata marginally more juvenile than rocks of the Sicker Group. to a newly recognized informal upper part of the Sicker Group. In all of these areas, the volcanic strata are directly overlain by limestones of the Mt. Mark Formation (Figure 2) [4, 14]. 2.3.4. Arc Migration in the Southern Wrangellia Arc System. Near Port Alberni, in the northernmost Cowichan uplift, An apparent northward migration of arc magmatism is felsic volcanic and tuffaceous rocks yield detrital zircon LA- observed from the pattern of U-Pb ages and geographic posi- ICPMS ages from 351:3±2:1 Ma to 334:9±2:2 Ma, and rhy- tions of igneous rock units on southern Vancouver Island olite porphyry bodies yield LA-ICPMS ages from 352:8±3:1 [14]. Igneous zircon U-Pb ages from rocks of the Sicker Ma to 335:0±2:9 Ma, with a single result of 293:6±1:6 Ma Group in the Cowichan uplift peak at an age of 358 Ma, (Figure 1) [14]. Felsic volcanic rocks associated with VMS whereas rocks of the Buttle Lake Group in the Bedingfield deposits in the Bedingfield uplift yield igneous zircon LA- uplift and Dragon property peak at an age of 305 Ma ICPMS ages of 312:4±3:5 Ma, 308:3±3:2 Ma, and 305:6± (Figures 1 and 5). A considerable lack of ages ranging from 4:1 Ma (Figure 1) [14]. Similar VMS-associated felsic rocks 334 to 314 Ma between southern and northern uplifts on in the Dragon property yield U-Pb zircon ages from 311 Ma Vancouver Island has been interpreted as a “magmatic gap” to 300 Ma and are overlain by felsic tuffs with ages of ca. in the arc system of southern Wrangellia [14].

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 8 Lithosphere

2.4. Mesozoic Rocks. Paleozoic rocks in northern and southern Vancouver Island (Figures 4(c) and 4(d)). The Comox portions of Wrangellia are overlain by ~4-7 km of the Triassic Formation is interpreted to consist of fluvial sediments Nikolai-Karmutsen flood basalts [1]. A younger, latest trapped by localized topographic highs during deposition, Triassic-Jurassic arc is represented by stratified volcanic- as well as deltaic deposits containing large incised fluvial sedimentary rocks of the Bonanza Group and associated channels [43, 44, 53, 57]. Island Intrusions in southern Wrangellia [21, 33–36]. Upper Previous samples analyzed by LA-ICPMS techniques Cretaceous sediments of the Nanaimo Group accumulated from the Comox Formation (samples CO1 and CO2 of on and adjacent to the older stratigraphic units of southern Matthews et al. [42]) from the east-central coast of Wrangellia in a convergent margin basin west of the Coast Vancouver Island (Figure 1) yield peak ages of 154-153 Ma Mountain Batholith [37–44]. and 94-91 Ma (Figure 6). Sample CO2 also yields 50 grains that range in age from 204 Ma to 163 Ma, with a peak age 2.4.1. Karmutsen Formation. Paleozoic rocks on Vancouver of 167 Ma (Figure 6). A source from the central Coast Moun- Island are overlain unconformably by ca. 3.5 km of ca. 230- tain Batholith is the preferred provenance interpretation by 225 Ma submarine flows and pillow basalts of the Karmutsen Matthews et al. [42] for zircon grains with ages from 150 to Formation (Figures 1 and 3) [2, 24, 45]. Gabbroic rocks 80 Ma. This interpretation is consistent with previous studies related to Karmutsen basalts yield U-Pb zircon ages of ca. indicating that Late Jurassic to mid-Cretaceous plutons of the 228-226 Ma from outcrops on Saltspring Island [26]. Coast Mountain Batholith contributed significant sediment to the Nanaimo Group [39, 40, 58, 59]. 2.4.2. Bonanza Group. The youngest widespread igneous assemblage on Vancouver Island is latest Triassic-Jurassic 3. Methods in age and consists of arc-type rocks of the Bonanza Group, Island Intrusions, and West Coast Crystalline Complex Fieldwork focused primarily on collecting samples for geochro- (WCC) (Figures 1 and 3) [33, 36]. The WCC contains gab- nologic analysis. The sampling strategy was targeted at precise broic to dioritic plutonic rocks, migmatites, amphibolites, localities of the Fourth Lake Formation containing fine- to and metasedimentary rocks and preserves the deepest section medium-grained sandstone fractions, guided by descriptions of the Bonanza arc. Rocks of the WCC have yielded U-Pb zir- and coordinates of Massey [5, 6]. Recent forest clear cuts pro- con ages ranging from ca.190 to 177 Ma [33, 35], a Rb-Sr age vided additional exposure of previously unknown outcrops of of 151 Ma [33], and K-Ar ages of 172 and 163 Ma [33]. The the Fourth Lake Formation (Figure 4(b)). Isolating millimeter midcrustal section is preserved in batholiths and felsic intru- to centimeter thick sandstone layers within the Fourth Lake sions of the Island Intrusive Suite (Figure 1) [35, 36], which Formation proved difficult considering the interbedded nature have yielded U-Pb zircon crystallization ages from ca.175 to of these samples, which lead to the collection and processing of 168 Ma, 40Ar-39Ar cooling ages of ca.176 and 166 Ma, and a entire chunks of interbedded sandstone, mudstone, and/or K-Ar age range of 181-152 Ma [46, 47]. The Bonanza Group chert (Figures 4(a) and 4(e)–4(g)). During collection of the consists of ~2,500 meters of lava flows, pyroclastic flows, thin Fourth Lake Formation, in the southern Cowichan uplift, field interbedded sedimentary units, and minor low-grade meta- relations of the Comox Formation unconformably overlaying morphic rock that are interpreted to be the volcanic equiva- the Fourth Lake Formation and observations of clasts resem- lents of the Island Intrusions (Figure 1) [34, 36, 48–50]. bling the Fourth Lake Formation prompted sample collection These rocks have yielded a U-Pb zircon age range from ca. (Figures 4(c) and 4(d)). These observations lead to the hypoth- 202 Ma to 165 Ma [35, 51]. esis that the Comox Formation was likely to yield detrital zir- Whole-rock geochemical data from volcanic rocks of the con grains that would provide insights into the Paleozoic and Bonanza Group indicate that melts were primarily mantle Early Mesozoic evolution of southern Wrangellia. derived, but that moderately intermediate to juvenile geo- Zircon extraction was performed at the Arizona Laser- chemical signature record varying contributions from older Chron Center (http://www.laserchron.org) using methods rocks on Vancouver Island (e.g., Sicker Group, Buttle Lake described by Gehrels et al. [60], Gehrels and Pecha, [61], Group, and Karmutsen flood basalts) [36]. and Pullen et al. [62]. Primary steps included crushing/pul- The equivalent latest Triassic-Jurassic arc sequence in verizing, usage of a Wilfley table, Frantz magnetic separator, southern Alaska is referred to as the Talkeetna arc system [52]. and heavy liquids. Grains from each sample were poured in a 1-inch epoxy mount alongside fragments of U-Pb zircon 2.4.3. Nanaimo Group. Paleozoic through Jurassic rocks on standards (FC-1, SL2, and R33) and Hf zircon standards Vancouver Island are overlain unconformably by ca. 5 km (Mud Tank, Temora-2, FC-1, 91500, Plesovice, R33, and of Upper Cretaceous nonmarine to deep-marine basinal SL2). Mounts were polished with fine sandpaper and finished strata of the Nanaimo Group (Figures 1 and 3) [38–42, 44, with a 1 μm diamond polish. All sample mounts were imaged 53]. The basal unit of the Nanaimo Group is the Comox For- using cathodoluminesence (CL) and backscatter electron mation, which contains various macrofossil assemblages that (BSE) methods. Samples were cleaned with a 2% HNO3 indicate a Turonian to Coniacian depositional age [54–57]. A and 1% HCL solution prior to isotopic analysis. CL and varying thickness of 0-350 m is reported for the Comox For- BSE images were utilized to select analytical points, avoiding mation, which consists primarily of poorly bedded pebble to complex internal structures and nonzircon grains. boulder conglomerate. Clasts are angular to rounded, and U-Pb analyses were conducted by laser ablation induc- compositions indicate derivation from older strata on tively coupled plasma mass spectrometry (LA-ICPMS) using

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 Lithosphere 9

20 Fourth Lake Formation DM (n = 1422) 15

5 337 Ma CHUR Crustal evolution

0 Epsilon Hf Epsilon –5

–10

–15 Crustal evolution (x5 VVE) –20 275 300 325 350 375 400 800 1200 1600 2000 2400 2800 Detrital zircon age (Ma)

Ruks, (2015) 17AVI03 Nd converted 17AVI08 17AVI04 17AVI05 17AVI02

Figure 6: U-Pb and ƐHf ðtÞ values from the Fourth Lake Formation. Converted ƐNDðtÞ values from Ruks [14] are displayed as triangles. Lower curve is a cumulative probability curve for the Fourth Lake Formation samples from this study. Ages > 400 Ma are vertically exaggerated by a factor of ﬁve relative to <400 Ma grains. In order to compare ƐNdðtÞ and ƐHf ðtÞ values, a conversion to ƐHf ðtÞ is used (ƐHf ðtÞ =1:36 × ƐNdðtÞ +2:95), as described by Vervoort et al. [85]. Upper plot shows ƐHf ðtÞ values for individual samples. Reference lines on the Hf plot are as follows: DM—depleted mantle, calculated using 176Hf/177Hf0 = 0:283225 and 176Lu/177Hf0 = 0:038512 [86]; CHUR—chondritic uniform reservoir, calculated using 176Hf/177Hf = 0:282785 and 176Lu/177Hf = 0:0336 [87]. Black arrows show interpreted crustal evolution trajectories assuming present-day 176Lu/177Hf = 0:0093 [64, 85, 88].

a Photon Machines G2 excimer laser connected to an Ele- 4. U-Pb Geochronologic and Hf Isotope ment2 single-collector ICPMS. One spot was analyzed per Geochemical Results grain using a laser beam diameter of 15 or 20 microns. The number of analyses per mount ranged from 109 to 380, 4.1. Fourth Lake Formation. Sample 17AVI02, an interbed- depending on the abundance of zircon extracted from each ded chert and ﬁne-grained sandstone (Figure 4(a)), yields sample. Results of U-Pb zircon analyses are reported in DR 313 U-Pb ages with a single group ranging from 365 to Table 3. 310 Ma and a peak age of 339 Ma (Figure 5). ƐHf ðtÞ values Lu-Hf analyses were conducted using a Photon from 36 analyses yield a dominantly juvenile signature, rang- Machines G2 excimer laser connected to an NU Plasma ing from +17 to +5, with two analyses yielding more evolved HR multicollector ICPMS, utilizing analytical methods values of -15 and -11 (Figure 6). described by Cecil et al. [63] and Gehrels and Pecha [61]. Sample 17AVI03 (see Figure 4 caption and Table 1 for Hf analyses were conducted with a 40-micron spot cen- sample details) yields 380 U-Pb ages with a dominant range tered over the preexisting pit evacuated by U-Pb analyses. of ages from 371 to 313 Ma and a peak age of 344 Ma Hf analyses were conducted on grains representative of (Figure 5). This sample also contains older grains displaying each age group in each sample. The average uncertainty minor age peaks of 1202, 1015, 691, 668, 597, 572, 536, 442, for analyses is 1.4 epsilon units (2σ). External precision, and 402 Ma (Figure 5). 17AVI03 also yielded several single ages based on analysis of standards from the same mounts as of 2679, 2133, 1665, and 1465 Ma (Figure 5). Grains represen- unknowns, is 1.2 epsilon units (2σ). Analytical parameters tative of the 344 Ma peak age yield predominantly juvenile Ɛ are reported in DR Table 2. To assist in interpretations of Hf data, ƐHf ðtÞ values (at the time of crystallization) Hf ðtÞ values, ranging from +16 to +5, whereas a single grain Ɛ Ɛ within 5 units of depleted mantle (DM) are referred to as yields an evolved Hf ðtÞ value of -17. Hf ðtÞ values for eight juvenile in composition, values between 5 and 12 units zircon grains that fall within an age range of 714 to 442 Ma below DM are considered moderately juvenile are intermediate to evolved, with values ranging from +2 to (intermediate), and values >12 units below DM are -19. Six zircon grains older than 1000 Ma yield intermediate considered evolved (following Bahlburg et al. [64]). to evolved ƐHf ðtÞ values ranging from +6 to -7 (Figure 6).

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 10 Lithosphere

Table 1: Sample descriptions, peak ages, # of analyses done, and locations from the Fourth Lake Formation and Comox Formation on southern Vancouver Island.

# of U-Pb #ofHf Sample Unit Rock type Prominent U − Pb age peaks < 400 Ma Lat/long analysis analysis Fourth Lake Interbedded chert and fine- 49.006, 17AVI02 339 Ma 313 36 Formation grained sandstone -124.42 Fourth Lake Interbedded chert and fine- 48.906, 17AVI03 344 Ma 380 58 Formation grained sandstone -124.1937 Fourth Lake Cobbles with interbedded 48.908, 17AVI04 336 Ma 313 36 Formation mudstone and sandstone -124.1937 Fourth Lake Sandstone-siltstone argillite 48.873, 17AVI05 331 Ma 109 32 Formation interbedded with chert -124.072 Fourth Lake Interbedded chert and fine- to 49.0672, 17AVI08 317 Ma 307 41 Formation coarse-grained sandstone -124.3869 Comox 48.873, 17AVI06 Pebble conglomerate 201 Ma, 342 Ma 118 26 Formation -124.077 Comox 49.044, 17AVI07A Pebble to cobble conglomerate 166 Ma, 194 Ma 314 33 Formation -124.431 Comox 49.043, 17AVI07B Pebble to cobble conglomerate 167 Ma, 197 Ma 314 39 Formation -124.426 Comox 49.2067, 17AVI09 Pebble to cobble conglomerate 87 Ma, 152 Ma, 165 Ma, 356 Ma 312 142 Formation -124.035

Sample 17AVI04, an interbedded mudstone and sand- of +17 to +9. Grains older than 371 Ma comprise 3.4% of the stone cobble collected near sample 17AVI03, yields 313 U- total grains analyzed. ƐHf ðtÞ values of these older grains Pb ages with a dominant range of 364-316 Ma and a peak range from intermediate to highly evolved, with a range of age of 336 Ma (Figure 5). This sample yields four older grains +5 to -20 (Figure 6). with ages of 2802, 941, 492, and 374 Ma (Figure 5). Grains within the main age group yield juvenile to intermediate Ɛ 4.2. Comox Formation. Sample 17AVI06 was collected near Hf ðtÞ values of +15 to +6. Older grains, from oldest to youn- the unconformable contact between the Fourth Lake and Ɛ gest, yield Hf ðtÞ values of -1, 0, +4, and +15 (Figure 6). Comox Formations and contains large pebble-sized clasts of Sample 17AVI05 (see Figure 4 caption and Table 1 for the Fourth Lake Formation within a medium-grained sandy sample details) yields 109 U-Pb ages with a dominant group matrix (Figure 4(c)). Processing of entire conglomerate of 359-311 Ma and a peak age of 331 Ma (Figure 5). Two chunks yielded 118 detrital zircon U-Pb ages with prominent grains older than 359 Ma yield ages of 2749 Ma and peak ages at 341 Ma and 202 Ma and subordinate peak ages 1612 Ma (Figure 5). ƐHf ðtÞ values for zircon grains from of 196 and 159 Ma (Figure 7). Single-grain ages are 263, 223, the main cluster yield primarily juvenile values ranging from 128, and 87 Ma. Grains within the main age groups yield juve- Ɛ +15 to +5. The two older grains yield ƐHf ðtÞ values of +7 and nile Hf ðtÞ values ranging from +15 to +6 (Figures 8 and 9). +3, respectively (Figure 6). Three Hf analyses were conducted on grains within the subor- Sample 17AVI08 was collected furthest to the north in dinate younger age group and yield ƐHf ðtÞ values ranging the southern Cowichan uplift and consisted of interbedded from +13 to +9. Two single grains with ages of 263 Ma and chert and fine-grained sandstone with a distinctive black 87 Ma yield ƐHf ðtÞ values of +10 and +13, respectively color (Figure 4(g)). This sample yields 307 U-Pb ages with (Figure 9). a prominent group of 361-281 Ma and a peak age of Sample 17AVI07A, exposed due to recent clear cuts in 317 Ma. There are three single-grain ages of 1025, 273, and the southern Cowichan uplift, had relatively smaller pebble- 270 Ma (Figure 5). Hf isotope results from 41 zircon crystals sized clasts of the Fourth Lake Formation within a finer Ɛ of 361-281 Ma yield juvenile to intermediate Hf ðtÞ values of sandy matrix (Figure 4(d)). This sample yielded 314 U-Pb +17 to +2. The single older grain of 1025 Ma yields a ƐHf ðtÞ ages that belong to two different groups. The older group value of +5 (Figure 6). ranges from 217 to 184 Ma, with a peak age of 194 Ma, As shown in Figure 5, our five samples from the Fourth whereas the younger group yields an age range of 180- Lake Formation yield peak ages of 344, 339, 336, 331, and 159 Ma, with a peak age of 166 Ma (Figure 7). Two single 317 Ma, that young northward, and produce a combined grains yielded ages of 363 Ma and 360 Ma. ƐHf ðtÞ values from peak age of 337 Ma. The Hf isotopic compositions of these 33 zircon grains of 217-159 Ma range from +13 to +5, Paleozoic zircon grains record derivation from juvenile to whereas the two older grains observed in this sample yield intermediate signatures, with most ƐHf ðtÞ values in the range ƐHf ðtÞ values of +13 and +10 (Figures 8 and 9).

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 Lithosphere 11

94 Matthews et al., 2017 153 167 CO2 (n = 207)

91 Matthews et al., 2017 154 CO1 (n = 248)

367 (10x VE) 07JBM06 (n = 460) 155 196 (10x VE) 152 17AVI09 87 (n = 312) 165 356

167 197 17AVI07B (n = 314)

166 17AVI07A (n = 314)

194

341 202 17AVI06 (n = 118)

167 (10x VE) 17AVI Comox 195 (n = 1055) 86 152 341

166 (10x VE) Comox Formation 91 153 195 n 365 ( = 1970) 0 50 100 150 200 250 300 350 400 450 500 550 600 .9 1.2 1.5 1.8 2.1 2.4 2.7 3.0 Detrital zircon age (Ma)

Figure 7: Normalized age distribution diagram for detrital zircons from the Comox Formation compared to Comox Formation samples from Matthews et al. [42] (CO1 and CO2). Main peaks are noted for each sample in millions of years. Proportions of >600 Ma grains have been vertically exaggerated by a factor of ten relative to <600 Ma grains (VE). Number of analyses in each sample denoted with “n.”

Sample 17AVI07B, collected at the top of the outcrop several single-grain ages between 308 Ma and 09 Ma, and where sample 17AVI07A was taken, yields 314 U-Pb ages four single grains of 1712, 529, 396, and 381 Ma (Figure 7). with two equally prominent age peaks of 197 Ma and Of the 142 grains analyzed for Hf isotopic data, ca. 152 Ma 167 Ma and a minor peak age of 83 Ma (Figure 7). A scatter- grains yield ƐHf ðtÞ values ranging from +16 to +6 (with a sin- ing of single grains yields older ages of 365, 356, 349, and gle grain yielding an evolved value of -11), ca. 165 Ma grains 321 Ma (Figure 7). 39 ƐHf ðtÞ values from sample 17AVI07B yield ƐHf ðtÞ values of +14 to +6 (with a single grain yielding a were generated, with ƐHf values for the two prominent (t) value of -6), and ca. 87 Ma grains yield ƐHf ðtÞ values ranging age groups yielding ƐHf ðtÞ values ranging from +14 to +7 from +13 to +8 (Figures 8 and 9). Zircons within the age (Figures 8 and 9). ƐHf ðtÞ analyses of the younger ca. 83 Ma range of 360-340 Ma yield ƐHf ðtÞ values ranging from +11 grains yield values of +7 (Figures 8 and 9). Older ages, to +5 (Figures 8 and 9), whereas 308-209 Ma grains have Ɛ encompassing grains from 365 to 321 Ma, yield ƐHf ðtÞ values Hf ðtÞ values ranging from +15 to +3 (Figures 8 and 9). Single ranging from +13 to +9 (Figures 8 and 9). Ɛ Sample 17AVI09, a poorly sorted pebble to cobble con- grains with ages of 529 Ma and 396 Ma yield positive Hf ðtÞ glomerate was sampled from a road cut near the city of values of +1 and+9, respectively (Figures 8 and 9). Nanaimo (Figure 1 and Table 1), yields 312 U-Pb ages with A single sample (07JBM06) of the Comox Formation dominant peak ages of 152 Ma and 87 Ma (Figure 7). A outcropping on Saltspring Island (Figure 1) was collected minor peak of 165 Ma is observed on the old side of the dom- by Brian Mahoney and is included in this study. This sample inant 152 Ma age group (Figure 7). Older ages include a sub- yields dominant age peaks of 367, 196, and 155 Ma (Figure 7). ordinate group of 360-340 Ma with a peak age of 356 Ma, Six individual grain ages older than 1000 Ma are 2678, 1828,

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 12 Lithosphere

20 DM 15

CHUR

0 Epsilon Hf Epsilon Crustal evolution –5

–10

–15

–20 167 Comox Formation (n = 1055)

195

86 152 341

337 Fourth Lake Formation (n = 1422)

0 100 200 300 400 500 600 Detrital zircon age (Ma) Fourth Lake Formation Comox Formation

Figure 8: U-Pb and ƐHf ðtÞ values from the Fourth Lake Formation and the Comox Formation analyzed in this study. Lower curves are cumulative normalized probability curves for each formation. Upper plot shows ƐHf ðtÞ values for all samples analyzed. Reference lines on the Hf plot are as follows: DM—depleted mantle, calculated using 176Hf/177Hf0 = 0:283225 and 176Lu/177Hf0 = 0:038512 [85]; CHUR—chondritic uniform reservoir, calculated using 176Hf/177Hf = 0:282785 and 176Lu/177Hf = 0:0336 [87]. Black arrows show interpreted crustal evolution trajectories assuming present-day 176Lu/177Hf = 0:0093 [64, 85, 88].

1134, 1072, 1040, and 1037 Ma (Figure 7). Younger grains tion (ﬁve from this study and two from Ruks [14]) show con- that are less than 110 Ma yield ages of 108 and 87, with two siderable overlap with the igneous U-Pb ages from rocks of ages at 81 Ma (Figure 7). Hf isotope data are not available the Sicker Group and Buttle Lake Group. Additionally, detri- from this sample. tal zircon peak ages and cumulative peak ages of igneous zir- Collectively, Comox samples yield peak ages of 365, 195, con ages are apparently young northward. Maximum peak 166, 153, and 91 Ma (Figure 7). Most ƐHf ðtÞ values from ages in the range of ca. 360-340 Ma are associated with Sicker these grains range from +15 to +5 (Figures 8 and 9). Group magmatism, predominately observed in the Cowi- chan uplift, whereas 320-300 Ma peak ages are sourced from 5. Provenance of the Fourth Lake Formation Buttle Lake Group igneous rocks commonly found in the Bedingﬁeld uplift and Dragon property to the north 5.1. Paleozoic Zircons. As shown in Figure 5, most of the (Figure 1). This younging trend is most likely due to the U-Pb ages of detrital zircons from the Fourth Lake Forma- progressive rifting of the arc shifting magmatic centers

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 Lithosphere 13

20 DM 15

Crustal evolution CHUR 0

Epsilon Hf

–5 166

195 –10 91 153 Comox Formation 365 n = 1970 –15 Central Coast Mountain Batholith n = 3015

Southern Coast Mountain Batholith n = 2822

0 50 100 150 200 250 300 350 400 Age (Ma)

Central CMB Magmatic fux curve 17AVI09 Age-distribution 17AVI07B Magmatic fux curve 17AVI07A Age-distribution 17AVI06

Figure 9: U-Pb and ƐHf ðtÞ values from the Comox Formation and central Coast Mountain Batholith ƐHf ðtÞ values are from Cecil et al. [63]. Filled grey curve is the cumulative age distribution for the Comox Formation, which includes samples from Matthews et al. [42]. The pink age-distribution curve for the Central Coast Mountains is from Gehrels et al. [68] and Cecil et al. [63]. Filled purple age-distribution curve for the Southern Coast Mountains is from Cecil et al. [89]. Dashed magmatic ﬂux curves for the Central Coast Mountains [68] and Southern Coast Mountain Batholith [89] are represented with dashed lines. Reference lines on the Hf plot are as follows: DM—depleted mantle, calculated using 176Hf/177Hf0 = 0:283225 and 176Lu/177Hf0 = 0:038512 [86]; CHUR—chondritic uniform reservoir, calculated using 176Hf/177Hf = 0:282785 and 176Lu/177Hf = 0:0336 [87]. Black arrows show interpreted crustal evolution trajectories assuming present-day 176Lu/177Hf = 0:0093 [64, 85, 88].

northward along with northerly propagation of depositional of Sicker Group and Buttle Lake Group aged zircons, these centers. ages must have been derived from igneous rocks of southern One of the main diﬀerences between Paleozoic igneous Wrangellia. This suggests that southern Wrangellia was highly versus Paleozoic detrital records of southern Wrangellia is active throughout the Carboniferous, contradicting the previ- the abundance of detrital zircon ages within the 334-314 Ma ously interpreted magmatic gap of 334–314 Ma. A proposed gap observed in igneous ages reported by Ruks [14]. Three location for these igneous bodies would be somewhere samples contain dominant peak ages of 331 Ma (sample between the Buttle Lake and Cowichan uplifts based on the 17AVI05), 320 Ma (08TR017, from Ruks [14]), and 317 Ma perceived trajectory of arc migration to the north as proposed (17AVI08) (Figure 5), and all of these samples were collected above (Figure 1). However, these rocks are most likely covered further north in the Cowichan uplift and near the city of beneath the widespread Karmutsen basalts (Figure 1). Nanaimo (Figure 1). Considering the low-energy depositional The Hf isotope data acquired in this study and the Nd environment of the Fourth Lake Formation and the presence isotope data from Ruks [14] can also be used to compare

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 14 Lithosphere

the detrital and igneous records. Converted ƐNdðtÞ to ƐHf ðtÞ in these strata were derived directly from the Sicker and But- values from Ruks [14] yield values for ca. 360 Ma samples tle Lake Groups or recycled from the Fourth Lake Formation. ranging from +11 to +9, +11 for a single rock with an age The latter interpretation is preferred for detrital zircons ana- of ca. 317 Ma, and +12 to +8.7 for rocks with ages of ca. lyzed in this study based on the occurrence of clasts derived 300 Ma (Figure 6). As shown on Figure 6, ƐHf ðtÞ values for from the Fourth Lake Formation observed in these samples detrital zircons from the Fourth Lake Formation yield a sim- (Figures 4(c) and 4(d)). Sample 07JBM06 yields the oldest ilar age range of +17 to +10 (Figures 5 and 6). In addition, Paleozoic peak age at 367 Ma, which is consistent with deri- both igneous and detrital zircons are interpreted to show a vation from nearby rocks of the Sicker Group and the Salt- spring Intrusive Suite, which is of similar age and decrease in the abundance of intermediate ƐHf ðtÞ values representative of first cycle deposition of zircons into the from 325 Na to 290 Ma (Figure 6). The occurrence of similar Comox Formation (Figures 1 and 7). Conglomeratic samples U-Pb ages and Hf isotope values in both data sets supports of the Comox Formation located southeast of Vancouver the interpretation that strata of the Fourth Lake Formation Island in the San Juan Islands similarly yield several Sicker were sourced from the Sicker Group and igneous rocks of Group aged detrital zircons and a minor amount of 400- the Buttle Lake Group (Figure 6). 500 Ma ages [53]. As shown in Figure 8, Paleozoic zircons from the Comox Formation yield ƐHf ðtÞ values of +15 to 5.2. Neoproterozoic to Paleozoic Zircons. The Fourth Lake Formation contains 49 grains (3.4% of the total grains ana- +5, which directly overlap the juvenile values of strata of the lyzed) with ages between 2802 Ma and 442 Ma and with cor- Fourth Lake Formation and igneous rocks of the Sicker arc. The similarity of U-Pb ages, ƐHf ðtÞ values, and previous responding ƐHf ðtÞ values ranging from +13 to -20 (Figure 6). These grains, plus the occurrence of some 360-300 Ma grains detrital zircon interpretations of the basal Comox Formation in the San Juan Islands suggests that Paleozoic zircons were with more evolved ƐHf ðtÞ signatures, suggest that sources for derived in large part from Paleozoic rocks of the Sicker arc strata of the Fourth Lake Formation also included pre-mid- and overlying Fourth Lake Formation. Paleozoic rocks or sediments with pre-mid-Paleozoic grains. The lack of evidence for such components of >380 Ma in 6.2. Early to Late Triassic Zircons. The Comox Formation samples from the Sicker Group and Buttle Lake Group [14] yields 12 single-grain ages ranging from 277 to 223 Ma raises the possibility that some grains in the Fourth Lake For- (Figure 7) with corresponding ƐHf ðtÞ values ranging from mation were sourced from rocks that are not currently part of +12 to +3 (Figure 8). The source of these grains 277- southern Wrangellia. 223 Ma is uncertain given the lull in magmatism in southern One potential source area for these older grains is the Wrangellia during this time period. The only magmatism Alexander terrane, given that previous workers have sug- recorded in southern Wrangellia during this age range is gested that the Wrangellia and Alexander terranes were near, the eruption of the Karmutsen flood basalts, which occurred or even built partially on, each other during late Paleozoic from 232 to 225 Ma [2], but zircon is rare in basaltic rocks. time [12, 13, 15], and that some portions of the Alexander terrane contain older and more evolved crustal components 6.3. Late Triassic to Middle Jurassic Zircons. The abundance [17, 18, 65]. Three distinct portions of the Alexander terrane of 217-160 Ma detrital zircons (Figure 7) indicates that strata have yielded U-Pb ages and Hf isotope data, including the of the Comox Formation may have been shed from latest Saint Elias Mountain region [17, 18], Southeast Alaska [63, Triassic-Jurassic igneous rocks of the West Coast Crystalline 66], and the Banks Island assemblage [65] (Figure 1). Of Complex, Island Intrusive Suite, and/or Bonanza Group, these regions, connections with the Banks Island assemblage which are widespread on Vancouver Island (Figure 1). These and the Saint Elias Mountain region are considered most rocks yield U-Pb ages of 202-165 Ma [33, 35, 51] and juvenile likely given the geologic, geochronologic, isotopic, and paleo- to intermediate isotopic signatures [36], which are quite magnetic evidence that Vancouver Island was within prox- similar to Hf isotope values from the Comox Formation imity of northern Wrangellia, thusly within proximity of (Figures 8 and 9). the Alexander terrane, prior to Early Cretaceous time [65– 67]. Although Precambian–early Paleozoic grains within 6.4. Late Jurassic to Cretaceous Zircons. Previous researchers Wrangellia cannot be traced to one individual entity of the have suggested that Late Jurassic and Cretaceous detrital zir- Alexander terrane due to the relatively low percentage of cons from the Comox Formation were shed from the Coast these grains in our samples, they do fit well within detrital Mountain Batholith [39, 40, 42–44, 58, 59]. More specifically, zircon populations observed in the Alexander terrane. Matthews et al. [42] suggested derivation from the central CMB because of consistencies in timing and volume of mag- 6. Provenance of the Comox Formation matic flux events in the central CMB (160-140 Ma, 120- 78 Ma, and 55-48 Ma; Gehrels et al. [68]) correlating to peak 6.1. Paleozoic Zircons. The Comox Formation yields Paleo- ages. In contrast, Huang [43] has suggested that Nanaimo zoic peak ages of 367, 356, and 341 Ma (Figure 7). With the detrital zircons were shed from the southern CMB given that inclusion of Comox samples from Matthews et al. [42], a the observed ages are more similar. cumulative Paleozoic peak age for the Comox Formation is A comparison of our age and Hf isotopic data with the 364 Ma with a range of ages encompassing rocks of the Sicker information available from the southern and central CMB and Buttle Lake Groups (Figure 7). This suggests that zircons is shown in Figure 9. We conclude that the age distributions

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 Lithosphere 15

from strata of the Comox Formation are not an excellent fit isotopes [29, 30], depletion of radiogenic Nd in felsic volcanic with either portion of the CMB, especially the lower and intrusive rocks of the Cowichan uplift, and the occur- proportion of 130-100 Ma ages and higher proportion of rence of peraluminous granitoids within the Saltspring Intru- ~170-160 Ma ages in Comox samples. ƐHf ðtÞ values of sive Suite [14]. Comox zircons are also significantly more juvenile than zir- There are two scenarios in which intermediate early arc cons from the central CMB (Figure 9). We accordingly agree magmas were generated in southern Wrangellia. (1) with previous workers that the Late Jurassic and Cretaceous Subduction-related magmatism was contaminated with sed- zircons in the Comox Formation were likely shed from the iment derived from the nearby Alexander terrane from the CMB, but conclude that the available data do not support upper crust and was consumed during subduction-related firm connections with either the southern or northern por- arc magmatism. (2) In the scenario we prefer, a transitional tions of the batholith. portion of the Alexander terrane crust underlays the nascent arc of southern Wrangellia and was since buried due to 7. Discussion extensive arc construction. This would provide the older, more evolved crust required to produce rocks with interme- 7.1. Tectonic Implications for Paleozoic Wrangellia. Previous diate to evolved geochemical signatures observed in the ear- workers have suggested that Wrangellia originated at the liest phase of the Sicker arc. The presence of Alexander-like margin of the Paleo-Pacific and Paleo-Arctic realms, devel- older crust would also provide a source of pre-Devonian zir- oping within close proximity to the Alexander terrane from cons deposited into the Fourth Lake Formation throughout Devonian through Permian time [13, 15–18, 69–77]. the Carboniferous. Previously reported connections between the Alexander and Wrangellia terranes [13–15, 74, 77, 78] combined with 7.1.2. Mississippian (354-323 Ma; Figure 10(b)). In southern our new data lead to the following integrated tectonic model Wrangellia, slab rollback rifted the preexisting Sicker arc for Paleozoic development of northern and southern Wran- and formed a back-arc spreading center in its place, focused gellia (Figures 10(a)–10(d)). in the Alberni area between the Buttle Lake and south Cow- ichan areas [14]. Within the new back-arc rift region, local 7.1.1. Late Devonian to Earliest Mississippian (370-355 Ma; bimodal magmatism accompanied the emplacement of Figure 10(a)). According to Colpron and Nelson [74] and VMS-type deposits [14]. The earliest deposits of the Buttle Nelson et al. [77], the Alexander terrane was extruded west- Lake Group, ribbon cherts associated with the Fourth Lake ward out of the Paleo-Arctic into the NE Paleo-Pacific realm. Formation, are deposited in the new back-arc basin and on An east-dipping subduction zone was established at this time portions of the rifted fragments of the Sicker arc [14]. A mod- outboard of western Laurentia, with the overriding plate ern analog to Paleozoic southern Wrangellia’s rifted arc-type bearing the Alexander terrane [13]. This subduction zone setting is the Taupo-Tonga-Kermadec arc system in the was proposed to extend into an intraoceanic setting beyond southwest Pacific Ocean [15, 79]. Maximum arc activity in the terrane boundary of Alexander, where it gave rise to the southern Wrangellia occurred in the Middle to Late Missis- Skolai arc of northern Wrangellia and Sicker arc of southern sippian, as shown by a cumulative detrital zircon age peak Wrangellia [13, 15]. of 337 Ma. Igneous rocks with these ages have not been pre- In northern Wrangellia, initiation of back-arc rifting is viously recognized in Wrangellia. On Vancouver Island, they inferred from the presence of coeval nonarc gabbros in are probably covered by extensive younger formations. The Wrangellia (Steele Creek complex) and in the Craig subter- southern Wrangellia arc was highly active throughout the rane of the Alexander terrane (Constantine complex and Carboniferous base on the abundance of detrital zircons related gabbro dikes) yielding ages of ca. 363 Ma [13]. These within this age range. Corresponding Hf isotope data are gabbros suggest a close connection between the Alexander highly juvenile compared to earlier Wrangellian sources. terrane and northern Wrangellia. The shift to exceptionally juvenile Hf values in the Sicker Although such gabbros are not recognized in southern arc may reflect the progressive rifting of the arc northward Wrangellia’s Sicker arc, Late Devonian nonarc basalt in the coupled with a process similarly invoked to explain juvenile lower Duck Lake Formation may similarly represent the εNd signatures in volcaniclastic rocks of the Klinkit Group onset of back-arc rifting during subduction-initiated arc (late Paleozoic) that were deposited on older pericratonic development of characteristic IAT, L-IAT, and minor E- strata of the Yukon-Tanana composite terrane in northern MORB-type rocks of the Sicker Group [15, 20–22]. The British Columbia and Yukon [80]. They proposed that the Sicker Group preserves the beginning of the southern Wran- voluminous asthenospheric melts that fed Klinkit volcanism gellia arc, ca. 370-365 Ma, on top of oceanic to transitional were rapidly and repeatedly emplaced along coated conduits, crust bordering the Alexander terrane. Relatively intermedi- insulated from sources of crustal contamination [80]. ate Hf isotope values from Late Devonian to Early Mississip- In contrast with the arc-rift to back-arc environment rep- pian detrital zircons from the Fourth Lake Formation suggest resented by Mississippian Sicker Group volcanic rocks, in that evolved older crust must have been present in order to northern Wrangellia, basalts and basaltic andesites of the Sta- influence the chemistry of the early stages of the Sicker arc. tion Creek Formation were more likely the products of arc- Further geochemical indicators for an evolved crustal com- axial magmatism. They range as old as ca. 352 Ma, as shown ponent in early arc construction include VMS deposits in by a U-Pb age determination on a rhyolite near the exposed the Myra Falls area yielding elevated levels of radiogenic Pb bottom of the section. At present, there is no other absolute

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 16 Lithosphere

Late Devonian - earliest Miss. Mississippian Late Miss. - Early Pennsylvanian Late Pennsylvanian - Early Permian (370-355 Ma) N (354-323 Ma) (322-315 Ma) (314-280 Ma)

nWR BGP nWR nWR CRP nWR LGP SC SGP sWR AX MC AX BVMS DGB AX sWR AX sWR AVMS AXa sWR sWR sWR sWR

= Nascent nWR arc = Old spreading center = Fourth Lake and Telwood Frm. = Transient AX crust = Fourth Lake Formation (a) (b) (c) (d)

Figure 10: Plate-tectonic reconstruction of Paleozoic Wrangellia during the (a) Late Devonian to earliest Mississippian (370-355 Ma), (b) Mississippian (354-323 Ma), (c) Early Pennsylvanian (322-315 Ma), and (d) Late Pennsylvanian to Early Permian (314-280 Ma). Abbreviations: AX: the Alexander terrane; AXa: Admiralty subterrane of Alexander terrane; nWR: northern Wrangellia; sWR: southern Wrangellia; AVMS: Alberni VMS deposits; BVMS: Bedingﬁeld uplift and Dragon property VMS deposits; BGP: Barnard Glacier porphyry; CRP: Centennial Ridge pluton; LGP: Logan Glacier porphyry; DGB: Donjek Glacier Batholith; SGP: Steel Glacier pluton.

age constraint on Station Creek volcanism. It may well have within both terranes and across their boundary [15]. The persisted through the Middle Mississippian interval docu- Early Permian Hasen Creek Formation was deposited as a mented in this study. clastic wedge atop the now-defunct Station Creek arc; Hasen Creek conglomerates were also deposited over the exhumed 7.1.3. Early Pennsylvanian (322-315 Ma; Figure 10(c)). Dur- gabbros [13]. ing this time interval, Sicker Group magmatic centers con- Arc-related magmatism continued in southern Wrangel- tinue to migrate towards the northwest (modern lia, with a northwesterly migration of locally bimodal mag- coordinates) [14]. Deposition of low-energy marine sedi- matism into the Bedingfield uplift and at the Dragon ments of the Fourth Lake and Thelwood Formation contin- property (northwestern Vancouver Island; Figure 1), where ued. Northward-younging detrital zircon peak ages in the 312-300 Ma volcanic rocks and associated VMS deposits are Fourth Lake Formation also reflect the direction of magmatic recognized [14]. The volcanic succession is overlain by and depositional migration. ~290 Ma tuffs, the final volcanic deposits recorded in south- Volcanism of the Skolai arc continued through this time ern Wrangellia [14]. Deposition of the Fourth Lake and Thel- interval, preserved in arc-type rocks of the upper Station wood formations continued in a back-arc setting. The lack of Creek Formation [15]. During this time, the arc system in deformation and postcollisional plutons suggest that south- northern Wrangellia is thought to have undergone a subduc- ern Wrangellia was largely unaffected by the collisional tion reversal, with evidence in changing chemistry from the events of northern Wrangellia and the Alexander terrane. lower to upper Station Creek igneous rocks [13, 15]. A rever- After collision, the subduction zone shifted behind the sal in subduction is also inferred in southern Wrangellia, Alexander terrane, as observed in emplacement of the Don- marked by the migration of volcanism from the Alberni area jek Glacier Batholith (286-284 Ma), the Steele Glacier pluton in the Mississippian to northwestern Vancouver Island by ca. (291 Ma), and other granitoids with ages from 290 to 280 Ma 312 Ma [14]. This change in subduction broke along the in the Alexander terrane [15]. In Early Permian time extinct Late Devonian spreading center and initiated the (~280 Ma), a fragment from the Alexander terrane, the encroachment of the Alexander terrane towards northern Admiralty subterrane, clogged the subduction zone, which Wrangellia [13, 15]. finally shuts off the arc magmatism within the Alexander terrane [15, 80]. Limestones of the Mt. Mark Formation on 7.1.4. Late Pennsylvanian to Early Permian (314-280 Ma; Vancouver Island record the cessation of arc-related volca- Figure 10(d)). Northern Wrangellia and Alexander collided, nism of southern Wrangellia. Fossiliferous Lower Permian as the Alexander block entered the Wrangellia subduction limestones of the Pybus Formation overlap the suture zone [15]. The collision led to the exhumation of basement between the Craig and Admiralty subterranes in southeast- gabbros and stitched the northern portion of Wrangellia ern Alaska [81]. These coeval limestone formations in south- and Alexander together, indicated by emplacement of the ern Wrangellia and Alexander terrane indicate that they Barnard Glacier porphyry (ca. 307 Ma), Centennial Ridge shared a quiescent postcollisional tectonic regime and depo- pluton (304 Ma), and the Logan Glacier porphyry (307 Ma) sitional environment.

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 Lithosphere 17

8. Conclusions Acknowledgments

This study presents the first robust detrital zircon analysis of Field work for this project was supported by NSF EAR Paleozoic strata from southern Wrangellia. Detrital zircons 1347375. Laboratory analyses were supported by NSF EAR in sedimentary rocks of the Fourth Lake Formation primar- 1649254 to the Arizona LaserChron Center. Thanks are due ily record penecontemporaneous magmatic activity of the to Graham Nixon for presubmission comments and edits. Sicker arc. Cumulative detrital zircon U-Pb ages peak at Lithosphere reviewers Grant Lowey and Sarah Roeske made 337 Ma. This broad Middle Mississippian peak lies between many helpful comments that improved the manuscript. We fi U-Pb igneous age clusters of 352-335 Ma and 312-300 Ma thank Rebecca Alberts for assistance in the eld and Mark identified by Ruks [14]. The combined data represents a Pecha, Nicky Giesler, Chelsi White, Kojo Plange, Clay Kelty, series of closely spaced magmatic episodes in the Sicker arc Gayland Simpson, Heather Alvarez, and Ken Kanipe for system between Late Devonian and Early Permian time. assistance with the laboratory analyses and logistics. Each episode, illustrated in our updated tectonic model, reflects a significant change in the arc and back-arc configu- ration in southern Wrangellia, which in turn were caused by Supplementary Materials changes in subduction geometry and eventually subduction Supplementary Materials include two data tables. Data table polarity. Progressively younger peak ages are observed from 1 (DR.1) contains U-Pb zircon data from this study and Ruks sample locations from south to north, suggesting volcanic [14]. Data table 2 (DR. 2) contains Lu-Hf isotopic data from centers in southern Wrangellia migrated northward during this study and Sm-Nd data from Ruks [14]. (Supplementary back-arc basin development. Materials) Mississippian and younger detrital zircons yield predominately juvenile Hf isotopic signatures, whereas interme- References diate Hf signatures are observed in Late Devonian zircons derived from the nascent Sicker arc. We propose that Late [1] D. L. Jones, N. J. Silberling, and J. Hillhouse, “Wrangellia – a Ɛ ðtÞ Devonian zircons with intermediate Hf values formed displaced terrane in northwestern North America,” Canadian in Sicker arc-related magmas that assimilated minor Journal of Earth Sciences, vol. 14, no. 11, pp. 2565–2577, 1977. amounts of crust from an Alexander terrane-like source. [2] A. R. Greene, J. S. Scoates, D. Weis, E. C. Katvala, S. Israel, and Incorporation of an older, more evolved component from G. T. Nixon, “The architecture of oceanic plateaus revealed by the Alexander terrane is further supported by the presence the volcanic stratigraphy of the accreted Wrangellia oceanic of detrital zircon grains older than 380 Ma with highly vari- plateau,” Geosphere, vol. 6, no. 1, pp. 47–73, 2010. able Hf signatures. This suggests the presence of a Precam- [3] A. Sutherland Brown, C. J. Yorath, R. G. Anderson, and brian to lower Paleozoic crustal component, similar in K. Dom, “Geological Maps of Southern Vancouver Island, profile to detrital zircon populations in the Alexander ter- LITHOPROBE 1,” Geological Survey of Canada, Open File rane. A key result of this study is that although southern 1272, 10 Sheets, 1986. Wrangellia is not connected to the Alexander terrane by [4] N. W. D. Massey, “Geology and mineral resources of the common geological elements and stitching plutons as is Alberni-Nanaimo Lakes sheet, Vancouver Island 92F/1W, ” northern Wrangellia [13, 15], the detrital zircon ages and 92F/2E, and part of 92F/7E, B.C. Ministry of Energy, Mines Hf isotopic signatures presented here provide a strong case and Petroleum Resources. Paper 1992-2, 1995. “ for early Alexander terrane linkages. Thus, all of Wrangellia [5] N. W. D. Massey, Geology and mineral resources of the Cow- ” and the Alexander terrane can be considered as parts of a ichan Lake sheet, Vancouver Island 92C/16, B.C. Ministry of coevolving tectonic system, beginning with the birth of Energy, Mines and Petroleum Resources. Paper 1992-3, 1995. [6] N. W. D. Massey, “Geology and mineral resources of the Dun- Wrangellia marked by the inception of the Sicker and Skolai ” arcs. We recommend a study similar to this one in northern can sheet, Vancouver Island 92B/13, B.C. Ministry of Energy, Mines and Petroleum Resources. Paper 1992-4, 1995. Wrangellia which is aimed at U-Pb age distributions and Hf isotopic data, which could be compared to results from [7] C. J. Yorath, A. Sutherland Brown, and N. W. D. Massey, “LITHOPROBE, southern Vancouver Island, British Colum- southern Wrangellia to provide an integrated analysis of ” – bia: geological survey of Canada, Bulletin, vol. 498, p. 145, the combined Sicker-Skolai arc back-arc system. 1999. [8] J. G. Smith and E. M. McKevett, “The Skolai Group in the McCarthy B-4, C4 and C-5 quadrangles, Wrangell Mountains, Data Availability Alaska. U.S,” Geological Survey Bulletin, vol. 1274-Q, pp. 1–26, 1970. Data tables of data generated from this study are [9] G. C. Bond, “A Late Paleozoic volcanic arc in the eastern included in supplementary files. Ruks [14] data is also Alaska Range, Alaska,” Journal of Geology, vol. 81, no. 5, included. pp. 557–575, 1973. [10] D. H. Richter and D. L. Jones, “Structure and stratigraphy of the eastern Alaska Range,” AAPG Memoir, Arctic Geology, Conflicts of Interest vol. 19, pp. 408–420, 1973. [11] P. B. Read and J. W. H. Monger, “Pre-Cenozoic volcanic The authors declare that they have no conflicts of interest. assemblages of the Kluane and Alsek ranges, southwestern

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 18 Lithosphere

Yukon Territory, Geological Survey of Canada,” Open File, the Paleozoic Sicker Group and associated intrusive rocks, vol. 381, p. 96, 1976. Vancouver Island, British Columbia,” Current Research, Geo- [12] M. C. Gardner, S. C. Bergman, G. W. Cushing et al., “Pennsyl- logical Survey of Canada, Paper 86-1A, 1986. vanian pluton stitching of Wrangellia and the Alexander Ter- [26] C. L. Sluggett, “Uranium-lead age and geochemical constraints rane, Wrangell Mountains, Alaska,” Geology, vol. 16, no. 11, on Paleozoic and Early Mesozoic magmatism in Wrangellia pp. 967–971, 1988. Terrane, Saltspring Island, British Columbia,” B.Sc. thesis, [13] S. Israel, L. Beranek, R. M. Friedman, and J. L. Crowley, “New University of British Columbia, Vancouver, B.C., 2003. ties between the Alexander terrane and Wrangellia and impli- [27] M. D. Hannington, C. E. J. de Ronde, and S. Peterson, “Sea- cations for North America Cordilleran evolution,” Litho- floor tectonics and submarine hydrothermal systems,” in Eco- sphere, vol. 6, no. 4, pp. 270–276, 2014. nomic Geology One Hundredth Anniversary Volume,J.W. [14] T. W. Ruks, “Stratigraphic and Paleotectonic Studies of Paleo- Hedenquist, J. F. H. Thompson, R. J. Goldfarb, and J. P. – zoic Wrangellia and it’s Contained Volcanogenic Massive Sul- Richards, Eds., pp. 111 141, Society of Economic Geologists, fide (VMS) Occurrences, Vancouver Island, British Columbia, 2005. Canada,” Dissertation, University of British Columbia, 2015. [28] A. Andrew, R. L. Armstrong, and D. Runkle, “Neodymium– – [15] L. Beranek, C. Van Staal, W. McClelland, N. Joyce, and strontium lead isotope study of Vancouver Island igneous “ rocks,” Canadian Journal of Earth Sciences, vol. 28, no. 11, S. Israel, Late Paleozoic assembly of the Alexander- – Wrangellia-Peninsular composite terrane, Canadian and Alas- pp. 1744 1752, 1991. “ kan Cordillera,” Geological Society of America Bulletin, [29] A. Andrew and C. I. Godwin, Lead- and strontium-isotope vol. 126, no. 11-12, pp. 1531–1550, 2014. geochemistry of Paleozoic Sicker Group and Jurassic Bonanza Group volcanic rocks and Island Intrusions, Vancouver Island, [16] L. P. Beranek, C. R. van Staal, S. M. Gordee, W. C. British Columbia,” Canadian Journal of Earth Sciences, vol. 26, McClelland, S. Israel, and M. G. Mihalynuk, “Tectonic sig- pp. 894–907, 1989. nificance of Upper Cambrian–Middle Ordovician maficvolcanic “ rocks on the Alexander terrane, Saint Elias Mountains, north- [30] C. I. Godwin and M. S. Robinson, Galena lead isotopes, Buttle western Canada,” Journal of Geology,vol.120,pp.293–314, 2013. Lake mining camp, Vancouver Island, British Columbia, Can- ada, Economic,” Geology, vol. 91, pp. 549–562, 1996. [17] L. P. Beranek, C. R. van Staal, W. C. McClelland, S. Israel, and “ [31] S. D. Samson, P. J. Patchett, G. E. Gehrels, and R. G. Anderson, M. G. Mihalynuk, Baltican crustal provenance for Cambrian- “ Ordovician sandstones of the Alexander terrane, North Amer- Nd and Sr isotopic characterization of the Wrangellia terrane and implications for crustal growth of the Canadian Cordil- ican Cordillera: evidence from detrital zircon U-Pb geochro- ” – nology and Hf isotope geochemistry,” Journal of the lera, The Journal of Geology, vol. 98, no. 5, pp. 749 762, 1990. Geological Society of London, vol. 170, no. 1, pp. 7–18, 2013. [32] G. D. Webster, J. W. Haggart, C. Saxifrage, B. Saxifrage, “ fi [18] L. P. Beranek, C. R. van Staal, W. C. McClelland, S. Israel, and C. Gronau, and A. Douglas, Globally signi cant Early Perm- M. G. Mihalynuk, “Detrital zircon Hf isotopic compositions ian crinoids from the Mount Mark Formation in Strathcona Provincial Park, Vancouver Island, British Columbia — pre- indicate a northern Caledonian connection for the Alexander ” terrane,” Lithosphere, vol. 5, no. 2, pp. 163–168, 2013. liminary analysis of a disappearing fauna, Canadian Journal of Earth Sciences, vol. 46, no. 9, pp. 663–674, 2009. [19] A. R. Greene, J. S. Scoates, D. Weis, and S. Israel, “Geochemis- [33] C. E. Isachsen, “Geology, geochemistry, and cooling history of try of Triassic flood basalts from the Yukon (Canada) segment the West Coast Crystalline Complex and related rocks, Meares of the accreted Wrangellia oceanic plateau,” Lithos, vol. 110, Island and vicinity, Vancouver Island, British Columbia,” no. 1-4, pp. 1–19, 2009. vol. 24, Tech. Rep. 10, Canadian Journal of Earth Sciences, “ [20] N. W. D. Massey and S. J. Friday, Geology of the Cowichan 1987. Lake Area, Vancouver Island (92C/16),” in Geological Field- [34] G. T. Nixon, J. L. Hammack, V. M. Koyanagi et al., “Prelimi- work 1986, pp. 223–229, B.C. Ministry of Energy, Mines and nary Geology of the Quatsino - PortMcNeill Map Areas, Petroleum Resources Paper 1987-1, 1987. Northern Vancouver Island (92L/12, 11),” in Geological Field- “ [21] N. W. D. Massey and S. J. Friday, Geology of the Chemainus work 1993, B. Grant and J. M. Newell, Eds., pp. 63–85, B.C. ” River-Duncan Area, Vancouver Island (92C/16; 92B/13), in Ministry of Energy, Mines and Petroleum Resources, Paper – Geological Fieldwork 1987, pp. 81 91, B.C. Ministry of Energy, 1994-1, 1994. Mines and Petroleum Resources Paper 1988-1, 1988. [35] S. M. DeBari, R. G. Anderson, and J. K. Mortensen, “Correla- [22] N. W. D. Massey and S. J. Friday, “Geology of the Alberni- tion among lower to upper crustal components in an island Nanaimo Lakes Area, Vancouver Island (91F/1W, 92F/2E arc: the Jurassic Bonanza arc, Vancouver Island, Canada,” and part of 92F/7),” in Geological Fieldwork 1987, pp. 61–74, Canadian Journal of Earth Sciences, vol. 36, pp. 1371–1413, B.C. Ministry of Energy, Mines and Petroleum Resources 1999. Paper 1989-1, 1989. [36] B. D. Paulson, “Magmatic processes in the Jurassic Bonanza [23] S. Juras, “Geology of the polymetallic volcanogenic Buttle Lake arc: insights from the Alberni region of Vancouver Island, camp, with emphasis on the Price-Hillside, central Vancouver Canada,” unpublished Ph.D. thesis, Western Washington Uni- Island, British Columbia,” Ph.D. thesis, University of British versity Graduate School Collection, 2010. Columbia, 1987. [37] J. E. Muller and J. A. Jeletzky, “Geology of the Upper Creta- [24] R. R. Parrish and V. J. McNicoll, “U-Pb age determinations ceous Nanaimo Group, Vancouver Island and Gulf Islands, from the southern Vancouver Island, area, British Columbia, British Columbia: Geological Survey of Canada, Paper 69- Radiogenic Age and Isotope Studies: Report 5,” Geological 25,” p. 77, 1970. Survey of Canada, Paper 91-2, 1992. [38] T. D. J. England, “Late Cretaceous to Paleogene structural evo- [25] M. T. Brandon, M. J. Orchard, R. R. Parrish, A. Sutherland lution of the Georgia basin, southwestern British Columbia,” Brown, and C. J. Yorath, “Fossil ages and isotopic dates from Ph.D. thesis, Memorial University, Newfoundland, 1990.

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 Lithosphere 19

[39] P. S. Mustard, “The Upper Cretaceous Nanaimo Group, Geor- [53] E. H. Brown, “Obducted nappe sequence in the San Juan gia basin,” in Geology and Geological Hazards of the Vancouver Islands–northwest Cascades thrust system, Washington and Region, Southwestern British Columbia, J. W. H. Monger, Ed., British Columbia,” Canadian Journal of Earth Sciences, vol. 481, pp. 27–95, Geological Survey of Canada Bulletin, vol. 49, no. 7, pp. 796–817, 2012. 1994. [54] J. W. Haggart, “A synthesis of Cretaceous stratigraphy, Queen [40] D. C. Katnick and P. S. Mustard, “Geology of Denman and Charlotte Islands, British Columbia,” in Evolution and Hydro- Hornby islands, British Columbia: implications for Nanaimo carbon Potential of the Queen Charlotte Basin, British Colum- basin evolution and formal definition of the Geoffrey and bia, G. J. Woodsworth, Ed., pp. 253–277, Geological Survey of Spray formations, Upper Cretaceous Nanaimo Group,” Cana- Canada Paper 90-10, 1991. dian Journal of Earth Sciences, vol. 40, no. 3, pp. 375–393, [55] J. W. Haggart and E. S. Carter, “Biogeography of latest Jurassic 2003. and Cretaceous mollusc and radiolarian faunas of the Insular [41] H. Bain and S. Hubbard, “Stratigraphic evolution of a long- Belt, British Columbia, suggest minimal northward displace- lived submarine channel system in the late Cretaceous ment,” Geological Society of America Abstracts with Programs, Nanaimo Group, British Columbia, Canada,” Journal of Sedi- vol. 26, no. 7, article A148, 1994. mentary Geology, vol. 337, pp. 113–132, 2016. [56] J. W. Haggart, P. D. Ward, and W. Orr, “Turonian (Upper [42] W. A. Matthews, B. Guest, D. Coutts, H. Bain, and S. Hubbard, Cretaceous) lithostratigraphy and biochronology, southern “Detrital zircons from the Nanaimo basin, Vancouver Island, Gulf Islands, British Columbia, and northern San Juan Islands, British Columbia: an independent test of Late Cretaceous to Washington State,” Canadian Journal of Earth Sciences, Cenozoic northward translation,” Tectonics, vol. 36, no. 5, vol. 42, no. 11, pp. 2001–2020, 2005. pp. 854–876, 2017. [57] P. D. Johnstone, P. S. Mustard, and J. A. MacEachern, “The basal [43] C. Huang, “Refining the chronostratigraphy of the Lower unconformity of the Nanaimo Group, southwestern British Nanaimo Group, Vancouver Island, Canada, using detrital zir- Columbia: a Late Cretaceous storm-swept rocky shoreline,” con geochronologyB.Sc. Thesis, Simon Fraser University,” Canadian Journal of Earth Sciences,vol.43,pp.1165–1181, 2006. 2018. [58] P. S. Mustard, R. R. Parrish, and V. McNicoll, “Provenance of [44] C. Huang, S. E. Dashtgard, B. A. P. Kent, H. D. Gibson, and the Upper Cretaceous Nanaimo Group, British Columbia: evi- W. A. Matthews, “Resolving the architecture and early evolu- dence from U-Pb analyses of detrital zircons,” vol. 52, Strati- tion of a forearc basin (Georgia Basin, Canada) using detrital graphic Development in Foreland Basins, SEPM Special zircon,” Science Report, vol. 9, no. 1, p. 15360, 2019. Publication,, 1995. [45] J. C. Lassiter, “Geochemical investigations of plume related [59] J. B. Mahoney, P. S. Mustard, J. W. Haggart, R. M. Friedman, lavas: constraints on the structure of mantle plumes and the C. M. Fanning, and V. J. McNicoll, “Archean zircons in Creta- nature of plume/lithosphere interactions [Ph.D. thesis],” Uni- ceous strata of the western Canadian Cordillera: The “Baja versity of California, Berkeley, 1995. B.C.” hypothesis fails a “crucial test”,” Geology, vol. 27, no. 3, [46] D. J. T. Carson, “Petrography, chemistry, age and emplace- pp. 195–198, 1999. ment of plutonic rocks of Vancouver Island: Geological Survey [60] G. E. Gehrels, V. A. Valencia, and J. Ruiz, “Enhanced preci- of Canada Paper 72-44,” 1973. sion, accuracy, efficiency, and spatial resolution of U-Pb ages [47] J. E. Muller, “Geology of Vancouver Island: Geological Survey by laser ablation-multicollector-inductively coupled plasma of Canada, Open File 463; 1:250 000 Scale,” 1977. mass spectrometry,” Geochemistry Geophysics Geosystems, [48] J. E. Muller, R. K. Wanless, and W. D. Loveridge, “A Paleozoic vol. 9, no. Q03017, p. 13, 2008. zircon age of the Westcoast Crystalline Complex of Vancouver [61] G. Gehrels and M. Pecha, “Detrital zircon U-Pb geochronol- Island, British Columbia,” Canadian Journal of Earth Sciences, ogy and Hf isotope geochemistry of Paleozoic and Triassic vol. 11, pp. 1717–1722, 1974. passive margin strata of western North America,” Geosphere, [49] G. T. Nixon and A. J. Orr, “Recent revisions to the Early Meso- vol. 10, no. 1, pp. 49–65, 2014. zoic stratigraphy of northern Vancouver Island (NTS 102I; [62] A. Pullen, M. Ibanez-Mejia, G. Gehrels, D. Giesler, and 092L) and metallogenic implications, British Columbia,” in M. Pecha, “Optimization of a laser ablation-single collector- Geologic Fieldwork 2006, pp. 163–178, British Columbia Min- inductively coupled plasma-mass spectrometer (Thermo Ele- istry of Energy, Mines and Petroleum Resources, Paper 2007- ment 2) for accurate, precise, and efficient zircon U-Th-Pb 1, 2006. geochronology,” Geochemistry, Geophysics, Geosystems, [50] G. T. Nixon, J. L. Hammack, V. M. Koyanagi et al., Geology, vol. 19, no. 10, pp. 3689–3705, 2018. Geochronology, Lithogeochemistry and Metamorphism of the [63] M. R. Cecil, G. E. Gehrels, M. N. Ducea, and P. J. Patchett, “U- Mahatta Creek Areas, Northern Vancouver Island (NTS Pb-Hf characterization of the central Coast Mountains batho- 92L/11, and Parts of 92L/05, 12, and 13), British Columbia lith: implications for petrogenesis and crustal architecture,” Geomap 2011-02, 1:50,000 Scale, 2011. Lithosphere, vol. 3, no. 4, pp. 247–260, 2011. [51] R. M. Friedman and G. T. Nixon, “U–Pb zircon dating of [64] H. Bahlburg, J. D. Vervoort, S. A. DuFrane, V. Carlotto, Jurassic porphyry Cu(–Au) and associated acid sulphate sys- C. Reimann, and J. Cardenas, “The U–Pb and Hf isotope evi- tems, northern Vancouver Island, British Columbia,” Geolog- dence of detrital zircons of the Ordovician Ollantaytambo For- ical Association of Canada – Mineralogical Association of mation, southern Peru, and the Ordovician provenance and Canada Annual Meeting, Victoria, B.C., 1995. paleogeography of southern Peru and northern Bolivia,” Jour- [52] P. D. Clift, T. Pavlis, S. M. DeBari, A. E. Draut, M. Rioux, and nal of South American Earth Sciences, vol. 32, no. 3, pp. 196– P. B. Kelemen, “Subduction erosion of the Jurassic Talkeetna- 209, 2011. Bonanza arc and the Mesozoic accretionary tectonics of west- [65] C. J. Tochilin, G. E. Gehrels, J. L. Nelson, and J. B. Mahoney, ern North America,” Geology, vol. 33, no. 11, pp. 881–884, “U-Pb and Hf isotope analysis of detrital zircons from the 2005. Banks Island assemblage (coastal British Columbia) and

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021 20 Lithosphere

southern Alexander terrane (southeast Alaska),” Lithosphere, [79] N. Mortimer, P. B. Gans, J. M. Palin, S. Meffre, R. H. Herzer, vol. 6, no. 3, pp. 200–215, 2014. and D. N. B. Skinner, “Location and migration of Miocene- [66] C. White, G. E. Gehrels, M. Pecha et al., “U-Pb and Hf isotope Quaternary volcanic arcs in the SW Pacific region,” Journal analysis of detrital zircons from Paleozoic strata of the south- of Volcanology and Geothermal Research, vol. 190, pp. 1–10, ern Alexander terrane (southeast Alaska),” Lithosphere, 2010. vol. 8, no. 1, pp. 83–96, 2016. [80] R.-L. Simard, J. Dostal, and C. F. Roots, “Development of late [67] G. Plafker, W. J. Nokleberg, and J. S. Lull, “Bedrock geology Paleozoic volcanic arcs in the Canadian Cordillera: an example and tectonic evolution of the Wrangellia, Peninsular and Chu- from the Klinkit Group, northern British Columbia and south- gach terranes along the Trans-Alaska crustal transect in the ern Yukon,” Canadian Journal of Earth Sciences, vol. 40, no. 7, Chugach Mountains and southern Copper River Basin, pp. 907–924, 2003. Alaska,” Journal of Geophysical Research, vol. 94, no. B4, [81] S. M. Karl, P. W. Layer, A. G. Harris, P. J. Haeussler, and B. L. pp. 4255–4295, 1989. Murchey, “The Cannery Formation - Devonian to Early Perm- [68] G. Gehrels, M. Rusmore, G. Woodsworth et al., “U-Th-Pb geo- ian arc-marginal deposits within the Alexander terrane, south- chronology of the Coast Mountains batholith in north-coastal eastern Alaska,” Studies by the U.S. Geological Survey in British Columbia: constraints on age and tectonic evolution,” Alaska, 2008–2009: U.S. Geological Survey Professional Paper Geological Society of America Bulletin, vol. 121, pp. 1341– 1776-B, 2010. 1361, 2009. [82] N. W. D. Massey, D. G. MacIntyre, P. J. Desjardins, and R. T. [69] C. M. Soja, “Significance of Silurian stromatolite sphinctozoan Cooney, “Digital geology map of British Columbia: Tile NM9 reefs,” Geology, vol. 22, pp. 355–358, 1994. Mid Coast, B.C,” B.C. Ministry of Energy and Mines Geofile [70] D. R. Bazard, R. F. Butler, G. E. Gehrels, and C. M. Soja, 2005-2, 2005. “Paleomagnetism of the Early Devonian Karheen Formation, [83] N. W. D. Massey, P. J. Desjardins, and R. T. Cooney, “Digital southeast Alaska: implications for Alexander terrane paleoge- map of British Columbia: Tile NM10 Southwest B.C,” B.C. ography,” Geology, vol. 23, pp. 707–710, 1995. Ministry of Energy, Mines and Petroleum Resources, GeoFile [71] G. E. Gehrels, R. F. Butler, and D. R. Bazard, “Detrital zircon 2005-3, 2005. ” geochronology of the Alexander terrane, southeastern Alaska, [84] K. M. Cohen, S. C. Finney, P. L. Gibbard, and J.-X. Fan, Geological Society of America Bulletin, vol. 108, no. 6, Updated 2020, vol. 36, The ICS international chrono- – pp. 0722 0734, 1996. stratigraphic chart, Episodes, 2013. “ [72] R. F. Butler, G. E. Gehrels, and D. R. Bazard, Paleomagnetism [85] J. D. Vervoort, P. J. Patchett, J. Blichert-Toft, and F. Albarede, of Paleozoic strata of the Alexander terrane, southeastern “Relationships between Lu-Hf and Sm-Nd isotopic systems in Alaska,” Geological Society of America Bulletin, vol. 109, ” – the global sedimentary system, Earth and Planetary Science no. 10, pp. 1372 1388, 1997. Letters, vol. 168, pp. 79–99, 1999. [73] C. M. Soja and L. Krutikov, “Provenance, depositional setting, [86] J. D. Vervoort and J. Blichert-Toft, “Evolution of the depleted and tectonic implications of Silurian polymictic conglomerate mantle: Hf isotope evidence from juvenile rocks through in Alaska’s Alexander terrane,” in The Terrane Puzzle: New time,” Geochimica et Cosmochimica Acta, vol. 63, pp. 533– Perspectives on Paleontology and Stratigraphy from the North 556, 1999. American Cordillera, R. B. Blodgett and G. D. Stanley Jr., “ Eds., pp. 63–75, Geological Society of America Special Paper [87] A. Bouvier, J. D. Vervoort, and J. D. Patchett, The Lu-Hf and 442, 2008. Sm-Nd isotopic composition of CHUR: constraints from unequilibrated chondrites and implications for the bulk com- [74] M. Colpron and J. L. Nelson, A Palaeozoic Northwest Passage: position of terrestrial planets,” Earth and Planetary Science incursion of Caledonian, Baltican and Siberian terranes into Letters, vol. 273, no. 1-2, pp. 48–57, 2008. eastern Panthalassa, and the early evolution of the North Amer- “ ican Cordillera, vol. 318, no. 1, 2009Geological Society of Lon- [88] J. D. Vervoort and P. J. Patchett, Behavior of hafnium and don, Special Publications, 2009. neodymium isotopes in the crust: constraints from crustally “ derived granites,” Geochimica et Cosmochimica Acta, vol. 60, [75] M. Colpron and J. L. Nelson, A Palaeozoic Northwest Passage – and the Timanian, Caledonian and Uralian connections of no. 19, pp. 3717 3733, 1996. some exotic terranes in the North American Cordillera, Chap- [89] M. R. Cecil, M. E. Rusmore, G. E. Gehrels et al., “Along-strike ter 31,” Arctic Petroleum Geology: Geological Society of London, variation in the magmatic tempo of the Coast Mountains Memoirs, A. M. Spencer, A. Embry, D. Gautier, A. Stoupakova, batholith, British Columbia, and implications for processes and K. Sorensen, Eds., vol. 35, pp. 463–484, 2011. controlling episodicity in arcs,” Geochemistry, Geophysics, [76] E. L. Miller, N. Kuznetsov, A. Soboleva, O. Udoratina, M. J. Geosystems, vol. 19, p. 16, 2018. Grove, and G. E. Gehrels, “Baltica in the Cordillera?,” Geology, vol. 39, no. 8, pp. 791–794, 2011. [77] J. L. Nelson, M. Colpron, and S. Israel, “The Cordillera of British Columbia, Yukon, and Alaska: Tectonics and Metallo- geny,” in Tectonics, Metallogeny, and Discovery, M. Colpron, T. Bissig, B. G. Rusk, and J. F. H. Thompson, Eds., vol. 17, no. - ch.3pp. 53–109, The North American Cordillera and Similar Accretionary Settings: Society of Economic Geologists Special Publication, 2013. [78] M. Colpron, J. L. Nelson, and D. C. Murphy, “Northern Cor- dilleran terranes and their interactions through time,” GSA Today, vol. 17, no. 4, pp. 4–10, 2007.

Downloaded from http://pubs.geoscienceworld.org/gsa/lithosphere/article-pdf/doi/10.2113/2021/7866944/5292498/7866944.pdf by guest on 23 September 2021