USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Contaminant Hydrology”

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,659 records · Page 70Linked to original sources
USGS2026

GST-1: A high-resolution global sediment thickness model

Global Sediment Thickness 1 (GST-1) is a high-resolution sedimentary thickness model calculated on a 0.125° x 0.125° grid. It modifies the sediment thickness of the 1° x 1° Earth Crustal Model 1 (ECM1) by means of 3D inversions of free air gravity anomalies. GST-1 is calculated by performing structural inversions on high-density contrasts across two crustal boundaries: the sediment – basement interface and the crystalline crust – upper mantle interface. The inversions are calculated in each of ten overlapping 3D models that span the globe. These ten models are merged to obtain the GST-1 global model, providing an eight-fold increase in lateral spatial resolution in comparison with ECM1 and CRUST 1.0. Our sediment thickness model exploits the nearly continuous sampling of gravity data when compared to the irregular, sparse sampling of seismic refraction data. Sediment thickness values in GST-1 are in excellent agreement with independently derived cross sections from well-studied sedimentary basins, and within expected resolution limits of seismic refraction data. GST-1 offers a robust, high resolution global model of sedimentary thickness to support studies of sedimentary basins.

Tectonophysics↗
USGS2025

Haunted Summerville: Ghostly lights or earthquake lights?

Among the colorful local lore in the Charleston, South Carolina, area, are a number of ghost stories, shared not only over campfires but also in published books. Among the most well-known of the stories is the tale of the Summerville Light. Local lore holds that a strange light sometimes seen in a remote area is a lantern carried by the ghost of a local woman who once waited hours for her husband, who turned out to have been decapitated earlier that day in a train accident (DePoppe, 2023). Extant sources suggest the ghost stories began to circulate in the 1950s to 1960s. So pervasive was the lore that (Old) Sheep Island Road became known among local residents as Light Road, with a local stretch of road known today as Old Light Road. Reviewing the location where the lights appear as well as the nature of accounts, I suggest that many if not all of the anecdotal observations can be most readily attributed to natural phenomena, including earthquake lights from earthquakes that were too small to be felt. Accounts of lights near Summerville cluster in proximity to the generally accepted epicenter of the 1886 Charleston, South Carolina, earthquake, where foreshocks to the 1886 mainshock were apparently concentrated, and within a few kilometers of three M3.5 – 4.4 earthquakes in 1959 and 1960.

South Carolina↗
USGS2026

Open-source gravity reduction workflows for geothermal resource assessment

Potential-field geophysical data such as gravity can enhance understanding of geothermal resources at all stages of the resource life cycle, including assessment, exploration, development, and monitoring, and at multiple scales, from the reservoir scale to regional scale. However, to make gravity data useful for geothermal resource characterization, several processing steps are required to isolate the effects of density variations in the Earth’s crust to enable the identification of structural features associated with geothermal resources. Although this process is well-established, standard computational implementations for processing gravity data that are FAIR (Findable, Accessible, Interoperable, and Reproduceable) are still lacking. This paper details ongoing efforts at the U.S. Geological Survey (USGS) to develop a standard set of open-source Python tools for gravity data reduction that align with the FAIR principles. This workflow makes use of existing open-source tools for geophysical data processing with the goal of maximizing opportunities for rapid improvements, interoperability, and adaptability to other types of geophysical data.

Conference Paper↗
USGS2024

Monitoring marine eruptions

