High-Precision CA-ID-TIMS U-Pb Zircon Geochronology of Felsic High-Precision CA-ID-TIMS U-Pb Zircon Geochronology of Felsic Rocks in the Finlayson Lake VMS District, Yukon: Linking Rocks in the Finlayson Lake VMS District, Yukon: Linking Paleozoic Basin-Scale Accumulation Rates to the Occurrence of Paleozoic Basin-Scale Accumulation Rates to the Occurrence of Subseafloor Replacement-Style Mineralization Subseafloor Replacement-Style Mineralization

Felsic igneous complexes and associated volcano-sedimentary rocks in continental back-arc environments host large-tonnage and/or high-grade volcanogenic massive sulfide (VMS) deposits. The emplacement mechanisms, style, and preservation of these deposits is thought to be partially dependent on depositional rates of the host lithofacies (i.e., discrete volcanic eruptions) relative to the setting of massive sulfide genesis on the seafloor as mounds and/or via subseafloor replacement of existing strata. The localization and occurrence of subsea-floor replacement-style VMS deposits is therefore strongly influenced by the characteristics of the volcano-sedimentary facies in the hosting basin and the rates of their emplacement; the latter are poorly constrained in the literature due to the difficulty of obtaining high-precision dates that make this possible in Phanerozoic and older rocks. New high-resolution U-Pb geochronology and detailed regional stratigraphic investigation indicate that Devonian-Mississippian volcanic rocks and associated VMS mineralization in the Yukon-Tanana terrane in the Finlayson Lake district, Yukon, Canada, were erupted or emplaced during distinct time periods (ca. 363.3, 362.8, and 355.2 Ma) in two discrete submarine basins: the Kudz Ze Kayah formation and the Wolverine Lake group. The VMS deposits in both settings are contained within intrabasinal rocks that accumulated at rapid rates of ~350 to 2,000 m/m.y. over 0.6 to 1.4 m.y. Locally, these rates reach peak rates up to 7,500 m/m.y. in the Wolverine Lake group, which are interpreted to reflect facies deposition by mass transport complexes or turbidity currents. These new dates indicate that rapid accumulation of volcanic rocks in the back-arc basins was critical for localizing subseafloor replacement-style mineralization and the development of the Zn-enriched GP4F, Kudz Ze Kayah, and Wolverine VMS deposits. Rapid depositional processes observed in these deposits and their host basins are interpreted to have an important role in developing highly porous and permeable, water-saturated lithofacies that provide optimal conditions for enhancing zone refining processes and subsequent preservation of massive sulfide mineralization, which are key in the development of high-grade and large-tonnage VMS deposits. It is herein suggested that quantitative basin-scale accumulation rates, as a result of new U-Pb geochronological methods and increased precision combined with detailed stratigraphic and facies analysis, may provide important perspectives on the formation of continental back-arc basins and the localization of VMS deposits in other continental margin environments globally.


Introduction
Silicic magmatism at convergent margins contributes to the growth of continental crust and drives the formation of hydrothermal ore deposits in both the arc (e.g., porphyry Cu deposits; Sillitoe, 2010) and back-arc environments (e.g., volcanogenic massive sulfide [VMS] deposits; Franklin et al., 2005;Piercey, 2011).In the back-arc environment, slab rollback, crustal extension, and upwelling of asthenosphere-derived melts elevate the geothermal gradient, increase heat flow, and enhance fluid circulation in the upper crust (McKenzie and Bickle, 1988;Heuret and Lallemand, 2005).Extensional faults within the back-arc rifts facilitate magma transport in the crust and the formation of subvolcanic intrusions, volcanic flows, or explosive pyroclastic eruptions on the seafloor (Gibson et al., 1999).Rapid accumulation of volcanic-subvolcanic facies (± epiclastic sedimentary rocks) in the subaqueous environment of continental back-arc basins provides insight into the nature of volcanism, sedimentation, and subsidence in ancient seafloor environments.Collectively, volcanic and volcano-sedimentary rock emplacement style can also influence the styles of VMS mineralization.Flow-dominated environments commonly lead to exhalative mound development (i.e., black smoker chimneys), whereas in more clastic-dominated environments, replacement-style deposits are common where sulfides replace porous and, to a lesser extent, coherent volcanic and sedimentary facies (e.g., Franklin et al., 2005;Gibson, 2005).Rapid accumulation of volcaniclastic rocks in VMS settings has been proposed as an important aspect of forming large and/or high-grade VMS deposits in the subseafloor environment (e.g., Doyle and Allen, 2003).However, despite the potential tectonic and economic importance of emplacement rates at which volcanic facies are accumulated, the rates, how they are related to the formation and localization of VMS de-posits, and their role in facilitating subseafloor replacement mineralization remain poorly understood.
Deciphering rates of volcanic and sediment accumulation in the ancient rock record is challenging due to the paucity of robust geochronological data in ancient assemblages.Further, most historical geochronological data, or dates obtained by less precise in situ methods (e.g., laser ablation-inductively coupled plasma-mass spectrometry [LA-ICP-MS] or secondary ion mass spectrometry [SIMS]), are not sufficient for robust calculations of ancient emplacement rates at the resolution required to investigate VMS formation.Thus, determining volcanic emplacement rates requires the application of modern, high-precision geochronological methods, such as chemical abrasion-isotope dilution-thermal ionization mass spectrometry (CA-ID-TIMS) U-Pb zircon geochronology (Mattinson, 2005) and precise controls on stratigraphic and structural relationships for the VMS-hosting rock packages (e.g., Bleeker et al., 1999;Ross et al., 2014;Oliver et al., 2021).Current U-Pb zircon geochronological methods allow for a precision of <0.1% on single zircon crystals with weighted mean dates that can decipher distinct Paleozoic magmatic events less than 100 k.y.apart (Schoene and Baxter, 2017).Utilizing modern CA-ID-TIMS methods, we present the first set of new high-precision U-Pb zircon dates for mid-Paleozoic back-arc related rocks in the Kudz Ze Kayah formation and Wolverine Lake group in the Finlayson Lake VMS district, Yukon, Canada, which together contain >25 Mt of polymetallic base metal sulfide resources (Regan, 2007;Peter et al., 2007;van Olden et al., 2019).The rock packages in this district are interpreted to be deposited in relatively small (<600 km 2 ), restricted continental back-arc basins that formed along the western margin of ancestral North America in the Late Devonian to Early Mississippian (e.g., Bradshaw et al., 2008).We integrate stratigraphic data and, to our knowledge, the highest-resolution U-Pb results ever reported for any VMS district in the ancient geologic record to calculate linear and volumetric, time-resolved accumulation rates for volcanic and sedimentary rock assemblages within the Finlayson Lake VMS district.These calculated rates indicate rapid accumulation of large volumes of subaqueous volcanic and sedimentary facies during the early back-arc rifting stage.The results may indicate a possible link between rapid, basin-scale volcanosedimentary processes and the localization and genesis of Zn-rich, replacement-style VMS deposits in the Finlayson Lake district.

