Abstract

Continent-continent collisional orogens are the hallmark of modern plate tectonics. The scarcity of well-preserved high-pressure granulite facies terranes minimally obscured by later tectonic events has limited our ability to understand how closely Archean tectonic processes may have resembled better-understood modern processes. Here we describe Neoarchean gneisses in the Teton Range of Wyoming, USA, that record 2.70 Ga high-pressure granulite facies metamorphism, followed by juxtaposition of gneisses with different protoliths, and then by intrusion of leucogranites generated through decompression melting in response to post-collisional uplift. This evidence is best explained as the result of a 2.70–2.68 Ga Himalayan-style orogeny, and suggests that, although subduction may have been occurring earlier in the Archean, doubling of continental thickness by continent-continent collisions may date back to at least 2.7 Ga.

INTRODUCTION

The nature of tectonic processes affecting the early Earth is an enduring and controversial topic (Van Kranendonk, 2010; Hamilton, 2011; O’Neill et al., 2016; Roberts et al., 2015; Smart et al., 2016; Tang et al., 2016). One of the most significant obstacles to unraveling the early history of tectonic processes on Earth is the scarcity in both time and space of representative rock samples; the situation becomes more limiting as time before present increases. Given this fundamental fact, studies of Archean rocks produced by tectonic processes of that era and minimally affected by later tectonism and metamorphism are of primary interest in the study of these phenomena.

Because tonalite-trondhjemite-granodiorite (TTG) and granite-greenstone terrains are distinctive Archean assemblages, they have been the subject of considerable study. However, their tectonic interpretation remains unsettled, with some workers concluding that they may form in a static lid or plume-dominated regime whereas others call on plate tectonics to supply the water needed for their formation (i.e., Condie, 2016). In this study, we look instead at Archean high-pressure granulites as indicators of tectonic regime. In the Phanerozoic, granulite-facies crustal rocks that preserve a record of transport from the surface to depth along a clockwise P-T-t path are diagnostic of collisional plate tectonic processes that double the thickness of the crust. However, high-pressure granulites are uncommon in the Archean. A compilation by Brown and Johnson (2018) of more than 450 localities with well-determined metamorphic conditions and ages identifies very few occurrences of Archean rocks that record pressures of 10 kbar or greater (Table 1). Most are from relatively small outcrops and commonly they have been overprinted by subsequent metamorphic events, making their tectonic significance difficult to establish. All but two are younger than 3 Ga; in one of these two, the Barberton terrane, temperatures do not reach granulite facies.

TABLE 1.

SUMMARY OF HIGH-PRESSURE ARCHEAN GRANULITE TERRAINS (AFTER BROWN AND JOHNSON, 2018)

In this contribution we evaluate the geological record of an Archean high-pressure granulite terrain in the Wyoming Province, North America. No Proterozoic or Phanerozoic metamorphism or deformation has affected these rocks, making the Teton Range a suitable place to identify Archean tectonic processes. Archean rocks of the northern Teton Range include metapelitic gneisses that record 2.70 Ga high-pressure granulite facies metamorphism followed by tectonic juxtaposition with gneisses of different protoliths and geologic histories, and then by intrusion of leucogranites. This evidence is best explained by doubling of crustal thickness through continent-continent collision followed by exhumation, decompression melting, and intrusion of 2.68 Ga leucogranitic sills.

Field and structural relations, petrology and geochemistry, Nd isotopic compositions, and U-Pb zircon geochronology allow us to describe the origin and evolution of this well-preserved assemblage. This terrain preserves evidence of both crustal thickening (based on high-pressure granulite facies metamorphic conditions) and exhumation as evidenced by later, amphibolite facies metamorphism. Zircon geochronology brackets the time for the whole cycle to occur. Although the terrain is only 50 × 15 km in size, it is an excellent place to study Archean tectonic processes because there are rocks that are sensitive recorders of metamorphic conditions and there are no significant younger tectonic overprints to complicate the interpretations.

GEOLOGIC BACKGROUND

The Teton Range, located in northwestern Wyoming, USA, exposes some of the westernmost outcrops of the Archean Wyoming Province (Fig. 1). The northern portion of the Teton Range is underlain by a series of variably deformed paragneisses and orthogneisses (Reed, 1973) that are intruded by dikes of undeformed 2547 ± 3 Ma peraluminous leucogranite known as the Mt. Owen Quartz Monzonite (Zartman and Reed, 1998). The metamorphism and deformation described in this contribution are thus required to be older than the 2547 Ma Mt. Owen Quartz Monzonite. The geologic map of the northern Teton Range (Fig. 2) is derived from the map of Reed (1973). Reed recognized three gneiss units: the Layered Gneiss, the leucogranitic Webb Canyon Gneiss, and Augen Gneiss. A portion of the Layered Gneiss of Reed consists of mafic gneiss and pelitic gneiss, both of which show evidence of high-pressure granulite metamorphism. Because these rocks are so distinctive, we separate them out as the Moose Basin gneiss. (We do not capitalize “gneiss” here because this unit has not been formally recognized on U.S. Geological Survey maps [Frost et al., 2006]).

Figure 1.

Precambrian geology of Wyoming and adjacent states showing the location of the Teton Range and its relation to other Archean exposures in the Wyoming province. Precambrian rocks are exposed in the cores of Laramide uplifts (patterned). The southwestern extent of the province is uncertain. The northeastern boundary of the Wyoming province is after Worthington et al. (2015).

Figure 1.

Precambrian geology of Wyoming and adjacent states showing the location of the Teton Range and its relation to other Archean exposures in the Wyoming province. Precambrian rocks are exposed in the cores of Laramide uplifts (patterned). The southwestern extent of the province is uncertain. The northeastern boundary of the Wyoming province is after Worthington et al. (2015).

Figure 2.

Geologic map across the northern Teton Range showing the relations between the high-pressure granulites (Moose Basin gneiss), leucogranites (undifferentiated Webb Canyon and Bitch Creek gneisses of Frost et al., 2016) and Layered Gneiss. Modified after Love et al. (1992). Numbers identify the samples cited in this paper. Two samples of Moose Basin gneiss analyzed for Sm-Nd isotopic compositions are located south of the map area (08T10 and 08T25).

Figure 2.

Geologic map across the northern Teton Range showing the relations between the high-pressure granulites (Moose Basin gneiss), leucogranites (undifferentiated Webb Canyon and Bitch Creek gneisses of Frost et al., 2016) and Layered Gneiss. Modified after Love et al. (1992). Numbers identify the samples cited in this paper. Two samples of Moose Basin gneiss analyzed for Sm-Nd isotopic compositions are located south of the map area (08T10 and 08T25).

A characteristic feature of the Moose Basin gneiss is large kyanite blades in the pelitic gneiss, which are commonly more than a cm in length (Fig. 3A). Structurally beneath the pelitic gneiss is a mafic gneiss that locally retains granulite assemblages (Fig. 3B). In a few places, the mafic gneiss is cut by garnet-bearing leucosomes (Fig. 3C). The areal extent of the Moose Basin gneiss is reasonably well constrained. Kyanite is found as far south as where upper Bitch Creek crosses the Teton gneisses near the southern margin of Figure 2. Highly retrograded pelitic gneisses occur along the Teton divide at the head of Leigh Creek, ∼5 km south of the southern margin of Figure 2. These occurrences also may be part of the Moose Basin gneiss. Isolated, highly retrograded pelites occur in Osprey Canyon, on the eastern slope of the range (sample 07T32). The eastern margin of the Moose Basin gneiss is a complex contact zone ∼100 m wide consisting of interlayered amphibolite and Webb Canyon Gneiss. Following Reed (1973), we include this contact zone within the Moose Basin gneiss on Figure 2.

Figure 3.

Photos of field relations in the Moose Basin and Layered gneisses. Moose Basin gneiss: (A) kyanite-bearing pelitic gneiss with leucosome interpreted as the result of partial melting during granulite-facies metamorphism; (B) granulite-facies assemblages surrounded by amphibolite in mafic gneiss; (C) partially melted mafic gneiss that contains garnet-bearing leucosomes. Layered Gneiss: (D) folded paragneiss with layers of leucosome.

Figure 3.

Photos of field relations in the Moose Basin and Layered gneisses. Moose Basin gneiss: (A) kyanite-bearing pelitic gneiss with leucosome interpreted as the result of partial melting during granulite-facies metamorphism; (B) granulite-facies assemblages surrounded by amphibolite in mafic gneiss; (C) partially melted mafic gneiss that contains garnet-bearing leucosomes. Layered Gneiss: (D) folded paragneiss with layers of leucosome.

The Layered Gneiss consists of pelitic and psammitic paragneisses that are interlayered with quartzofeldspathic orthogneiss, amphibolite, metagabbro, and metaperidotite (Reed, 1963; Miller et al., 1986); pelitic gneisses are rare. Locally the paragneiss is migmatitic (Fig. 3D).

Intrusive into the Moose Basin and Layered gneisses is foliated leucogranitic Webb Canyon Gneiss of Reed (1973). On the basis of major element geochemistry, Frost et al. (2016) subdivided this unit into the Webb Canyon Gneiss and the Bitch Creek gneiss. The Webb Canyon Gneiss is the major leucogranite gneiss in the area. It forms the large tabular exposures shown on Fig. 2. In contrast, the Bitch Creek gneiss occurs in small sills, dikes, and plutons that cannot be delineated on Figure 2.

Reed (1973) mapped Augen Gneiss near the western limits of Precambrian outcrop. Because Miller et al. (1986) and Frost et al. (2016) have interpreted this unit as a porphyritic portion of the Webb Canyon Gneiss, we have shown it as part of the Webb Canyon Gneiss on Figure 2. The major ferromagnesian mineral is biotite with accessory hornblende. It is moderately well-foliated and lineated, with little or no compositional layering.

Titanite-bearing amphibolites with well-developed fabric occur as generally concordant but locally cross-cutting bodies within both the Moose Basin gneiss and the Layered Gneiss. These rocks are described in more detail below under “Thermobarometry.”

METHODS

Structural measurements made in the field were plotted and interpreted using Stereonet 9 software described by Cardozo and Allmendinger (2013). A compilation of the structural measurements is available as Table S4 (footnote 1).

U-Pb zircon geochronology was undertaken using the sensitive high-resolution ion microprobe with reverse geometry (SHRIMP-RG) at Stanford University. Zircon grains were mounted in epoxy. After polishing, cathodoluminescence (CL) scanning electron microscope (SEM) images were taken for all zircon grains. Isotopic ratios and U, Th, and Pb concentrations were standardized against VP-10 zircons. The data were reduced following the methods of Williams (1998) and using the SQUID Excel macro of Ludwig (2000). Uncertainties given for individual analyses (ratios and ages) are at the 1σ level, with correction for common lead made using the 204Pb lead and the model Pb evolution curve of Stacey and Kramers (1975). Concordia plots and linear discordia regression fits were carried out using Isoplot (Ludwig, 2003) and uncertainties are reported at the 95% confidence level. U-Pb data are reported in Supplemental Tables S1 and S21.

