Canadian Journal of Earth Sciences

Cambrian detrital zircon signatures of the northern Canadian Cordilleran passive margin, Liard area, Canada: evidence of sediment recycling, non-Laurentian ultimate sources and basement denudation.

Journal: Canadian Journal of Earth Sciences

Manuscript ID cjes-2016-0127.R1

Manuscript Type: Article

Date Submitted by the Author: 21-Feb-2017 Complete List of Authors: McMechan,Draft Margaret; Geological Survey of Canada, Currie, Lisel; Geological Survey of Canada – Calgary Ferri, Filippo; BC Ministry of Natural Gas Development, Tenure and Geoscience Branch Matthews, William; University of Calgary, Department of Geosciences O'Sullivan, Paul; Apatite To Zircon, Inc.; Geosep Services

detrital zircon, , sediment recycling, non-Laurentian sources, Keyword: Cordilleran passive margin

Note: The following files were submitted by the author for peer review, but cannot be converted to PDF. You must view these files (e.g. movies) online. cjes-2016-0127.R1-suplemental-2.xlsm

https://mc06.manuscriptcentral.com/cjes-pubs Page 1 of 57 Canadian Journal of Earth Sciences

Cambrian detrital zircon signatures of the northern Cordilleran passive margin, Liard area, Canada: evidence of sediment recycling, non-Laurentian ultimate sources and basement denudation

Margot McMechan 1, Lisel Currie 2, Filippo Ferri 3, William Matthews 4 and Paul O’Sullivan 5

1 Geological Survey of Canada – Calgary, 3303 33 rd St. NW, Calgary, AB T2L 2A7;

[email protected]

2 Geological Survey of Canada – Calgary, 3303 33 rd St. NW, Calgary, AB T2L 2A7;

[email protected]

3 Tenure and Geoscience Branch, British Columbia Ministry of Natural Gas Development

Victoria, BC; [email protected] Draft

5 Department of Geoscience, University of Calgary, 2500 University Drive NW, Calgary, AB

T2N 1N4; [email protected]

5 GeoSep Services, 1521 Pine Cone Rd, Moscow, Idaho 83843, USA;

[email protected]

Author responsible for correspondence:

Margot McMechan Geological Survey of Canada – Calgary, 3303 33 rd St. NW, Calgary, AB

T2L 2A7; [email protected] ; phone: 403-292-7154, fax 403-292-5377

1

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 2 of 57

Cambrian detrital zircon signatures of the northern Cordilleran passive margin, Liard area,

Canada: evidence of sediment recycling, non-Laurentian ultimate sources and basement denudation.

Margot McMechan, Lisel Currie, Filippo Ferri, William Matthews and Paul O’Sullivan

ABSTRACT Detrital zircon U-Pb age probability distributionsDraft for the Cambrian Vizer formation (informal) and Mount Roosevelt Formation (middle member) of the northern Canadian Cordilleran passive margin indicate extensive recycling from ~1.7 to 1.6 Ga Paleoproterozoic sandstones and

Proterozoic and Lower Cambrian strata, respectively. The units have minor or no first cycle input from Laurentian basement. The lower part of the Vizer formation contains North American

Magmatic Gap (1610 to1490 Ma) detrital zircons and lacks ultimate Grenvillian sourced grains, indicating that the grains were likely sourced from a nearby Mesoproterozoic basin and have an ultimate non-Laurentian source. Detrital zircon U-Pb ages of 670 to 640 Ma from the middle member of the Mount Roosevelt Formation indicate associated volcanic clasts were locally sourced, and are not of syn-sedimentary Middle Cambrian age. Provenance of these units was indirectly impacted by the Liard Line basement feature.

Detrital zircon U-Pb age probability distributions from the northern Canadian Cordilleran passive margin indicate sediments were sourced from the east in the Early Cambrian

(Terreneuvian; Vizer formation and correlatives) and the northeast during Early Cambrian

2

https://mc06.manuscriptcentral.com/cjes-pubs Page 3 of 57 Canadian Journal of Earth Sciences

(Series 2) deposition of Sekwi Formation and correlative strata. In the early Middle Cambrian

the middle member of the Mount Roosevelt Formation was primarily locally sourced, whereas

the upper member was derived from Laurentian basement to the east and southeast. The change

from reworked Paleoproterozoic cover in the Terrenuvian to primary basement sources in the

Middle Cambrian suggests significant denudation of the basement occurred southeast of the

Liard Line.

Keywords: detrital zircon, Cambrian, sediment recycling, non-Laurentian sources,

Cordilleran passive margin

Draft

3

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 4 of 57

INTRODUCTION

The area now forming the western margin of Laurentia has had a protracted history of extension, sediment accumulation and shortening. Rifting events in the middle and late Neoproterozoic

(Cryogenian , Ediacaran) established the margin (e.g. Stewart 1972; Arnott and Hein 1986; Ross et al. 1995) and resulted in the deposition of two prominent sedimentary successions: the

Neoproterozoic Windermere Assemblage and the Ediacaran to Middle Cordilleran passive margin. Several northeast-trending, basement-controlled fault zones affected sedimentation within the Laurentian continental margin succession of the Canadian Cordillera.

The most prominent occurs near 60 o N Draftat a feature termed the Liard Line, which has been interpreted as a Neoproterozoic and lower Paleozoic transfer fault zone separating a wide lower plate margin on the north from a narrow upper plate margin on the south (Fig. 1; Cecile et al.

1997). This paper describes the detrital zircon U-Pb age signatures of the pre- quartz arenite of theVizer formation (informal) and Middle Cambrian volcanic clast-bearing conglomerates of the Mount Roosevelt formation exposed near the Liard Line and uses this information together with published detrital zircon age information from the various units described in the geological framework below to constrain their provenance and outline the major changes in detrital zircon source areas during the Early and Middle Cambrian. Our results indicate: 1) extensive recycling from the ~ 1.70 to 1.6 Ga Paleoproterozoic sedimentary blanket deposited over northern Laurentia and a complete absence of ultimate Grenville sourced grains in the Vizer formation and middle member of the Mount Roosevelt formation; 2) the component of North American magmatic gap (1610 to 1490 Ma) aged zircons found in the lower part of the

4

https://mc06.manuscriptcentral.com/cjes-pubs Page 5 of 57 Canadian Journal of Earth Sciences

Vizer Formation were ultimately derived from non-Laurentian sources and most likely recycled

from a nearby Mesoproterozoic basin; 3) the Vizer formation is most likely Lower Cambrian

Terreneuvian in age; 4) the major change to first-cycle basement sources in the upper member of

the Mount Roosevelt Formation likely reflects denudation of the Laurentian basement; and 6)

stratigraphic changes associated with the Liard Line had an impact on the provenance of the

Vizer and Mount Roosevelt formations.

GEOLOGICAL FRAMEWORK

The ancestral North American continent Laurentia, comprising Archean to Early Proterozoic

microcontinents separated by collisional and magmatic orogenic belts, formed between 2.0 and

1.8 Ga (Hoffman 1989; Whitmeyer andDraft Karlstom 2007) as part of a larger continental

amalgamation termed Columbia (Rogers and Santosh 2002). Between 1.7-1.68 Ga and 1.65-1.60

Ma juvenile orogens formed along the southeastern side of Laurentia (Whitmeyer and Karlstrom

2007) and thick skinned deformation of the Forward Orogeny occurred inboard of the northwest

Laurentia margin prior to 1.663 Ga (Cook and MacLean 1995). Rifting along the northwestern

margin of Laurentia around 1.65 Ga resulted in rapid deposition of the 13 km thick Wernecke

Supergroup (Furlanetto et al. 2016) and most likely the 6 km thick Muskwa assemblage (Fig. 2).

Detrital zircon studies record a predominant Laurentian provenance for both Wernecke

(Furlanetto et al. 2009, 2016) and Muskwa (Ross et al. 2001) strata. Paleocurrent indicators and

age constraints suggest they are the distal equivalents of a sandstone package derived from

westward flowing river systems that blanketed at least the northern part of Laurentia in the late

Paleoproterozoic. Eastern remnants of this package are now preserved in the Athabasca, Thelon,

and Hornby Bay basins (Fig. 1; Young 1978; Ross et al. 2001; Rainbird et al. 2005). Similar

5

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 6 of 57

detrital zircon age distributions from the Athabasca and Muskwa basins support their correlation, whereas, detrital zircon age distributions from the Dubawnt Supergroup of the Thelon Basin are dissimilar and indicate a different provenance (Rainbird et al. 2005; Rainbird and Davies 2007).

Local intense deformation of the northwestern edge of Laurentia during the Racklan Orogeny occurred after deposition of the Wernecke Supergroup and before intrusion of the Wernecke breccia ca. 1595 Ma (Thorkelson et al. 2005; Furlanetto et al. 2013). By 1.5 Ga north Australia, south Australia and eastern Antartica were firmly attached to the western margin of Laurentia

(Karlestrom et al. 2001; Thorkelson et al. 2005; Ross and Villeneuve 2003; Stewart et al. 2010;

Furlanetto et al. 2013; Medig et al. 2014).

Extension along the western margin of DraftLaurentia between about 1490 to 1400 Ma did not produce full continental separation. It is recorded in western Canada by the PR1 (basal Fifteen

Mile Group, Medig et al. 2014) and Belt-Purcell basins (Fig. 1; Ross and Villeneuve 2003).

Sometime after 1380 Ma, renewed subsidence of the northwestern margin of Laurentia initiated deposition of the Pinguicula Group (Medig et al. 2010) and the roughly correlative Dismal Lakes

Group (Cook and MacLean, 1995). Extension along the northern margin of Laurentia resulted in intrusion of the Mackenzie dyke swarm (1267 ± 2 Ma, LeCheminant and Heaman 1989) and the extrusion of the Coppermine River basalts above the Dismal Lake Group. After a period of planation, predominantly platformal strata of the Mackenzie Mountains Supergroup and Shaler

Group (Fig. 1) were deposited within an intracratonic basin developed along the northwest edge of Laurentia between about 1000 Ma and 775 to 720 Ma (Heaman et al. 1992; Rainbird et al.

1996; Milton 2015). Thick sandstones in the lower part of the succession with 1600 to 1000 Ma detrital zircons derived predominantly from the Grenville orogenic belt of eastern Laurentia were deposited by west-northwest flowing rivers (Rainbird et al. 1997, 2012).

6

https://mc06.manuscriptcentral.com/cjes-pubs Page 7 of 57 Canadian Journal of Earth Sciences

The local extrusion of mantle sourced volcanic rocks at 775 Ma (Little Dal Basalts) in the

Mackenzie Mountains correlates with the Gunbarrel magmatic event and signals the initiation of

continental extension along this part of the western edge of Laurentia that would ultimately lead

to the deposition of the Windermere Supergroup (Milton, 2015). In northwest Laurentia, the first

narrow rift basin formed in the Mackenzie Mountains and was filled by the Coates Lake Group

whose uppermost strata have yielded a 732.2 ± 3.9 Ma Re-Os age (Rooney et al. 2014). The

onset of continental rifting along the length of the western margin of Laurentia to form either the

proto-Pacific Ocean (e.g. Stewart 1972; Arnott and Hein 1986; Ross et al. 1989; Ross et al.

1995) or a large intercontinental rift basin (Colpron et al. 2002) resulted in deposition of the

Windermere Supergroup. Along the segment north of 57 o N sedimentation began sometime after

728 +8/-7 Ma (Evenchick et al. 1984) andDraft before 717.43 ± 0.14 Ma (Macdonald et al. 2010). In

this same area local rift-related volcanism is recorded by the Gataga volcanics (689.1 ± 4.6 Ma;

Ferri et al. 1999).

In the late Ediacaran a major thermal and rifting event initiated deposition of the Ediacaran to

Middle Devonian Cordilleran passive margin succession (e.g. Bond and Kominz 1984), from

which our samples were obtained. In most areas along the eastern edge of this passive margin an

Ediacaran to Lower Cambrian sandstone package forms the basal unit (Gog Group and

correlatives). This facies appears to be offset 250 km to the northeast across the Liard Line

(Aitken 1993). Changes in older strata suggest this feature may have been established during

Cryogenian (Aitken 1993) or Paleoproterozoic(?) rifting events and reactivated in the Ediacaran

(Aitken, 1993). Local depositional thickening and volcanism within the Canadian Cordilleran

passive margin indicate widespread episodic extension beginning in the Middle Cambrian

(Goodfellow et al. 1995). In the Liard area this resulted in the development of the northwest-

7

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 8 of 57

trending Roosevelt graben system above the western Liard Line and the local deposition of coarse clastic strata of the Mount Roosevelt Formation (Taylor and Stott 1973; Fritz 1979, 1991;

Post and Long 2008).

SAMPLES ANALYSED

Samples were collected from two Lower Paleozoic horizons along the western part of the Liard

Line: the Vizer and Mount Roosevelt formations (Fig. 1).

Vizer formation (informal)