Geologic and Stratigraphic Setting
The Finlayson Lake VMS district, Yukon, Canada, is a faultbounded portion of the Yukon-Tanana and Slide Mountain terranes that formed adjacent to the peri-Laurentian continental margin in the mid-to late-Paleozoic (Fig. 1; Tempelman-Kluit, 1979;Mortensen and Jilson, 1985;Murphy et al., 2006).Rocks that compose the Yukon-Tanana terrane arc and back-arc assemblages are volcanic, plutonic, and sedimen-tary rocks that were deposited above or intruded a pre-Late Devonian basement; these rocks are variably deformed and metamorphosed but locally retain primary geologic and geochemical features (Colpron et al., 2006;Murphy et al., 2006;Piercey et al., 2006;Piercey and Colpron, 2009).The Jules Creek transform fault juxtaposes rocks of the Yukon-Tanana terrane with ophiolitic rocks of the Slide Mountain terrane (Figs. 1, 2; Murphy et al., 2006), which were then together thrust above North American platformal strata along the Inconnu thrust in the Late Jurassic (Murphy et al., 2002).In the Eocene, Yukon-Tanana terrane rocks were displaced ~430 km along the Tintina strike-slip dextral fault system, resulting in the Finlayson Lake district (Fig. 1; Gabrielse et al., 2006).
Yukon-Tanana terrane rocks within the Finlayson Lake VMS district are hosted in three distinct structurally imbricated panels (from deepest to shallowest): the Big Campbell, Money Creek, and Cleaver Lake thrust sheets (Figs. 1, 2; Murphy et al., 2006).We will only discuss the geology of the Big Campbell thrust sheet in the scope of this paper, as the Money Creek and Cleaver Lake thrust sheets contain no known VMS mineralization; we refer the reader to Murphy et al. (2006) for descriptions of these rock packages.The Big Campbell thrust sheet is bounded below by the post-Late Triassic Big Campbell thrust fault and above by the Early Permian Money Creek thrust fault (Figs. 1,2).Rocks in the Big Campbell thrust sheet include the basal Upper Devonian metasedimentary rocks of the North River formation overlain by Upper Devonian mafic and felsic metavolcanic and metasedimentary rocks of the Grass Lakes group, which include the Fire Lake, Kudz Ze Kayah, and Wind Lake formations (Figs. 1,2).Late Devonian gabbroic to dioritic intrusions of the North Lakes intrusion and granitoid intrusions of the Grass Lakes plutonic suite cut the Grass Lakes group.The Grass Lakes group is unconformably overlain by Lower Mississippian metaclastic and mafic to felsic metavolcanic rocks of the Wolverine Lake group .Metasedimentary rocks of the Lower Permian Money Creek formation are interpreted to unconformably overlie the Wolverine Lake group.Bimodal felsic and felsic-siliciclastic Zn-enriched VMS mineralization in the Big Campbell thrust sheet is present throughout the Grass Lakes group stratigraphy (GP4F and Kudz Ze Kayah deposits; Piercey et al., 2001;Sebert et al., 2004;Peter et al., 2007) and in the Wolverine Lake group (Wolverine deposit; Bradshaw et al., 2008;Figs. 1-4); mafic-siliciclastic VMS mineralization is present in the Fire Lake formation (Kona Cu-Co deposit; Sebert et al., 2004;Peter et al., 2007).Rocks and associated VMS mineralization in the Grass Lakes and Wolverine Lake groups are interpreted to have formed in an evolving continental arc to back-arc basin environment that represents the earliest stages of rifting on the western peri-Laurentian continental margin (Piercey et al., 2001(Piercey et al., , 2002(Piercey et al., , 2004(Piercey et al., , 2006;;Murphy et al., 2006).
Our discussion below focuses on the felsic VMS-bearing assemblages in the Big Campbell thrust sheet that have been Fig. 1.Regional geologic setting of the Finlayson Lake district, Yukon-Tanana terrane, Yukon (modified after Murphy et al., 2006 andColpron andNelson, 2011).Numbers indicate locations of VMS deposits in the region.Abbreviations: BCT = Big Campbell thrust, CLT = Cleaver Lake thrust, JCF = Jules Creek fault, MCT = Money Creek thrust, NRF = North River fault, VMS = volcanogenic massive sulfide.Inset map of Yukon shows the distribution of Paleozoic terranes (YTT = Yukon-Tanana terrane; SMT = Slide Mountain terrane) relative to the North American continental margin.variably metamorphosed to greenschist facies (Murphy et al., 2006).However, primary textures are preserved throughout the rock units in the Grass Lakes and Wolverine Lake groups (e.g., lapilli, laminations), and as such, we omit the prefix "meta" from all rock descriptions below.Descriptive terminology is used for coherent and volcaniclastic lithofacies (i.e., tuff, lapilli tuff, aphyric) that follows the non-genetic, size-based nomenclature of Fisher (1966), recently updated by White and Houghton (2006).

Map Elements
Fault, normal

Unconformity
Stream/river Lake/pond  (Piercey et al., 2001), were deposited in two distinct continental margin basins, and both contain multiple occurrences of sulfide mineralization (e.g., Figs. 1-3;Peter et al., 2007).Geologic mapping, cross sections, and stratigraphic reconstructions indicate that both stratigraphic units are characterized by lateral variations in thickness over minimum strike lengths of ~15 km (Kudz Ze Kayah formation) and ~23 km (Wolverine Lake group; Figs.3-8).The lower part of the Kudz Ze Kayah stratigraphy consists of siltstone and carbonaceous argillite (~130-500 m thick), but volcaniclastic rocks (i.e., tuffs, lapilli tuffs, and crystal tuffs), aphanitic rhyolitic intrusive rocks, and minor mudstone and mafic intrusive rocks compose the remaining ~400 to 1,300 m of stratigraphy throughout the district (Figs. 4,6,8).The volcaniclastic rocks (RHYv; Fig. 9A-D) are volumetrically dominant and contain facies variations that are indistinguishable at the regional scale; these variations occur both along strike and up-section but are mapped as one regional  .Geologic map of the Wolverine Lake group, merged with regional mapping (after Yukon Geological Survey, 2018).Symbols and labels as in Figures 3 and 4  geologic unit (Fig. 4).Coherent rhyolites (RHYc; Fig. 9E, F) are present as discrete units that crop out over a 6-km strike length and have been transposed subparallel to the regional S1 cleavage (Fig. 4).These rocks are spatially associated with aphyric rhyolites that contain minor quartz eyes (RHYi), and both RHYc and RHYi units exhibit features that are indicative of being stratigraphically concordant; however, surface expressions of the RHYi unit are limited to between the B-B' and C-C' cross sections (<2 km along strike; e.g., Figs. 4,6).Tuffaceous facies are commonly poorly sorted and, where lapilli and/or crystals are present, are matrix-to clast-supported (e.g., Fig. 9B, C) with random grain sizes and modal distributions; rare instances of graded bedding within the tuffs are present.In the Kudz Ze Kayah formation, massive sulfides are present in the felsic-siliciclastic GP4F deposit (1.5 Mt at 6.4% Zn, 3.1% Pb, 0.1% Cu, 2.0 g/t Au, and 81.7 g/t Ag; MacRobbie and Holroyd, 2000), located ~70 m above the siltstonevolcaniclastic facies transition in the lower Kudz Ze Kayah formation, and in the bimodal-felsic Kudz Ze Kayah deposit (18.1 Mt at 6.6% Zn, 2.0% Pb, 0.9% Cu, 1.4 g/t Au, and 156 g/t Ag; van Olden et al., 2019)-consisting of the ABM and Krakatoa zones-located higher in the Kudz Ze Kayah stratigraphy, ~500 to 600 m above GP4F (Fig. 4; Peter et al., 2007).

Structure
Rocks in the Grass Lakes and Wolverine Lake groups have been metamorphosed to greenschist facies and exhibit four phases of deformation (e.g., Fig. 10; Murphy and Piercey, 1999;Murphy et al., 2006).In both areas, local occurrences of fine laminations in argillite define primary bedding planes Fig. 7. Cross sections for lines F-F', G-G', and H-H' in the Wolverine Lake group, drafted with no vertical exaggeration.All aspects of cross sections the same as in Fig. 6.Abbreviations: CRF = Campbell Range formation, MCf = Money Creek formation, WLg = Wolverine Lake group.(S0; e.g., Fig. 10C).The first deformation (D1) is represented by penetrative foliation in the Grass Lakes group (S1-GLg) that is typically subparallel to the S0 bedding planes and is defined by oriented micas, flattened lapilli and feldspar crystals, and fissile argillite (e.g., Figs. 9, 10).The D1 structures are interpreted to result from the regional deformation event at ca. 358 Ma (Murphy and Piercey, 1999;Murphy et al., 2006).A second phase of deformation (D2) is recognized by the following: 1) local crenulation fabrics (S2-GLg) in the Grass Lakes group that are oblique to the S1-GLg in rocks with low competency (e.g., argillite and sericite-altered, well-foliated volcaniclastic rocks; Fig. 10B); 2) rare, mesoscale folds that contain the S0 and S1-GLg foliation (F1-GLg); 3) penetrative foliation in the Wolverine Lake group (S1-WLg) that is subparallel to bedding in argillite (S0), represented by similar oriented mica, feldspar, and lapilli as in the Grass Lakes group; and 4) a weakly foliated, discordant granitic intrusion that cuts previously deformed (S1-GLg) argillites in the Wind Lake formation (Murphy et al., 2006).The D1 and D2 deformations are further supported by distinct S1-GLg (~20°-30° dip to the northnortheast) and S1-WLg (~20°-30° dip to the east) that define the dominant penetrative foliation in the Grass Lakes group and Wolverine Lake group, respectively.The third deformation phase (D3) is present primarily as localized, southwestvergent isoclinal folds (F1-WLg) in the upper Wolverine Lake group footwall and hanging wall (mm-to m-scale) that fold the earlier S 1-WLg fabric.The orientation of these F1-WLg folds together with E-dipping S1-WLg cleavage supports earlier hypotheses that the Wolverine Lake group stratigraphy makes up part of the eastern limb of an open, upright, N-NE-trending, regional anticline (Murphy and Piercey, 1999;Bradshaw et al., 2008).Stratigraphically above the Wolverine Lake group in the Campbell Range formation basalts, tight isoclinal anticline-syncline pairs have been interpreted as the areas of highest strain during the D3 deformation which are considered to be near a synclinal hinge of the regional anticline (Murphy and Piercey, 1999).Finally, the latest deformation (D4) is represented by late, brittle faults in two orientations: 1) NE-E-trending structures that cut both the Grass Lakes and Wolverine Lake groups with ~320 m of mappable offset (e.g., East and Fault Creek faults; Figs. 4, 5) and 2) N-NWtrending structures in the Grass Lakes group (Fig. 3; Murphy and Piercey, 1999).

Unit thickness corrections and assumptions
The unit thicknesses for lithofacies in the Kudz Ze Kayah formation and Wolverine Lake group can be estimated despite deformation in the Big Campbell thrust sheet.We base this interpretation on the following arguments.First, argillites in both the Grass Lakes and Wolverine Lake groups contain local laminae that signify S0 bedding planes (Fig. 10C).Second, primary features are recognized in volcaniclastic rocks despite vertical compressional strain (i.e., flattening), such as lapilli and crystal fragments that are variably oriented with the S1 fabric (Figs.9-11).Bradshaw et al. (2008) note that stratigraphy in the vicinity of the Wolverine deposit is deformed and significantly flattened, but primary stratigraphic and deposit morphologies can still be deciphered.This indicates that although the rocks are metamorphosed to phyllite and schist, conditions were not conducive to full recrystallization and fold-induced thickening of the strata in the study areas.Third, although folding of early penetrative foliation is relatively common throughout the Wolverine Lake group, we interpret the footwall stratigraphy to be situated on the eastern limb of an open N-NE-trending anticline that was not affected by intense strain of a hinge region (e.g., Murphy and Piercey, 1999).Strain intensity increases above the Wolverine Lake group hanging wall in the overlying Money Creek (YTT) and Campbell Range formations (SMT); however, folding is interpreted to not impact the stratigraphic thickness of the Wolverine Lake group footwall stratigraphy (Figs. 5, 7;Murphy and Piercey, 1999;Murphy et al., 2006).Furthermore, we are confident that thickness estimates are reasonable despite widespread strain.
In light of the above points, cross sections were constructed in orientations that are roughly orthogonal to the strike of the penetrative foliation in each area (S1-GLg and S1-WLg; Figs.4-7).For the Kudz Ze Kayah formation, Late Devonian porphyritic granites of the Grass Lakes plutonic suite intrude the entire Grass Lakes group and are thus omitted from thickness reconstructions (Figs. 1,3,4).On the surface, these porphyritic rocks are ~40-to 60-m-thick sills in the middle part of the Kudz Ze Kayah formation, whereas units in the western Kudz Ze Kayah formation are thicker (~200-250 m; Fig. 4).In the Wolverine Lake group, porphyritic rocks are interpreted as synvolcanic intrusive rocks (e.g., Piercey et al., 2008), Creek faults (D4) and are corrected for a maximum of ~320 m of normal displacement; the D4 faults in the Wolverine Lake group, however, are roughly parallel to and therefore do not intersect the cross sections.To avoid assumptions related to the intensity of strain throughout the study areas, the degree of compression is left uncorrected, and the final thicknesses and associated rates presented are thus minimum estimates.Finally, due to the lack of high-precision geochronological constraints above the Kudz Ze Kayah formation and Wolverine Lake group, the Wind Lake formation and Wolverine Lake group hanging-wall strata, respectively, are not included in rate calculations below.

Analytical Methods
Sample preparation and mineral separation procedures were carried out at Memorial University of Newfoundland (MUN).Zircon was extracted from rock samples following a series of crushing, grinding, splitting, sieving (to <500 μm to >63 μm), heavy liquid (bromoform, methylene iodide), and magnetic separation methods.The zircon grains were then hand-picked in ethanol, loaded into quartz crucibles, and placed into a furnace for annealing at 900°C for 60 hours to treat minor radiation damage; annealing enhances cathodoluminescence (CL) emission (Nasdala et al., 2002), promotes more reproducible interelement fractionation during laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS; Allen and Campbell, 2012), and prepares the crystals for subsequent chemical abrasion (Mattinson, 2005).Annealed grains were mounted in epoxy and ground and polished with 12-μm Al oxide grit, 3-μm diamond paper, and 1-μm Al oxide powder in steps until the grains were sectioned approximately halfway through.Samples were then carbon coated and imaged with cathodoluminscence (CL), backscattered electron (BSE), and secondary electron (SE) techniques at voltages of 10kV (CL) and 15kV (BSE, SE) on a JEOL JSM 7100F scanning electron microscope equipped with a Schottky-type Field Emission gun and a Deben Centaurus CL detector.Zircon CL images are presented in Figure 12 and Appendix Figures A3 and A4.LA-ICP-MS and CA-ID-TIMS analyses were performed at the Isotope Geology Laboratory at Boise State University, Idaho (BSU).The LA-ICP-MS analyses were completed prior to CA-ID-TIMS to target and omit zircon grains with components of Pb loss or inheritance that would otherwise  A-D) and the Wolverine Lake group (E-I).A) Felsic crystal tuff with blue quartz eyes located ~5 m below the Wind Lake-Kudz Ze Kayah formation contact (17MM-002).B) Felsic tuff to lapilli tuff in the hanging wall of the ABM zone (17MM-031).C) Felsic tuff with minor lapilli in the footwall of the Krakatoa horizon (18MM-133).D) Felsic crystal tuff with feldspar and blue quartz in the footwall of the GP4F deposit (17MM-074).E) Feldspar porphyritic intrusive rock in the Fisher zone (18MM-102).F) Felsic crystal tuff with feldspar and quartz fragments from the footwall of the Wolverine/Lynx zone (P99-WV-4K).G) Feldspar porphyritic intrusive rock from the footwall of the Wolverine/Lynx zone (P00-WV-1C).H) Quartz-feldspar porphyritic intrusive rock from the Sable zone (P00-WV-12).I) Felsic crystal tuff with quartz and feldspar in the lower Wolverine Lake group footwall (18MM-114).Mineral abbreviations: fsp = feldspar, py = pyrite, qtz = quartz.
be unsuitable for obtaining high-precision CA-ID-TIMS analyses.Laser ablation was performed on each grain using a New Wave Research UP-213 Nd:YAG UV laser (213 nm) using fluence and pulse rates of ~5 J/cm 2 and 10 Hz, respectively, during a 45-second analysis (15-sec gas blank, 30-sec ablation).Spot sizes are 25 μm in diameter and, where necessary, raster lines of 10-μm diameter and ~50-mm length were used to avoid burning through grains.Ablated material was carried by a 1.2-L/minHe gas stream to the nebulizer flow of the plasma and analyzed in a ThermoElectron X-Series II quadrupole ICP-MS.Quadrupole dwell times were 5 ms for Si and Zr, 200 ms for 49 Ti and 207 Pb, 80 ms for 206 Pb, 40 ms for 202 Hg, 204 Pb, 208 Pb, 232 Th, and 238 U, and 10 ms for all other high field strength elements (HFSEs) and rare earth elements (REEs); total sweep duration is 950 ms.Background count rates for each analyte were obtained prior to each spot analysis and subtracted from the raw count rate for each analyte.For concentration calculations, background-subtracted count rates for each analyte were internally normalized to 29 Si and calibrated with respect to NIST SRM-610 and -612 glasses as the primary standards.Standards (Plešovice, Seiland, and Zirconia) were used throughout the run, four for every ten unknowns, to monitor for instrument drift.Data processing was completed with a BSU in-house Excel-based program.
CA-ID-TIMS U-Pb geochronology methods for isotope dilution thermal ionization mass spectrometry follow those previously published by Davydov et al. (2010) and Schmitz and Davydov (2012).Zircon crystals were plucked from the epoxy mount following LA-ICP-MS analysis and subjected to a modified version of the chemical abrasion method of Mattinson ( 2005), whereby single crystal fragments plucked from grain mounts were individually abraded in a single step with concentrated HF at 190°C for 12 h.The CA treatment effectively eliminates zircon domains that have experienced postcrystallization Pb loss, thereby producing more concordant and accurate U-Pb dates.Zircon fragments were dissolved in Parr bombs at 220°C for 48 h.Dissolved zircon solutions were subsequently dried down, redissolved in 100 μl of 6N HCl and converted to chlorides in Parr bombs at 180°C for 12 h, after which solutions were dried again and brought up in 50 μl of 3N HCl.Uranium and Pb were isolated by anion exchange column chromatography using 50-μl columns and AG-1 X8 resin (200-400 mesh, chloride form [Eichrom]; Krogh, 1973).
The U-Pb aliquot was loaded in a silica gel emitter (Gerstenberger and Haase, 1997) to an outgassed, zone-refined Re filament.Isotopic determinations were performed using an IsotopX PhoeniX-62 TIMS.A correction for mass-dependent Pb fractionation was applied based on repeated measurements of NBS 982 (Catanzaro et al., 1968) Pb (on both the Daly ion counter [0.16 ± 0.03 %] amu -1 ; 1s) and the Faraday cups (0.10 × [1 ± 0.02%] amu -1 ; 1s).Uranium was run as an oxide (UO2) and measured in static mode on Faraday detectors equipped with 10 12 Ω resistors.The U mass fractionation for the same analyses was calculated using the 233 U/ 235 U ratio of the EARTHTIME ET535 mixed 233-235 U-205 Pb double spike solution (0.99506 ± 0.01%, 1s).The EARTHTIME ET2535 mixed 233-235 U- 202-205 Pb double spike solution (0.99923 ± 0.027%, 1s) was also utilized for certain samples (denoted by * beside the sample number in Table 1) and, where applicable, the mass-dependent Pb fractionation was corrected for during the run.
Raw U and Pb data were filtered using the Tripoli software program (Bowring et al., 2011), and the U-Pb dates and uncertainties for each analysis were calculated using the algorithms of Schmitz and Schoene (2007), the U decay constants of Jaffey et al. (1971), and a value of 238 U/ 235 U = 137.88.Uranium oxide measurements were corrected for isobaric interferences using an 18 O/ 16 O value of 0.00206.Uncertainties are based upon non-systematic analytical errors, including counting statistics, instrumental fractionation, tracer subtraction, and blank subtraction.All non-radiogenic Pb was attributed to laboratory blank with a mean isotopic composition determined by total procedural blank measurements.These error estimates should be considered when comparing the 206 Pb/ 238 U dates Fig. 12. Representative cathodoluminescence images for zircon selected for chemical abrasion-isotope dilution-thermal ionization mass spectrometer (CA-ID-TIMS) analysis in A) the Kudz Ze Kayah formation and B) the Wolverine Lake group footwall.Each grain corresponds to an individual CA-ID-TIMS analysis that corresponds with Table 1.The 50-micron scale bar is applicable to all grains in the figure.with those from other laboratories that used tracer solutions calibrated against the EARTHTIME gravimetric standards.When comparing our dates with those derived from other isotopic decay schemes (e.g., 40 Ar/ 39 Ar, 187 Re-187 Os), the uncertainties in tracer calibration (0.05%; Condon et al., 2015;McLean et al., 2015) and U decay constants (0.108%; Jaffey et al., 1971) should be added to the internal error in quadrature.
Quoted errors for calculated weighted means are thus of the form ±X/Y/Z (±2s), where X is solely analytical uncertainty, Y is the combined analytical and tracer uncertainty, and Z is the combined analytical, tracer and 238 U decay constant uncertainty (Schoene et al., 2006).Concordia diagrams, regression intercepts, and weighted averages were produced with Isoplot 4.1 (Ludwig, 2003) A1 and A2, and Figures A5 and A6.

