Search USGS⌕ Search

SEARCH · Search USGS

Results for “Professional Paper”

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.

188 records · Page 6Linked to original sources
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

Development of USGS NSHMs: Do small changes in hazard imply small changes in risk?

One of the flagship products from the U.S. Geological Survey (USGS) is the National Seismic Hazard Model (NSHM). Since 1976, the NSHM has been periodically updated to reflect newly published earthquake science and provide probabilistic estimates of seismic hazard for the United States. During each update cycle, alternative models are deliberated, analyzed, and documented through logic trees and their corresponding logic tree branch weights. For example, the decision to modify a logic tree branch weight may be influenced by sensitivity analyses of the logic tree branches in their effects on the mean hazard. However, do small changes in traditional measures of hazard imply small changes in risk? In this study, we make use of two update cycles of the USGS NSHMs and a National Bridge Inventory (NBI) from the Federal Highway Administration (FHWA) to explore the preceding question. Specifically, we first identify geographic locations in the conterminous United States in which the change in hazard from one cycle to another is relatively small. Next, we model the seismic risk to highway bridges for these locations and for each update cycle, while simultaneously distinguishing low hazard environments from high hazard environments. These data enable quantitative analysis of how much changes in risk can be expected given small changes in hazard, investigating the importance of risk models in decision-making during development of the USGS NSHMs.

Conference Paper↗
USGS2025

Evidence for offset of Cretaceous plutons by the Tintina fault in eastern Alaska: Implications for regional metallogeny

Cretaceous magmatism in eastern interior Alaska is voluminous, but temporally and spatially diverse – suggestive of varying sources and drivers. More than 150 new U-Pb zircon and more than 500 geochemical analyses of Cretaceous plutonic units allow for the grouping of distinct plutonic suites. Magmatism was continuous from 120-66 Ma but can be grouped into temporally distinct pulses from ca. 115-100 Ma, 100-90 Ma, and 75-66 Ma. Geochemical diversity occurs during each pulse, further distinguishing multiple suites. Diverse metallogenic epochs are strongly correlated to pluton chemistry. Mineralization is largely absent prior to 108 Ma. From 108-100 Ma, plutonism is coeval with sparse, but notable Au-quartz veins with variable Bi, As, W, and Mo. From 100-90 Ma, intrusion-related mineralization zones from Au-Cu(-Bi) and U-Th in the northwest to central Au- Bi-As-Te(-W), and Mo-W to the southeast. Porphyry style Cu-Mo(-Au) occurrences occur with the latest Cretaceous plutons emplaced from 75-66 Ma . Restoration of ~450 km of dextral movement on the Tintina fault and comparison of metallogenic and geochemical characteristics of Alaska plutons suggest 100-90 Ma plutons may be the continuation of the metallogenically significant Tombstone, Mayo, and Tungsten suites from the Yukon.

Alaska↗
USGS2026

Birds, breakpoints, and baselines: How citizen science data can reveal ecological boundaries in Kenya’s Upper Tana watershed

Tropical watersheds are increasingly threatened by climate change, land-use conversion, and resource extraction, yet conventional biodiversity monitoring in these systems is often spatially and temporally limited. Citizen science offers a complementary approach, enabling biodiversity data collection over large areas that can supplement professional scientific surveys. We analyzed 10 years (2012–2022) of Upper Tana Watershed bird data from the Kenya Bird Map project, covering 114 pentads (9 × 9 km) within a 17,000 km 2 watershed to assess patterns of bird community composition and distribution across this watershed and to also evaluate the effects of environmental variables, seasonality, and sampling effort to help inform improvements in future citizen science projects. Citizen (or community) scientists recorded 575 species (>50% of Kenya’s total avifauna) in 74 families. Asymptotic species accumulation indicates that most probable species present in the watershed were detected. Threshold indicator taxa analysis revealed distinct ecological boundaries along elevation (∼1,500 m), precipitation (∼1,100 mm), and mean temperature (∼19°C) gradients, corresponding to a turnover from xeric savanna to mesic montane forest assemblages. Notably, bird communities showed little seasonal differentiation between wet and dry periods, consistent with dominance by resident year-round species. Data limitations including uneven survey distribution and frequency, absence of abundance metrics, and coarse representation of local environmental conditions that likely reduced our ability to detect fine-scale species–habitat relationships. Addressing these gaps through spatially balanced sampling at greater resolution, greater survey frequency in underrepresented areas, and improved capture of habitat metrics could strengthen the use of citizen-science bird data for watershed bioassessment. Our findings demonstrate that structured citizen-science initiatives can identify ecological boundaries and inform adaptive management of tropical socio-ecological systems under rapid environmental change.

