SEARCH · USGS Science

Results for “Journal of Veterinary Diagnostic Investigation”

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 32Linked to original sources

USGS2024

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

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

Conference Paper↗

USGS2024

A great tsunami earthquake component of the 1957 Aleutian Islands earthquake

The great 1957 Aleutian Islands earthquake ruptured ∼1200 km of the plate boundary along the Aleutian subduction zone and produced a destructive tsunami across Hawaiʻi. Early seismic and tsunami analyses indicated that large megathrust fault slip was concentrated in the western Aleutian Islands, but tsunami waves generated by slip in the west cannot explain the large observed runup in Hawaiʻi far to the southeast. Recently mapped 1957 geologic deposits on eastern Aleutian Islands suggest occurrence of very large nearby slip. Jointly modeling tsunami runup along the eastern Aleutian and Hawaiian Islands together with tide gauge recordings across the Pacific resolves 12-26 m shallow slip along 600 km of the eastern Aleutian Islands in addition to modest, deeper western slip inferred from seismic records. The eastern near-trench slip results in an M W 8.3-8.6 tsunami earthquake component of the M W 8.6-8.8 rupture, comparable in size to the adjacent 1946 Aleutian tsunami earthquake to the east. The reexamination of the 1957 rupture confirms the tsunami hazards posed by the eastern Aleutian subduction zone to Hawaiʻi and lays the groundwork for investigation of large prehistoric earthquakes through modeling tsunami runup inferred from stratigraphic observations to constrain their rupture processes.

Alaska↗

USGS2026

Tectonic controls on volcanism and associated hydrothermal activity in a sediment-dominated mid-ocean ridge; Escanaba Trough

Mid-ocean ridges, the Earth's most extensive volcanic system, exhibit unique characteristics in sediment-dominant environments. Thick sediment cover insulates the crust and channels fluid along pathways that can lead to the formation of distinct crustal alteration patterns, exceptionally large mineral deposits, and specialized chemosynthetic ecosystems. This study presents an interdisciplinary investigation into the tectonics of the Escanaba Trough, a heavily sedimented axial valley at the southern Gorda Ridge in the Northeast Pacific Ocean. A primary challenge in such environments is overcoming the masking effect of thick sediments on basement structures that control magmatic and hydrothermal activity. We address this by employing three-dimensional (3D) magnetic modeling of high-resolution near-seafloor magnetic data collected by an autonomous underwater vehicle (AUV). The 2022 surveys with AUV Sentry provided data for 3D magnetic susceptibility models, refining our understanding of the geometry of sub-sediment laccoliths/saucer-shaped sills and hydrothermal alteration. In conjunction with a new 1:100,000 scale lithostratigraphic map, we outline the tectonic controls on the emplacement of Escanaba Trough's three main volcanic centers, characterize the geometry of its spreading segments, and provide volumetric data on the distribution of sub-sediment volcanism in the southern Gorda Ridge.

Escanaba Trough, Pacific Ocean↗

USGS2025

Geochemical processes related to mined, milled, or natural metal deposits in a rapidly changing global environment

The demand for metals and raw materials, such as nickel and copper, has been projected to expand in the coming decades, driven by the global energy transition, the need for green technologies, and expanding infrastructure. Consequently, the increasing extraction and production of mining waste can have adverse impacts on surrounding environments and human health. The aim of this thematic collection is to fill critical knowledge gaps in the present-day cycles of metal(loid)s from source to larger sinks, and the effect of environmental management, anthropogenic development, and climate change. Altogether, the studies have been conducted in different natural settings around the world and comprise investigations in laterites, a soil-medicinal plant system, watersheds, and banded iron formations, among others. The geochemical applications in tracing mineralization, its secondary products, and/or potential impact on the immediate environment are highly diverse with applied tools ranging from isotope tracers to major and trace element systematics. Particularly the use of rare earth elements, their patterns and anomalies are methods employed by several studies in this collection. We summarize the findings to offer a potential future direction for the use of geochemical tracing techniques in resource exploration in the context of climate change and environmental challenges.

Geochemistry: Exploration, Environment, Analysis↗