Introduction Submarine volcanoes produce much of the same seismicity and eruptive activity as subaerial volcanoes and can pose hazards to society. Although they can be monitored with similar techniques and methods as described in other chapters of this volume, their submerged location brings unique challenges. This chapter addresses these challenges and provides recommendations for monitoring volcanoes fully or partly in marine environments to meet the capabilities described in other chapters of this volume. The United States and its territories host dozens of submarine volcanoes with most (around 60) in the Commonwealth of the Northern Mariana Islands. Approximately 20 of the Northern Mariana Islands submarine volcanoes are known to be hydrothermally active, and 10 have confirmed eruptions since the 1950s (for example, Baker and others, 2008; Tepp and others, 2019a). Nine of those volcanoes were considered by the National Volcanic Threat Assessment (Ewert and others, 2018) to have a combination of eruptive type and summit depth that poses a higher risk of hazardous eruptions, although only one was listed as a moderate (level 3) threat. Other notable submarine volcanoes of interest to the United States that have historically erupted are Axial Seamount off the Washington State coast, Kamaʻehuakanaloa in Hawaiʻi, and Vailuluʻu seamount in American Samoa. All of these, however, have a low risk of hazards because of their depth (greater than 600 meters below sea level) and eruptive type and so are not included in the National Volcanic Threat Assessment. In addition to submarine volcanoes, the submerged flanks of island volcanoes can also be a source of hazardous submarine eruptions—for example, the 1877 eruption of Mauna Loa, Hawai‘i, in Kealakekua Bay (Wanless and others, 2006). The most notable submarine eruption in recent times was the 2022 eruption of Hunga Tonga–Hunga Haʻapai in Tonga, which was one of the largest eruptions on Earth in the past 100 years. It created a massive volcanic plume, unprecedented shock waves, and far-reaching tsunami (Lynett and others, 2022). Other recent submarine eruptions in the Pacific Ocean Basin have produced subaerial plumes that reached aircraft heights (Carey and others, 2014) and large pumice rafts that can affect marine traffic and harbors (for example, Jutzeler and others, 2014; Kornei, 2019). These examples illustrate the potential hazards of major submarine eruptions. Yet, submarine volcanoes are largely unmonitored, and many eruptions occur that are unnoticed or only identified hours or days afterward. Within U.S. territory, submarine volcanoes in the Northern Mariana Islands have been known to produce eruptive activity that can affect society. Reports from fishermen and other marine vessels in the Northern Mariana Islands have noted underwater explosions, sea-surface discoloration, and bubbling water, all of which are known to be signs of submarine volcanic activity. South Sarigan seamount, located about 160 kilometers (km) north of Saipan, erupted in 2010 from greater than 150 meters below the sea surface, resulting in a gas and ash plume that reached more than 11.9 km into the atmosphere (for example, Searcy, 2013; Embley and others, 2014), high enough to affect international air traffic. Precursory and co-eruptive seismicity was detected on the regional Northern Mariana Islands seismic network (Searcy, 2013) and on global monitoring instruments (Green and others, 2013). Monitoring of submarine volcanoes is best accomplished with marine-based instrumentation, which is also useful for monitoring small island volcanoes that may not have the land area necessary for comprehensive subaerial monitoring. The primary marine-based instrumentation used for submarine volcanoes includes ocean-bottom pressure sensors to assess sea-floor deformation, ocean-bottom seismometers (OBSs) to detect seismicity, and both moored and ocean-bottom hydrophones to detect submarine explosions. Other sensors offer important monitoring data, such as turbidity, temperature, and chemistry of hydrothermal emissions. Marine-based instruments are typically deployed in campaign-style networks with no real-time telemetry owing to cost considerations and technical limitations. However, when necessary, marine instruments can be operated in real time using cables to transmit data to land-based facilities; other technologies for this purpose are in use or in development, such as acoustic transmission from the instrument to a moored buoy (Matsumoto and others, 2016) and a winch-based system with a satellite antenna that is part of the instrument mooring (Matsumoto and others, 2019). Emerging technologies for marine-based monitoring may be considered as part of a future monitoring plan. These technologies include ocean gliders and floats with on-board hydrophones that have been used to record earthquakes and submarine eruptions (for example, Matsumoto and others, 2013; Sukhovich and others, 2015) and fiber-optic cables that have been used as strainmeters to detect earthquakes (for example, Marra and others, 2018; Lindsey and others, 2019). Land-based instruments and satellites can also provide some capability for monitoring submarine volcanoes, but they provide more limited observations than marine-based instrumentation.

Scientific Investigations Report↗
USGS2026

Evaluation of stream capture related to groundwater pumping, middle Humboldt River Basin, Nevada

