Geology topics

Ian W. Hillenbrand

Publications and source records attributed to Ian W. Hillenbrand.

At least 19 recordsLinked to original sources

USGS2026

USGS Geochron—A database of geochronological and thermochronological dates and data—Technical documentation

Geochronological and thermochronological data are essential for constraining the timing and rates of geological processes, supporting geologic mapping, natural hazard assessment, and resource exploration. The U.S. Geological Survey (USGS) Geochron Database is a centralized, relational database that integrates USGS and State geological survey data in accordance with the National Geologic Mapping Act of 1992 and its reauthorizations. This report documents the structure and development of the database, including its standardized schema, controlled vocabularies, and protocols for compiling legacy and newly published data. The database is designed to adhere to FAIR (Findable, Accessible, Interoperable, and Reusable) data principles and currently (2026) includes data from a wide range of geochronological and thermochronological methods. Public access is provided through the USGS Geochron Database Explorer, a browser-based geographic information system (GIS) interface, and through versioned data releases available on ScienceBase, which provide sample-level summary data and detailed analytical information. The USGS Geochron Database is structured to be extensible, supporting the inclusion of additional methods and data types as they are compiled. This format ensures long-term utility and aligns with the Mapping Act’s directive to create a national archive of geochronological information that adds interpretive value to geologic map data. Its design reflects a commitment to data transparency, scientific reproducibility, and national geoscience priorities.

Data Report↗

USGS2026

Zircon petrochronology fingerprints mantle sources of magmatic rare earth element deposits

The mantle sources of carbonatite and alkaline magmatic systems that host rare earth element (REE) deposits are difficult to constrain because whole-rock compositions are commonly altered. Alteration-resistant zircon trace-element compositions are used to fingerprint mantle sources at Mountain Pass, California, the largest REE deposit in the United States, and extend this approach to a global carbonatite−alkaline rock dataset. Autocrystic zircon samples from the Mountain Pass intrusive suite record subduction-like U-Sc-Nb-Yb compositions, oxidized fO2 (0.8−2.4 relative to fayalite-magnetite-quartz buffer [ΔFMQ]), and Ti-in-zircon temperatures up to ∼970 °C, indicating a lithospheric mantle source enriched and oxidized by Paleoproterozoic subduction and subsequently tapped during Mesoproterozoic postcollisional lithospheric thinning. This contrasts with the ocean island basalt−like enriched mantle recorded by Mesoproterozoic zircon from the Bayan Obo carbonatite complex, China, demonstrating that world-class REE deposits of comparable age can reflect fundamentally different mantle sources. New classifiers separate subduction-metasomatized from non-subduction-enriched mantle domains with ∼90% accuracy, providing a zircon-based tool for fingerprinting the mantle source regions of magmatic REE deposits.

Geology↗

USGS2026

USGS Geochron Database

Introduction Geochronology helps us understand Earth’s history by determining when important events, like volcanic eruptions, the rise of mountains, the formation of mineral resources, and changes in the landscape, happened. Geochronological data directly support geologic mapping and can inform decisions about geologic hazard mitigation, natural resource management, and infrastructure resilience. The U.S. Geological Survey (USGS) Geochron database provides access to more than 300,000 published, publicly available age measurements from more than 40,000 geological samples. This database is the result of a collaborative effort with State geological surveys and geoscientists from across the globe. The USGS Geochron database is the most comprehensive collection of geochronological data available for the United States. Users can view data through an interactive map explorer, download datasets, and integrate data into geospatial software or other analysis tools.

Fact Sheet↗

USGS2026

The use of fluorite geochemistry and machine learning to identify critical mineral systems