CA-ID-TIMS and LA-ICP-MS U-Pb Zircon Geochronology
Integrated field mapping and drill core logging were utilized for precise stratigraphic sampling in the Upper Devonian Kudz Ze Kayah formation and the Early Mississippian Wolverine Lake group (Figs.3-5).Four CA-ID-TIMS samples were collected in the Kudz Ze Kayah formation, and five samples were obtained from the Wolverine Lake group for the purpose of precisely defining the ages of important VMS deposits in the Finlayson Lake district (i.e., Kudz Ze Kayah, GP4F, and Wolverine) and constraining ages of host stratigraphy in each rock package (Fig. 11).Detailed rock descriptions are presented in the Appendix.

Reporting of CA-ID-TIMS U-Pb zircon dates
The dispersion of individual CA-ID-TIMS 206 Pb/ 238 U dates indicates that protracted crystallization of zircon occurred during the formation of rocks in the Finlayson Lake district (Table 1; Figs. 12, 13).Protracted crystallization suggests that the oldest zircon date represents the time of zircon saturation, whereas the youngest date represents the final stage of zircon crystallization at the solidus temperature (e.g., Wotzlaw et al., 2013;Rioux et al., 2012;Samperton et al., 2015); the established array is thus defined by a total duration (Dt) of zircon within the host magma that includes the 2s uncertainty (95% confidence interval).However, for each individual sample presented for the Kudz Ze Kayah formation and Wolverine Lake group, the individual zircon dates record uncertainties that are in many cases larger than the dispersion between each zircon date (e.g., Fig. 13).Therefore, we report a single weighted mean 206 Pb/ 238 U date for each sample population with a 2s uncertainty that directly overlaps within error of the youngest individual zircon date.
Krakatoa footwall (18MM-133): Zircon grains from sample 18MM-133 are euhedral and ~80 to 120 mm in length with aspect ratios near ~2:1.The grains have oscillatory zoning with alternating CL-bright and CL-dark zones, and minor sector zoning (Fig. 12A, App.Fig. A3).The CA-ID-TIMS results show seven concordant zircon fractions with high Th/U (0.99-1.2, n = 6), where one fraction yielded a low Th/U of 0.20 (z5).These seven grains gave a weighted mean 206 Pb/ 238 U date of 362.847 ± 0.099 Ma (2s; MSWD = 0.47; Table 1; Fig. 13C), which is equivalent within uncertainty to the ABM hanging wall (17MM-031); this date is interpreted as the eruption age of this rock and approximate age of the Krakatoa deposit.One zircon fraction (z10) overlaps with the weighted mean U-Pb date but is not included in the weighted mean calculation because of its high uncertainty that increases the MSWD.The LA-ICP-MS results are concordant with variable Th/U (0.42-3.8) and a weighted mean 206 Pb/ 238 U date of 360.0 ± 1.6 Ma (2s; MSWD = 0.70; n = 34), which is younger outside of 2s uncertainty than the CA-ID-TIMS date.One grain gave an older 207 Pb/ 206 Pb date of ca.1130 Ma and is interpreted as a xenocryst.
Wolverine/Lynx feldspar porphyry (P00-WV-1C): A feldspar porphyritic rock (P00-WV-1C) in the Wolverine/Lynx deposit footwall yielded zircon that is euhedral and relatively small (~60-90 mm in length) with low aspect ratios (1:1 to 2:1).The zircon grains display oscillatory zoning and sector zoning shown by CL-bright zones (e.g., z2 and z3; Fig. 12B, App.Fig. A4).Results from CA-ID-TIMS analyses reveal zircon grains that are concordant with high Th/U (0.90-1.1) and a weighted mean 206 Pb/ 238 U date of 355.262 ± 0.066 Ma (2s; MSWD = 0.99; Table 1; Fig. 13G); this date is inferred as the crystallization age of the rock and the age of mineralization in the Wolverine/Lynx deposit.One zircon fraction (z6) is inferred to be a xenocryst with a concordant U-Pb date of ~363.3Ma.The LA-ICP-MS analytical results are concordant and yield variable Th/U (0.14-3.7) with a weighted mean 206 Pb/ 238 U date Fig. 14.Summary of chemical abrasion-thermal ionization mass spectrometer (CA-TIMS) U-Pb geochronological results in ranked order.A) Ranked U-Pb zircon CA-isotope dilution (ID)-TIMS geochronology results for the Kudz Ze Kayah formation (blue) and the Wolverine Lake group (green).Each vertical bar represents a single 206 Pb/ 238 U date with maximum 2s extents (Table 1 of 352.8 ± 2.2 Ma (2s; MSWD = 1.07; n = 25); this date is slightly younger outside of uncertainty than the CA-ID-TIMS date.Three older grains gave dates of ca.370 to 371 Ma (n = 2) and a 207 Pb/ 206 Pb date of 1105 Ma (n = 1), which are all interpreted to be xenocrystic zircon.
Sable quartz-feldspar porphyry (P00-WV-12): Zircon grains in a quartz-feldspar porphyritic rock (P00-WV-12) in the Sable zone are euhedral and prismatic with grain sizes between 100 and 250 mm in length.The grains have the highest aspect ratios of samples in the Wolverine Lake group, ranging between ~2:1 and 4:1, and show oscillatory zoning with CL-bright cores and CL-dark rims (Fig. 12B, App.Fig. A4).The CA-ID-TIMS data show five concordant zircon fractions with relatively low Th/U (0.45-0.70) and a weighted mean 206 Pb/ 238 U date of 356.094 ± 0.059 Ma (2s; MSWD = 0.67; Table 1; Fig. 13H), interpreted as the crystallization age of the rock.Two zircon fractions (z5 and z7) have high Th/U (1.15-1.25)and concordant U-Pb dates of ~363.4 and ~363.5 Ma, which are interpreted as xenocrysts.The LA-ICP-MS data give concordant results with variable Th/U (0.15-2.2) and a weighted mean 206 Pb/ 238 U date of 354.6 ± 2.1 Ma (2s; MSWD = 1.05; n = 22), which overlaps within error of the CA-ID-TIMS date.Five zircon grains have older dates between ca.366 and 369 Ma and are interpreted as xenocrysts.

Volcano-Sedimentary Accumulation Rates and Volume Estimates
Rates of volcano-sedimentary accumulation in submarine settings are dependent on the character of volcanism, tectonic and depositional setting, and proximity to the continental slope and/or arc edifice (e.g., Karig and Moore, 1975;Carey and Sigurdsson, 1984;Díaz-Asencio et al., 2019).Quantitative estimates of volcaniclastic sedimentation rates on the modern seafloor are typically between 300 and 1000 m/m.y. for Eocene continental margin basins and younger early intra-arc rift settings (Karig and Moore, 1975;Clift, 1995;Marsaglia et al., 1995;Carey, 1999).Rates for background sedimentation devoid of volcanism are typically lower and in continental shelf environments, such as in the Gulf of Mexico, where rates decrease from the slope to the deeper abyssal plain areas (~130-180 m/m.y. to ~30-50 m/m.y.; Santschi and Rowe, 2008;Díaz-Asencio et al., 2019).The Middle Valley sedimented rift in the Juan de Fuca Ridge contains background Holocene sedimentation rates of 55 to 65 m/m.y., but much higher rates related to Pleistocene turbidite influx have been also recorded at this site (400 m/m.y.; Goodfellow and Franklin, 1993).Volumetric rates for subaerial, silicic volcanic eruptions in continental and oceanic tectonic settings range from 2 to 11,550 km 3 /m.y.
(median = 921 km 3 /m.y.; n = 79; White et al., 2006, and references therein).Moreover, the volcano-sedimentary accumulation rates in the literature vary significantly between 30 and 1,000 m/m.y. and are dependent on the sediment source and the association with active submarine volcanism (e.g., Karig and Moore, 1975;Díaz-Asencio et al., 2019).Stratigraphic reconstructions and new CA-ID-TIMS U-Pb results provide constraints for estimations of volcano-sedimentary accumulation rates and volume estimates for intrabasinal facies in the Kudz Ze Kayah formation and Wolverine Lake group.Due to the lack of drilling constraints and minimal outcrop exposure in low-lying areas, we have simplified the volcanic facies in the entire Kudz Ze Kayah formation and Wolverine Lake group footwall into discrete units, although both are composed of numerous individual eruptions that were likely deposited at much faster rates.In the case of the Kudz Ze Kayah formation and Wolverine Lake group, the major map units are relatively continuous along strike, and are estimated to have strike lengths of at least 15 and 23 km, respectively (Figs. 4,5).The lateral continuity of these basins is similar in basin structure to various parts of the Okinawa Trough (e.g., Arai et al., 2017;Minami et al., 2020), which has been interpreted as the closest modern analogue to the Finlayson Lake district (Piercey et al., 2001).The stratigraphic locations of U-Pb samples are thus extrapolated to similar locations in all cross sections and rates are calculated based on the thickness of the column divided by the total Dt for the entire unit, unless otherwise noted.For each case, the minimum Dt is the smallest difference between the weighted mean 2s envelopes and the maximum Dt is the total range of single crystal dates, including the 2s error.The Wolverine Lake group has been attributed to forming in a basin of at least ~14 by 10 km (Bradshaw et al., 2008); however, we recognize a much greater strike length of ~23 km (Fig. 3).Volumetric estimates for the regional units are calculated using a modified method of Kokelaar and Busby (1992) and Hudak et al. (2003) to account for a basin shape that is assumed to have a length greater than its width.Thus, we use the formula for an elliptical cylinder (V = pabh), where a is half the strike length, b is 70% of a, and h is the dip-corrected unit thickness (Fig. 8).Moreover, these calculations represent minimum estimates for both rates and volumes due to post-depositional lithification, deformational compression, and erosion (e.g., Kokelaar and Busby, 1992;Hudak et al., 2003).
The U-Pb dates for the Kudz Ze Kayah formation indicate that volcanism lasted <1 m.y. and accumulated with rates between 350 and 1,900 m/m.y.(Table 2).Taking only the facies thickness of ~600 m and Dt ~400 k.y. between the GP4F and Kudz Ze Kayah deposits, higher calculated rates between 1,000 and 2,900 m/m.y. are observed.Repeated occurrences of thick (>10-to 50-m), unsorted to poorly sorted crystal and lapilli-rich tuffs at and in between the two VMS deposits may indicate that these units were deposited as volcaniclastic mass-flows, which have been interpreted to represent rapid emplacement and burial rates during subaqueous volcanic eruptions (e.g., App.Fig. A1; McPhie et al., 1993;White, 2000).Estimated volumes of facies along the ~15-km strike length yield an average of ~150 km3 of basin fill, which correspond to time-resolved rates between 110 and 180 km 3 /m.y.In the Wolverine Lake group, the footwall volcano-sedimentary facies reflect up to 1.5 m.y. of deposition that accounts for average accumulation rates of ~800 to 1,400 m/m.y., similar to those in the Kudz Ze Kayah formation (Table 2).Facies below the stratigraphic level of the Sable zone (i.e., sample P00-WV-12), not including the lowermost quartzofeldspathic meta-sandstone, were deposited at a rapid rate of between 2,300 and 7,800 m/m.y.(Dt = ~150-400 k.y.) and are interpreted to be a result of volcanic and sediment influx by mass transport complexes or turbidity currents similar to those observed in the hanging wall of the Wolverine deposit (Piercey et al., 2016).The longer strike length (~23 km) of the Wolverine Lake group strata alone indicates larger volumes of basin fill than in the Kudz Ze Kayah formation, yielding a minimum volume of ~440 km 3 with rates between 280 and 360 km 3 /m.y.(Table 2).The total volcanic contributions in the Wolverine Lake group are approximately 1.8 times greater than argillite, whereas in the Kudz Ze Kayah formation, volcanic facies are much more abundant (~3.8 times) than siltstones and argillites.

Discussion
The new lithostratigraphic and U-Pb geochronological constraints for the Finlayson Lake VMS district provide insights into the timing and duration of VMS mineralization and rates of basin-scale volcano-sedimentary depositional processes in an ancient continental margin setting.Our results provide the first robust dates, absent of Pb loss effects and inheritance, for all known felsic-hosted, Zn-enriched VMS deposits in the Finlayson Lake district and provide insight into depositional rates at a resolution that has not been previously documented in the ancient geologic record.These results provide improved stratigraphic constraints for the Kudz Ze Kayah formation and Wolverine Lake group and directly impact regional correlations throughout the Yukon-Tanana terrane.The discussion below highlights the importance of integrating detailed stratigraphic and high-precision CA-ID-TIMS geochronological results to establish links between depositional rates and replacement-style mineralization in seafloor and subseafloor environments.

Timing and duration of felsic-hosted VMS mineralization
The age and duration of mineralizing processes in VMS districts is critical to understanding the development of hydrothermal systems in the upper crust and identifying prospective areas for future exploration; however, depositional rates are challenging to determine in the ancient record, and this limits our ability to compare rates of deposit formation between modern seafloor massive sulfide (SMS) and ancient VMS systems.Current U-Pb zircon geochronology methods, coupled with precise stratigraphic controls, can constrain the lower and upper age extents of VMS mineralization within a district (e.g., Bleeker and Parrish, 1996;Bleeker et al., 1999;Ross et al., 2014;Mortensen et al., 2015;Oliver et al., 2021), but do not allow for direct dating of hydrothermal mineralization except in the exceptionally rare cases where hydrothermal zircon is observed (e.g., Iberian Pyrite Belt; Nesbitt et al., 1999).Therefore, quantifying the age and duration of mineralization has historically been challenging due to the uncertainties on U-Pb dates (>0.2% at best on individual zircon crystals before 2010; Schoene, 2014) that often are much greater than the lifespan of hydrothermal systems (~10 k.y. to 1 m.y.; e.g., Cathles et al., 1997;Barrie et al., 1999;Hannington et al., 2005).With modern CA-ID-TIMS U-Pb zircon geochronology, however, the uncertainties are now typically <0.1% (Schoene and Baxter, 2017) and stratigraphically bracketed U-Pb dates in a VMS deposit (i.e., stratigraphically above and below, or below and cross cutting the VMS deposit) are at much higher resolution than dates that have been obtained from direct dating of the mineralization event using U-Pb or other isotopic methods (e.g., Mathur et al., 1999;Nesbitt et al., 1999;Nozaki et al., 2014).Zircon from nine U-Pb samples from the Kudz Ze Kayah formation and Wolverine Lake group were used to define the lower and upper temporal bounds of each rock succession and precise ages of VMS deposits in each unit .Zircon grains dated from four samples of volcaniclastic rocks from the Kudz Ze Kayah formation yield tightly constrained weighted mean 206 Pb/ 238 U dates (ca. 363.3-362.4Ma) that distinguish the ages of VMS mineralization in the GP4F, ABM, and Krakatoa zones, and the absolute duration of volcanic activity (Dt = 0.65-1.0m.y.; Fig. 14). At lowe stratigraphic levels, the earliest-known volcanic phase in the Kudz Ze Kayah formation is 363.254 ± 0.098 Ma (17MM-074), interpreted to be the age of the GP4F deposit.Approximately 600 m stratigraphically above GP4F, both the Krakatoa and ABM zones of the Kudz Ze Kayah deposit yield nearly identical dates of 362.847 ± 0.099 Ma (18MM-133) and 362.82 ± 0.12 Ma (17MM-031), respectively.The last phase of felsic volcanism in the Kudz Ze Kayah formation occurred at 362.404 ± 0.098 Ma (17MM-002), which was determined from zircon in a felsic crystal tuff that is ~250 m stratigraphically above the Kudz Ze Kayah deposit.Similarly, the new U-Pb zircon results for five rocks in the Wolverine Lake group (ca.357.5-354.9Ma) constrain the timing of mineralization for the Wolverine/Lynx, Sable, and Fisher zones, and define the duration of volcanosedimentary activity for the footwall facies (Dt = 1.2 to 1.5 m.y.; Fig. 13).Following a period of sedimentation at ca. 357.5 Ma (Fig. 15; Murphy et al., 2006), felsic volcanism in the lower footwall unit began at 356.371 ± 0.091 .Approximately 1,000 m of crystal tuffs and argillites separate the lower footwall from the Sable zone, which was cut by a quartz-feldspar porphyritic intrusion at 356.094 ± 0.059 Ma (P00-WV-12).The Wolverine/Lynx zone is ~2 km northwest along strike from the Sable zone, where overlapping dates for a crystal tuff (355.236± 0.058 Ma; P99-WV-4K) and a feldspar-porphyritic rhyolite (355.262± 0.066 Ma; P00-WV-1C) define the age of the Wolverine deposit footwall and the maximum age of the VMS mineralization (Piercey et al., 2008).Stratigraphy in the Fisher zone, approximately 9 km along strike to the northwest of the Wolverine/Lynx zone, was intruded by a large feldspar-porphyritic stock at 354.982 ± 0.068 Ma and defines the latest phase of volcanism in the Wolverine Lake group footwall unit, interpreted to be the minimum estimate for VMS mineralization.
The strong lateral continuity of lithofacies in the Kudz Ze Kayah formation and Wolverine Lake group, combined with new high-precision U-Pb dates for VMS-bearing stratigraphy, is interpreted to suggest that volcanism and associated VMS mineralization occurred at distinct periods during the life of the back-arc basin (Table 1; Fig. 14).For example, CA-ID-TIMS dates for the GP4F (sample 17MM-074) and Kudz Ze Kayah deposits (17MM-031 and 18MM-133) show that the maximum extent of 2s errors for weighted mean 206 Pb/ 238 U dates are ~200 and ~250 k.y., respectively, indicating that deposition of mineralization occurred over 200 to 250 k.y.(Fig. 14).In the Wolverine Lake group, porphyritic rocks have been shown to be either pre-mineralization (quartz-feldspar porphyritic rocks) or syn-to post-mineralization (feldspar porphyritic rocks; Piercey et al., 2008), thus providing maximum and minimum ages for mineralization.The quartz-feldspar porphyritic rocks intruded at ca. 356.1 Ma, whereas the feldspar-porphyritic rocks intruded between ca.355.2 and 355.0 Ma; each sample accounts for a maximum of ~120 to 140 k.y. of 2s uncertainty on the weighted mean date.These results further indicate that deposition of VMS mineralization for the Wolverine/Lynx (P00-WV-1C and P99-WV-4K) and Fisher zones (18MM-102) occurred over ~150 and 140 k.y., respectively (Fig. 14).
Accumulation rates for massive sulfide horizons in ancient VMS deposits in the literature are exceptionally rare (e.g., Bleeker and Parrish, 1996) but critical to facilitate comparisons to modern SMS deposits (e.g., Jamieson et al., 2014) and to provide a broader understanding of the rates of SMS/VMS formation in modern and ancient settings.The calculated rates were determined as follows: 1) the maximum and minimum age ranges for a given deposit based on stratigraphic and cross cutting relationships were determined; and 2) the massive sulfide tonnage was then divided by these age ranges to determine the minimum and maximum accumulation rates.Estimates for the Kudz Ze Kayah deposit (~19 Mt; van Olden et al., 2019) suggest mass accumulation rates of ~75 to 800 t/ yr, and much lower rates of ~8 t/yr for GP4F (~1.5 Mt; Peter et al., 2007).Estimated depositional rates for the Wolverine deposit (~5.2 Mt; Regan, 2007) are between ~40 and 200 t/yr.These rates are comparable to estimates for the Kidd Creek deposit, where average rates were suggested to be between 10 and 100 t/yr and peak rates reached as high as 1,000 t/ yr (Bleeker and Parrish, 1996), and for modern hydrothermal deposits where mass accumulation rates lie between 1 and 794 t/yr (e.g., Jamieson et al., 2014).Our data further indicate the following: 1) it is possible to compare ancient VMS accumulation rates with modern SMS accumulation rates using highresolution U-Pb geochronology; 2) rates are similar in ancient deposits when compared to modern deposits; however, the ancient deposits are limited to the Finlayson Lake district (this study) and Kidd Creek (Bleeker and Parrish, 1996); and 3) rate resolution between modern SMS and ancient VMS can only be determined using high-precision geochronology methods, like the CA-ID-TIMS data presented herein.

