SEARCH · USGS Science

Results for “Rapid Communications in Mass Spectrometry”

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,602 records · Page 59Linked to original sources

USGS2024

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Earth and Planetary Science Letters↗

USGS2025

Cooling perspectives on the risk of pathogenic viruses from thawing permafrost

Climate change is inducing wide-scale permafrost thaw in the Arctic and subarctic, triggering concerns that long-dormant pathogens could reemerge from the thawing ground and initiate epidemics or pandemics. Viruses, as opposed to bacterial pathogens, garner particular interest because outbreaks cannot be controlled with antibiotics, though the effects can be mitigated by vaccines and newer antiviral drugs. To evaluate the potential hazards posed by viral pathogens emerging from thawing permafrost, we review information from a diverse range of disciplines. This includes efforts to recover infectious virus from human remains, studies on disease occurrence in polar animal populations, investigations into viral persistence and infectivity in permafrost, and assessments of human exposure to the enormous viral diversity present in the environment. Based on currently available knowledge, we conclude that the risk posed by viruses from thawing permafrost is no greater than viruses in other environments such as temperate soils and aquatic systems.

mSystems↗

USGS2024

Seismicity zoning at Coso geothermal field and stress changes from fluid production and migration

The Coso geothermal field is a major geothermal power production site in the western United States. It has been observed that low-magnitude seismic events ( M < 3.71) are unevenly distributed in three distinct zones, namely, nearfield (<3 km), midfield (3–6 km), and farfield (> 6 km) from the Coso geothermal plant. These zones exhibit distinct changes in earthquake location before and during geothermal production episodes that began in 1986. After 1986, the midfield region of the main flank experiences a significantly lower seismicity rate than the surrounding areas before production episodes. During 2014–2019, the farfield earthquakes cluster in the eastern and western parts of the greater Coso area, which is discernably different from how those pre-production earthquake events were distributed along the conjugate NW-SE and SW-NW trending structures across the main flank. Here, we analyze the stage of stress with finite-element-based poroelastic simulations to illustrate how the spatiotemporal evolution of the seismicity is associated with the pattern of stress perturbations caused by fluid migration amid the operations of geothermal power plants. Generally, ∼70% of co-production seismicity is found in zones of increased Coulomb stress between 2014 and 2019 at >99% confidence. Meanwhile, the midfield zone of seismic paucity overlaps with the zone of decreasing pore-fluid pressure. Overall, the results provide a physical explanation of how decadal geothermal operations at Coso have perturbed stress-field changes and contributed to the evolving characteristic seismic pattern, shedding insights into assessing the seismic hazard in other geothermal settings.

California↗

USGS2026

Top Elevation of Glacial Till and Thickness of the Big Sioux Aquifer Delineated From Electrical Resistivity Tomography Surveys Near Sioux Falls, South Dakota, 2022 and 2025

The City of Sioux Falls, South Dakota, requested the U.S. Geological Survey perform electrical resistivity surveys on three parcels of land north of the city. Electrical resistivity data were collected along a total of 22 transects during March 14–18, 2022, and November 17–21, 2025. Results from electrical resistivity surveys were used to delineate the top of glacial till deposits for the purpose of characterizing the Big Sioux aquifer near the city. Delineating geologic contacts provides important information on groundwater storage, flow dynamics, well design and placement, contaminant transport, groundwater–surface-water interactions, and regional water modeling. The top elevation of glacial till and the thickness of the Big Sioux aquifer varied among the three survey areas. The interpreted top elevation of glacial till in the North survey area decreases from east to west toward a slough, with elevations ranging from 1,403 to 1,418 feet (ft). The estimated thickness of the Big Sioux aquifer in the North survey area increased from east to west, with thicknesses ranging from 23 to 38 ft. The top elevation of glacial till in the Well 72 survey area generally decreases from northwest to southeast. Top elevations of the glacial till in the Well 72 survey area ranged from 1,400 to 1,409 ft along the southern end of transect W72_2. The estimated thickness of the Big Sioux aquifer in the Well 72 survey area was greatest along a southeast to northwest trending channel, with thicknesses ranging from 28 to 40 ft. The top elevation of glacial till in the Nose survey area generally decreases west toward the Big Sioux River. Top elevations of the glacial till in the Nose survey area ranged from 1,362 to 1,395 ft. The estimated thickness of the Big Sioux aquifer in the Nose survey area ranged from 33 to 70 ft.

South Dakota↗

USGS2023