Sm-Nd isotopic data were obtained at the University of Wyoming on a VG Sector 54 thermal ionization multiple collector mass spectrometer. Samples were dissolved in HF and HNO3, then converted to chlorides. One aliquot of the sample was spiked with 149Sm and 146Nd. Rare-earth elements (REE) were separated using conventional cation exchange chromatography, and Sm and Nd were further separated by eluting HCl through di(2-ethylhexyl) orthophosphoric acid columns. An average 143Nd/144Nd ratio of 0.511846 ± 11 (2σ) was measured for the La Jolla Nd standard, normalized to 146Nd/144Nd = 0.7219. Uncertainties in Nd isotopic ratio measurements are ± 0.00001 (2 standard errors). Uncertainties on Sm and Nd concentrations are ± 2% of the measured values. Sm-Nd isotopic data are given in Table 2.

TABLE 2.

Sm-Nd ISOTOPIC DATA FOR LAYERED AND MOOSE BASIN GNEISSES OF THE NORTHERN TETON RANGE

Mineral compositional data were used in conjunction with whole-rock chemical composition data for thermobarometric analysis and for construction of pseudosections. Whole-rock geochemical analyses used for pseudosection calculations were obtained either from Frost et al. (2018) or by X-ray fluorescence (XRF) using a Panalytical Axios 4 KW wavelength dispersive XRF instrument at the University of Wyoming. Whole-rock compositional data used in calculation of pseudosections are given in Table 3.

TABLE 3.

WHOLE-ROCK GEOCHEMICAL DATA FOR PELITES FROM MOOSE BASIN GNEISS ANALYZED WITH PSEUDOSECTIONS

Electron probe microanalysis of rutile, garnet, pyroxene, and hornblende were conducted at the University of Wyoming using the JEOL 8900 microprobe (Tables 4A4C).

TABLE 4A.

GARNET ANALYSES, NORTHERN TETON PELITES

TABLE 4B.

RUTILE IN MOOSE BASIN PELITES

TABLE 4C.

PYROXENE ANALYSES FROM MOOSE BASIN GNEISS MAFIC GRANULITES

Pseudosections were constructed using Perple_X_6.7.7 (Connolly 2009) and fluid equation of state for H2O and thermodynamic data set from Holland and Powell (1991, 1998, 2011). In all samples, it was assumed that there was a saturated component (SiO2). The thermodynamic components selected were Na2O, MgO, Al2O3, K2O, CaO, TiO2, FeOt and H2O. Solution models employed in the calculation of both pseudosections include: Gt(W), Chl(W), St(W), Crd(W), Ctd(W), Bi(W), Mica(W), Ilm(W), and melt(W) (White et al., 2014) and feldspar (Fuhrman and Lindsley, 1988). Assumptions for fluid composition are given in the text for each pseudosection.

Mineral abbreviations used in the text are after Whitney and Evans (2010). Mineral abbreviations used in the pseudosection plots are those generated by the Perplex program, correlate with the solution models listed above, and adhere to the following conventions. If an abbreviation begins with an upper case letter, a solution model was used for that mineral in generating the pseudosection. If an abbreviation begins with a lower case letter, a pure endmember was used in generating the pseudosection. Mineral abbreviations beginning with upper case letters and shown in italics are subsets of the solution model used in the calculation of the pseudosection. For example, Pl refers to plagioclase and Afs refers to alkali feldspar produced when the feldspar (ternary feldspars) solution model was employed.

Locations and brief descriptions for all samples referenced in this contribution are given in Table S3 (footnote 1).

RESULTS

Structure

The older gneisses in the northern Teton Range have been affected by three phases of folding. The first, F1, is manifested by rarely preserved refolded isoclinal fold hinges within the pelitic gneiss in the Moose Basin gneiss and in the Layered Gneiss (Fig. 4). The second, F2, produced decameter-scale isoclinal folds in the Moose Basin gneiss, the limbs of which define the major foliation in the region (A in Fig. 4). In Moose Basin the mafic gneiss forms the cores of the F2 antiforms. Adjacent to the pelitic gneiss the mafic gneiss is amphibolite, locally with garnet. Distal from the contact with the pelitic gneiss the mafic gneiss contains granulite facies assemblages. The contact between the Moose Basin gneiss and the leucogranites is conformable to the limbs of the F2 folds (B in Fig. 4).

Figure 4.

Block diagram showing relations among the major units in the northern Teton Range. Key features: (A) The Moose Basin gneiss consists of pelitic schist and mafic gneiss with the mafic gneiss forming the cores of the F2 antiforms. The margins of the mafic gneiss bodies are always foliated amphibolites but the core of the antiform may contain granulite assemblages. The F1 folds occur only as isolated fold axes in the pelitic schist and in the Layered Gneiss. Major foliation in the area is parallel to the limbs of the F2 folds. (B) Contacts between the Moose Basin gneiss, Webb Canyon Gneiss, and Layered Gneiss lie parallel to the major foliation in the area. (C) Lineation in the Webb Canyon Gneiss is parallel to the F2 fold axes in the Moose Basin gneiss. (D) In most areas, the Webb Canyon Gneiss shows no folding; F2 folds are found only locally. Late amphibolite dikes have locally been boudinaged. All units have been folded into large F3 folds.

Figure 4.

Block diagram showing relations among the major units in the northern Teton Range. Key features: (A) The Moose Basin gneiss consists of pelitic schist and mafic gneiss with the mafic gneiss forming the cores of the F2 antiforms. The margins of the mafic gneiss bodies are always foliated amphibolites but the core of the antiform may contain granulite assemblages. The F1 folds occur only as isolated fold axes in the pelitic schist and in the Layered Gneiss. Major foliation in the area is parallel to the limbs of the F2 folds. (B) Contacts between the Moose Basin gneiss, Webb Canyon Gneiss, and Layered Gneiss lie parallel to the major foliation in the area. (C) Lineation in the Webb Canyon Gneiss is parallel to the F2 fold axes in the Moose Basin gneiss. (D) In most areas, the Webb Canyon Gneiss shows no folding; F2 folds are found only locally. Late amphibolite dikes have locally been boudinaged. All units have been folded into large F3 folds.

The Webb Canyon and Bitch Creek gneisses do not contain F1 elements, but they have a lineation that is parallel to the trend and plunge of the F2 fold axes of the Webb Canyon Gneiss (C in Fig. 4); F2 fold hinges are seen only locally in the leucogranites (D in Fig. 4), and F2 folds occur in some of the migmatites formed during peak metamorphism. The leucogranites are intruded by amphibolite dikes that locally show extensional sense-of-shear (Frost et al., 2016).

The F3 event is seen as a broad open fold that exhibits no axial planar schistosity and involves all metamorphic rocks in the area (Fig. 4). The Moose Basin gneiss occurs in the core of a major F3 synform and hence lies structurally above the Layered Gneiss, although the direct contact is nowhere preserved. Poles to the F2 foliations measured in a small region around Moose Basin show a strong maximum in the southeast quadrant (Fig. 5A). Lineations in the same area lie ∼90° away in the northeast quadrant (Fig. 5B). The poles of the foliation for all the gneisses across the northern Teton Range define a fold, the axis of which trends 20° and plunges 24° (Fig. 5C). We interpret this as the axial orientation of the F3 folding. Lineations across the northern part of the range form a strong maximum that lies in the northeast quadrant, with the same orientation as the axial plunge of the F3 fold. The fact that the lineations in the Moose Basin area (Fig. 5B) have the same trend as all the lineations across the northern part of the range (Fig. 5D) indicates that the F2 and F3 folds were nearly coaxial and suggests that they formed in the same deformation event.

Figure 5.

Structural data for the northern Teton Range. (A) Poles to the F2 foliations in the Moose Basin area. (B) Lineations (mostly L2 lineations) in the Moose Basin area. (C) Poles to foliations across the whole northern Teton Range. Cylindrical fit to the foliations defines the F3 fold, which has a trend of 20° and a plunge of 24°. (D) Lineations (mostly L2) from across the northern Teton Range. The lineations have a strong maximum near the axial trend of the F3 fold axes. Diagrams were calculated from program Stereonet (Cardozo and Allmendinger, 2013). Contour intervals are 2σ above standard deviation (Kamb, 1959).

Figure 5.

Structural data for the northern Teton Range. (A) Poles to the F2 foliations in the Moose Basin area. (B) Lineations (mostly L2 lineations) in the Moose Basin area. (C) Poles to foliations across the whole northern Teton Range. Cylindrical fit to the foliations defines the F3 fold, which has a trend of 20° and a plunge of 24°. (D) Lineations (mostly L2) from across the northern Teton Range. The lineations have a strong maximum near the axial trend of the F3 fold axes. Diagrams were calculated from program Stereonet (Cardozo and Allmendinger, 2013). Contour intervals are 2σ above standard deviation (Kamb, 1959).

GEOCHRONOLOGY

Two leucosome rocks contain zircon that formed during crystallization of the partial melts, and U-Pb dating of zircon from these rocks constrains the timing of metamorphism in the northern Teton Range (Tables S1 and S2 [footnote 1]). One is a garnet-bearing leucosome that formed from its host mafic granulite in the Moose Basin gneiss (06T16, Fig. 3B). The zircon grains in this sample range in shape from rounded to euhedral and contain inherited cores that are resorbed, zoned, and relatively bright in cathodoluminescence (Fig. 6A). The second sample is a leucosome from partially melted Layered Gneiss (06T9, Fig. 3C). The zircon grains from this sample are euhedral, show primarily oscillatory growth zoning under cathodoluminescence, and lack recognizable older cores (Fig. 6B).

Figure 6.

Cathodoluminescence image of some of the zircon grains from 06T16 and 06T9. Grains from 06T16 show bright low-U cores that are considerably older than the margins. Grains from 06T9 lack recognizable older cores and yield concordant ages from all positions in the grains. Ages are 204Pb corrected 207Pb/206Pb ages with 1 σ errors.

Figure 6.

Cathodoluminescence image of some of the zircon grains from 06T16 and 06T9. Grains from 06T16 show bright low-U cores that are considerably older than the margins. Grains from 06T9 lack recognizable older cores and yield concordant ages from all positions in the grains. Ages are 204Pb corrected 207Pb/206Pb ages with 1 σ errors.

