SEARCH · USGS Science

Results for “Environmental Science: Processes and Impacts”

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

SourceResearch connectionFrom yearThrough year

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

1,662 records · Page 38Linked to original sources

USGS2025

Sand provenance boundary in the Mu Us Sandy Land of northern China

Desert dunes are often assumed to have uniform mineral compositions due to extensive mixing during lateral transport, which complicates provenance studies. The Mu Us Sandy Land in north-central China, near the East Asian summer monsoon precipitation boundary, experiences a wetter climate than most deserts. Climate wetting as a result of a warming climate, and the ‘Sand Control Project’ implemented by the Chinese government to decrease the lateral movement of sand dunes in this area provide an opportunity to study surface processes of sand production and transport. Previous studies using zircon U Pb geochronology and heavy mineral composition indicate distinct sand sources for the Mu Us Sandy Land: local basement-derived middle Yellow River sediments and recycled dried-up lacustrine sediments for the eastern part, and northeastern Tibetan Plateau-derived upper Yellow River sediments for the western part. However, zircons and heavy minerals only represent trace amounts of the bulk mineralogy within a sand dune, so broader provenance analysis targeting common minerals is essential. We focus on the area near the proposed provenance difference boundary between the western and eastern Mu Us Sandy Land, using comprehensive sampling and provenance techniques to confirm distinct provenance characteristics and delineate the provenance boundary. Our findings reveal that sand from most of the Mu Us Sandy Land originated from the erosion of local basement by the middle Yellow River and recycled local dried-up lacustrine sediments, whereas the southwestern corner and the neighboring western-central Chinese Loess Plateau received sediments from the distal northeastern Tibetan Plateau.

Mu Us Sandy Land↗

USGS2026

Estimation of parental abundance using hierarchical Bayesian modeling with data augmentation

Pedigree-based estimation methods leverage the fact that each offspring in a cohort is genotypically “marked” by its parents and represent a recent and promising toolset for estimating population dynamics. This includes pedigree accumulation estimators that model the “accumulation” of inferred unique parents within a given cohort to estimate parental abundance. Unlike close-kin mark-recapture approaches, which rely on intercohort comparisons, pedigree accumulation modeling can be completed solely using intracohort samples. This is particularly advantageous for semelparous species, where intercohort pairs are impossible and adult life stages can be difficult to sample without affecting their likelihood of successfully reproducing. Previous work has evaluated a range of estimators for such datasets, concluding that the non-parametric Chao estimator provides the most accurate and precise estimates for feasible levels of sampling effort. We used simulated data to evaluate an alternative estimator based on hierarchical modeling and data augmentation in a Bayesian framework. Results indicate that estimates from the hierarchical Bayesian estimator had comparable accuracy and better precision than both the previously tested Chao1 estimator and the improved iChao formulation across a range of sample sizes and sex ratios. Furthermore, the Bayesian estimator was far more robust to simulated errors in pedigree reconstruction, especially the presence of false negatives. Hierarchical Bayesian pedigree accumulation models can also provide additional insight into underlying reproductive ecology through their use of an explicit observation process, allowing for the incorporation or estimation of species- and population-specific reproductive dynamics. More broadly, the parametric nature of these models offers opportunities to efficiently pool information among datasets as well as to propagate uncertainty within more complex models.

Ecology and Evolution↗

USGS1995

Evaluation of ground-water flow and hydrologic budget for Lake Five-O, a seepage lake in northwestern Florida

Temporal and spatial distributions of ground-water inflow to, and leakage from Lake Five-O, a softwater, seepage lake in northwestern Florida, were evaluated using hydrologic data and simulation models of the shallow ground-water system adjacent to the lake. The simulation models indicate that ground-water inflow to the lake and leakage from the lake to the ground-water system are the dominant components in the total inflow (precipitation plus ground-water inflow) and total outflow (evaporation plus leakage) budgets of Lake Five-O. Simlulated ground-water inflow and leakage were approximately 4 and 5 times larger than precipitation inputs and evaporative losses, respectively, during calendar years 1989-90. Exchanges of water between Lake Five-O and the ground-water system were consistently larger than atmospheric-lake exchanges. A consistent pattern of shallow ground-water inflow and deep leakage was also evident throughout the study period. The mean time of travel from ground-water that discharges at Lake Five-O (time from recharge at the water table to discharge at the lake) was estimated to be within a range of 3 to 6 years. Flow-path evaluations indicated that the intermediate confining unit probably has a negligible influence on the geochemistry of ground-water inflow to Lake Five-O. The hydrologic budgets and flow-path evaluations provide critical information for developing geochemical budgets for Lake Five-O and for improving the understanding of the relative importance of various processes that regulate the acid-neutralizing capacity of softwater seepage lakes in Florida.