Relationship between Kudz Ze Kayah and Wolverine Lake basins
The Kudz Ze Kayah formation and Wolverine Lake group contain volcano-sedimentary facies that are similar in texture and geochemical affinity; however, the units represent two distinct periods of magmatism and deposition along the western Laurentian continental margin (Piercey et al., 2001;Murphy et al., 2006).Our U-Pb results show that the Kudz Ze Kayah formation and Wolverine Lake group formed between ca.363.3 and 362.4 and ca. 357.5 and 354.9 Ma,respectively (Figs. 13,14).Further, Proterozoic zircon xenocrysts, as shown by LA-ICP-MS results, are present in both Kudz Ze Kayah formation and Wolverine Lake group footwall samples and are interpreted to reflect melting of a crustal substrate with Laurentian affinity (e.g., Snowcap assemblage; Piercey and Colpron, 2009;Fig. 15; App.Fig. A6, Table A2); these Fig. 15.Schematic diagrams illustrating the evolution of Upper Devonian and Early Mississippian submarine back-arc basins in the Finlayson Lake volcanogenic massive sulfide (VMS) district, Yukon-Tanana terrane.A) Early-stage rifting resulting in the deposition of the Kudz Ze Kayah formation (KZKf), where felsic volcaniclastic rocks were deposited in association with hydrothermal upwelling to form the GP4F deposit.B) Upwelling basalt and accelerated extension resulting in tectonic subsidence provide favorable conditions for rapid volcaniclastic accumulation, hydrothermal activity, and subseafloor replacement-style mineralization in the Kudz Ze Kayah deposit (ABM and KKT zones) at ca. 362.8 Ma.C) Waning of the rifting leads to eruption and deposition of alkalic basalts and mudstones in the Wind Lake formation (WLf) above the Kudz Ze Kayah formation, which are then cut by plutons and dikes of the Grass Lakes plutonic suite (GLPS).D) A period of D1 deformation causes penetrative S1-GLg fabrics in the Grass Lakes group (GLg).The rocks are uplifted, tilted, and eroded into a new, incipient back-arc basin (Wolverine Lake group; WLg), transporting zircon with ages ca.363.5 to 358 Ma into the new basin contained in unit Wcl (Murphy and Piercey, 1999;Murphy et al., 2006).E) Rifting resumes and volcaniclastic rocks (unit XLT) are deposited in thick accumulations, partly on top of the older rift stratigraphy.Thick argillite packages are deposited from mass transport complexes and/or turbidity currents from the continental slope as a result of rapid extension and steep basin topography.These strata are intruded by quartz-feldspar porphyritic sills and dikes at ca. 356.1 Ma.F) Felsic volcaniclastic rocks (unit FLT) are deposited and intercalated with argillite.Rapid deposition of footwall tuffs and feldspar porphyritic intrusive rocks between ca.355.2 and 355.0 Ma, coupled with hydrothermal upwelling, facilitated subseafloor and intermittent seafloor mineralization in the Wolverine deposit and Fisher zone.Yellow symbols and numbers correspond to U-Pb samples and stratigraphic locations as in Figures 4 through 8 and 14.For simplicity in panels (E) and (F), units in the Grass Lakes group are shaded to a single gray color.results also concur with previous geochronology work in the Finlayson Lake district (Grant, 1997;Murphy et al., 2006;Piercey et al., 2008).From CA-ID-TIMS results, xenocrystic zircon in three of five rocks from the Wolverine Lake group indicate that Early Mississippian magmatic rocks in the group record assimilation or mixing with sources that contain zircon with dates coeval with the stratigraphically lowest parts of the Kudz Ze Kayah formation (ca. 363.5-363.3Ma; n = 3; Table 1; Figs. 13,14).Similarly, one xenocrystic zircon (ca.361.9 Ma) was dated in sample 18MM-114 from the lower Wolverine Lake group footwall (unit XLT; Fig. 5) and has an age similar to those for the Grass Lakes plutonic suite that intrudes the Grass Lakes group (Murphy et al., 2006).These observations are important for three reasons.First, the stratigraphic position of the crystal tuff unit (XLT) in the lower Wolverine Lake group, which contains the ca.361.9 Ma xenocryst, suggests that erosion or mixing of the underlying Grass Lakes group units was more widespread and was not just recorded in the basal sedimentary rocks (i.e., Wcl unit) of the group (e.g., Murphy and Piercey, 1999); however, it is also possible that this grain was inherited by the magmas that were responsible for generating the crystal tuffs (Fig. 3).Second, older ca.363.5-363.3Ma xenocrysts in feldspar porphyritic, quartz-feldspar porphyritic, and tuffaceous units at the highest levels of Wolverine Lake group stratigraphy indicate that zircon was either entrained by magmas through crustal assimilation or, in the case of the felsic crystal tuff (P99-WV-4K), mixed by sedimentary processes either by passing through zones of the lower Kudz Ze Kayah formation or zones of the basal sedimentary rock unit (Wcl) that contained older zircon grains from the Kudz Ze Kayah formation (Fig. 15).Third, the oldest xenocrysts from CA-ID-TIMS results are among the oldest Late Devonian dates presented in this study and could indicate that there are older Late Devonian rocks (>363.5 Ma) beneath the GP4F deposit (ca.363.3 Ma) that have yet to be recognized.These results illustrate that there may be more important basement connections and crustal inheritance between the two basins at depth that, until now, have not been recognized.