There are two age populations of zircon from 06T16 (Fig. 7). The cores of the zircon grains, some of which are relatively bright under cathodoluminescence, give ages around 3.1 Ga (Fig. 7A). We interpret this as reflecting the age of the protolith. The rims of the grains, which are generally dark under cathodoluminescence, yield a U-Pb age of 2694.8 ± 6.7 Ma (MSWD = 1.0) (Fig. 7B; one highly discordant analysis was excluded from the calculation). We interpret these data as indicating the age of peak metamorphism and the partial melting of the mafic rocks. The zircon from 06T9 yields U-Pb ages of 2684.8 ± 4.9 Ma (MSWD = 0.3), which we interpret to represent the age of zircon growth from partial melts associated with the F2 deformation event and the M2 metamorphic event (Fig. 7C). The formation of leucosome in 06T9 was coeval with the intrusion of Webb Canyon and Bitch Creek gneisses between 2686 and 2675 Ma (Frost et al., 2016).

Figure 7.

U-Pb concordia diagrams for 06T16 and 06T9. The low-uranium relict cores in 06T16 yield Mesoarchean ages (A), while the higher-uranium rims and unzoned crystals yield younger ages (B). Zircon from 06T9 lacks recognizable older cores and yields Neoarchean ages from all locations in the crystals (C). Complete analytical data are included in Tables S1 and S2 (see footnote 1). MSWD—mean square of weighted deviates.

Figure 7.

U-Pb concordia diagrams for 06T16 and 06T9. The low-uranium relict cores in 06T16 yield Mesoarchean ages (A), while the higher-uranium rims and unzoned crystals yield younger ages (B). Zircon from 06T9 lacks recognizable older cores and yields Neoarchean ages from all locations in the crystals (C). Complete analytical data are included in Tables S1 and S2 (see footnote 1). MSWD—mean square of weighted deviates.

ISOTOPE GEOCHEMISTRY

Sm-Nd isotopic data were obtained to help constrain the protoliths of the Moose Basin and Layered gneisses. We report data for 20 felsic samples (8 samples of the Moose Basin gneiss, 12 samples of the Layered Gneiss), and 4 samples of mafic rocks within the Layered Gneiss (Table 2; Fig. 8). In general, Sm and Nd contents and Sm/Nd ratios are higher in the Layered Gneiss than in the Moose Basin gneiss. These characteristics reflect the high REE contents and lack of light REE (LREE)-enrichment typical of the Layered Gneiss compared to the generally lower REE contents and LREE enrichment of the Moose Basin gneiss (Frost et al., 2018). Mafic rocks interlayered within the Layered Gneiss have variable Sm and Nd contents and Sm/Nd that reflect flat to LREE-enriched REE patterns (Frost et al., 2018). Amphibolite dike samples (07T9 and 07T27) have higher Sm/Nd and positive initial εNd, whereas gabbro samples (07T13 and T35E) have lower Sm/Nd and slightly negative initial εNd.

Figure 8.

Diagram comparing initial εNd calculated at 2685 Ma for the Moose Basin gneiss with that for the Layered Gneiss. Data including sample numbers and sources, Sm and Nd abundances, and all ratios are included in Table 2.

Figure 8.

Diagram comparing initial εNd calculated at 2685 Ma for the Moose Basin gneiss with that for the Layered Gneiss. Data including sample numbers and sources, Sm and Nd abundances, and all ratios are included in Table 2.

Both the Moose Basin and Layered gneisses are characterized by variation in initial εNd, but there are important differences between the two units. With one exception, the Moose Basin gneisses have initial εNd less than –0.8, whereas with one exception the felsic and mafic Layered Gneiss samples all have initial εNd –0.8 or greater (Fig. 8). These differences in initial Nd isotopic composition of the Moose Basin and Layered gneisses indicate that the two paragneiss units were composed of clastic detritus that was derived from different sources: those of the Moose Basin gneiss are on average older crustal sources, and those of the Layered Gneiss are more juvenile. This conclusion is also supported by the presence of >3 Ga zircon cores in the Moose Basin gneiss and the absence of these ancient cores in all zircons analyzed from the Layered Gneiss.

THERMOBAROMETRY

Rocks with low variance assemblages were selected for thermobarometric analysis. Mineral and whole-rock chemical data were evaluated in the context of petrographic observations described above to deduce the P-T-t path inferred below.

Moose Basin Gneiss Pelitic Rocks

The pelitic rocks in the Moose Basin gneiss are interlayered with leucosomes, indicating significant melting at peak metamorphic conditions (Fig. 3A). The least retrograded Moose Basin gneiss pelites preserve evidence of an early, high-pressure granulite facies metamorphic event predating or coeval with penetrative deformation and a later amphibolite facies event that is largely unaffected by penetrative deformation. Garnets with rare kyanite, sillimanite, zircon, pyrite, and rutile inclusions record the peak metamorphic conditions experienced by these rocks. While sparse biotite, plagioclase, and ilmenite after rutile occur along fractures within the garnets and along the garnet margins, none of these phases was observed to occur as inclusions within homogeneous portions of the garnets suggesting they were not stable during early garnet growth at peak metamorphic conditions (Figs. 9A and 9B).

Figure 9.

Backscattered electron and visible light micrographs illustrating textures in Moose Basin gneiss pelites and in Layered Gneiss pelites. Moose Basin gneiss images from sample 04T30: (A) Large, partially resorbed anhedral to subhedral garnets with very few inclusions, Fe and Mn-enriched margins, and cordierite (Crd) coronas. Py—pyrite. (B) Anhedral garnet (Gt) with a single, pristine rutile (Rt) inclusion that is not associated with ilmenite. (C) Minimally deformed partial replacements of cm-scale kyanite crystals by sillimanite. (D) Staurolite (St) being replaced by cordierite and sillimanite (Sil). (E) Rutile in matrix partially replaced by ilmenite (Ilm). Layered Gneiss images from sample 06T10: (F) Large subhedral garnet with abundant inclusions (quartz, plagioclase, biotite, rutile + chlorite, chlorite, zircon and monazite). Detail for red box is shown in H. (G) Typical myrmekitic intergrowth of sericite (Ser) and quartz (Qz). (H) Rutile and chlorite (Chl) inferred to have formed from oxidation of primary ilmenite. (I) Typical intergrowth of rutile and chlorite associated with matrix biotite; the rutile + chlorite is inferred to have formed from oxidation of primary ilmenite.

Figure 9.

Backscattered electron and visible light micrographs illustrating textures in Moose Basin gneiss pelites and in Layered Gneiss pelites. Moose Basin gneiss images from sample 04T30: (A) Large, partially resorbed anhedral to subhedral garnets with very few inclusions, Fe and Mn-enriched margins, and cordierite (Crd) coronas. Py—pyrite. (B) Anhedral garnet (Gt) with a single, pristine rutile (Rt) inclusion that is not associated with ilmenite. (C) Minimally deformed partial replacements of cm-scale kyanite crystals by sillimanite. (D) Staurolite (St) being replaced by cordierite and sillimanite (Sil). (E) Rutile in matrix partially replaced by ilmenite (Ilm). Layered Gneiss images from sample 06T10: (F) Large subhedral garnet with abundant inclusions (quartz, plagioclase, biotite, rutile + chlorite, chlorite, zircon and monazite). Detail for red box is shown in H. (G) Typical myrmekitic intergrowth of sericite (Ser) and quartz (Qz). (H) Rutile and chlorite (Chl) inferred to have formed from oxidation of primary ilmenite. (I) Typical intergrowth of rutile and chlorite associated with matrix biotite; the rutile + chlorite is inferred to have formed from oxidation of primary ilmenite.

The later amphibolite facies reaction assemblage (in the same rocks) is Gt–St–Crd–Bt–Sil–Pl–Ilm–Zrn–Qz. Retrograde chlorite after biotite and garnet and intergranular sericite are also present. The garnets are rimmed by cordierite, are typically irregular in shape and have likely experienced significant resorption, consistent with enrichment of Mn and Fe around current grain boundaries (Table 4A). Because of extensive resorption and generally anhedral shapes of garnets, relative timing of crystallization of sillimanite and kyanite is not constrained by spatial distribution of the aluminosilicate inclusions in garnets, but cm-sized kyanite blades in the matrix are partially replaced by sillimanite, indicating that the kyanite predates sillimanite in these samples (Fig. 9C). Staurolite is locally reacting to cordierite + sillimanite (Fig. 9D). Ilmenite occurs as mantles on rutile within the matrix (Fig. 9E), but is absent on rutile inclusions in unfractured garnet (Fig. 9B). Textural relations clearly indicate that the Gt–St–Crd–Bt–Sil–Pl–Ilm–Zrn–Qz assemblage post-dates the Gt–Ky–Sil–Zrn–Qz assemblage, and the delicate reaction textures lack evidence of significant deformation after reaction. Finally, alkali feldspar is restricted to the leucosomes and is completely absent in the pelitic rocks.

Temperatures for both assemblages are constrained by the Zr-in-rutile thermometer (Tomkins et al., 2007). Rutile occurring as isolated inclusions in garnet has Zr as high as 3800 ppm, whereas rutile in the matrix and near fractures in the garnet ranges from 260 to 400 ppm Zr (Table 4B). The Zr-in-rutile thermometer yields temperatures of 893 °C (8 kbar) to 920 °C (12 kbar) for the rutile inclusions in garnet, and temperatures of 630 °C to 690 °C (depending on the SiO2 polymorph and the pressure) for the rutile in the matrix. Because rutile loses Zr during re-equilibration at lower temperature, the maximum temperature was at least 900 °C.

A pseudosection was calculated for one sample (04T30) assuming SiO2 saturation, fixed 4% H2O by weight, and components Na2O, MgO, Al2O3, K2O, CaO, TiO2, and FeO (Fig. 10). Details for solution models used in the calculation are given in the “Methods” section. We did not include MnO in this calculation because of its low concentration in the whole-rock chemical analysis. We did not include Fe2O3 because the ilmenite is very nearly stoichiometric and pyrite is present in both the matrix and as inclusions in the garnets in this sample. We based our estimate of H2O on loss on ignition from preparation of the XRF borosilicate beads; this estimate is supported by the fact that significantly lower estimates for H2O (2.5%) yield large alkali feldspar fields near the independently estimated T and P, and significantly higher estimates for H2O (6% and H2O saturation) truncate the kyanite/sillimanite reaction curve at a melt field above 700 °C.

Figure 10.

Pseudosection for migmatitic pelite from Moose Basin gneiss (sample 04T30) showing coexisting high-pressure granulite facies assemblage (M1, 900 °C and 12 kbar) and younger amphibolite facies assemblage (M2, 690 °C and 7 kbar). The T/P conditions for the granulite facies assemblage are constrained by very high Zr-in-rutile occurring as inclusions in undamaged garnet, the presence of both sillimanite and kyanite and the absence of ilmenite as inclusions in the garnets, the general scarcity of inclusions and absence of feldspar and biotite inclusions in the garnet, the evidence for pervasive melting, and the absence of alkali feldspar in the samples. The T/P conditions for the amphibolite facies assemblage are constrained by the St–Crd–Sil–Rt/Ilm–Zrn–Qz assemblage and are consistent with the temperature independently predicted by the much lower Zr-in-rutile occurring in the matrix and near fractures in the garnets. Mineral abbreviations used in the pseudosection plots are those generated by the Perplex program (see text for details).