Florida↗

USGS2025

Hydrologic response of groundwater and streamflow to natural and anthropogenic drivers of change in headwaters of the upper Colorado River basin during recent wet (1982–1999) and drought (2000–2022) conditions

Study region: Headwaters of the upper Colorado River basin (UCOL), USA Study focus: Surface-water and groundwater numerical models incorporating water-use information were used to investigate changes in climate, water use, and simulated hydrologic responses of snow processes, evapotranspiration, groundwater, and streamflow during recent wet (1982–1999) and drought (2000–2022) periods in the headwater subregions of the upper Colorado River basin. New hydrologic insights for the region: Decreases in average streamflow between wet and drought periods ranged from 20 % in the Colorado River headwaters subregion to 23 % in the Gunnison River headwaters subregion. Like streamflow, average surface runoff was statistically less during the drought than the wet period, with decreases from 24–31 % in the headwaters. On a volume basis, runoff decreases were greater than streamflow decreases in both the Colorado River and Gunnison River headwaters. Although the amount of water-year groundwater discharge to streams remained nearly the same between the wet and drought periods, groundwater as a percentage of streamflow increased between the wet and drought periods, highlighting the importance of groundwater in sustaining streamflow during drought conditions. Multiple linear regression analyses revealed that snowmelt-only models were better than the best precipitation and temperature models at explaining streamflow variability from all headwater subregions for both the wet and drought periods.

Arizona, Colorado, New Mexico, Utah, Wyoming↗

USGS2025

Laboratory assessment for recovery of porcine circovirus 2 and porcine reproductive and respiratory syndrome virus using two types of commercially available hollow-fiber ultrafilters

Groundwater near swine farms is an uninvestigated reservoir for porcine reproductive and respiratory syndrome virus (PRRSV) and porcine circoviruses (PCVs). Enteric microorganisms are often collected from groundwater via dead-end ultrafiltration, but recovery of PRRSV and PCV with this method has not been assessed. We recovered PRRSV2 and PCV2 by dead-end ultrafiltration followed by polyethylene glycol (PEG) precipitation, nucleic acid extraction, and reverse-transcription quantitative real-time PCR. We also compared 2 commercial hemodialysis ultrafilters (Asahi Kasei Rexeed-25A, Nipro Elisio-25H) and compared PRRSV2 recovery in these filters to other waterborne microorganisms. On average, 8 ± 1% of PRRSV2 was recovered by dead-end ultrafiltration and PEG precipitation, compared to 25 ± 6% for adenovirus 41. Full-process recovery of bacteria in the same filters was 5–15%; Cryptosporidium parvum recovery was 42 ± 12%. PCV2 was detected in 4 of 12 replicate filters, but low stock concentrations precluded quantitative recovery estimates. Elisio-25H ultrafilters performed similarly to Rexeed-25A filters for all organisms tested and is an effective replacement for the Rexeed-25A, which is no longer available in the United States. Our recovery of PRRSV2 and PCV2 by dead-end ultrafiltration in the laboratory suggests that PRRSV2 detection limits are as low as 3–50 genomic copies/L in sample volumes of 100–1,500 L. Based on quantitative microbial risk assessment, these concentrations are relevant to PRRSV2 infection rates in the U.S. swine herd.

Journal of Veterinary Diagnostic Investigation↗

USGS2026

Transcriptional changes in wild Yukon River Chinook Salmon associated with Ichthyophonus infections