A quartz arenite unit exposed in the Cariboo Range immediately south of the B.C. – Yukon border (Locality 3b in Fig. 1, 2) has beenDraft informally called the Vizer formation by Benjamin et al. (2011), who studied the unit in detail as a potential source for fine-grained, pure (>97%) quartz sands to be used by the petroleum industry in the fracking process. The Vizer formation is comprised mainly of pale grey to white, fine- to medium-grained, thin- to thick- bedded, sub- angular to subrounded quartz arenite (Benjamin et al. 2011). It is unconformably overlain by

Silurian carbonate and its base is not exposed because it lies in the immediate hangingwall of the

Larsen Fault (McMechan et al. 2012). One of our samples is from the lowermost exposures and the other comes from the middle part. Taylor and Stott (1980) suggested a Lower and (?) Middle

Cambrian age and correlated these strata mainly with strata now called the Mount Roosevelt

Formation, whereas McMechan et al. (2012) suggested an Ediacaran to Lower Cambrian age based on lithologic similarity with the undated quartz arenite unit underlying the Mount

Roosevelt Formation in the southwest Liard area, as shown in Figure 2. Recently, Fallas et al.

(2014) suggested these strata should be assigned to the uppermost Cambrian to Crow

Formation. The Crow Formation comprises variably subarkosic sandstone, conglomerate,

8

https://mc06.manuscriptcentral.com/cjes-pubs Page 9 of 57 Canadian Journal of Earth Sciences

maroon siltstone, minor carbonate and local volcanics in exposures immediately northwest of the

Cariboo Range (Pigage 2009; Fig. 1, localities 2, 3a). From a sedimentology and compositional

perspective it seems unlikely that the nearly pure quartz arenite succession exposed in the

Cariboo Range is part of the variably subfeldspathic package comprising the Crow Formation,

particularly since crossbeds in the Crow Formation indicate a west to southwest transport

direction (Pigage, 2009). Detrital zircon U-Pb age results for 6 samples of the Crow Formation

reported by Pigage (2009) are compared with detrital zircon signatures from the Mount

Roosevelt and Vizer formations below.

Mount Roosevelt Formation

The lower Middle Cambrian Mount RooseveltDraft Formation (Fritz 1979; 1991) represents the progradation of an alluvial fan delta into a marine environment (Post and Long 2008). Three

members comprise the up to 1500 m thick formation. The lower member is characterized by

ooid-bearing siltstone and sandstone, interbedded with carbonate and hematitic conglomerate.

The thick middle member comprises polymict, cobbly pebble conglomerate and the upper

member is comprised of karstified dolostone, calcareous cemented conglomerate and sandstone

and limestone. Deposition of the Mount Roosevelt Formation occurred in fault-controlled basins

with the thickness and distribution of the formation controlled by the location of active faults

(e.g. Fritz 1991; Ferri et al. 1999). The thickest development of the Mount Roosevelt Formation

occurs in the hanging wall of the Forcier fault immediately northwest of Muncho Lake (Post and

Long, 2008). Our sample (locality 3B in Fig. 1, 2) comes from a volcanic clast-bearing polymict

conglomerate in the middle member of the formation.

9

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 10 of 57

Detrital zircons from a second sample of the Mount Roosevelt Formation were analysed by

Gehrels and Ross (1998) and reanalysed by Gehrels and Pecha (2014), who included this sample in the Atan Group following the now discontinued terminology used by Taylor and Stott (1973).

This sample is from a coarse-grained to pebbly, carbonate-cemented quartz arenite horizon in a clastic unit that is overlain by shale and limestone (Gehrels and Ross 1998). It comes from the upper member of the Mount Roosevelt Formation exposed 14 km southeast of the middle member sample location. In the Liard area, the Middle Cambrian Mount Roosevelt Formation unconformably overlies Paleoproterozoic strata of the Muskwa assemblage in the east and up to

300 m of fine- to coarse-grained quartz arenite that likely correlates with the Vizer formation of the Cariboo Range in the west (Fig. 2, see above).

METHODS Draft

Detrital zircons from three new samples, two from the Vizer formation quartz arenite and one from conglomerate of the Middle Cambrian Mount Roosevelt Formation, were processed and analysed by AtoZ Incorporated, Bosie, Idaho using the LA- ICP-MS at the Geoanalytical

Laboratory, Washington State University, Pullman, Washington. Full details of the analytical approach, grain descriptions and the data can be found in the Supplementary Information.

For grains that yield ages older than ca. 1000 Ma, the 207 Pb/ 206 Pb age is used, whereas for younger grains the 206 Pb/ 238 U age is used. Analyses were filtered following the criteria of

Gehrels and Pecha (2014). Analyses with >10% uncertainty (1 σ) in 206 Pb/ 238 U or the 206 Pb/ 207 Pb age are not included in our analysis. Concordance is based on 206 Pb/ 238 U age/ 206 Pb/ 207 Pb age, with 100% considered concordant. Analyses with >20% discordance (<80% concordance) or

>5% reverse discordance (>105% concordance) are not included.

10

https://mc06.manuscriptcentral.com/cjes-pubs Page 11 of 57 Canadian Journal of Earth Sciences

The resulting detrital zircon U-Pb age distributions for our samples and the Mount Roosevelt

Formation sample analysed by Gehrels and Pecha (2014) are compared with: 1) the U-Pb zircon

age distributions for other late Ediacaran and Cambrian aged strata; and 2) potential intermediate

sources (i.e. Proterozoic sedimentary strata through which grains may have been recycled) using

direct comparison of normalized probability plots (calculated using the Excel™ plug-in Isoplot

v. 4.15; Ludwig 1998, Ludwig 2012, Berkeley Geochronolgy Center, 2016) and

multidimentional scaling (MDS). MDS is a method of objectively interpreting large datasets.

Using the IsoplotR program of Vermeesch (2013; 2017) a ‘map’ of data points is produced on

which ‘dissimilar’ or unrelated samples plot far apart and ‘similar’ samples, such as samples

with similar provenance, or a sample and its ultimate or intermediate provenance, cluster closely

together. The ‘map’ is based on the D valuesDraft (the maximum vertical separations between two

cumulative density functions) that are produced as part of the Kolmogorov-Smirnov test

(Vermeesch 2013). This study uses data from numerous sources including a few with small

sample sizes (n <50). Recent work by Pullen et al. (2014) suggests that the number of concordant

grains analysed for each sample needs to approach 300 for more meaningful statistical analysis.

Analyses from these other publications were also filtered using the criteria of Gehrels and Pecha

(2014).

DETRITAL ZIRCON GEOCHRONOLOGY

Results

Detrital zircons from two samples of the Vizer formation were analysed. Sample FF10-1 from

the lowest exposures of the unit contains a dominant 1900 to 1500 Ma age fraction with a lesser

2900 to 2500 Ma population. A distinct “shoulder” between 1700 and 1500 Ma is formed in part

11

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 12 of 57

by 5 concordant grains with ages between 1605 and 1538 Ma that fall into a hiatus in Laurentian magmatism between 1610 and 1490 Ma termed the North American magmatic gap (NAMG;

Ross and Villenueve 2003; Fig. 3). Sample 11EPF55A from the middle part of the formation has a preponderance of 2000 to 1600 Ma grains with less common 3000 to 2200 Ma grains. Both samples have a significant peak between 1900 and 1750 Ma. Neither contains zircons younger than 1500 Ma and thus the zircon ages provide no useful information on the depositional age of the Vizer formation. The samples are similar and have therefore been grouped together on the normalized probability and MDS plots (Fig. 4, 5).

One sample from the middle member of the Middle Cambrian Mount Roos evelt Formation (11EPF67A) yielded dominantly 2100 toDraft 1600 Ma and 2800 to 2300 Ma grains, with two concordant 700 to 600 Ma grains and one reliable ca. 3200 Ma grain. The probability plot has a prominent peak between 1900 and 1700 Ma and a lesser peak between 2650 and 2500 Ma (Fig.

3). Zircons older than 2200 Ma make up a greater proportion of the zircons in this sample than in the Vizer formation quartz arenites. There were no Early or Middle Cambrian (541 to 497 Ma) grains.

DISCUSSION

Comparison with Cordilleran Passive Margin Cambrian to Ordovician Units

Detrital zircon results from our samples and the sample from the upper member of the Mount

Roosevelt Formation (Gehrels and Pecha, 2014) were compared with those of other Late

Ediacaran to Cambro-Ordovician units in the Canadian Cordilleran passive margin to determine which units could have had similar detrital zircon provenances.

12

https://mc06.manuscriptcentral.com/cjes-pubs Page 13 of 57 Canadian Journal of Earth Sciences

Direct comparison of normalized probability plots (Fig. 4) for these units suggest the Vizer

formation has a different detrital zircon provenance from 1) the Franklin Mountain and Sekwi

formations and the Sekwi correlative unit because of the abundance of 1490 to 1000 Ma grains

and relative lack of 2000 to 1750 Ma grains in those samples; 2) the upper member of the Mount

Roosevelt Formation sample because of the abundance of 2500 to 2300 Ma grains and absence

of a strong 1880 to 1750 Ma peak in that sample; 3) the Crow Formation Northwest Liard

locality samples because of the presence of Ordovician aged zircon and the absence of 1750 to

1500 Ma grains in those samples; 4) the Yanks Peak Formation because of the lack of 1700 to

1600 Ma grains and presence of 1400 to 1000 Ma grains in that sample; and likely 5) the Crow

Formation Coal River locality samples because of the subdued 2000 to 1750 Ma peak and

component of 1400 to 600 Ma grains (8%)Draft in those samples.

The MDS plot for comparison of detrital zircon signatures of Ediacaran to Cambro-Ordovician

units along the Canadian Cordilleran passive margin (Fig. 5) suggests the detrital zircon

signature of the Vizer formation is most similar to those of the McNaughton, Hamill, Vampire-

Narchilla units and the middle member of the Mount Roosevelt Formation. The similar detrital

zircon profiles of the Vizer formation with Lower Cambrian (Terrenevian) aged samples from

the McNaughton and Hamill formations suggests that the Vizer is most likely part of the

regionally extensive transgressive sandstone package at the base of the Cordilleran passive

margin succession. In the Mackenzie Mountains stratigraphic relationships (Aitken 1993;

MacNaughton et al. 2008) show that the Vampire interval is the basinal equivalent of nearshore

sandstones of the Backbone Ranges Formation. Comparable detrital zircon profiles of the

Vampire-Narchilla and Vizer formations in the Liard area suggests that a similar facies

relationship likely occurs between the Lower Cambrian (Terrenuevian; Pigage et al. 2015)

13

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 14 of 57

Vampire-Narchilla formation and the nearshore sandstones of the Vizer formation. This would also indicate a Terrenuevian Early Cambrian age for the quartz arenites of the Vizer formation.

Direct comparison of normalized probability plots of Ediacaran to Cambro-Ordovician units along the Canadian Cordilleran passive margin (Fig. 4) suggest that our sample from the middle member of the Mount Roosevelt Formation 1) is most similar in provenance to the Vizer formation samples because of the width of the 2000 to 1600 Ma peak and subdued presence of

2500 to 2000 Ma grains, 2) has a different provenance from the Franklin Mountain and Sekwi formations and Sekwi correlative unit because of the abundance of 1490 to 1000 Ma ultimate

Grenville sourced grains in those samples and 3) has a different provenance from the upper member of the Mount Roosevelt FormationDraft sample because of the relative abundance of 2400 to 2300 Ma aged grains and lack of a 1900 to 1700 Ma peak in that sample. The MDS plot (Fig. 5) suggests the detrital zircon signature of the middle member of the Mount Roosevelt Formation sample is most similar to the McNaughton, Hamill and Vampire-Narchilla units. Both direct comparison of normalized probablility plots (Fig. 4) and the MDS plot (Fig. 5) suggest the detrital zircon signature of the upper member of the Mount Roosevelt Formation differs from all the other Ediacaran to Cambro-Ordovician units sampled to date in the Canadian Cordilleran passive margin succession.

Comparison with Potential Intermediate Detrital Source Units

Due to its mechanical durability, zircon is commonly recycled from one sedimentary unit to another without impacting the reliability of the U-Pb age determinations (e.g. Dickenson et al.,

2009). The significance of sedimentary recycling for Neoproterozoic and Cambrian strata in northwestern Laurentia has been documented by Hadlari et al. (2012) and Lane and Gehrels

14

https://mc06.manuscriptcentral.com/cjes-pubs Page 15 of 57 Canadian Journal of Earth Sciences

(2014) who concluded most Neoproterozoic and Cambrian units in the Mackenzie – Wernecke

Mountains and adjacent Richardson Trough were recycled from intermediate (Proterozoic

sedimentary) sources. Almost all the Proterozoic successions mentioned previously in the

Geological Framework were exposed during the Cambrian (Fig. 1), making them potential

intermediate sources of detrital zircons found in the Vizer and Mount Roosevelt formations.

Direct comparison of normalized probability plots for the Vizer and Mount Roosevelt formation

with those from potential Proterozoic intermediate sources (Fig. 6) rules out the Neoproterozoic

Shaler, Mackenzie Mountains and Windermere supergroup units exposed in the Wernecke and

Mackenzie mountains and northwest Liard area (Fig. 1, 2, localities 1, 3b) as intermediate source