Figure 10.

Pseudosection for migmatitic pelite from Moose Basin gneiss (sample 04T30) showing coexisting high-pressure granulite facies assemblage (M1, 900 °C and 12 kbar) and younger amphibolite facies assemblage (M2, 690 °C and 7 kbar). The T/P conditions for the granulite facies assemblage are constrained by very high Zr-in-rutile occurring as inclusions in undamaged garnet, the presence of both sillimanite and kyanite and the absence of ilmenite as inclusions in the garnets, the general scarcity of inclusions and absence of feldspar and biotite inclusions in the garnet, the evidence for pervasive melting, and the absence of alkali feldspar in the samples. The T/P conditions for the amphibolite facies assemblage are constrained by the St–Crd–Sil–Rt/Ilm–Zrn–Qz assemblage and are consistent with the temperature independently predicted by the much lower Zr-in-rutile occurring in the matrix and near fractures in the garnets. Mineral abbreviations used in the pseudosection plots are those generated by the Perplex program (see text for details).

The 900 °C Zr-in-rutile temperature intersects the kyanite/sillimanite phase boundary at 12 kbar, constraining the peak granulite facies conditions. These results are independent of assumptions about H2O, MnO and Fe2O3 abundances or activities. These conditions are well above the wet melting curve for biotite, consistent with the presence of significant amounts of leucocratic material in the rocks. Biotite, plagioclase, and alkali feldspar are not stable in this sample at these conditions, consistent with their absence as inclusions in the garnet. Interestingly, garnet core compositions (Table 4A) combined with 900 °C estimate for temperature also yield a minimum GRAIL (garnet–rutile–Al2SiO5 polymorph–ilmenite–quartz) pressure of equilibration for these rocks of 12 kbar (Bohlen et al., 1983). The Gt–St–Crd–Bt–Sil–Pl–Ilm–Zrn–Qz amphibolite facies assemblage occurs at 675 °C < T < 710 °C and 6.4 kbar < P < 7.5 kbar. This result is completely independent of but consistent with the Zr-in-rutile temperature based on rutile occurring outside of garnets (and rimmed by ilmenite) in this sample.

Moose Basin Gneiss Mafic Granulites

Like the pelites, the least retrograded Moose Basin gneiss mafic granulites preserve evidence of an early, granulite facies metamorphic event and a later amphibolite facies event. Mafic rocks in the cores of F2 folds preserve the granulite facies assemblage Opx–Cpx–Pl–Ilm–(Gt)±Qz, and garnet-bearing leucosomes from partial melting of these rocks occur locally (Fig 3B, C). The presence of garnet in the pre-amphibolite facies mafic granulites is inferred based on occurrence of likely pseudomorphs of cummingtonite/gedrite/plagioclase after garnet in the rocks and on the presence of garnet in the associated leucosomes; no garnet has been identified in the mafic granulites. The later amphibolite-facies assemblage (in the same rocks) is Hbl–Pl–Ilm–Qz, where hornblende occurs exclusively as mantles on ilmenite and relict orthopyroxene and clinopyroxene formed during the earlier granulite facies metamorphism. These rocks are generally massive with only weakly developed fabric. The absence of surviving garnet and the strong dependence of P, T estimates for these rocks on PH2O, however, limits their value in establishing quantitative estimates of metamorphic conditions for either the granulite facies or the amphibolite facies metamorphic events.

The two-pyroxene mafic granulite assemblage does provide temperature constraints for these rocks but not pressure constraints. Both pyroxenes are moderately exsolved and have outer rims depleted in exsolution lamellae, indicating that some re-equilibration of peak compositions has occurred. An estimate of maximum T can be obtained using cation exchange thermometry for the two-pyroxene assemblage, which is not dependent on PH2O. Using integrated compositions from cores of orthopyroxene and clinopyroxene grains (Table 4C), the QUILF program of Andersen et al. (1993) yields temperatures in the range of 816 °C to 832 °C. The presence of garnet-bearing leucosomes in the mafic granulites suggests temperatures closer to 900 °C (Fig. 3C; Pattison et al., 2003). Without preserved garnet these rocks provide no constraints on pressure of equilibration. However, they yield temperatures comparable to the more reliable pelite-based temperatures described above, and their granulite-facies assemblage and close proximity to those high-pressure pelites suggests that they probably equilibrated at pressures comparable to those above.

Layered Gneiss Migmatitic Pelites

The Layered Gneiss includes interlayered quartzofeldspathic, migmatitic, and amphibolitic rocks. The most common assemblage in the Layered Gneiss is Qz–Bt–Pl±Afs±Gt, and rocks with assemblages suitable for thermobarometry are rare. We did not recognize high-pressure granulite facies assemblages in the Layered Gneiss, but rare migmatitic pelitic samples yield amphibolite facies metamorphic conditions (Fig. 3D). Sample 06T10 is a pelitic enclave in a leucosome (06T9) from which the zircons yielded an age of 2685 Ma (Fig. 7). This rock contains Gt–St–Bt–Sil–Pl–Rt–Zrn–Qz. The biotite is partially replaced by chlorite and the plagioclase is partially replaced by sericite. Plagioclase laths >5 cm in length also occur, suggesting the presence of melt zones. The garnets are subhedral, fractured, and contain abundant inclusions of quartz, plagioclase, biotite, rutile, zircon, and rare, minute grains of ilmenite (Fig. 9F). Patches of sericite and quartz with myrmekitic texture are common (Fig.9G), indicating significant hydrothermal alteration and possibly potassium metasomatism. The rutile invariably occurs intergrown with chlorite and Zr concentration in rutile occurring within the matrix and as inclusions in garnets is generally below detection limits (Figs. 9H and 9I). For these reasons, we interpret the rutile as a retrograde product forming after ilmenite, and we believe the principal Ti phase at peak metamorphic conditions was ilmenite.

Because of evidence for late re-equilibration of mineral compositions, we restrict our analysis of thermobarometric conditions of equilibration for this sample to inferences based on assemblage alone (and not on individual mineral compositions). Garnet and staurolite are texturally early, suggesting that 06T10 records a melting event in which the pelitic rocks contained staurolite. We calculated a pseudosection for this sample assuming 4% H2O, SiO2 saturation, and no Fe2O3 (Fig.11). This pseudosection requires the Gt–St–Bt–Sil–Pl–Ilm–Zrn–Qz assemblage to coexist along a reaction curve from 665 °C, 6 kbar to 690 °C and 7.3 kbar; the absence of kyanite and primary rutile restricts this assemblage to P < 7.3 kbar (Fig. 11). Allowing some Fe2O3 would lower the activity of ilmenite and extend its stability to higher pressure. Based on this calculation, coexistence of this assemblage and melt restricts the assemblage to P > ∼7.2 kbar. The assemblage Gt–St–Sil–Pl–Qz has been shown experimentally to coexist with a water-saturated melt in the limited P-T range 665 °C < T < 680 °C and 6.6 kbar < P < 7.5 kbar (shaded area in Fig. 11; Garcia-Casco et al., 2003). The results of the pseudosection calculation are within 10 °C and 100 bars of the experimental results for the same assemblage. These relations suggest that the melting of 06T10 occurred at ∼680–700 °C and 6.5–7.5 kbar.

Figure 11.

Pseudosection for a migmatitic pelitic sample from the Layered Gneiss (sample 06T10) showing the amphibolite facies assemblage (M2, ∼7.3 kbar and ∼690 °C). This sample includes no evidence for an earlier, granulite-facies metamorphic event. The T/P conditions for the amphibolite facies assemblage are constrained by the assemblage Gt–Bt–Pl–St–Sil–(Ilm)–Zrn–Qz and by the petrographically indicated coexistence of melt with that assemblage (darkly shaded triangle, Garcia-Casco et al., 2003). Ilm is inferred to have been a part of the assemblage before hydration and retrograde processes oxidized it to produce Rt based on the fact that Rt is everywhere intimately intergrown with chlorite, the Zr concentrations in the Rt are everywhere near or below detection limits in this sample, and there is no stability field with coexisting Rt + Sil + St for this sample. Mineral abbreviations used in the pseudosection plots are those generated by the Perplex program (see text for details). WM—white mica.

Figure 11.

Pseudosection for a migmatitic pelitic sample from the Layered Gneiss (sample 06T10) showing the amphibolite facies assemblage (M2, ∼7.3 kbar and ∼690 °C). This sample includes no evidence for an earlier, granulite-facies metamorphic event. The T/P conditions for the amphibolite facies assemblage are constrained by the assemblage Gt–Bt–Pl–St–Sil–(Ilm)–Zrn–Qz and by the petrographically indicated coexistence of melt with that assemblage (darkly shaded triangle, Garcia-Casco et al., 2003). Ilm is inferred to have been a part of the assemblage before hydration and retrograde processes oxidized it to produce Rt based on the fact that Rt is everywhere intimately intergrown with chlorite, the Zr concentrations in the Rt are everywhere near or below detection limits in this sample, and there is no stability field with coexisting Rt + Sil + St for this sample. Mineral abbreviations used in the pseudosection plots are those generated by the Perplex program (see text for details). WM—white mica.

Garnet Amphibolites, Moose Basin Gneiss, and Layered Gneiss

Two compositionally distinct groups of amphibolites occur in the northern Teton Range with the assemblage Hbl–Pl–Bt–Ilm–Qz±Gt±Ttn±Cum. One set has high Fe/(Fe+Mg) whole-rock and mineral compositions, comparatively sodic plagioclase, and garnet is present in virtually all of these samples. The second set has lower Fe/(Fe+Mg) ratios and more calcic plagioclase and only a few samples in this category contain garnet (Fitz-Gerald, 2008). Both varieties occur in both the Moose Basin gneiss and in the Layered Gneiss. The amphibolites form generally concordant but locally cross-cutting bodies within the host gneisses and leucogranites; they all have well-developed fabrics. Fitz-Gerald evaluated pressures and temperatures for these rocks using the barometer of Kohn and Spear (1990) and the thermometer of Dale et al. (2000). Mineral chemistry appropriate for this barometer occurred only in the lower Fe/(Fe+Mg) amphibolites. Two garnet-bearing samples within the prescribed composition range for the barometer yield 540 °C to 600 °C and 4.0 ± 1 kbar and there is no systematic difference between amphibolites from the Layered Gneiss and the Moose Basin gneiss (Fitz-Gerald, 2008). These results are consistent with the presence of titanite in the assemblage and the absence of garnet in all but the most iron-rich amphibolites in the region.

DISCUSSION

