Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Bulletin of the Seismological Society of America”

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

Analyzing multi-year nitrate concentration evolution in Alabama aquatic systems using a machine learning model

Rising nitrate contamination in water systems poses significant risks to public health and ecosystem stability, necessitating advanced modeling to understand nitrate dynamics more accurately. This study applies the long short-term memory (LSTM) modeling to investigate the hydrologic and environmental factors influencing nitrate concentration dynamics in rivers and aquifers across the state of Alabama in the southeast of the United States. By integrating dynamic data such as streamflow and groundwater levels with static catchment attributes, the machine learning model identifies primary drivers of nitrate fluctuations, offering detailed insights into the complex interactions affecting multi-year nitrate concentrations in natural aquatic systems. In addition, a novel LSTM-based approach utilizes synthetic surface water nitrate data to predict groundwater nitrate levels, helping to address monitoring gaps in aquifers connected to these rivers. This method reveals potential correlations between surface water and groundwater nitrate dynamics, which is particularly meaningful given the lack of water quality observations in many aquifers. Field applications further show that, while the LSTM model effectively captures seasonal trends, limitations in representing extreme nitrate events suggest areas for further refinement. These findings contribute to data-driven water quality management, enhancing understanding of nitrate behavior in interconnected water systems.

Alabama↗
USGS2026

Predictions of anthropogenic background PFAS concentrations in soil and relation to bedrock lithology and groundwater quality

Detectable concentrations of per- and polyfluoroalkyl substances (PFAS) have been observed in soils in remote areas and presumably originate from atmospheric deposition. These anthropogenic background concentrations may enable some PFAS to leach to groundwater at levels that exceed regulatory criteria for drinking water. However, anthropogenic background soil concentrations and their connection to groundwater are not well characterized. We developed a boosted regression tree model to predict perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) concentrations in shallow soils across northern New England. Low soil pH was the most important predictor of elevated anthropogenic PFOS and PFOA concentrations in background soils, rather than potential PFAS sources, land use, or population density. Total organic carbon (TOC) was also an important predictor for PFOS soil concentrations. Model predictions indicate that 73% of the shallow soils within Maine, Vermont, and New Hampshire exceed New Hampshire’s Soil Remediation Standard for PFOS (0.5 ng/g) and 41% exceed the PFOA standard (0.4 ng/g). Analysis of soil model results and groundwater data suggests that areas with high soil pH are associated with higher groundwater detection frequencies, illustrating how areas with less retention in soil are, conversely, also areas with potentially greater groundwater vulnerability. Further analysis indicates that groundwater may be more vulnerable in calcareous lithologies.

Maine, New Hampshire, Vermont↗
USGS2025

Flow heterogeneity controls dissolution dynamics in topologically complex rocks

Rock dissolution is a common subsurface geochemical reaction affecting pore space properties, crucial for reservoir stimulation, carbon storage, and geothermal energy. Predictive models for dissolution remain limited due to incomplete understanding of the mechanisms involved. We examine the influence of flow, transport, and reaction regimes on mineral dissolution using 29 time-resolved data from 3D rocks. We find that initial pore structure significantly influences the dissolution pattern, with reaction rates up to two orders of magnitude lower than batch conditions, given solute and fluid-solid boundary constraints. Flow unevenness determines the location and rate of dissolution. We propose two models describing expected dissolution patterns and effective reaction rates based on dimensionless metrics for flow, transport, and reaction. Finally, we analyze feedback between evolving flow and pore structure to understand conditions that regulate/reinforce dissolution hotspots. Our findings underscore the major impact of flow arrangement on reaction-front propagation and provide a foundation for controlling dissolution hotspots.

Geophysical Research Letters↗
USGS2025

Metal-rich lacustrine sediments from legacy mining perpetuate copper exposure to aquatic-riparian food webs