Upper Tana watershed↗
USGS2025

Avak Creek oil occurrence, North Slope, Alaska: Newly discovered oil seep on Native lands, near village of Utqiagvik

An unknown occurrence of oil was detected near Avak Creek on Native lands on the North Slope of Alaska. Determining the source of oil was imperative for allowing stakeholders (Federal, State, and local government agencies and the landowner, an Alaska Native corporation) to make timely and informed decisions and mount a mitigation response, if required. The regional and local geological framework of the Avak Creek site was constructed using seismic surveys, well data, and basin modeling results, to identify local petroleum systems, map structural geometry and faults, define source rock thermal maturity distributions, and infer likely oil-migration pathways. Molecular hydrocarbon fingerprints (biomarkers, diamondoids, compound-specific isotopes) of the oil were compared to those of local and regional oil seeps, exploration well tests, and produced oils. Biomarker acid distributions characterized the history and extent of petroleum biodegradation. Integrating subsurface and geochemical parameters, the oil is interpreted to be a natural seep generated locally, predominantly from the Brookian Lower Cretaceous Hue Shale/gamma-ray zone, rather than an anthropogenic source of pollution. Results highlight sophisticated analytical technologies used to characterize complex, compositionally altered hydrocarbons. Results also advance our understanding of Brookian source rock distribution, subsurface petroleum migration pathways, and Arctic Alaska petroleum systems.

Alaska↗
USGS2026

In-situ geochronology using LA-ICP-MS/MS: Application of the Lu-Hf system in carbonate, apatite and fluorite

Geochronology is a fundamental tool for reconstructing earth history and constraining the timing of ore deposit formation. Traditionally, U-rich minerals like zircon, titanite, monazite and apatite have been analysed to determine the timing of mineralization and alteration events. However, not all deposits or alteration assemblages contain these minerals, and the U-Pb system can be disturbed due to open system behaviour. We present a methodology for obtaining in-situ Lu-Hf dates from various minerals to demonstrate the potential for dating ore deposits that lack U-rich minerals. Carbonate, apatite and fluorite minerals from different tectonic settings and geological environments were analysed. We acquired Lu-Hf and REE (rare earth element) data from carbonate minerals from a carbonatite orebody and dike samples from the Mountain Pass REE carbonatite deposit in California, USA. We also tested a ‘double-dating’ method, where U-Th-Pb dates and Lu-Hf dates were obtained from the same analysis of apatite from rocks in the Michigan’s Great Lakes Tectonic Zone. Preliminary results for Lu-Hf dating of fluorite from carbonatite rocks in the Lemitar Mountains of New Mexico are also included. The analytical methodology shows promising results in each of these cases where there is sufficient parent Lu and U and measurable daughter Hf and Pb respectively.

Conference Paper↗
USGS2025

Tracing fluid evolution during batholith and deposit formation: In-situ and bulk rock Mo isotopes from the Questa porphyry Mo deposit, New Mexico, USA

The Oligocene Questa porphyry Mo deposit, located in New Mexico, USA, is a well-characterized system that formed over 400 kyr, punctuating the 8 Myr formation of the Questa-Latir batholith. Detailed sampling of the intrusions and minerals associated with magmatichydrothermal mineralization, as well as bookending barren intrusions, allows for examination of the Mo isotope record of the deposit and its host batholith. Through this established chronologic framework, such data can yield insights into the metallogenesis of deposit formation and long-term magma genesis evolution for upper crustal magmatism. The Mo isotopic compositions of barren intrusions yield a younging trend towards lighter δ 98 Mo. However, intrusions associated with Mo deposit formation yield heavier δ 98 Mo relative to barren magmatism, by up to 2.0‰. This indicates that they were affected by an isotopically heavy and Mo-rich fluid. The in-situ data from molybdenite covers the range observed in mineralized bulk rocks. However, individual samples consistently yield Mo isotopic variations of 0.4‰ to 0.5‰. The consistency of the δ 98 Mo variation in in-situ analyses of molybdenite within hand samples and individual mineral grains is likely the result of MoS 2 crystallization in a closed system and reflects an almost complete precipitation of Mo.

New Mexico↗
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

Simulated ground-motion records for the seismic assessment of monumental masonry structures

