Geology topics

Michael L. Williams

Publications and source records attributed to Michael L. Williams.

10 recordsLinked to original sources

USGS2025

Anatectic origin of Mississippian spodumene-bearing pegmatites in western Maine during orogenic plateau collapse

Spodumene pegmatites are an important lithium source, but the processes and tectonic settings in which they form are poorly understood. The Rumford pegmatite district surrounding Plumbago Mountain, western Maine, is host to numerous spodumene pegmatites, including the Plumbago North pegmatite (a world-class spodumene resource). Competing petrogenetic models for these spodumene pegmatites include (1) highly fractionated melts of the Mooselookmeguntic igneous complex and (2) anatexis. We tested these hypotheses by constraining the geologic, magmatic, metamorphic, and tectonic history of the Plumbago Mountain area with detailed geologic mapping and U-(Th)-Pb geochronology. The Silurian Rangeley Formation records initial isoclinal folding prior to, and contact-related metamorphism synchronous with, the intrusion of the 417 ± 4 Ma Plumbago Mountain pluton. Peak amphibolite facies metamorphism and crustal melting occurred during the ca. 410 to 400 Ma Acadian orogeny. Pulsed emplacement of the Mooselookmeguntic igneous complex occurred between ca. 389 and 356 Ma. Cassiterite U-Pb dates of spodumene pegmatites (333–327 Ma) are ≥23 m.y. younger than nearby granitic plutons, strongly arguing against the fractional crystallization model. Metamorphic monazite and xenotime (346–328 Ma) and 330 to 308 Ma 40 Ar/ 39 Ar hornblende dates indicate metamorphism coeval with spodumene pegmatite emplacement, supporting anatectic models. Reheating, anatexis, and spodumene pegmatite emplacement occurred during collapse of the 380 to 330 Ma Acadian orogenic plateau. Lithium enrichment may be linked to one or more stages of partial melting of metasedimentary and plutonic rocks during the formation, tenure, and collapse of the Acadian altiplano and emphasizes the role of anatexis in producing spodumene pegmatites of economic significance.

Maine↗

USGS2025

Re-evaluating the tectonic affinity of Proterozoic crustal provinces in the Southwest USA: Detrital zircon evidence for a Laurentian source for the Yavapai and Mojave Provinces

Models for crustal growth commonly involve the accretion of dominantly juvenile crust to continental margins. However, tracking the provenance and tectonic affinity of dominantly juvenile crustal provinces is challenging. This difficulty is highlighted by uncertainty over whether the Yavapai and Mojave Provinces, part of the >1300-km-wide system of Proterozoic orogens in southwestern Laurentia, (1) have similar crustal and tectonic histories and (2) if they formed on or near Laurentian, Australian, or Antarctic cratons. Here, we contribute new large-n detrital zircon U-Pb geochronology and Sm-Nd whole-rock isotope geochemistry to help constrain the provenance of the Yavapai Province and address these broader questions. Yavapai Province metasedimentary rocks from central Colorado in the southwestern USA have abundant pre-1.80 Ga detrital zircon grains, with ca. 1.85 Ga, 2.30 Ga, and 2.50−2.70 Ga peaks, and variable amounts of 1.79−1.78 Ga grains. Evolved whole-rock Sm-Nd isotopic compositions from these rocks, including 2.36−2.08 Ga model ages, also suggest mixing between 1.79 Ga and 1.78 Ga Yavapai Province arcs and early Proterozoic to Archean sources. Nearly identical pre-1.8 Ga detrital and inherited zircon age distributions suggest that the Yavapai and Mojave Provinces formed on and/or incorporated similar material. The Trans-Hudson orogen, and to a slightly lesser extent the Penokean orogen, provide the closest matches to the pre-1.80 Ga material in the Yavapai and Mojave Provinces. This similarity, coupled with a weaker resemblance to Australian and Antarctic sources, support a Laurentian affinity for the Yavapai and Mojave Provinces. We envision Paleoproterozoic arc formation on both oceanic crust and material of Laurentian affinity and multiple phases of arc-back-arc genesis, closure, and accretionary tectonism along the long-lived margin of the supercontinent Columbia (Nuna).

Arizona, Colorado, Idaho, New Mexico, Utah, Wyomin↗