USGS2025

Earthquake stress-drop values delineate spatial variations in maximum shear stress in the Japanese forearc lithosphere

Earthquake stress drop (Δσ) may increase with depth and stress in the brittle lithosphere. However, the range of uncertainty in Δσ and the lack of constraints on absolute stress make it difficult to establish whether they are correlated. Here, we investigate Δσ dependence on depth and maximum shear stress ( τ max ) based on ~11 years of seismicity in the northeastern Japanese forearc following the 2011 Tohoku-Oki megathrust earthquake. We interpret Δσ estimates computed using both individual spectra and spectral-ratio methods and find that Δσ exhibits a clear depth dependence within the seismically active upper ~60 km of the forearc lithosphere ( ~ 0.8 MPa per 10 km). We further compare Δσ values with quantitative τ max estimates from finite-element models of force balance. We find that median Δσ values increase with τ max in the brittle forearc lithosphere and that earthquake stress release is proportional to τ max . The dependence of Δσ on τ max explains the apparent depth dependence of Δσ and suggests that average Δσ values provide a relative measure of the stress at failure. In the northeastern Japanese forearc, Δσ values remained roughly constant in the decade following the Tohoku-Oki earthquake, suggesting negligible changes in failure stress in the forearc since the mainshock.

Communications Earth and Environment↗

USGS2023

Reconnaissance basement geology and tectonics of North Zealandia

New rock dredge samples supply key information to establish the tectonic and geological framework of the northern two-thirds of the 95% submerged Zealandia continent. The R/V Investigator voyage IN2016T01 to the Fairway Ridge, Coral Sea, obtained poorly sorted poly-lithologic pebbly to cobbly sandstones, well sorted fine grained sandstones, mudstones, bioclastic limestones, and basaltic lavas. Post-cruise analytical work comprised petrography, whole rock geochemical and Sr and Nd isotopic analyses, and U-Pb zircon, Rb-Sr, and Ar-Ar geochronology. A Fairway Ridge cobbly sandstone has a ∼95 Ma (early Late Cretaceous) depositional age; two biotite granite cobbles are 111 ± 1 and 128 ± 1 Ma in age, and some volcanic pebbles are also likely Early Cretaceous. Fairway Ridge basalts have intraplate alkaline chemistry and are of Late Eocene age (∼40–36 Ma). By analogy with South Zealandia, we interpret strong positive continental magnetic anomalies of North Zealandia to mainly result from Late Cretaceous to Cenozoic intraplate basalts, many of them rift-related lavas. A new basement geological map of North Zealandia shows the position of the Mesozoic Gondwana magmatic arc axis (Median Batholith) and other major geological units. This study completes onland and offshore reconnaissance geological mapping of the entire 5 Mkm 2 Zealandia continent.

Tectonics↗

USGS2024

Compositional and structural mapping of Northwest Africa 15507 angrite

Angrite meteorites represent interesting sampling of planetary crustal environments. Quench-textured angrites with strong crystal zoning originated from the shallow surface region, with evidence of reducing conditions during solidification. Plutonic angrites have more coarse-grained igneous and metamorphic textures with comparatively less zoning and are interpreted as having equilibrated at greater depth. Plutonic angrites contain the minerals magnetite Fe 3 O 4 comprised of Fe 2+ and Fe 3+ , and rhönite, where Fe 3+ is required by inspection of mineral stoichiometry. NWA 15507 is a plutonic angrite with a microgabbroic texture (mean grainsize ∼1.4 mm) composed of zoned Al-Ti-augite, Ca-bearing olivine, anorthite, with accessory kirschsteinite, rhönite, hercynite, low-Ni kamacite, merrillite, Ti-free magnetite and troilite [ 1 ]. Rhönite analyzed by electron-probe microanalysis (EPMA) has the formula Ca 2.04 (Mg 0.32 Fe 2+ 4.25 Fe 3+ 0.47 Ti 0.33 Al 0.61 )(Si 3.74 Al 2.26 )O 20 , where Fe 3+ was estimated by stoichiometric analysis. During preliminary cathodoluminescence (CL) analysis, complex subgrain and oscillatory zoning was observed in the anorthite. In this study we use a combined approach of electron probe microanalysis (EPMA), cathodoluminescence (CL), electron-backscatter diffraction (EBSD), and laser ablation inductively-coupled plasma mass spectrometry (LA-ICP-MS) to further investigate the anorthite and distribution of Fe 3+ in NWA 15507.