Fluorite (CaF 2 ) is a potential pathfinder to critical mineral and rare earth element (REE) deposits but its application has been limited to a narrow range of mineralization types. I show that fluorite is a robust recorder of mineralization fertility by applying statistical and machine-learning methods to a new global fluorite geochemical database. Distinct median rare earth and trace element patterns are observed among deposit types and genetic environments. Fluorite associated with carbonatites and REE deposits are relatively enriched in Sr and have minimal Eu anomalies. These characteristics define new bivariate discrimination diagrams that correctly identify 78% of carbonatite-related fluorite and 88% of fluorite from REE deposits. Random forest classifiers were developed for a wide range of mineralization types and genetic settings. Trained solely on rare earth element patterns, these models achieve accuracies of 77–79%. Higher classification accuracies (up to 88–96%) are obtained when including elements such as Sr, highlighting the significance of trace elements for optimal fluorite classification. The recognition of diagnostic fluorite compositional fingerprints, particularly in REE-fertile systems, underscores its potential as a pathfinder and indicator for critical mineral exploration in F-bearing environments.

Mineralium Deposita↗

USGS2026

Tectonically driven integration of the 4.8 Ma Colorado River USA tracked with detrital sanidine and fish genetics

The development of the continental-scale Colorado River system, western USA, from 8 to 4.8 Ma, is revealed using 60-40 Ma detrital sanidine tracer grains and fish phylogeny. Here we show that precursor paleoriver segments became integrated north to south as traced by 60-40 Ma sand grains that were derived from the north and sequentially appeared in the 25-8 Ma Browns Park Formation of Utah, 7-6 Ma upper Bidahochi Formation of Arizona, and 4.8 Ma Bouse Formation of the lower Colorado River and proto Gulf of California. This timing is mimicked by molecular clock estimates of divergence times among fish lineages. River integration was a response to headwater uplifts in the Yellowstone hotspot track and Rocky Mountains. 40 Ar/ 39 Ar ages refine the timing for mantle-drips that caused subsidence, then uplift, of depositional basins that influenced the integration pathway and tempo. The ~ 3-million-year timescale suggests that multiscale mantle-driven uplift, rather than lake spillover, was the primary driver for integration of the proto-Colorado River through Grand Canyon.

Colorado River↗

USGS2026

Unraveling protracted modification of Archean and Paleoproterozoic crust in central Laurentia, Penokean orogen, with garnet and accessory mineral geochronology and microstructural analysis

Proterozoic metamorphism and deformation of the southern margin of the Superior craton in the Lake Superior region is attributed to the Penokean orogeny (1890−1830 Ma). This model includes a period of crustal inversion in which Archean basement blocks were exhumed through overlying Paleoproterozoic strata, producing the corridor of gneiss domes that parallels the trend of the Penokean orogen across the northern Midcontinent, USA. However, recent geologic mapping and 40 Ar/ 39 Ar geochronology challenge this interpretation, suggesting instead that the gneiss dome structures reflect younger episodes of tectonic activity along the southern margin of Laurentia. In absence of integrated pressure-temperature-time-deformation constraints for these rocks, interpretations are largely limited to their final cooling history, making it difficult to both identify the tectonic forces that shaped the architecture of the Penokean orogenic belt and assess the extent to which later Proterozoic tectonism modified the southern Superior craton. We address this problem with an approach joining thermodynamic modeling, garnet and accessory mineral geochronology, and microstructural analysis for several metamorphic rocks across the gneiss dome corridor. The U-Pb ages of titanite reveal that the Proterozoic geometries of exhumed basement gneiss domes are governed by preexisting Archean structures. Garnet Lu-Hf geochronology constrains the timing of prograde-to-peak metamorphism in the Penokean orogenic belt. Granulite facies metamorphism is related to the final stages of the Penokean orogeny at 1837 Ma and localized in a belt of high-grade rocks near a major Penokean suture. Garnet Lu-Hf ages of samples adjacent to gneiss domes reflect regional metamorphism following the accretionary phase of the Penokean orogeny, between 1825 Ma and 1782 Ma, which we suggest reflects continued crustal thickening related to convergence farther south during this time interval. Combination of garnet microstructures and Sm-Nd ages reflects later exhumation of gneiss domes and buried metasedimentary rocks by ca. 1750 Ma, consistent with previously published 40 Ar/ 39 Ar cooling ages across the region. Reset Lu-Hf and Sm-Nd garnet ages and U-Pb ages of syn-kinematic titanite reflect reactivation of primary Penokean structures during this period of basement uplift. These data document significant modification of the Penokean orogen and the Archean crust of the southern Superior province between 1800 Ma and 1700 Ma. Tectonic activity during this interval coincides with collisional events recognized in western Laurentia, suggesting that the period immediately following the Penokean orogeny may be a broadly important time for crustal growth and modification in proto-North America.