Key Features of the Northern Teton Range

Based on the information presented above, the following observations must be taken into account in a tectonic interpretation of the Archean rocks of the northern Teton Range.

  • 1. Two contrasting gneiss units are present in the northern Teton Range:

    • a. Moose Basin gneiss, composed of metapelitic and mafic rocks. These have reached granulite facies and have partially melted. Peak pressures and temperatures are in excess of 12 kbar and 900 °C. The pelitic rocks contain zircon cores ca. 3.1 Ga and their Nd isotopic compositions are consistent with derivation from older continental crustal sources. The Moose Basin gneiss preserves evidence of both F1 and F2 deformation.

    • b. The Layered Gneiss, composed of quartzofeldspathic paragneiss, interlayered amphibolite and areas of peridotite and gabbro. These quartzofeldspathic gneisses are migmatitic but apparently did not exceed amphibolite facies conditions of ∼7.5 kbar and 700 °C. The Layered Gneiss has relatively radiogenic Nd isotopic compositions consistent with derivation from juvenile sources. The Layered Gneiss preserves evidence of both F1 and F2 deformation.

    • c. Amphibolites, some of which contain garnet, occur in both the Moose Basin gneiss and in the Layered Gneiss. Amphibolite dikes are generally concordant but locally cut the leucogranites and migmatites, suggesting that they were emplaced either during or after the emplacement of the leucogranites. Titanite-bearing garnet amphibolites in both gneisses record pressures and temperatures of 4 ± 1 kbar and 540 °C to 600 °C.

  • 2. These units are intruded by the Webb Canyon and Bitch Creek leucogranites. These leucogranites did not experience the earliest folding event, F1, which is preserved in both the Moose Basin gneiss and Layered Gneiss.

  • 3. Peak granulite facies metamorphism is dated at 2695 ± 7 Ma by zircon growth in leucosome associated with mafic granulite in the Moose Basin gneiss. The date of partial melting of the Layered Gneiss is 2685 ± 5 Ma, within the range of crystallization ages for the Webb Canyon and Bitch Creek leucogranites (2675–2686 Ma).

P-T-t Path of the Northern Teton Range

The thermobarometry shows that the rocks from the northern Teton Range record three distinct thermal regimes (Fig. 12). The Moose Basin gneiss records a high-pressure granulite-facies event of T ∼ 900 °C and P > 12 kbar (M1) with a cooling trend to 690 °C and P ∼7 kbar (M2). The Layered Gneiss records peak conditions of only amphibolite-facies, with a pressure of ∼7.3 kbar and temperature of ∼690 °C. Significant melting accompanied both metamorphic events, and some unknown amount of melt was probably lost from the affected units. Garnet amphibolites that occur within the Moose Basin gneiss and the Layered Gneiss record pressures lower than 5 kbar and temperatures less than 600 °C (M3, significantly lower grade conditions than M2).

Figure 12.

P-T-t path inferred for the northern Teton Range.

Figure 12.

P-T-t path inferred for the northern Teton Range.

From the P-T relations and the geochronology we infer that the Moose Basin gneiss underwent high-pressure granulite metamorphism and melting at 2694.8 ± 6.7 Ma; these rocks were then tectonically juxtaposed with the Layered Gneiss at conditions near those of Barrovian metamorphism at 2684.8 ± 4.9 Ma (Fig. 12). A final lower grade amphibolite facies event is recorded in both the Moose Basin gneiss and the Layered Gneiss garnet amphibolites at 540 °C < T < 600 °C and P 4 ± 1 kbar (M3). This event also produced titanite in amphibolite and gabbro throughout the Teton Range. This event is interpreted as the time of accretion of crustal blocks, including the northern Teton Range, to the Wyoming craton (Frost et al., 2018) and is unrelated to the 2.70–2.68 Ga collisional orogeny described in this paper. This event may also be responsible for some of the extensive retrogression of the higher grade assemblages observed in the Layered Gneiss and in the Moose Basin gneiss.

These observations require a clockwise P-T path for the Teton Range (Fig. 12). In Moose Basin gneiss pelites, both kyanite and sillimanite occur as inclusions in garnets, but sillimanite clearly replaced kyanite as metamorphism and deformation progressed. Rutile inclusions in garnet record temperatures of 900 °C and garnet-bearing leucosomes formed in mafic granulites (M1). Approximately 10 m.y. later, amphibolite facies assemblages replace the granulite facies assemblages in Moose Basin gneiss pelites and mafic granulites, and rutile everywhere except as inclusions in garnet lost most of its Zr in response to falling temperatures. Migmatitic Layered Gneiss lacks evidence of the early granulite facies metamorphic event (M1), but clearly records evidence of the later amphibolite facies event (M2). This clockwise P-T path is likely the product of rapid burial by thrusting, heating, uplift, and exhumation, all characteristic of continent-continent collisions in Phanerozoic time (England and Thompson, 1984; Thompson and England, 1984).

The boundary between the Moose Basin and Layered Gneiss was intruded by leucogranites between 2674 ± 5 Ma and 2686 ± 5 Ma; the dates for the Webb Canyon and Bitch Creek leucogranites are indistinguishable (Fig. 13; Frost et al., 2016). These dates indicate that the high-pressure metamorphism occurred around 10 m.y. before tectonic assembly, although within error it could have been nearly coeval. Extensive emplacement of leucogranites could have followed the 2685 Ma tectonic assembly for up to 10 m.y. The ages of the leucogranites fall within error of each other, but because some leucogranites show F2 folding and some do not, Frost et al. (2016) proposed that the emplacement must have occurred in at least two stages.

Figure 13.

Summary of geochronology of the northern Teton Range. The age of high-pressure metamorphism (M1) is interpreted from U-Pb zircon ages in the leucosome of partially melted mafic granulite (sample 06T16). The timing of tectonic assembly and M2 amphibolite grade metamorphism is provided by the age of leucosome in partially melted Layered Gneiss (06T9). Leucogranitic gneisses were intruded coincident with tectonic assembly; samples of two distinct leucogranitic suites, the Bitch Creek gneiss and the Webb Canyon Gneiss, are described in detail by Frost et al. (2016) and are arranged in order of their location from east to west.

Figure 13.

Summary of geochronology of the northern Teton Range. The age of high-pressure metamorphism (M1) is interpreted from U-Pb zircon ages in the leucosome of partially melted mafic granulite (sample 06T16). The timing of tectonic assembly and M2 amphibolite grade metamorphism is provided by the age of leucosome in partially melted Layered Gneiss (06T9). Leucogranitic gneisses were intruded coincident with tectonic assembly; samples of two distinct leucogranitic suites, the Bitch Creek gneiss and the Webb Canyon Gneiss, are described in detail by Frost et al. (2016) and are arranged in order of their location from east to west.

The Webb Canyon and Bitch Creek leucogranites, although both trondhjemitic, have distinct geochemical characteristics. The Webb Canyon Gneiss is strongly ferroan, comparatively low in alumina, and is characterized by high Zr and Y, low Sr, and high REE contents that define “seagull” shaped patterns. The Bitch Creek gneiss is lower in Zr, Y, and REE and alumina and Sr are higher than in the Webb Canyon Gneiss. Frost et al. (2016) concluded that these differences reflect two different petrogenetic processes, both of which are typical of collisional environments. The geochemical characteristics of the Bitch Creek gneiss are consistent with water-excess melting in a collision-related overthrust, where a relatively cool, hydrous lower plate releases water into a hotter upper plate. The Webb Canyon magmas formed by dehydration melting caused when dramatically thickened crust undergoes gravitational collapse and tectonic extension. Because the slope of the dehydration reaction is positive, melts formed by this process can migrate to shallower levels without intersecting their solidus. The sheet-like form and F2 fabrics of the Webb Canyon Gneiss are consistent with layer-parallel magma migration during orogenic collapse. Both water-excess and dehydration melting have been called upon to explain leucogranite of the Himalaya (Le Fort et al., 1987; Reichardt and Weinberg, 2012; Searle et al., 2009; Visona’ et al., 2012).

Similarity between Northern Tetons and Phanerozoic Continental Collisions

Five features in the gneisses of the northern Teton Range are characteristic of modern continent-continent collisional orogens.

  • (1) High-pressure granulite-facies metamorphism is followed by clockwise cooling.

  • (2) High-pressure metamorphism precedes the tectonic assembly of the gneisses by ∼10 m.y. This is typical of high-pressure rocks from the Himalaya and Alps, where rocks subducted to great depths are exhumed as the subduction zone migrates prior to the final collision (Rubatto et al., 1999; Babist et al., 2006).

  • (3) Leucogranites are emplaced in an extensional environment immediately following tectonic assembly (Frost et al., 2016). Such a time scale is typical of the Himalayan orogeny (Searle et al., 2003). Our results suggest that leucogranite emplacement accompanied tectonic assembly at 2685 Ma, and continued for as long as 10 m.y.

  • (4) Initial εNd of the metasedimentary rocks in Moose Basin gneiss and those in the Layered Gneiss are distinct, which suggests the two sedimentary packages were derived from different sources and were deposited in separate basins. The Nd isotopic evidence thus supports the interpretation that the Moose Basin gneiss represents detritus derived from a different, more ancient continental block than the Layered Gneiss, which is derived from relatively juvenile crustal sources, and that these were juxtaposed during collision at high-pressure granulite facies conditions. Ultramafic rocks within the Moose Basin and Layered gneisses may represent remnants of oceanic crust caught up with the two sedimentary packages during collision.

  • (5) The association of metaperidotite, metagabbro, and juvenile metasedimentary rocks in the Layered Gneiss is an assemblage that is found in Phanerozoic mountain belts where an ocean basin has been consumed during continental collision. An example is the Tectonic Accretion Channel in the central Alps (Engi et al., 2001).

Archean High-P Granulites

Korenaga (2013) observes that for high-pressure granulites to be observed at the Earth’s surface, two processes must occur. First, crustal rocks must be taken down to depth where they are metamorphosed under high pressure and temperatures conditions. Second, they must then be exhumed rather rapidly in order to preserve peak metamorphic conditions. One might predict that in a hotter Archean Earth, hotter slabs would be weaker and less likely to subduct. Korenaga’s analysis of the strength of Archean slabs suggests that to the contrary, they are likely to be stronger. A hotter mantle early in Earth history is likely to have experienced a greater degree of partial melting and partial melting also may have taken place at greater depths. Because water in the source mantle is transferred almost completely into the melt phase, a thicker dehydrated lithosphere is formed, one which is stiffer and stronger than hydrated mantle. As these thick plates become more dense with age they will reach neutral buoyancy, and then subduct.

