USGS ScienceSearch

SEARCH · USGS Science

Results for “Integrated Environmental Assessment and Management”

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,661 records · Page 21Linked to original sources

Synthesizing a twelve-year sediment trap time series of planktic foraminiferal flux in the Gulf of America (Mexico)

Sediment trap time series provide powerful frameworks for testing hypotheses about planktic foraminiferal assemblage composition, seasonality, and geochemical responses to environmental variability, all of which are central to improving paleoceanographic reconstructions. We present high resolution foraminiferal assemblage data from a long-running sediment trap (2008–2020) in the northern Gulf of America (Mexico). This study summarizes the species composition, seasonality, and size distribution of the fifteen most abundant species of planktic foraminifera, which account for 98% of total flux. Foraminiferal flux peaks in winter and reaches a minimum in summer, following the seasonal pattern of primary production in the northern Gulf. Winter assemblages are dominated by non-spinose taxa, whereas spinose taxa prevail during summer. Across nearly all species, average monthly test size covaries with temperature, independent of seasonal flux trends. Notably, Trilobatus sacculifer and Neogloboquadrina dutertrei show a significant decline in relative abundance and flux during 2017–2020, nearly disappearing from the assemblage.

Journal of Foraminiferal Research

Modeling legacy nitrogen transport under instantaneous, steady-state, and transient groundwater flow conditions

In hydrologic settings where groundwater discharge contributes substantially to surface waters, legacy nitrogen in groundwater can confound surface water nitrogen loads estimated exclusively from current terrestrial sources. Additionally, legacy nitrogen in groundwater can contribute to lagged responses to nitrogen management efforts. Some methods of estimating groundwater contributions to surface water nitrogen loads account for legacy nitrogen, while others do not. The resulting differences are rarely quantified. We used a numerical modeling framework to compare three methods of estimating time-varying annual groundwater nitrogen loads to surface water receptors on eastern Long Island, New York. The instantaneous load method used steady-state contributing areas and includes no temporal groundwater lag. The second method used numerical simulations of nitrogen loads under steady-state flow, which captures groundwater transport lags but omits the annual variability in transient hydrologic stresses. The third method numerically simulated both transient groundwater flow and nitrogen transport to explicitly capture the effects of legacy nitrogen in groundwater. Depending on antecedent nitrogen and hydrologic conditions, historical nitrogen loads estimated from the numerical simulations were sometimes similar (<10% difference) and other times substantially different (±100%) from the instantaneous load estimates. Additionally, simulated future surface water nitrogen loads responded asymptotically over several decades following reductions in terrestrial nitrogen sources, further highlighting the effect of groundwater transport lag times. The comparison of the three methods, quantification of historical interannual variability, and prediction of lagged responses to nitrogen source reductions provide important context for decision makers using estimated groundwater nitrogen loads to help evaluate nitrogen management efficacy.

New York

Utilizing downhole datasets for modelling the aeromagnetic signature of the Iron Creek Co-Cu deposit in the Idaho Cobalt Belt

The Idaho Cobalt Belt in east-central Idaho is host to some of the largest domestic Co resources, including the informal Iron Creek deposit. The two main ore zones of this deposit, the Iron Creek and the Ruby, are hosted in greenschist-grade interbedded argillite/siltstone and quartz-rich units of the Mesoproterozoic Apple Creek Formation of Lemhi Group. The primary ore mineral is cobalt-bearing pyrite, which occurs with pyrrhotite, chalcopyrite, and magnetite. This study integrates mineralogical and core-scale geophysical measurements with regional-scale aeromagnetic surveys. The high magnetite content within the Ruby zone produces elevated magnetic susceptibility, but the relatively limited spatial footprint of the ore zone would produce a small-scale anomaly that may be overlooked in regional surveys. The low magnetite content in the Iron Creek zone results in low magnetic susceptibility, creating a relatively low amplitude geophysical response. By characterizing the magnetic properties and mineralogy of these ore zones, this study enhances the interpretation of aeromagnetic data, enabling the identification of small or faint anomalies as potential Co targets. These findings improve can improve exploration strategies, both within the Idaho Cobalt Belt as well as for similar deposit types globally.

