Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Water Process Engineering”

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,711 records · Page 92Linked to original sources
USGS2025

Fingerprinting magmatic REE deposit sources with zircon petrochronology

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

Conference Paper↗
USGS2024

Uncertainty and spatial correlation in station measurements for mb magnitude estimation

The body‐wave magnitude (⁠⁠) is a long‐standing network‐averaged, amplitude‐based magnitude used to estimate the magnitude of seismic sources from teleseismic observations. The U.S. Geological Survey National Earthquake Information Center (NEIC) relies on in its global real‐time earthquake monitoring mission. Although waveform modeling‐based moment magnitudes are the modern standard to characterize earthquake size, is important because (1) in many cases, waveform modeling is not possible (e.g., low signal‐to‐noise events), (2) is applicable over a broad range of magnitudes, ∼M 4–7, and (3) there is a many decades‐long history of estimating magnitudes. We use the NEIC Preliminary Determination of Epicenters earthquake catalog to investigate the uncertainty in NEIC station measurements. We show that measurements are spatially correlated, which can bias event ⁠, and we describe an empirical relation between this spatial correlation and station‐to‐station distance. We further describe an approach to mitigate bias from the spatial correlation. Accounting for the spatial covariance of observations can change the event from −0.15 to 0.07 units (10th to 90th percentile) for smaller events (⁠⁠). These smaller events have the largest standard deviations ranging from 0.05 to 0.15 units (10th to 90th percentile).

The Seismic Record↗
USGS2026

SlideDetect: Spatio-temporal landslide detection using a three-dimensional convolutional neural network

Landslides pose a serious and ongoing threat to both human lives and infrastructure worldwide; therefore, it is of interest to predict where and when landslides are likely to occur. Advances in machine learning techniques have spurred numerous studies aimed at estimating relative landslide propensity, but are limited to spatial (as opposed to temporal) prediction due to the sparsity of landslide timing data. We address this data gap by training SlideDetect, a 3-dimensional convolutional neural network (3D CNN), to identify landslides based on their spatial and temporal occurrence within multitemporal image stacks. We use an inventory of landsides triggered by the 2018 Hokkaido earthquake and two years of monthly composite optical imagery spanning this event. The model can identify not only landslide location but also landslide date with an area under the precision-recall curve (PR-AUC) of 0.84. We further present a new standard for presenting PR curve results that explicitly compares model performance at different confidence thresholds, allowing for clearer model evaluation and comparison. Our new approach to constraining landslide timing paired with this more consistent and objective method for evaluating model performance shows considerable promise, and with further application and testing, SlideDetect could enhance the data availability and tools needed to advance landslide hazard and risk assessments.

JGR Machine Learning and Computation↗
USGS2026

Integrating detrital magnetite geochemistry and (U-Th)/He chronometry as a sediment provenance tool in geologic and metallogenic terranes

Magnetite is ubiquitous in porphyry Cu systems and in sediment sourcing both barren and mineralized regions, with potential as an indicator mineral in concealed and coarsely-mapped terranes. We develop and test a workflow for integrated geochemistry and (U–Th)/He (He) dating for inferring detrital magnetite (DMt) provenance in these settings. The ca. 70 Ma Taurus porphyry Cu–Mo(–Au) district in eastern interior Alaska serves as a test case. DMt from streams draining porphyry-related mineralization was characterized by geochemistry and mineral inclusion and microstructure observations, complemented by similar data for potential porphyry and host rock sources. Principal component analysis and clustering of DMt geochemical data resolve multiple populations in our samples geochemically and texturally compatible with derivation from metamorphic, porphyry-related hydrothermal, and igneous sources. Hydrothermal magnetite comprises ∼16–50% of DMt nearest porphyry mineralization but diminishes to ∼4% ∼15 km downstream. Subsampled grains within populations yield ∼160–110 Ma, ∼70 Ma, ∼55 Ma, and ∼20 Ma magnetite He date modes. Combined with provenance, He dates capture Early Cretaceous regional exhumation of metamorphic host rock and Late Cretaceous porphyry Cu mineralization. DMt grains showing partial hematite replacement yield ca. 55–20 Ma dates regardless of source, overlapping regional warm/wet climatic intervals. We interpret Cenozoic dates to reflect exhumation to near-surface oxidizing conditions and(or) supergene weathering. Magnetite is thus a promising target phase for (1) tracking the spatiotemporal distribution of porphyry systems, and (2) linking the formation and exhumation of these systems to a regional geologic history, both in Alaska and globally.