Historic copper mining left a legacy of metal-rich tailings resulting in ecological impacts along and within Torch Lake, an area of concern in the Keweenaw Peninsula, Michigan, USA. Given the toxicity of copper to invertebrates, this study assessed the influence of this legacy on present day nearshore aquatic and terrestrial ecosystems. We measured the metal (Co, Cu, Ni, Zn, Cd) and metalloid (As) concentrations in sediment, pore water, surface water, larval and adult insects, and two riparian spider taxa collected from Torch Lake and a nearby reference lake. Overall, elevated metal and metalloid concentrations, particularly Cu, were measured in all sediment samples and some surface and pore water samples collected from Torch Lake. For instance, Cu concentrations in the Torch Lake sediment were ∼200% higher than the reference lake and all measured concentrations exceeded predicted effects concentrations by at least ninefold. Within larval insect tissues, we observed 160% higher Cu concentrations than measured in the reference lake, and Cu was the only measured element above predicted effects concentrations in Torch Lake. Adult insects collected at both lakes had similar metal concentrations irrespective of exposure levels. Yet we found 100% higher copper concentrations in Torch Lake riparian spiders, demonstrating elevated exposure risk to insectivores across the aquatic-terrestrial boundary. Our results highlight that other metals in the mixture may not be as concerning to adjacent riparian ecosystems, but copper remains a contaminant of concern in Torch Lake 60 years after mining ceased.

Michigan↗
USGS2025

Methods for estimating selected low-flow statistics at gaged and ungaged stream sites in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Conservation and Recreation, Office of Water Resources, computed selected at-site streamflow statistics at U.S. Geological Survey streamgages in and near Massachusetts and developed regional regression equations for estimating selected streamflows at ungaged stream sites in Massachusetts. Two sets of regional regression equations were developed: (1) the “mainland” equations, for mainland Massachusetts excluding the area covered by the second set, and (2) the “southeastern” equations, for the Plymouth-Carver-Kingston-Duxbury aquifer area in southeastern Massachusetts and for Cape Cod. The regression equations and at-site statistics may be used by Federal, State, and local water managers in addressing water-resources issues relevant in Massachusetts. Regional regression analyses for the mainland equations were developed to estimate the following 27 streamflow statistics: 99-, 98-, 95-, 90-, 85-, 80-, 75-, 70-, 60-, and 50-percent flow durations; monthly June, July, August, and September 90- and 50-percent flow durations; February, June, and August median of the monthly means; harmonic mean; and medians of the following annual low-flow frequency statistics: 7-day; 7-day, 2-year; 7-day, 10-year; 30-day, 2-year; and 30-day, 10-year. The analyses used 81 streamgages with minimal to no regulations in and near Massachusetts. The regression analyses determined that four basin characteristics—drainage area, combined hydrologic soils A and B, streamflow variability index, and annual mean temperature—were the only significant explanatory variables for the different mainland equations. Regional regression equations were developed for the Plymouth-Carver-Kingston-Duxbury aquifer area in southeastern Massachusetts and Cape Cod, because surface-water drainage areas and groundwater contributing areas do not always coincide in this area of the State. The regression analyses to estimate 10 flow durations from the 99th to 50th percentiles used 18 streamflow sites with some occasional minor regulations—because there are few unregulated streams in southeastern Massachusetts. The analyses determined that groundwater contributing area and storage (combined water bodies and wetlands) were the only significant explanatory variables in the southeastern equations.

Massachusetts↗
USGS2024

The EnMAP spaceborne imaging spectroscopy mission: Initial scientific results two years after launch