Objective We compared differentially expressed genes in Chinook Salmon Oncorhynchus tshawytscha with three divergent Ichthyophonus statuses (undetected, subclinical infections, or clinical disease; n = 100) to investigate associated transcriptomic responses. Disease associated with the fish parasite Ichthyophonus sp. was first diagnosed in adult Chinook Salmon from the Yukon River in the late 1980s and has subsequently been implicated in premature host mortality. Methods Chinook Salmon tissue sample collections and Ichthyophonus infection data were leveraged from a multi-agency collaboration during summer 2022 at three locations along the main-stem Yukon River that spanned nearly 2,000 km of freshwater migration (lower, middle, and upper river). We sequenced the transcriptome and compared this to infection status based on routine diagnostic procedures. Results Among the 17,569 genes for which messenger RNA was detected, we identified a transcription signature in the skeletal muscle that was associated with Ichthyophonus infections and included 53 differentially expressed genes. The differentially expressed genes and their pathways included those known for involvement in immune functions, energy synthesis, cellular breakdown, and reproduction—all processes that are known to be influenced by senescence during spawning migrations. Conclusions Results demonstrate a clear transcriptional difference between diseased fish (clinical disease group) and those in which Ichthyophonus was undetected, including identifying candidate markers for infection in this population. These results provide a foundation for development of nonlethal biomarkers to evaluate potential Ichthyophonus infections in Chinook Salmon based on gene transcription, protein products, or gene variants (e.g., polymorphisms).

Journal of Aquatic Animal Health↗

USGS2025

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↗

USGS2024

Are field observations of surface rupture useful? An example from the 2023 Mw 7.8 Pazarcık, Turkey (Türkiye), earthquake

Field investigations have long been an important component of the scientific response to surface‐faulting earthquakes. However, in light of advances in remote data and models, the question arises whether field‐based observations of surface rupture remain useful for understanding rupture processes and seismic hazards. We approach this question using a field‐based study of the central 2023 M w 7.8 Pazarcık, Turkey (Türkiye), earthquake rupture, at the intersection of the east Anatolian fault (EAF) and Narlı fault. Our field observations include the surface rupture expression and extent of the central EAF and northernmost Narlı fault in generally forested and steep terrain and 68 measurements of left‐lateral surface displacement. These data improve our understanding of the Pazarcık rupture complexity, resolve the surface geometry of the Narlı fault–EAF connection, and exhibit a clear (>2 m) change in surface displacement across this intersection zone that confirms remote‐based coseismic slip models. Our study shows that focusing field efforts in areas of obscured or low‐resolution remote data can yield essential data for refining rupture extent, documenting perishable on‐fault displacement, and improving postearthquake situational awareness. A comparison of similarly large‐magnitude continental surface‐rupturing earthquakes indicates that displacement uncertainties relate to a complex set of factors, including measurement methods, rupture complexity, and displacement magnitude. Our study validates the need for postearthquake field observations, which, when driven by clear motivating questions and knowledge of methodological strengths and limitations, provide high‐resolution rupture data that complement remote‐based models.

Pazarcık↗

USGS2025

Using peak geometry and shifts in the x-ray spectrum of carbon from electron probe microanalysis to determine thermal maturity of organic matter

During the burial of mudstones, the associated organic matter undergoes gradual thermal maturation, a key process that can influence the reactivity of organic matter during catagenesis, the formation of hydrocarbon deposits and the chemical weathering of mudstones. Conventional methods for assessing the thermal maturity of organic matter often fail to reflect the geochemical heterogeneity between individual organic phases in mudstone samples. Here, we report an alternative, non-destructive, surficial and micro-scale (analytical spot size of ~ 300 nm with about 4 μm diffusion depth for micrometre-size organic grains) method to evaluate the thermal maturity of organic matter in mudstones using the carbon K α X-ray spectrum measured by field emission-electron probe microanalyser (FE-EPMA). Using this method, we observed correlations between parameter values derived from FE-EPMA spectra, including the peak position, the peak area and the intra-sample heterogeneity of these measurements, and independently measured vitrinite/solid bitumen reflectance for a suite of mudstones, representing different age, geological context and burial depth. With the increased values in peak area and position, we identified an increase in the carbon mass fraction of organic matter and the mean nominal oxidation state of carbon approaching zero. These trends, which are consistent with aromatisation and graphitisation, provide the rationale for using FE-EPMA to estimate the thermal maturity of organic matter. To explore some of these trends in more detail, we employed time-of-flight secondary ionisation mass spectrometry, X-ray photoelectron spectroscopy and optical reflectance measurements on a subset of samples.