Alaska, Yukon↗
USGS2025

Fossil diatoms in Arctic marine surface sediments

Diatoms are one of the main constituents of marine phytoplankton in the Arctic, and thanks to their siliceous skeletons, diatom fossils are relatively well preserved in sediments. Due to their species-specific sensitivity to different ocean conditions, their abundance and assemblages in sediments are routinely used by paleoceanographers to reconstruct the state of the surface ocean in the past (Koç 2007). By using statistical methods, we can obtain valuable knowledge about their ecological preferences (Oksman et al. 2019) and generate quantitative reconstructions of various parameters, such as sea-surface temperature and sea-ice concentration through time (Krawczyk et al. 2021; Sha et al. 2014). The Marine Arctic Diatoms (MARDI) working group (WG) ( pastglobalchanges.org/mardi ) aims to advance knowledge on marine-diatom ecology and diatom-based reconstructions by compiling and harmonizing data from surface-sediment samples across the Arctic.

Pages Magazine↗
USGS2026

Early Miocene volcanic rocks and associated tectonics, Lava Hills and southern Bristol Mountains, California

Volcanic rocks of latest Oligocene to early Miocene age form an east-west belt across part of the central eastern Mojave Desert from the Whipple Mountains on the east to the Rosamond Hills on the west. We term this the central belt because it is separated from northern and southern belts by swaths with no volcanic rocks. Limited geochronologic data indicate that much of the belt is latest Oligocene and early Miocene in age, about 24 to 19 Ma, a finding that is consistent with these rocks being overlain by the 18.8 Ma Peach Spring Tuff in many places. We describe Miocene geology in a central area of the belt, in the Lava Hills, southern Bristol Mountains, and southern Old Dad Mountains. Sedimentary basins formed coeval with early andesite to rhyolite volcanism, progressing from fluvial and lacustrine tuffaceous sandstone to volcanic lavas, tuffs, and breccias, indicating that early basins formed proximal to volcanic edifices. Higher strata are fluvial and lacustrine with lavas punctuating the sequence. Although basins may partly have been formed within topographic lows bounded by volcanic domes, plateaus, and stratovolcanoes, consistent stratigraphic sections over wide areas indicate that tectonic basin evolution affected broad areas. The volcanic section is capped by local basalt flows and the regional Peach Spring Tuff. Limited data on normal faults support interpretations of early extensional basin development caused by northeast-southwest oriented stretching. Later extension caused stratal rotations, tilting early deposits down to the southwest. This tilted and subsequently beveled basin architecture was overlain by the youngest volcanic deposits, primarily rhyolite and basalt. The Peach Spring Tuff, 18.8 Ma, lies within this upper unit. Similar stratigraphic and structural relations are exposed in the nearby Marble Mountains and Van Winkle Mountain sections, reinforcing that a broad area underwent similar volcanism and tectonism. In our study area the upper unit is only very gently tilted except near dextral strike-slip faults of the eastern California shear zone. These late Miocene to Recent faults are represented as four main fault zones spaced about 5 km apart, representing distributed shear north of the Bristol Lake basin.

California↗
USGS2022

South Carolina and Landsat

