SHRIMP U-Pb and 40Ar/39Ar age constraints for relating plutonism and mineralization in the Boulder batholith region, Montana

Karen Lund; John N. Aleinikoff; Michael J. Kunk; Daniel M. Unruh; G. D. Zeihen; W. C. Hodges; Edward A. du Bray; J. Michael O’Neill

doi:10.2113/97.2.241

USGS · 70024823

SHRIMP U-Pb and 40Ar/39Ar age constraints for relating plutonism and mineralization in the Boulder batholith region, Montana

Karen Lund·John N. Aleinikoff·Michael J. Kunk·Daniel M. Unruh·G. D. Zeihen·W. C. Hodges·Edward A. du Bray·J. Michael O’Neill

View original source ↗Export metadata ↓

Abstract

The composite Boulder batholith, Montana, hosts a variety of mineral deposit types, including important silver-rich polymetallic quartz vein districts in the northern part of the batholith and the giant Butte porphyry copper-molybdenum pre-Main Stage system and crosscutting copper-rich Main Stage vein system in the southern part of the batholith. Previous dating studies have identified ambiguous relationships among igneous and mineralizing events. Mineralizing hydrothermal fluids for these types of deposits and magma for quartz porphyry dikes at Butte have all been considered to be late-stage differentiates of the Boulder batholith. However, previous dating studies indicated that the Boulder batholith plutons cooled from about 78 to 72 Ma, whereas copper-rich Main Stage veins at Butte were dated at about 61 Ma. Recent efforts to date the porphyry copper-molybdenum pre-Main Stage deposits at Butte resulted in conflicting estimates of both 64 and 76 Ma for the mineralizing events. Silver-rich polymetallic quartz vein deposits elsewhere in the batholith have not been dated previously. To resolve this controversy, we used the U.S. Geological Survey, Stanford, SHRIMP RG ion microprobe to date single-age domains within zircons from plutonic rock samples and 40 Ar/ 39 Ar geochronology to date white mica, biotite, and K-feldspar from mineral deposits. U-Pb zircon ages are Rader Creek Granodiorite, 80.4 ± 1.2 Ma; Unionville Granodiorite, 78.2 ± 0.8 Ma; Pulpit Rock granite, 76.5 ± 0.8 Ma; Butte Granite, 74.5 ± 0.9 Ma; altered Steward-type quartz porphyry dike (I-15 roadcut), 66.5 ± 1.0 Ma; altered Steward-type quartz porphyry dike (Continental pit), 65.7 ± 0.9 Ma; and quartz monzodiorite of Boulder Baldy (Big Belt Mountains), 66.2 ± 0.9 Ma. Zircons from Rader Creek Granodiorite and quartz porphyry dike samples contain Archean inheritance. The 40 Ar/ 39 Ar ages are muscovite, silver-rich polymetallic quartz vein (Basin district), 74.4 ± 0.3 Ma; muscovite, silver-rich polymetallic quartz vein (Boulder district), 74.4 ± 1.2 Ma; muscovite, early dark micaceous vein (Continental pit), 63.6 ± 0.2 Ma; biotite, early dark micaceous vein (Continental pit), 63.6 ± 0.2 Ma; potassium feldspar, early dark micaceous vein (Continental pit), 63 to 59 Ma; and biotite, biotite breccia dike (Continental pit), 63.6 ± 0.2 Ma. Outlying silver-rich polymetallic quartz veins of the Basin and Boulder mining districts probably are directly related to the 74.5 Ma Butte Granite, whereas Steward-type east-west quartz porphyry dikes and Butte pre-Main Stage deposits are parts of a 66 to 64 Ma magmatic-mineralization system unrelated to emplacement of the Boulder batholith. The age of the crosscutting Main Stage veins may be about 61 Ma as originally reported but can only be bracketed as younger than the 64 Ma pre-Main Stage mineralization and older than the about 50 Ma Eocene Lowland Creek intrusions. The 66 Ma age for the quartz monzodiorite of Boulder Baldy and consideration of previous dating studies in the region indicate that small ca. 66 Ma plutonic systems may be common in the Boulder batholith region and especially to the east. The approximately 64 Ma porphyry copper systems at Butte and gold mineralization at Miller Mountain are indicative of regionally important mineralizing systems of this age in the Boulder batholith region. Resolution of the age and probable magmatic source of the Butte pre-Main Stage porphyry copper-molybdenum system and of the silver-rich polymetallic quartz vein systems in the northern part of the Boulder batholith documents that these deposits formed from two discrete periods of hydrothermal mineralization related to two discrete magmatic events.