Historical, future, and potential stream capture from groundwater pumping in the middle Humboldt River Basin (MHRB), Nevada, is estimated using a calibrated numerical groundwater flow model. The model was developed to estimate (1) stream capture, which is the change in flux between the groundwater system and the Humboldt River and tributaries, and (2) change in streamflow, which is the change in streamflow estimated for the Imlay gage on the Humboldt River (U.S. Geological Survey streamgage 10333000). Historical stream capture for water years (WYs) 1961–2015 is estimated using recorded and estimated groundwater pumping during that period. Future (predictive) stream capture was based on historical stresses (WYs 1961–2015) using a scenario that simulated non-mine pumping from WY 2015 at a uniform rate for 100 years into the future. Potential stream capture throughout the middle Humboldt River Basin from groundwater pumping during varying durations of time are presented in a series of capture maps. Maps also are presented that show the potential to capture from groundwater evapotranspiration, as well as the storage changes for pumping duration of 100 years. Estimates of historical stream capture from the mainstem Humboldt River during the early 1960s are less than 400 acre-feet per year (acre-ft/yr) when groundwater withdrawals and pumping rates were relatively small compared to more recent times. In the late 1980s and early 1990s, groundwater withdrawals increased and estimated historical stream capture also increased from about 4,000 acre-ft/yr in the late 1980s and early 1990s to as much as 18,800 acre-feet (acre-ft) in WY 1998. In WY 2015, estimated historical stream capture declined to about 13,000 acre-ft because of decreasing groundwater withdrawals and lower streamflow during the drought of WYs 2012–15, resulting in less stream water available for capture. Stream capture was estimated for 100 years into the future based on WY 2015 non-mine pumping rates and mine-dewatering activity through WY 2015. Stream capture is forecast to increase to about 23,000 acre-ft/yr, and streamflow in the Humboldt River could decrease by as much as 19,000 acre-ft/yr. Pumping for mine-dewatering and the associated discharge of that water affects streamflow in the Humboldt River at Imlay, Nevada (U.S. Geological Survey streamgage 10333000). Historically, from WYs 1991 to 2015, streamflow was greater at Imlay gage during active mine-dewatering from mine-water discharge operations and increased by as much as 105,000 acre-ft in WY 1998. The increase was attributed mostly to the discharge of groundwater from mine-related dewatering operations directly into the mainstem Humboldt River or its tributaries, with some of this increase associated with return flows from discharge to rapid infiltration basins. Results indicate that streamflow at Imlay gage is expected to decrease by as much as 1,600 acre-ft/yr 30 years after mine-related pumping and discharge are discontinued. The streamflow reductions at the Imlay gage are expected to then decrease to around 500 acre-ft/yr, 100 years after mine-related pumping and discharge are discontinued. Potential capture maps were produced for pumping durations of 10, 25, 50, and 100 years. Capture map results indicate that areas of greater potential stream capture occur adjacent to the Humboldt River and for upstream tributaries areas north of the Humboldt River.

Nevada↗
USGS2025

Protected from Pterygoplichthys? Predicting thermal habitat suitability for nonnative armored catfish in the Suwannee River

Objective Nonnative fishes can modify ecosystems and harm economies when they are introduced to new environments. Climate change is likely to assist the spread and establishment of some nonnative fishes (e.g., warmwater species), but spatiotemporal gaps in water temperature monitoring and modeling may prevent ecologists and managers from forecasting thermal habitat suitability for these taxa. The purpose of this study was to develop a predictive model of winter water temperatures and thermal habitat suitability for two priority nonnative armored catfish, Vermiculated Sailfin Catfish Pterygoplichthys disjunctivus and Orinoco Sailfin Catfish P. multiradiatus , in the Suwannee River, Florida and Georgia. Methods Precipitation- and groundwater-corrected air–water temperature models were developed and evaluated using a model selection procedure to predict water temperatures at four sites in the Suwannee River. These models were chosen because they blend the simplicity of air–water temperature models with the accuracy of hydrometeorological models to create an efficient, economical, management-relevant approach for analyzing and forecasting water temperature. Results Most of the top-performing water temperature models (92%) had precipitation or groundwater corrections to air–water temperature formulations. Projected mean and maximum water temperatures increased as simulated climate change intensified. All four Suwannee River sites studied were projected to be thermally hospitable to the survival of Vermiculated Sailfin Catfish. Lower river sites, noticeably warmer than upper river sites, were conducive to the survival of Orinoco Sailfin Catfish throughout the winter months. The upper river sites were too cold for Orinoco Sailfin Catfish survival in some climate-change scenarios, but the Suwannee River has an abundance of constant-temperature springs that are likely hospitable to Vermiculated Sailfin Catfish and Orinoco Sailfin Catfish throughout the year. Conclusions The findings suggest that winter water temperatures will likely not be a barrier to the survival of Pterygoplichthys catfish in the Suwannee River, amplifying the importance of conservation and management approaches to inhibit their spread and establishment. If the Pterygoplichthys population remains small and isolated and decision makers are able to devote required staff time and resources to managing these species, removal and eradication at local if not broader scales may be reasonable goals. This study provides a water temperature modeling approach that can aid ecologists and managers in prioritizing sites to prevent the introduction, slow the dispersal, eradicate, and control Pterygoplichthys catfish and other nonnative fishes in the Suwannee River and beyond.

