USGS Science⌕ Search

SEARCH · USGS Science

Results for “Frontiers in Ecology and Environment”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,673 records · Page 76Linked to original sources
USGS2025

The effects of line simplification on planform geometry

Data on maps should retain accuracy regardless of scale. Yet, as cartographic lines are generalized, there can be impacts on properties such as topology, density, and planform geometry. Here, we investigate the use of the Scale Specific Sinuosity (S3) metric (Stanislawski et al., 2023) to evaluate the effects of line simplification on planform geometry, which is the bends of streams in map view. We employ an open-source Python S3 workflow to characterize the geometry of five diverse stream channels in the United States. The original data are extracted from the U.S. Geological Survey National Hydrography Dataset 1:24,000-scale vector data (U.S. Geological Survey, 2000) (Table 1), and the simplification is done using the Visvalingam and Whyatt method (2017) with a simplification tolerance of 0.5, 1.0, 1.5, and 2.0 km. The S3 analysis is calculated at each level of simplification and S3 derivatives are generated. Derivatives include measures of sinuosity, fractal dimension, and the dominant bend wavelength. The findings show that the change in planform geometry is scale-dependent, though simplification will have little effect on straighter lines. The change becomes more apparent in complex lines as the degree of simplification aligns with the scale of the dominant bend geometries. These logical conclusions are evidence that the S3 is a useful metric for automated characterization of bend geometry regardless of line complexity.

Conference Paper↗
USGS2024

Deep syntectonic burial of the Anthracite belt, Eastern Pennsylvania

Fluid inclusion microthermometry and Raman spectroscopy of fluid inclusions in quartz veins from the Pennsylvanian rocks of the Anthracite belt, eastern Pennsylvania support a deep burial model of coalification in favor of focused orogenic hot fluid flow. High-temperature (250 to 255 °C) trapping of CH 4 ± CO 2 saturated aqueous fluids and CH 4 ± CO 2 inclusions indicate fluid trapping at depths of 11.5 to 13.4 km under a cover of Pennsylvanian to Permian(?) syntectonic load. In the folded rocks to the south of the Anthracite belt, CH 4 ± CO 2 fluid inclusions indicate a sediment load that was up to 16.3 km thick. Re -equilibrated aqueous fluid inclusions from veins in Silurian through Devonian rocks give the same range of trapping conditions but a wide range of fluid salinities suggesting that folding, fracturing, and meteoric recharge resulted in the intermixing of fluids from throughout the stratigraphic succession.

Pennsylvania↗
USGS2025

Stratigraphic notes—Volume 2, 2025