Microscopy and Microanalysis↗

USGS2025

Characterization and validation of tidally calibrated strains from the Alto Tiberina Near Fault Observatory Strainmeter Array (TABOO-NFO-STAR)

Six horizontal borehole tensor strainmeters (TSM1-6) installed from Fall 2021 to Spring 2022 comprise the Alto Tiberina Near Fault Observatory Strainmeter Array (STAR), providing an unprecedented opportunity to investigate seismic and aseismic deformation from hazardous high- and low-angle normal faults in Italy. Prior to use in tectonic applications, they require in-situ calibration and correction for non-tectonic signals. We tidally calibrate the instruments, characterize the calibration uncertainty, and test the results against environmental and earthquake signals originating from local to teleseismic distances. The STAR sites demonstrably deviate from assumptions common to the standard manufacturer's calibrations, including negative areal coupling at TSM3-6. While the tidally calibrated strains have ~3-56% uncertainty, the calibrated dynamic strains show interstation precision and accuracy to nanostrain levels, and static coseismic offsets in the array footprint are within uncertainty. TSM3 records a complex series of strains that may arise from dynamically triggered near-borehole fracture slip and fluid flow that does not appear to affect its sensitivity to lower strain rate deformation. Future calibration improvement may be afforded with longer stable timeseries, particularly for TSM4. Overall, our analyses demonstrate expanded geodetic capability for detecting deformation in the Alto Tiberina Near Fault Observatory.

Alto Tiberina Near Fault Observatory↗

USGS2024

On the provenance of field reports of the 1886 Charleston, South Carolina, earthquake: A seismo-historical whodunnit

Much of what is known about the effects of the 1886 Charleston, South Carolina, earthquake throughout the epicentral region can be attributed to meticulous field investigations by an individual with training in geology and engineering, Earle Sloan ( Clendenin, 1926 ). In a recent study, Bilham and Hough (2024) undertook a detailed analysis of the effects of the earthquake on railroads in the Charleston region, drawing heavily from Sloan’s reports. This exercise identified several inconsistencies in Sloan’s field reports, including understandable measurement imprecision, inferred data entry mistakes, and transcription errors. The study also begged the question, where was Sloan at the time of the mainshock and over the following week? And to what extent did he draw from secondhand information in compiling his reports? On this question Sloan’s reports were sometimes enigmatic, lending themselves to misinterpretation in contemporaneous as well as modern interpretations. Beyond the details that were germane for, and briefly summarized by, the studies of Bilham and Hough (2023 , 2024) , in this report we don our historical seismologist caps to chronicle Sloan’s activities following the earthquake. We summarize our inferences here for the benefit of future scholars who might attempt to retrace either Sloan’s footsteps or our own. This study also serves to highlight Sloan’s singular contributions to earthquake science, which were never published separately.

South Carolina↗

USGS2026

Groundwater dependency and hydroclimatic influences on riparian and upland vegetation productivity, Upper San Pedro, Arizona, United States

In arid and semi-arid regions, groundwater sustains vegetation through subsurface water access, yet the responses of groundwater-dependent ecosystems (GDEs) to changing hydroclimate and groundwater availability are relatively understudied. This study investigates seasonal and spatial patterns in vegetation greenness using Landsat Enhanced Vegetation Index (EVI) values across riparian and upland zones in the semi-arid Upper San Pedro (USP) watershed, southern Arizona, which experiences a bimodal precipitation regime. We paired 25 years (2000–2024) of EVI and depth to groundwater (DTG) data from 89 wells and climate metrics (precipitation and vapour pressure deficit) to quantify the sensitivity of vegetation to subsurface moisture as well as atmospheric moisture supply and demand. Vegetation at wells near the USP riparian area showed strong associations between EVI and DTG anomalies during the monsoon season, indicating sustained groundwater use even during this wet period when summer precipitation is abundant. In contrast, upland vegetation that lacked access to groundwater showed minimal sensitivity in EVI to DTG and was generally less responsive to vapour pressure deficit. Interestingly, the riparian GDEs were not decoupled from precipitation and climate variability. These results underscore the importance of groundwater for maintaining riparian productivity and highlight the utility of remote sensing in identifying vegetation-climate-groundwater linkages across heterogeneous dryland landscapes.