Idaho

Tracking diagenetic alteration of magnetic susceptibility in thrust ridge and slope basin sediments of the Cascadia margin (ODP Sites 1249 and 1252; IODP Site U1325)

We investigated sediment core records from the Cascadia Margin (Ocean Drilling Program Sites 1249 and 1252 at Hydrate Ridge; Integrated Ocean Drilling Program Site U1325 offshore Vancouver Island) using a Zr/Rb heavy mineral proxy from X-ray fluorescence (XRF) core scanning to identify intervals of primary detrital magnetic susceptibility (κ) and predict intervals where diagenesis caused magnetite dissolution by hydrogen sulfide. We also measured total sulfur (TS) content, grain size distributions, total organic carbon (TOC) content, and the magnetic mineral assemblage to further constrain the role of diagenesis on κ. Understanding how κ can be used to better characterize the varied effects of detrital and diagenetic signals in marine settings is important for understanding biogeochemical cycling and records of paleoenvironmental change. The upper 100 m of slope basin Site 1252 contains multiple intervals (> 90 m total) of decreased κ correlated with elevated TS content, consistent with dissolution of magnetite and precipitation of pyrite, iron monosulfides, and/or elemental sulfur. Similarly at the other slope basin site, Site U1325, κ is lower and TS is elevated in the interval between 24 and 51 mbsf, due to sulfide formation. At both slope basin sites, these low κ intervals correspond with high TOC, suggesting the possibility that organoclastic sulfate reduction (OSR) is likely a major driver of diagenetic alteration of κ at these sites. High TS:TOC ratios at Site U1325 suggests anaerobic oxidation of methane (AOM) during sulfate-methane transition zone (SMTZ) migration may have contributed to alteration of κ. In contrast, within the upper 90 m of Site 1249, a methane seep site at the summit of Hydrate Ridge, κ is almost entirely altered by diagenetic processes, with much of the low κ explained by a high degree of iron sulfide formation, while some intervals are affected by precipitation of magnetic iron sulfides that maintain or even increase κ. The presence of abundant methane seepage and gas hydrate as well as chemosynthetic seafloor fauna at this site, suggests that sulfide is released to the water column and AOM, rather than OSR, drives diagenetic alteration of κ at this site. Overall, the slope basin sites show episodic variation of κ that is influenced by TOC content, likely driven by changes in marine primary productivity and sedimentation rate, while the seep site shows consistently altered κ with lower TS content and no correlation with TOC. Methane seep environments likely experience loss of hydrogen sulfide to the water column and oxidation of hydrogen sulfide by seafloor seep fauna, which limits the amount of solid phase sulfur (pyrite, iron monosulfides, elemental sulfur) that can be precipitated within the sediments. In contrast, the migration of a buried SMTZ at slope basin sites results in enhanced sulfur precipitation within the sediments. This integrated magnetic and geochemical approach reveals the diagenetic production pathway and residence time of sulfide with the sediment column ultimately controls the style and degree of diagenetic loss of κ in marine sediments. This approach works best in environments with unaltered reference intervals/sites, consistent magnetic mineralogy, and clay-to-silt grain sizes.

Cascadia margin, Pacific Ocean

A crosswalk of the 2015 World Terrestrial Ecosystems to the International Union for the Conservation of Nature Global Ecosystem Typology Framework

To support ecosystem mapping and accounting applications, we aligned the 2015 U.S. Geological Survey/Esri/The Nature Conservancy-World Terrestrial Ecosystems (WTEs) with the International Union for Conservation of Nature Global Ecosystem Typology (GET) framework. This process, known as “crosswalking,” enabled the development of a global map of GET level 3 Ecosystem Functional Groups (EFGs) at a 250-meter spatial resolution. Crosswalking involved manually assigning 1,781 biogeographically stratified WTEs to their most probable EFG based on similarities in climate, terrain, vegetation, and geographic distribution. We compared attributes of the WTE dataset with summary characteristics of the EFGs. The resulting crosswalked global map of International Union for Conservation of Nature GET ecosystems is intended to be useful for standardizing ecosystem classification and reporting under frameworks such as the Kunming-Montreal Global Biodiversity Framework and the United Nations System of Environmental-Economic Accounting. We discuss key challenges in reconciling non-identical classifications, such as many-to-one relationships and variation in data quality.