units for both the Mount Roosevelt and Vizer formations because they contain abundant 1490 to

1000 Ma (Grenville source) aged grainsDraft that the Vizer and Mount Roosevelt formations lack.

Unit Ps is not considered an important intermediate source for the Vizer formation because of the

abundant 2.0 to 1.9 Ga aged detrital zircons in that unit. Paleoproterozoic sandstones of the

Athabasca, Muskwa and Wernecke successions and Neoproterozoic deep water turbidites of the

Horsethief Creek Group and Ediacaran Yusezyu Formation are potential sources for the Vizer

formation. However, none of the potential intermediate source units contain the robust

contribution 2.5 to 2.2 Ga contribution found in the upper Mount Roosevelt Formation (Fig. 6).

Comparison with published probability density plots for samples from the Mesoproterozoic

Fifteen Mile and Pinguicula groups in the Hart Creek Inlier (Medig et al. 2012, 2014) suggest

that these would be unlikely intermediate sources for the Vizer and Mount Roosevelt formations

because they contain abundant detrital zircons with ages less than or equal to 1500 Ma (Fig. 7).

Whereas, the basal unit of the Pinguicula Group in the Wernecke inlier (Fig. 1, location 1a) has

very few <1500 Ma detrital zircon grains and could be a potential intermediate source (Fig. 7 A-

15

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 16 of 57

C). Two of the three samples from this unit contained 1610 to 1490 Ma aged NAMG detrital zircons (Fig. 7; Medig et al. 2012).

The MDS plot for comparison of detrital zircon signatures of the Vizer and Mount Roosevelt formations with potential Proterozoic intermediate sources (Fig. 8) suggests that 1) the detrital zircon signature of the Vizer formation is most similar to that of the Ediacaran Yusezyu

Formation; 2) the detrital zircon signature of the upper member of the Mount Roosevelt

Formation differs significantly from all the potential intermediate sources; and 3) the detrital zircon signature of the middle member of the Mount Roosevelt Formation is similar to those of the Paleoproterozoic Muskwa succession, Neoproterozoic Horsethief Creek Group and the Ediacaran Yusezyu Formation. Data forDraft the Pinguicula Group is not yet available for MDS comparison.

The close similarity of detrital zircon signature for the Vizer formation to those of Proterozoic intermediate sources suggests there is little if any first cycle Laurentian basement sourced material in this unit. The abundance of 2500 to 2300 Ma zircons in the upper member of the

Mount Roosevelt Formation (25%, Fig. 4) as compared to the Muskwa (5%) implies a direct cratonic source (Gehrels and Ross 1998). Magmatism between 2500 and 2300 Ma has a limited distribution in western Canada (Fig. 1, Canadian Geochronology Knowledgebase 2013). The most likely ultimate source for the 2500 to 2390 Ma grains are basement rocks in the western

Trans-Hudson Orogen. Other potential source areas include the La Ronge arc, Taltson magmatic belt, Rae and Hearne provinces and eastern Trans-Hudson Orogen (Fig. 1). Abundant 2380 to

2300 Ma grains in the upper Mount Roosevelt Formation could have come from the Buffalo

Head, Thorsby or Wabamun terranes (Gehrels and Ross 1998) or rocks exposed near the western margin of the Rae Province, north of the Athabasca basin (Fig. 1). Basement rocks in the eastern

16

https://mc06.manuscriptcentral.com/cjes-pubs Page 17 of 57 Canadian Journal of Earth Sciences

Trans-Hudson Orogen are considered a less likely source for 2500 to 2300 Ma grains because the

Trans-Hudson Orogen would provide a significant component of 1900 to 1800 Ma detrital

zircons, which are not observed in the upper Mount Roosevelt Formation sample. Recycled

zircons with 2500 to 2300 Ma ages occur in the northern Hottah Terrane (Davis et al. 2015). This

is considered to be an unlikely source since Cambrian sandstone samples unconformably

overlying the northern Hottah terrain at Great Bear Lake contain only a small component (~10%)

of these grains and a large component of 2000 to 1750 Ma grains (Hadlari, 2012). The 2050 to

1800 Ma grains in the upper Mount Roosevelt sample match well with plutons in several nearby

terranes to the east. The 2900 to 2500 Ma grains could have come from the Rae, Hearn or Slave

provinces (Fig. 1). The middle member of the Mount Roosevelt Formation contains a noticeable

component of 2500 to 2400 Ma aged detritalDraft zircon grains that were ultimately derived from the

Laurentian basement areas mentioned above. They may indicate that a small component of first

cycle Laurentian basement material made its way into the Roosevelt graben or a lack of

sufficient sampling of zircon from feldspathic strata at the top of the underlying Muskwa

Assemblage (Tuchodi Formation, n = 19; Ross et al. 2001). Two concordant first-cycle

approximately 670 and 640 Ma detrital zircon grains occur in the middle Mount Roosevelt

Formation sample. Magmatism of this age is known from only two localities in western Canada

(Canadian Geochronology Knowledgebase, 2013; Yukon Geochronology Database, 2015). The

Gataga volcanics (689.1 ± 4.6 Ma; Ferri et al. 1999), which outcrop 75 km southwest of the

middle Mount Roosevelt Formation sample location and the Pool Creek Syenite (Fig. 2; 650 to

640 Ma, Pigage and Mortensen, 2004), located 150 km north of the middle Mount Roosevelt

sample locality. The older grain (669 ± 14 Ma) is intermediate in age while the younger (640 ±

10 Ma) detrital zircon grain is the same age as the Pool Creek Syenite. The erosion of volcanic

17

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 18 of 57

and/or plutonic rocks associated with the Gataga Volcanics and the Pool Creek Syenite intrusion form the likely sources. Occasional younger Neoproterozoic grains (690 to 640 Ma) are also found in samples from the lower Cambrian Vampire-Narchilla and Cambro-Ordovician Crow formations at the Coal River locality (Pigage 2009; Pigage et al. 2015). There was no detrital zircon record of Middle Cambrian volcanism in the middle member of the Mount Roosevelt

Formation sample.

Potential Sources 1610 to 1490 Ma “North American Magmatic Gap” Grains, lower Vizer formation

The detrital zircon age profile presented here for the lower Vizer formation is are the first reported occurrence of NAMG detrital Draftzircons (5% of accepted analyses) from the Cambrian Laurentian margin succession not associated with Grenville aged grains (Fig. 3). This absence of

Grenville aged grains places serious doubt on the possibility that the NAMG zircons were derived from eastern Laurentia, despite the occurance of 1520 and 1460 Ma magmatism in the

Pinware terrane (Fig. 1 inset map) of the easternmost Grenville orogen (Gower and Krogh 2002;

Heaman et al. 2004). The allochtonous Pinware terrane was accreted to the Archean and

Paleoproterozoic core of Laurentia during the (1250 to 980 Ma) Grenville Orogeny.

Consequently, any Pinwarian aged detrital zircon grains would be associated with younger (1350 to 1000 Ma) Grenville aged detrital zircon grains, as is found in the Shaler Group (Rayner and

Rainbird 2013). The complete absence of Grenville aged detrital zircon grains combined with the older age range of the Vizer detrital zircon grains (1608 to 1532 Ma) as compared to the age of

Pinwarian magmatism (1520 to 1460 Ma) rule out Laurentia as a potential ultimate source.

18

https://mc06.manuscriptcentral.com/cjes-pubs Page 19 of 57 Canadian Journal of Earth Sciences

NAMG aged detrital zircons derived from non-Laurentian sources have been found from

northcentral Yukon to southern Arizona in Mesoproterozoic basins formed along the western

margin of Laurentia (Medig et al. 2014 and references therein). The ultimate source of the 1610

to 1490 Ma detrital zircon grains is thought to be the Mawson and North Australia cratons

(Stewart et al., 2010; Medig et al., 2014). Direct recycling from the three exposed Canadian

Laurentian margin Mesoproterozoic successions into the Cambrian Vizer Formation is unlikely.

First, the PR1 basin (Fig. 1; basal Fifteen Mile Group) located 900 km northwest of the Liard

area is characterized by a unimodal 1500 Ma detrital zircon population (Fig. 7D; Medig et al.

2014). These strata are directly overlain by younger Mesoproterozoic sediments of the Fifteen

Mile Group and were not available to provide material directly to Cambrian strata (Medig et al.

2014). Second, the basal unit of the PinguiculaDraft Group (unit A) in the Wernecke Mountains 700

km northwest of the Liard area contain a small component of NAMG grains not associated with

younger grains (Fig. 7A-C; Medig et al. 2012). However, the Pinguicula Group in these areas is

unconformably overlain by the Neoproterozoic Mackenzie Mountains Supergroup (Medig et al.

2012) and was not available to provide material directly to Cambrian strata. Finally, the Belt-

Purcell basin, located 1200 km south of the Liard area (Fig. 1), is unconformably overlain by

Early Cambrian strata in northeastern exposures and potentially a source. The rarity of NAMG

grains in Ediacaran and Early Cambrian strata near the Belt-Purcell basin and to the north of this

basin (Gehrels and Pecha 2014; Matthews 2016) suggests the NAMG grains in the Vizer

formation were likely recycled from a more proximal source.

The absence of an exposed and therefore geologically viable source for the recycled NAMG

detrital zircon grains in the Vizer formation indicates that these grains most likely came from a

unit that has either been completely removed by lower Paleozoic erosion, is preserved only in the

19

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 20 of 57

subsurface or has not yet been sampled. Unit Ps in the northwest Liard area provides evidence for such a unit. Pigage (2009) suggested Proterozoic strata of unit Ps (Fig. 2) most likely correlate with the Mesoproterozoic Pinguicula Group on the basis of a detrital zircon signature that 1) lacks 1490 to 1000 Ma Grenville-sourced detrital zircons characteristic of the

Neoproterozoic Mackenzie Mountains and Shaler supergroups and eastern Windermere

Supergroup successions and 2) contains a strong 2000 to 1900 Ma peak that is absent in the

Paleoproteroic Wernecke, Athabasca and Muskwa successions but present in some Pinguicula

Group samples from the Wernecke Mountains (Fig. 6, 7). Although the exposed portion of unit

Ps could not be a source of NAMG grains in the lower part of the Cambrian Vizer formation because it is overlain by Neoproterozoic strata, its presence does indicate that a Mesoproterozoic succession younger than the Muskwa assemblageDraft (Paleoproterozoic) and older than the

Mackenzie Mountain Supergroup (Neoproterozoic –Tonian) occurs in the subsurface of the northwest Liard area. By analogy with Mesoproterozoic successions exposed along the western edge of Laurentia, this succession probably contained a component of NAMG aged detrital zircons.

During deposition of the lower part of the Vizer formation the postulated Mesoproterozoic succession would have been progressively truncated eastward and available to provide detritus to the basal Cambrian sandstone as it was progressively overlapped. Erosion and overlap of the postulated Mesoproterozoic succession would also have happened in the southern part of the

Liard area as Cambrian or Silurian strata directly overlie the Muskwa assemblage. Progressive removal of a Mesoproterozoic succession by sub-Cambrian erosion in the eastern and southern part of the Liard area forms the most geologically feasible source for the NAMG detrital zircon grains in the lower Vizer formation.

20

https://mc06.manuscriptcentral.com/cjes-pubs Page 21 of 57 Canadian Journal of Earth Sciences

Paleogeographic implications

New detrital zircon results from the Vizer and Mount Roosevelt formations combined with

previously published results show important changes in detrital zircon provenance occurred in

the northern Canadian Cordilleran passive margin succession during the Early and Middle

Cambrian. These changes are illustrated in Figure 9 using hypothetical drainage pathways for

four time intervals corresponding to deposition of the Vampire-Narchilla and Vizer formations

(Early Cambrian Terreneuvian series; Figure 9A); Sekwi Formation and correlatives (Early

Cambrian, Series 2; Figure 9B); middle member of the Mount Roosevelt Formation (early

Middle Cambrian; Figure 9C); and upper member of the Mount Roosevelt Formation (early

Middle Cambrian; Figure 9D) in the Coal River and Liard areas. Draft In the Canadian Cordilleran passive margin succession the Early Cambrian Terreneuvian and

Series 2 strata represent a major transgression that reached its maximum before the end of the

Early Cambrian. This was followed by a major regression represented by a disconformity at the

base of the Middle Cambrian and then a major Middle Cambrian transgression that progressively

onlapped the Laurentian craton (Aitken 1989; Gabrielse 1991). South of 55 o N most of the

western Canadian Laurentian basement east of the passive margin was directly overlapped by

Middle Cambrian strata, whereas north of the Liard Line the Cambrian overlapped large regions

of Proterozoic sedimentary cover east of the passive margin (Fig. 1).

Detrital zircon age distributions for the basinal Vampire-Narchilla and the nearshore Vizer

formations (Fig. 4 -6; Terreneuvian) indicate extensive recycling from Paleoproterozoic

Athabasca-Muskwa sandstones exposed to the east and southeast, no recycling from the

Neoproterozoic Mackenzie Mountain/Shaler and Windermere successions exposed to the