Geostandards and Geoanalytical Research↗

USGS2026

Earthquake-hazard exposure of residents with potential access and functional needs in the United States

Earthquake response plans and earthquake early warning (EEW) systems designed for general populations may not consider potential access and functional needs (AFN) of individuals with physical, sensory, cognitive, or social limitations. Previous efforts to map the distribution of these populations have focused on social-vulnerability indices that ignore or oversimply these limitations. The descriptive and exploratory analysis summarized in this United States (U.S.) case study addresses this gap by identifying and integrating spatially explicit data for AFN-related residential populations, earthquake hazards, and county and county equivalents for the conterminous U.S., Alaska, Puerto Rico, and Hawaii. We focus on 13 AFN-related attributes that relate to an individual's ability to access information contained in an EEW alert, to understand and process earthquake information or observed ground shaking, and to take self-protective actions based on this information and physical cues of an earthquake. Depending on the demographic attribute, there are millions to tens of millions of U.S. residents with AFN-related attributes in areas considered to have varying likelihoods (2%, 10%, and 50%) of exceedance of a damaging earthquake in the next 50 years. Although these amounts represent low percentages at the national level, the percentage of individuals with AFN-related attributes in many counties and county equivalents substantially exceeds national percentages. No one county, county equivalent, U.S. state, or U.S. territory has the highest percentage of individuals in all AFN-related attributes; therefore, future efforts to increase individual resilience to earthquakes may benefit from understanding the local context of individuals with potential access and functional needs.

International Journal of Disaster Risk Reduction↗

USGS2024

Volcanic gas monitoring

Introduction As magma rises through the crust, decreasing pressure conditions allow volatiles to exsolve from the magma. These volatiles then migrate upward through the crust, where they can be stored at shallower levels or escape to the atmosphere. Rising magma also heats rock masses beneath volcanic centers, causing water in shallow aquifers and hydrothermal systems to boil and release additional gases and steam (see chapter F , this volume; Ingebritsen and Hurwitz, 2024). The chemistry and quantity of gases that reach the surface during periods of quiescence or volcanic unrest can reveal that gas-rich magma is ascending, crystallizing, or alternatively stalling, with important implications for volcanic hazard (for example, Sutton and others, 1992; Aiuppa and others, 2007, 2021; Werner and others, 2009, 2011, 2012; Moretti and others, 2013; de Moor and others, 2016; Lewicki and others, 2019; Edmonds and others, 2022; Kern and others, 2022; Kunrat and others, 2022). Most volcanoes in Alaska and the western United States are characterized by weak degassing, with one or more low-temperature fumaroles (typically near the local boiling temperature of water) and connect to a deeper and sometimes extensive hydrothermal system (for example, McGee and others, 2001; Symonds and others, 2003a, b). Hydrothermal systems will affect the chemistry of rising gases exsolved from deeper magma (Symonds and others, 2001), including sulfur dioxide (SO 2 ), hydrogen chloride (HCl), and water vapor (for example, Doukas and Gerlach, 1995; Gerlach and others, 1998, 2008; Symonds and others, 2001; Werner and others, 2013). As an example, depending on factors such as temperature, pressure, and oxidation state, rising SO 2 will react with groundwater to form hydrogen sulfide (H 2 S) gas, dissolved sulfate (SO 4 2− ), or elemental sulfur (Christenson, 2000; Symonds and others, 2001; Werner and others, 2008). The reaction and dissolution of SO 2 into shallow groundwater is commonly referred to as scrubbing, and can reduce the likelihood that ascending, degassing magma can be detected. Carbon dioxide, however, in addition to exsolving from magma early in the ascent process, is not easily removed by hydrothermal fluids (Lowenstern, 2001). As scrubbing and other processes take place, the SO 2 /H 2 S, CO 2 /SO 2 , and CO 2 /H 2 S ratios may change. High rates of SO 2 emission indicate that magma has moved to relatively shallow levels in the volcano and that the system has heated up enough to establish dry pathways from depth to the surface. Monitoring multiple gas species and the total output of those species is thereby useful for volcano monitoring during both periods of quiescence, to establish background degassing conditions, and during unrest, when gas geochemistry and emission rates can provide information on changing conditions, such as magma ascent. To provide context for multidisciplinary volcano forecasts, we focus on the following two key required capabilities: (1) characterizing baseline geochemistry and gas discharge from volcanoes and volcanic regions and (2) monitoring changes in gas geochemistry and discharge to inform forecasts of volcanic eruptions and their effects. Sufficient baseline data must be collected to identify and interpret anomalous degassing associated with volcanic unrest (for example, Sorey and others, 1998; Rouwet and others, 2014). Differences in volcano type, baseline degassing rates, local hydrology, and geography (for example, high versus low latitude) will result in a different baseline for each volcano. Volcanoes of any threat level that exhibit one or more degassing phenomena would ideally be monitored by techniques needed to establish baseline degassing data, with the sampling frequency of baseline data dictated by the threat level (table E1). Additional monitoring techniques become necessary during periods of unrest. In general, three of the most important techniques for gas monitoring are (1) direct sampling of fumarole, spring, and soil gases for laboratory geochemical measurements, (2) measurements of the chemical composition of the volcanic plume and emission rates of major gas species (for example, H 2 O, CO 2 , SO 2 , and H 2 S) by satellite, airborne, or ground-based techniques, and (3) measurements of diffuse emissions of CO 2 and other gases through soils. Various methods and instruments may be useful both for baseline studies and during unrest.