USGS2025

Alteration mapping in granitic gneiss using handheld geophysical and geochemical instruments: Implications for iron oxide-apatite and rare earth elements exploration

The Adirondack Mountains of New York, U.S.A. contain iron oxide-apatite (IOA) mineral deposits with variable concentrations of rare earth elements (REE). The IOA mineral deposits are typically hosted in the Lyon Mountain Granite Gneiss and are spatially correlated with extensive Na metasomatism (albitization) of the surrounding country rocks, although some mineral deposits also occur in metagabbro, paragneiss, and anorthosite hosts. The location of albitization zones is key to finding new IOA mineral deposits and better understanding the mineralization processes associated with their genesis. However, the Na alteration zones are generally not visibly identifiable in outcrop or hand specimens because the color and textural changes are subtle and are thus difficult to map. Here, we discuss the results from testing two spectroscopic techniques (handheld, field gamma-ray spectroscopy, and portable X-ray fluorescence (pXRF)) to characterize albitization zones surrounding IOA mineral deposits. The gamma-ray and pXRF spectrometers successfully distinguished the relative intensity of alteration in the Lyon Mountain Granite Gneiss based on the K abundance. The measured K content decreases towards the IOA mineral deposits, and the magnitude and width of the gradient are similar along the strike of each mineral deposit. Elevated Th and U values are present in host rocks adjacent to REE-bearing IOA mineral deposits. The pXRF and gamma-ray spectrometer K values are remarkably consistent with laboratory-based, whole-rock XRF compositional data and, therefore, useful for semi-quantitative analysis. Notably, albitization aureoles are consistent around REE-rich and REE-poor mineral deposits, suggesting that saline fluids are essential for the petrogenesis of IOA mineral deposits but may not be critical to REE mineralization. Ultimately, this study demonstrates the utility of handheld gamma-ray and pXRF spectrometry for identifying otherwise cryptic albitization gradients associated with IOA mineral deposits in granitic gneiss.

New York↗

USGS2024

Paleoproterozoic reworking of Archean crust and extreme back-arc metamorphism in the enigmatic southern Trans-Hudson orogen

The crustal evolution of the southernmost ∼2000–1800 Ma Trans-Hudson orogen (THO) is enigmatic due to burial by Phanerozoic sediments. We provide new insights through petrochronologic analysis of a paragneiss drill core sample. Detrital zircon age peaks at 2625, 2340, and 1880 Ma and Hf isotopes suggest Paleoproterozoic arc development proximal to Archean source(s). Phase equilibria modeling and ternary feldspar thermometry suggest peak conditions of ≥1 GPa, ≥900°C, the first recognition of extreme, ultra-high temperature metamorphism in the THO. The largely isobaric P-T path, rapid heating rate, and ∼20 Myr duration (1872–1850 Ma) of peak conditions suggest that this metamorphism occurred in a back-arc tectonic setting. The sample records post-peak (1850–1815 Ma) mid-crustal residence, slow cooling, and exhumation. Further retrogression occurred during Proterozoic regional exhumation (1630–1470 Ma) and Phanerozoic (360–220 Ma) reheating and/or fluid influx. Evidence for Paleoproterozoic arc(s) supports geophysical data for Archean cratonic and Paleoproterozoic arc crust in this region.

Geophysical Research Letters↗

USGS2024

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↗

USGS2023

Monazite and xenotime petrochronologic constraints on four Proterozoic tectonic episodes and ca. 1705 Ma age of the Uncompahgre Formation, southwestern Colorado, USA