Florida, Georgia↗
USGS2025

Characterizing Meteor Crater impact melts through geochemistry and textural analysis

The U.S. Geological Survey Astrogeology Science Center houses the Meteor Crater sample collection, an assemblage of over 2,500 meters of cuttings from 161 drill holes into Meteor Crater’s rim, flanks, and ejecta blanket. We have utilized this unique collection to study the composition and spatial distribution of impact-generated materials from within the ejecta blanket. Meteor Crater has historically been known to have generated only a relatively small amount of impact melt compared to other terrestrial craters of similar size. A detailed compositional and textural dataset of impact-derived melts from this impact can therefore be a useful asset in improving our understanding of crater formation, and in particular impact melt formation. We have characterized 42 impact-melt particles from Meteor Crater using a scanning electron microscope and an electron microprobe for textural and compositional analysis. We analyzed samples from six drill holes in the ejecta blanket, situated to the northwest, southeast, south, and southwest of the crater (ejecta northeast of the crater is devoid of impact melts). Impact melts were collected from drill cuttings at various depths within the ejecta blanket, ranging from a few centimeters below the surface down to ~6.5 meters. Backscattered electron (BSE) images were acquired for each analyzed impact-melt particle. To characterize the various textures and phases present in each impact melt, we also took many detailed BSE images. Our geochemical analyses include full spectral profiles using energy dispersive X-ray spectrometry and well-calibrated wavelength dispersive spectrometry for a number of phases, including minerals (olivine, pyroxene, and so on), pristine glass, and metallic inclusions. The full dataset is available in ScienceBase as a data release (Gullikson and others, 2024), accessible at https://doi.org/10.5066/P9OGAJ8P . Our goal for this Open-File Report is to provide a summary of this immense dataset, details on data collection, descriptions of the different phases observed within impact-melt particles (both geochemically and texturally), and observable trends.

Arizona↗
USGS2025

Ultramafic float rocks at Jezero crater (Mars): Excavation of lower crustal rocks or mantle peridotites by impact cratering?

Based on observation and data from meteorites and in situ scientific missions, experiments as well as models, the Martian mantle is assumed to share some compositional and mineralogical affinity with the terrestrial mantle. However, there might be subtle differences like the Martian mantle being more ferroan. Yet, we do not have any direct analysis of a Martian mantle rock to confirm this assumption. NASA’s Perseverance rover found olivine-rich boulder-sized float rocks on the upper Jezero fan (Mars). These boulders have an ultramafic composition and their mineralogy is dominantly composed of Fo 73±3 olivine with high-Mg orthopyroxene, Cr-rich Ti-Fe oxides and minor plagioclase and high-Ca pyroxene. Microtextural and petrological analysis reveals that these minerals crystallized at equilibrium. In addition, these boulders are different from all the bedrocks analyzed by Perseverance along its traverse which are crustal igneous rocks and sediments. Comparing our data to Martian meteorites and available Mars bulk silicate models (BSM), we discuss that these boulders could represent primitive melts and/or lower crustal material, and we specifically hypothesize that they could be mantle peridotites. We propose that these putative mantle rocks could have been excavated by the succession of impacts from the shallow mantle or lower crust in the Isidis region where Jezero crater is located. These olivine-rich boulders could thereby constitute the first direct analysis of a Martian mantle rock.

Earth and Planetary Science Letters↗
USGS2024

Nitrogen deposition weakens soil carbon control of nitrogen dynamics across the contiguous United States