This is the second volume in the U.S. Geological Survey (USGS) series of reports on stratigraphy entitled “Stratigraphic Notes,” which consists of short papers that highlight stratigraphic studies, changes in stratigraphic nomenclature, and explanations of stratigraphic names and concepts used on published geologic maps. “Stratigraphic Notes” is a long-term (multiyear), multivolume publication containing articles that address updates or revisions to stratigraphic nomenclature (and whose content ultimately will be incorporated by National Geologic Map Database personnel into Geolex, https://ngmdb.usgs.gov/Geolex/ ). We welcome papers for the “Stratigraphic Notes” series from geoscientists of the USGS, of State Geological Surveys, and from academicians. Papers can be submitted for publication in “Stratigraphic Notes” by contacting the USGS Geologic Names Committee ( gnc@usgs.gov ). As new “Stratigraphic Notes” volumes are published, links to the volumes will be posted at https://doi.org/10.3133/pp1879 .

Professional Paper↗
USGS2026

Stratigraphic notes—Volume 3, 2026

This is the third volume in the U.S. Geological Survey (USGS) series of reports on stratigraphy entitled “Stratigraphic Notes,” which consists of short papers that highlight stratigraphic studies, changes in stratigraphic nomenclature, and explanations of stratigraphic names and concepts used on published geologic maps. “Stratigraphic Notes” is a long-term (multiyear), multivolume publication containing articles that address updates or revisions to stratigraphic nomenclature (and whose content ultimately will be incorporated by National Geologic Map Database personnel into Geolex, https://ngmdb.usgs.gov/Geolex/ ). We welcome papers for the “Stratigraphic Notes” series from geoscientists of the USGS, of State Geological Surveys, and from academicians. Papers can be submitted for publication in “Stratigraphic Notes” by contacting the USGS Geologic Names Committee ( gnc@usgs.gov ). As new “Stratigraphic Notes” volumes are published, links to the volumes will be posted at https://doi.org/10.3133/pp1879 .

Professional Paper↗
USGS2024

Aftershock forecasting

Aftershocks can compound the impacts of a major earthquake, disrupting recovery efforts and potentially further damaging weakened buildings and infrastructure. Forecasts of the probability of aftershocks can therefore aid decision-making during earthquake response and recovery. Several countries issue authoritative aftershock forecasts. Most aftershock forecasts are based on simple statistical models that were first developed in the 1980s and remain the best available models. We review these statistical models, and the wide-ranging research to advance aftershock forecasting through better statistical, physical, and machine learning methods. Physics-based forecasts based on mainshock stress changes can sometimes match the statistical models in testing, but don’t yet outperform them. Physical models are also hampered by unsolved problems such as the mechanics of dynamic triggering and the influence of background conditions. Initial work on machine learning forecasts shows promise, and new machine learning earthquake catalogs provide an opportunity to advance all types of aftershock forecasts.

Annual Review of Earth and Planetary Sciences↗
USGS2025

The feasibility of using lidar-derived digital elevation models for gravity data reduction

Gravity data require submeter elevation accuracy for data processing, and differential global navigation satellite system (dGNSS) equipment is commonly used to acquire three-dimensional positional data to achieve such accuracy. However, lidar (light detection and ranging) data are commonly used to develop digital elevation models (DEMs) of Earth’s surface. Therefore, using elevations from lidar-derived DEMs for gravity-data acquisition and reduction may improve field efficiency and reduce cost. This study examines the feasibility of using DEMs for gravity-data reduction by comparing dGNSS elevation data from 435 gravity stations in Michigan, Wyoming, and Colorado with their respective DEM elevations. The results show that the average difference between DEM and dGNSS elevations is 13 centimeters (cm) and that 93 percent of those differences are less than 50 cm, even in areas with steep terrain. Because an elevation discrepancy of 50 cm corresponds to an error of roughly 0.1 milligals (mGal) in the simple Bouguer gravity anomaly, the results suggest that lidar-derived DEMs are a viable source for acquiring the elevation data needed to process gravity data, thus improving both the cost and efficiency of data collection for regional surveys where an accuracy of less than 1.0 mGal is desired.

Open-File Report↗
USGS2023

Cursed? Why one does not simply add new data sets to supervised geothermal machine learning models

Recent advances in machine learning (ML) identifying areas favorable to hydrothermal systems indicate that the resolution of feature data remains a subject of necessary improvement before ML can reliably produce better models. Herein, we consider the value of adding new features or replacing other, low-value features with new input features in existing ML pipelines. Our previous work identified stress and seismicity as having less value than the other feature types (i.e., heat flow, distance to faults, and distance to magmatic activity) for the 2008 USGS hydrothermal energy assessment; hence, a fundamental question regards if the addition of new but partially correlated features will improve resulting models for hydrothermal favorability. Therefore, we add new maps for shear strain rate and dilation strain rate to fit logistic regression and XGBoost models, resulting in new 7-feature models that are compared to the old 5-feature models. Because these new features share a degree of correlation with the original relatively uninformative stress and seismicity features, we also consider replacement of the two lower-value features with the two new features, creating new 5-feature models. Adding the new features improves the predictive skill of the new 7-feature model over that of the old 5-feature model; albeit, that improvement is not statistically significant because the new features are correlated with the old features and, consequently, the new features do not present considerable new information. However, the new 5-feature XGBoost model has a statistically significant increase in predictive skill for known positives over the old 5-feature model at p = 0.06. This improved performance is due to the lower-dimensional feature space of the former than that of the latter. In higher-dimensional feature space, relationships between features and the presence or absence of hydrothermal systems are harder to discern (i.e., the 7-feature model likely suffers from the “curse of dimensionality”).

Geothermal Resources Council Transactions↗
USGS2024

Constraining mean landslide occurrence rates for non-temporal landslide inventories using high-resolution elevation data

Constraining landslide occurrence rates can help to generate landslide hazard models that predict the spatial and temporal occurrence of landslides. However, most landslide inventories do not include any temporal data due to the difficulties of dating landslide deposits. Here we introduce a method for estimating the mean landslide occurrence rate of deep-seated rotational and translational slides derived solely from high-resolution (≤3 m) elevation data and globally available estimates of the diffusion coefficient for sediment flux. The method applies a linear diffusion model to the roughest landslide deposits until they reach a representative non-landslide roughness distribution. This estimates the time for a landslide deposit to be unrecognizable in high-resolution digital elevation data, which we term the mean lifetime of the landslide. Using the mean lifetime and number of landslides within an area of interest, we can estimate the mean occurrence rate of landslides over that domain. We validate this approach using a comprehensive temporal inventory of landslides in western Oregon created using age-roughness curves that are calibrated with high-resolution elevation data and radiocarbon data. We find good agreement between our diffusion method and the existing age-roughness-derived estimates, producing mean lifetimes of 4500 and 5200 years (4% difference), respectively. Hazard maps produced using the two methodologies generally agree, with the maximum differences in landslide probability reaching 0.1. Due to the relative abundance of high-resolution elevation data compared with age-dated landslides, our method could help constrain landslide occurrence rates in areas previously considered unfeasible.

Oregon↗
USGS1976

Ground magnetic and VLF studies at Midnite uranium mine, Stevens County, Washington

Ground magnetic and E-mode VLF measurements were made at more than 200 stations at the Midnite uranium Mine, Stevens County, Washington. The equipment was found to be light and foolproof, and each set of data was collected in one day, exclusive of surveying and marking the station points. The maps of contoured magnetic and VLF data resembled geologic maps of the area, showing many common features, as well as a few unique features on each map. Each showed some evidence of a small normal fault cutting the region; taken together, the evidence for this fault was judged to be much stronger, so that further structural details were inferred. Because of soil cover and inadequate outcrops, the presence of this fault had hitherto been unsuspected. A target area for further uranium mineralization north of the fault was indicated by this quick, inexpensive geophysical survey. The drilling program which Midnite Mine management had underway at that time had already outlined a mineralized zone in the indicated area.

Washington↗
USGS2024

Geophysical modeling of a possible blind geothermal system near Battle Mountain, NV

The northeastern portion of the Reese River basin in north-central Nevada is the focus of detailed geophysical and geological studies as part of the INGENIOUS project, which aims to identify new, commercially viable hidden geothermal systems in the Great Basin region of the western U.S. This location, herein referred to as Argenta Rise, occupies a broad (~15km wide) left-step between major range-front fault systems along the northwestern edge of the Shoshone Range and Argenta Rim, with numerous ENE-striking intra-basin faults presumably accommodating sinistral-normal oblique slip across the step-over. Four discrete regions have been identified within the study area that have favorable structural settings for hosting a blind hydrothermal system. However, with no definitive or extensive surface manifestations of an active hydrothermal system (e.g., geysers, steam vents, sinter, etc.), detailed geophysical studies are necessary to resolve subsurface geology and structure, and identify zones of enhanced structural complexity that may promote hydrothermal fluid flow. Hence, we collected high-resolution gravity, MT, and rock property data (density, magnetic susceptibility), and analyzed the recently acquired GeoDAWN aeromagnetic data to characterize potential geothermal resources in this region. Using the new geophysical datasets, we jointly modeled gravity and magnetic data along a series of intersecting 2D profiles that integrated information from recent, local-scale fault mapping. Rock property measurements performed on outcrops and hand samples throughout the study area constrained the models. The MT data were used to construct a 3D resistivity model that highlights the location of inferred alteration and fluids in the subsurface. Combined MT and potential field results reveal which structures may be most important for controlling hydrothermal fluid migration, as well as which geologic units may host hydrothermal fluids. Our gravity derived depth to basement surface coincides well with the base of shallow conductive anomalies, suggesting hydrothermal fluids may be confined to basin fill sediments and volcanics. This work supports our development of 3D geophysical and geologic models that are focused along the western flank of the northern Shoshone Range and aids the process of selecting sites for temperature gradient drilling.

Nevada↗
USGS2025

New constraints on location and timing of the Great Lakes tectonic zone, central Upper Peninsula, Michigan, USA

The Great Lakes tectonic zone (GLTZ) forms the boundary between the Wawa–Abitibi and Minnesota River Valley subprovinces within the Archean Superior Province. The GLTZ is concealed for all of its 1100 km length, except for a segment in the central Upper Peninsula of Michigan. There, it is exposed as a northwest-striking mylonite zone along a 11 km segment, extending to the onlap of Paleozoic rocks to the east. Farther east, its location has been unknown. Here, we use aeromagnetic and gravity data to develop interpretations of the expression of the GLTZ and to define its extent under cover. Aeromagnetic gradients over the mylonite zone are interpreted to be produced by structurally juxtaposed rocks with varying magnetizations. Gravity data show a regional gradient along the GLTZ, produced by the juxtaposition of a dense greenstone belt on the north against lower-density gneisses and granites on the south. The GLTZ is interpreted to extend ∼55 km under cover to the east. The GLTZ is terminated on the east by the buried eastern arm of the ca. 1100 Ma Midcontinent Rift. An undeformed granitic dike that cuts the mylonitic foliation produces a U–Pb apatite age of 2523 ± 33 Ma, implying no major post-Archean shearing occurred, and is at odds with previous interpretations of major Proterozoic reactivation. A granite intrusion in the Minnesota River Valley subprovince produces a Pb–Pb zircon age of 2606.9 ± 3.6/7.4 Ma. This suggests that magmatism related to the Sacred Heart orogeny, previously known in Minnesota, extended to Michigan.

Michigan↗
USGS2025

New paleomagnetic constraints on the eruption timing, stratigraphy, and post-emplacement deformation of the Picture Gorge Basalt within the Columbia River Basalt Group

The Picture Gorge Basalt (PGB) is part of the Miocene Columbia River Basalt Group (CRBG). The PGB, which outcrops in eastern Oregon, is considered coincident in time with the voluminous Grande Ronde Basalt. New radiometric ages have expanded the age‐range of the PGB, suggesting it began erupting prior to the Steens Basalt to its south and continued until after cessation of the Grande Ronde Basalt eruptions, an interval of 1.5 Ma. However, the existing paleomagnetism of the PGB implies this eruption timeline is an overestimate. To reconcile the radiometric and paleomagnetic timescales for the PGB, we conducted a paleomagnetic study on sections of the PGB to construct a detailed, high‐quality magnetostratigraphy. Our data indicate the stratigraphically lowest lava flows in the PGB are of reversed polarity, revealing a new paleomagnetic transition with the PGB and a reversed (R)–normal (N)–reversed (R) sequence. This suggests one of two timeline possibilities for PGB volcanism: (a) eruptions began and during through CRBG polarity chrons R0–N0–R1, penecontemporaneous with Steens Basalt, or (b) eruptions began and persisted during CRBG polarity chrons R1–N1–R2. Our work supports a longer interval of PGB volcanism than was suggested by previous paleomagnetic data but is at odds with the suggestion that PGB eruptions lasted through the entire main CRBG. We favor a scenario wherein PGB eruptions begin with R0 and continue into the R1 paleomagnetic interval. The paleomagnetic results also record a ∼18° vertical‐axis rotation of east‐central Oregon after ∼16 Ma with respect to stable North America.

Oregon↗
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↗
USGS2026

The Great Lakes Geologic Mapping Coalition—Working collaboratively to understand the geology of the Great Lakes Region

Introduction The Great Lakes Geologic Mapping Coalition (GLGMC), commonly referred to as the “Coalition,” is a partnership between the U.S. Geological Survey (USGS), the U.S. States of Illinois, Indiana, Michigan, Minnesota, New York, Ohio, Pennsylvania, and Wisconsin and the Canadian province of Ontario. The member States receive funding for geologic mapping work from the USGS National Cooperative Geologic Mapping Program (NCGMP), whereas Ontario participates as a nonfunded partner. The mission of the GLGMC is to produce three-dimensional (3D) geologic maps that depict unconsolidated sediments and near-surface bedrock in the Great Lakes region of North America. Geologic maps are the basis of most earth science investigations and help support resource exploration (energy, minerals, groundwater), natural hazard mitigation, infrastructure development, and land-use planning, all of which can be used to advance economic development and strengthen national security in the Great Lakes region. During the last few million years, the Great Lakes region has experienced repeated glacial advances and retreats, leaving behind extensive sediments, abundant natural resources, and widespread effects on the underlying bedrock geology (Swezey and others, 2022). Linked by shared histories of past glaciations, industrial agriculture, and legacy automotive, coal, steel, and manufacturing industries, the GLGMC member States collaborate to improve the understanding of the 3D distribution of the sediments overlying the region’s bedrock (fig. 1). Developing a comprehensive subsurface 3D framework of this glaciated terrain can provide earth science data to policymakers at all levels. These insights facilitate informed decisions on the exploration, use, and protection of vital resources, such as critical minerals, industrial materials, and aquifers, thereby supporting economic prosperity and the well-being of the citizens of this region. Since its inception in 1998, the Coalition has completed more than 100 geologic mapping projects across the Great Lakes region. Each project aims to deliver geologic maps, 3D datasets, and other information that improves understanding of the geology of the Great Lakes region, with an emphasis on economic and water resources. Key deliverables include 3D geologic maps and models typically portraying sediment thickness, often derived from top-of-bedrock and borehole data. These products are developed through a combination of fieldwork, subsurface modeling, and the collection and analysis of rock and sediment cores. To support Coalition goals, member States collaborate with scientists working on related STATEMAP, EDMAP, and FEDMAP projects. Coalition scientists also engage with Tribal Nations in the Great Lakes region to ensure that Tribal interests pertaining to Coalition work are addressed. Through this collaboration, the Coalition unites the efforts of State, Federal, and Tribal Nation stakeholders to advance geologic data production and enhance understanding of the geologic resources of the Great Lakes region.

Illinois, Indiana, Michigan, Minnesota, New York, ↗
USGS2025

Overview of The SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest earthquake sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress‐drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude‐spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress‐drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two‐week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite‐fault modeling, ground‐motion‐based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between‐method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress‐drop values.

California↗
USGS2025

Surficial geology and Quaternary tectonics of the Madison Valley and fault zone, Madison, Gallatin, and Beaverhead Counties, southwest Montana

The north-northwest-striking Madison fault is approximately 95 kilometers in length, lying at the confluence of the northeastern Basin and Range province and the Yellowstone tectonic parabola. The fault zone consists primarily of west-dipping normal faults that have east-dipping antithetic faults, which create the Madison Valley graben and several northeast-trending intrabasin faults. The Madison fault and associated sections discussed herein refer to the main west-dipping, range-bounding fault along the eastern side of the valley. Detailed geologic mapping (1:12,000 scale) of the entire fault zone and fault scarp profiling (total of 102 profiles) of the Madison fault reveal greater late Quaternary paleoseismic activity towards the south, including at least three paleoevents along the southern part of the fault that postdate Pinedale glaciation. Early to middle Holocene alluvial fans have vertical surface offsets that average between 2.0 and 3.0 meters and define the characteristic single-event surface offset. Pinedale lateral moraines have vertical surface offsets as great as 12.0 meters. Late Pleistocene to Holocene multiple-event fault scarps show little evidence of beveling, suggesting short seismic recurrence intervals and potential late Pleistocene and Holocene temporal clustering. Long-term average tectonic activity rates indicate slip rates ranging from 0.18–0.6 millimeters per year. Based on a comparison of fault-scarp height versus maximum slope angle of known regression lines developed from other paleoseismic investigations, the most recent event ranges from 5–1 ka. The northern section of the fault zone is defined by multiple normal faults, which detached the hanging walls of Laramide thrust faults within the Paleozoic and Mesozoic strata. This resulted in the partitioning of extension along multiple preexisting structures and less displacement along individual normal fault strands. Structural controls on lateral propagation of individual paleoevents involve the position of lateral ramps along preexisting Laramide contractional faults. This resulted in greater displacement within the larger basement-cored structures along the southern section of the fault zone, where extension is accommodated by one inferred principal basement-involved normal fault. Inferred east-northeast trending, intrabasin, normal faults within the southern half of the fault zone have no late Pleistocene displacement.

Idaho, Montana, Wyoming↗
USGS2025

Ultramafic float rocks at Jezero crater (Mars): Excavation of lower crustal rocks or mantle peridotites by impact cratering?

Based on observation and data from meteorites and in situ scientific missions, experiments as well as models, the Martian mantle is assumed to share some compositional and mineralogical affinity with the terrestrial mantle. However, there might be subtle differences like the Martian mantle being more ferroan. Yet, we do not have any direct analysis of a Martian mantle rock to confirm this assumption. NASA’s Perseverance rover found olivine-rich boulder-sized float rocks on the upper Jezero fan (Mars). These boulders have an ultramafic composition and their mineralogy is dominantly composed of Fo 73±3 olivine with high-Mg orthopyroxene, Cr-rich Ti-Fe oxides and minor plagioclase and high-Ca pyroxene. Microtextural and petrological analysis reveals that these minerals crystallized at equilibrium. In addition, these boulders are different from all the bedrocks analyzed by Perseverance along its traverse which are crustal igneous rocks and sediments. Comparing our data to Martian meteorites and available Mars bulk silicate models (BSM), we discuss that these boulders could represent primitive melts and/or lower crustal material, and we specifically hypothesize that they could be mantle peridotites. We propose that these putative mantle rocks could have been excavated by the succession of impacts from the shallow mantle or lower crust in the Isidis region where Jezero crater is located. These olivine-rich boulders could thereby constitute the first direct analysis of a Martian mantle rock.

Earth and Planetary Science Letters↗
USGS2025

Comparing magmatism and hydrothermal alteration using magnetic modelling and stable isotopes at the Clementine porphyry copper prospect, Montana, USA

Recent mapping of the underexplored Clementine prospect in southwestern Montana has revealed evidence of hydrothermal alteration and mineralized breccia vein gossan interpreted to represent the upper expression of a Cretaceous, sediment-hosted copper porphyry system. The prospect is at the nexus of several Cretaceous granites, including the pre-mineralization Butte Granite and Big Hole Canyon plutons and local pre- to post-mineralization granite intrusions. Here, 3D magnetic inversions and stable isotope data are used to evaluate the spatial and genetic relationship of mineralization and Cretaceous magmatism. Magnetic inversions reveal a zone of high magnetic susceptibility beneath the prospect that is related to an exposure of unaltered, likely post-mineralization granite. The granite appears to connect to the Big Hole Canyon pluton at depth, suggesting it is a late-stage differentiate of the cooling, pre-mineralization pluton that may have been one of several potential heat sources for the mineral system at Clementine. This is supported by carbon isotope values of graphite that indicate organic material in shale of the Cretaceous Kootenai Formation was subject to metamorphism. Preliminary light hydrogen isotope values of actinolite hydrothermal orbs suggest formation following magmatic outgassing during thermal advection from an underlying intrusion.

Montana↗