The Proterozoic tectonic evolution of the south-western USA remains incompletely understood due to limited constraints on the timing and conditions of the tectono-metamorphic phases and depositional age of metasedimentary successions. We integrated multi-scale compositional mapping, petrologic modeling, and in situ geochronology to constrain pressure-temperature-time paths from samples of Paleoproterozoic basement gneisses and overlying quartzites in southwestern Colorado, USA. Basement gneiss from the western Needle Mountains records metamorphic conditions of 600 °C at 0.75 GPa at 1764 ± 9 Ma and ~575 °C at 1741 ± 10 Ma. Gneiss sampled from drill core near Pagosa Springs, Colorado, records conditions of 700 °C at 1748 ± 9 Ma, 800 °C at 1.1 GPa at 1650 ± 40 Ma, 540 °C at 1570 ± 36 Ma, and 440 °C at 1424 ± 12 Ma. The Uncompahgre Formation was deposited at ca. 1705 Ma, as constrained by detrital monazite (1707 ± 8 Ma) and xenotime (1692 ± 40, 1725 ± 50 Ma), metamorphic xenotime (1650 ± 10 Ma), and published 40 Ar/ 39 Ar and detrital zircon data. Compositions of ca. 1705 Ma detrital monazite and xenotime are consistent with derivation from a garnet-bearing source in the Yavapai orogenic hinterland. The Vallecito Conglomerate and Uncompahgre Formation record macroscopic folding and greenschist-facies metamorphism at 1650 ± 10 Ma and temperatures of 270 °C to >570 °C at 1470–1400 Ma. Laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) zircon geochronology yielded dates of 1775 ± 18 Ma from the Twilight Gneiss and 1696 ± 7 Ma from the Bakers Bridge Granite, supporting previous isotope dilution–thermal ionization mass spectrometry (ID-TIMS) dates. The Eolus Granite yielded a date of 1463 ± 6 Ma, which is older than previous 1.44–1.43 Ga ID-TIMS dates. The newly dated granite of Cataract Gulch is 1421 ± 12 Ma. In situ analysis of detrital and metamorphic monazite and xenotime, igneous zircon, and quantitative thermobarometry, integrated with previously published constraints, indicate multiple tectonic episodes after the emplacement of 1800–1760 Ma arc-related rocks. The region experienced greenschist- to amphibolite-facies metamorphism (M1) from 1760 Ma to 1740 Ma, which was followed by the intrusion of granites at 1730–1695 Ma and deposition of the Uncompahgre Formation at ca. 1705 Ma, contemporaneous with the Yavapai orogeny. Metamorphism at 1680–1600 Ma was characterized by greenschist-facies conditions near Ouray, Colorado, and granulite-facies conditions near Pagosa Springs (M2) during the Mazatzal orogeny. From 1470 Ma to 1400 Ma, greenschist- to amphibolite-facies metamorphism (M3) and largely granitic plutonism occurred during the protracted Picuris orogeny. These results demonstrate the power of monazite and xenotime analyses to constrain depositional ages, provenance, and pressure-temperature-time ( P - T - t ) paths to resolve the compound orogenic history that is characteristic of most mountain belts.

Colorado↗

USGS2023

Geochemical evidence for evolving Proterozoic crustal thickness and orogenic styles in southwestern Laurentia

It has long been challenging for researchers to track the crustal thickness and mode(s) of crustal modification in ancient convergent margins, limiting evaluation of the tectonic styles and processes that modify continental crust during orogenesis. We present trace element igneous geochemical crustal thickness proxies that quantitatively track the crustal thickness evolution of the long-lived Proterozoic active margin in the southwestern U.S.A. We integrate these results with geobarometric data to constrain the mode of crustal modification. The data indicate a complex record of crustal thickness change in space and time and evolving orogenic styles. Geochemical proxies at 1.84–1.72 Ga are consistent with 20–40 km thick magmatic arcs that were locally thickened to ∼50 km during ∼1.75 Ga tectonism. During the Yavapai orogeny, 1.72–1.69 Ga, a ∼200-km-wide belt of 50-60 km thick crust extended from southern California to northern Colorado and was rapidly thinned and exhumed by ∼1.68 Ga. Crustal thickening and thinning during the Yavapai orogeny largely occurred by shortening and exhumation, respectively, in the upper 25 km of the crust. Subsequent 1.68–1.60 Ga tectonism involved crustal growth, local crustal thickening, and low-P, high-T metamorphism, consistent with extensional accretionary orogenesis. The 1.47–1.37 Ga Picuris orogeny was associated with 50–60 km thick crust across much of the Southwest and involved crustal shortening with ∼10 km of magmatic underplating. Advective heat from the emplacement of ferroan granites in the mid-crust likely contributed to elevated geothermal gradients and rheologically weakened the crust. Our results suggest evolving orogenic styles in the Southwest from 1.75–1.69 Ga short-lived crustal thickening associated with terrane accretion to 1.69–1.60 Ga largely extensional accretionary orogenesis, and regional, long-lived crustal thickening at 1.47–1.37 Ga involving extensive basaltic underplating. Contrasting with some recent hypotheses, our data document a complex middle Proterozoic record for the Southwest that was not orogenically quiescent or tectonically stagnant but involved complex mountain building styles.