Imaging spectroscopy has been a recognized and established remote sensing technology since the 1980s, mainly using airborne and field-based platforms to identify and quantify key bio- and geo-chemical surface and atmospheric compounds, based on characteristic spectral reflectance features in the visible-near infrared (VNIR) and short-wave infrared (SWIR). Spaceborne missions, a leap in technology, were sparse, starting with the CHRIS/PROBA and EO1/Hyperion missions in the early 2000s, and providing spectroscopy data with limited spectral coverage and/or low data quality in the SWIR. Since 2019, several countries and agencies have successfully launched a number of spaceborne imaging spectroscopy systems into orbit or deployed them on the International Space Station (ISS) such as DESIS, PRISMA, HISUI, GF-5, EnMAP and EMIT. Among these recent missions, the German Environmental Mapping and Analysis Program (EnMAP) stands for its long-term development, sophisticated design with on-board calibration, high data quality requirements, and extensive accompanying science program. EnMAP was launched in April 2022 and, following a successful commissioning phase, started its operational activities in November 2022. The EnMAP mission encompasses global coverage from 80° N to 80° S through on-demand data acquisitions. Data are free and open access with 30 m spatial resolution, a high spectral resolution with a spectral sampling distance of 6.5 nm and 10 nm in the VNIR and SWIR regions respectively, and a high signal-to-noise ratio. In this paper, we aim to present the mission's current status, coverage, science capabilities and performance two years after launch. We show the potential of EnMAP for space-based imaging spectroscopy to operate in various environments, including high and low light levels, dense forests, Antarctic glaciers, and arid agricultural areas. EnMAP enables various applications in fields such as agriculture and forestry, soil compositional, raw materials, and methane mapping, as well as water quality assessment, and snow and ice properties. The results show that EnMAP's performance exceeds the mission requirements, and highlights the significant potential for contribution to scientific exploitation in various geo- and biochemical sciences. EnMAP is also expected to serve as a key tool for the development and testing of data processing algorithms for upcoming global operational missions.

Remote Sensing of Environment↗
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↗
USGS2025

Aftershocks in stress shadows are inconsistent with modeled static Coulomb stress changes

Aftershock triggering is commonly attributed to increases in static Coulomb stress. In some areas, termed "stress shadows", a decrease in Coulomb stress is predicted to suppress earthquake occurrence. However, aftershocks are often observed in the modeled stress shadows. We examine several hypotheses that attempt to reconcile these shadow aftershocks with the static Coulomb stress change model: (1) they appear to be in shadows because of inaccuracy in the stress change calculations, (2) they occur on faults of unusual orientation which actually experienced increased Coulomb stress, (3) they occur on faults with different frictional properties, not modeled well by Coulomb stress, and (4) they are secondary aftershocks triggered by prior aftershocks or afterslip. When tested on the 2016 Mw7.0 Kumamoto, Japan, and 2019 Mw7.1 Ridgecrest, California, aftershock sequences, none of these hypotheses can explain the majority of the shadow aftershocks, and taken together these hypotheses can explain only about half of these aftershocks. This implies that Coulomb stress modeling that lacks small-scale fault zone heterogeneity might be inadequate to fully capture the true static stress changes and/or that other physical triggering models are needed, for example transient processes such as delayed triggering by dynamic stress changes from the passing seismic waves.

California↗
USGS2025

ARCHI: A new R package for automated imputation of regionally correlated hydrologic records

Missing data in hydrological records can limit resource assessment, process understanding, and predictive modeling. Here, we present ARCHI (Automated Regional Correlation Analysis for Hydrologic Record Imputation), a new, open-source software package in R designed to aggregate, impute, cluster, and visualize regionally correlated hydrologic records. ARCHI imputes missing data in “target” records by linear regression using more complete “reference” records as predictors. Automated imputation is implemented using a novel, iterative algorithm that allows each site to be considered a target or reference for regression, growing the pool of complete references with each imputed record until viable gap-filling ceases. Users can limit artifacts from spurious correlations by specifying model-acceptance criteria and applying geospatial, correlation, and group-based filters to control reference selection. ARCHI provides additional functions for visualizing results, clustering records with similar correlation structures, evaluating holdout data, and interactive parameterization with an accessible and intuitive graphical user interface (GUI). This methods brief provides an overview of the ARCHI package, modeling guidelines, and benchmarking on two regional groundwater-level datasets from the Central Valley, CA and Long Island, NY. We evaluate ARCHI alongside widely used multivariate imputation software to highlight and contextualize its computational efficiency, imputation accuracy, and model transparency when applied to large, groundwater-level datasets.

