The Medicine Hat Block and the Early Paleoproterozoic Assembly of Western Laurentia

: The accretion of the Wyoming, Hearne, and Superior Provinces to form the Archean core of western Laurentia occurred rapidly in the Paleoproterozoic. Missing from Hoffman’s (1988) original rapid aggregation model was the Medicine Hat block (MHB). The MHB is a structurally distinct, complex block of Precambrian crystalline crust located between the Archean Wyoming Craton and the Archean Hearne Province and overlain by an extensive Phanerozoic cover. It is distinguished on the basis of geophysical evidence and limited geochemical data from crustal xenoliths and drill core. New U ‐ Pb ages and Lu ‐ Hf data from zircons reveal protolith crystallization ages from 2.50 to 3.28 Ga, magmatism/metamorphism at 1.76 to 1.81 Ga, and ε Hf T values from − 23.3 to 8.5 in the Archean and Proterozoic rocks of the MHB. These data suggest that the MHB played a pivotal role in the complex assembly of western Laurentia in the Paleoproterozoic as a conjugate or extension to the Montana Metasedimentary Terrane (MMT) of the northwestern Wyoming Province. This MMT–MHB connection likely existed in the Mesoarchean, but it was broken sometime during the earliest Paleoproterozoic with the formation and closure of a small ocean basin. Closure of the ocean led to formation of the Little Belt arc along the southern margin of the MHB beginning at approximately 1.9 Ga. The MHB and MMT re ‐ joined at this time as they amalgamated into the supercontinent Laurentia during the Great Falls orogeny (1.7–1.9 Ga), which formed the Great Falls tectonic zone (GFTZ). The GFTZ developed in the same timeframe as the better ‐ known Trans ‐ Hudson orogen to the east that marks the merger of the Wyoming, Hearne, and Superior Provinces, which along with the MHB, formed the Archean core of western Laurentia. ‐ 206 22 grains 207 206 Pb 2.84 ± Ga Younger, concordant ages from both cores and rims The error ‐ weighted mean age of 2.84 Ga is interpreted to be the best estimate of the minimum crystallization age of the protolith.


Introduction
The Precambrian Medicine Hat Block (MHB) lies south of the Archean Hearne Province ( Figure  1), and it was originally defined by structural features in overlying Paleozoic sedimentary rocks [1] and subsequently by potential field and seismic surveys [2,3]. Although its history is poorly constrained due to extensive Phanerozoic sedimentary cover (Figure 1), an Archean age for the MHB was first proposed on the basis of multi-grain U-Pb zircon ages acquired by ID-TIMS (isotope dilution-thermal ionization mass spectrometry) from five samples of drill core recovered from oil and gas exploration wells [4]. U-Pb ages of these samples range from approximately 2.6 to 3.3 Ga [4]. With its long history, the MHB clearly occupied an important position between the Wyoming Craton and the Hearne Province during the formation of the supercontinent Laurentia.
In this paper, we focus on refining the evolution of the northwestern crystalline basement of the MHB and western Canada sedimentary basin and relate that history to that of other terranes that accreted to Laurentia prior to 1.7 Ga. We report results of 484 new single zircon U-Pb ages obtained by laser ablation multi-collector ICP-MS from drill core and xenoliths previously described and analyzed by Villeneuve et al. [4] and Davis et al. [5] in order to extract the maximum information from these unique samples, including all new in situ Lu-Hf analyses by the same method. The xenoliths are from Eocene minettes of the Sweetgrass Hills area of the Montana alkali province (MAP) located within southern Alberta and northern Montana (Figure 1; [6]). These alkaline volcanicentrained xenoliths provide a unique opportunity to study the age and history of the enigmatic continental crust of the northwestern MHB (e.g., [5,7]) and the buried southern (Alberta) part of the western Canada sedimentary basin [8][9][10][11][12][13]. Constraining the relationships of the Medicine Hat block to the adjacent Paleoproterozoic belts (e.g., Great Falls Tectonic Zone (GFTZ), Trans-Hudson orogen (THO), and Vulcan low) is important for understanding the geodynamics of the rapid, Paleoproterozoic assembly of Laurentia proposed by Hoffman [14].