Earthquakes are natural disasters that can cause widespread devastation and loss of life. Simulated ground-motion records can be useful in regions with limited seismic stations or a history of damaging but infrequent earthquakes. This is especially true in areas with a high concentration of heritage masonry structures, which are especially susceptible to damage, as simulated records can be crucial in predicting their seismic response and protecting these buildings from seismic damage. Despite the importance of simulated earthquakes, few studies have investigated how effective they are compared to real earthquakes when assessing the structural response of heritage buildings. To address this knowledge gap, we employ two different simulation methods of the Mw 6.2 Faial earthquake, which occurred on July 9th, 1998, in the Azores, to replicate the recorded time-series at four available stations within an epicentral distance of 150 km. The study has two objectives: first, to validate the simulated records of the 1998 Faial earthquake using alternative stochastic ground-motion simulation approaches, and second, to determine how these approaches impact the seismic assessment of historic masonry structures. To accomplish these objectives, this study uses real and simulated ground-motion datasets to conduct non-linear response history analyses of the São Francisco Church, a monumental structure in Horta that sustained damage during the Faial earthquake. The results show that both simulation approaches yield structural responses similar to the observed records.

Conference Paper↗
USGS2025

Detecting hidden sedimentary geothermal systems in the Upper Colorado River Basin

Geothermal resources exist in sedimentary rock where circulation of water for efficient extraction or storage of heat is possible. Except in rare instances where hot water is expressed at the land surface, sedimentary geothermal resources are hidden, so the identification of these systems is optimally accomplished using predictive subsurface modeling. An integrated approach using detailed paleogeographic interpretations, subsurface geologic mapping, and numerical modeling has produced regional geologic and temperature models for the Upper Colorado River Basin, a large watershed in central North America that contains many sedimentary basins. These models identify areas of hidden sedimentary geothermal resource potential in low temperature (<90°C), moderate temperature (90–150°C), and high temperature (>150°C) fairways across the study area. These models incorporate maps of key horizons in outcrop and the subsurface to create a robust structural framework that can be used to target favorable geology for natural or engineered permeability. This framework is populated with lithologies derived from detailed palaeogeographical maps and over 40,000 bottom hole temperature (BHT) values were used to create a calibrated three-dimensional (3D) temperature model across the region. The resulting maps serve as a regional sedimentary geothermal play fairway screening tool for evaluating different grades of sedimentary geothermal resources and for identifying areas of interest where more detailed, prospect-scale studies can be undertaken.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2025

The role of alkali bicarbonate-sulfate brines in the genesis of carbonatite REE resources at the Bear Lodge Alkaline Complex, Wyoming

Rare-earth element (REE) resources in the Bear Lodge Alkaline Complex, Wyoming, are hosted in a variably leached carbonatite dike swarm spatially related to bodies of diatreme breccia. This study examines fluid inclusions in carbonatite dikes, peripheral fluorite breccias, and smoky quartz veins to reconstruct the physiochemical conditions of REE mineralization. Results reveal a multi-stage evolution of carbonatitic fluids: 1) high-temperature (330–432°C), CO 2 -rich magmatic fluids, 2) REE-enriched alkali bicarbonate-sulfate brines, and 3) low-temperature meteoric (117–182°C) or diluted magmatic fluids. Multiphase inclusions (stage-2) contain burbankite, nahcolite, and alkali sulfate daughter crystals, linking the brines to early burbankite [(Na,Ca) 3 (Sr,Ba,REE) 3 (CO 3 ) 5 ] mineralization in carbonatite. REE mobility has been linked to the stability of aqueous complexes with ligands such as Cl-, SO 4 2- , and HCO 3 - . However, light REEs were not initially transported significant distances, rather they crystallized as burbankite from brines or brine-melts associated with the carbonatite magma. Crystallization of burbankite removed Na from the system, leaving residual fluids with high K/Na ratios and promoting potassic fenitization throughout the complex. Fractionation of light and heavy REEs and involvement of low temperature (<300°C), high K/Na bicarbonate brines may have facilitated transport and formation of peripheral occurrences that are more enriched in heavy REEs at Bear Lodge.

Wyoming↗
USGS2024

Recent applications of the USGS National Crustal Model for Seismic Hazard Studies

