USGS ScienceSearch

USGS · 70211684

Mapping the 3-D extent of the Stillwater Complex, Montana—Implications for potential platinum group element exploration and development

Abstract

Geophysical models characterize the exposed and interpreted buried extent of the Stillwater Complex, critical for understanding the origin of the layered mafic intrusion and its associated high-grade platinum group element resources. The 3D models, constrained by gravity, magnetic, xenolith, seismic, borehole, and rock property data indicate that the likely maximum extent of the Stillwater Complex beneath Phanerozoic cover is ~10 times greater than its outcrop, ~2240 km 2 . The thickness values are poorly constrained but vary from ~7000 to 12,000 m, depending on crustal and mantle density variations and depths to the top of the lower crust and mantle. This thickness may include dense metasedimentary units of the basin into which the Stillwater Complex intruded. Using the modeled thickness results in a volume estimate of ~24,700 km 3 , albeit poorly constrained. New analyses of xenoliths from the Cretaceous Sliderock and Suzie Peak intrusions produce ages of 2706–2716 Ma, corresponding to the age of the Stillwater Complex, and 2813 Ma, corresponding to the age of Archean gneissic basement. Seismic reflectors in inferred Archean crystalline basement, possibly including the Stillwater Complex, dip ~25–30° north, with segments dipping as much as 70° north. Layered reflectors beneath the Phanerozoic sedimentary section and above the inferred Archean crystalline basement may represent metasedimentary units, perhaps a southern extension of the Mesoproterozoic Belt Basin. The potential field models and the seismic reflection data suggest that the Stillwater Complex was dipping northward prior to deposition of Cambrian strata, perhaps uplifted in the late Archean or Proterozoic as previously proposed, and that during Laramide times, faulting and intrusions highly disrupted the complex. Temperature measurements from boreholes help constrain the depths of feasible mining of the complex.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 45.0502402697946° to 45.98169518512228° latitude; -111.70898437499999° to -108.7646484375° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Carol A. Finn, Michael L. Zientek, Heather L. Parks, Dana E. Peterson. 2020. Mapping the 3-D extent of the Stillwater Complex, Montana—Implications for potential platinum group element exploration and development. https://doi.org/10.1016/j.precamres.2020.105860

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

KEEP EXPLORING

Related USGS reports

Petrogenesis and mineralization potential of spinifex komatiitic basalts in the Bradley Peak greenstone terrane, Wyoming Province

Komatiitic volcanic rocks are important hosts of Ni sulfide mineralization and record early Earth evolution; however, those in the well-studied Archean Wyoming Province have received little attention. Here, we elucidate the timing and petrogenesis of the Bradley Peak komatiitic volcanic rocks using field and textural observations, geochronology, and geochemistry. Detrital and igneous zircon U-Pb ages for two samples from previously undated units support published age determinations, placing the eruption age at 2.72 Ga. Stratigraphy of the volcanic flows was mapped and 36 samples including cumulates, greenschists, and spinifex-textured rocks were collected. Whole-rock geochemistry was used to classify the spinifex-textured samples as Al-undepleted komatiitic basalts (11–17 wt% MgO). Platinum-group element concentrations (n = 25) are like those in global Al-undepleted komatiitic basalts, and PGE/Ti ratios do not indicate the volcanic flows likely host sulfide mineralization. Initial ε Nd values of −0.5 to +4.7 (n = 16), indicate that these lavas were derived from a depleted mantle source and have negligible evolved crust contamination. The primary magma to the komatiitic basalt flows is estimated to have had 19 wt% MgO and be derived from ∼15 to 25 % mantle partial melting at 3–4 GPa. Trace element chemistry and thermodynamic modeling suggest the primary melt assimilated local banded iron formation. Although the Bradley Peak komatiitic basalts do not contain positive evidence of magmatic sulfide deposits, depleted Au in the flows suggests they could be source rocks for nearby orogenic gold deposits.

Wyoming

Heterogeneous multi-stage accretionary orogenesis — Evidence from the Gunnison block in the Yavapai Province, southwest USA

Proterozoic rocks exposed in the southwestern U.S.A. represent one of the best examples of crustal growth by arc-related magmatism and accretionary orogenesis. Within the Southwest the 1.8–1.7 Ga Yavapai Province is widely regarded as a classic example of juvenile arc crust, however 1.8–2.5 Ga inherited zircon and Nd and Hf model ages have been recognized near Gunnison in central Colorado. These data have led to questions regarding the extent and nature of pre-1.8 Ga crustal material and the genesis of the Yavapai Province. We present evidence for a geochemically distinct, spatially restricted crustal block underlain by pre-1.8 Ga crust material (referred to here as the Gunnison block) in central to western Colorado within the Yavapai Province. The Gunnison block is characterized by 1.8–1.9 and 2.4–2.6 Ga inherited zircon, Pb isotopic systematics (μ = 9.8 ± 0.1, κ = 3.7 ± 0.1) elevated relative to 1.8 Ga depleted mantle values, 1.8–2.5 Ga Nd and Hf model ages, and a distinct pressure-temperature-time history. The geochemical data are consistent with mixing between juvenile 1.8 Ga and pre-1.8 Ga sources. The older crustal component is most similar to the isotopically enriched Mojave Province of eastern California and western Arizona, suggesting greater similarities between these provinces than previously recognized. Monazite and xenotime petrochronology indicate ca. 1.75–1.74, 1.72–1.69, 1.67, and 1.47–1.38 Ga tectono-metamorphic events. These data suggest that the Gunnison block accreted to other components of the Yavapai Province outboard of Laurentia at 1.75–1.74 Ga. The composite Yavapai Province was accreted to the margin of Laurentia during the 1.72–1.69 Ga Yavapai orogeny. Later overprinting is associated with the ∼1.68–1.60 Ga Mazatzal and ∼1.47–1.37 Ga Picuris orogenies. Identification of distinct crustal terranes within the Yavapai Province supports models involving multiple arcs and back-arcs that were progressively assembled prior to their accretion to Laurentia, perhaps akin to the present-day Banda Sea in Indonesia.

Yavapai Province

Deep magmatic staging chambers for crustal layered mafic intrusions: An example from the Bushveld Complex of southern Africa

The deep mafic magmatic staging chambers of layered mafic intrusions have been conjectured but not imaged. Their existence has long been postulated from geochemical models which require multiple magma injections from staging chambers to account for their multi-scale igneous layering and variations in sources and degrees of crustal contamination. For the Bushveld Complex of southern Africa, the world’s largest layered mafic intrusion, seismic receiver functions identify a diffuse crust-mantle transition beneath the Complex that suggests high velocity lower crust and/or uppermost lithospheric mantle. Here we present 3D gravity modelling of the Bushveld Complex that includes dense material at the crust-mantle boundary, imaging for the first time, remnants of magma staging chambers. They underlie the whole Bushveld Complex and extend westwards to the Molopo Farms Complex in Botswana. Feeders to the Bushveld Complex coincide with intersections of major faults like the Thabazimbi-Murchison Lineament and Sugarbush Fault with the staging chamber. This identification of magmatic staging chambers beneath the Bushveld relates to similar geophysically imaged lower crustal features beneath the Duluth and Stillwater Complexes. Comparison of seismic and gravity data with geochemical models from other complexes aid in development of models for the magmatic architecture of layered mafic intrusions in general.

Precambrian Research