Arizona↗

USGS2025

Depth of magma crystallization and fluid exsolution beneath the porphyry-skarn Cu deposits at Santa Rita and Hanover-Fierro, New Mexico, USA

The depth level at which porphyry Cu–forming magmas fractionated and exsolved mineralizing fluids is actively debated. In the classic model, extensive magma fractionation occurs in large, upper crustal magma chambers, and concomitant fluid exsolution leads to forceful expulsion of residual magmas in the form of porphyry dikes, stocks, and breccia pipes, which subsequently serve as pathways for the mineralizing fluids. In contrast, some recent studies highlighting the role of deep crustal magma fractionation in the production of fertile magmas essentially deny the existence of upper crustal magma chambers at the time of mineralization. To address this, we conducted a detailed thermobarometric investigation of 13 intermediate to felsic, porphyritic intrusive rocks related to porphyry-skarn Cu mineralization at Santa Rita and Hanover-Fierro, New Mexico, United States, representing two premineralization magmas (61–60 Ma), seven synmineralization magmas (60–58 Ma), and four late- to postmineralization magmas (58–57 Ma). For each sample, the pressure of last magma crystallization before final magma ascent to the current exposure level was reconstructed based on Al-in-hornblende barometry of small hornblende inclusions trapped within quartz phenocrysts and through titanium-in-quartz (TitaniQ) thermobarometry of the host quartz phenocrysts themselves. Since quartz is one of the last crystallizing magmatic minerals, and no significant phenocryst growth could have occurred in small dikes and stocks after final magma emplacement, quartz phenocrysts and their contained hornblende inclusions record the depth of last magma crystallization before final magma ascent. When present, hornblende phenocrysts and hornblende inclusions within other major phenocrysts were also analyzed. Both quartz and hornblende barometers return consistent average pressures of 3.2 ± 0.4 kbar for the entire suite of pre- to postmineralization magmas, corresponding to depths of 11 to 14 km. The synmineralization magmas return even more consistent average pressures of 3.1 ± 0.2 kbar, corresponding to a depth of 12 ± 1 km. The volume of the mineralizing porphyry dikes and stocks at the emplacement level is far too small to have provided all the fluids and metals required to form the observed ore deposits. Therefore, the majority of the ore-forming fluids must have originated from the magmas that crystallized at 12 ± 1 km depth. The ore deposits, conversely, formed at ~5-km paleodepth. This implies that most of the mineralizing fluids traveled an average vertical distance of ~7 km from their magmatic source to the eventual site of ore precipitation. The relatively unaltered nature and low veining degree of deeper parts of mineralized porphyry dikes and stocks suggest that the fluid transport through these intrusive bodies occurred mostly at near-solidus conditions by means of fluid percolation along grain boundaries. In summary, our results suggest that (1) a large, upper crustal pluton exists ~7 km beneath the Santa Rita and Hanover-Fierro deposits; (2) abundant phenocryst crystallization occurred at this depth level; and (3) this pluton was the main source for the exsolution of ore-forming fluids. However, the investigated rocks have elevated whole-rock Sr/Y ratios, indicating magma fractionation at deep crustal levels. As a result, our preferred model is a combination of the two end-member models introduced above, with most magma fractionation having occurred in the deep crust and with residual, intermediate to felsic melts having ascended and accumulated at 11 to 14 km paleodepth, where they continued to crystallize with comparatively little crystal-liquid separation, before some of these magmas ascended further to shallow levels and quenched to porphyries.

New Mexico↗

USGS2025

The effects of line simplification on planform geometry