The U.S. Geological Survey is developing the National Crustal Model (NCM) for seismic hazard studies to facilitate modeling site, path, and source components of seismic hazard across the conterminous United States. The NCM is composed of a 1km grid of geophysical profiles, extending from the Earth’s surface into the upper mantle. It is constructed from a threedimensional (3D) geologic framework and geophysical rules that use (1) a petrologic and mineral physics database; (2) a 3D temperature model; and (3) a calibrated rock type- and age-dependent porosity model. Parameters needed to estimate site response for existing ground motion models (GMMs), including the time-averaged velocity in the upper 30 meters (VS30), the depths to 1.0 and 2.5 km/s shear-wave velocity (Z1.0 and Z2.5), and sediment thickness, can be computed from the NCM. As GMMs continue to improve in the future, other metrics could also be extracted or derived from the NCM, such as fundamental period, site attenuation (ko), a fully frequency-dependent site response function, or 3D geophysical volumes for wavefield simulations. Application of the NCM may also benefit other aspects of seismic hazard analysis, including better accounting for path-dependent attenuation and geometric spreading, more accurate estimation of earthquake source properties such as hypocentral location and stress drop, and calculation of crustal strength profiles that inform estimates of the base of seismicity.

conterminous United States↗
USGS2025

Constraining large magnitude event source and path effects using ground motion simulations

The purpose of this study is to use ground motion simulations to investigate ways in which source and path effects for large magnitude events can be represented in non-ergodic GMMs. While we initially developed computation techniques using CyberShake simulations, the range of magnitudes and source-site combinations is not adequate to replicate what is observed empirically. We therefore designed a new ground motion simulation study, which includes earthquakes with a large range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of simulations (M4-M7), we then develop a non-ergodic GMM with the simulation data. We find that the within-site residuals are dominated by the radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within-site residuals among multiple source realizations. Finally, after removing the source effects from the within-site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site with the large events, while the second approach considers all small events on the fault plane. The results indicate that the path effects of large events cannot be satisfactorily approximated with that of small events using either approach.

Conference Paper↗
USGS2024

Groundwater-Surface water interactions research: Past trends and future directions

Interactions between groundwater and surface water sustain groundwater-dependent ecosystems and regulate river temperature and biogeochemical cycles, amongst many other processes. These interactions occur in freshwater environments including rivers, springs, lakes, and wetlands, and in coastal environments via tidal pumping, submarine groundwater discharge, and seawater intrusion. Here, we explore groundwater-surface water interactions research using bibliometric analyses of titles, abstracts, and keywords from 20,275 journal papers published between 1970 and 2023 extracted from Scopus. Analyses show that research into groundwater-surface water interactions is highly multi-disciplinary, with growing contributions from the social and biological sciences. The number of groundwater-surface water interactions papers is rapidly increasing with over 1200 papers published per year since 2020. Drawing on our data-driven approach and expert knowledge, we synthesise current research trends and identify critical future research directions. Despite the thousands of papers on groundwater-surface water interactions, important processes are still difficult to quantify or predict at meaningful spatial scales to inform water-resources management. We see benefits in future groundwater-surface water interactions research focusing on: (1) using new technologies including internet-of-things-based sensors, uncrewed vehicles, and remote-sensing approaches for data collection to inform groundwater-surface water interactions at large scales, (2) seeking approaches to upscale site-specific findings to better inform management, and (3) continuing the movement towards multi-disciplinary investigations to better inform the understanding of groundwater-surface water interactions and processes that will enable better management outcomes.

Journal of Hydrology↗
USGS2025

Quantitative mineral resource assessment of lithium pegmatite deposits in the Appalachian Orogen, USA

Lithium is classified as a U.S. critical mineral commodity, and its demand is projected to drastically increase through 2040, driven by electric vehicle production and energy storage applications (IEA 2021).Most global lithium production is not in the United States increasing vulnerability to a supply disruption. The U.S. Geological Survey is actively assessing domestic lithium deposits including lithium-bearing pegmatites in the Appalachian orogen. Permissive tracts for lithium pegmatite deposits were delineated by integrating lithological, tectonic, geochemical, geophysical, and mineral occurrence data. The geospatial data and permissive tracts were used to estimate the number of undiscovered lithium pegmatite deposits. Estimates were then integrated into probabilistic simulations along with a new global lithium pegmatite grade and tonnage dataset to quantify potential contained undiscovered lithium resources. An economic filter was used to estimate the amount of potentially recoverable undiscovered resources. Preliminary computations for the northern Appalachians, including application of the economic filter to the median recoverable contained resource, yields 900,000 metric tons of Li 2 O that correspond to enough Li 2 O to replace 127 years of import reliance at the current rate (7,100 t Li 2 O/yr; USGS, 2025). For the southern Appalachians, preliminary computations yielded 1,430,000 metric tons of Li 2 O, which corresponds to 201 years of import reliance.