Basin-scale accumulation rates and comparisons to modern environments
Our U-Pb results are the first to constrain the rates of magmatism and volcano-sedimentary deposition in an ancient Paleozoic seafloor VMS environment, with accumulation rates of ~350 to >2,000 m/m.y. in the Finlayson Lake district.These estimates are greater than those observed on the modern seafloor, such as in continental margin basins and early intra-arc rift settings (300-1,000 m/m.y.; Karig and Moore, 1975;Clift, 1995;Marsaglia et al., 1995;Carey, 1999).Volcano-sedimentary accumulation rates in deep marine back-arc basins, the interpreted setting of the Finlayson Lake district, are dependent on the volcanic style (e.g., deep-vs.shallowlevel eruptions, degree of explosivity; Gibson et al., 1999;Allen and Stewart, 2003;Cas and Simmons, 2018) and the stage of back-arc development (e.g., Carey and Sigurdsson, 1984).Furthermore, the source and eruptive style of volcanoes can greatly influence the quantity of particles available for deposition, especially in submarine environments where hydrostatic pressure and the presence of water affect the separation and distribution of facies generated by pyroclastic eruptions (Cas and Simmons, 2018).Subaerial systems are known to facilitate large-volume silicic eruptions (>10 3 km 3 ) that are observed to deliver volcanic ash to the water column by ash plumes (e.g., Mt.Pinatubo; Bryan et al., 2010, and references therein) or by pyroclastic eruptions that breach the seawater interface (e.g., Soufrière Hills volcano, Montserrat; le Friant et al., 2009).We do not envision subaerial arc eruptions to have had a significant impact on volcanic facies in the units studied here for three primary reasons.First, there is little evidence for shallow water textural features in the volcanic facies (e.g., tractional current-related bedforms, welding) as would be expected of ash columns entering a back-arc basin, suggesting that rocks herein were formed below the storm wave base (Gibson et al., 1999).Second, the rocks do not preserve features consistent with eruption in a subaerial environment, such as accretionary lapilli (Piercey et al., 2001;2016).Lastly, the district contains local-to map-scale features in volcanic and sedimentary facies that indicate facies deposition from turbidity currents and/or debris flows.Examples include normal grading of sand-to mud-rich facies (i.e., Kudz Ze Kayah-Wind Lake formation; Fig. 4), tuff or mudstone bedforms in the Wolverine Lake group footwall (Fig. 5), and alternating tuff-argillite sequences attributed to turbidity currents in the Wolverine Lake group hanging wall (Piercey et al., 2016).These geologic features are much more consistent with having formed during the early-rift stages of back-arc development where steep topography was created on the seafloor, leading to the development of a basin architecture that facilitated abundant submarine eruptions and mass transport complexes or turbidity current-related deposits (Fig. 15; Carey and Sigurdsson, 1984;Clift, 1995).Our new U-Pb dates and associated accumulation rates, together with chemical and regional facies associations, also suggest that the Kudz Ze Kayah formation and Wolverine Lake group were formed in incipient back-arc basins with an early-rift architecture, which then facilitated rapid, submarine volcanism and debris flows and/or turbidity currents, likely coupled with rapid subsidence (Fig. 15).Together, these characteristics suggest that the facies that directly host abundant VMS mineralization were deposited rapidly in an evolving back-arc basin.
Rapid emplacement of volcanic facies on the seafloor facilitates the formation of a water-saturated interface where unconsolidated sediments and rock fragments dominate (Cas and Wright, 1987;Doyle and Allen, 2003).This water-saturated interface facilitates a chaotic volcanic environment if subsequent lavas and/or subvolcanic intrusions are emplaced, which would then be subjected to quench fragmentation and creation of autoclastic and fractured coherent facies (e.g., Mc-Phie and Allen, 1992;Cas and Simmons, 2018).This layer of unconsolidated sediment should continue to be present with further rapid volcanic and volcano-sedimentary emplacement, especially if large volumes of sediment are supplied to the basin in rapid succession, which appears to have been the case in both the Kudz Ze Kayah and Wolverine basins based on our depositional rate calculations.Contemporaneous hydrothermal upwelling within a seafloor-subseafloor environment like the above leads to conditions that are optimal for deposition of sulfides in both the seafloor and subseafloor setting as the metal-rich hydrothermal fluids mix with seawa-ter in the bottom waters of the ocean, as well as with fluids trapped in pore spaces and/or fractures in clastic and coherent facies, respectively (Doyle and Allen, 2003;Piercey, 2015).We have demonstrated that rapid volcanic emplacement processes have occurred in the Finlayson Lake district and are associated with VMS deposits; however, it is uncertain what their role is in creating replacement-style mineralization as predicted above.Below, we evaluate the role of emplacement rate and sulfide emplacement style using deposits in the Finlayson Lake district.