Data on maps should retain accuracy regardless of scale. Yet, as cartographic lines are generalized, there can be impacts on properties such as topology, density, and planform geometry. Here, we investigate the use of the Scale Specific Sinuosity (S3) metric (Stanislawski et al., 2023) to evaluate the effects of line simplification on planform geometry, which is the bends of streams in map view. We employ an open-source Python S3 workflow to characterize the geometry of five diverse stream channels in the United States. The original data are extracted from the U.S. Geological Survey National Hydrography Dataset 1:24,000-scale vector data (U.S. Geological Survey, 2000) (Table 1), and the simplification is done using the Visvalingam and Whyatt method (2017) with a simplification tolerance of 0.5, 1.0, 1.5, and 2.0 km. The S3 analysis is calculated at each level of simplification and S3 derivatives are generated. Derivatives include measures of sinuosity, fractal dimension, and the dominant bend wavelength. The findings show that the change in planform geometry is scale-dependent, though simplification will have little effect on straighter lines. The change becomes more apparent in complex lines as the degree of simplification aligns with the scale of the dominant bend geometries. These logical conclusions are evidence that the S3 is a useful metric for automated characterization of bend geometry regardless of line complexity.

Conference Paper↗

USGS2024

Deep lithospheric controls on surface deformation and seismicity around the East Anatolian Fault Zone and A3 Triple Junction

The East Anatolian Fault Zone (EAFZ) is a plate-bounding strike-slip fault capable of hosting large earthquakes, as demonstrated by the extremely damaging February 2023 M w 7.8 and M w 7.7 mainshocks of the Kahramanmaraş earthquake sequence. Deformation related to this boundary, part of the Anatolia-Arabia-Africa (A 3 ) Triple Junction, is diffuse, as was shown by part of this earthquake sequence occurring on a northern splay of the EAFZ (the Sürgü-Çardak Fault Zone; SCFZ). Controls on surface deformation are commonly linked to stress in the brittle upper crust, but the complex deformation and seismicity patterns in this region may also reflect deeper processes, such as variations in the location and extent of the strong Arabian Plate lithospheric mantle. Seismic tomography indicates that the Arabian Plate underthrusts Anatolia as far north as the SCFZ and extends as far west as the central Adana Basin, coincident with a zone of relatively deep (>30 km) strike-slip seismogenesis that has produced M w >6 earthquakes. By investigating the relationship between deformation since the inception of the EAFZ (ca. 5 Ma), seismic structure, and seismicity, we infer that the SCFZ will become the future SE boundary of the Anatolian Plate as part of the evolving A 3 Triple Junction.

GSA Today↗

USGS2024

Integrating Sr isotopes, microchemistry, and genetics to reconstruct Salmonidae species and life history

Recent approaches to fisheries research emphasize the importance of the coproduction of knowledge in building resilient and culturally mindful fisheries management frameworks. Despite widespread recognition of the need for Indigenous knowledge and historical reference points as baseline data, archaeological data are rarely included in conservation biology research designs. Here we propose a novel multiproxy method to learn from former fisheries stewards by generating archaeological data on past salmonid population parameters. We used a newly developed, high throughput qPCR (HT-qPCR) chip, originally designed for environmental DNA (eDNA), for species identification of archaeological salmonid vertebrae. We combine this with the laser ablation split-stream (LASS) approach to identify ocean-migration versus freshwater residency. We test this multidisciplinary approach using both contemporary and archaeological salmonid samples and new radiocarbon dates from the Tronsdal Site on the Skagit River, Washington State, USA. This is a useful approach for extracting information about Salmonidae species and life history diversity from archaeological remains to reconstruct historic baselines for several population parameters in anadromous species with long periods of freshwater residency. The approach outlined in this paper may be particularly useful for research investigating past fisheries dynamics, offering hundreds to thousands of years of temporal depth for modern fisheries management, harvest policies, restoration ecology, and conservation biology.

Idaho, Oregon, Washington↗

USGS2024

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

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

Conference Paper↗

USGS2025

Examining the role of elevated and sustained strain in dynamically triggering earthquakes on the Anza section of the San Jacinto fault