California, New York↗
USGS2025

Linking stream-reach nitrogen loads and groundwater “reachsheds” to inform wastewater-nitrogen management actions, Cape Cod, Massachusetts

Study Region Cape Cod, Massachusetts, U.S.A. Study Focus Anthropogenic nitrogen (N) is a key factor in degrading groundwater and surface-water quality, particularly in coastal New England where onsite wastewater systems are prevalent. This study evaluated whether direct N-load measurements in streams on Cape Cod, Massachusetts, coupled with flow-path information from groundwater-flow models, can effectively identify potential land areas where nitrogen mitigation could substantially reduce loads to receiving waters. Nitrogen fluxes were measured along 63 stream reaches during winter and summer and paired with simulated groundwater recharge areas to identify and rank potential areas for reduction of nonpoint-source N inputs. New Hydrologic Insights for the Region Reach-scale nitrate-N loads ranged from −39.1–1182 kg-N/yr per 100 m of stream, indicating spatially variable groundwater inputs across seasons. “Reachsheds” — areas contributing groundwater recharge to specific stream reaches — were delineated using a regional groundwater-flow model. Strong correlations were found between observed N loads and land-use characteristics, especially the number of septic systems and total N inputs from the sum of considered sources. Observed N loads were moderately correlated with recharge area size and wastewater flow estimates. Correlating reach-specific groundwater N loads with land use and parcel-scale nitrogen-yield data identified reachsheds with the highest potential for N load reduction. This approach enables targeted implementation of restoration efforts to optimize nutrient management and support regional load reduction.

Massachusetts↗
USGS2024

New developments at the Center for Engineering Strong-Motion Data (CESMD)

The Center for Engineering Strong-Motion Data (CESMD), an internationally utilized joint center of the U.S. Geological Survey (USGS) and the California Geological Survey (CGS), provides a single access point for earthquake strong-motion records and station metadata from the CGS California Strong-Motion Instrumentation Program (CSMIP), the USGS National Strong-Motion Project (NSMP), the USGS Advanced National Seismic System, and other affiliates. The CESMD has been continuously improving its webtools to facilitate the access of strong-motion data and metadata for use in post-earthquake response and for scientific and engineering research applications. The Center provides raw and processed strong-motion data via the Engineering Data Center (EDC) and the Virtual Data Center (VDC) web portals. This paper focuses on the strong-motion products provided by the EDC where more than 48,000 records with peak ground accelerations greater than 0.1% g from over 2400 earthquakes are currently hosted. and on the ongoing efforts to develop data access tools and applications. The new developments and ongoing efforts in the EDC include: 1) enhancements to the CESMD webservices to facilitate access to station metadata, earthquake information, and strong motion records 2) new features to the interactive map interface, improving the visualization and access to earthquake, station, and record information, 3) efforts to develop a new web application tool for data format conversion from a number of data formats, 4) efforts to unify varying waveform data formats into a consistent format, 5) ongoing efforts to compile seismic station site geology, measured or inferred Vs30 values, shear-wave profiles, NEHRP site class, and available structural instrument deployment schematics, and 6) a special studies pages for research topic-specific ground motion datasets that offer uniform processing of records from a variety of sources.

Conference Paper↗
USGS2024

Separating signals in elevation data improves supervised machine learning predictions for hydrothermal favorability