Geologic Background and Previous Work
The boundary separating the MHB and the Hearne Province to the north is the Vulcan structure, which is a prominent, linear feature characterized by both gravity and magnetic lows (Figure 1; [9]). The northern limit of the MHB in southern Alberta and northern Montana, including the Vulcan structure, is entirely buried beneath the sediments filling the Western Canada sedimentary basin, which is a Cordilleran foreland basin. On its eastern margin, the MHB is bounded by Paleoproterozoic rocks of the Trans-Hudson orogen (Figure 1; [20][21][22]). To the south, at least parts of the MHB lie within the Great Falls tectonic zone (GFTZ) that formed during the Great Falls orogeny (1.7-1.9 Ga; [23]). The GFTZ also marks the northern limit of the Archean Beartooth-Bighorn magmatic zone of the Wyoming Province [24,25]. Proposals for the origin of the GFTZ have ranged from Archean and Paleoproterozoic suture zones to an intra-continental shear zone [2,[26][27][28].

U-Pb Geochronology of Zircon
Zircon grains provided by the Geological Survey of Canada were prepared using traditional methods of crushing and grinding, followed by density separation by panning, high-density liquids, and by magnetic susceptibility using a Frantz magnetic separator. Hand-picked zircon grains were embedded in a 25 mm epoxy disk, together with fragments of a natural zircon standard FC-1 [70][71][72]. The mounts were polished to approximately half thickness of the individual zircon grains, and imaged using cathodoluminescence (CL). In addition to the 484 U-Pb analyses, 106 Lu-Hf analyses were obtained from selected zircon domains identified in the CL images for each zircon. The zircon mounts were first analyzed for U-Pb, followed by Lu-Hf analysis in the same CL domain as the U-Pb analysis, but conducted in a separate session.
The U-Pb ages were determined using the Nu Plasma laser ablation-multi-collector-inductively coupled plasma mass spectrometry (LA-MC-ICP-MS) at the University of Florida (Gainesville, FL, USA) following methods in Mueller et al. [73]. The analyses utilized a New Wave 213 nm laser and a 20 μm spot diameter. Each individual U-Pb analysis consisted of 30 seconds of measurement during ablation in order to minimize ablation pit depth and elemental fractionation. Isotopic data were acquired using the Nu-Instruments Time Resolved Analysis software, which allowed isotopic ratios to be calculated from data collected simultaneously in individual faraday and ion-counting channels during the 30 second ablation. Fractionation and drift corrections were calibrated against multiple ablations of the FC-1 natural zircon standard (Duluth Gabbro; 1098 Ma; [70][71][72]); each batch of 10 unknowns was bracketed between two pairs of FC-1 analyses. The zircon isotopic data were reduced with the in-house program hosted in Microsoft Excel (CALAMARI, © P.A. Mueller), with ages and concordia diagrams calculated using Isoplot v. 4.0 [74,75]. Representative age errors based on the long-term reproducibility of FC-1 are 2% for 206 Pb/ 238 Pb (standard error of the mean (s.e.m.)) and 1% for 207 Pb/ 206 Pb (s.e.m.). U-Pb data, including for FC-1, are provided in the Supplementary Data (Tables S1,S3).

Lu-Hf Isotopic Analysis of Zircon
Lu-Hf analyses were also conducted by laser ablation with the MC-ICP-MS in static mode following the methods in Mueller et al. [73]. Faraday collectors were used for the simultaneous measurement of 180 Hf, 178 Hf, 177 Hf, 176 Hf, 175 Lu, 174 Hf, and 172 Yb [73]. The analyses were performed with online Lu and Yb isobaric interference corrections, using 176 Lu/ 175 Lu = 0.02653 and 176 Yb/ 172 Yb = 0.5870, both of which are within the range of published values [76,77]. The non-radiogenic isotopic ratio of 180 Hf/ 177 Hf was measured to monitor mass bias corrections determined using 177 Hf/ 179 Hf and determined to be within the error of the true value (1.88710) for all unknowns and standards. Analyses that did not meet this criterion were discarded. Corrections for Lu and Yb interferences (peak stripping) were done on-line and based on mass bias-corrected ratios using the mass bias factors derived from Hf isotopes. This method was chosen because of the low reproducibility of Yb isotope ratios. Although mass bias for Yb, Lu, and Hf may have an elemental component to them, we do not utilize 176 Hf/ 177 Hf data that have combined Yb + Lu corrections in excess of 25%. Analyses with Yb + Lu corrections exceeding 20% are not used to define the limits of model ages or initial εHf for any suite of zircons. All Hf, Lu, and Yb isotopic ratios were corrected for mass bias using 178 Hf/ 177 Hf = 1.46718. Each analysis consisted of a 40 μm spot size and 20-second background measurement, up to one minute of measurement during ablation, and a 30-second purging period between analyses. Multiple analyses of FC-1 (Duluth Gabbro zircon standard; [71,78]) yielded 176 Hf/ 177 Hf = 0.28217 (± 0.00002, 2σ, n = 9), which was within the error of solution analyses of this standard ( 176 Hf/ 177 Hf = 0.282174 ± 0.000013, 2σ; [73]) as well as data published by Woodhead and Hergt (2005;176 Hf/ 177 Hf = 0.282172 ± 0.000042, 2σ). The reproducibility of 176 Hf/ 177 Hf for FC-1 analyses in epsilon notation without age correction was -21.8 +/−1.1 εHf (2σ), which is within the error of reported values and the Florida lab's measurements of FC-1 by wet plasma [73,78]. The standard data were obtained during the same session as the unknown samples (1 standard, 15 unknowns, 1 standard). Measured and online mass-bias-corrected present-day 176 Lu/ 177 Hf ratios were utilized to calculate initial 176 Hf/ 177 Hf ratios using the decay constant of Lu (λ = 1.867e −11 ; [79,80]) and the 207 Pb/ 206 Pb age of the analytical spot, following established procedures [81,82]. Overall, because of the very low Lu/Hf ratios (FC-1: 176 Lu/ 177 Hf = 0.001329 ± 0.000091, 2σ, n = 9), the difference between the present-day measured and calculated initial 176 Hf/ 177 Hf ratios in most cases is < 1 ε unit (e.g., [73]).
Zircons were selected for Lu-Hf analysis on the basis of low discordance (<10% 206 Pb/ 238 U versus 207 Pb/ 206 Pb), zonation patterns, and grain size in order for the Lu-Hf analyses (40 μm spot) to be in the same domain as the U-Pb analysis (20 μm spot). Measuring Hf isotopes on the same zircons that are dated has the advantage of estimating the likelihood of open system behavior: Zircons that have experienced no, or only very little, Pb loss (i.e., concordant or nearly concordant) are the least likely to have suffered open system behavior of either Lu or Hf. The 176 Hf/ 177 Hf ratio is commonly represented as εHf, which is the 176 Hf/ 177 Hf in the unknown sample minus the 176 Hf/ 177 Hf in chondritic uniform reservoir, and then that value divided by 176 Hf/ 177 Hf in the chondritic uniform reservoir times 10,000 [83]. The chondritic uniform reservoir (CHUR) values are those of Bouvier et al. [84]: 176 Lu/ 177 Hf = 0.0336 and 176 Hf/ 177 Hf = 0.282785. Initial ratios are given as 176 Hf/ 177 Hf(T) and εHf(T) for all analyzed samples and were calculated using the 207 Pb/ 206 Pb crystallization age of the zircon or its host rock and the measured Lu/Hf ratio, (e.g., initial Hf at 2.4 Ga = εHf(2.4 Ga)); propagated net uncertainties for initial ratios are approximately 1.  [85] and Sanfilippo et al. [86]. The mean values for 176 Hf/ 177 Hf(T) and εHf(T) (and related errors (± 2σ)) are summarized in Table 1; details of the specific grain analyses are available in Supplementary Table S2. FC-1 standard Lu-Hf isotopic data are available in Supplementary Data Table S4. Note that some grains were too small for both U-Pb and Lu-Hf analyses; only age data are reported for these grains.
In addition to providing insight into the role of mantle and crustal reservoirs in the development of the magmas from which the zircons crystallized, the range of initial εHf values can also be used to identify zircons that may not have crystallized from that magma, i.e., xenocrysts. Although there is no fixed range of initial εHf that can be used to constrain a felsic magma generated in a terrane with very old crust, the viability of the age can be enhanced if those ages lead to initial εHf values that cover a limited range and that range of values is comparable to the εHf values measured in the zircons today. The limits are unavoidably ad hoc to some degree, but comparisons with the ranges of initial εHf values determined for rocks of similar age in the region provide useful benchmarks. For example, the approximately 2.8 Ga Long Lake magmatic complex of the Beartooth Mountains has a range of initial εHf of 15 epsilon units [39,87], and the Paleoproterozoic igneous rocks of the Little Belt Mountains range over 23 epsilon units [88]. Measurement errors (i.e., ± 2σ) are inconsequential compared to the ranges of initial ratios.  Grains analyzed by ICP-MS for this study were selected from the zircon concentrates provided by the Geological Survey of Canada and previously described as part of earlier multi-grain ID-TIMS studies [4,5]. Here, individual analyses were first evaluated as acceptable if they met two criteria: (1) The 207 Pb/ 206 Pb age was < 5% discordant compared to the 206 Pb/ 238 U age, and (2) common Pb (i.e., mass 204 (Pb + Hg)) cps were less than the measured cps at mass 204 for the FC-1 standard. The values for 204 Pb cps and 204 Pb/ 206 Pb for the analyses and standards are available in the supplementary data. CA-TIMS analyses [72] showed that typical FC-1 zircons have 204 Pb/ 206 Pb < 0.001. All analyses were plotted on a Tera-Wasserburg (TW) concordia diagram using the Isoplot program of Ludwig [74,75] to examine trends in discordance (Figures 3-5). The probability distribution function in Isoplot was used to identify age populations. Plots of 207 Pb/ 206 Pb ages versus percent discordance provided further insight into the behavior of age populations (Supplementary Table S1). These plots assisted in selecting grains for weighted average 207 Pb/ 206 Pb age calculations by clarifying relationships between age populations and discordance. Grains were selected based on their near concordance and limited range of 207 Pb/ 206 Pb ages as demonstrated within the plot. The calculation of a useful concordia intercept age was hindered in some samples by multi-stage Pb-loss. In these cases, we conducted a second evaluation (described with each sample) and used error-weighted mean calculations (Isoplot v. 4.3, [74,75]) of low-discordance 207 Pb/ 206 Pb ages as the best minimum age estimate for the crystallization of the protolith. The external reproducibility of each individual zircon analysis ( 207 Pb/ 206 Pb age) is approximately 1% (2σ) on the basis of reproducibility of the standard (FC-1; See supplementary data). Error error-weighted mean calculations [74,75] U-Pb age estimates (where applicable) are reported in Table 2.

Home Pacific Knappen
Cathodoluminescence (CL) imaging revealed oscillatory zoned zircon grains typical of magmatic crystallization, with some grains exhibiting core-and-rim structures and fractures ( Figure  2). Grains range in size from approximately 30 μm to 400 μm (long axis). Of the 90 analyses, discordance ranged to 90% (Supplementary Table S1). Seven low discordance grains (<5%) yielded an error-weighted mean age of 3.28 ± 0.01 Ga (MSWD = 2.3), which is interpreted as the age of protolith crystallization (Figure 3, Table 2).  [73]. e εHf values were calculated using the individual U-Pb age for each grain.

Imperial Calstan Lake Newell
Zircon grains in this sample were mostly fragmented and dark in CL images, with limited observable zonation. Grains range in size from approximately 40 μm to 110 μm (long axis). The probability distribution function revealed two age populations in the total of 69 analyses: An older population of 31 grains at ca. 2.74 Ga, and a population of grains ≤ 2.7 Ga. The youngest two grains yielded ages of 1.76 ± 0.01 Ga and 1.91 ± 0.01 Ga, which may reflect zircon growth during later tectonothermal events [56]. The ca. 2.74 Ga analyses were further evaluated using the plots comparing percent discordance versus apparent 207 Pb/ 206 Pb ages grains, revealing a near concordant cluster of grains with discordance values between -2.1% and +2.2%, and 207 Pb/ 206 Pb ages from 2735 to 2755 Ma. This population of 20 grains yields an error-weighted mean age of 2.74 ± 0.03 (MSWD = 1.7; Figure 3, Table 2), which was interpreted as a crystallization age.

PCA Calstan Parkland
Most zircon grains in this sample are fragmented, but singly and doubly terminated grains are present as well. Grains range in size from approximately 60 to 200 μm (long axis). CL images reveal many grains with oscillatory growth zoning, in some cases also showing thin overgrowths or dark cores ( Figure 2). Data were omitted based on the filtering criteria and occurrence along a Pb-loss trajectory observed on the concordia plot ( Figure 3). The data were further evaluated in age versus percent discordance space, which revealed a population of grains with percent discordance values from −1.2 to + 0.9, and 207 Pb/ 206 Pb ages from 2634 to 2656 Ma. This population of 16 grains yield an error-weighted mean age of 2.64 ± 0.02 (MSWD = 0.88; Figure 3, Table 2). This age is interpreted as the minimum crystallization age, with the < 2.624 Ma ages possibly reflecting ancient Pb loss or total to partial metamorphic recrystallization.

PCP Medicine Hat
Zircon grains in this sample are often fragmented, but some elongate, rounded grains are present, ranging in size from approximately 50 to 170 μm (long axis). All grains were dark in CL images, with limited observable zonation (Figure 2). Out of the 70 individual zircon grains analyzed, only 11 pass the filtering criteria (Figure 4), and they range in age from approximately 2.70 to 2.78 Ga. This age distribution is interpreted to represent xenocrystic zircon in the protolith described as a granodiorite gneiss, with possible later overgrowths too small to be analyzed.

DRA-93-26
Zircon grains from this xenolith include stubby, rounded morphologies and some blockier grain fragments, ranging in size from approximately 70 to 160 μm (long axis). CL images reveal apparent oscillatory zoned cores and unzoned dark CL overgrowths in some grains (Figure 2). Other grains display simpler oscillatory growth zoning (Figure 2). Many analyses that pass the filtering criteria lie on a discordia with intercepts at approximately 1.79 Ga and 2.86 Ga; however, no grains with 207 Pb/ 206 Pb ages older than 2.84 Ga are observed ( Figure 5). The TW plot of this sample reveals an older age cluster at ca. 2.83 Ga, with a young tail (discordance >5%) trending down to approximately 2.45 Ga ( Figure 5; Supplementary Table S1). We further evaluated analyses for an error-weighted mean age, favoring grains yielding percent discordance values between −1.6 and +2.5, and apparent 207 Pb/ 206 Pb ages ca. 2.82-2.84 Ga. The 22 grains passing these criteria yielded a 207 Pb/ 206 Pb age of 2.84 ± 0.02 Ga (MSWD = 1.2; Figure 5, Table 2). Younger, concordant ages from both cores and rims likely reflect Pb loss due to Paleoproterozoic tectonothermal event(s). The error-weighted mean age of 2.84 Ga is interpreted to be the best estimate of the minimum crystallization age of the protolith.

DRA-93-01
Zircon grains from this xenolith are stubby to elongate, often with missing or rounded tips, ranging in size from approximately 50 to 180 μm (long axis). CL images reveal oscillatory growth zoning that, in some cases, is overgrown by more homogenous rims (Figure 2). Thirteen of 59 analyses were rejected based on the initial filtering criteria. The three oldest concordant analyses from zircon cores yielding ages of ca. 2.69 Ga are interpreted as xenocrystic. The Paleoproterozoic aged analyses were further evaluated using plots comparing percent discordance versus apparent 207 Pb/ 206 Pb ages. We favored analyses yielding percent discordance values between −0.6 and + 2.5. Sixteen analyses from zircon cores passed these secondary criteria, and 15 grains yielded an error-weighted mean age of 1.81 ± 0.05 Ga (MSWD = 2.1; Figure 5, Table 2). Petrographic analysis indicates a high-grade mineral assemblage (i.e., garnet-clinopyroxene-pigeonite-rutile +/− hornblende +/− quartz assemblage; [5,7,89]), suggesting that the original magma was emplaced into the lower crust (i.e., a charnockite).

DRA-93-20
Zircon grains from this xenolith are generally small and round to oval in shape, ranging in size from approximately 40 to 140 μm (long axis). CL images reveal small, oscillatory zoned cores in some grains, with thicker and rounder overgrowths. Other grains tend to show patchy zonation or little visible zonation. Four analyses were omitted based on the basic filtering criteria ( Figure 5). Analyses were further evaluated using plots comparing percent discordance versus age (Supplementary Table  S1). We favored analyses yielding percent discordance values between −1.2 and +2.5, and apparent 207 Pb/ 206 Pb ages ca. 1.75-1.78 Ga. The youngest coherent population of eight analyses passing these criteria were used to determine an error-weighted mean age of 1.76 ± 0.06 Ga (MSWD = 2.1; Figure 5, Table 2). This age is interpreted as the crystallization age for this sample.

PCP Medicine Hat (PCPMH)
The sample yielded εHf(0) values ranging from −54.6 to −60.9 (n = 15; 6.3 ε-unit range) ( Table 2). Due to the lack of a coherent crystallization age for this sample, εHf(T) values are calculated to the individual zircon age. These values range from −2.9 to 5.3 (8.2 ε-unit range).

Discussion
The new, single grain U-Pb zircon ages presented here largely support previous multi-grain analyses, mostly concordia upper intercept ages, reported by Villeneuve et al. [4] and Davis et al. [5] for the same zircon concentrates. The Home Pacific Knappen core yields zircons with reliable ages of approximately 3.28 Ga, while the remainder of the samples yielded slightly older ages than previously cited (Tables 1 and 2), with the exception of DRA-93-20, which yielded a slightly younger age than previously published (Tables 1 and 2; [4,5]). Collectively, the Archean crystallization ages presented here and by Davis et al. [5] and in other samples analyzed by Villeneuve et al. [4] confirm the widespread presence of Mesoarchean crust in the MHB and southern Hearne Province.
Whole-rock Sm-Nd isotopes of two samples discussed in Villeneuve et al. [4] yielded initial εNd values of + 0.6 and −0.1 (PCP Travers and HPK respectively), crustal residence, or mantle extraction ages of 2.86 and 3.48 Ga, respectively, and 147 Sm/ 144 Nd ratios not typical of continental crust ( 147 Sm/ 144 Nd = 0.1532; [90]). The Sm-Nd isotopic data of Villeneuve et al. [4] are similar to the data of Frost and Burwash [91] for the MHB and southern Hearne Province with crustal residence ages ranging from approximately 2.6 Ga to 3.  Figure 6; current ratios, as well as depleted mantle model ages, are shown in Table 2. There is evidence in the Lu/Hf ratios of the Mesoarchean zircons, albeit limited, of Archean crust with components >3.2 Ga within the MHB (Figure 6; this study; [4,5,56,57]), which indicates that there might be parts of the MHB not yet discovered that record an older history. These possible exceptions aside, the data summarized here strongly suggest that the MHB is similar to other Archean provinces around the world in containing a large amount of crust created between approximately 2.8 and 2.5 Ga [4,5,66]. This conclusion is compatible with interpretations of extant whole rock Sm-Nd isotopic data that yield Archean crustal residence ages for all of the xenoliths and basement cores discussed here and elsewhere [4,5,91]. Similarly, εHf values of zircons for the drill cores and xenoliths range from DM to more negative epsilon values, indicating the mixing and reworking of older crustal materials. The Lu-Hf TDM of approximately 3.6 Ga for zircons from the HPK core (Table 2) indicates that older material equivalent to the MMT [32,[92][93][94] was likely present in the MHB during development of the Neoarchean and Paleoproterozoic crust. This buried Meso-to Paleoarchean crust is a good candidate for the older contaminating component that contributed to lowering the εHf values of the Neoarchean samples (e.g., Figure 6 lower crustal evolution path). Similarly, the Sm-Nd isotopic values for the Mesoarchean MMT crust also involve the mixing of juvenile and older crustal material, although showing a larger range in ε-units (approximately 8.5 Nd ε-units; + 4.9 εNd to −4.3 εNdT values [49,95]).
Although the extent of the Wyoming province, Medicine Hat block, and Hearne provinces and their temporal relationships to the Superior Province across the Trans-Hudson orogen are still debated, the similarity of ages and isotopic compositions between the MHB and the northern Wyoming Province, particularly the MMT, suggests that these two terranes may have shared at least a partial history ( Figure 7A). If so, the MMT-MHB connection likely existed as far back as the Mesoarchean, strengthening proposals for a connection between the MHB and Wyoming Province [6]. The connection was likely broken sometime during the earliest Paleoproterozoic with the formation of a small ocean basin (Medicine Hat ocean) that separated the MMT and other parts of the northern Wyoming Province from the MHB and possibly related terranes such as the Clearwater and Priest River (e.g., [96][97][98][99]). Closure of this ocean by north dipping subduction [28] led to the formation of the Little Belt continental arc on MHB crust beginning at approximately 1.9 Ga [58]. The MHB and MMT re-joined at this time as they amalgamated into the supercontinent Laurentia during the Great Falls orogeny (1.7-1.9 Ga), which formed the GFTZ [1,28,[56][57][58][59]. The GFTZ developed in the same timeframe as the better known Trans-Hudson orogen [60][61][62][63][64][65][66][67] to the east that marks the merger of the Wyoming, Hearne, and Superior Provinces, which along with the MHB formed the Archean core of western Laurentia. The impetus for closure of the Medicine Hat ocean, formation of the GFTZ, and amalgamation of the Superior, Hearne, and Wyoming Provinces remains unclear, but the presence of the approximately 1.8 Ga Rimbey arc immediately northwest of the MHB/Vulcan zone and its proposed relation to the collision of Australia and Laurentia (e.g., [100,101]) coincides with ages presented here and high-grade 1.75-1.85 metamorphism of the western Hearne Province [102]. This collisional scenario is complex and associated with the formation of Nuna-Columbia. Although there are several extant proposals for the sequence of events that led to the formation of a combined Nuna-Columbia-Laurentia scale supercontinent in the Paleoproterozoic, many authors have proposed that a firm connection was established between Laurentia and Australia in the interval 1.8-1.6 Ga (e.g., [103][104][105]). Establishing this connection involved at least one continent-continent or craton-craton collision, which would be of the appropriate scale to induce a wide area of modern northwestern North America characterized by 1.7-1.8 Ga magmatic crystallization ages and metamorphic cooling ages.  [4,5,[56][57][58][59]95,98,99], all of which is being reported here. For reference, ages are plotted compared to U-Pb ages from this study. Ages interpreted as times of igneous crystallization are represented as squares, and potential metamorphic ages are represented as circles. Ages from this study are solid symbols, and ages from the literature are open symbols. (B) The approximate locations and relative positions of the Wyoming Province, MHB, Hearne Province, GFTZ (tightly dashed lines), Trans-Hudson, and Vulcan low (represented by the widely spaced waved lines) [14,24,25,91,109] with corresponding ages cited above.

Conclusions
The data presented here help expand our understanding of the age and tectonic evolution of the largely concealed Medicine Hat block. The new U-Pb ages presented generally agree with previous Archean and Paleoproterozoic age determinations for the region [4,5]; however, they help resolve complexities caused by discordance only hinted at by previous analyses of these samples. The ages generated in this study tend to be slightly older than the earlier multi-grain analyses, which may have been more affected by Pb loss or the incorporation of later metamorphic zircon growth in multi-grain dissolutions. Lu-Hf isotopic analyses paired with single-grain U-Pb data help further constrain the crustal evolution of the MHB, reinforcing that it contains a large amount of crust created between approximately 2.8 and 2.5 Ga. The initial εHf values for the majority of the samples from this study range from DM, indicating juvenile materials, mixed with reworked older Meso-to Paleoarchean crustal materials. The samples reveal similarities to the ages and isotopic compositions of the MMT, suggesting a connection between the MHB and the MMT in their history. This connection was likely broken sometime in the earliest Paleoproterozoic, leading to the formation of a small ocean basin that later closed, resulting in the Little Belt continental arc at approximately 1.9 Ga. The prevalence of Paleoproterozoic ages across the MHB, GFTZ, and MMT (approximately 1.9 Ga to 1.7 Ga) record the entire Great Falls Orogeny, and time that the combined MHB-WY-Hearne crust was amalgamated with the Superior Province to form the Archean core of Laurentia.