northeast and no significant first cycle basement material. The Liard Line most likely impacted

21

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 22 of 57

the detrital zircon provenance of the Vizer formation indirectly by controlling the southward extent of the Mackenzie Mountains Supergroup (Aitken 1993). The minor component of NAMG aged detrital zircon in the lower part of the Vizer formation was most likely recycled from

Mesoproterozoic strata exposed during the Terreneuvian in the eastern and southern Liard area

(discussed above). Further north in the Mackenzie Mountains, detrital zircon age distributions for the lower and upper Backbone Ranges Formation (Fig. 4) indicate derivation mainly from cratonic basement sources to the east, including the Hottah terrane, Great Bear magmatic zone,

Coronation Margin and the Slave craton (Fig. 9A; Leslie 2009; Lane and Gehrels 2014). The detrital zircon provenance of all three formations resulted from a westward paleoflow. The significant difference in detrital zircon signatures between the Vampire-Narchilla-Vizer and the

Backbone Ranges successions most likelyDraft reflects the greater amount of basement exposure east of the Backbone Ranges localities in the Terreneuvian.

The detrital zircon age distribution for the Lower Cambrian Series 2 Sekwi Formation in the

Mackenzie Mountains indicates extensive recycling of ultimate Grenville orogeny sourced grains. These most likely came from the Mackenzie Mountain Supergroup exposed to their northeast in the Mackenzie Arch (Leslie 2009; Lane and Gehrels 2014). Northeastern source areas are also indicated by the detrital zircon age distribution for the Sekwi correlative unit in the

Coal River area (Fig. 4, 9B), which has a broad spectrum of subdued peaks extending from 2800 to 1200 Ma that suggests derivation predominantly from the Mackenzie Mountain/Shaler supergroup (1500 to 1200 Ma grains) and cratonic sources including the Slave craton, Hottah terrane and Great Bear Magmatic zone (Pigage et al. 2015). Over most of the Liard area, Early

Cambrian Series 2 time is represented by a disconformity and no detrital zircon data is available.

22

https://mc06.manuscriptcentral.com/cjes-pubs Page 23 of 57 Canadian Journal of Earth Sciences

In the Liard area, Middle Cambrian strata are restricted to the syn-sedimentary Mount Roosevelt

graben. No detrital zircon data is available for Middle Cambrian deposits elsewhere along the

northern Canadian Cordilleran passive margin. The detrital zircon age distribution for the

conglomeratic middle member of the Mount Roosevelt Formation indicates a local provenance

consistent with the alluvial fan deposition model proposed by Post and Long (2008) with a minor

contribution from 2500 to 2400 Ma basement sources. Most detrital zircon grains were recycled

from the underlying Muskwa assemblage either directly or indirectly through Lower Cambrian

sandstones, and their provenance was indirectly influenced by the Proterozoic stratigraphic

changes associated with the Liard Line. The concordant 670 and 640 Ma detrital zircon grains

strongly suggest that volcanic clasts in the conglomerate were also derived locally as the only

magmatic rocks of this age near the CanadianDraft Cordilleran passive margin are in the Liard area

(discussed above; Fig. 9C).

In contrast, the detrital zircon age distribution for the sandstone-dominated upper member of the

Mount Roosevelt Formation suggests derivation mainly from more distant first cycle cratonic

sources to the east or southeast (Fig. 9D). The upper member of the Mount Roosevelt Formation

contains a significant component of 2500 to 2300 Ma Laurentian basement sourced detrital

zircon not present in any of the intermediate sources. Magmatism of this age has a limited

distribution in western Canada (Figure 1, Canadian Geochronology Knowledgebase 2013). As

discussed above the most likely potential source areas for the 2500 to 2390 Ma grains are in or

near the western Trans-Hudson orogen (Fig. 1). The Buffalo Head, Thorsby or Wabamun

terranes and western edge of the Rae Province form the most likely source areas for 2380 to

2300 Ma grains (Fig. 1; Gehrels and Ross 1998) (Fig. 1). The younger 1825 to 1770 Ma grains

also indicate source areas to the east and southeast as intrusions of this age only occur southeast

23

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 24 of 57

of the Snowbird tectonic zone in the Hearne Province, Trans-Hudson orogen and intervening thrust-fold belt (Fig. 9D; Canadian Geochronology Knowledgebase 2013). The 2050 to 1830 Ma grains match well with plutons in several nearby arcs and terranes and Archean grains could have been derived from the Hearne, Slave or Rae provinces (Fig. 1; Gehrels and Ross 1998).

Both the Vizer and Vampire-Narchilla intervals and the upper member of the Mount Roosevelt

Formation were sourced from the east or southeast. The dramatic change in detrital zircon provenance from recycled Paleoproterozoic blanket sandstone in the Terreneuvian to first cycle

Laurentian basement in the early Middle Cambrian strongly suggests significant erosion of the

Paleoproterozoic sedimentary cover and denudation of Laurentian basement occurred east of the Liard area during late Early Cambrian andDraft Middle Cambrian regression. Further stratigraphic and detrital zircon provenance studies of Ediacaran and Cambrian strata along the northern

Canadian Cordilleran passive margin should be able to refine the evolution of drainage divides and basement denudation.

CONCLUSIONS

The Vizer formation was derived primarily from Paleoproterozoic intermediate sources, with little or no first cycle Laurentian basement input. The conglomeratic middle member of the

Mount Roosevelt Formation was derived primarily from Proterozoic and Lower Cambrian intermediate sources, possibly with a minor component of first cycle Laurentian basement material.

The North American Magmatic Gap aged detrital zircons in the lower part of the Vizer formation were likely derived from a local Mesoproterozoic succession that was completely removed by lower Paleozoic erosion or is preserved in the subsurface of the north central Liard area.

24

https://mc06.manuscriptcentral.com/cjes-pubs Page 25 of 57 Canadian Journal of Earth Sciences

An Early Cambrian Terreneuvian age is suggested for the nearshore Vizer formation based on

similar detrital zircon profiles with Terreneuvian transgressive sandstones to the south, the facies

relationships of Terreneuvian strata in the Mackenzie Mountains and the similarity of the Vizer

formation detrital zircon age distribution with that of the nearby basinal Vampire-Narchilla

Formation.

Volcanic clasts in the middle member of the Mount Roosevelt Formation were likely locally

sourced from the 689.1 ± 4.6 Ma Gataga volcanics (Ferri et al. 1999) and volcanic and/or

intrusive rocks related to the 650 to 640 Ma Pool Creek Syenite (Pigage and Mortensen 2004),

based on the presence of concordant 670 and 640 Ma detrital zircon grains. There is no detrital zircon evidence to support syn-extensionalDraft volcanism associated with the Middle Cambrian development of the Mount Roosevelt Graben.

The detrital zircon provenance of Early and Middle Cambrian strata deposited along the northern

Canadian Cordilleran passive margin indicates major variations in source area and paleoflow

occurred with time. Source areas were to the east when the Vizer, Vampire-Narchilla and

Backbone Ranges formations were deposited (Early Cambrian, Terreneuvian series) and changed

to the northeast as Sekwi Formation and Sekwi correlative strata formed (Early Cambrian, Series

2). The middle member of the Mount Roosevelt Formation (Middle Cambrian) in the Roosevelt

graben was primarily locally sourced, whereas the detritus that formed the upper member was

derived directly from Laurentian basement exposed over a large region to the east or southeast.

The change from recycled Paleoproterozoic sedimentary cover to primary basement sources

suggests that significant denudation of the cratonic basement occurred east of the Liard area in

the Early Cambrian and Middle Cambrian.

25

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 26 of 57

Proterozoic stratigraphic changes associated with the Liard Line indirectly affected the provenance of the Vizer formation and middle member of the Mount Roosevelt Formation.

ACKNOWLEDGEMENTS

This study was funded by the Geological Survey of Canada and the B.C. Ministry of Natural Gas

Development. Discussions with Thomas Hadlari, Tony Hamblin and Larry Lane are greatly appreciated. Wing Chan and Kelsey Mah are thanked for their assistance with data manipulation.

This manuscript benefitted from the constructive comments and reviews of Thomas Hadlari, associate editor Michael Rygel and two anonymous reviewers. Natural Resources Canada ESS contribution number 0160262 . Draft

26

https://mc06.manuscriptcentral.com/cjes-pubs Page 27 of 57 Canadian Journal of Earth Sciences

REFERENCES CITED

Aitken, J,D, 1989. The Sauk sequence – Cambrian to Lower Ordovician miogeocline and

platform. In The Western Canadian Sedimentary Basin. Edited by B.D. Ricketts. Canadian

Society of Petroleum Geologists, pp. 105- 119.

Aitken, J.D. 1993. Cambrian and Lower Ordovician – Sauk Sequence. Sedimentary cover of the

craton in Canada. Edited by D.F. Stott and J.D. Aitken. Geological Survey of Canada,

Geology of Canada Series no. 5, pp. 96-124.

Arnott, R.W.C. and Hein, F.J. 1986. Submarine canyon fills of the Hector Formation, Lake

Louise, Alberta: Late Precambrian, syn-rift deposits of the proto-Pacific miogeocline. Bulletin of Canadian Petroleum Geology,Draft 34 : 395-407. Benjamin, V., Triebel, K, Grewal, I., Hope, R., Kaminski, K., Wilkins, G., Haldane-Wilsone, R.,

Meyer, T., Adanjo, R., Ramsey, D. and Heidari, R. 2011. Technical report and preliminary

economic assessment of the Nonda project, northern BC. Available from

http://www.stikineenergy.com in 2012.

Berkeley Geochronology Center, 2016. Isoplot 3.7 and 4.1:

http://www.bgc.org/isoplot_etc/isoplot.html

Bond, G. C. and Kominz, M. A. 1984. Construction of tectonic subsidence curves for the early

Paleozoic miogeocline, southern Canadian Rocky Mountains--implications for subsidence

mechanisms, age of breakup, and crustal thinning: Geological Society of America,

Bulletin, 95 : 155-173

Campbell, R.B., Mountjoy, E.W. and Struik, L.C. 1982. Structural cross-section through south-

central Rocky and Cariboo mountains to the coast range: Geological Survey of Canada

Open-File Report 831, 1 sheet.

27

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 28 of 57

Canadian Geochronology Knowledgebase, 2013. Geological Survey of Canada, Earth Science

Sector, Natural Resources Canada. http://www.nrcan.gc.ca/earth-sciences/geography/atlas-

Canada/geochron/17012#cgkb

Cecile, M.P., Morrow, D.W. and Williams, G.K. 1997. Early Paleozoic (Cambrian to Early

Devonian) tectonic framework, Canadian Cordillera. .Bulletin of Canadian Petroleum

Geology 45 : 54–74.

Colpron, M., Logan, J.M. and Mortensen, J.K. 2002. U-Pb zircon age constraint for late

Neoproterozoic rifting and initiation of the lower Paleozoic passive margin of western

Laurentia. Canadian Journal of Earth Sciences, 39 : 133-143. Cook, D.G. and MacLean, B.C. 1995. TheDraft intracratonic Paleoproterozoic Forward Orogeny, and implications for regional correlations, Northwest Territories, Canada. Canadian

Journal of Earth Sciences, 32 : 1991–2008.

Cook, D.G. and MacLean, B.C. 2004. Subsurface stratigraphy and tectonics of the western plains

of the Northwest Territories: Geological Survey of Canada, Bulletin 575, 26 PDF sheets,

92 p. text.

Cook, F.A., van der Velden, A.J., Hall, K.W. and Roberts, B.J. 1999. Frozen subduction in

Canada’s Northwest Territories: Lithoprobe deep lithospheric reflection profiling of the

western Canadian Shield. Tectonics, 18 : 1-24.

Davis, W.J., Ootes, Luke, Newton, L., Jackson, V. and Stern, R.A. 2015. Characterization of the

Paleoproterozoic Hottah terrane, WopmayOrogen using multi-isotopic (U-Pb, Hf and O)

detrital zircon analyses: An evaluation of linkages to northwest Laurentian

Paleoproterozoic domains. Precambrian Research, 269 : 296-310.

28

https://mc06.manuscriptcentral.com/cjes-pubs Page 29 of 57 Canadian Journal of Earth Sciences

Dickinson, W.R., Lawton T.F. and Gehrels, G.E. 2009. Recycling detrital zircons: A case study

from the Bisbee Group of southern Arizona. Geology. 37: 503–506.

Evenchick, C.A., Parrish, R.R. and Gabrielse, H. 1984. Precambrian gneiss and late Proterozoic

sedimentation in north-central British Columbia. Geology, 12 : 233-237.

Fallas, K.M., Lane, L.S., and Pigage, L.C., 2014. Geology, La Biche River Yukon-Northwest

Territories; Geological Survey of Canada, Canadian Geoscience Map 144, scale 1:250 00.

doi:10.4095/294606

Fermor, P.R. and Moffat, I.W. 1992. Tectonics and structure of the Western Canada foreland

basin. I n Foreland basins and fold belts. Edited by R.W. Macqueen and D.A. Leckie.