Scientific Investigations Report↗

USGS2024

Peri-Gondwanan sediment in the Arkoma Basin derived from the north: The detrital zircon record of a uniquely concentrated non-Laurentian source signal in the late Paleozoic

During the assembly of Pangea, peri-Gondwanan terranes collided with the eastern and southern margins of Laurentia and brought with them unique detrital zircon U-Pb signatures. Discriminating between individual peri-Gondwanan terranes in the detrital record is difficult due to their similar geologic histories. However, characterization of this provenance is critical for understanding late Paleozoic sediment routing during development of Pangea. Along southeastern Laurentia, in the Arkoma Basin (present-day Arkansas and eastern Oklahoma, southeastern United States), we identified Middle Pennsylvanian (Desmoinesian) strata that exhibit a concentrated peri-Gondwanan detrital zircon signature (e.g., ca. 800–550 Ma). Although several southern peri-Gondwanan terranes (e.g., Maya, Suwannee) are closer to the Arkoma Basin, geologic data, such as predominantly north-to-south paleocurrents and proximal-to-distal facies relationships in these Desmoinesian strata, support a northern source (e.g., Ganderia, Avalonia, Meguma). Further evidence of a northern source comes from detrital zircon source mapping, which reveals the persistence of this peri-Gondwanan signal in depocenters to the north of the basin after the signal had diminished in the Arkoma Basin. To this end, bottom-up detrital zircon source modeling, source mapping, regional stratigraphy, paleocurrent data, and sandstone petrography allow us to reconstruct the evolution of this Middle Pennsylvanian (Desmoinesian) sediment pathway in the context of intraplate and plate-margin tectonic activity. This reconstruction documents processes affecting Earth’s surface (e.g., tectonics, climate) during the assembly of Pangea and describes in detail part of a dynamic continental-scale drainage system.

Geosphere↗

USGS2025

Land application of drill waste: A scope analysis

Drilling fluid waste land application, a process where drilling wastes are spread and tilled into the land surface, has become common in some petroleum-producing states, however, the potential benefits and risks of this practice are not well studied. Drilling fluids can be water- or oil-based and can have high concentrations of total soluble salts and total petroleum hydrocarbons. Comprehensive chemical characterization of these fluids is not well documented in the literature, and the extent of land application is largely unknown. We hypothesized that the land application of drill waste would fluctuate over time due to economic factors. To begin to understand the extent of historical and potential future land application, we analyzed data from over 5,800 drilling fluid land application permits collected by the Oklahoma Corporation Commission for years 2000, 2005, 2010, and 2015–2020. During the years studied, drilling fluid wastes were applied to more than 250,000 acres in Oklahoma, with over 54,000 thousand barrels (Mbbl) of liquids and nearly 21,000 Mbbl of solids applied. Land application is widespread (occurring in 59/77 counties), however recent drilling activity, land availability, and the economics of transportation have created conditions favorable for land application specifically in the Anadarko Basin. Land application can co-occur with sensitive areas, such as important groundwater and surface-water drinking sources and agricultural fields used for subsistence or feed crop production. Our approach for quantifying the extent of land application, along with further chemical characterization studies, can aid operators and land managers who are considering this practice in assessing the associated benefits and risks.