Explore related subjects

Economic Geology Article

Source-reported bounding extent: 45.57944511437787° to 46.837649560937464° latitude; -113.345947265625° to -111.6595458984375° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Karen Lund, John N. Aleinikoff, Michael J. Kunk, Daniel M. Unruh, G. D. Zeihen, W. C. Hodges, Edward A. du Bray, J. Michael O’Neill. 2002. SHRIMP U-Pb and 40Ar/39Ar age constraints for relating plutonism and mineralization in the Boulder batholith region, Montana. https://doi.org/10.2113/97.2.241

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Origin of the Pd/Pt ratio of the J-M Reef, Stillwater Complex, Montana, USA

The J-M reef of the Stillwater Complex is characterized by a high Pd/Pt ratio (mean ~3.8 with a standard error of 0.03) with a homogeneous geospatial distribution at the deposit scale. In this contribution, we demonstrate that the Pd/Pt ratio of the reef is the product of equilibration of an immiscible sulfide liquid with a silicate melt rich in Pd relative to Pt. Despite the high tenors of the J-M reef sulfides (avg 2,700 ppm Pt and 770 ppm Pt), numerical modeling shows that the parental melts did not have extraordinary Pd and Pt concentrations. Instead, the initial composition of a plausible parental silicate melt can have Pd and Pt contents well within the expected range of a normal, mantle-derived partial melt (i.e., ~10–20 ppb for both Pd and Pt with Pd/Pt of ~1). The relative differences in the partitioning behavior of Pt and Pd between sulfide liquid and silicate melt are unlikely to produce a consistent Pd/Pt ratio across a wide range of silicate melt to sulfide liquid mass ratios (i.e., R factors). Instead, the pre-emplacement fractionation of Pt alloy from S-undersaturated silicate magma accounts for the homogeneous and high Pd/Pt ratio of the J-M reef. We show that batch equilibration of sulfide liquid with silicate melt can produce the high Pd/Pt ratios of the reef if the partition coefficients between sulfide liquid and silicate melt for Pd and Pt are extremely high (>10 6 ). In an alternative model, Pd enrichment could be achieved by sulfide upgrading in resident footwall mush even if the partition coefficients between sulfide liquid and silicate melt are relatively small (between 10 4 and 10 6 ) because the instantaneous mass ratio of silicate melt to sulfide liquid is small (R ≈ 100–700), so the partitioning behavior of Pt and Pd has little impact on the composition of sulfide liquid.

Montana↗

USGS2025

Depth of magma crystallization and fluid exsolution beneath the porphyry-skarn Cu deposits at Santa Rita and Hanover-Fierro, New Mexico, USA