American Association of Petroleum Geologists, Memoir 55, pp. 81-105. Draft Ferri, F., Rees, C., Nelson, J. and Legun, A. 1999. Geology and mineral deposits of the northern

Keckika Trough between Gataga River and the 60 th parallel. British Columbia Ministry of

Energy, Mines and Petroleum Resources, Bulletin 107, 122 p.

Fritz, W.H. 1979. Cambrian stratigraphy in the northern Rocky Mountains, British Columbia.

Geological Survey of Canada, Current Research, Part B, Paper 77-1B, pp. 99-109.

Fritz, W.H. 1991. Cambrian assemblages. In Geological Survey of Canada, Geology of Canada,

no. 4. Edited by H. Gabrielse and C.J. Yorath,.pp. 155-184.

Furlanetto, F., Thorkelson, D., Davis, W.J., Gibson, H.D., Rainbird, R.H. and Marshall,

D.D. 2009. Preliminary results of detrital zircon geochronology, Wernecke

Supergroup, Yukon. In Yukon Exploration and Geology 2008. Edited by L.H. Weston,

L.R., Blackburn and L.L. Lewis. Yukon Geological Survey, pp. 125–135.

29

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 30 of 57

Furlanetto, F., Thorkelson, D.J., Daniel, G.H., Marshall, D.D., Rainbird, R.H., Davis,W.J.,

Crowley, J.L. and Vervoort, J.D. 2013. Late Paleoproterozoic terrane accre-tion in

northwestern Canada and the case for circum-Columbian orogenesis.Precambrian

Research, 224 : 512–528.

Gabrielse .H. 1991. Summary, Cambrian to Middle Devonian assemblages. In Geological Survey

of Canada, Geology of Canada, no. 4. Edited by H. Gabrielse and C.J. Yorath, pp. 153-155.

Gabrielse, H. and Taylor, G.C. 1982. Geological maps and cross sections of the northern

Canadian Cordillera from southwest of Fort Nelson, British Columbia to Gravina island,

southeastern Alaska: Geological Survey of Canada Open-File Report 864, 1 sheet, 15 p. text. Draft Gehrels, G.E. and Pecha, M. 2014. Detrital zircon U-Pb geochronology and Hf isotope

geochemistry of Paleozoic and passive margin strata of western North America:

Geosphere, 10: 49-65, doi:10.1130/GES00889.

Gehrels, G.E., and Ross, G.M. 1998. Detrital zircon geochronology of Neoproterozoic to

Permian miogeoclinal strata in British Columbia and Alberta: Canadian Journal of Earth

Sciences, 35 , 1380–1401, doi:10.1139/e98-071.

Goodfellow, W.D., Cecile, M.P., and Leybourne, M.I. 1995. Geochemistry, petrogenesis, and

tectonic setting of lower Paleozoic alkalic and potassic volcanic rocks, northern Canadian

Cordilleran miogeocline. Canadian Journal of Earth Sciences, 32 : 1236–1254.

Goodge, J.W., Fanning, C.M., Brecke, D.M., Licht, K.J. and Palmer, E.F. 2010. Continua-tion of

the Laurentian Grenville Province across the Ross Sea Margin of EastAntarctica. Journal

of Geology, 118 : 601–619.

30

https://mc06.manuscriptcentral.com/cjes-pubs Page 31 of 57 Canadian Journal of Earth Sciences

Gordey, S.P. 2013. Evolution of the Selwyn Basin region, Sheldon Lake and Tay River map

areas, central Yukon: Geological Survey of Canada, Bulletin 599, 5 sheets, 176 p.

Gower, C.F., and Krogh, T.E. 2002. A U–Pb geochronological review of the Proterozoic history

of the eastern Grenville Province. Canadian Journal of Earth Sciences, 39 : 795–829.

Hadlari, T., Davis, W.J., Dewing, K., Heaman, L.M., Lemieux, Y., Ootes, L., Pratt, B.R., and

Pyle, L.J.,

2012. Two detrital zircon signatures for the Cambrian passive margin of northern Laurentia

highlighted by new U-Pb results from northern Canada: Geological Society of America

Bulletin, 124 , 1155–1168, doi:10.1130 /B30530.1. Heaman, L.M., LeCheminant, A.N. andDraft Rainbird, R.H. 1992. Nature and timing of Franklin igneous events, Canada: Implications for a late Proterozoic mantle plume and the break-up

of Laurentia: Earth and Planetary Science Letters, 109: 117-131.

Heaman, L.M., Gower, C.F. and Perreault, S. 2004. The timing of Proterozoic magmatism in the

Pinware terrane of southeast Labrador, easternmost Quebec and northwest Newfoundland.

Canadian Journal of Earth Sciences, 41: 127-150.

Hoffman, P.F. 1989. Precambrian geology and tectonic history of North America. In Bally,

A.W., Palmer, A.R. (Eds.), The Geology of North America — an overview. Edited by

A.W. Bally and A.R. Palmer. Geology of North America, A. Geological Society of

America, pp. 447–512.

Karlstrom, K.E., Åhäll, K.-I., Harlan, S.S., Williams, M.L., McLelland, J. and Geissman,J.W.

2001. Long-lived (1.8–1.0 Ga) convergent orogen in southern Laurentia, its extensions to

Australia and Baltica, and implications for refining Rodinia.Precambrian Research. 111: 5–

30.

31

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 32 of 57

Lane, L.S. and Gehrels, G.E. 2014. Detrital zircon lineages of late Neoproterozoic and Cambrian

strata, NW Laurentia: Geological Society of America Bulletin, 126: 398–414, doi:

10.1130/B30848.1.

Lechimant, A.N. and Heaman, L.M. 1989. Mackenzie igneous events, Canada: Middle

Proterozoic hotspot magmatism associated with ocean opening. Earth and Planetary

Science Letters, 96 : 38-48.

Leslie, C.D. 2009. Detrital zircon geochronology and rift-related magmatism: central Mackenzie

Mountains, Northwest Territories. M.Sc. thesis: University of British Columbia, 236 p.

Macdonald, F. A., Schmitz, M. D., Crowley, J. L., Roots, C. F., Jones, D. S., Maloof, A. C.,

Strauss, J. V., Cohen, P. A., Johnston, D. T. and Schrag, D. P. 2010. Calibrating the

Cryogenian. Science, 327: 1241-1243.Draft

MacNaughton, R.B., Roots, C.F. and Martel, E. 2008. Neoproterozoic-(?)Cambrian

lithostratigraphy, northeast Sekwi Mountain map area, Mackenzie Mountains, Northwest

Territories: new data from measured sections. Geological Survey of Canada, Current

Research (Online) no. 2008-16, 17 p.

Mair, J.L., Hart, C.J.R. and Stephens, J.R. 2006. Deformation history of the northwestern Selwyn

Basin, Yukon, Canada: Implications for orogen evolution and mid-Cretaceous magmatism:

Geological Society of America Bulletin, 118: 304-323.

Matthews, W.A. Detrital zircon geochronology with applications to North American tectonics.

Ph.D. thesis: University of Calgary, 329 p.

Medig, K.P.R., Thorkelson, D.J. and Dunlop, R.L. 2010. The Proterozoic Pinguicula Group:

Stratigraphy, contact relationships and possible correlations. In Yukon Exploration and

32

https://mc06.manuscriptcentral.com/cjes-pubs Page 33 of 57 Canadian Journal of Earth Sciences

Geology 2009. Edited by K.E. MacFarlane, L.H. Weston and L.R. Blackburn. Yukon

Geological Survey, pp. 265-278.

Medig, K.P.R., Thorkelson, D.J., Turner, E.C., Davis, Gibson, H.D., Rainbird, R.H., and

Marshall, D.D. 2012. The Proterozoic Pinguicula Group, Wernecke Mountains, Yukon: A

siliciclastic and carbonate slope to basin succession with local and exotic sediment

provenance. In Yukon Exploration and Geology 2011. Edited by K.E. MacFarlane and P.J.

Sack. Yukon Geological Survey, pp. 129-149.

Medig, K.P.R., Thorkelson, D.J., Davis, Rainbird, R.H., Gibson, H.D., Turner, E.C., and

Marshall, D.D. 2014. Pinning northeastern Australia to northwestern Laurentia in the

Mesoproterozoic: Precambrian Research, 249 : 88-99.

McMechan, M.E. 1994. Geology, DawsonDraft Creek, west of the sixth meridian, British Columbia;

Geological Survey of Canada, Map 1858A, scale 1:250,000, with balanced structural cross-

section, 1 sheet.

McMechan, M.E. 2000. Walker Creek Fault Zone, central Rocky Mountains, British Columbia -

southern continuation of the Northern Rocky Mountain Trench Fault Zone: Canadian

Journal of Earth Sciences, 37: 1259-1273.

McMechan, M.E. 2002. Structural geometry in the Carbon Creek area of the Rocky Mountain

Fold and Thrust Belt: Bulletin of Canadian Petroleum Geology, 50: 407- 418.

McMechan, M.E. 2013. Geology, Kakwa area, British Columbia – Alberta: Geological Survey of

Canada Map 2065A, scale 1:125,000, with balanced structural cross-sections, 3 sheets,

doi:10.4095/293062

33

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 34 of 57

McMechan, M., Ferri, F. and MacDonald, L. 2012. Geology of the Toad River area (NTS 094N),

northeast British Columbia. In Geoscience Reports 2012. British Columbia Ministry of

Energy and Mines, pp. 17–39.

Milton, J.E. 2015. Sedimentary rock hosted copper mineralization in the Neoproterozoic

Redstone Copperbelt, Mackenzie Mountains, Northwest Territories, Canada. Ph.D. thesis,

Department of Geological Siences, University of British Coulumbia, Vancouver. 396 p.

Nelson, J.L., Colpron, M. and Israel, S, 2013. The Cordillera of British Columbia, Yukon and

Alaska : tectonics and metallogeny : Society of Economic Geologists, Special Publication

17, p. 53-109.Omma, J.E., Pease, V., and Scott, R.A., 2011, U–Pb SIMS zircon geochronology of Triassic and JurassicDraft sandstones on northwestern Axel Heiberg Island, northern Sverdrup Basin, Arctic Canada. In Arctic Petroleum Geology. Edited by A. M.

Spencer, , A. F. Embry, , D. L. Gautier, , A. V. Stoupakova , and K. Sørensen. Geological

Society, London, Memoir, 35, pp. 559–566.

Ootes, L. Davis, W. J., Jackson, V. A. and van Breemen O. 2015. Chronostratigraphy of the

Hottah terrane and Great Bear magmatic zone of Wopmay Orogen, Canada, and

exploration of a terrane translation model. Canadian Journal of Earth Sciences, 52 : 1062-

1092, 10.1139/cjes-2015-0026

Pigage, L.C. 2009. Bedrock geology of NTS 95C/5 (Pool Creek) and NTS 95D/8 map sheets,

southeast Yukon. Yukon Geological Survey, Bulletin 16, 150 p.

Pigage, L.C. and Mortensen, J.K. 2004. Superimposed Neoproterozoic and early Tertiary

alkaline magmatism in the La Biche River area, southeast Yukon Territory. Bulletin of

Canadian Petroleum Geology, 52 : 325-342.

34

https://mc06.manuscriptcentral.com/cjes-pubs Page 35 of 57 Canadian Journal of Earth Sciences

Pigage, L.C., Roots, C.F. and Abbott, J.G. 2015. Regional bedrock geology for Coal River map

area (NTS 95D), southeast Yukon: Yukon Geological Survey, Bulletin 17, 155 p.

Post, R.T. and Long, D.G.F. (2008). The Middle Cambrian Mount Roosevelt Formation (new) of

northeastern British Columbia: evidence for rifting and development of the Kechika graben

system; Canadian Journal of Earth Sciences, 45 : 483–498.

Price, R.A. and Fermor, P.R. 1984. Structure section of the cordilleran foreland thrust and fold

belt west of Calgary, Alberta: Geological Survey of Canada, Paper 84-14, 1 sheet.

Pugh, D.C. 1975. Cambrian stratigraphy from western Alberta to northeastern British Columbia:

Geological Survey of Canada, Paper 74-37, 10 sheets, 31 p.

Pullen, A., Ibáñez-Mejía, M., Gehrels, G.E.,Draft Ibáñez-Mejíacd, J.C. and Pecha, M. 2014. What happens when n= 1000? Creating large-n geochronological datasets with LA-ICP-MS for

geologic investigations. Journal of Analytical Atomic Spectrometry, 29: 971-980. DOI:

10.1039/C4JA00024B

Rainbird, R.H. and Davis, W.J. 2007. U-Pb detrital zircon geochronology and provenance of the

Late Paleoproterozoic Dubawnt Supergroup: linking sedimentation with tectonic

reworking of the western Churchill Province, Canada: Geological Society of America

Bulletin, 119 , 314-328.

Rainbird, R.H., Jefferson, C.W. and Young, G.M. 1996. The early Neoproterozoic sedimentary

Succession B of northwestern Laurentia: Correlations and paleogeographic significance:

Geological Society of America Bulletin, 108 : 454-470.