South Carolina, the eighth State admitted to the union, transcends its size with its deep, rich history; striking beauty; vast natural resources; and extensive cultural diversity. Home to part of the Blue Ridge Mountains of the Central Appalachians, the Upstate is graced with more than 100 waterfalls, while the Lowcountry borders the Atlantic Ocean with 187 miles of coastline and 35 barrier islands. Forests cover two-thirds of the State, and forestry and agriculture together, as agribusiness, make up South Carolina’s leading industry. Two historic crops—cotton and tobacco—still rank in the top 10 commodities, though corn and soybeans now rank higher. Poultry, cattle, peanuts, and flowers also make the list. South Carolina’s population totals more than five million. Other residents include a variety of wildlife, bird, reptile, and fish species, including Ursus americanus (black bears), Alligator mississippiensis (American alligators), and Tursiops truncatus (bottlenose dolphins). More than 100 tree species also reside in South Carolina, which pays homage to one with its “The Palmetto State” nickname. South Carolina’s subtropical climate, long coastline, and lower elevations make it highly susceptible to tornado and hurricane activity and coastal flooding. Projected sea-level rise is a growing concern. A view from space can help monitor and manage natural resources on the land and in rivers, marshes, and the coast. Landsat reveals not just what an area looks like now, but also insights from decades ago.

South Carolina↗
USGS2025

Drone-based radiometric surveys provide high-resolution mine waste characterization

Airborne radiometric surveys use passive geophysical techniques to characterize geochemical variations at or near earth’s surface. These methods have been used for a variety of mapping applications, including mineral resource evaluation. However, detailed characterization of smaller geologic targets, including mine waste features, requires flying at lower altitudes and with tighter line spacing than is feasible with traditional aircraft. Here, a small uncrewed aircraft system (sUAS) equipped with a radiometric sensor was used to acquire high-resolution gamma-spectrometry over small mine waste features and a low-grade stockpile in southwestern New Mexico. The sUAS radiometric system mapped local variability within each survey area and revealed ~2–10 m wide zones where radioelements K, Th, and U may be elevated 2–10× the surrounding material. Additionally, the sUAS radiometric data revealed radioelement variability across survey sites, which correlated reasonably well with variability seen in geochemical samples at each survey site, even though samples collected from individual sites showed high local variability. The sUAS data characterized local heterogeneity within mine waste and other small geologic targets at scales of a few meters to tens of meters, which is not possible with traditional crewed aircraft, and with continuity of coverage that is not possible with ground surveys, thus filling a key gap in geophysical survey spatial resolution.

New Mexico↗
USGS2025

Estimating earthquake source depth using teleseismic broadband waveform modeling at the USGS National Earthquake Information Center

The U.S. Geologic Survey National Earthquake Information Center (NEIC) monitors global seismicity, producing a catalog of earthquake source parameters in near-real-time to provide information that can help mitigate the societal impact of earthquakes. The NEIC commonly relies on teleseismic observations to constrain earthquake source parameters (e.g., location, depth, magnitude, and mechanism) due to a lack of local and regional observations. For these ‘teleseismic’ events, depth phase (i.e., pP , sP ) arrival time observations provide the best estimate on source depth. However, depth phases are often difficult to accurately identify and/or pick. Therefore, NEIC relies on waveform modeling, such as those determined from W-phase ( M ww ), body wave ( M wb ), and regional ( M wr ) moment tensor estimations, to provide constraints on source depth. While depth estimates from these approaches are informative, higher frequency observations provide more precise estimates because depth phases are more prominently observed at higher frequencies. Here, we present NEIC’s relatively high-frequency (~0.04 to 1 Hz) teleseismic waveform modeling approach, termed Synthetic Depth Phase Modeling (SynDepth), for determining source depth. SynDepth was developed to provide NEIC with a tool that enables rapid, accurate, and quantifiable estimates of earthquake source depth in cases where locator depths are not reliable. This relatively simple and fast procedure searches over 1 km-incremented source depths and an expanding triangular source-time function to find the best-fitting solution. We compare automatic SynDepth solutions for a dataset of 1,216 earthquakes (M5.5-M7.6) between 2017 and 2021 to NEIC-derived depth estimates from other methods. Our approach provides a robust depth estimate for earthquakes lacking local arrival time data, and it minimizes the need for analyst review of depth-phase picks ( pP , sP ) or using predefined ‘fixed’ depths.

