ABSTRACT A large portion of the Belt-Purcell Supergroup is well exposed in the vicinity of Glacier and Waterton Lakes national parks of northwestern Montana, USA, and southwestern Alberta, Canada. These strata were deposited in the northeastern part of the Mesoproterozoic Belt Basin. The dramatic rate of subsidence combined with dominantly fine-grained sediment influx produced thick units of broadly uniform lithology, which constitute the spectacular and unusually colorful mountain scenery of this region. Seemingly fairly simple at first glance, in detail these rocks exhibit a great deal of facies heterogeneity and a number of unusual attributes. This has resulted in contrasting and controversial interpretations of sedimentary features, depositional dynamics, sedimentary environments, and consequently the overall understanding of the entire basin. The Belt Basin reveals itself to be a unique setting in many respects, but ideas stemming from these rocks have implications for other strata, not just those of pre-Cambrian age, but for the entire Phanerozoic as well. The Belt Supergroup is therefore a particularly stimulating field-trip destination that challenges textbook interpretations.

Crinkle cracks are sand-filled cracks up to 5 mm wide in plan view that pinch at their ends. In cross section, they are canted and crinkled. They cut mudstone beds that underlie hummocky cross-laminated sandstone lenses. They are here described from the Piegan Group, Proterozoic Belt Supergroup, but they are widespread in Proterozoic and Phanerozoic marine and lacustrine rocks. However, they represent a distinctive, descriptive style of mudcracks, not attributed to inferred syneresis processes, although they have been commonly attributed to syneresis. In plan view, crinkle cracks closely resemble cracks formed where oscillatory waves striking viscous mud banks are transformed into fluid solitary-like waves that open surface cracks on their trailing limbs and close the cracks on their leading limbs as they pass through the viscous mud. Crinkle cracks preserved in rocks are hypothetically attributed to oscillatory waves moving sand over viscous mud. The oscillatory waves are transformed into solitary-like waves as they pass down into the mud, forming the cracks. The surface sand falls down into the cracks, preserving them. With burial, the water escapes, and the viscous mud compacts, crinkling the sand-filled cracks.

A recent 1:24,000 scale mapping project within the northern Beaverhead Mountains along the Idaho-Montana border has resulted in a reinterpretation of both the Mesoproterozoic stratigraphy and the regional structural framework. A 15-km-thick stratigraphic section of the Mesoproterozoic Lemhi subbasin was initially deformed by northeast-southwest shortening into giant northwest-striking, northeast-verging folds, probably during Cretaceous Sevier orogenesis. These initial folds were then dissected by a system of subparallel and anastomosing, oblique-slip reverse, thrust, and normal faults that generally strike northwest, but that exhibit east-west–oriented lineations, suggesting components of strike-slip displacement. Contractional faulting appears to have been followed by Eocene to Miocene extensional faulting, with many normal faults following the preexisting fabrics. Extension opened Tertiary basins along some of these faults, including the Salmon Basin along the southwestern side of the Beaverhead Range. Subparallel faults in the surrounding region appear to have a similar complex history, and all appear to be part of a major northwest-striking Cretaceous fold-and-thrust belt that was later dissected by Tertiary extension. Although the faults of the Beaverhead Mountains are significant and long-lived, they are not terrane-bounding structures separating the Belt and Lemhi sedimentary sequences. Instead, Lemhi strata extend across the range and northward to Missoula, where they grade into correlative Missoula Group strata.