Anthropogenic nitrogen (N) deposition is unequally distributed across space and time, with inputs to terrestrial ecosystems impacted by industry regulations and variations in human activity. Soil carbon (C) content normally controls the fraction of mineralized N that is nitrified ( ƒ nitrified ), affecting N bioavailability for plants and microbes. However, it is unknown whether N deposition has modified the relationships among soil C, net N mineralization, and net nitrification. To test whether N deposition alters the relationship between soil C and net N transformations, we collected soils from coniferous and deciduous forests, grasslands, and residential yards in 14 regions across the contiguous United States that vary in N deposition rates. We quantified rates of net nitrification and N mineralization, soil chemistry (soil C, N, and pH), and microbial biomass and function (as beta-glucosidase (BG) and N -acetylglucosaminidase (NAG) activity) across these regions. Following expectations, soil C was a driver of ƒ nitrified across regions, whereby increasing soil C resulted in a decline in net nitrification and ƒ nitrified . The ƒ nitrified value increased with lower microbial enzymatic investment in N acquisition (increasing BG:NAG ratio) and lower active microbial biomass, providing some evidence that heterotrophic microbial N demand controls the ammonium pool for nitrifiers. However, higher total N deposition increased ƒ nitrified , including for high soil C sites predicted to have low ƒ nitrified , which decreased the role of soil C as a predictor of ƒ nitrified . Notably, the drop in contemporary atmospheric N deposition rates during the 2020 COVID-19 pandemic did not weaken the effect of N deposition on relationships between soil C and ƒ nitrified . Our results suggest that N deposition can disrupt the relationship between soil C and net N transformations, with this change potentially explained by weaker microbial competition for N. Therefore, past N inputs and soil C should be used together to predict N dynamics across terrestrial ecosystems.

contiguous United States↗
USGS2024

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri River between Kansas City and St. Louis, Missouri, May 19–26, 2021

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near nine bridges at eight highway crossings of the Missouri River between Kansas City and St. Louis, Missouri, from May 19 to 26, 2021. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches about 1,640 to 1,840 feet (ft) longitudinally and generally extending laterally across the active channel from bank to bank during low to moderate flood-flow conditions. These surveys provided channel geometry and hydraulic conditions at the time of the surveys and provided characteristics of scour holes that may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. These data also may be useful to the Missouri Department of Transportation as a low to moderate flood-flow assessment of the bridges for stability and integrity issues with respect to bridge scour during floods. Bathymetric data were collected around every in-channel pier. Scour holes were present at most piers for which bathymetry could be obtained, except those on banks or surrounded by riprap. Occasionally, scour holes were minor and difficult to discern from nearby dunes and ripples. All the bridge sites in this study were previously surveyed and documented in previous studies. Comparisons between bathymetric surfaces from the previous surveys and those of the current (2021) study do not indicate any consistent correlation between channel-bed elevations and streamflow conditions. The average difference between the bathymetric surfaces varied from 1.59 ft higher to 0.95 ft lower in 2021 than 2017, which corresponds to a gain of 100,200 cubic yards and a loss of 55,800 cubic yards, respectively. The average difference between the bathymetric surfaces varied from 2.74 ft higher to 3.05 ft lower in 2021 than 2013, which corresponds to a gain of 111,500 cubic yards and a loss of 169,200 cubic yards, respectively. The average difference between the bathymetric surfaces varied from 4.52 ft higher to 1.38 ft lower in 2021 than 2011, which corresponds to a gain of 221,100 cubic yards and a loss of 90,300 cubic yards, respectively. The most substantial overall net gain was 221,100 cubic yards between 2011 and 2021 at structures L0550 and A4497 at Jefferson City (site 20). The large net gain likely results from a combination of the mitigation of the scour holes near pier 4 of both bridges and the substantially lower flow in 2021 than in 2011. Alternatively, the most substantial overall net loss was 169,200 cubic yards between 2013 and 2021 at structure A6288 at Hermann (site 21), despite comparable streamflows. Pier size, nose shape, and skew to approach flow had a substantial effect on the size of the scour hole observed at a given pier. Larger and deeper scour holes were present at piers with wide or blunt noses caused by exposed footings or caissons. When a pier was skewed to primary approach flow, the scour hole was generally deeper and larger than at a similar pier without skew; furthermore, the shape of the scour hole near skewed piers in this study generally was longer and deeper on the side with impinging flow. At structure A6288 at Hermann (site 21), the scour hole near pier 5 was difficult to discern from nearby dunes and ripples, whereas the upstream edge of the footing was visible at pier 4, which likely contributes to the larger scour hole near that pier; the top of the footing may blunt the horseshoe vortex at pier 5, but the exposed front of the footing may exacerbate the vortex at pier 4.