Rainbird, R.H., McNicoll, V.J., Thériault, R.J., Heaman, L.M., Abbott, J.G., Long, D.G.F. and

Thorkelson, D.J. 1997. Pan-continental river system draining Grenville orogen recorded by

35

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 36 of 57

U-Pb and Sm-Nd Geochronology of Neoproterozoic quartzarenites and mudrocks,

northwestern Canada: The Journal of Geology, 105 : 1-17.

Rainbird, R.H., Stern, R.A., Rayner, N. and Jefferson, C.W. 2005. Age, provenance, and regional

correlation of the Athabasca Group, Saskatchewan and Alberta, constrained by igneous and

detrital zircon geochronology. In EXTECH IV: Geology and Uranium EXploration

TECHnology of the Proterozoic Athabasca Basin, Saskatchewan and Alberta. Edited by

C.W. Jefferson and G. Delaney. Geological Survey of Canada, Bulletin 588, pp. 1-16.

Rainbird, R.H., Cawood, P.A. and Gehrels, G.E. 2012. The Great Grenvillian sedimentation

episode: record of Supercontinent Rodinia’s assembly. In Recent Advances in Tectonics of

Sedimentary Basins. Edited by C. Busby and A. Azor. Blackwell Publishing Ltd., pp. 583-

601. Draft

Rayner, N.M. and Rainbird, R.H. 2013. U-Pb geochronology of the Shaler Supergroup, Victoria

Island, northwest Canada: 2009-2013. Geological Survey of Canada, Open File 7419, 62 p.

doi:10.4095/292694

Rayner, N., Stern, R.A. and Rainbird, R.H. 2003. SHRIMP U-Pb detrital zircon geochronology

of Athabasca Group sandstones, northern Saskatchewan and Alberta; Geological Survey of

Canada, Current Research 2003-F2, 20 p.

Rogers, J. and Santosh, M. 2002. Configuration of Columbia, a Mesoproterozoic super-

continent. Gondwana Research. 5: 5–22.

Rooney, A. D., Macdonald, F. A., Strauss, J. V., Dudás, F. Ö., Hallmann, C. and Selby, D. 2014.

Re-Os geochronology and coupled Os-Sr isotope constraints on the Sturtian snowball

Earth: Proceedings of the National Academy of Sciences, 111 : 51-56.

36

https://mc06.manuscriptcentral.com/cjes-pubs Page 37 of 57 Canadian Journal of Earth Sciences

Ross, G.M., McMechan, M.E. and Hein, F.J. 1989. Proterozoic history: Birth of the miogeocline.

In The Western Canadian Sedimentary Basin. Edited by B.D. Ricketts. Canadian Society

of Petroleum Geologists, pp. 79-104.

Ross, G.M., Bloch, J.D. and Krouse, H.R. 1995. Neoproterozoic strata of the southern Canadian

Cordillera and the isotopic evolution of seawater sulfate. Precambrian Research, 73 : 71-99.

Ross, G.M., Villeneuve, M.E. and Theriault, R.J. 2001. Isotopic provenance of the lower

Muskwa assemblage (Mesoproterozoic, Rocky Mountains, British Columbia); new clues to

correlation and source areas: Precambrian Research, 111 : 57-77.

Ross, G. M. and Villeneuve, M. 2003. Provenance of the Mesoproterozoic (1.45 Ga) Belt basin

(western North America): AnotherDraft piece in the pre-Rodinia paleogeographic puzzle:

Geological Society of America Bulletin, 115 : 1191-1217.

Slind , O. L., Andrews, G. D., Murray, D. L., Norford, B. S., Paterson, D. F., Salas, C. J.,

Tawadros, E. E. and Aitken, J. D. 1994. Chapter 8 – Middle Cambrian to Lower

Ordovician strata of the Western Canada Sedimentary Basin. In Geological atlas of the

Western Canada sedimentary basin. Compiled by G.D. Mossop and I Shetson. Canadian

Society of Petroleum Geologists and Alberta Research Council, pp. 87-108.

Stewart, E.D., Link, P.K., Fanning, C.M., Frost, C.D. and McCurry, M. 2010. Paleogeographic

implications of non-North American sediment in the Mesoproterozoic upper Belt

Supergroup and Lemhi Group, Idaho and Montana. Geology, 38: 927-930, doi:

10.1130/G31194.1

Stewart, J.H. 1972. Initial deposits in the Cordilleran geosyncline; Evidence of Late Precambrian

37

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 38 of 57

(<850m.y.) Continental Separation. Geological Society of America Bulletin, 83 : 1345-

1360.

Taylor, G.C. and Stott. D.F. 1973. Tuchodi Lakes map-area, British Columbia. Geological

Survey of Canada, Memoir 373, 37 p..

Taylor, G.C. and Stott, D.F. (1980). Geology, Toad River, British Columbia. Geological Survey

of Canada, Open File 674, scale 1:250 000.

Thorkelson, D.J., Abbott, J.G., Mortensen, J.K., Creaser, R.A., Villeneuve, M.E., McNicoll,V.J.

and Layer, P.W. 2005. Early and middle Proterozoic evolution of Yukon, Canada.

Canadian Journal of Earth Sciences, 42 : 1045–1071.

Vermeesch, P. 2013. Multi-sample comparison of detrital age distributions. Chemical Geology

341: 140-146, doi:10.1016/j.chemgeo.2013.01.010Draft

Vermeesch, P. 2017. IsoplotR, http://isoplotr.london-geochron.com/

Whitmeyer, S.J., Karlstrom, K.E. 2007. Tectonic model for the Proterozoic growth of North

America. Geosphere, 3: 220-259.

Young, G. M. 1978. Proterozoic (< 1.7 b.y.) stratigraphy, paleocurrents and orogeny in North

America. Egyptian Journal of Geology, 22: 45-64.

Yukon Geochronology Database, 2015. Geochronology (Rock Ages) for the Yukon.

http://www.geology.gov.yk.ca/databases_gis.html#Databases

38

https://mc06.manuscriptcentral.com/cjes-pubs Page 39 of 57 Canadian Journal of Earth Sciences

Potential Intermediate Reference Sources Athabasca Group Rayner et al. (2003) Horsethief Creek Group Gehrels and Pecha (2014) Mackenzie Mountain Leslie (2009) Supergroup (Katherine +Little Dal Grps) Keele Fm. (L+G) Lane and Gehrels (2014) Keele Fm. (L) Leslie (2009) Muskwa Group Ross et al. (2001) Pa Pigage (2009) Ps Pigage (2009) Pinguicula Group Medig et al. (2012) Rapitan Group (L+G) Lane and Gehrels (2014) Toobally Formation Pigage (2009) Wernecke Supergroup Furlanetto et al. (2016) Yusezyu Formation Pigage et al. (2015)

Cambrian Units for Comparison Backbone Ranges Fm. Leslie (2009) Crow Formation Pigage (2009) Franklin Mountains Fm. Leslie (2009) Draft Hamill Group Gehrels and Pecha (2014) McNaughton Formation Matthews (2016) Mt. Roosevelt Fm. (middle) This study Mt. Roosevelt Fm. (upper) Gehrels and Pecha (2014) Sekwi Fm. correlative Pigage et al. (2015) Sekwi Formation Leslie (2009) Vampire-Narchilla unit Pigage et al. (2015) Vizer formation This study Yanks Peak Formation Matthews (2016)

McMechan et al Cambrian DZ Table 1

Table 1. References for detrital zircon age data used in Figures 4 - 6, 8. Fm. – formation, Grps. –

Groups.

39

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 40 of 57

FIGURE CAPTIONS

Figure 1. Subcrop map at the base of the Paleozoic for northwest Laurentia showing distribution of Cambrian sandstone, and Cambrian and Proterozoic detrital zircon localities palispastically restored for Cordilleran deformation. Numbers correspond to sections in Figure 2. Sources of information: Aitken 1993, Cook et al. 1999; Canadian Geochronology Knowledgebase 2013;

Cook and MacLean 2004; McMechan 2002; Ootes et al. 2015. Pugh 1975; Ross et al. 2001;

Slind et al. 1994; unpublished B.C. well logs. Palinspastic restoration based on: Campbell et al.

1982; K.M. Fallas, pers. comm. 2015; Fermor and Moffat 1992; Gabrielse and Taylor 1982;

Gordey 2013; Mair 2006; McMechan 1994, 2000, 2013, unpublished; Pigage et al. 2015; Price and Fermor 1984. BP – Belt-Purcell basin;Draft C – Coronation Margin; G – Great Bear magmatic zone; H – Hamill Formation and Horsethief Creek Group; HBB – Horby Bay Basin; HCI – Hart

Creek Inlier; MM – Mackenzie Mountain Supergroup; Mn – McNaughton Formation; P –

Pinware terrane (inset map); PR1 = PR1 Basin; S-l – upper Coppermine River Group and lower

Shaler Supergroup; S-u – upper Shaler Supergroup; T – Thorsby terrane; W – Wabamun terrane;

Yp – Yanks Peak Formation.

Figure 2. Proterozoic to Middle Ordovician stratigraphic correlations and tectonic events along the northern Canadian passive margin showing stratigraphic positions of detrital zircons samples from this study or used for comparison purposes in Figures 4 to 8. See Table 1 for sources of detrital zircon U-Pb age data. FO – Forward Orogeny; GB – Gunbarrel magmatic event; RO –

Racklan Orogeny.

Figure 3. U-Pb zircon probability density distributions and histograms for Vizer and Mount

Roosevelt formation samples.

40

https://mc06.manuscriptcentral.com/cjes-pubs Page 41 of 57 Canadian Journal of Earth Sciences

Figure 4. Normalized U-Pb zircon age probability plots for Ediacaran to Cambro-Ordovician

units in the Canadian Cordilleran passive margin succession. See Table 1 for sources of detrital

zircon U-Pb age data from other studies. Grey bars highlight age range of ultimate Grenville

Orogeny source area (Gren), NAMG (N) and 2500 to 2300 Ma. Edi – Ediacaran; eC2 – Early

Cambrian – series 2; eCT – Early Cambrian - Terrenuevian; lC- eO – Late Cambrian to Early

Ordovician; lC- mO – Late Cambrian to Middle Ordovician; mC – Middle Cambrian.

Figure 5. MDS plots (Vermeesh 2017) showing relative statistical similarity of detrital zircon age

distributions (based on the Kolmogorov-Smirnov test D value) for Ediacaran and Cambro-

Ordovician units in the Canadian Cordilleran passive margin succession. See Table 1 for sources of detrital zircon U-Pb age data. BBl – Draftlower Backbone Ranges Formation, BBu – upper Backbone Ranges Formation. The thick pale grey dashed line separates samples with ultimate

Grenville sourced grains (to the right of the line) from those lacking or having very few

Grenville sourced zircons (to the left of the line).

Figure 6. Normalized U-Pb zircon age probability plots for potential sources of zircons found in

Mount Roosevelt and Vizer formation samples. See Table 1 for sources of detrital zircon U-Pb

age data. Grey bars highlight age range of ultimate Grenville Orogeny source area (Gren),

NAMG (N) and 2500 to 2300 Ma. eCT – Early Cambrian – Terrenuevian; mC – Middle

Cambrian; mP – Mesoproterozoic; nP – Neoproterozoic; pP – Paleoproterozoic. MM –

Mackenzie Mountain Supergroup; l Shaler – upper Coppermine River Group and lower Shaler

Supergroup; L – Leslie (2009); LG – Lane and Gehrels (2014).

41

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 42 of 57

Figure 7. U-Pb zircon age probability density distributions from samples of the basal Pinguicula

Group (unit A) exposed in the Wernecke Inlier, Wernecke Mountains, Yukon (A – C; Medig et al. 2012) and the PR1 basin (D; Medig et al. 2014). Grey bar indicates age range of NAMG.

Figure 8. MDS plots (Versmeech, 2017) showing relative statistical similarity of detrital zircon age distributions (based on the Kolmogorov-Smirnov test D value) for potential intermediate sources of detrital zircons found in the Vizer and Mount Roosevelt formation. Samples from this study are labelled in italics. See Table 1 for sources of detrital zircon U-Pb age data. HTC –

Horsethief Creek Group; L – Leslie (2009); LG – Lane and Gehrels (2014); MM – Mackenzie

Mountain Supergroup; Shaler- l – upper Coppermine River Group and lower Shaler Supergroup; Shaler- u – Shaler Supergroup, upper part.Draft Note: the lower part of the Shaler Supergroup includes one sample from sandstone at top of Coppermine River Group.

Figure 9. Hypothetical drainage pathways suggested by observed detrital zircon signatures superimposed on sub-Paleozoic subcrop map (see Figure 1 for details). A) Vampire-Narchilla

(VN), Vizer (V) and Backbone Ranges (BR) formations (Early Cambrian – Terreneuvian); B)

Sekwi Formation (S) and Sekwi correlative unit (SC) (Early Cambrian – Series 2); C) middle member Mount Roosevelt Formation (Middle Cambrian – Series 3); D) upper member Mount

Roosevelt Formation (Middle Cambrian – Series 3). BHT – Buffalo Hills terrane; C –