If Archean lithosphere is strong enough to subduct, why are Archean high-pressure granulites so rare? Korenaga speculates that the paucity of high-pressure rocks early in Earth’s history may reflect some difficulty in the exhumation stage. Thermomechanical numerical models indicate that the ultimate driver is the ratio of the intrinsic buoyancy of the subducted continental crust to side traction forces in the conduit (Butler et al., 2014). Models by Husson et al. (2009) show that slab rollback within subduction zones induces conditions where material in the subduction channel flows upwards. Moreover, when rollback is associated with a decrease in slab dip, the exhumation process becomes more efficient. These studies suggest that exhumation is not expected to accompany vertical tectonic processes that have been called on to explain many Archean terrains. High-pressure granulites may instead be restricted to those areas where subduction enables exhumation rather than where sagduction or crustal overturn processes occur.

We have argued above that the Archean rocks of the northern Teton Range are best interpreted in terms of a collisional orogeny in which fine-grained sediments and ultramafic rocks are transported along a clockwise P-T path to depths of 35 km or greater, then exhumed during orogenic collapse and intrusion of leucogranites. This distinctly plate tectonic interpretation has not been universally invoked for the other Archean high-pressure granulite terrains shown on Table 1 and Figure 14. In many cases, incomplete preservation or a subsequent metamorphic overprint precludes definitive identification of tectonic process. For the oldest high-pressure rocks, both plate tectonic and pre-plate vertical tectonic interpretations have been proposed. The oldest high-pressure granulites occur within the 3.8 Ga Itsaq gneiss. Brown and Johnson (2018) estimate that the assemblage Gt–Cpx–Hbl–Qz formed at ∼12 kbar and 870 °C. Some authors have interpreted these rocks to record a collisional orogeny (Nutman et al., 2013, 2015), others propose that stagnant-lid tectonics were dominant at this time and that the Itsaq rock record reflects local, episodic subduction (e.g., Brown, 2016). Likewise, the high-pressure amphibolites of the Inyoni shear zone, Barberton granite-greenstone terrane have yielded multiple tectonic interpretations. Moyen et al. (2006) suggested that the high-grade rocks were subducted to depths of at least 35 km before being exhumed along the shear zone in a cycle that took place between 3.23 and 3.22 Ga. They interpreted synkinematic trondhjemites as the result of decompression melting during return flow. Alternatively, Van Kranendonk et al. (2014) propose that the greenstones foundered to depths of 12–15 kbar during partial convective overturn at around 3.23 Ga. They suggest that the greenstones were exhumed along a mylonite zone associated with an intrusive event at 3.11 Ga (Van Kranendonk et al., 2014). If so, then exhumation was much later than expected by analogy with modern continent-continent collisions.

Figure 14.

Diagram showing pressure conditions for Precambrian high-pressure granulites as a function of age. Plotted are localities that record pressures in excess of 10 kbar from Brown and Johnson (2018). Kola and Barberton did not reach granulite facies temperatures. High-pressure granulites are rare prior to the Neoarchean. High-pressure granulites of the Teton Range record a 20 million year cycle of subduction and exhumation, a history that is not overprinted by post-Archean metamorphism and deformation. VCM—Voronezh Crystalline Massif.

Figure 14.

Diagram showing pressure conditions for Precambrian high-pressure granulites as a function of age. Plotted are localities that record pressures in excess of 10 kbar from Brown and Johnson (2018). Kola and Barberton did not reach granulite facies temperatures. High-pressure granulites are rare prior to the Neoarchean. High-pressure granulites of the Teton Range record a 20 million year cycle of subduction and exhumation, a history that is not overprinted by post-Archean metamorphism and deformation. VCM—Voronezh Crystalline Massif.

Most Archean high-pressure granulites are Neoarchean, and one of the best studied also has disparate tectonic interpretations. The granulites of the Assynt block, Lewisian Gneiss Complex, record multiple periods of deformation and metamorphism, making the Archean high-pressure event difficult to interpret. Johnson et al. (2016) conclude that a subduction-related origin is possible in part because many of the rocks have trace element signatures characteristics of modern arcs. However, the authors propose that a process of sagduction is also possible, in which mafic and ultramafic rocks sank into the deep crust due to their greater density compared to underlying, partially molten felsic orthogneisses. The authors suggest that downward flow would be arrested by increasing stiffness of the orthogneiss residua as partial melt was extracted.

In contrast to Archean terranes that have been affected by Proterozoic metamorphism and deformation, the last such events in the Wyoming province were Neoarchean. The lack of subsequent tectonism greatly simplifies the interpretation of Archean events in the northern Teton Range, allowing the peak conditions and clockwise P-T-t path to be identified and the 20 million year cycle of subduction and exhumation to be quantified.

CONCLUSIONS

We suggest that the Archean rocks of the northern Teton Range formed by a process involving lateral, collisional orogeny. The metapelitic rocks of the Moose Basin gneiss derived from a ca. 3.1 Ga continental source were tectonically buried to depths corresponding to >12 kbar at 2695 Ma. This material was exhumed and juxtaposed with the Layered Gneiss, composed of juvenile metasediments and interlayered amphibolites, gabbros, and peridotite, at depths equivalent to ∼7 kbar at 2685 Ma. The mafic and ultramafic rocks are interpreted as the remnants of ocean crust consumed prior to collision. Shortly thereafter and for as long as 10 million years following, leucogranites formed by water excess melting and dehydration melting were intruded. This history is analogous to Cenozoic continental collisions such as the Alps (Engi et al., 2001; Babist et al., 2006; Searle et al., 2009).

A number of authors, including Voice et al. (2011), Condie and Aster (2010), and Condie et al. (2011) have identified major zircon U-Pb age peaks at 2.7–2.5 Ga, 2.0–1.7 Ga, 1.6 Ga, 1.2–1.0 Ga, and 0.7–0.5 Ga. These age peaks correlate with times of supercontinent formation. Brown and Johnson (2018) noted that ages of peak metamorphism also cluster at these times, and suggest that metamorphism records the amalgamation of continental fragments into supercontinents. Neoarchean continent assembly includes the construction of Superior, Sclavia, and Vaalbara (Bleeker, 2003). Collisional orogeny recorded in the Teton Range supports the hypothesis that plate tectonic processes were involved in the amagalmation of crustal blocks into larger continental masses by the Neoarchean.

ACKNOWLEDGMENTS

The authors acknowledge financial support provided by National Science Foundation grant EAR 0537670 to B.R. Frost, C.D. Frost, and S.M. Swapp and by University of Wyoming/Grand Teton National Park grant 49182 to B.R. Frost and C.D. Frost. We extend sincere thanks to Dave Mogk and to one anonymous reviewer for extremely thorough and very helpful reviews; the paper is greatly improved thanks to their comments and suggestions. We also thank Al McGrew for his help as the managing editor for this submission. Dr. Joe Wooden is thanked for hosting our analytical sessions on the SHRIMP-RG instruments at Stanford. This manuscript was completed while C.D. Frost was serving at the National Science Foundation.

1 Supplemental Tables. Table S1: Zircon analyses from leucosomes in partially melted mafic granulite (06T16). Table S2: Zircon analyses from leucosomes in partially melted layered gneiss (06T9) Table S3: Locations and brief descriptions for all samples referenced in this contribution. Table S4: Compilation of the structural measurements made in the field, plotted and interpreted using Stereonet 9 software described by Cardozo and Allmendinger (2013). Please visit https://doi.org/10.1130/GES01607.S1 or the full-text article on www.gsapubs.org to view the Supplemental Tables.
Science Editor: Shanaka de Silva
Guest Associate Editor: Allen J. McGrew

REFERENCES CITED

Andersen
,
D.J.
,
Lindsley
,
D.H.
, and
Davidson
,
P.M.
,
1993
,
QUILF: A PASCAL program to assess equilibria among Fe-Mg-Mn-Ti oxides, pyroxenes, olivine, and quartz
:
Computers & Geosciences
 , v.
19
, p.
1333
1350
, https://doi.org/10.1016/0098-3004(93)90033-2.
Babist
,
J.
,
Handy
,
M.R.
,
Konrad-Schmolke
,
M.
, and
Hammerschmidt
,
K.
,
2006
,
Precollisional, multistage exhumation of subducted continental crust: The Sesia Zone, western Alps
:
Tectonics
 , v.
25
,
TC6008
, https://doi.org/10.1029/2005TC001927.
Bleeker
,
W.
,
2003
,
The late Archean record: a puzzle in ca. 35 pieces
:
Lithos
 , v.
71
, p.
99
134
, https://doi.org/10.1016/j.lithos.2003.07.003.
Bohlen
,
S.R.
,
Wall
,
V.J.
, and
Boettcher
,
A.L.
,
1983
,
Experimental investigations and geological applications of equilibria in the system FeO-TiO2-Al2O3-SiO2-H2O
:
The American Mineralogist
 , v.
68
, p.
1049
1058
.
Brown
,
M.
,
2016
,
Archean metamorphism and geodynamics: from episodic and local to continuous and global subduction and mobile-lid plate tectonics
:
International Geological Congress Abstracts
 , v.
35
, abstract 3493.
Brown
,
M.
, and
Johnson
,
T.
,
2018
,
Secular change in metamorphism and the onset of global plate tectonics
:
The American Mineralogist
 , v.
103
, p.
181
196
, https://doi.org/10.2138/am-2018-6166.
Butler
,
J.P.
,
Beaumont
,
C.
, and
Jamieson
,
R.A.
,
2014
,
The Alps 2: Controls on crustal subduction and (ultra)high-pressure rock exhumation in Alpine-type orogens
:
Journal of Geophysical Research
 , v.
119
, p.
5987
6022
.
Cardozo
,
N.
, and
Allmendinger
,
R.W.
,
2013
,
Spherical projections with OSXStereonet
:
Computers & Geosciences
 , v.
51
, p.
193
205
, https://doi.org/10.1016/j.cageo.2012.07.021.
Condie
,
K.C.
,
2016
,
Crust and mantle evolution
, in
Condie
,
K.C.
,
Earth as an Evolving Planetary System
  (
third edition
):
London
,
Academic Press
, p.
147
199
.
Condie
,
K.C.
, and
Aster
,
R.D.
,
2010
,
Episodic zircon age spectra of orogenic granitoids: The supercontinent connection and continental growth
:
Precambrian Research
 , v.
180
, p.
227
236
, https://doi.org/10.1016/j.precamres.2010.03.008.
Condie
,
K.C.
,
Bickford
,
M.E.
,
Aster
,
R.C.
,
Belousova
,
E.
, and
Scholl
,
D.W.
,
2011
,
Episodic zircon ages, Hf isotopic composition, and the preservation rate of continental crust
:
Geological Society of America Bulletin
 , v.
123
, p.
951
957
, https://doi.org/10.1130/B30344.1.
Connolly
,
J.A.D.
,
2009
,
The geodynamic equation of state: What and how
:
Geochemistry, Geophysics, Geosystems
 , v.