Missouri↗
USGS2025

Peak streamflow trends in Montana and northern Wyoming and their relation to changes in climate, water years 1921–2020

Frequency analysis on annual peak streamflow (hereinafter, peak flow) is essential to water-resources management applications, including critical structure design (for example, bridges and culverts) and floodplain mapping. Nonstationarity is a statistical property of a peak-flow series such that the distributional properties (the mean, variance, or skew) change either gradually (monotonic trend) or abruptly (shift, step change or change point) through time. Not incorporating or accounting for observed nonstationarity into peak-flow frequency analysis might result in a poor representation of the true probability of large floods and thus misrepresent the actual flood risks to life and property. This report summarizes how hydroclimatic variability might affect the temporal and spatial distributions of peak-flow data in the State of Montana (and northern Wyoming) and is part of a larger study to document peak-flow nonstationarity and hydroclimatic changes across a nine-State region consisting of Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flows. This study was completed in cooperation with the Montana Department of Natural Resources and Conservation. The purpose of this report is to characterize temporal and spatial patterns of nonstationarity in peak flows and hydroclimatology in Montana and northern Wyoming. In this evaluation, peak-flow, daily streamflow, and model-simulated gridded climatic data were examined for monotonic trends, change points, and other statistical properties that might indicate changing climatic and environmental conditions. This report includes background information on the study area, the history of U.S. Geological Survey peak-flow data collection and frequency analysis in Montana, and the review of research relating to hydroclimatic variability and change in Montana. This study might help provide a framework for addressing potential nonstationarity issues in peak-flow frequency updates that commonly are completed by the U.S. Geological Survey in cooperation with other agencies throughout the Nation. The analytical structure of this study includes analyses of monotonic trends and change points in numerous hydroclimatic variables in assigned 30-, 50-, 75-, and 100-year analysis periods. For Montana and part of Wyoming, the 30-, 50-, 75, and 100-year analyses included 157, 70, 48, and 12 streamgages, respectively. For those streamgages, nonstationarities were analyzed in the following variables: (1) climatic variables, including annual and seasonal (winter, spring, summer, and fall) temperature and precipitation; (2) daily streamflow variables, including the annual center of volume duration, annual center of volume median, and peaks over threshold with a mean of four events per year; and (3) annual peak-flow variables, including peak-flow timing and magnitude. A likelihood approach was used to express statistical confidence and assign the nonstationarity results as likely upward or downward (highest statistical confidence), somewhat likely upward or downward (less statistical confidence), or about as likely as not (little statistical confidence; hereinafter, neutral). For the nonstationarity analyses of the climatic, daily streamflow, and peak-flow variables, the results are presented in detail and discussed with respect to statewide patterns and geographic variability. For each of the 30-, 50-, and 75-year analyses, peak-flow change-point and monotonic trend analyses were compiled for streamgages classified with likely downward or likely upward trends. For those streamgages, the associated basin characteristics and nonstationarity results for peak-flow timing, daily streamflow, and climatic variables were investigated and statistically compared to discern associations among other variables that might contribute to the peak-flow nonstationarity results. The 50- and 75-year peak-flow nonstationarities identified in this study are mostly downward, in association with mostly upward temperature and potential evapotranspiration:precipitation monotonic trends. For the 50-, 75-, and 100-year analyses, the peak-flow change points are predominantly downward and are concentrated in the 1970s and 1980s, which indicates general consistency among the longer trend periods. These findings are in association with substantial research documenting globally rising temperature and atmospheric greenhouse gas concentrations that might be largely attributed to anthropogenic activities. Anthropogenic effects might represent long-term (on the order of several decades to more than a century) climate changes that might happen within highly variable natural climate fluctuations. Several paleo studies in the north-central United States have indicated that hydroclimatic extremes (that is, low- and high-streamflow conditions) before European settlement have been outside of extremes since the 1900s. Depending on the interactions of anthropogenic effects and natural climate variability, extreme high-streamflow conditions might occur in the future, even in the presence of long-term downward peak-flow trends.