Microearthquakes can be dynamically triggered in southern California by remote earthquakes. However, directly connecting dynamic triggering mechanisms with observational data remains challenging. One proposed failure mechanism suggests that both the amplitude and duration of cyclic fatigue caused by the passing seismic wave contribute to triggering occurrence. Here, we measure dynamic strains recorded by borehole strainmeters in the Anza section of the San Jacinto fault zone from 710 earthquakes that occurred over 300 km away between 2008 and 2017 to systematically investigate the role of elevated and sustained strain in controlling dynamic triggering. We design a suite of tests to evaluate whether specific amplitude thresholds and durations of strain can predict dynamic triggering cases. We further test whether the peak dynamic strain (PDS) can predict triggering occurrence in combination with the strain amplitude and duration. Based on these tests, there is no strain amplitude–duration threshold that can distinguish triggering occurrence in Anza. Dynamic triggering is more likely to occur if a remote earthquake causes a PDS above 100 nanostrain, though many cases were triggered at smaller PDSs. The lack of clear correlation between triggering and characteristics of the dynamic strain field suggests that the tested features of the incoming waves do not determine triggering occurrence and local fault conditions and slip processes are more important in controlling dynamic triggering in Anza.

California↗

USGS2025

Review and synthesis of the applications of machine learning to coalbed methane recovery

Over the last 30 years, a substantial literature has evolved on the use of machine learning (ML) to assess, predict, and improve the efficiency of coalbed methane (CBM) recovery. In the United States, the production of CBM declined as shale gas production matured, but CBM continues to be an important energy resource in other parts of the world. ML applications that have the potential to improve CBM reservoir management and production forecasts, and to increase exploration and operational efficiency, are still of significant interest. The integration of geostatistical techniques into the CBM ML applications has been largely absent but represents an opportunity for improvement. The literature demonstrates the widespread interest in, and applicability of, ML algorithms applied to CBM problems, and that they continue to result in improvements in predictive performance. However, (1) much of the research is more academic than operational, (2) many results are based on simulations, or small or proprietary datasets, (3) ML performance information can be inconsistent and sometimes entirely omitted, (4) most methodologies are unique to the specific CBM situation and likely not generalizable, (5) no standard data repositories are available to directly compare the performance of competing algorithms, and (6) the spatial component is often omitted. Finally, relatively new ML protocols involving causality analysis and reinforced learning, as well as hybrid workflows combining both supervised and unsupervised learning, are anticipated to dominate the future investigations. Integration of geostatistical and geospatial analysis with ML should enhance performance.

Book chapter↗

USGS2026

Decoding the oxidative digestion mechanism for polystyrene nanoplastic detection in the Great Lakes using a customizable Raman spectral processing algorithm

Despite the concerns surging around nanoplastics (NPs) regarding their prevalence and bioavailability in freshwater systems, robust detection of NPs in complex environmental matrices is hindered by the lack of standardized sample pretreatment and a mechanistic understanding of oxidative digestion. Here, we systematically investigate the interaction between hydrogen peroxide (H 2 O 2 ) and polystyrene (PS) NPs during digestion in deionized (DI) water and four environmental matrices from in and around the Great Lakes Basin. To facilitate high-throughput analysis, we develop Pre_peak, a customizable Raman spectral processing algorithm that achieves >99% accuracy for both NP identification and interference rejection, allowing reliable NP quantification via pixel counting and systematic decoding of the oxidative digestion mechanisms. In DI water, varying H 2 O 2 doses from 0 to 30% has negligible effects on the recovery and Raman signal intensity of PS NPs over 24 hours of digestion. However, morphological changes and aggregation of PS NPs are observed when the H 2 O 2 dose exceeds 20%. Prolonged digestion further leads to progressive NP loss. In natural waters, the optimal dosage and digestion duration depend on matrix characteristics, including dissolved organic matter (DOM) and ion composition. This study provides mechanistic insights into NP–oxidant interactions and underscores the need for matrix-tailored digestion protocols to advance standardized NP detection in freshwater environments.

Great Lakes↗