Michigan, Wisconsin↗

USGS2026

Reconstructing ancient sedimentary source-to-sink systems – Examples from southern Laurentia’s Proterozoic accretionary orogens

Provenance analysis is a powerful tool for investigating sediment delivery networks, constraining magmatic histories, and reconstructing the tectonic evolution of orogenic belts and basins. Basin analysis studies increasingly use detrital zircon (DZ) U-Pb forward mixture modeling to enhance provenance interpretations by quantifying the relative contributions of different sources. Forward mixture modeling requires significant a priori knowledge that limits deep-time applications. This challenge is overcome with an inverse mixture modeling approach non-negative matrix factorization to reconstruct the number and age distributions of paleo-source regions of Proterozoic metasedimentary rocks in the southwestern United States. This analysis indicates eight reconstructed end-member distributions representing unique sediment sources: two multi-modal end members characterized by ages older than ca. 1.8 Ga from cratonic Laurentia, five unimodal age distributions between ca. 1.80 Ga and 1.65 Ga consistent with Paleoproterozoic arc magmatic sources, and a ca. 1.6−1.5 Ga end member likely derived from exotic cratons in supercontinent Nuna (Columbia). Sediments deposited between ca. 1.80 Ga and 1.73 Ga yield heterogeneous age distributions suggesting multiple arc-backarc systems and several phases of slab roll back, contraction, and accretionary orogenesis, including input from pre−1.8 Ga Laurentian cratons. Homogenization of DZ signatures during the Yavapai orogeny (ca. 1.72−1.68 Ga) reflect crustal assembly as well as the uplift of Paleoproterozoic arcs in the orogenic hinterland. Detrital zircon age distributions from strata deposited during the Mazatzal orogeny (ca. 1.65−1.60 Ga) suggest the Mazatzal Province is a continental arc constructed on older crust. Mesoproterozoic samples are consistent with multiple basins derived from local recycling and long-distance sediment transport. Collectively, these data record the tectonic transition from the episodic accretion of disparate crustal domains to an increasingly integrated continental margin. These results provide new insights into the Proterozoic tectonic and paleogeographic evolution of the southwestern United States at basin to orogen scales and highlight the power of inverse DZ modeling to extract geologically meaningful quantitative mixture models from sedimentary records alone, offering a powerful tool for deep-time tectonic and basin analysis.

Arizona, California, Colorado, New Mexico, Nevada,↗

USGS2025

Finding the hidden orogeny – The Proterozoic polymetamorphic history of northern New Mexico

Pressure–temperature–time-deformation histories provide key constraints on orogenic processes but can be affected by later overprinting. This is exemplified in the Proterozoic orogenic belts of southwestern Laurentia where competing tectonic models involve either a single progressive Mesoproterozoic event, the Picuris orogeny, or a polyorogenic history that also includes the ~1.65 Ga Mazatzal orogeny. We address this controversy with structural analysis and petrochronology in the type locality of the Picuris orogeny. Xenotime and monazite domains associated with the greenschist-facies axial planar fabric of early F 1 folds yield 1644 ± 11 Ma (xenotime) and 1641 ± 15 Ma (monazite) ages and 450°C–482°C temperatures recording metamorphism and shortening associated with the Mazatzal orogeny. This Palaeoproterozoic greenschist-facies assemblage was overprinted by higher grade, protracted (1470–1350 Ma) tectonism that included the Picuris orogeny. Our results document the complex polyphase crustal assembly of Laurentia and highlight how petrochronology can effectively see through higher grade overprints to identify a more complete orogenic evolution.

New Mexico↗

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

Geochemistry and radiogenic isotopes constrain the mantle source region of the Mountain Pass Intrusive Suite, California