Montana, Wyoming↗
USGS2024

Using the horizontal-to-vertical spectral ratio method to estimate thickness of the Barry Arm landslide, Prince William Sound, Alaska

Conducting detailed investigations of large landslides is difficult, especially in the subsurface, largely due to environmental factors such as steep slopes, difficult access, and numerous objective hazards. These factors have made it challenging to accurately estimate the depth to the failure surface of the Barry Arm landslide, a large (roughly 10 8 cubic meters), deep-seated bedrock landslide in Prince William Sound, Alaska, recognized in 2019. The landslide has exhibited accelerated movement in recent years and poses a potential tsunamigenic hazard if rapid failure occurs. Failure surface depth, equivalent to landslide thickness, is a necessary metric for landslide-volume calculations and associated tsunami wave models. In this report, we used seismic noise recorded by a seismometer located on the Barry Arm landslide in Alaska to calculate the horizontal-to-vertical spectral ratio (HVSR) to investigate the site fundamental frequency ( f 0 ) and depth of the failure surface. To ensure that observed peak frequencies in the spectral ratio were related to the underlying stratigraphy (and not caused by other noise sources like nearby glaciers, topographic resonance, weather, or human activities), we also calculated HVSRs using earthquake signals, HVSRs at other seismic stations within a 2.5-kilometer radius, and a standard spectral ratio between the landslide station and other sites. We observed multiple peaks in the landslide HVSR curves at 1.5 hertz (Hz), 4–5 Hz, and 7–11 Hz. The frequencies of these peaks were consistent at the landslide site through time and across methods and were dissimilar to those identified at other seismic stations in the area, making it unlikely the peaks were caused by local noise. Directional HVSRs calculated at 15-degree intervals showed amplification of the higher frequency peaks in the direction parallel to slip, indicating two-dimensional site effects. We used the distinct frequency peaks in the seismic record to develop a 4-layer conceptual model of the landslide wherein the top of the deepest layer represents the primary failure surface, or the boundary between damaged (mobile) and undamaged material. We inverted Rayleigh wave ellipticity curves within this 4-layer configuration with constraints on S-wave velocity and layer thickness based on analogous material properties identified in the literature. This was necessary absent any site-specific subsurface S-wave velocity data. The best-fitting models indicate a mean slope-normal depth to the failure surface of 188 (±9) meters (m), with additional stratigraphic boundaries at 4 and 20 m below ground surface, potentially representing layered motion. These results agree with and improve upon ranges estimated by previous studies and can support future modeling and assessment efforts at Barry Arm.

Alaska↗
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↗
USGS2026

Ground-motion aleatory-variability models for Puerto Rico and the U.S. Virgin Islands

I develop independent logic trees for aleatory variability for crustal and subduction-zone (interface and intraslab) earthquakes for seismic hazards analyses in Puerto Rico and the U.S. Virgin Islands (PRVI) from existing suites of ground-motion models (GMMs) and from ground-motion datasets, including a regional PRVI dataset. The aleatory variability models are parameterized using a partially nonergodic partitioning of standard deviation that consists of independently developed between-event ( ), site-to-site ( ), and event-corrected single-station ( ) standard deviation components. The effects of nonlinear site response on aleatory variability are incorporated through additional terms that modify the standard deviation components. Because one goal of this work is to develop independent logic trees for aleatory variability that synthesize the aleatory variability models from GMMs, I make use of the functional forms of the input GMMs. The PRVI dataset contains a limited number of stations with high-quality site metadata and does not contain records from earthquakes with magnitudes greater than 6.1, so I choose not to develop the aleatory variability models from the regional dataset alone. Instead, the standard deviation components from regional ground-motion data are evaluated against the components derived from GMMs and from available global datasets, and regionalized standard deviation components are incorporated where there is evidence that regional effects exhibit substantial differences. The resulting logic trees for aleatory variability consist of models of and that are consistent with semiempirical GMMs for active crustal and subduction-zone regimes, and two alternative models of , including one model that exhibits site-to-site variability informed by PRVI data, with values that exceed global models. The aleatory variability models may be considered in future hazards assessments in PRVI to simplify the hazard calculations, to incorporate regional ground-motion variability effects, and to enable direct logic-tree weighs of aleatory variability.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