Oklahoma↗

USGS2025

Rainfall thresholds for postfire debris-flow initiation vary with short-duration rainfall climatology

The size, frequency, and geographic scope of severe wildfires are expanding across the globe, including in the Western United States. Recently burned steeplands have an increased likelihood of debris flows, which pose hazards to downstream communities. The conditions for postfire debris-flow initiation are commonly expressed as rainfall intensity-duration thresholds, which can be estimated given sufficient observational history. However, the spread of wildfire across diverse climates poses a challenge for accurate threshold prediction in areas with limited observations. Studies of mass-movement processes in unburned areas indicate that thresholds vary with local climate, such that higher rainfall rates are required for initiation in climates characterized by frequent intense rainfall. Here, we use three independent methods to test whether initiation of postfire runoff-generated debris flows across the Western United States varies similarly with climate. Through the compilation of observed thresholds at various fires, analysis of the spatial density of observed debris flows, and quantification of feature importance at different spatial scales, we show that postfire debris-flow initiation thresholds vary systematically with short-duration rainfall-intensity climatology. The predictive power of climatological data sets that are readily available before a fire occurs offers a much-needed tool for hazard management in regions that are facing increased wildfire activity, have sparse observational history, and/or have limited resources for field-based hazard assessment. Furthermore, if the observed variation in thresholds reflects long-term adjustment of the landscape to local climate, rapid shifts in rainfall intensity related to climate change will likely induce spatially variable shifts in postfire debris-flow likelihood.

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

USGS2026

Localization of spatiotemporally heterogeneous subsurface flows using autoencoder-based deep learning framework for time-lapse self-potential tomography

Self-potential (SP) monitoring has emerged as a valuable method for characterizing subsurface hydrogeological features and processes due to its sensitivity to fluid-induced electrokinetic effects. Despite advancements in SP inversion, challenges remain in imaging groundwater dynamics from SP activities due to complex hydrological settings and transient noise. In this study, a deep learning autoencoder (AE)-based framework is proposed for the spatiotemporal localization of subsurface fluid movement from time-lapse SP tomography. Temporal segments of time-lapse numerical inversions were first derived from long-term SP monitoring conducted from a floodplain site in Oak Ridge, Tennessee, known for active hyporheic exchange. Subsequently, AE models based on vision transformer (ViT), convolutional long short-term memory (ConvLSTM), convolutional neural network, and temporal convolutional network were individually trained and compared on the SP tomography segments for reconstruction performance. Finally, the reconstruction error over time serves as an anomaly score to identify moments of active SP variation, whereas spatial distributions of errors within these moments are analyzed to image and localize regions associated with anomalous subsurface fluid movement. The results demonstrate that ConvLSTM- and ViT-AE are most capable for the localization task with contrasting error distributions and consistent delineation of anomalies. Applying the method to both SP arrays parallel and perpendicular to the stream produced consistent anomaly zones near a fault or karst feature, validating the robustness and generalization of the approach. These results demonstrate the potential of the proposed framework as a scalable and interpretable tool for spatiotemporal analysis of subsurface flow dynamics in complex hydrogeological systems.

Tennessee↗

USGS2024

Earthquake-triggered ground-failure inventory associated with the M7.1 2018 Southcentral Alaska earthquake