Subseafloor replacement-style mineralization and potential relationships to accumulation rates
Subseafloor replacement is recognized as an important process in the formation of many VMS deposits globally (Doyle and Allen, 2003); however, in the ancient record, the geologic information required to document these processes can often be obscured due to post-depositional tectonic events, such as metamorphism and deformation (Allen et al., 2002;Piercey, 2015;Tornos et al., 2015).Replacement textures in VMS systems are typically observed as sulfide infiltration into either porous volcano-sedimentary facies or coherent and/or autoclastic lavas or subvolcanic intrusions (Doyle and Allen, 2003).The relatively intact stratigraphy in the Finlayson Lake district hosts at least three Zn-enriched VMS deposits that exhibit various primary textural and facies characteristics of subseafloor replacement-style mineralization processes (e.g., Fig. 16; Peter et al., 2007;Bradshaw et al., 2008).Previous workers have suggested subseafloor replacement was an important process in the genesis of some of the Finlayson Lake deposits (Peter et al., 2007;Bradshaw et al., 2008), and we provide additional evidence herein.For example, at Wolverine, the Zn-Pb-Fe-rich sulfides have textures that resemble partial or complete replacement of carbonaceous argillite and minor felsic tuff beds (Fig. 16A, B), filling fractures in coherent porphyritic intrusions (Fig. 16C), and microscale bud-bybud replacement of framboidal sulfides with metal-rich sulfides (Peter et al., 2007;Bradshaw et al., 2008;Piercey, 2015).Replacement textures in the GP4F deposit include sulfide deposition subparallel to bedding in volcaniclastic rocks and relict lapilli and/or crystal fragments in the sulfide mineralization (Fig. 16D-F; Peter et al., 2007).The Kudz Ze Kayah deposit exhibits many of the same replacement characteristics, including relicts of rhyolite fragments or lapilli in sulfide interpreted to reflect partial to complete replacement of primary volcaniclastic facies by sulfide minerals (Fig. 16G, I-K), fracture infill of sphalerite in coherent rhyolite (Fig. 16H), and the presence of sulfides that have contacts discordant to bedding surfaces (Fig. 16L).Additionally, similar alteration assemblages occur in both hanging-wall and footwall facies of all three deposits (e.g., Peter et al., 2007;Bradshaw et al., 2008;Denisová, unpub. data).
Subseafloor replacement is a key process in the formation of large and high-grade VMS deposits (e.g., Rosebery, Tasmania and Kidd Creek, Ontario;Hannington et al., 1999;Martin, 2004), and recently SMS deposits (e.g., Hakurei Site, middle Okinawa Trough;Nozaki et al., 2021), because it enhances the amount of sulfide precipitated to increase the deposit tonnage, promotes higher degrees of preservation with increased burial in the ancient record, and facilitates zone refining processes during long-lived hydrothermal systems (Doyle and Allen, 2003;Piercey, 2015;Nozaki et al., 2021).All felsic-hosted VMS deposits in the Finlayson Lake district show abundant evidence for subseafloor replacement in both unconsolidated and coherent volcanic facies (e.g., Fig. 16), and our new U-Pb results, along with textural observations, indicate that they were hosted in facies sequences that were emplaced rapidly at both local and basin scales (>1,000 m/m.y.; e.g., Table 2).Rapid deposition of host lithofacies has been proposed as a diagnostic factor for facilitating replacement-style mineralization (Doyle and Allen, 2003); however, this attribute has been identified primarily by textural facies analysis and not with quantitative rates.There is significant variability in timing relative to volcano-sedimentary deposition for replacement-style deposits, where hydrothermal activity can be pre-, syn-, and post-deposition relative to the volcano-sedimentary sequence (e.g., Goodfellow and Franklin, 1993;Allen et al., 1996a, b;Hannington et al., 1999;Thurston et al., 2008;Belford et al., 2015).Subseafloor replacement is recognized at Middle Valley on Juan de Fuca Ridge, where sediment thicknesses and sedimentation rates are well-defined and indicate that hydrothermal sediment is absent below ~80 cm depth; this is interpreted to represent high sedimentation rates and subsidence prior to hydrothermal activity in the active vent sites (Goodfellow and Franklin, 1993).Further, Goodfellow and Franklin (1993) identified slower sediment accumulation rates at bathymetric highs compared to lows, which they suggest is due to shedding of sulfide mounds and clastic sediments into depressions because of tectonic and/or paleoslope instability related to movement on active rift-related normal faults.This relationship of sediment thickness, rate, and mineralization is also evident at a larger scale within the Kudz Ze Kayah formation despite a much larger volcanic component, where the Kudz Ze Kayah and GP4F deposits are spatially associated with the thickest and highest accumulation rates of volcaniclastic facies in the basin (Figs. 4,8).Therefore, this comparison shows that localization of replacement-style deposits is related to both the thickest and most rapid rates of volcano-sedimentary accumulation.
We argue that rapid volcano-sedimentary accumulation rates are required to deposit volcano-sedimentary lithofacies with the infrastructure (i.e., porosity and permeability) for seawater (+hydrothermal fluid) infiltration that can then sustain the formation of large-tonnage and/or high-grade VMS deposits through zone refining processes and enhanced preservation.Continued rapid emplacement of lithofacies, particularly poorly sorted volcaniclastic mass-flows, as observed in the Kudz Ze Kayah and Wolverine deposits and their host stratigraphy, will both bury and preserve existing subseafloor mineralization and provide new substrate that has the potential to be mineralized.Rapid facies emplacement in these basins was likely an important factor for consistently developing thick accumulations of fertile lithofacies at the seafloor interface that are porous, permeable, and water saturated.Therefore, accumulation rates influenced the scale of hydrothermal mixing with seawater and resulted in the increased tonnages and high Zn grades found in all three felsic-hosted VMS deposits in the Finlayson Lake district.
Rapidly emplaced volcanic and volcano-sedimentary facies were previously argued to be important for facilitating subseafloor replacement-style VMS mineralization (Doyle and Allen, 2003); however, documentation of the actual emplacement rates of associated volcanic sequences has heretofore been inadequate and unquantified.Our approach underscores the significance of modern high-precision U-Pb geochronology to resolve rates and timescales of emplacement for VMS-related lithofacies in the ancient record.It is only with this enhanced high-resolution geochronology, coupled with detailed stratigraphic reconstructions and facies analysis, that we can more precisely unravel the governing factors between facies deposition, rates of emplacement, and mineralization styles that define high-value VMS deposits (i.e., high grade and/or tonnage) and most prospective basins in convergent margins globally.