Open-File Report

Widespread abyssal turbidites record megathrust earthquake-triggered landslides and coseismic deformation in the Cascadia subduction zone

Abyssal marine turbidites provide some of the longest and most spatially extensive records of subduction zone earthquake recurrence globally; however, correlation of these deposits over long distances and interpretation of synchronous emplacement requires both an understanding of the turbidite generating systems and precise dating. Here, we present an integrated suite of high-resolution bathymetry, subbottom profiles, and sediment cores from combined autonomous underwater vehicle, remotely operated vehicle, and ship-based studies at a key paleoseismic site in the southern Cascadia subduction zone. We demonstrate how widespread, earthquake-triggered landslides on the lower slope deposit discrete, proximal mass transport deposits (MTDs) that grade offshore into complex, interfingered abyssal turbidites, which correspond to records of megathrust earthquake history. We propose accretion and oversteepening of thrust folds on the lower slope both preconditions the slope to fail and provides a perpetual source of unstable material to fail during every earthquake cycle. Furthermore, we suggest the periodic and pervasive landsliding indicates coseismic deformation of the outer accretionary wedge during megathrust rupture.

California, Oregon

The influence of nodule versus crust morphology on the composition of seamount-hosted ferromanganese minerals

Ferromanganese mineral precipitation in the global oceans is ubiquitous, occurring in the form of both crusts and nodules at a broad range of depths and seafloor terrains. Although ferromanganese crusts and nodules are both composed of ferromanganese minerals, mineralogy and mean element concentrations compiled for regional crust versus nodule occurrences differ. Notably, most published compilations compare nodules from abyssal plain sediments to crusts forming on seamounts, and do not address the question as to whether mineral morphology may affect the composition of crusts and nodules in addition to environmental factors. To address this, we have compiled a dataset of co-located hydrogenetic (seawater sourced) open-ocean ferromanganese crusts and seamount-hosted nodules. The preliminary data presented here support our hypothesis that co-located crusts and nodules exhibit similar compositions, both in terms of their mineralogy and mean element concentrations, as well as overall variance in concentrations for any given element. These data indicate that local processes influence sample-to-sample variation in both crusts and nodules, in contrast to broader ferromanganese variation reported between open-ocean seamounts- hosted crusts and abyssal plain-hosted nodules.

Conference Paper

Foundational uncertainties in terminal Ediacaran chronostratigraphy revealed by high-precision zircon U-Pb geochronology of the Nama Group, Namibia

The Nama Group of southern Namibia and northwestern South Africa hosts the best-dated mixed carbonate-siliciclastic foreland basin succession of the terminal Ediacaran [ca. 551 million years (Ma) ago to <538 Ma] and is key for resolving the chronology of early metazoan evolution. Numerous silicified volcanic tuff interbeds are present, but differing interpretations regarding the fidelity of their ages lead to different regional stratigraphic correlations, especially for the Urusis Formation of the Schwarzrand Subgroup. An expanded record of the Urusis Formation is found in the Swartpunt area of southern Namibia, which has yielded an important metazoan biota. But the succession in this area is preserved as a series of thrusts at the leading edge of the Gariep orogenic belt and zircon U-Pb data show systematic age repetition. We use regional stratigraphic and structural mapping, integrated with carbonate carbon isotope (δ 13 C carb ) chemostratigraphy and high-precision radioisotope U-Pb zircon geochronology from outcrop and recently acquired drill core to develop a temporally calibrated basin-wide depositional model. This integrated dataset either reflects complex zircon reworking, inheritance, or potential analytical biases (Scenario 1) or the presence of a Gariep-related cryptic décollement within the Spitskop Member that has resulted in stratigraphic repetition (Scenario 2). We investigate the evidence for and against both scenarios and consider their implications for stratigraphic and δ 13 C carb correlations between the Swartpunt area and coeval exposures along the Orange River border with South Africa. Given that these issues are in an area that hosts numerous silicified ash beds and extensive exposure, an inability to confidently discount either scenario highlights a level of compounding uncertainty in zircon U-Pb geochronology that must be considered when attempting to build global chronostratigraphic frameworks. Scenario 1 implies that some of the weighted mean ages and Bayesian eruption ages from the Swartpunt area may be >1 Myr older than the depositional age of their respective ash beds when assuming existing stratigraphic correlations. If this scenario is preferred, then a cautious approach would be to consider all weighted mean zircon U-Pb ages from ash beds to reflect maximum depositional ages. Both scenarios support deposition of the Huns Member >540 Ma in the Swartpunt area if the oldest weighted mean age reported here represents a near-depositional age, which has significant implications for the temporal calibration of important terminal Ediacaran ichnofossil assemblages and future cyclostratigraphic studies. Stratigraphic correlations common to both scenarios allow us to temporally calibrate a basin evolution model for the Nama Group. Temporal trends in initial hafnium isotope (εHf) compositions of zircon grains from ash beds throughout the succession may support progressive crustal thickening associated with underplating of the Damara orogenic belt along the northern periphery of the Kalahari craton from ca. 547 Ma to ca. 538 Ma. The compilation of new and published zircon U-Pb ages may also imply that the locus of carbonate platform development migrated from north to south (present co-ordinates), tracking the migration of foredeep subsidence.