10
,
Q10014
, https://doi.org/10.1029/2009GC002540.
Crowley
,
Q.C.
,
Key
,
R.
, and
Noble
,
S.R.
,
2015
,
High-precision U-Pb dating of complex zircon from the Lewisian Gneiss Complex of Scotland using an incremental CA-ID-TIMS approach
:
Gondwana Research
 , v.
27
, p.
1381
1391
, https://doi.org/10.1016/j.gr.2014.04.001.
Dale
,
J.
,
Holland
,
T.J.B.
, and
Powell
,
R.
,
2000
,
Hornblende-garnet-plagioclase thermobarometry; a natural assemblage calibration of the thermodynamics of hornblende
:
Contributions to Mineralogy and Petrology
 , v.
140
, p.
353
362
, https://doi.org/10.1007/s004100000187.
Dumond
,
G.
,
Goncalves
,
P.
,
Williams
,
M.L.
, and
Jercinovic
,
M.J.
,
2015
,
Monazite as a monitor of melting, garnet growth and feldspar recrystallization in continental lower crust
:
Journal of Metamorphic Geology
 , v.
33
, p.
735
762
, https://doi.org/10.1111/jmg.12150.
Dyck
,
B.
,
Reno
,
B.L.
, and
Kokfelt
,
T.F.
,
2015
,
The Majorqaq Belt: A record of Neoarchaean orogenesis during final assembly of the North Atlantic Craton, southern West Greenland
:
Lithos
 , v.
220–223
, p.
253
271
, https://doi.org/10.1016/j.lithos.2015.01.024.
Engi
,
M.
,
Berger
,
A.
, and
Roselle
,
G.T.
,
2001
,
Role of the tectonic accretion channel in collisional orogeny
:
Geology
 , v.
29
, p.
1143
1146
, https://doi.org/10.1130/0091-7613(2001)029<1143:ROTTAC>2.0.CO;2.
England
,
P.C.
, and
Thompson
,
A.B.
,
1984
,
Pressure-temperature-time paths of regional metamorphism I. Heat transfer during the evolution of regions of thickened continental crust
:
Journal of Petrology
 , v.
25
, p.
894
928
, https://doi.org/10.1093/petrology/25.4.894.
Finley-Blasi
,
L.
,
2009
,
The significance of the Webb Canyon Gneiss in the Archean rocks of the northern Teton Range, Wyoming [M.Sc. thesis]
:
University of Wyoming
,
Laramie, Wyoming
.
Fitz-Gerald
,
D.B.
,
2008
,
Evidence for an Archean Himalayan-style orogenic event in the northern Teton Range, Wyoming [M.Sc. thesis]
:
University of Wyoming
,
Laramie, Wyoming
.
Flowers
,
R.M.
,
Bowring
,
S.M.
,
Mahan
,
K.H.
,
Williams
,
M.L.
, and
Williams
,
I.S.
,
2008
,
Stabilization and reactivation of cratonic lithosphere from the lower crustal record in the western Canadian shield
:
Contributions to Mineralogy and Petrology
 , v.
156
, p.
529
549
, https://doi.org/10.1007/s00410-008-0301-5.
Fonarev
,
V.I.
,
Pilugin
,
S.M.
,
Savko
,
K.A.
, and
Novikova
,
M.A.
,
2006
,
Exsolution textures of orthopyroxene and clinopyroxene in high-grade BIF of the Voronezh Crystalline Massif: evidence of ultrahigh-temperature metamorphism
:
Journal of Metamorphic Geology
 , v.
24
, p.
135
151
, https://doi.org/10.1111/j.1525-1314.2006.00630.x.
Frost
,
B.R.
,
Frost
,
C.D.
,
Cornia
,
M.
,
Chamberlain
,
K.R.
, and
Kirkwood
,
R.
,
2006
,
The Teton-Wind River domain, a 2.68–2.67 Ga active margin in the western Wyoming Province
:
Canadian Journal of Earth Sciences
 , v.
43
, p.
1489
1510
, https://doi.org/10.1139/e06-102.
Frost
,
B.R.
,
Swapp
,
S.M.
,
Frost
,
C.D.
,
Bagdonas
,
D.A.
, and
Chamberlain
,
K.R.
,
2018
,
Neoarchean tectonic history of the Teton Range: Record of accretion against the present-day western margin of the Wyoming Province
:
Geosphere
 , v.
14
, 10.1130/GES01559.1.
Frost
,
C.D.
,
Swapp
,
S.M.
,
Frost
,
B.R.
,
Finley-Blasi
,
L.
, and
Fitz-Gerald
,
D.B.
,
2016
,
Leucogranites of the Teton Range, Wyoming: A record of Archean collisional orogeny
:
Geochimica et Cosmochimica Acta
 , v.
185
, p.
528
549
, https://doi.org/10.1016/j.gca.2015.12.015.
Fuhrman
,
M.L.
, and
Lindsley
,
D.H.
,
1988
,
Ternary-feldspar modeling and thermometry
:
The American Mineralogist
 , v.
73
, p.
201
215
.
Garcia-Casco
,
A.
,
Haissen
,
F.
,
Castro
,
A.
,
El-Hmidi
,
H.
,
Torres-Roldan
,
R.L.
, and
Millan
,
G.
,
2003
,
Synthesis of staurolite in melting experiments of natural metapelite: Consequences for the phase relations in low-temperature politic migmatites
:
Journal of Petrology
 , v.
44
, p.
1727
1757
, https://doi.org/10.1093/petrology/egg056.
Hamilton
,
W.B.
,
2011
,
Plate tectonics began in Neoproterozoic time, and plumes from deep mantle never operated
:
Lithos
 , v.
123
, p.
1
20
, https://doi.org/10.1016/j.lithos.2010.12.007.
Holland
,
T.J.B.
, and
Powell
,
R.
,
1991
,
A compensated-Redlich-Kwong (CORK) equation for volumes and fugacities of CO2 and H2O in the range 1 bar to 50 kbar and 100–1600 C
:
Contributions to Mineralogy and Petrology
 , v.
109
, p.
265
273
, https://doi.org/10.1007/BF00306484.
Holland
,
T.J.B.
, and
Powell
,
R.
,
1998
,
An internally consistent thermodynamic data set for phases of petrological interest
:
Journal of Metamorphic Geology
 , v.
16
, p.
309
343
, https://doi.org/10.1111/j.1525-1314.1998.00140.x.
Holland
,
T.J.B.
, and
Powell
,
R.
,
2011
,
An improved and extended internally consistent thermodynamic dataset for phases of petrological interest, involving a new equation of state for solids
:
Journal of Metamorphic Geology
 , v.
29
, p.
333
383
, https://doi.org/10.1111/j.1525-1314.2010.00923.x.
Husson
,
L.
,
Brun
,
J.-P.
,
Yamato
,
P.
, and
Faccenna
,
C.
,
2009
,
Episodic slab rollback fosters exhumation of HP-UHP rocks
:
Geophysical Journal International
 , v.
179
, p.
1292
1300
, https://doi.org/10.1111/j.1365-246X.2009.04372.x.
Johnson
,
T.E.
,
Brown
,
M.
,
Goodenough
,
K.M.
,
Clark
,
C.
,
Kinny
,
P.D.
, and
White
,
R.W.
,
2016
,
Subduction or sagduction? Ambiguity in constraining the origin of ultramafic-mafic bodies in the Archean crust of NW Scotland
:
Precambrian Research
 , v.
283
, p.
89
105
, https://doi.org/10.1016/j.precamres.2016.07.013.
Kamb
,
W.B.
,
1959
,
Ice petrofabric observations from Blue Glacier, Washington in relation to theory and experiment
:
Journal of Geophysical Research
 , v.
64
, p.
1891
1909
, https://doi.org/10.1029/JZ064i011p01891.
Karmakar
,
S.
, and
Schenk
,
V.
,
2015
,
Neoarchean UHT metamorphism and Paleoproterozoic UHT reworking at Uweinat in the East Sahara Ghost Craton, SW Egypt: evidence from petrology and texturally controlled in situ monazite dating
:
Journal of Petrology
 , v.
56
, p.
1703
1742
, https://doi.org/10.1093/petrology/egv051.
Kohn
,
M.J.
, and
Spear
,
F.S.
,
1990
,
Two new geobarometers for garnet amphibolites, with applications to southeastern Vermont
:
The American Mineralogist
 , v.
75
, p.
89
96
.
Korenaga
,
J.
,
2013
,
Initiation and evolution of plate tectonics on Earth: theories and observations
:
Annual Review of Earth and Planetary Sciences
 , v.
41
, p.
117
151
, https://doi.org/10.1146/annurev-earth-050212-124208.
Le Fort
,
P.
,
Cuney
,
M.
,
Deniel
,
C.
,
France-Lanord
,
C.
,
Sheppard
,
S.M.F.
,
Upreti
,
B.N.
, and
Vidal
,
P.
,
1987
,
Crustal generation of Himalayan leucogranites
:
Tectonophysics
 , v.
134
, p.
39
57
, https://doi.org/10.1016/0040-1951(87)90248-4.
Li
,
X.
,
Zhang
,
L.
,
Wei
,
C.
, and
Slabunov
,
A.I.
,
2015
,
Metamorphic PT path and zircon U-Pb dating of Archean eclogite association in Gridino complex, Belomorian province, Russia
:
Precambrian Research
 , v.
268
, p.
74
96
, https://doi.org/10.1016/j.precamres.2015.07.009.
Love
,
J.D.
,
Reed
,
J.C.
, Jr
., and
Christiansen
,
A.C.
,
1992
,
Geologic map of Grand Teton National Park, Teton County, Wyoming
:
U.S. Geological Survey Miscellaneous Investigations Series Map I-2031, scale 1:62, 500
 ,
17
p. text.
Ludwig
,
K.R.
,
2000
,
SQUID 1.02, A User’s Manual
:
Berkeley, California
,
USA, Berkeley Geochronology Center Special Publication
.
Ludwig
,
K.R.
,
2003
,
ISOPLOT: A geochronological toolkit for Microsoft EXCEL
:
Berkeley, California, USA, Berkeley Geochronology Center Special Publication No. 4
 ,
70
p.
Miller
,
S.H.
,
Hildebrandt
,
P.K.
,
Erslev
,
E.A.
, and
Reed
,
J.C.
, Jr
.,
1986
,
Metamorphic and deformation history of the gneiss complex in the northern Teton Range, Wyoming
, in
Garrison
,
P.B.
, ed.,
Geology of the Beartooth Uplift and Adjacent Basins, United States (USA): Montana Geological Society and Yellowstone Bighorn Research Association Joint Field Conference and Symposium
 , p.