The Blackbird cobalt-copper (Co-Cu) district in the Salmon River Mountains of east-central Idaho occupies the central part of the Idaho cobalt belt—a northwest-elongate, 55-km-long belt of Co-Cu occurrences, hosted in grayish siliciclastic metasedimentary strata of the Lemhi subbasin (of the Mesoproterozoic Belt-Purcell Basin). The Blackbird district contains at least eight stratabound ore zones and many discordant lodes, mostly in the upper part of the banded siltite unit of the Apple Creek Formation of Yellow Lake, which generally consists of interbedded siltite and argillite. In the Blackbird mine area, argillite beds in six stratigraphic intervals are altered to biotitite containing over 75 vol% of greenish hydrothermal biotite, which is preferentially mineralized. Past production and currently estimated resources of the Blackbird district total ~17 Mt of ore, averaging 0.74% Co, 1.4% Cu, and 1.0 ppm Au (not including downdip projections of ore zones that are open downward). A compilation of relative-age relationships and isotopic age determinations indicates that most cobalt mineralization occurred in Mesoproterozoic time, whereas most copper mineralization occurred in Cretaceous time. Mesoproterozoic cobaltite mineralization accompanied and followed dynamothermal metamorphism and bimodal plutonism during the Middle Mesoproterozoic East Kootenay orogeny (ca. 1379–1325 Ma), and also accompanied Grenvilleage (Late Mesoproterozoic) thermal metamorphism (ca. 1200–1000 Ma). Stratabound cobaltite-biotite ore zones typically contain cobaltite 1 in a matrix of biotitite ± tourmaline ± minor xenotime (ca. 1370–1320 Ma) ± minor chalcopyrite ± sparse allanite ± sparse microscopic native gold in cobaltite. Such cobaltite-biotite lodes are locally folded into tight F 2 folds with axial-planar S 2 cleavage and schistosity. Discordant replacement-style lodes of cobaltite 2 -biotite ore ± xenotime 2 (ca. 1320–1270 Ma) commonly follow S 2 fractures and fabrics. Discordant quartz-biotite and quartz-tourmaline breccias, and veins contain cobaltite 3 ± xenotime 3 (ca. 1058–990 Ma). Mesoproterozoic cobaltite deposition was followed by: (1) within-plate plutonism (530–485 Ma) and emplacement of mafic dikes (which cut cobaltite lodes but are cut by quartz-Fe-Cu-sulfide veins); (2) garnet-grade metamorphism (ca. 151–93 Ma); (3) Fe-Cu-sulfide mineralization (ca. 110–92 Ma); and (4) minor quartz ± Au-Ag ± Bi mineralization (ca. 92–83 Ma). Cretaceous Fe-Cu-sulfide vein, breccia, and replacement-style deposits contain various combinations of chalcopyrite ± pyrrhotite ± pyrite ± cobaltian arsenopyrite (not cobaltite) ± arsenopyrite ± quartz ± siderite ± monazite (ca. 144–88 Ma but mostly 110–92 Ma) ± xenotime (104–93 Ma). Highly radiogenic Pb (in these sulfides) and Sr (in siderite) indicate that these elements resided in Mesoproterozoic source rocks until they were mobilized after ca. 100 Ma. Fe-Cu-sulfide veins, breccias, and replacement deposits appear relatively undeformed and generally lack metamorphic fabrics. Composite Co-Cu-Au ore contains early cobaltite-biotite lodes, cut by Fe-Cu-sulfide veins and breccias, or overprinted by Fe-Cu-sulfide replacement-style deposits, and locally cut by quartz veinlets ± Au-Ag ± Bi minerals.

The redox state of the mid-Proterozoic oceans, lakes, and atmospheres is still debated, but it is vital for understanding the emergence and rise of macroscopic organisms and eukaryotes. The Appekunny Formation, Belt Supergroup, Montana, contains some of these early macrofossils dated between 1.47 Ga and 1.40 Ga and provides a well-preserved record of paleoenvironmental conditions. We analyzed the iron chemistry and mineralogy in samples from Glacier National Park, Montana, by pairing bulk rock magnetic techniques with textural techniques, including light microscopy, scanning electron microscopy, and synchrotron-based X-ray absorption spectroscopy. Field observations of the Appekunny Formation combined with mineralogical information allowed revised correlations of stratigraphic members across the park. However, late diagenetic and/or metasomatic fluids affected primary iron phases, as evidenced by prevalent postdepositional phases including base-metal sulfides. On the west side of the park, pyrrhotite and chlorite rims formed during burial metamorphism in at least two recrystallization events. These complex postdepositional transformations could affect bulk proxies for paleoredox. By pairing bulk and textural techniques, we show primary records of redox chemistry were preserved in early diagenetic and often recrystallized framboidal pyrite, submicron magnetite grains interpreted to be detrital in origin, and red-bed laminae interpreted to record primary detrital oxides. Based on these observations, we hypothesize that the shallow waters of the mid-Proterozoic Belt Basin were similar to those in modern marine and lacustrine waters: fully oxygenated, with detrital reactive iron fluxes that mineralized pyrite during organic diagenesis in suboxic, anoxic, and sulfidic conditions in sedimentary pore waters.