Earth-Science Reviews

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi

Testing characteristic magnitude distributions in modern PSHA models

The characteristic magnitude distribution hypothesis predicts a higher rate of large earthquakes than a Gutenberg–Richter extrapolation of the small‐earthquake rate would imply. Characteristic magnitude distributions have been commonly applied to faults in probabilistic seismic hazard analysis (PSHA), and in modern models they can emerge from the way short‐term seismicity constraints are combined with long‐term geologic and geodetic constraints. We test the characteristic magnitude distribution hypothesis by comparing the fault‐based magnitude distributions from the 2023 update to the National Seismic Hazard Model (NSHM23) in the Western United States with observed seismicity over the past 93 yr. We find that observed magnitude distributions fall outside the model‐predicted confidence bounds in regions where NSHM23 produces characteristic magnitude distributions: in these regions, the model predicts higher rates of large earthquakes than are observed. An analysis of the earlier California model (Uniform California Earthquake Rupture Forecast, version 3) also reveals discrepancies between the modeled and observed magnitude distributions. In addition, we find that observed magnitude distributions near modeled faults are not significantly different from those in background regions. These results challenge the prevalence of characteristic magnitude distributions in fault‐based seismic hazard models and call for a reassessment of how disparate data sets are integrated in PSHA.

western United States

Self-guided decision support groundwater modelling with Python

The GMDSI tutorial notebooks repository provides learners with a comprehensive set of tutorials for self-guided training on decision-support groundwater modelling using Python-based tools. Although targeted at groundwater modelling, they are based around model-agnostic tools and readily transferable to other environmental modelling workflows. The tutorials are divided into three parts. The first covers fundamental theoretical concepts. These are intended as background reading for reference on an as-needed basis. Tutorials in the second part introduce learners to some of the core concepts parameter estimation in a groundwater modelling context, as well as providing a gentle introduction to the PEST, PEST++ and pyEMU software. Lastly, the third part demonstrates how to implement highly-parameterized applied decision-support modelling workflows. The tutorials aim to provide examples of both “how to use” the software as well as “how to think” about using the software. A key advantage to using notebooks in this context is that the workflows described run the same code as practitioners would run on a large-scale real- world application. Using a small synthetic model facilitates rapid progression through the workflow.

Journal of Open Source Education

Classification of lakebed geologic substrate in autonomously collected benthic imagery using machine learning