The Mountain Pass carbonatite stock is the largest rare earth element (REE) deposit and only active REE mine in the United States. The carbonatite intrusion and spatially associated alkaline silicate intrusions constitute the Mountain Pass Intrusive Suite, which is located within the Mojave Province in California. Both the carbonatite and the alkaline silicate rocks are enriched in large ion lithophile elements and light REEs, and less enriched to depleted in high field strength elements, indicating the mantle source region was metasomatically enriched in incompatible trace elements. The cause of this metasomatic mantle enrichment and the genetic relationship between the carbonatite and the alkaline silicate stocks are poorly understood. In this study, major and trace element geochemical data and isotopic (Rb-Sr, Sm-Nd, and Lu-Hf) data are presented to constrain genesis of the Mountain Pass Intrusive Suite, from mantle source region to the intrusion of the stocks. Our geochemical data are consistent with derivation of the alkaline silicate and carbonatite melts through partial melting from a shared mantle source region rather than through liquid immiscibility or fractional crystallization and separation of a carbothermal fluid. Although the Rb-Sr isotopic system in the Mountain Pass Intrusive Suite is disturbed at the whole-rock scale, the isotopic systems for whole-rock Sm-Nd (εNd i = ‐2.2 ± 0.8) and zircon Lu-Hf (εHf i = 0.1 ± 1.1) are robust and support mantle derivation of the magmas. Geochemical modeling using experimentally derived partition coefficients was used to identify possible causes of enrichment in incompatible elements through metasomatism in the mantle source region. Modeling of metasomatism by melts derived by partial melting of deeply subducted carbonated sediments approximates observed Mountain Pass Intrusive Suite trace element chemistry. Scattered εHf i in inherited zircon (2.8 ± 2.6) is consistent with derivation from an arc-related environment with substantial crustal contamination. Paleotectonic studies in the Mojave Province indicate that regional subduction preceded emplacement of the Mountain Pass Intrusive Suite by ∼300 Ma. Melting of the Mountain Pass source region may have been caused by post-collisional thermal relaxation and extension.

California↗

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↗

USGS2025

Multi-stage crustal thickening, surface uplift, and collision in the western Himalaya-Karakoram-Tibet orogen revealed by chemical mohometry

The Himalaya-Karakoram-Tibet (HKT) orogen provides an unrivaled opportunity to study the dynamic linkages between deep and surface processes during collisional orogenesis. However, these efforts are hindered by conflicting interpretations on the number and timing of collisional events, and the timing of crustal thickening and associated surface uplift. Here, we resolve this with quantitative paleo-crustal thickness estimates in the northwestern HKT orogen. We show that: (a) the paleo-Asian margin had thick crust (50–60 km) at least 65 Ma prior to terminal collision, consistent with a continental arc setting, (b) crustal thickening to 60 km or more occurred at ca. 60–50 Ma in the Kohistan-Ladakh arc and by 40–25 Ma in the paleo-Asian margin, indicating a multi-stage Himalayan collision, and (c) modern crustal thicknesses in the northwestern HKT have been sustained since ca. 40–25 Ma suggesting an orogenic steady-state in which crustal thickening, crustal flow, and surface uplift have been balanced by erosion.

Geophysical Research Letters↗

USGS2025

Zircon as a pathfinder to REE mineralization

Carbonatites and alkaline silicate rocks are major primary sources of the rare earth elements (REE) and other critical metals, such as Nb. Despite the economic significance of these rocks, their formation and the processes of REE enrichment are poorly understood. Here, statistical analysis of a global dataset demonstrates that zircon geochemistry is a powerful recorder of REE metallogenesis and a potential pathfinder for REE deposits. Zircons from REE and Nb fertile intrusions lack Eu anomalies and have elevated Gd/Yb and Th/Yb, indicating they crystallised from magmas that originated from deep, oxidised and enriched mantle sources. Complexes with Nb enrichment have low U/Nb, reflecting an enriched mantle source, whereas high U/Nb in REE-only fertile intrusions suggest a subduction-metasomatised mantle source. Machine learning models demonstrate high accuracy in classifying zircon from barren and fertile deposits. Classification of detrital zircons shows that REE-enriched deposits correlate with supercontinent assembly, whereas Nb fertile complexes are associated with supercontinent breakup. This approach offers a new, mineral to global scale, petrologic and exploration tool that enhances understanding of REE metallogenesis.