Earth and Planetary Science Letters↗

USGS2023

The Arctic Rivers Project: Using an equitable co-production framework for integrating meaningful community engagement and science to understand climate impacts

As the Arctic and its rivers continue to warm, a better understanding of the possible future impacts on people would benefit from close partnership with Indigenous communities and scientists from diverse fields of study. We present efforts by the Arctic Rivers Project to conduct community-engaged research to increase collective understanding of the historical and potential future impacts of climate change on rivers, fish, and Indigenous communities. Working in central to northern Alaska and the Yukon Territory in Canada, the project seeks to engage with Indigenous communities in ethical and equitable ways to produces science that is useful, useable, and used that may serve as an example for future research efforts. Toward this goal, we formed an Indigenous Advisory Council and together developed project-specific knowledge co-production protocols. This paper provides a novel model of design and implementation to co-produce knowledge with communities across a large study domain.

Alaska↗

USGS2022

Geochemical evidence for diachronous uplift and synchronous collapse of the high elevation Variscan hinterland

Competing end-member models for the late Paleozoic Variscan orogeny (ca. 360-290 Ma) alternatively suggest moderate 2-3 km elevations underlain by relatively thin crust (<50 km) or a thick crust (>55 km) that supported high 4-5 km elevations. We tested these models and quantified the crustal thickness and elevation evolution of the Variscan orogeny using igneous trace element geochemical proxies. The data suggest that thick crust (55-70 km) capable of supporting 3-5 km elevations developed diachronously from east to west between ca. 350 and 315 Ma. Crustal thinning occurred from ca. 315 Ma to 290 Ma across the orogen. Crustal thickness and elevation changes at ca. 340-325 Ma and 315-290 Ma correspond with increases in silicate weathering recorded by Sr and Li isotopes, consistent with models in which silicate weathering of the Variscan orogen contributed to global cooling associated with the late Paleozoic ice age.

Black Forest↗

USGS2019

Syn-collisional exhumation of hot middle crust in the Adirondack Mountains (New York, USA): Implications for extensional orogenesis in the southern Grenville province

Extensional deformation in the lower to middle continental crust is increasingly recognized and shown to have significant impact on crustal architecture, magma emplacement, fluid flow, and ore deposits. Application of the concept of extensional strain to ancient orogenic systems, like the Grenville province of eastern North America, has helped decipher the structural evolution of these regions. The Marcy massif is a ~3000 km2 Mesoproterozoic anorthosite batholith in the Adirondack Mountains (New York, USA) of the southern Grenville province. Bedrock geology mapping at 1:24,000 scale paired with characterization of bedrock exposed by recent landslides provides a glimpse into the structural architecture of the massif and its margin. New data demonstrate granulite- to amphibolite-facies deformational fabrics parallel the margin of the batholith, and that the Marcy massif is draped by a southeast-directed detachment zone. Within the massif, strain is localized into mutually offsetting conjugate shear zones with antithetic kinematic indicators. These relationships indicate that strain was coaxial within the Marcy massif, and that subsimple shear components of strain were partitioned along its margin. In situ U–Th–total Pb monazite analysis shows that deformation around and over the Marcy massif occurred from 1070 to 1060 Ma during granulite-facies metamorphism, and monazite from all samples record evidence for fluid-mediated dissolution reprecipitation from 1050 to 980 Ma. We interpret that rocks cooled isobarically after accretionary orogenesis and emplacement of the anorthosite- mangerite- charnockitegranite plutonic suite at ca. 1160–1140 Ma. Gravitational collapse during the Ottawan phase of the Grenville orogeny initiated along a southeast-directed detachment zone (Marcy massif detachment zone), which accommodated intrusion of the Lyon Mountain Granite Gneiss, and facilitated substantial fluid flow that catalyzed the formation of major ore deposits in the Adirondack Highlands.

New York↗