Conclusions
New high-precision U-Pb geochronological constraints on felsic volcanic rocks in the Finlayson VMS district provide well-constrained limits on the age of volcanism and associated VMS mineralization.The Kudz Ze Kayah formation is dominated by volcaniclastic and lesser subvolcanic and sedimentary facies that were erupted or deposited between 363.3 and 362.4 Ma (Dt = 0.6-1.0m.y.).In contrast, the Wolverine Lake group began sedimentary deposition at ca. 357.5 Ma, which then transitioned into a period of voluminous interbedded volcanic and argillite facies between 356.4 and 354.9 Ma (Dt = 1.2-1.5 m.y.).These results indicate that the facies comprising the Kudz Ze Kayah formation and Wolverine Lake group were formed at rapid, time-averaged accumulation rates of ~300 to 1,900 m/m.y.(up to 2,900 m/m.y.) and ~800 to 1,400 m/m.y.(up to 7,800 m/m.y.), respectively, where the highest rates are interpreted to result from mass transport complexor turbidity current-related deposition.Numerous subseafloor replacement-style VMS deposits are observed and documented in each of the geologic units and occur at distinct time intervals (GP4F = 363.3Ma; Kudz Ze Kayah = 362.8Ma; Wolverine = 355.2Ma).We suggest that the thick, rapidly emplaced volcano-sedimentary facies in the Finlayson Lake district created a fertile, water-saturated environment that facilitated subseafloor replacement-style mineralization, rapidly covered the deposits, and thus enhanced their preservation and enabled zone refining in these subseafloor systems, all of which led to the development of VMS deposits with elevated Zn grades and tonnages.Together, high-precision geochronology, stratigraphic reconstructions, and detailed facies analysis can greatly improve our understanding of the rates, timescales, and characteristics of volcanism and associated VMS mineralization that produced the world-class VMS deposits in ancient convergent margins worldwide.