Coronation Margin; G – Great Bear magmatic zone; H – Hottah terrane; LR – La Ronge arc; MA

– Mackenzie arch; MM- Mackenzie Mountains Supergroup; PRA – Peace River Arch; STZ –

Snowbird Tectonic zone; T – Thorsby terrane; TA – Taltson magmatic belt; W – Wabamun terrane.

42

https://mc06.manuscriptcentral.com/cjes-pubs Page 43 of 57 Canadian Journal of Earth Sciences

‘DISMAL SHALER LAKES’ S-l QUEEN ArcticS-u Islands MAUD 1 ‘HORNBY 10 BLOCK PR1 HBB BAY’ 2.1-1.85 Ga o 2.85-2.00 Ga Liard Area W P WERNECKE C o N 1a G 65 p HCI DUBAWNT a MM 1.87-1.85 Ga (THELON ra SLAVE Laurentia -a BASIN) u 1b PROVINCE to 2.5-1.89 Ga c h to n W o o u s ? RAE t 140 e PROVINCE r r a n e s AH 2 o 3a 60 N HOTT 3b LIARD LINE 1.94-1.90 Ga TALTSON ATHABASCA 4 Draft Liard (BASIN) MUSKWA BUFFALO Area HEAD TERRANE 2.00-1.92 Ga 2.32-1.99 Ga

1.96-1.90 Ga o N 55 TRANS- HUDSON Ea st er OROGEN n lim T it o f 2.32 Ga W M <2.3 Ga id d 2.3-1.91 Ga l YP e C Mn a m b 1.85-1.78 Ga HEARNE ri an PROVINCE s tra ta 50o N

1

10 W o H BP o W

W

300 km o 130

120

CORDILLERAN PASSIVE MARGIN 1.9-1.8 Ga juvenile crust with 2.32-1.90 Ga accreted Proterozoic crust and Lower Cambrian sandstones ‘basement’ inliers Neoproterozoic - Earliest Cambrian 2.0-1.78 Ga continental magmatic arcs Archean (>2.5 Ga) crust Windermere Supergroup, Hyland Group Neoproterozoic Tonian 1.9-1.80 Ga thrust-fold belts Cambrian detrital zircon sample locality, Mackenzie Mountains/ Shaler Supergroups this study; other studies Mesoproterozoic Sedimentary Sequences 2.3 Ga magmatic crust Proterozoic detrital zircon sample locality other studies

Paleoproterozoic 1760 - 1600 Ma Eastern limit Cordilleran deformation Magmatism: 2.4 to 2.3 Ga, 2.5 to 2.4 Ga Sedimentary Sequences

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 44 of 57

N. Mackenzie - Coal River NW Liard area SW Liard area Tectonic Events Wernecke area 450 1a Mountains 1b 2 3a 3b 4 area 1 areas 2-4 Sunblood Sunblood Fm. Sunblood Crow Fm. Franklin Rabbitkettle Crow Fm. Kechika Fm. Ordovician Mountain Fm. Fm. u Otter Ck Fm. 500 m Mt. Roosevelt

Sekwi Fm. Sekwi - correlative ? Cambrian l Vam- Backbone unnamed formation Cordilleran pire Vizer fm. quartzite Ranges Fm. Vampire - passive margin Narchilla unit 550 Risky Fm. base not exposed

Blueflower Fm. Hyland Group Draft Yusezyu Fm. Gametrail Fm. Ediacaran

Sheepbed Fm. base not exposed Toobally Ice Brook Fm. Fm. 650 Keele Fm. + + 650 Ma ? Pool Creek Windermere Supergroup Twitya Fm. unit Pa Syenite formation Cryogenian ? Rapitan Gp proto-Pacific margin? Neoproterozoic

750 CL GB

Mackenzie Mountains onian

T Supergroup 1000 base not exposed ? Pinguicula unit Ps Group ? base not 1350 exposed PR1 Mesoproterozoic

1600 ? RO WSG Muskwa Assemb. base not ? base not exposed FO 1750 exposed

Ma Paleo- proterozoic

Unit with samples analysed for Other studies for comparison, see Extension, rifting U-Pb detrital zircon ages and Figs. 4-6, 8 lacking Grenville age zircons See Medig et al. (2012, 2014) and This study Folding, reverse faulting Fig. 7

https://mc06.manuscriptcentral.com/cjes-pubs Page 45 of 57 Canadian Journal of Earth Sciences

40 FF10-1 Vizer formation 35 11EPF55A Vizer formation 11EPF67A Mount Roosevelt 25 35 Cariboo Range Cariboo Range (n = 99) (n = 94) (n = 96) Formation Relative probability 30 (middle mbr.) 30 20 25 25 15 20 20 15 15 10 Number of Grains 10 10 5 5 5 0 0 0 0 500 1000 1500 2000 2500 3000 3500 0 500 1000 1500 2000 2500 3000 3500 0 500 1000 1500 2000 2500 3000 3500 Age (Ma) Age (Ma) Age (Ma)

Draft

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 46 of 57

Unit THIS STUDY Mt. Roosevelt Fm. (middle) Age n=95 mC THIS STUDY Vizer formation n=194 ? Crow Formation NW Liard loc.(3A) n=264 lC- mO Crow Formation Coal River loc. n = 271 lC- mO Franklin Mtn. Fm. uC- n=45 Draft eO Mt. Roosevelt Formation (upper) n = 201 mC Sekwi Formation n=60 eC2 Sekwi Correlative n=74 eC2 upper Backbone Range Formation n=57 eCT lower Backbone Range Formation n=48 eCT Vampire - Narchilla unit n = 230 eCT McNaughton Formation n=222 eCT

Hamill Group n=192 eCT

Yanks Peak Formation n=112 Edi GREN N 500 1000 1500 2000 2500 3000 3500 4000 Ma

https://mc06.manuscriptcentral.com/cjes-pubs Page 47 of 57 Canadian Journal of Earth Sciences

Mount Roosevelt (upper) This Study 0.2 Crow (NW Liard) BBl BBu Franklin Mountain Crow (Coal River) 0.0 V-N Sekwi correlative Sekwi Mt. Roosevelt H (middle) McNaughton Unlikely correlatives of Vizer and Mount -0.2 Roosevelt formations: Vizer units contain abundant DraftGrenville age zircons Yanks Peak

-0.2 0 0.2 0.4 0.6

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 48 of 57

Unit Age upper Mt. Roosevelt n = 201 mC middle Mt. Roosevelt THIS STUDY n = 95 mC

THIS STUDY Vizer formation n = 194 eCT

Yusezyu n = 114 nP

Muskwa n = 102 mP Ps Draft n = 229 pP

Athabasca n = 278 pP Wernecke n = 275 pP

Horsethief Creek n = 388 nP Toolbally n = 63 nP Keele (L) n = 55 nP TION SOURCE Keele (LG) n = 99

Wiindermere Supergroup nP Pa n = 98 nP Rapitan n = 96 nP upper Shaler

n = 330 TIBLE VIZER FORMA

nP A MM n = 101 nP l Shaler n = 257 INCOMP nP-mP GREN N 500 1000 1500 2000 2500 3000 3500

https://mc06.manuscriptcentral.com/cjes-pubs Page 49 of 57 Canadian Journal of Earth Sciences

A KM09-1-4-1C 0.005 Pinguicula Grp. Wernecke inlier n= 53/77 0.003 Probability

0.001

900 1300 1700 2100 2500 2900 3300 Age (Ma) B KM09-10-3-1C Pinguicula Grp. 0.003 Wernecke inlier n= 66/82

Probability 0.001

900 1300 1700 2100 2500 2900 3300 Age (Ma) Draft C KM09-12-2-2C 0.005 Pinguicula Grp. Wernecke inlier n= 31/47 0.003 Probability

0.001

900 1300 1700 2100 2500 2900 3300 Age (Ma)

D Basal Fifteen Mile Group PR1 basin n= 83 Relative Probability

900 1300 1700 2100 2500 2900 3300 Age (Ma)

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 50 of 57

0.4

Mt. Roosevelt (upper) Likely Unlikely Sources: 0.2 Sources units contain abundant Grenville Wernecke age zircons Rapitan Athabasca Shaler- l Keele (LG) Ps 0 Muskwa Pa Mt.Roosevelt Shaler-u MM HTC (middle) Keele (L) Toobally Yusezyu -0.2 Viser

-0.2 0 0.2 0.4 Cambrian Draftunit Potential intermediate source Cambrian unit This Study

https://mc06.manuscriptcentral.com/cjes-pubs Page 51 of 57 Canadian Journal of Earth Sciences

SHALER SHALER

C o o 65 N MA C 65 N MM SLAVE H G SLAVE BR PROVINCE S MM PROVINCE W W o H G o ? 140 140 VN Sc

LIARD LINE LIARD LINE V ATHABASCA ATHABASCA MUSKW BASIN BASIN 300 km 300 km A

PRA PRA

W W

o o

130 130 o o 50 N 50 N

1 1

10 10

Vizer formation o Sekwi Formation o A and correlatives W B and correlatives W

SHALER SHALER

o o 65 N 65 N MM Draft MM SLAVE SLAVE PROVINCE PROVINCE W W o o ? ? RAE 140 140 PROVINCE

LIARD LINE LIARD LINE ATHABASCA TA ATHABASCA BASIN BASIN 300 km 300 km BHT STZ PRA PRA TRANS- HUDSON W OROGEN

W W

o o T HEARNE PROVINCE

130 130 o o 50 N 50 N

1 1

10 10

Middle Member Mount o Upper Member Mount o C Roosevelt Formation W D Roosevelt Formation W

Neoproterozoic - Earliest Cambrian Mesoproterozoic Sedimentary Inferred extent of Paleoproterozoic Windermere Supergroup, Hyland Group Sequences sedimentary cover

Neoproterozoic Tonian Paleoproterozoic 1760 - 1600 Magmatism: 2.4 to 2.3 Ga, 2.5 to 2.4 Ga Mackenzie Mountains/ Shaler Supergroups Ma Sedimentary Sequences Detrital zircon sample

https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 52 of 57

1

Zircon U-Pb Dating Methodology

Published Methods

Detailed descriptions of the methods followed by Geosep Services to produce and process their ZrnUPb data have been presented in numerous peerreviewed manuscripts. These include Bradley et al. (2009), Hults et al. (2013), and Moore et al. (2015).

Sample Preparation

Zircon grains were isolated and prepared for LAICPMS analysis using standard procedures combined with specific customized procedures described by Donelick et al. (2005). These customized procedures were designed to maximize recovery of: 1) all possible grain sizes present within a sample by minimizing the potential loss of smaller grain sizes through the use of watertable devices, and 2) complete grains with as close to full terminations by minimizing grain breakage and/or fracturing inherent with the standard procedures typically used to separate individual grains from the original rock material. Use of these procedures results in a significantly greater range of recovered grain sizes, as well as a higher percentage of “complete” grains being retained during the mineral separation process (Fig. 1). Draft

Fig. 1. Example zircon separate from GSS Project 025, showing a wide range in grainshapes and grainsizes after sample processing by GSS mineral separation procedures. Importantly, a significant number of tiny grains, as well as complete grains were recovered.

Whole rock samples were first run multiple (minimum = 3) times through a Chipmunk brand jaw crusher with the minimum jaw separation set to 23 mm. The crushed material was then sieved through 300 µm nylon mesh, and the <300 µm size fraction washed with tap water and allowed to dry at room temperature. Zircon grains were separated from other mineral species using a combination of lithium metatungstate (density ~2.9 g/cm 3), Frantz magnetic separator,

 2016, GeoSep Services, LLC https://mc06.manuscriptcentral.com/cjes-pubs Page 53 of 57 Canadian Journal of Earth Sciences

2

diiodomethane (density ~3.3 g/cm 3), and handpanning separation procedures. Epoxy wafers (~1 cm x 1 cm) containing zircon grains for LAICPMS were polished manually using 3.0 µm and 0.3 µm Al 2O3 slurries to expose internal zircon grain surfaces. The polished zircon grain surfaces were washed in 5.5 M HNO 3 for 20 s at 21ºC in order to clean the grain surfaces prior to introduction into the laser system sample cell.

LA-ICP-MS Session Details

LAICPMS data collection was performed at the Geoanalytical Laboratory, Washington State University, Pullman, Washington, U.S.A following conditions and parameters presented in Table 1. Individual zircon grains were targeted for data collection using a New Wave YP213 213 nm solid state laser ablation system using a 20 µm diameter laser spot size, 5 Hz laser firing rate, and ultra high purity He as the carrier gas. Isotopic analyses of the ablated zircon material were performed using a ThermoScientific Element2 magnetic sector mass spectrometer using high purity Ar as the plasma gas. The following masses (in amu) were monitored for 0.005 s each in pulse detection mode(Pb, Th, and U isotopes): 202, 204, 206, 207, 208, 232, 235, and 238.

Table 1. ICPMS and laser ablations system operating conditions and data acquisition parameters