The depth level at which porphyry Cu–forming magmas fractionated and exsolved mineralizing fluids is actively debated. In the classic model, extensive magma fractionation occurs in large, upper crustal magma chambers, and concomitant fluid exsolution leads to forceful expulsion of residual magmas in the form of porphyry dikes, stocks, and breccia pipes, which subsequently serve as pathways for the mineralizing fluids. In contrast, some recent studies highlighting the role of deep crustal magma fractionation in the production of fertile magmas essentially deny the existence of upper crustal magma chambers at the time of mineralization. To address this, we conducted a detailed thermobarometric investigation of 13 intermediate to felsic, porphyritic intrusive rocks related to porphyry-skarn Cu mineralization at Santa Rita and Hanover-Fierro, New Mexico, United States, representing two premineralization magmas (61–60 Ma), seven synmineralization magmas (60–58 Ma), and four late- to postmineralization magmas (58–57 Ma). For each sample, the pressure of last magma crystallization before final magma ascent to the current exposure level was reconstructed based on Al-in-hornblende barometry of small hornblende inclusions trapped within quartz phenocrysts and through titanium-in-quartz (TitaniQ) thermobarometry of the host quartz phenocrysts themselves. Since quartz is one of the last crystallizing magmatic minerals, and no significant phenocryst growth could have occurred in small dikes and stocks after final magma emplacement, quartz phenocrysts and their contained hornblende inclusions record the depth of last magma crystallization before final magma ascent. When present, hornblende phenocrysts and hornblende inclusions within other major phenocrysts were also analyzed. Both quartz and hornblende barometers return consistent average pressures of 3.2 ± 0.4 kbar for the entire suite of pre- to postmineralization magmas, corresponding to depths of 11 to 14 km. The synmineralization magmas return even more consistent average pressures of 3.1 ± 0.2 kbar, corresponding to a depth of 12 ± 1 km. The volume of the mineralizing porphyry dikes and stocks at the emplacement level is far too small to have provided all the fluids and metals required to form the observed ore deposits. Therefore, the majority of the ore-forming fluids must have originated from the magmas that crystallized at 12 ± 1 km depth. The ore deposits, conversely, formed at ~5-km paleodepth. This implies that most of the mineralizing fluids traveled an average vertical distance of ~7 km from their magmatic source to the eventual site of ore precipitation. The relatively unaltered nature and low veining degree of deeper parts of mineralized porphyry dikes and stocks suggest that the fluid transport through these intrusive bodies occurred mostly at near-solidus conditions by means of fluid percolation along grain boundaries. In summary, our results suggest that (1) a large, upper crustal pluton exists ~7 km beneath the Santa Rita and Hanover-Fierro deposits; (2) abundant phenocryst crystallization occurred at this depth level; and (3) this pluton was the main source for the exsolution of ore-forming fluids. However, the investigated rocks have elevated whole-rock Sr/Y ratios, indicating magma fractionation at deep crustal levels. As a result, our preferred model is a combination of the two end-member models introduced above, with most magma fractionation having occurred in the deep crust and with residual, intermediate to felsic melts having ascended and accumulated at 11 to 14 km paleodepth, where they continued to crystallize with comparatively little crystal-liquid separation, before some of these magmas ascended further to shallow levels and quenched to porphyries.

New Mexico↗

USGS2025

Assimilation of reduced carbon triggers platinum alloy saturation in mafic and ultramafic magmas

It is generally observed that magmatic sulfide ores have higher ratios of Pd/Pt than the mantle-like values of their parental magmas. This discrepancy has defied simple explanation because the partitioning behavior of both elements between sulfide and silicate liquids is very similar. Assimilation of sulfur- and carbon-rich country rocks by mafic and ultramafic magmas is considered a critical, if not essential, step in the formation of magmatic base metal sulfide deposits. Although there is general consensus that the assimilation of external sulfur and carbon promotes sulfide saturation, the effect of carbon assimilation on the solubilities of platinum-group elements in natural S-bearing silicate melt has been overlooked. In this study, we investigate the variations of platinum and palladium solubilities during assimilation of graphite and methane through thermodynamic modeling, in comparison with data from an array of highly distinctive magmatic sulfide ore systems representing ages from Archean to Paleozoic, melt compositions from komatiite to basalt, and magmatic settings including lavas, hypabyssal intrusions, plutonic continental arc roots, and plutonic layered intrusions, namely: Raglan, Norilsk-Talnakh, Lac des Iles, and the J-M Reef of the Stillwater Complex. We model assimilation-fractional crystallization processes to estimate the reduction of oxygen fugacity ( ⁠fO2⁠ ) of the melt due to incorporation of graphite and methane. The simulations show that although Pd remains highly soluble during the progressive assimilation of reduced carbon, Pt solubility decreases significantly as the silicate melt becomes increasingly reduced. With less than 8 % of sediment assimilation, Pt alloy may saturate and then deviate from sulfide-undersaturated silicate melts, concomitantly increasing the Pd/Pt value of the remaining melts of the Raglan and Norilsk-Talnakh systems. For the Lac des Iles and Stillwater systems, a higher extent of assimilation is needed to reach Pt saturation because of the relatively carbon-poor nature of the lower crustal rocks. The assimilation of methane volatiles is shown to be more effective than graphite assimilation, and it provides a pathway to Pt alloy fractionation in the absence of detectable amounts of bulk host-rock assimilation. High Pd/Pt values have been documented in many world-class magmatic sulfide deposits whose parental magmas have demonstrably experienced crustal contamination. Our model suggests that although anomalous Pd/Pt values may be explained by other mechanisms such as incongruent melting of preexisting sulfide or differences in the diffusivities of the metals within achieving equilibration, the assimilation of graphite or methane may play an important role in the global occurrence of magmatic sulfide ores with elevated Pd/Pt values.

Economic Geology↗