The 30 November 2018, magnitude (Mw) 7.1 earthquake in Southcentral Alaska triggered substantial landslides, liquefaction, and ground cracking throughout the region, resulting in widespread geotechnical damage to buildings and infrastructure. Despite a challenging reconnaissance and remote-sensing environment, we constructed a detailed digital inventory of ground failure associated with the event from several sources. Sources included information derived from remotely sensed data, and data compiled from literature, social media postings, and earthquake damage information compiled by local, state, and federal agencies. Each instance of ground failure within the inventory contains information on the location and type of observed ground failure, and the methods and data used to document the occurrence. Where high-quality data, such as LIDAR or satellite imagery, were available and showed the ground-failure instance clearly, the extent is mapped as a polygon or polyline. All other locations are mapped as points. There are a total of 886 ground-failure instances documented within the inventory (400 landslides, 286 liquefaction features, and 200 features unattributed to specific processes). A semi-quantitative confidence scheme is used to describe mapping certainty associated with each ground-failure feature. This inventory represents a relatively moderate ground-failure-triggering event that occurred in a subarctic environment. This data paper describes the content within the inventory, the inventory data collection procedures, and limitations of the data. Events of this type are not often documented in detail; thus, adding the inventory data to the US Geological Survey Open Repository of Earthquake-Triggered Ground-Failure Inventories further diversifies the datasets available to the scientific community to be used to better understand and model earthquake-triggered ground failure.

Alaska↗

USGS2026

Syn-magmatic subsidence during the early stages of continental rifting in the Mesoproterozoic—A reanalysis of legacy data for the Midcontinent Rift, western Lake Superior

The Midcontinent Rift system (ca. 1.1 Ga) is a 2000-km-long series of elongated volcanic and sedimentary troughs and associated intrusive centers exposed chiefly in the Lake Superior region of North America. The rift system represents a long history of intense magmatism and subsequent sedimentation that was arrested by far-field tectonic events before sea-floor spreading was established. The premature cessation preserved a record of processes related to the beginning of continental rifting. The rift system under Lake Superior has been long studied using seismic-reflection data collected as part of the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE). We reexamine GLIMPCE Line C by developing a detailed velocity model for time to depth conversion constrained by other legacy data. We corroborate the model and develop a geologic interpretation using gravity and magnetic modeling and ties to geology mapped onshore. We recognize superposed subsiding sedimentary and volcanic basins for the southern half of the Line C depth section. This interpretation differs from previous paradigms that show major crustal faults that bound half-grabens or full grabens. We conclude that high-velocity (6.9 km/s) intrusive zones rather than major crustal faults border the sides of the basins. We speculate that the volcanic basin represents the initiation of seaward dipping reflectors. The syn-magmatic subsidence can be explained by dike injection and volcanic loading. Discrete lava basins throughout the region likely subsided at different times in a disorganized manner along the rift trend, raising questions about the long-term role of lithospheric thinning and melt generation.

Michigan, Minnesota, Ontario, Wisconsin↗

USGS2025

Anatectic origin of Mississippian spodumene-bearing pegmatites in western Maine during orogenic plateau collapse

Spodumene pegmatites are an important lithium source, but the processes and tectonic settings in which they form are poorly understood. The Rumford pegmatite district surrounding Plumbago Mountain, western Maine, is host to numerous spodumene pegmatites, including the Plumbago North pegmatite (a world-class spodumene resource). Competing petrogenetic models for these spodumene pegmatites include (1) highly fractionated melts of the Mooselookmeguntic igneous complex and (2) anatexis. We tested these hypotheses by constraining the geologic, magmatic, metamorphic, and tectonic history of the Plumbago Mountain area with detailed geologic mapping and U-(Th)-Pb geochronology. The Silurian Rangeley Formation records initial isoclinal folding prior to, and contact-related metamorphism synchronous with, the intrusion of the 417 ± 4 Ma Plumbago Mountain pluton. Peak amphibolite facies metamorphism and crustal melting occurred during the ca. 410 to 400 Ma Acadian orogeny. Pulsed emplacement of the Mooselookmeguntic igneous complex occurred between ca. 389 and 356 Ma. Cassiterite U-Pb dates of spodumene pegmatites (333–327 Ma) are ≥23 m.y. younger than nearby granitic plutons, strongly arguing against the fractional crystallization model. Metamorphic monazite and xenotime (346–328 Ma) and 330 to 308 Ma 40 Ar/ 39 Ar hornblende dates indicate metamorphism coeval with spodumene pegmatite emplacement, supporting anatectic models. Reheating, anatexis, and spodumene pegmatite emplacement occurred during collapse of the 380 to 330 Ma Acadian orogenic plateau. Lithium enrichment may be linked to one or more stages of partial melting of metasedimentary and plutonic rocks during the formation, tenure, and collapse of the Acadian altiplano and emphasizes the role of anatexis in producing spodumene pegmatites of economic significance.

Maine↗