Seismological Research Letters↗
USGS2024

The 2023 Alaska National Seismic Hazard Model

US Geological Survey (USGS) National Seismic Hazard Models (NSHMs) are used extensively for seismic design regulations in the United States and earthquake scenario development, as well as risk assessment and mitigation for both buildings and infrastructure. This 2023 update of the long-term, time-independent Alaska NSHM includes substantial changes to both the earthquake rupture forecast (ERF) and ground motion models (GMMs). The ERF includes numerous additions to the finite-fault model, considers two deformation models, and introduces updated declustering and smoothing algorithms in the gridded background seismicity model. For the Alaska–Aleutian subduction zone, megathrust earthquakes occur on an updated structural and segmentation model, and the moment magnitude (M) 8+ rupture and rate model include a logic tree branch that considers slip rates derived from geodetic models of interface coupling. The megathrust model considers multiple models of down-dip width, and magnitudes are computed using newly developed scaling relations. For subduction intraslab events and subduction interface events with M < 7, the 2023 update uses a smoothed seismicity model with rupture depths derived from Slab2. The 2023 model updates GMMs in all tectonic settings using the recently published Next Generation Attenuation Subduction (NGA-Sub) GMMs for subduction interface and intraslab events, and the NGA-West2 GMMs for active crustal settings. Collectively, additions and updates to the Alaska NSHM result in hazard increases across most of south-central Alaska relative to the previous model, published in 2007. These changes are primarily due to the adoption of updated rate models for the large-magnitude interface events and the NGA-Sub GMMs that have much higher aleatory variability (sigma), consistent with global observations, and that include models of epistemic uncertainty.

Alaska↗
USGS2025

Unveiling coseismic deformation from differenced legacy aerial photography and modern lidar topography: The 1983 M6.9 Borah Peak earthquake, Idaho, USA

The 1983 M6.9 Borah Peak, Idaho, earthquake is one of the largest historical normal fault earthquakes in the western United States. We quantified meter-scale vertical change along the 35 km-long rupture using topographic differencing of 1966 aerial imagery and 2019 lidar-derived data. The initial differencing results are largely obscured by horizontal and vertical georeferencing errors and flight-line stripes. Our error corrections are designed to be insensitive to the coseismic deformation and reduced error by 50%. We calculated vertical separation and resolved a maximum of 2.02 ± 0.46 m at Doublespring Pass. Our vertical separation measurements are generally consistent with those from prior studies using field data and post-earthquake topographic data. However, the differencing measurements are a few decimeters lower than these prior measurements, indicating that differencing can isolate historical from prehistoric earthquake deformation. Our study demonstrates that revisiting historical earthquakes can provide new insights into the magnitude and patterns of coseismic deformation.

Idaho↗
USGS2025

Mapping a Carrington storm

A map is presented of median 1-min-resolution peak geoelectric-field strength across the United States as would be induced by magnetic storms as intense as the 2 September 1859 Carrington storm. The map is constructed from two data sets: Magnetometer time series from 22 ground-based observatories recording 40 magnetic storms, and surface impedance tensors derived from magnetotelluric measurements acquired at 1616 survey sites across the contiguous United States. Carrington-class storm geoelectric fields are likely to be very strong in the United States East and Midwest; > 5.00 V/km at many places. In Virginia, strengths would likely range from 30.30 V/km, with a 68% confidence interval of [19.44,47.20] V/km, to as low as 0.05 [0.03,0.07] V/km. Comparison of model geopotentials with those measured on 30 long lines, indicates errors of about 18%. A Carrington-class storm would likely induce geoelectric fields with strengths 55% greater than for the 13–14 March 1989 storm.