Mapping benthic habitats with bathymetric, acoustic, and spectral data requires georeferenced ground-truth information about habitat types and characteristics. New technologies like autonomous underwater vehicles (AUVs) collect tens of thousands of images per mission making image-based ground truthing particularly attractive. Two types of machine learning (ML) models, random forest (RF) and deep neural network (DNN), were tested to determine whether ML models could serve as an accurate substitute for manual classification of AUV images for substrate type interpretation. RF models were trained to predict substrate class as a function of texture, edge, and intensity metrics (i.e., features) calculated for each image. Models were tested using a manually classified image dataset with 9-, 6-, and 2-class schemes based on the Coastal and Marine Ecological Classification Standard (CMECS). Results suggest that both RF and DNN models achieve comparable accuracies, with the 9-class models being least accurate (~73–78%) and the 2-class models being the most accurate (~95–96%). However, the DNN models were more efficient to train and apply because they did not require feature estimation before training or classification. Integrating ML models into benthic habitat mapping process can improve our ability to efficiently and accurately ground-truth large areas of benthic habitat using AUV or similar images.

Michigan, Wisconsin

Effect of mineral deposit data on predictions from the three-part approach to quantitative mineral resource assessment—A study of 16 previous U.S. Geological Survey assessments

The three-part approach to quantitative mineral resource assessment requires information about the properties of undiscovered mineral deposits in an assessment area. These properties are unknown, so the properties of discovered mineral deposits of the same mineral deposit type are used instead. In the three-part approach, these discovered mineral deposits come from around the world, and their properties constitute the pooled data for that mineral deposit type. Alternatively, these discovered mineral deposits could come from the assessment area, and their properties constitute the tract data for that mineral deposit type. Tract data may be more representative of the undiscovered mineral deposits in the assessment area than the pooled data. The goal of this study was to determine whether resource predictions using pooled data are equivalent to resource predictions using tract data. To this end, 16 previous U.S Geological Survey assessments were studied. For each assessment, resources were predicted for one undiscovered mineral deposit in the assessment area. One set of predictions used pooled data, and another used tract data. The two sets of predictions were compared with an equivalence test, using the six assessment statistics that are commonly reported for mineral resource assessments. Practical equivalence is the condition that two corresponding assessment statistics are within a factor of 1.5 of one another. For each of 2 assessments, all 6 assessment statistics were practically equivalent. For both assessments, the assessment statistics from the pooled data, relative to the corresponding assessment statistics from the tract data, ranged from 1.30 times smaller to 1.03 times larger. For each of 14 assessments, 1 or more of the 6 assessment statistics were not practically equivalent. The assessment statistics from the pooled data, relative to the corresponding assessment statistics from the tract data, ranged from 26.6 times smaller to 5.53 times larger. The use of pooled data has been a standard procedure in the three-part approach since at least 1986. The 16 assessments in this study are not a representative sample of those prior assessments that used pooled data. So, it is inappropriate to use the study results to infer whether pooled data affected the resource predictions for those prior assessments.

Scientific Investigations Report

High-resolution landscape-scale maps of nutrients, phytoplankton, and related water-quality constituents in Suisun Bay and the Sacramento–San Joaquin Delta, California, during 2018—Results from high-resolution underway surveys

We examined the abundance and distribution of nutrients and phytoplankton in the tidal aquatic environments of the Sacramento–San Joaquin Delta (Delta) and Suisun, Grizzly, and Honker Bays in California by completing three spatial surveys that used continuous underway high frequency (1 hertz) measurements and sampling onboard a high-speed boat. Surveys were conducted in May, July, and October 2018. Water was sampled continuously during surveys to simultaneously collect information about the concentration and spatial distribution of all major nutrient forms together with information about the major classes of phytoplankton and associated standard field measurements of water quality, such as dissolved oxygen, dissolved organic matter, pH, salinity, specific conductance, turbidity, and water temperature. The results show a greater than 50-fold variation in nutrient and phytoplankton concentrations across space and time, providing evidence of the dynamic environmental processes that shape the ways nutrients interact with and affect Delta aquatic habitats.

Data Report

Using a time-of-travel sampling approach to quantify per- and polyfluoroalkyl substances (PFAS) stream loading and source inputs in a mixed-source, urban catchment