Geochemical Perspectives Letters↗

USGS2025

Fingerprinting magmatic REE deposit sources with zircon petrochronology

Carbonatites and associated alkaline silicate rocks are of considerable economic interest due to their enrichments in rare earth elements. The petrogenesis and source(s) of these complexes, however, are poorly understood. Models propose either mantle plume-derived carbon-rich melts or a mantle source enriched by subduction-related metasomatism. We use zircon trace elements to fingerprint carbonatite-alkaline complex sources, focusing on the economically significant 1.4 Ga Mountain Pass intrusive suite (MPIS). Autocrystic zircon from MPIS alkaline silicate rocks are enriched in Th and U relative to Nb suggesting a subduction influenced, oxidized source region. Lower Sc/Yb, higher Ti concentrations, and the absence of Eu anomalies in autocrystic zircon suggest derivation from a less hydrous, hotter, and deeper mantle source relative to arc-related 1.8–1.6 Ga inherited zircon. These data are best explained by syn- to post-collisional Mesoproterozoic reactivation of a subduction-metasomatized mantle source. The source and tectonic setting of the MPIS contrasts with Mesoproterozoic Gifford Creek and Bayan Obo carbonatites which exhibit plume-like compositions and are associated with rifts, suggesting different sources and tectonic settings for these economically significant deposits. We demonstrate that zircon petrochronology is a robust method for distinguishing carbonatite sources and can inform more targeted exploration strategies for critical mineral resources.

Conference Paper↗

USGS2025

Enhancing mineral systems exploration through geochronology, thermochronology, and isotope analysis: USGS Geochron and USGS Isotope databases

A mineral systems approach to mineral exploration provides a comprehensive framework for understanding ore deposit formation by examining the geodynamic, magmatic, hydrothermal, and sedimentary processes responsible for mineralization, alteration, and remobilization of economic mineral deposits. Temporal and thermal constraints on ore genesis are crucial for refining mineral system models and guiding predictive exploration strategies. Geochronology and thermochronology offer invaluable insights into the timing and thermal evolution of ore-forming processes, whereas isotopic analyses provide critical information on the source and geochemical history of ore-forming fluids. Combining these methodologies have proven highly effective for mineral exploration in regions like Australia, however, their combined application has been limited in the United States. To apply these tools to mineral systems-based exploration, the U.S. Geological Survey (USGS) has developed two products: (1) The USGS Geochron Database, and (2) the USGS Isotope Database. These databases provide centralized repositories of geo/thermochronological dates and data (Geochron Database) and both radiogenic and stable isotope data (Isotope Database) generated by the USGS and partners over the past decades. Integrating these datasets together and with traditional exploration approaches provides the mineral exploration community with powerful tools for determining the temporal and thermal histories of ore systems and identifying metallogenic source provinces.

Continental United States↗

USGS2024

Lead isotopes constrain Precambrian crustal architecture, thermal history, and lithospheric foundering in Laurentia

Laurentia (ancestral North America) records nearly 4 billion years of crustal evolution. Here, a newly compiled continental-scale Pb isotopic database is used to evaluate the Precambrian crustal evolution of Laurentia. Pb model ages yield a 2.7 Ga peak, a 2.5–1.8 Ga minimum and 1.8–0.9 Ga continuum. Pb model ages yield thermochronometric data and track crustal growth via arc-related magmatism and accretionary orogenesis. Model 232 Th/ 204 Pb and 238 U/ 204 Pb broadly correlate with mapped crustal domains. More homogeneous and less radiogenic 238 U/ 204 Pb and 232 Th/ 238 U after 2.7 Ga suggests a shift to more juvenile sources, loss of early isotopic reservoirs and greater crustal reworking. U and Th are fractionated from Pb in Proterozoic orogens with abundant ferroan and anorthosite–mangerite–charnockite–granite(AMCG)-suite magmatism. This fractionation suggests the removal of Pb-rich lower crust, supporting petrogenetic models involving lithospheric foundering and magmatic underplating. Lithospheric thinning and associated magmatism may have contributed to high middle Proterozoic geothermal gradients.

Terra Nova↗

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↗