A recent study identified topography (land surface elevation above sea level) as an important input dataset (feature) for predicting the location of hydrothermal systems in the Great Basin in Nevada. Yet, topography is generally a result of more than one geological process and may consequently contain multiple distinct signals. For example, the geologic evolution of the Great Basin has produced both crustal thickening (i.e., regional-scale trends in elevation) and thinning via Basin and Range extensional faulting (i.e., valley-scale topographic relief). We postulate that these geologic processes may affect the occurrence of hydrothermal systems differently. Therefore, we separate the regional trend from the valley-scale signal in the Great Basin, and then use them separately to evaluate the importance of each as predictors for hydrothermal favorability. Our prior work applying supervised machine learning (ML) using the data from the Nevada Machine Learning Project demonstrated that employing a training strategy that randomly selects negative training sites produces better performing models for predicting hydrothermal favorability than a training strategy that uses expert-selected negatives. The models created using both training strategies exhibited a west-east geographic trend in the predictions for the favorability of hydrothermal resources. These models generally predicted higher favorability in western Nevada and lower favorability in eastern Nevada. This west-east trend in predicted favorability correlates with elevation across the Great Basin, which trends higher from west to east. By separating the original elevation feature into distinct features for elevation trend (i.e., regional-scale topography) and detrended elevation (i.e., valley-scale or local relative topography), we find that models using the separated topographic signals consistently outperform competing models that use the original elevation feature. Although western Nevada still exhibits higher favorability than eastern Nevada, using separated signals for regional elevation and local structure reduces the west-east prediction trend in the region and emphasizes structures associated with hydrothermal upflow. This work emphasizes how carefully engineering features to represent geological conditions relevant to hydrothermal systems allows ML algorithms to detect important patterns for predicting hydrothermal resource favorability and leads to better model performance.

Conference Paper↗
USGS2025

Evidence of mineral alteration in a salt marsh subterranean estuary: Implications for carbon and trace element cycling

Subterranean estuaries (STE) in salt marshes are biogeochemically active zones where interactions between terrestrial groundwater and seawater drive complex cycling of carbon and trace elements, influenced by mineral dissolution. These systems, characterized by fine-grained organic-rich peat overlying permeable coastal aquifers, play a crucial role as a blue carbon sink, yet their geochemical dynamics remain poorly understood. We investigated dissolved trace elements, carbon, silica, and radium isotopes in a salt marsh STE (Sage Lot Pond, Waquoit Bay, MA) over seasonal and annual cycles. Our results reveal that groundwater and estuarine water circulation through marsh peat and aquifer sediments leads to enrichments of dissolved organic and inorganic carbon (DOC and DIC), Si, Ba, and Mn, with variable source/sink behavior of Fe and net removal of U. Submarine groundwater discharge dominated Ba fluxes, whereas pore water drainage from marsh peat acted as the main sink for U and source of Si. Fe cycling was variable, with terrestrial Fe largely removed as groundwater passed through the STE, consistent with Fe-sulfide and amorphous phase formation. Radium isotope ratios identified two distinct subsurface flow pathways, influenced by metal-oxide cycling and organic matter breakdown. Si production was decoupled from DIC, suggesting Si originates from mineral alteration, whereas DIC results from both mineral weathering and microbial respiration. Silicate mineral alteration, coupled with marsh pore water drainage, accounts for up to 16% of annual DIC exports (66 g C m −2 y −1 ), highlighting the importance of STEs in coastal carbon and trace element cycling, especially as marshes face environmental change.

Massachusetts↗
USGS2025

Towards mobile wind measurements using joust configured ultrasonic anemometer for applications in gas flux quantification

Small uncrewed aerial systems (sUASs) can be used to quantify emissions of greenhouse and other gases, providing flexibility in quantifying these emissions from a multitude of sources, including oil and gas infrastructure, volcano plumes, wildfire emissions, and natural sources. However, sUAS-based emission estimates are sensitive to the accuracy of wind speed and direction measurements. In this study, we examined how filtering and correcting sUAS-based wind measurements affects data accuracy by comparing data from a miniature ultrasonic anemometer mounted on a sUAS in a joust configuration to highly accurate wind data taken from a nearby eddy covariance flux tower (aka the Tower). These corrections had a small effect on wind speed error, but reduced wind direction errors from 50° to >120° to 20–30°. A concurrent experiment examining the amount of error due to the sUAS and the Tower not being co-located showed that the impact of this separation was 0.16–0.21 ms − 1 "> ms − 1 , a small influence on wind speed errors. Lower wind speed errors were correlated with lower turbulence intensity and higher relative wind speeds. There were also some loose trends in diminished wind direction errors at higher relative wind speeds. Therefore, to improve the quality of sUAS-based wind measurements, our study suggested that flight planning consider optimizing conditions that can lower turbulence intensity and maximize relative wind speeds as well as include post-flight corrections.