Explosive trachyte eruptions from the Al Efairia volcanic center in northern Harrat Rahat, Kingdom of Saudi Arabia

Harrat Rahat is an alkali basalt, continental, intraplate volcanic field located within the central-western part of the Kingdom of Saudi Arabia. The northern quarter of Harrat Rahat contains evolved volcanic products that achieve trachyte compositions (>60 weight percent SiO 2 ). Within the Al Efairia volcanic center, pyroclastic-flow and -surge deposits that reflect explosive trachyte volcanism (and minor exposed lava domes that reflect effusive trachyte eruptions) sit at the surface as the youngest expression of volcanic activity within this part of the Harrat Rahat volcanic field. Five trachyte deposits emplaced explosively have been identified within the Al Efairia volcanic center based on geologic mapping, petrography, geochemistry, and paleomagnetism. These units are the trachytes of Um Rgaibah, Gura 5, Gura 4, Al Efairia, and Al Qayf, in descending stratigraphic order. Here, we present 14 40 Ar/ 39 Ar analyses from four of these units, which yield eruption ages of 4.2±5.2 thousand years (ka) for the trachyte of Um Rgaibah, 79.7±1.6 ka for the trachyte of Gura 5, 84.3±1.6 ka for the trachyte of Gura 4, and 88.0±1.8 ka for the trachyte of Al Efairia. The eruption age of the trachyte of Al Qayf has been constrained to between 410.3±3.4 and 418.8±1.9 ka using paleomagnetic correlations and 40 Ar/ 39 Ar ages from overlying and underlying intermediate composition lava flows. Most of these trachytes have distinct geochemical compositions, petrographic characteristics, and directions of remanent magnetization. The exceptions are for the trachytes of Gura 4 and Gura 5, which overlap in their geochemical, petrographic, paleomagnetic, and geochronologic affinities. Based on these similarities, we interpret the trachytes of Gura 4 and Gura 5 to have erupted during a closely spaced (a few decades) time interval from the same magma batch but from craters that are >2 kilometers (km) apart. The eruption of the trachyte of Al Efairia at 88.0±1.8 ka is the result of a different magma batch that erupted a few thousand years prior to the trachytes of Gura 4 and Gura 5. The Al Efairia volcanic center is remarkably different from the Matan volcanic center located ~10 km to the north, which has also erupted young (<150 ka) trachytes. The Matan volcanic center has been shown to produce trachyte compositions only after eruption of basalt followed by intermediate lava flows, whereas only trachyte compositions have erupted within the Al Efairia volcanic center over this same time interval.

Al Efairia volcanic center, Harrat Rahat↗

USGS2025

Permafrost–wildfire interactions: active layer thickness estimates for paired burned and unburned sites in northern high latitudes