Proterozoic mafic magmatic rocks exposed along the western side of North America, or western Laurentia, from Kimberley, British Columbia, through to northwestern Wyoming have been previously divided into two large igneous provinces: the ca. 1460 Ma Moyie-Purcell and the ca. 780 Ma Gunbarrel large igneous provinces. New geochemical analysis from this study demonstrates that there are additional intraplate mafic magmatic rocks present. Distinguishable by variable normalized rare earth element patterns combined with differing slopes on a binary Ti versus V plot, there are 17 identifiable geochemical signatures in the 307 whole-rock and trace-element analyses from this study. Only seven of these signatures can be linked to the ca. 1460 Ma Moyie-Purcell large igneous province, and one signature to the 780 Ma Gunbarrel large igneous province. This study has identified two groups of intrusions with distinct geochemical signatures previously linked with the ca. 1460 Ma Moyie-Purcell large igneous province but now recognized to be separate events, a single unique geochemical signature with a U-Pb age correlative with the Moyie-Purcell large igneous province and seven other heretofore unidentified signatures interpreted to belong to additional undated events.

We investigated the crustal structure of the central Mesoproterozoic Belt Basin in northwestern Montana and northern Idaho using a crustal resistivity section derived from a transect of new short- and long-period magnetotelluric (MT) stations. Two- and three-dimensional resistivity models were generated from these data in combination with data collected previously along three parallel short-period MT profiles and from EarthScope MT stations. The models were interpreted together with coincident deep seismic-reflection data collected during the Consortium for Continental Reflection Profiling (COCORP) program. The upper-crustal portion of the resistivity model correlates well with the mapped surface geology and reveals a three-layer resistivity stratigraphy, best expressed beneath the axis of the Libby syncline. Prominent features in the resistivity models are thick conductive horizons that serve as markers in reconstructing the disrupted basin stratigraphy. The uppermost unit (up to 5 km thick), consisting of all of the Belt Supergroup above the Prichard Formation, is highly resistive (1000–10,000 Ω·m) and has relatively low seismic layer velocities. The intermediate unit (up to 7 km thick) consists of the exposed Prichard Formation and 3+ km of stratigraphy below the deepest stratigraphic exposures of the unit. The intermediate unit has low to moderate resistivity (30–200 Ω·m), relatively high seismic velocities, and high seismic reflectivity, with the latter two characteristics resulting from an abundance of thick syndepositional mafic sills. The lowest unit (5–10 km thick) is nowhere exposed but underlies the intermediate unit and has very high conductivity (4–8 Ω·m) and intermediate seismic velocities. This 17–22-km-thick three-layer stratigraphy is repeated below the Libby syncline, with a base at ~37 km depth. Seismic layer velocities indicate high mantle-like velocities below 37 km beneath the Libby syncline. The continuous high-conductivity layer in the lower repeated section is apparently displaced ~26 km to the east above a low-angle normal fault inferred to be the downdip continuation of the Eocene, east-dipping Purcell Trench detachment fault. Reversal of that and other Eocene displacements reveals a >50-km-thick crustal section at late Paleocene time. Further reversal of apparent thrust displacements of the three-layer stratigraphy along the Lewis, Pinkham, Libby, and Moyie thrusts allows construction of a restored cross section prior to the onset of Cordilleran thrusting in the Jurassic. A total of ~220 km of Jurassic–Paleocene shortening along these faults is indicated. The enhanced conductivity within the lowest (unexposed) Belt stratigraphic unit is primarily attributed to one or more horizons of laminated metallic sulfides; graphite, though not described within the Belt Supergroup, may also contribute to the enhanced conductivity of the lowest stratigraphic unit. A narrow conductive horizon observed within the Prichard Formation in the eastern part of the transect correlates with the stratigraphic position of the world-class Sullivan sedimentary exhalative massive sulfide deposit in southern British Columbia, and it may represent a distal sulfide blanket deposit broadly dispersed across the Belt Basin. By analogy, the thick conductive sub–Prichard Formation unit may represent repeated sulfide depositional events within the early rift history of the basin, potentially driven by hydrothermal fluids released during basaltic underplating of attenuated continental crust.