Understanding per- and polyfluoroalkyl substances (PFAS) mass distribution in surface and groundwater systems can support source prioritization, load reduction, and water management. Thirteen sites within an urban catchment were sampled utilizing a time-of-travel sampling approach to minimize the influence of subdaily fluctuations in mass from PFAS point sources and to quantify PFAS and ancillary chemical loads from various PFAS sources. A larger increase in perfluoroalkyl sulfonate (PFSA) loads (8 to 11 μg/s, up to 618%) than in perfluoroalkyl carboxylate (PFCA) loads (no change to 3.4 μg/s, up to 122%) was observed at sites below tributaries influenced by military bases with known groundwater discharge. Point discharges from two sewage treatment plants (STPs) resulted in increases in PFCA and PFSA loads that were similar (6 and 10 μg/s respectively) below the first STP and greater for PFCA compared to PFSA loads (23 and 13 μg/s respectively) below the second STP. Overall, percent increases in total PFAS load ranged from 20 to 277% for military base inputs and 44 to 77% for STP inputs. A focus catchment that represents only 14% (76.9 km 2 ) of the drainage area at the most downstream site (544 km 2 ) accounted for about 70% of PFSA and 40% of PFCA loads observed at the most downstream site. Results show that by using a time-of-travel sampling approach in mixed, urban settings with several PFAS sources, it is possible to quantify stream loads from individual PFAS sources, thereby improving source attribution and providing actionable data for water-resource managers.

ACS ES&T Water

A framework for understanding the effects of subsurface agricultural drainage on downstream flows

Understanding controls on streamflow volume and magnitude is important to water resource management applications, such as critical water and transportation structure design and floodplain mapping. Changes in land use and agricultural practices, such as subsurface agricultural drainage, may be contributing to changes in streamflow characteristics. Subsurface agricultural drainage, also known as tile drainage, is the practice of installing drains in the subsurface of agricultural fields to improve productivity. Because of the complex interactions between subsurface drainage systems, precipitation, local soil conditions, and land management practices, it is difficult to determine how subsurface agricultural drainage affects downstream flow. Previously developed subsurface agricultural drainage conceptual models under dry, saturated, and winter conditions are summarized, and current literature on the effects of subsurface agricultural drainage on downstream flows, focusing on peak flow, non-event flow, and total flow to develop frameworks for discussing these systems is compiled. The effects that subsurface drainage has on hydrologic systems are expected to vary by site and are seasonally based on system design, soil type, moisture conditions, precipitation characteristics, and land conditions. Subsurface drainage can affect the magnitude of peak flow by converting surface runoff from a storm event to subsurface runoff. By increasing hydrologic connectivity of a catchment, subsurface drainage can increase non-event flow or the flow between two storm events, typically dependent on lateral flow through the subsurface and groundwater. Theoretically, by diverting water from groundwater recharge or by reducing water available for evapotranspiration, subsurface drainage may increase the total volume of flow. Precipitation changes may increase infiltration, excess overland flow, and flood risk regardless of the presence or absence of subsurface drainage.

Illinois, Iowa, Michigan, Minnesota, Missouri, Mon

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

Neotectonic mapping of Puerto Rico

Puerto Rico is part of the Puerto Rico-Virgin Islands microplate, along the Caribbean–North American plate boundary between the Puerto Rico trench subduction zone and the Muertos Trough incipient subduction zone. Despite recent seismicity and geodetically constrained deformation of ~3 mm/yr of left-lateral shear across the island, Quaternary fault locations remain largely uncertain. Preservation of recent faulting in the landscape is masked by distributed faulting, high weathering rates of the tropical climate, steep topography, frequent landsliding, and extensive agriculture and urbanization along coastlines and river valleys. We present remote neotectonic mapping of historical aerial imagery from the 1930s-60s and <1-m lidar topography, integrated with field observations, to create an updated active fault map of Puerto Rico. We focus on faults that offset younger geomorphic surfaces, with ages estimated by geologic mapping, OSL and radiocarbon dating, and morphologic interpretations. We present new evidence for Quaternary activity on seven faults, including the South Lajas, Salinas, Punta Montalva, Great Southern Puerto Rico, Cerro Goden, Parguera, and San Marcos faults. We find that active faulting occurs preferentially along the southern and western coasts and does not always spatially coincide with preexisting bedrock faults, possibly reflecting partial reactivation of older faults in the modern strain field.

Puerto Rico