contiguous United States↗
USGS2026

Relating glassy rind thicknesses to ambient air temperatures at the Lost Jim flow field in the Imuruk Lake volcanic field, Alaska

The Lost Jim flow field, in the Imuruk Lake volcanic field, Alaska, extends west ~34 km from a single vent, crossing subarctic tundra and currently touches several lakes and streams. The weighted mean of five 36Cl cosmogenic exposure ages from the Lost Jim pāhoehoe flow is 7.73 ± 0.37 ka, indicating this eruption occurred substantially after the eruption of the underlying Camille flow, which was emplaced at 39.7 ± 1.3 ka. Paleoclimate records indicate the period when the Lost Jim flow field was emplaced was after deglaciation, and the climate was similar to today. We propose that the emplacement of lava in these cold subarctic conditions can lead to faster cooling of the lava surface compared to lava emplaced in warmer locations such as mid- latitude cold deserts. Glass abundance in the outermost rinds at the Lost Jim flow field was on average 74 % with 6.4 mm thick rims, compared to 60 % with 2.9 mm rims for cold mid-latitude desert samples. We interpret increased glass content as a proxy for rapid cooling likely occurring partly during winter. Glassiness values varied less across vent, margin, and mid-flow locations when compared to the mid-latitude flows suggesting the Lost Jim flow field was broadly impacted by the subarctic climate as opposed to responding to local microclimates. Our results indicate that lava glassiness may be a useful environmental indicator of cooler (in this case subarctic) conditions.

Alaska↗
USGS2025

The National Map Corps—Federal Emergency Management Agency and Oak Ridge National Laboratory pilot project report

This report provides an overview of the U.S. Geological Survey National Map Corps —Federal Emergency Management Agency and Oak Ridge National Laboratory pilot project in St. James Parish, Louisiana, that began in February 2024 and ended at the end of March 2024. The project used the power of The National Map Corps’ volunteer community to improve building classifications in the original Federal Emergency Management Agency’s U.S.A. Structures dataset. The report highlights the project’s completion and details the work and results achieved through a collaborative effort to enhance geospatial data quality and utility.

Open-File Report↗
USGS2026

Geothermal district energy systems coupled with seasonal underground thermal energy storage: A U.S. techno-economic screening by climate and geology

In the United States, cooling-dominated commercial building loads can cause geothermal heat pump-based district energy systems to accumulate a long-term subsurface thermal imbalance, motivating the incorporation of seasonal underground thermal energy storage. We developed a transferable workflow to evaluate geothermal district systems that pair ground heat exchangers with seasonal underground thermal energy storage. Using standardized hourly loads for seven commercial buildings and a uniform cost framework, we simulated ten U.S. cities with a physics-based ground heat exchanger model, subsurface storage simulations, and economic assessment to isolate the roles of climate and hydrogeology. In cooling-dominated cities, underground thermal energy storage supplied the majority of annual cooling, cutting electricity use and summer peaks substantially while achieving levelized costs comparable to or below conventional chiller-boiler plants. In cooler climates, the storage share shrunk, required borefield size and costs rose, and levelized cost of energy increased nearly linearly with declining underground thermal energy storage fraction, indicating storage fraction as the primary economic lever. Sensitivity analysis showed capital risk dominated by borefield drilling and surface heating, ventilation, and air-conditioning and piping, with underground thermal energy storage costs secondary. This workflow provides a transparent foundation for site-specific design and screening of next-generation geothermal district energy systems.

Arizona, Illinois, Maryland, Michigan, Mississippi↗
USGS2026

PFAS mixture composition and internal exposure profiles shape biological responses under field-realistic exposure