Alaska↗
USGS2025

Critical mineral inventory of select IOA-IOCG deposits, southwestern USA

Critical minerals are necessary for modern technology and strategic purposes. Their increasing importance requires finding new and nontraditional resources. Samples of ore, altered, and unaltered host rock were collected from 26 iron mines and prospects in California, Nevada, and Utah to assess the potential of these deposits to host economic quantities of different critical minerals. Geochemical analyses were conducted by 61 element ICP-OES-MS sodium peroxide fusion and major elements determined by WDXRF. These deposits concentrated many critical minerals beyond what is found in average upper crustal abundances, such as Sb, As, Bi, Co, Ga, Mg, Mn, Ni, Nb, Pd, REE, Sc, Te, Sn, Ti, W, V, and Zn. However, most of these are not concentrated enough in the ore to be considered as economic resources. Those critical minerals that are enriched enough in some of these deposits to possibly be considered as by-product commodities are Ni, REE, V, and potentially Co and Ga. These enrichments were not uniform, with REE more likely to be enriched in IOA deposits, whereas Co, Ga, Ni, and V could be found enriched in either IOA or IOCG deposits.

California, Nevada, Utah↗
USGS2025

Potential for continental scientific drilling to inform fault mechanics and earthquake science

Our understanding of fault mechanics and earthquake processes remains limited, largely due to minimal direct observations near active faults at seismogenic depths. This lack of data restricts our ability to accurately assess and mitigate both natural and human-induced seismic hazards. However, recent advancements in drilling capabilities and downhole sensing technologies offer an opportunity: the ability to observe the physical conditions within a volume near active fault zones. In this contribution, we highlight how scientific drilling can provide access to the near-fault environment, enabling measurements of the stress, temperature, fluid pressure, and rock properties at depths where ruptures initiate, propagate, and arrest. These observations are essential to refine models of earthquake nucleation and dynamic rupture, bridging gaps between laboratory experiments, numerical simulations, and surface observations. These insights can advance fundamental understanding in earthquake science but also support the development of more effective seismic hazard assessments and risk mitigation strategies.

Seismica↗
USGS2026

Monitoring changes in Landsat thermal features in urban and non-urban interfaces from 1986 to 2023 in two international urban centers: Implications for climate and global issues

Rapid urbanization is reshaping thermal environments worldwide, with the strongest impacts occurring at the interface between urban and non-urban areas. Impervious surfaces, as key indicators of urban expansion, are critical for monitoring urban growth and assessing surface urban heat island (SUHI) effects. Land use and land cover change (LULCC) provides an essential link between urban dynamics and their environmental and societal consequences. Here, we integrated the U.S. Geological Survey (USGS) Climate Global Issues (CGI) Land Cover Product with Landsat thermal time-series to investigate SUHI evolution in two contrasting metropolitan regions: Wuhan, China, and Brasília, Brazil. Using data spanning 1986–2023, we analyzed the relationships between land cover, Landsat-based land surface temperature (LST), and SUHI intensity, and identified persistent thermal hotspots. Results demonstrate that the land cover data utilized increases the accuracy of impervious surface mapping along urban–rural gradients. Average SUHI intensities were 3.4 °C in Wuhan and 3.3 °C in Brasília, with statistically significant warming trends of 0.04 °C/year and 0.01 °C/year, respectively. Maximum temperature proved to be a robust indicator of SUHI intensification, capturing long-term upward trends. Our findings highlight the important role of urban land cover dynamics in shaping temporal SUHI variability and hotspot emergence. This prototype framework demonstrates the scientific and policy value of combining long-term land cover monitoring information with satellite thermal monitoring to quantify and track SUHI at city scale, supporting sustainable urban planning and climate adaptation strategies.

Remote Sensing↗