ICP-MS: operating conditions Draft Instrument Finnigan Element II Magnetic Sector ICPMS Forward power 1.25 kW Reflected power <5 W Plasma gas Ar Coolant flow 15 l/min Carrier flow 1.0 l/min (Ar) 0.8 l/min (He) Auxiliary flow 1.0 l/min ICP-MS: data acquisition parameters Dwell time 24 milliseconds per peak point Points per peak 3 Mass window 5% Scans 200 Data acquisition time 29.5 sec Data acquisition mode E scanning Isotopes measured 43 Ca or 29 Si and 238 U Laser ablation system: operating conditions Laser type New Wave Neodymium: YAG Wavelength 213 nm Laser mode Q switched Laser output power 10 J/cm 2 Laser warm up time 6 sec Shot repetition rate 5 Hz Sampling scheme spot (20 m)

At time = 0.0 s, the mass spectrometer began monitoring signal intensities; at time = 6.0 s, the laser began ablating zircon material; at time = 30.0 s, the laser was turned off and the mass

 2016, GeoSep Services, LLC https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 54 of 57

3 spectrometer stopped monitoring signal intensities. A total of 200 data scans were collected for each zircon spot analyzed comprising: approximately 55 background scans; approximately 20 transitions scans between background and background+signal, approximately 125 background+signal scans. A scheme was developed to check whether mass 238 experienced a switch from pulse to analog mode during data collection and a correction procedure was employed to ensure the use of good quality intensity data for masses 235 and 238 when such a switch was observed.

UPb Data Analysis

Previous LAICPMS studies of UPb zircon dating used the socalled intercept method, which assumes that isotopic ratio varies linearly with scan number due solely to linearly varying isotopic fractionation (Chang et al., 2006; Gerhels et al., 2008). The data modeling approach favored here was the modeling of backgroundcorrected signal intensities for each isotope at each scan. Background intensity for each isotope was calculated using a fitted line (for decreasing background intensity) or using the arithmetic mean (for nondecreasing background intensity) at the global minimum of selected isotopes ( 206 Pb, 232 Th, and 238 U) for the spot. Background+signal intensity for each isotope at each scan was calculated using the median of fitted (2 nd order polynomial) intensity values for a moving window (7 scans wide here) that includes the scan. The precision of each backgroundcorrected signal intensity value was calculated from the precision of backgroundDraft intensity value and the precision of the background+signal intensity value.

Zircon UPb age standards used during analysis are summarized in Table 2, including the 1099±0.6 Ma FC zircon (FC1 of Paces and Miller, 1993) used here as the primary age standard. Isotopic data for FC were used to calculate Pb/U fractionation factors and their absolute errors for each FC data scan at each FC spot; these fractionation factors were smoothed sessionwide for each data scan using the median of fitted (1 st order polynomial) fractionation factor values for a moving window that includes the current FC spot and scan.

Table 2. Zircon age standards.

Standard Standard UPb age ( ±2σ) Reference FC Duluth complex 1099.0 ± 0.6 Ma Paces and Miller, 1993 F5 Duluth complex 1099.0 ± 0.6 Ma Paces and Miller, 1993 (assumed equal to FC1) IF Fish Canyon Tuff 28.201 ± 0.012 Ma Lanphere et al., 2001; Kuiper et al., 2008 MD Mount Dromedary 99.12 ± 0.14 Ma Renne et al., 1998 T2 Temora2, Middledale 416.78 ± 0.33 Ma Black et al., 2004 gabbroic diorite TR Tardree Rhyolite 61.23 ± 0.11 Ma Dave Chew, personal communication

Pb/U Fractionation Factor

Under the operating conditions of LAICPMS sessions, fractionation factors are occasionally found to vary strongly with scan number, decreasing with increasing scan number (presumably

 2016, GeoSep Services, LLC https://mc06.manuscriptcentral.com/cjes-pubs Page 55 of 57 Canadian Journal of Earth Sciences

4

due to increasing ablation pit depth and the effect this has on fractionation; e.g., Paton et al., 2010). The zircon crystal lattice is widely known to accumulate αradiation damage (e.g., Zhang et al., 2009 and references therein). It is assumed that increased αdamage in a zircon leads to a decrease in the hardness of the zircon; this in turn leads to a faster rate of laser penetration into the zircon during ablation leading to dependence of isotopic fractionation on the degree of zircon lattice radiation damage. Ages calculated for all zircon age standards, when those standards were treated as unknowns, were used to construct a fractionation factor correction curve (exponential form) in terms of accumulated radiation damage. The notion of matrixmatched zircon standard and zircon unknown has been proposed largely on the basis of trace element chemistry (e.g., Black et al., 2004). In this study, time and lattice damage, parameters invisible to instruments used to characterize trace element chemistry, were introduced and applied based on measured U and Th chemistries to effectively matrixmatch standard and unknown zircons.

Common Pb Correction

Common Pb was subtracted out using the Stacey and Kramer (1975) common Pb model for Earth. Ages and common Pb ratio were determined iteratively using a preset, sessionwide minimum common Pb age value (default for each session was the age of the oldest age standard which for both Ap and Zrn was 1099 Ma FC1 and/or FC5z).

Preferred Age Draft

Uranium decay constants and the 238 U/ 235 U isotopic ratio reported in Steiger and Yäger (1977) 207 235 235 238 206 238 207 206 were used in this study. Pb/ Uc ( Uc = 137.88 U), Pb/ U, and Pb/ Pb ages were calculated for each data scan and checked for concordance; concordance here was defined as overlap of all three ages at the 1 σ level (the use of 2 σ level was found to skew the results to include scans with significant common Pb). The backgroundcorrected isotopic sums of each isotope were calculated for all concordant scans. The precision of each isotopic ratio was calculated by using the background and signal errors for both isotopes. The fractionation factor for each data scan, corrected for the effect of accumulated αdamage, was weighted according to the 238 U or 232 Th signal value for that data scan; an overall weighted mean fractionation factor for all concordant data scans was used for final age calculation.

Ages and common Pb ratio were determined iteratively using a preset, sessionwide minimum common Pb age value (default for each session was the age of the oldest age standard which for both Ap and Zrn was 1099 Ma FC1 and/or FC5z).

If the number of concordant data scans for a spot was greater than zero, then for grains that yield ages older than ca. 1000 Ma, the 207 Pb/ 206 Pb age is used. For younger grains, the 206 Pb/ 238 U age is used.

Analyses were filtered following the criteria of Gehrels and Pecha (2014). Analyses with >10% uncertainty (1 σ) in 206 Pb/ 238 U or the 206 Pb/ 207 Pb age are not included in our analysis. Concordance is based on 206 Pb/ 238 U age/ 206 Pb/ 207 Pb age, with 100% = concordant. Analyses with >20% discordance (<80% concordance) or >5% reverse discordance (>105% concordance) are not included.

 2016, GeoSep Services, LLC https://mc06.manuscriptcentral.com/cjes-pubs Canadian Journal of Earth Sciences Page 56 of 57

5

Preferred Age Precision

207 235 235 238 206 238 207 206 Errors for the isotopic ratios Pb/ Uc ( Uc = 137.88 U), Pb/ U, and Pb/ Pb at each scan included errors from the backgroundcorrected signal values for each isotope, the fractionation factor error, and an additional relative error term required to force 95% of the FC 207 235 235 238 206 238 ages to be concordant. Errors for the isotopic ratios Pb/ Uc ( Uc = 137.88 U), Pb/ U, and 207 Pb/ 206 Pb at each scan included errors from the backgroundcorrected signal values for each isotope, the fractionation factor error, and an additional relative error term required to force 95% of the FC ages to be concordant. Asymmetrical negativedirection and positivedirection age errors were calculated by subtracting and adding, respectively, the isotopic ratio errors in the appropriate age equation (Chew and Donelick, 2012).

References Cited

Black, L.P., Kamo, S.L., Allen, C.M., Davis, D.W., Aleinikoff, J.N., Valley, J.W., Mundil, R., Campbell, I.H., Korsch, R.J., Williams, I.S., and Foudoulis, C. 2004. Improved 206 Pb/ 238 U microprobe geochronology by the monitoring of traceelementrelated matrix effect; SHRIMP, IDTIMS, ELAICPMS and oxygen isotope documentation for a series of zircon standards. Chemical Geology, 205, p.15140 Draft

Bradley, D., Haeussler, P., O’Sullivan, P., Friedman, R., Till, A., Bradley, D., and Trop, J., 2009, Detrital zircon geochronology of Cretaceous and Paleogene strata across the southcentral Alaskan convergent margin, in Haeussler, P.J., and Galloway, J.P., Studies by the U.S. Geological Survey in Alaska, 2007: U.S. Geological Survey Professional Paper 1760F, 36 p.

Chang, Z., Vervoort, J.D., McClelland, W.C., and Knaack, C., 2006. UPb dating of zircon by LAICPMS. Geochemistry, Geophysics, Geosystems, American Geophysical Union, v. 7, n. 5, 14 p.

Chew, D.M. and Donelick, R.A., 2012. Combined apatite fission track and UPb dating by LA ICPMS and its application in apatite provenance analysis. Mineralogical Association of Canada Short Course, v. 42, 219247.

Donelick, R. A., O'Sullivan, P. B., and Donelick, M. B., 2010, A DiscordiaBased Method of Zircon UPb Dating from LAICPMS Analysis of Single Spots. Smart Science for Exploration and Mining, v. 1 and 2, p. 276278.

Donelick, R.A, O’Sullivan, P.B., Ketcham, R.A., 2005. Apatite fissiontrack analysis. Reviews in Mineralogy and Geochemistry, Mineralogical Society of America, v. 58, p. 4994.

Gehrels, G.E. and Pecha, M. 2014. Detrital zircon UPb geochronology and Hf isotope geochemistry of Paleozoic and Triassic passive margin strata of western North America: Geosphere, 10: 4965, doi:10.1130/GES00889.

 2016, GeoSep Services, LLC https://mc06.manuscriptcentral.com/cjes-pubs Page 57 of 57 Canadian Journal of Earth Sciences

6

Gehrels, G.E., Valencia, V.A., and Ruiz, J., 2008. Enhanced precision, accuracy, efficiency, and spatial resolution of UPb ages by laser ablationmulticollectorinductively coupled plasmamass spectrometry. Geochemistry Geophysics Geosystems, American Geophysical Union, v. 9, 13 p.

Hults, C.P., Wilson, .H., Donelick , R.A., and O’Sullivan , P.B., 2013. Two flysch belts having distinctly different provenance suggest no stratigraphic link between the Wrangellia composite terrane and the paleoAlaskan margin Lithosphere, December 2013, v. 5, p. 575594, first published on November 15, 2013, doi:10.1130/L310.1.

Kuiper, K.F., Deino, A., Hilgen, P.J., Krijgsman, W., Renne, P.R., and Wijbrans, J.R., 2008. Synchronizing rock clocks of Earth history. Science, 320, p. 500504.

Lanphere, M.A. and Baadsraard, H., 2001. Precise KAr, 40 Ar/ 39 Ar, RbSr and UPb mineral ages from the 27.5 Ma Fish Canyon Tuff reference standard. Chemical Geology, v. 175, p. 653671.

Moore. T.E., O’Sullivan, P.B., Potter, C.J., and Donelick, R.A., 2015, Provenance and detrital zircon geochonologic evolution of lower Brookian foreland basin deposits of the western Brooks Range, Alaska, and implications for early Brookian tectonism: Geosphere, v. 11, p. 93–122, doi:10.1130/GES01043.1.

Paces, J.B. and Miller, J.D., 1993. PreciseDraft UPb ages of Duluth Complex and related mafic intrusions, northeastern Minnesota: Geochronological insights to physical, petrogenic, paleomagnetic, and tectonomagmatic processes associated with the 1.1 Ga Midcontinent Rift System. Journal of Geophysical Research, v. 98, no. B8, p. 1399714013.

Paton, C., Woodhead, J.D., Hellstrom, J.C., Hergt, J.M., Greig, A., and Maas, R., 2010, Improved laser ablation UPb zircon geochronology through robust downhole fractionation correction, Geochem. Geophys. Geosyst., 11, Q0AA06, doi:10.1029/2009GC002618.

Renne, P.R., Swisher, C.C., Deino, A.L., Karner, D.B., Owens, T.L., and DePaolo, D.J., 1998. Intercalibration of standards, absolute ages and uncertainties in 40 Ar/ 39 Ar dating. Chemical Geology, v. 45, p. 117152.

Stacey J. S. and Kramer J. D., 1975. Approximation of terrestrial lead isotope evolution by a twostage model. Earth and Planetary Science Letters, v. 26, pp. 207221

Steiger, R.H. and Jäger, E., 1977. Subcommission on geochronology: Convention on the use of decay constants in geo and cosmochronology. Earth and Planetary Science Letters, v. 36, p. 369371.

Zhang, M., Ewing, R.C., Boatner, L.A., Salje, E.K.H., Weber, W.J., Daniel, P., Zhang, Y., and * Farnan, I., 2009, Pb irradiation of synthetic zircon (ZrSiO 4): Infrared spectroscopic study – Reply. American Mineralogist, v. 94, pp. 856858.

 2016, GeoSep Services, LLC https://mc06.manuscriptcentral.com/cjes-pubs