From collision to subduction: Thermal-kinematic inversions constrain plate boundary structure and dip-slip activity in southwestern New Zealand

Constraints on three-dimensional fault system geometry and kinematics are required to understand both long-term tectonics and contemporary seismic hazard. In many cases, these parameters must be inferred from a combination of different data or analyses. Here, we invert low-temperature thermochronological data to quantitatively parameterize both the crustal-scale geometry and late Cenozoic dip-slip activity of the plate boundary system in southwest New Zealand. We find that the modern transpressive plate boundary system has accommodated the bulk of Australian-Pacific Plate convergence since establishment in the middle Miocene-early Pliocene. Shortening has occurred via (a) dextral-reverse slip on the moderate angle, listric central Alpine Fault décollement, and (b) reverse-slip on either the sharply curved, anti-listric northern Puysegur subduction interface or low-angle offshore thrust fault system that splays from the high-angle southern Alpine Fault at depth. We resolve late Cenozoic shortening rates on these structures that are indistinguishable, and interpret these to reflect the transfer of convergent deformation outboard of the predominantly strike-slip southern Alpine Fault. Our work helps elucidate the transition from strain localized oblique continental collision to strain partitioned subduction. We also provide an estimate of the across-strike width of the central-southern Alpine Fault section boundary at seismogenic depths. This is the first quantitative constraint on the geometry of a structural feature inferred to have arrested several major prehistoric earthquakes. This study furthers understanding of fault structure and dip-slip kinematics along the greatest sources of seismic hazard in and offshore of southwest New Zealand.

South Island↗
USGS2025

Using public participatory geographic information systems (PPGIS) to explore uses and values for Mojave Trails National Monument, California

Many people ascribe a variety of values to public lands and waters, but some values are more difficult to assess and quantify than others. Public participatory geographic information systems (PPGIS) are tools that have been used to help quantify and map the public’s diverse values for a landscape. This work describes the first known Office of Management and Budget–approved use of PPGIS by a Department of the Interior bureau. The U.S. Geological Survey developed an internet-based application to aid in gathering PPGIS data, called Values Mapping for Planning in Regional Ecosystems (VaMPIRE). Further, this work describes the first pilot of the VaMPIRE application in coordination with the Bureau of Land Management to collect spatial data and other survey data regarding the public’s uses of and values for locations within Mojave Trails National Monument. We emailed the link to the VaMPIRE application to an interested party email list in 2024 with 207 valid emails and received 74 responses; we also received 47 responses from members of an off-roading social media group. Of the list of 16 value options, recreation was the most popular value for the monument, followed by wilderness and inspirational. Over 1,000 points were placed throughout the monument, indicating locations people use or value, with the locations spread throughout the entire monument. Additionally, most survey respondents stated their ability to receive benefits in locations they mapped would not change in response to a hypothetical scenario related to recreational facility development. This report describes exploratory results from the first use of the VaMPIRE tool in Mojave Trails National Monument and includes reflections on how the process went and considerations for future use of VaMPIRE.

California↗
USGS2025

Grand Canyon landslide-dam and paleolake triggered by the Meteor Crater impact at 56 ka

This paper hypothesizes that the Meteor Crater impact in Arizona, USA, 56,000 years ago triggered landslides in Grand Canyon that dammed the Colorado River and formed Nankoweap paleolake. This is compatible with shock and earthquake physics for the impact that infer a M5.4 seismic event, attenuated to an effective magnitude of M3.5 at Grand Canyon. Results that support the hypothesis include radiocarbon dating of driftwood and luminescence dating of associated slack-water lake sediments that are preserved in caves up to 60 m above the modern Colorado River. Radiocarbon ages from two locations, including Stanton’s Cave, date the driftwood as 55.25 ± 2.44 ka (n = 4). Sediments associated with the driftwood gave a luminescence age of 56.00 ± 6.39 ka (n = 2). These six Grand Canyon dates, and three published ages for the Meteor Crater impact, show statistically indistinguishable results that support the hypothesis for a geologically instantaneous series of events with a mean age of 55.60 ± 1.30 ka. This work highlights the value of radiocarbon dating near the limits of the technique, integration of multiple dating methods, and seismic and landslide hazards associated with meteorite impacts in regions of extreme topography like Grand Canyon.

Arizona↗