Fig. 4 .
Fig. 4. Geologic map of the Kudz Ze Kayah formation (after van Olden et al., 2019) showing U-Pb geochronology sample locations (yellow circles and black numbers) and lines of cross section lines used to reconstruct stratigraphic sections.Structural data is shown for representative measurements used to construct cross sections for S0 and S1 cleavages only.Red boxes correspond to VMS zones.Unit labels and colors as in Figure 3.
Fig.5.Geologic map of the Wolverine Lake group, merged with regional mapping (after Yukon Geological Survey, 2018).Symbols and labels as in Figures3 and 4.
Fig.5.Geologic map of the Wolverine Lake group, merged with regional mapping (after Yukon Geological Survey, 2018).Symbols and labels as in Figures3 and 4.

Fig. 6 .
Fig. 6.Cross sections for lines A-A', B-B', C-C', D-D', and E-E' in the Kudz Ze Kayah formation, drafted with no vertical exaggeration.Sections correspond to plan view lines in Fig. 4. Unit colors and abbreviations as in the legend in Fig. 3. U-Pb samples are shown as yellow circles for reference and projected into the cross section where applicable (B-B' and C-C'; 17MM-002); numbers correspond to sample symbols in Figures 4 and 5. Thin, black dashed lines in the subsurface are interpreted cleavage plane variations based on outcrop measurements.Red dashed lines indicate approximate location of VMS lenses.Arrows on the surface indicate locations of outcrop measurements corrected as apparent dips.Abbreviations: Ef = East fault, FCf = Fault Creek fault, FLf = Fire Lake formation, GLPS = Grass Lakes plutonic suite, KKT = Krakatoa, KZKf = Kudz Ze Kayah formation, NRf = North River formation, WLf = Wind Lake formation.

Fig. 8 .
Fig.8.Regional stratigraphic columns for A) Kudz Ze Kayah formation, and B) the Wolverine Lake group.Columns show cumulative unit thicknesses, approximate locations of VMS deposits, and U-Pb sample locations (circles = this study; diamonds = data fromMurphy et al., 2006;Piercey et al., 2008).Unit thicknesses are measured from interpreted cross sections in Figures6 and 7. Thicknesses in panel (A) account for fault restoration and subtraction of Grass Lakes plutonic suite granites, as these rocks intruded the volcanic succession following the transition into the overlying Wind Lake formation.Thick dashed lines in the columns represent unconformities; solid lines are locations of facies changes.Bold numbers correspond to sample numbers in Figures4 and 5. Abbreviations: FW = footwall, HW = hanging wall, KKT = Krakatoa VMS zone.Unit colors correspond to Figure1.

Fig. 10 .
Fig. 10.Representative structural fabrics in the Kudz Ze Kayah formation (A-B) and the Wolverine Lake group footwall (C-D).Pervasive S1 cleavage is the dominant planar feature in both geologic units.A) Prominent S1 cleavage in an outcrop of felsic tuff, with a broadly west-east strike that dips to the north.B) Crenulation fabrics in a felsic tuff developed on an earlier S1 cleavage, with a S2 cleavage superimposed highly oblique to S1. C) Intercalated argillite (ARG) and felsic tuff (FLT) in the upper Wolverine Lake group footwall.Laminations of argillite beds represent primary bedding (S0), and S1 cleavage in the FLT unit is subparallel to bedding.D) Quartz-feldspar porphyritic intrusive rock with strong mineral alignment subparallel to the prominent S1 cleavage.

Fig. 13 .
Fig. 13.Concordia diagrams showing chemical abrasion-isotope dilution-thermal ionization mass spectrometery (CA-ID-TIMS) U-Pb geochronology results for zircon in rocks from the Kudz Ze Kayah formation (A-D) and the Wolverine Lake group (E-I).A) Felsic crystal tuff at the Kudz Ze Kayah-Wind Lake formation contact (17MM-002); B) felsic lapilli tuff in the ABM hanging wall (17MM-031); C) felsic tuff from the Krakatoa footwall (18MM-133); D) felsic crystal tuff in the GP4F footwall (17MM-074); E) feldspar-porphyritic rhyolite from the Fisher zone (FP; 18MM-102); F) felsic tuff in the Wolverine/Lynx footwall (P99-WV-4K); G) feldspar-porphyritic rhyolite from the Wolverine/Lynx zone (P00-WV-1C); H) quartzfeldspar-porphyritic rhyolite from the Sable zone (QFP; P00-WV-12); and I) felsic crystal tuff from the lower Wolverine Lake group footwall (18MM-114).Each individual ellipse encompasses the 2s error on a single zircon analytical result.The gray dashed lines represent the range of decay constant uncertainty on the concordia curve.The inset panels show individual grain analyses with their included 2s error (vertical bars); the solid black line indicates the calculated weighted mean 206 Pb/ 238 U date, and the larger gray region indicates the error associated with external reproducibility of the analyses.Ellipse and bar colors correspond to the ranked bars in Fig. 14.Abbreviations: FW = footwall, HW = hanging wall, MSWD = mean square of weighted deviates, Wtd.= weighted.

Fig. 16 .
Fig.16.Textural evidence for replacement-style mineralization in the Wolverine (A-C), GP4F (D-F), and Kudz Ze Kayah VMS deposits (G-L) in the Finlayson Lake VMS district.A) Primary shale fragments and beds that are replaced partially and fully to Zn-Pb sulfides, grading into massive sphalerite-pyrite-galena-sulfosalt assemblages.B) Massive Cu-Fe-rich sulfides that have replaced shale beds and tuff beds and fragments.C) Quartz-feldspar porphyritic rhyolite that has been fractured and replaced by sulfide veins.D) Chlorite-garnet-rich rock that is replaced by pyrite-pyrrhotite-sphalerite-galena assemblages, leaving rhyolite and lapilli clasts behind.E) Crystal tuff that has been variably replaced by sphalerite-pyrite ± pyrrhotite.F) Lapilli tuff with pyrite-pyrrhotite beds subparallel to tuff beds, where lapilli are preserved as fragments in the sulfides.G) Massive pyrite-pyrrhotite that replaces a coherent rhyolite flow.H) Sphalerite in fractures in a rhyolite flow.I) Remnant lapilli and feldspar crystals in massive pyrite-pyrrhotite-sphalerite beds.J) Massive pyrite ± sphalerite that contains remnant rhyolite clasts and tuff beds.K) Remnant lapilli fragments in pyrite-pyrrhotite-sphalerite-galena assemblage, grading down into chalcopyrite-rich section that replaces tuffs.L) Discordant replacement of coherent rhyolite with sphalerite-pyritepyrrhotite.Abbreviations: cpy = chalcopyrite, fsp = feldspar, gal = galena, grt = garnet, po = pyrrhotite, py = pyrite, qtz = quartz, rhy = rhyolite, sph = sphalerite, ss = sulfosalt.
. Results are presented as 206 Pb/ 238 U isotope ratios for all of the CA-ID-TIMS and LA-ICP-MS samples <1 Ga; those grains with dates >1 Ga are reported in the text as 207 Pb/ 206 Pb dates.Unless otherwise noted, grains with young LA-ICP-MS dates that are not included in the weighted mean are interpreted to reflect post-crystallization Pb loss.The CA-ID-TIMS results are reported in Table 1 and Figures 13 and 14; LA-ICP-MS results are presented in Appendix Tables . The CA-ID-TIMS analytical results for seven individual fractions are concordant, contain relatively low Th/U (0.16-0.54), and yield a weighted mean 206 Pb/ 238 U date of 362.404 ± 0.098 Ma (2s; MSWD = 0.55; Table1; Fig.13A).

WV-1C Wolverine/Lynx FP
Pb/ 238 U date for this zircon population is 363.254 ± 0.098 Ma (2s; MSWD = 0.57; Table1; Fig.13D), interpreted to be the eruption age of this rock, the approximate age of the GP4F deposit, and the oldest date in the Kudz Ze Kayah formation.The LA-ICP-MS data is concordant with a range of Th/U (0.13-4.3) and a weighted mean 206 Pb/ 238 U date of 363.1 ± 2.3 Ma (2s; MSWD = 0.67; n = 40); this date is within uncertainty of the CA-ID-TIMS date.Two grains have older dates of ca.383 and 1770 Ma, which are both interpreted as xenocrystic zircon.

Table 1 .
Mattinson (2005)n-ID-TIMS U-Pb Zircon Geochronology Results for Rocks in the Kudz Ze Kayah Formation and Wolverine Lake Group, Finlayson Lake District, Yukon Sample numbers denoted with a * indicate the use of the ET2535 tracer solution, whereas the ET535 tracer was used for all other samples; UTM coordinates in NAD83 Zone 9U; z1, z2 etc. are labels for single zircon grains or fragments annealed and chemically abraded afterMattinson (2005); bold text indicates results used in weighted mean calculations ). Horizontal black bars indicate the weighted mean 206 Pb/ 238 U date.Lighter vertical bars are not included in the weighted mean calculation for each sample and, except for z1 in 18MM-133, are interpreted as xenocrysts.Horizontal pink lines indicate distinct time intervals of VMS mineralization based on 2s uncertainties on the weighted mean date, indicated by the names of deposits from which rocks were sampled (red text).Broken Y-axis scale from 361 to 358 Ma indicates a period of deformation (D), represented by an angular unconformity between the two stratigraphic successions.Bold numbers correspond to sample numbers in Figures 4 through 8 and Table 1.B) Summary of age distributions (Dt) for each sample (blue = Kudz Ze Kayah formation; green = Wolverine Lake group) and the total minimum and maximum durations based on weighted mean 206 Pb/ 238 U dates (WM), 206 Pb/ 238 U dates, and their associated 2s uncertainties.Vertical boxes represent the Dt for each category.Durations given in k.y. Sable?

Table 2 .
Summary of Calculated Volcano-Sedimentary Accumulation Rates