Per- and polyfluoroalkyl substances (PFAS) occur as complex mixtures, yet mixture-dependent biological effects under environmentally realistic exposure conditions remain poorly understood. We conducted multiyear (2018, 2019, 2021) continuous-flow, field-based exposures of male fathead minnows ( Pimephales promelas ) using a low-PFAS reference well (REF; sum of measured PFAS (∑PFAS) 0.1–0.2 μg L –1 ) and PFAS-contaminated groundwater from a fire-training area (FTA) at Joint Base Cape Cod, Massachusetts. Dilution treatments enabled separation of concentration and mixture effects. Groundwater from well FTA1 was perfluorooctanesulfonate (PFOS) dominated (∑PFAS 10–31 μg L –1 ), whereas groundwater from well FTA2 had higher concentrations and was enriched in perfluorooctanoate (PFOA) and diverse precursors (∑PFAS ∼ 80 μg L –1 ). Despite comparable plasma ∑PFAS and PFOS in FTA1–100% and FTA2–50% on day-7, cumulative mortality reached approximately 35% in FTA2–50% and 17% in FTA2–100%, and secondary sex trait expression was reduced by approximately 65–80% relative to REF, whereas sperm motility effects were mixture- and time-dependent. Plasma PFAS profiles were dominated by perfluorohexanesulfonate (PFHxS), PFOS, and sulfonamide precursors. Liver transcriptomics from REF-acclimated fish revealed robust disruption of metabolic, mitochondrial, and endocrine pathways that link PFAS mixture chemistry and internal exposure profiles to organismal outcomes, and testis transcriptomics provided complementary insight into reproductive impairment. These results indicate that PFAS mixture composition and internal exposure profiles, including precursor-associated differences, are important determinants of ecotoxicological outcomes under field-realistic conditions.

Massachusetts↗
USGS2026

Season and antecedent conditions impact concentration-discharge relationships for dissolved organic carbon and alkalinity in southeast Alaskan watershed

Fluvial export of dissolved carbon plays an important role in watershed-scale biogeochemistry. Predicted changes in climate are expected to impact watershed hydrologic regimes, and in turn, the sources and export of dissolved carbon from watersheds. Here, we utilize high resolution measurements of discharge and dissolved carbon concentration to examine how concentration-discharge (CQ) relationships vary seasonally and during high flow events over the main runoff season (May–October) in a temperate forested watershed in Southeast Alaska. Concentration-discharge relationships for dissolved organic carbon (DOC) and alkalinity demonstrated strong seasonal patterns, with more linear relationships in May and June versus other months. Changing power law model slopes ( b values; the exponent in a power law regression between runoff and carbon yields) indicated potentially shifting watershed sources (biogenic vs. geologic) and contrasting dominant flowpaths (shallow vs. deeper groundwater) for DOC and alkalinity over the sampling period. During the largest storm event of the study, DOC and alkalinity b values shifted from an overall pattern of transport (mean b = 1.58 values >1.0 indicate transport limitation) and source limitation (mean b = 0.48, values <1.0 indicate source limitation) to chemostatic (DOC, b = 0.99; alkalinity, b = 1.019). In June through August, patterns in hysteresis index suggest that CQ relationships were altered when storms followed in close succession to each other. Together, these findings indicate that seasonal and antecedent flow conditions play a role in dissolved carbon export from forested watersheds. Understanding these dynamics, particularly during winter months, will become increasingly important as changes to hydroclimate impact riverine carbon export.

Alaska↗
USGS2026

pySATSI: A Python package for computing focal mechanism stress inversions

We introduce pySATSI, a Python package for computing earthquake focal mechanism stress inversions. This algorithm can handle a wide variety of types of stress inversion problems with a single script and can duplicate many capabilities of preceding methodologies. We also add new capabilities that include spatiotemporally variable inversion grids, damped stress estimates for clusters with few or no focal mechanisms, and variable fault‐plane ambiguities that the user can assign to individual events. In addition, we added the ability to use damped stress inversions with fault‐plane ambiguity probabilities that are weighted by fault instabilities. Our algorithm is computationally efficient with faster runtimes than previous algorithms, scales well for large datasets, and can be easily parallelized.

Seismological Research Letters↗