A succession of separate tectonic events affected Mesoproterozoic Belt Supergroup strata of NW Montana, just as in the better-displayed Coeur d'Alene Mining District of Idaho. A series of these established a tectonic zone historically known as the Lewis and Clark line, here re-identified as the Lewis and Clark tectonic zone, an apparent product of periodic reactivation of fundamental basement structures and physical constraint of a growth fault on developing folds. Six events identify a partial tectonic history of the west-central Belt Basin. The oldest produced growth faults concentrated along at least two structural lineaments, one of which, the Jocko line, substantially controlled the distribution of subsequent structures; the other, the north-trending Noxon line, is implicated in creation of metal-enriched source rock for Coeur d'Alene veins and provides a marker for right-lateral faulting within the Lewis and Clark tectonic zone. Subsequent deformation produced (1) west-northwest–trending folds, mostly confined to the Lewis and Clark tectonic zone and terminating northward against the Jocko line as the likely result of their having been compressed against this pre–Belt Basin structure; (2) north-trending regional folds, which extend southward from Canada and cross the ultimate Lewis and Clark tectonic zone; (3) foliated shear zones in the Lewis and Clark tectonic zone and associated Coeur d'Alene veins and reverse faults; (4) right-lateral, transcurrent faults, so identified by offsets of the Noxon line, north-trending regional folds, and the Montana overthrust belt and its associated foredeep basin; and, last, (5) Lewis and Clark tectonic zone normal faults and associated kink folds, which extensively reached gigantic, “megakink” proportions. These megakinks locally disrupted all prior structures, greatly confusing local structure; these need to be “unkinked,” so that structures resulting from prior tectonism may be fully recognized and correctly interpreted. Many faults of the Lewis and Clark tectonic zone trend southeasterly in its easterly part, tracking pre–Belt Basin structures separate from those associated with the Jocko line.

The bedding-plane anisotropy and structural configuration of the Mesoproterozoic Belt-Purcell Supergroup guided a narrow magmatic salient >350 km eastward from the Salmon River suture of Idaho to the foreland basin of central Montana, along a deep graben within the southern part of the Belt-Purcell Basin. The magma assimilated anatectic melt from the lower Belt-Purcell Supergroup in the western half of the graben, where the Lemhi subbasin had intersected and deepened the graben by several kilometers. The magma stepped across the stratigraphic section as it intruded eastward along the graben, spread laterally as it climbed into the overlying Paleozoic and Mesozoic strata, and eventually erupted into the foreland basin. This paper develops a model in which the magma formed a thick, east-tapering wedge beneath the Belt-Purcell carapace. The wedge elevated and tilted its lid, which failed along the trend of the graben to a terminus in the Crazy Mountains basin of the Helena structural salient, much like a tectonic-scale landslide. The carapace failed in two main phases between ca. 100 and 75 Ma. It slid ~100 km during the first failure phase, and ~40 km during the second, when the Boulder batholith and its volcanic cover filled a large pull-apart structure within the carapace. Slaty cleavage, tectonic slides that omit strata, and a nested series of hairpin-shaped allochthons characterize the failed carapace. Shear zones and nappes bound the carapace; the sinistral Lewis and Clark line bounds it on the north, and the dextral southwest Montana transverse zone bounds it on the south. The Lewis thrust fault and associated structures of the Rocky Mountain fold-and-thrust belt overprinted and displaced the magmatic salient and its carapace from ca. 74 to 59 Ma. The magma crystallized, cooled, and generated hydrothermal ore deposits from Late Cretaceous to middle Eocene time. Eocene extension overprinted the system from 53 to 39 Ma and exhumed its infrastructure in core complexes. Those exposures, together with regional structural tilt, enable reconstruction of a deep cross section of the magmatic wedge and its carapace.

The Mesoproterozoic Belt-Purcell Basin of the United States–Canadian Rocky Mountains formed in a complex intracontinental-rift system. The basin contained three main fault blocks: a northern half-graben, a central horst, and a southern graben. Each had distinct internal stratigraphy and mineralization that influenced Phanerozoic sedimentation; the northern half-graben and horst formed a platform with a condensed section, whereas the southern graben formed the subsiding Central Montana trough. They formed major crustal blocks that rotated clockwise during Cordilleran thrust displacement, with transpressional shear zones deforming their edges. The northern half-graben was deepest and filled with a structurally strong prism of quartz-rich sedimentary rocks and thick mafic sills that tapered toward the northeast from >15-km-thick near the basin-bounding fault. This strong, dense prism was driven into the foreland basin as a readymade, critically tapered tectonic wedge and was inverted into the Purcell anticlinorium. Erosion did not breech the Belt-Purcell Supergroup in this prism during thrusting. The southern graben was thinner, weaker, lacked mafic sills, and was engorged with sheets of granite during thrusting. It was internally deformed to achieve critical taper and shed thick deposits of syntectonic Belt-Purcell–clast conglomerate into the foreland basin. A palinspastic map of the basin combined with a detailed paleocontinental map that juxtaposes the northeastern corner of the Siberian craton against western North America indicates that the basin formed at the complicated junction of three continental-scale rift zones.