As the northern high-latitude permafrost zone experiences accelerated warming, permafrost has become vulnerable to widespread thaw. Simultaneously, wildfire activity across northern boreal forest and Arctic/subarctic tundra regions impacts permafrost stability through the combustion of insulating organic matter, vegetation, and post-fire changes in albedo. Efforts to synthesis the impacts of wildfire on permafrost are limited and are typically reliant on antecedent pre-fire conditions. To address this, we created the FireALT dataset by soliciting data contributions that included thaw depth measurements, site conditions, and fire event details with paired measurements at environmentally comparable burned and unburned sites. The solicitation resulted in 52 466 thaw depth measurements from 18 contributors across North America and Russia. Because thaw depths were taken at various times throughout the thawing season, we also estimated end-of-season active layer thickness (ALT) for each measurement using a modified version of the Stefan equation. Here, we describe our methods for collecting and quality-checking the data, estimating ALT, the data structure, strengths and limitations, and future research opportunities. The final dataset includes 48 669 ALT estimates with 32 attributes across 9446 plots and 157 burned–unburned pairs spanning Canada, Russia, and the United States. The data span fire events from 1900 to 2022 with measurements collected from 2001 to 2023. The time since fire ranges from 0 to 114 years. The FireALT dataset addresses a key challenge: the ability to assess impacts of wildfire on ALT when measurements are taken at various times throughout the thaw season depending on the time of field campaigns (typically June through August) by estimating ALT at the end-of-season maximum. This dataset can be used to address understudied research areas, particularly algorithm development, calibration, and validation for evolving process-based models as well as extrapolating across space and time, which could elucidate permafrost–wildfire interactions under accelerated warming across the high-northern-latitude permafrost zone. The FireALT dataset is available through the Arctic Data Center ( https://doi.org/10.18739/A2RN3092P , Talucci et al., 2024).

Arctic↗

USGS2025

3D Converted wave reverse time migration imaging

We describe a newly developed method for recovering high-resolution images of seismic discontinuities, such as subducting slabs, in 3D. Our method makes use of converted S → P or P → S waves observed by dense arrays of seismometers to infer the locations and relative strengths of seismic discontinuities at depth in a target region. Observed direct and converted waves are backpropagated to their times of origin. The time-reversed wavefield is then separated into its constituent P and S components via the Helmholtz decomposition, and those separated wavefields are used to compute imaging functions that characterize the locations and relative strengths of seismic discontinuities. Imaging functions may be designed to use either S → P or P → S waves, so that users can target those arrivals expected to be most dominant in a given dataset. We have previously demonstrated the efficacy of our method in two dimensions, and we now present a 3D implementation of our technique which addresses the significant computational challenges posed by the size of volumetric wavefield data in three dimensions. Through a series of synthetic examples, we demonstrate that our method is capable of recovering the fine scale structure of a subducting slab given realistic station coverage and earthquake sources. We investigate optimal seismic station geometries for our technique and explore image interpretability in regions with poor data coverage. We find that linear station geometries yield more optimal, interpretable imaging functions than collections of small arrays can. We also show that our method can successfully recover bothS → P or P → S images when realistic shear earthquake sources are used, and we explore the additional computational challenges presented by the high frequency content of S waves. Our results demonstrate the potential for our technique to recover high-resolution information about subducting slabs in real-world regions, given that relatively sparse seismic arrays with only approximately 100 stations are capable of recovering interpretable imaging functions from just a few realistic earthquake sources for multiple discontinuities at significant depth in an area of approximately 400~sq~km.

Seismica↗

USGS2025

U.S. Geological Survey—Department of the Interior, Region 11, Alaska—2023–24 biennial science report

Introduction U.S. Geological Survey (USGS) Mission—The USGS national mission is to monitor, analyze, and predict the current and evolving dynamics of complex human and natural Earth-system interactions and to deliver actionable information at scales and timeframes relevant to decision makers. Consistent with the national mission, the USGS in Alaska provides timely and objective scientific information to help address issues and inform management decisions across five interconnected focus areas: Energy and Minerals; Geospatial Mapping; Natural Hazards; Water Quality, Streamflow, and Ice Dynamics; and Ecosystems. The USGS in Alaska consists of approximately 350 scientists and support staff working in 3 Alaska-based science centers. USGS science activities are also initiated by the Cooperative Research Unit and USGS centers outside Alaska. In the last 5 years, USGS research in Alaska has produced many scientific benefits resulting from more than 900 publications. Publications relevant to Alaska can be conveniently searched by keyword through the USGS Publications Warehouse at https://pubs.usgs.gov/.

Alaska↗

USGS2025

Using subducting plate motion to constrain Cascadia slab geometry and interface strength

Subduction zones are home to multiple geohazards driven by the evolution of the regional tectonics, including earthquakes, volcanic eruptions and landslides. Past evolution builds the present-day structure of the margin, while the present-day configuration of the system determines the state-of-stress in which individual hazardous events manifest. Regional simulations of subduction zones provide a tool to synthesize the tectonic history of a region and investigate how geologic features lead to variations in the state of stress across the subduction system. However, it is challenging to design regional models that provide a force-balance that is consistent with the large-scale motion of surrounding tectonic plates while also not over-constraining the solution. Here, we present new models for the Cascadia subduction zone that meet these criteria and demonstrate how the motion of the subducting Juan de Fuca plate can be used to determine the along-strike variations in the viscous (long-term) coupling across the plate boundary. All successful models require lower viscous coupling in the northern section of the trench compared to the central and southern sections. However, due to uncertainties in the geometry of the Cascadia slab, we find that there is a trade-off between along-strike variation in viscous coupling and slab shape. Better constraints on the slab shape, and/or use of other observations are needed to resolve this trade-off. The approach presented here provides a framework for further exploring how geologic features in the overriding plate and the properties of the plate boundary region affect the state-of-stress across this and other subduction zones.

Cascadia subduction zone↗

USGS2025

Groundwater budget for the surficial aquifer surrounding Lake Nokomis, Minneapolis, Minnesota

During prolonged periods of above-average precipitation, rising groundwater levels have the potential to cause damage to and interfere with underground infrastructure and building foundations. To understand the relations between precipitation and groundwater in the vicinity of Lake Nokomis, the U.S. Geological Survey, in collaboration with the University of Minnesota, quantified five components of the groundwater budget: groundwater recharge, change in surficial aquifer storage, surficial aquifer groundwater discharge to Lake Nokomis, groundwater evapotranspiration, and groundwater discharge to underlying bedrock aquifers. Field data, geologic records, and empirical calculation methods were used to quantify groundwater budget components for April 2023 through April 2024. Lake water budget data indicate that Lake Nokomis is a flowthrough system during periods with no outflow through the weir, with groundwater inputs equal to outputs. Roughly 40 percent of precipitation that fell in the study area was added to the surficial aquifer as recharge. Uncertainty in the vertical hydraulic conductivity resulted in wide-ranging estimates (spanning three orders of magnitude) of water discharging from the surficial aquifer to the underlying bedrock aquifer. Drought conditions persisted for the duration of this study and were not representative of the conditions that motivated this study. This study is a start towards understanding relations between precipitation, Lake Nokomis levels, and groundwater levels that could affect local underground infrastructure.

Minnesota↗

USGS2024

StreamStats—A quarter century of delivering web-based geospatial and hydrologic information to the public, and lessons learned

StreamStats is a U.S. Geological Survey (USGS) web application that provides streamflow statistics, such as the 1-percent annual exceedance probability peak flow, the mean flow, and the 7-day, 10-year low flow, to the public through a map-based user interface. These statistics are used in many ways, such as in the design of roads, bridges, and other structures; in delineation of floodplains for land-use zoning and setting of insurance rates; for regulatory purposes, such as the permitting of wastewater discharges; and for hydrologic and climate change studies. StreamStats was first developed for Massachusetts and released in 2001. The application provided users with the ability to obtain streamflow statistics computed from data collected at USGS streamgages and to obtain estimates of streamflow statistics for user-selected ungaged sites. Massachusetts StreamStats used geographic information system software and digital mapping to compute drainage-basin characteristics, which were then used in statistical models to estimate streamflow statistics for the user-selected sites. The statistical models were in the form of equations that were developed through a process known as regression analysis. StreamStats was the first known web application with the ability to do interactive geoprocessing. The utility of Massachusetts StreamStats was instantly apparent, leading the USGS to develop a version of StreamStats that could be implemented nationally. USGS State offices normally were required to develop custom regression equations and prepare local digital mapping data needed for implementing StreamStats for their States. Funding needed to complete this work usually was provided through cooperative agreements between the USGS and State agencies. In 2004, Idaho became the first to be released in the national version of StreamStats. By 2023, 44 States were fully implemented and six were undergoing implementation. StreamStats has undergone many modifications over the years to keep up with changes to the underlying software and to add functionality. Customized functionality and separate linked StreamStats applications were developed for several States. Meeting the high demand for additions and improvements to StreamStats while also adhering to budgetary constraints has, at times, been challenging. The StreamStats development team has identified numerous additional improvements that could be made to provide better performance and more functionality. The lessons learned from the experience of building and operating StreamStats for nearly 25 years could be relevant to others interested in pursuing efforts of a similar scale.

Circular↗

USGS2026

Hydrogeologic conceptualization of the Arroyos subarea and Northeast Air Force Research Laboratory subarea of the Air Force Research Laboratory, Edwards Air Force Base, Antelope Valley, California, 2017–21

Research and testing within the Air Force Research Laboratory at Edwards Air Force Base in the southwestern Mojave Desert has been conducted since the 1950s. This testing requires large quantities of groundwater for cooling and operations involving hazardous constituents that have entered the groundwater. To better understand the potential movement of contaminants through the groundwater system, the U.S. Geological Survey entered into a cooperative agreement with the U.S. Air Force Civil Engineer Center with the purpose of (1) gaining a more detailed understanding of the hydraulic connection between the unconsolidated basin-fill sediment and the underlying weathered and crystalline bedrock and (2) investigating faults and their potential effect on groundwater flow. This study was specific to the Arroyos subarea and Northeast Air Force Research Laboratory subarea of the Air Force Research Laboratory. The refined hydrogeologic characterization was derived by evaluating surface geologic maps, a regional gravity model, new subsurface geophysical surveys, lithology data from boreholes, and groundwater-level data. Refinement of the hydrogeologic characterization of the study area involved developing a map of basin-fill sediment thickness, estimating the geometry and depth to bedrock, determining the presence of weathered bedrock material, and identifying the location of faults and subsurface structures. The thickness of basin-fill sediment was derived from a regional gravity model. Basin-fill sediment ranges from less than 20 feet (ft) thick to about 150 ft thick and increases substantially toward a deep structural basin in the northwestern part of the study area. The top of bedrock estimated from the new surface geophysical surveys is generally deeper than the top of bedrock from the regional gravity model. These differences in interpreted bedrock tops indicate that there may be a zone of weathered bedrock in the subsurface, overlying crystalline bedrock. Additional subsurface data evaluating the geometry and thickness of weathered bedrock are necessary for a comprehensive evaluation of weathered bedrock in the study area. Previously unmapped extensions of the Leuhman Fault, Spring Fault Zone, and Boron Gate Fault were identified in the new surface geophysical surveys. Evaluation of groundwater flow between bedrock and basin-fill sediment, and the effects of faults on groundwater flow, was done using groundwater-level data from 2020. The groundwater table was in crystalline bedrock throughout most of the Arroyos and Northeast Air Force Research Laboratory subareas, with a few exceptions where it was in basin-fill sediment and weathered bedrock. Basin-fill sediment, and likely weathered bedrock, was partially saturated in parts of the study area where basin-fill sediment was as thick as 150 ft. Weathered bedrock material may be present near the Spring Fault Zone, and if so, was likely saturated in 2020. The groundwater-level data near and within Spring Fault Zone indicate that the Spring Fault Zone acts as a vertical-along-strike conduit for, and as an across-strike impediment to, groundwater flow. In addition, steep groundwater-level gradients in the southeastern part of the study area indicate that a subsurface structure affecting groundwater flow could be related to the Leuhman Fault. Sufficient groundwater-level data were not available to conclusively determine if the Arroyos, Rich, and Boron Gate Faults, and their inferred extensions from the geophysical surveys, affect the movement or direction of groundwater flow.

California↗

USGS2025

Mechanics and statistics of postseismic shaking

Analysis of two weeks of continuous post-seismic shaking after the 2019 M7.1 Ridgecrest, CA earthquake sequence using 4 nearby borehole seismometers reveals that continuous ground motions decay as Omori’s law in time and follow the Gutenberg-Richter distribution in logarithmic amplitude. The measured temporal decay in amplitudes agrees with predictions of the rate-and-state framework and indicates shaking amplitudes are proportional to the velocity of afterslip. Our ground motion-based statistical framework provides a basis to forecast shaking intensity in the minutes to hours after a large earthquake.

California↗

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↗

USGS2000

Water quality of the Quaternary and Ada-Vamoosa aquifers on the Osage Reservation, Osage County, Oklahoma, 1997

The project was to provide information on the quality of ground water from rural-domestic-water wells within the Osage Reservation and compare the water-quality to proximity to oil wells. About 38,500 oil wells have been drilled in the Reservation since drilling began in 1896. About 1,480 square miles or 64 percent of the Reservation is within a quarter mile of an oil well. The unconfined Quaternary sand aquifer covers about 315 square miles or about 14 percent of the Reservation and the confined Ada-Vamoosa sandstone aquifer covers about 800 square miles or about 35 percent of the Reservation. Fifty-eight percent of the Quaternary aquifer and 69 percent of the outcrop area of the Ada-Vamoosa aquifer are within a quarter mile of an oil well. One hundred twenty domestic ground-water wells were sampled from the Quaternary and Ada-Vamoosa aquifers. Forty-nine percent of the Reservation is underlain by the aquifers. Ground-water quality is good on most of the Reservation, but the use of domestic-water-supply wells tend to minimize water-quality problems. Existing water-supply wells commonly are located in areas that produce usable volumes of potable water. Several constituents in samples from the Ada-Vamoosa aquifer within a quarter mile of an oil well were significantly greater than from the aquifer not near oil wells. The constituents include specific conductance, dissolved solids, sodium, sulfate, chloride, bromide, and silica. These ions are probably derived from brine water. In the Ada-Vamoosa aquifer subgroups, 57 percent of the samples near oil wells and 24 percent of the samples not near oil wells had dissolved-solids concentrations greater than 500 milligrams per liter. The water quality in the Quaternary and Ada-Vamoosa aquifers is similar in areas where no oil wells have been drilled but is significantly different for several constituents. Median concentrations of major constituents from the Ada-Vamoosa aquifer not near oil wells were less than or equal to values from the Quaternary aquifer. Sixty-four percent of the water-quality samples from the Quaternary and 51 percent from the Ada-Vamoosa aquifers have dissolved-solids concentrations less than the secondary drinking water regulations of 500 milligrams per liter. Fifty-nine percent of the aquifer samples in the Quaternary aquifer subgroups not near oil wells and 70 percent of the samples near oil wells had dissolved solids less than 500 milligrams per liter. Areas in the Ada-Vamoosa aquifer near Hominy, Pershing, and Hula Lake have dissolved-solids concentrations greater than the secondary drinking water regulations. Water-quality samples from the Quaternary aquifer in these areas also have dissolved-solids concentrations greater than 500 milligrams per liter.

Oklahoma↗

USGS2026

Earthquake ground-motion model adjustments for the San Francisco Bay area

We develop adjustments to ergodic ground‐motion models (GMMs) to improve their performance in the San Francisco Bay Area (SFBA). GMMs are widely used in hazard assessments to estimate characteristics of ground shaking based on known properties of the source, path, and site. Such models are often developed using datasets containing records from various regions, resulting in models that represent median ground‐motion behavior, which may not adequately represent ground motions within subregions. This is true for the SFBA, where ground motions attenuate more rapidly with distance than in many other parts of California that dominate GMM databases. To support improved seismic hazard estimates in the SFBA, we calculate regional constants and anelastic attenuation coefficient adjustments relative to two commonly used ergodic GMMs: BSSA14 ( Boore et al. , 2014 ) and ASK14 ( Abrahamson et al. , 2014 ). These adjustments are obtained for a suite of ground‐motion intensity measures (peak ground acceleration, peak ground velocity, and 5%‐damped pseudospectral acceleration at oscillator periods ranging from 0.075 to 10 s) using mixed‐effects regression. Use of the regionally adjusted models reduces the overall bias by up to 0.5 natural log units for BSSA14 and up to 0.6 natural log units for ASK14. We demonstrate one application of our attenuation adjustments and their implications in an earthquake early warning case study of the 2014 M 6.0 South Napa earthquake. The predicted extent of shaking using the adjusted models better matches observed shaking at large source‐to‐site distances, especially for lower shaking intensities, thus potentially reducing overalerting. We encourage the use of our model adjustments when ergodic models are considered for seismic hazard studies in the SFBA.

California↗

USGS2025

Migration of seismicity from the mantle to the upper crust beneath Harrat Lunayyir volcanic field, Saudi Arabia

Harrat Lunayyir is a volcanic field in Saudi Arabia that experienced a Mw~5.4 earthquake driven by an upper-crustal dike intrusion in May 2009. This volcanic field has exhibited numerous forms of volcanic seismicity both prior to and since the 2009 dike intrusion. Significantly, earthquakes within the lithospheric mantle and, rarely, the lower crust are present in the two-decade long seismicity catalog of Harrat Lunayyir. Here we analyze 24 years of volcanic seismicity at Harrat Lunayyir from 1998 to 2022. We find that: 1) precursory seismicity began at least eight years prior to the 2009 event, with a particularly notable seismic episode one year prior; 2) lithospheric mantle seismicity is highly localized in space and in time, largely occurring in discrete sequences lasting on the order of a few hours to a few days; 3) one seismic sequence clearly migrates upward from the lithospheric mantle to the upper crust, including seismicity within the nominally ductile lower crust; 4) crustal seismicity has been slowly declining over time; and 5) lithospheric-mantle seismicity does not show any apparent decline with time. From these observations we infer that the seismicity is driven by magmatic fluids or volatiles, and seismic monitoring of this volcanic field should continue into the future.

Harrat Lunayyir volcanic field↗

USGS2025

Landslide-channel feedbacks amplify channel widening during floods

Channel widening is a major hazard during floods, particularly in confined mountainous catchments. However, channel widening during floods is not well understood and not always explained by hydraulic variables alone. Floods in mountainous regions often coincide with landslides triggered by heavy rainfall, yet landslide-channel interactions during a flood event are not well known or documented. Here we demonstrate with an example from the Great Colorado Flood in 2013, a 1000 year precipitation event, how landslide-channel feedbacks can substantially amplify channel widening and flood risk. We use a combination of DEM differencing, field analysis, and multiphase flow modeling to document landslide-channel interaction during the flood event in which sediment delivered by landslides temporarily dammed the channel before failing and generating substantial channel widening. We propose that such landslide-flood interactions will become increasingly important to account for in flood hazard assessment as flooding and landsliding both increase with extreme rainfall under climate change.

Colorado↗