Alabama, Connecticut, Delaware, Georgia, Maine, Ma↗
USGS2026

Stepovers and beyond: Structural control of The Geysers geothermal system and the broader Clear Lake region

Fault geometry exerts a first-order control on geothermal systems by governing stress localization, fracture development, and permeability, yet in complex fault networks or broader shear zones, the relative influence of individual geometric features is often difficult to resolve. In the northern California Coast Ranges, The Geysers geothermal field is commonly interpreted to occur within a releasing stepover, although no single, clearly defined stepover is identified in published studies. To investigate the structural controls on The Geysers and the broader Clear Lake region, a two-dimensional elastic boundary element model is developed to evaluate spatial patterns of dilational strain associated with progressively more complete fault geometries. Model results show that dilation in the region is not controlled by a single structure but instead reflects the combined effects of multiple interacting fault elements. Three primary controls are identified: (1) opposing bends in the regional strike-slip fault system, including a releasing bend along the Maacama fault; (2) the southern fault tip of the Collayomi fault, which generates a prominent dilational lobe beneath the southern Geysers; and (3) a releasing stepover between the Collayomi fault and the Geyser Peak–Mercuryville–Big Sulphur Creek fault system, inferred to collectively behave as a right-lateral shear zone bounding the western margin of The Geysers. Predicted dilational strain magnitudes are sufficient to localize permeability between faults. These results highlight that incorporating complete fault networks and bedrock geological mapping can enhance geothermal assessments and provide a transferable framework for evaluating structurally controlled permeability in tectonically active regions.

California↗
USGS2025

Preliminary depth to basement modeling at Salton Sea, California

The San Andreas Fault – Imperial Fault (SAF-IF) transtensional step-over zone along the southern margin of the Salton Sea hosts substantial geothermal production and lithium brine resources. Recent volcanism at the Salton Buttes and active seismicity along the SAFIF fault system highlight active tectonic and magmatic processes that pose natural hazards and may impact energy and mineral production. Characterizing the subsurface architecture and extent of concealed alteration associated with this tectono-magmatic system enhances understanding of these active processes, associated hazards, and resources. We have compiled a gravity database, consisting of new and re-processed existing data, from which we have constructed a new isostatic residual gravity anomaly map of the Salton trough. We have used this new gravity dataset together with a compilation of publicly available borehole data to develop new depth to basement inversion models for the region. These depth to basement models help to constrain basin geometries, inform alteration mapping, and reveal variations in basement rocks. Due to the concealed nature of the complex tectonic framework at the Salton trough, it is necessary to utilize geophysical methods for subsurface characterization. These new depth to basement models are a first step toward constructing 2D and 3D geophysical and geologic models of the Imperial Valley and Salton Sea geothermal area. This analysis complements other geophysical initiatives, including magnetotelluric (MT) modeling (Tokmakoff et al., 2024), magnetic mapping (Glen and Earney, 2023, 2024) and potential field modeling, and seismic studies focused on hazard and resource investigations in the Imperial Valley.

California↗
USGS2024

Seismic response comparison of a historical masonry church subject to real and simulated ground motions

In recent years, advanced numerical models and high-performance computing have facilitated the utilization of ground motion time series in the assessment of the non-linear dynamic behavior of historic masonry structures. Since recorded accelerograms can be sparse for specific analysis conditions, stochastic ground motion simulations have become a viable alternative to overcome this limitation. This study simulates the recorded acceleration time series of the Central Italy 2016 earthquake event at the closest station to the town of Macerata using a site-based stochastic approach. The simulated motions are seismologically evaluated using a goodness-of-fit method in terms of various intensity measures. The simulated records, in conjunction with real records, are used to study the non-linear dynamic behavior of San Filippo Neri church located in Macerata. The church of San Filippo represents an important example of Baroque religious architecture in central Italy, which was damaged and closed off to the public after the 2016 earthquake events. The construction was investigated with a vast diagnostic campaign which included on-site testing and dynamic identification tests. The collected data is used to calibrate the dynamic response of a three-dimensional finite element model of the church. The model is finally used to compare the non-linear seismic responses under real and simulated ground motions with the site recorded damage. The results of structural responses demonstrate a strong agreement between the real and simulated records, providing evidence to support the validation of the site-based stochastic simulation.

Macerata↗