91
110
.
Mints
,
M.V.
,
Belousova
,
E.A.
,
Konilov
,
A.N.
,
Natapov
,
L.M.
,
Shchipansky
,
A.A.
,
Griffin
,
W.L.
,
O’Reilly
,
S.Y.
,
Dokukina
,
K.A.
, and
Kaulina
,
T.V.
,
2010
,
Meosarchean subduction processes: 2.87 eclogites from the Kola Peninsula, Russia
:
Geology
 , v.
38
, p.
739
742
, https://doi.org/10.1130/G31219.1.
Moyen
,
J.-F.
,
Stevens
,
G.
, and
Kisters
,
A.
,
2006
,
Record of mid-Archaean subduction from metamorphism in the Barberton terrain, South Africa
:
Nature
 , v.
442
, p.
559
562
, https://doi.org/10.1038/nature04972.
Nicoli
,
G.
,
Stevens
,
G.
,
Moyen
,
J.-F.
, and
Frei
,
D.
,
2015
,
Rapid evolution from sediment to anatectic granulite in an Archean continental collision zone: The example of the Bandelierkop Formation metapelites, South Marginal Zone, Limpopo Belt, South Africa
:
Journal of Metamorphic Geology
 , v.
33
, p.
177
202
, https://doi.org/10.1111/jmg.12116.
Nutman
,
A.P.
,
Bennett
,
V.C.
,
Friend
,
C.R.L.
,
Hidaka
,
H.
,
Yi
,
K.
,
Lee
,
S.R.
, and
Kamiichi
,
T.
,
2013
,
The Itsaq Gneiss Complex of Greenland: Episodic 3900 to 3660 Ma juvenile crust formation and recycling in the 3660 to 3600 Ma Isukasian orogeny
:
American Journal of Science
 , v.
313
, p.
877
911
, https://doi.org/10.2475/09.2013.03.
Nutman
,
A.P.
,
Bennett
,
V.C.
, and
Friend
,
C.R.L.
,
2015
,
Proposal for a continent ‘Itsaqia’ amalgamated at 3.66 Ga and rifted apart from 3.53 Ga: initiation of a Wilson cycle near the start of the rock record
:
American Journal of Science
 , v.
315
, p.
509
536
, https://doi.org/10.2475/06.2015.01.
O’Neill
,
C.
,
Lenardic
,
A.
,
Weller
,
M.
,
Moresi
,
L.
,
Quenette
,
S.
, and
Zhang
,
S.
,
2016
,
A window for plate tectonics in terrestrial planet evolution?
:
Physics of the Earth and Planetary Interiors
 , v.
255
, p.
80
92
, https://doi.org/10.1016/j.pepi.2016.04.002.
Pattison
,
D.R.M.
,
Chacko
,
T.
,
Farquhar
,
J.
, and
McFarlane
,
C.R.M.
,
2003
,
Temperatures of granulite-facies metamorphism: Constraints from experimental phase equilibria and thermobarometry corrected for retrograde exchange
:
Journal of Petrology
 , v.
44
, p.
867
900
, https://doi.org/10.1093/petrology/44.5.867.
Reed
,
J.C.
, Jr
.,
1963
,
Structure of Precambrian crystalline rocks in the northern part of Grand Teton National Park, Wyoming
:
U.S. Geological Survey Professional Paper
 , v.
475-C
, p.
C1
C6
.
Reed
,
J.C.
,
1973
,
Geologic map of the Teton Range, Wyoming
:
U.S. Geological Survey Open-File Map, scale 1:62,600
 .
Reichardt
,
H.
, and
Weinberg
,
R.F.
,
2012
,
The dike swarm of the Karakoram shear zone, Ladakh, NW India: linking granite source to batholith
:
Geological Society of America Bulletin
 , v.
124
, p.
89
103
, https://doi.org/10.1130/B30394.1.
Roberts
,
N.M.W.
,
van Kranendonk
,
M.J.
,
Parman
,
S.
, and
Clift
,
P.D.
,
2015
,
Continent formation through time
, in
Roberts
,
N.M.W.
, et al.  
, eds.,
Continent Formation Through Time: Geological Society of London Special Publication 389
 , p.
1
16
, https://doi.org/10.1144/SP389.13.
Rubatto
,
D.
,
Gebauer
,
D.
, and
Conpagnoni
,
R.
,
1999
,
Dating of eclogite-facies zircons: the age of Alpine metamorphism in the Sesia-Lanzo Zone (Western Alps)
:
Earth and Planetary Science Letters
 , v.
167
, p.
141
158
, https://doi.org/10.1016/S0012-821X(99)00031-X.
Savko
,
K.A.
,
Kotov
,
A.B.
,
Sal’nikova
,
E.B.
,
Korish
E.Kh.
,
Pilyugin
,
S.M.
,
Artemenko
,
G.V.
, and
Korikovskii
,
S.P.
,
2010
,
The age of metamorphism of granulite complexes of the Voronezh Crystalline Massif: the monazite U-Pb geochronology
:
Doklady Earth Sciences
 , v.
435
, p.
1575
1580
, https://doi.org/10.1134/S1028334X10120056.
Searle
,
M.P.
,
Simpson
,
R.L.
,
Law
,
R.D.
,
Parrish
,
R.R.
, and
Waters
,
D.J.
,
2003
,
Structural geometry, metamorphic and magmatic evolution of the Everest massif, High Himalaya of Nepal—South Tibet
:
Journal of the Geological Society
 , v.
160
, p.
345
366
, https://doi.org/10.1144/0016-764902-126.
Searle
,
M.P.
,
Cottle
,
J.M.
,
Streule
,
M.J.
, and
Waters
,
D.J.
,
2009
,
Crustal melt granites and migmatites along the Himalaya: melt source, segregation, transport and granite emplacement mechanisms
:
Earth and Environmental Science Transactions of the Royal Society of Edinburgh
 , v.
100
, p.
219
233
, https://doi.org/10.1017/S175569100901617X.
Smart
,
K.A.
,
Tappe
,
S.
,
Stern
,
R.A.
,
Webb
,
S.J.
, and
Ashwal
,
L.D.
,
2016
,
Early Archaean tectonics and mantle redox recorded in Witwatersrand diamonds
:
Nature Geoscience
 , v.
9
, p.
255
259
, https://doi.org/10.1038/ngeo2628.
Stacey
,
J.S.
, and
Kramers
,
J.D.
,
1975
,
Approximation of terrestrial lead isotope evolution by a two-stage model
:
Earth and Planetary Science Letters
 , v.
26
, p.
207
221
, https://doi.org/10.1016/0012-821X(75)90088-6.
Tang
,
M.
,
Chen
,
K.
, and
Rudnick
,
R.
,
2016
,
Archean upper crust transition from mafic to felsic marks the onset of plate tectonics
:
Science
 , v.
351
, p.
372
375
, https://doi.org/10.1126/science.aad5513.
Thompson
,
A.B.
, and
England
,
P.C.
,
1984
,
Pressure-temperature time paths of regional metamorphism II: Their inference and interpretation using mineral assemblages in metamorphic rocks
:
Journal of Petrology
 , v.
25
, p.
929
955
, https://doi.org/10.1093/petrology/25.4.929.
Tomkins
,
H.S.
,
Powell
,
R.
, and
Ellis
,
D.J.
,
2007
,
The pressure dependence of the zirconium-in-rutile thermometer
:
Journal of Metamorphic Geology
 , v.
25
, p.
703
713
, https://doi.org/10.1111/j.1525-1314.2007.00724.x.
Van Kranendonk
,
M.J.
,
2010
,
Two types of Archean continental crust: Plume and plate tectonics on early Earth
:
American Journal of Science
 , v.
310
, p.
1187
1209
, https://doi.org/10.2475/10.2010.01.
Van Kranendonk
,
M.J.
,
Kroner
,
A.
,
Hoffmann
,
J.E.
,
Nagel
,
T.
, and
Anhaeusser
,
C.R.
,
2014
,
Just another drip: re-analysis of a proposed Mesoarchean suture from the Barberton Mountain Land, South Africa
:
Precambrian Research
 , v.
254
, p.
19
35
, https://doi.org/10.1016/j.precamres.2014.07.022.
Visona’
,
D.
,
Carosi
,
R.
,
Montomoli
,
C.
,
Tiepolo
,
M.
, and
Peruzzo
,
L.
,
2012
,
Miocene andalusite leucogranite in central-east Himalaya (Everest-Masang Kang area): Low-pressure melting during heating
:
Lithos
 , v.
144–145
, p.
194
208
, https://doi.org/10.1016/j.lithos.2012.04.012.
Voice
,
P.J.
,
Kowalewski
,
M.
, and
Eriksson
,
K.A.
,
2011
,
Quantifying the timing and rate of crustal evolution: global compilation of radiometrically dated detrital zircon grains
:
The Journal of Geology
 , v.
119
, p.
109
126
, https://doi.org/10.1086/658295.
White
,
R.W.
,
Powell
,
R.
,
Holland
,
T.J.B.
,
Johnson
,
T.E.
, and
Green
,
E.C.R.
,
2014
,
New mineral activity-composition relations for thermodynamic calculations in metapelitic systems
:
Journal of Metamorphic Geology
 , v.
32
, p.
261
286
, https://doi.org/10.1111/jmg.12071.
Whitney
,
D.L.
, and
Evans
,
B.W.
,
2010
,
Abbreviations for names of rock-forming minerals
:
The American Mineralogist
 , v.
95
, p.
185
187
, https://doi.org/10.2138/am.2010.3371.
Williams
,
I.S.
,
1998
,
U-Th-Pb geochronology by ion microprobe
:
Reviews in Economic Geology
 , v.
7
, p.
1
35
.
Worthington
,
L.L.
,
Miller
,
K.C.
,
Chamberlain
,
K.R.
,
Erslev
,
E.
,
Sheehan
,
A.
, et al.  
.,
2015
,
Building on the legacy of deep probe and CD-ROM: Crustal structure of the Bighorn Mountains, Wyoming and its relationship to Precambrian and Laramide tectonics
:
Geological Society of America Abstracts with Programs
 , v.
47
, no.
7
, p.
463
.
Zartman
,
R.E.
, and
Reed
,
J.C.
,
1998
,
Zircon geochronology of the Webb Canyon Gneiss and the Mount Owen Quartz Monzonite, Teton Range Wyoming: significance to dating late Archean metamorphism in the Wyoming craton
:
The Mountain Geologist
 , v.
35
, p.
71
77
.
Zirkler
,
A.
,
Johnson
,
T.E.
,
White
,
R.W.
, and
Zach
,
T.
,
2012
,
Polymetamorphism in the mainland Lewisian complex, NW Scotland—phase equilibria and geochronological constraints from the Cnoc an t’Sidhean suite
:
Journal of Metamorphic Geology
 , v.
30
, p.
865
885
, https://doi.org/10.1111/j.1525-1314.2012.01003.x.
Gold Open Access: This paper is published under the terms of the CC-BY-NC license.

Supplementary data