USGS Science⌕ Search

SEARCH · USGS Science

Results for “Science of the Total Environment”

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

Multi-scale predictors of Northern Long-eared Bat (Myotis septentrionalis) occupancy in the United States

Historically, Myotis septentrionalis (Northern Long eared Bat) was among the most common forest-interior species in North America. Largely due to high mortality from white-nose syndrome, this species has experienced severe population declines across its range. To create an updated species distribution map representing summer occupancy probabilities from 2017 to 2022, we integrated stationary acoustic data with live-capture data from the database of the North American Bat Monitoring Program into a multi-scale, multi-method occupancy modeling framework. Our results provide data-driven predictions with quantified uncertainty for summer occupancy probabilities for Northern Long-eared Bats at 2 spatial scales across the range of the species, while also accounting for inherent observation biases (e.g., imperfect detection).

Journal of North American Bat Research↗
USGS2024

Layered intrusions in the Precambrian: Observations and perspectives

Layered intrusions are plutonic bodies of cumulates that form by the crystallization of mantle-derived melts. These intrusions are characterized by igneous layering distinguishable by shifts in mineralogy, texture, or composition. Layered intrusions have been fundamental to our understanding of igneous petrology; however, it is their status as important repositories of critical metals – such as platinum-group elements, chromium, and vanadium – that has predominantly driven associated research in recent decades. Many layered intrusions were emplaced during the Precambrian, predominantly at the margins of ancient cratons during intervals of supercontinent accretion and destruction. It appears that large, layered intrusions require rigid crust to ensure their preservation, and their geometry and layering is primarily controlled by the nature of melt emplacement. Layered intrusions are best investigated by integrating observations from various length-scales. At the macroscale, intrusion geometries can be discerned, and their presence understood in the context of the regional geology. At the mesoscale, the layering of an intrusion may be characterized, intrusion-host rock contact relationships studied, and the nature of stratiform mineral occurrences described. At the microscale, the mineralogy and texture of cumulate rocks and any mineralization are elucidated, particularly when novel microtextural and mineral chemical datasets are integrated. For example, here we demonstrate how mesoscale observations and microscale datasets can be combined to understand the petrogenesis of the perplexing snowball oiks outcrop located in the Upper Banded Series of the Stillwater Complex. Our data suggest that the orthopyroxene oikocrysts did not form in their present location, but rather formed in a dynamic magma chamber where crystals were transported either by convective currents or within crystal-rich slurries. Critical metals may be transported to the level of a nascent intrusion as dissolved components in the melt. Alternatively, ore minerals are entrained from elsewhere in a plumbing system, potentially facilitated by volatile-rich phases. There are many ore-forming processes propounded by researchers to occur at the level of emplacement; however, each must address the arrival of the ore mineral, its concentration of metals, and its accumulation into orebodies. In this contribution, several of these processes are described as well as our perspectives on the future of layered intrusion research.

Precambrian Research↗
USGS2026

Baseflow and snowmelt sustained streamflow in the Upper Colorado River Basin, 1986-2020

The Upper Colorado River Basin (UCRB) faces substantial water availability limitations. Although most streamflow originates as snowmelt, the partitioning of snowmelt between surface runoff and groundwater recharge and subsequent groundwater discharge to streams is highly uncertain. On average, over half of the streamflow in the UCRB is estimated to originate from groundwater discharge to streams, highlighting the importance of baseflow in sustaining surface water. However, the historical patterns of baseflow and streamflow, along with their variability over space and time and their specific sources, remain unknown at the basin scale. This study addresses those gaps by characterizing the sources and transport pathways of both baseflow and streamflow in the UCRB at a seasonal timestep from 1986 to 2020, including the lagged delivery of subsurface water to streams beyond the current season, using coupled models of baseflow and streamflow. Between 1986 and 2020, on average 63% of UCRB streamflow originated from baseflow. About half of this baseflow took longer than one season to reach streams, and outside the snowmelt season, baseflow was the dominant source of streamflow. Snowmelt was a key source of both baseflow and streamflow. Current season snowmelt contributed 33% of streamflow via runoff, and 22% of the 29% of streamflow that originated as current season baseflow via subsurface flow to streams. Over the study period, baseflow index (BFI) declined in headwaters and increased at mid-elevations. Springtime increases in BFI demonstrate the increasingly important role baseflow plays in water supply. Identifying the sources, locations, and timing of water that contributed to the UCRB outlet can inform management of water resources in the basin.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2024

Testing rate‐and‐state predictions of aftershock decay with distance

We analyze aftershocks of the 2019 M 7.1 Ridgecrest mainshock and isolated M 5–6 mainshocks in southern California to test predictions made by the rate‐and‐state friction model of Dieterich (1994) . Rate‐and‐state friction predicts that the seismicity rate after a stress step follows Omori decay, where the Omori c ‐value, which is the saturation in aftershock rate observed at small times, is larger for smaller stress steps. Put in the context of an aftershock sequence, this predicts that the Omori c ‐value will be systematically larger at greater distances from the mainshock. To our knowledge, this predicted effect has not been observed. In part this may be because the Omori c ‐value is difficult to measure because it often reflects short‐term catalog incompleteness rather than a true saturation in aftershock rate. We explore the dependence of the Omori c ‐value on the distance to the mainshock by applying the “ a ‐positive” method ( van der Elst and Page, 2023 ). This method is insensitive to short‐term aftershock incompleteness and allows resolution of the true aftershock rate deep into the mainshock coda. For aftershocks of the Ridgecrest mainshock and stacked M 5–6 mainshocks, we observe systematic differences in early aftershock rates, relative to mainshock distance, consistent with the predictions of rate‐and‐state friction. Furthermore, for the larger Ridgecrest dataset, we observe that aftershocks nearer to the mainshock start earlier, and we resolve a flattening of the Omori curve consistent with a larger Omori c ‐value for the farthest aftershocks, as predicted by Dieterich (1994) .

Seismological Research Letters↗
USGS1998

Assessment of the hydraulic connection between ground water and the Peace River, west-central Florida

The hydraulic connection between the Peace River and the underlying aquifers along the length of the Peace River from Bartow to Arcadia was assessed to evaluate flow exchanges between these hydrologic systems. Methods included an evaluation of hydrologic and geologic records and seismic-reflection profiles, seepage investigations, and thermal infrared imagery interpretation. Along the upper Peace River, a progressive long-term decline in streamflow has occurred since 1931 due to a lowering of the potentiometric surface of the Upper Floridan aquifer by as much as 60 feet because of intensive ground-water withdrawals for phosphate mining and agriculture. Another effect from lowering the potentiometric surface has been the cessation of flow at several springs located near and within the Peace River channel, including Kissengen Spring, that once averaged a flow of about 19 million gallons a day. The lowering of ground-water head resulted in flow reversals at locations where streamflow enters sinkholes along the streambed and floodplain. Hydrogeologic conditions along the Peace River vary from Bartow to Arcadia. Three distinctive hydrogeologic areas along the Peace River were delineated: (1) the upper Peace River near Bartow, where ground-water recharge occurs; (2) the middle Peace River near Bowling Green, where reversals of hydraulic gradients occur; and (3) the lower Peace River near Arcadia, where ground-water discharge occurs. Seismic-reflection data were used to identify geologic features that could serve as potential conduits for surface-water and ground-water exchange. Depending on the hydrologic regime, this exchange could be recharge of surface water into the aquifer system or discharge of ground water into the stream channel. Geologic features that would provide pathways for water movement were identified in the seismic record; they varied from buried irregular surfaces to large-scale subsidence flexures and vertical fractures or enlarged solution conduits. Generally, the upper Peace River is characterized by a shallow, buried irregular top of rock, numerous observed sinkholes, and subsidence depressions. The downward head gradient provides potential for the Peace River to lose water to the ground-water system. Along the middle Peace River area, head gradients alternate between downward and upward, creating both recharging and discharging ground-water conditions. Seismic records show that buried, laterally continuous reflectors in the lower Peace River pinch out in the middle Peace River streambed. Small springs have been observed along the streambed where these units pinch out. This area corresponds to the region where highest ground-water seepage volumes were measured during this study. Further south, along the lower Peace River, upward head gradients provide conditions for ground-water discharge into the Peace River. Generally, confinement between the surficial aquifer and the confined ground-water systems in this area is better than to the north. However, localized avenues for surface-water and ground-water interactions may exist along discontinuities observed in seismic reflectors associated with large-scale flexures or subsidence features. Ground-water seepage gains or losses along the Peace River were quantified by making three seepage runs during periods of: (1) low base flow, (2) high base flow, and (3) high flow. Low and high base-flow seepage runs were performed along a 74-mile length of the Peace River, between Bartow and Nocatee. Maximum losses of 17.3 cubic feet per second (11.2 million gallons per day) were measured along a 3.2-mile reach of the upper Peace River. The high-flow seepage run was conducted to quantify losses in the Peace River channel and floodplain between Bartow and Fort Meade. Seepage losses calculated during high-flow along a 7.2-mile reach of the Peace River, from the Clear Springs Mine bridge to the Mobil Mine bridge, were approximately 10 percent of the river flow, or 118 cubic feet per second. Calculated seepages along the Peace River in Hardee and De Soto Counties were inconclusive, because most seepages were within the range of discharge measurement error. Two continuous aerial thermal infrared imagery surveys were conducted to locate sites of ground-water discharge along the Peace River. Although temperature and hydrologic conditions were ideal to observe spring flow using thermal infrared imaging techniques, no sources of ground-water discharge were identified using this method. Diffuse ground-water seepage may, however, provide significant ground-water discharge.

Florida↗
USGS2026

Influence of new residential construction varying in housing density on bird species, human tolerance guilds, and communities

Human population growth and changing settlement patterns fuel the development of urban fringe lands worldwide, with implications for biodiversity. We conducted a 12-year study of birds in the fast-developing urban fringe lands of the central Puget Sound region, Washington, USA, to examine the effect of development configuration on birds. We hypothesized that lower-intensity conservation developments, compared to higher-intensity planned community developments, would benefit the overall bird community, as well as native forest birds (avoiders of human development) and avian generalist species (adapters to human development), but that higher-intensity planned community developments would benefit synanthropic species (exploiters of human development). We fit single-species and multi-species occupancy models to test these hypotheses. Consistent with our hypotheses, we found that a greater proportion of the overall community, avoiders, and adapters occupied lower-intensity conservation developments compared to higher-intensity planned community developments. However, we did not detect an effect of development type on the exploiter guild, and we found that species in the exploiter guild are variable in their response to the configuration of suburban developments. We also hypothesized that human tolerance guilds would be a useful predictor of individual species responses to development type. This hypothesis was somewhat supported: we found that, for avoiders, 87% of species in the guild had the same response to development type as the overall guild; for adapters, 63% had the same response as the overall guild, and for exploiters, only 44% had the same response as the overall guild. Our results indicate that the configuration of suburban developments can have a meaningful impact on bird communities, particularly on those species that are most sensitive to any level of development.

Washington↗
USGS2025

Footprints of past mining in Alaska (USA) derived from high-resolution satellite imagery

Mapping the land area used for mining in the past is essential for guiding the remediation of affected landscapes and assessing the resource potential of related waste products. Despite significant recent progress delineating footprints of active and inactive mining globally, the known inventory of such mine lands remains incomplete. Here, I describe a new map dataset of footprints of land surface disturbance and waste at sites of past mining in Alaska (USA) based on visual interpretation of satellite imagery. This dataset maps 6–14 times the area of previous regional and global mine footprint maps in Alaska and is the first in the region to explicitly delineate mine waste landforms (e.g., tailings piles). The data are publicly available from the U.S. Geological Survey under a “no rights reserved” Creative Commons (CC0) license agreement.

Alaska↗
USGS2026

Temporal and spatial changes in seismic attenuation associated with inferred fluid migration in the 2016 central Apennines earthquake sequence

Prior work suggests that high‐frequency seismic attenuation acts as a highly sensitive proxy for crustal permeability and fluid mobility in fractured media. We test the hypothesis that the fault system responsible for the 2016–2017 Amatrice–Visso–Norcia–Capitignano sequence acted as an impermeable seal, compartmentalizing pressurized fluids until dynamic rupture triggered widespread fluid diffusion. By tracking across the sequence the spatiotemporal evolution of the S ‐wave anelastic attenuation parameter, we identify large, positive low‐frequency attenuation anomalies emerging within the hanging wall following the Amatrice mainshock and strictly preceding subsequent large ruptures. Conversely, we observe weaker, negative anomalies in the footwall, anticorrelated in time with those of the hanging wall, revealing a massive asymmetry in fluid redistribution and permeability evolution across the fault system. Furthermore, aftershock migration rates reveal distinct linear alignments in a distance‐reduced time space, allowing us to explicitly track and quantify episodes of lateral and upward fluid migration. These physically consistent patterns suggest that stress‐driven fluid diffusion directly weakens adjacent fault patches, dictating the spatiotemporal migration of seismicity. We conclude that near‐real‐time monitoring of seismic attenuation may help detect fluid redistribution in active fault systems and may provide useful information for time‐dependent seismic hazard assessment.

central Apennines↗
USGS2026

Lithium brine potential of North Dakota: Results from the 2025 sampling program

Lithium (Li), classified as a critical mineral and key in energy storage applications because of its high energy density, faces unprecedented demand. This has driven interest in unconventional sources of Li, especially brines from sedimentary basins. Accordingly, the spatial distribution of basinal brinehosted Li resources across North America is becoming better defined, including the Williston Basin underlying Saskatchewan, Canada. This has led to growing investment and exploration by industry throughout the 2020s; however, despite this interest, the Li potential of basinal brines of the Williston Basin in North Dakota remain underexplored. This work presents the results from the 2025 joint brine sampling campaign between the North Dakota Geological Survey, the United States Geological Survey (USGS), and The University of Texas at Austin. Twenty-seven fluid samples were collected from producing oil and water source wells from eight stratigraphic intervals across western North Dakota and analyzed for major cation and anion concentrations. Similar to southeastern Saskatchewan, the highest Li concentrations were identified in the Devonian Duperow Formation, with concentrations up to 179 mg/L being observed in a well commingled between the Duperow and Red River Formations. Some samples from the overlying Birdbear Formation exceeded 75 mg/L, whereas Li concentrations were lower in the Madison Group (33-44 mg/L) and the Three Forks (50 mg/L), Dawson Bay (43 mg/L), Winnipegosis (23 mg/L), and Red River (40.0 mg/L) Formations. This work indicates Li resource potential exists across the Williston Basin and advances our understanding of the Li distribution in basinal brines which could provide a new domestic source of this critical mineral.

North Dakota↗
USGS2025

Near-fault amplification and ground motion variability during the 2019 Ridgecrest, California sequence

We estimate ground-motion variability near the 2019 M 7.1 Ridgecrest earthquake sequence. Accurate seismic hazard estimation requires understanding ground-motion spatial correlations, yet many studies lack the dense station coverage needed to resolve small-scale variability. The 2019 M 7.1 Ridgecrest earthquake sequence presents a unique opportunity to examine ground motions and their spatial correlations at a range of interstation distances. The permanent seismic network was augmented with hundreds of temporary stations including several fault-crossing nodal arrays. We compute the event ( δE i ) and within-event ( δW ij ) residuals from the observed peak ground velocity and peak ground acceleration data to isolate potential sources of ground-motion variability. We then compare δW ij between station pairs that record an event to understand the semivariance of the ground motion versus interstation distance. By fitting an exponential model to the semivariances, we determine a correlation range of 25 km for the Ridgecrest region. Although the exponential model fits the broad-scale increase of semivariance with interstation distance, we also observe smaller-scale trends. We find that ground motions are less correlated for station pairs that are near or across faults that ruptured during the 2019 Ridgecrest sequence. We also find large, positive median δW ij with relative values 2–3 times larger than nearby stations for individual stations’ near-fault traces. Near-fault amplification and greater ground-motion variability can delineate fault zones and may locally increase the seismic hazard.

California↗
USGS2025

Flood-Inundation Maps of the Current and Jacks Fork Rivers including the Ozark National Scenic Riverways, Southeast Missouri, 2023

Digital flood-inundation maps for a 131.8-mile reach of the Current River and a 44.6-mile reach of the Jacks Fork River, in southeast Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the Ozark Foothills Regional Planning Commission and the South Central Ozark Council of Governments. The maps also encompass the 134 miles of the Current and Jacks Fork Rivers within the Ozark National Scenic Riverways, which is the first national park area to protect a river system. The flood-inundation maps, which can be accessed through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding corresponding to selected water levels (hereafter referred to as “stages”) at eight reference USGS streamgages—five on the Current River (USGS station numbers 07064440, 07064533, 07066510, 07067000, and 07068000) and three on the Jacks Fork River (USGS station numbers 07065200, 07065495, and 07066000). Near-real-time stages at these streamgages may be obtained from the USGS National Water Information System at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at four of these sites (USGS station numbers 07067000, 07068000, 07065495, and 07066000). Flood profiles were computed for seven of the eight map reaches by means of two-dimensional hydraulic models and the remaining reach by a one-dimensional hydraulic model. The models were calibrated by using stage-streamflow relations or streamflow measurements at the USGS streamgages and from high-flow stage measurements from water-level loggers distributed throughout the reaches. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “action stage” or near bankfull, to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root mean square error) to delineate the area flooded at each water level and the associated water depths. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that is critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Missouri↗
USGS2025

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

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

Mu Us Sandy Land↗
USGS2026

Surface rupture and slip distribution of the 2025 Mw7.7 Mandalay earthquake and updated length scaling of supershear earthquakes

The 2025 M w 7.7 Mandalay, Burma (Myanmar), earthquake ruptured 475 km of the central Sagaing fault and is the longest continental strike-slip rupture on record. The observed rupture length is 1.6–4.7 times the value expected (100–300 km) from existing length-magnitude scaling relations for strike-slip earthquakes. The earthquake resulted from shallow dextral faulting and ruptured bilaterally with supershear speeds south of the epicenter, rupturing close to three major cities in Myanmar and exposing over six million people to violent or extreme shaking. We report on the surface rupture character, length, and slip distribution based on sub-pixel correlation of Sentinel-2 (10 m) and Planet Dove (3 m) optical images and visual analysis of SkySat and WorldView (0.3–0.5 m) optical images. The earthquake had moderate surface slip (average = 3.3 m, maximum = 5.6 m, 25–75% range = 3.0–4.0 m), narrow deformation zone width (1–10 pixels in sub-pixel correlation and up to 190 meters for the detailed surface rupture mapping), and simple fault geometry (no stepovers or large changes in strike, 87% of the rupture that was mapped in detail is single-stranded). We attribute the extreme length of the Mandalay earthquake to supershear rupture speed, simple fault geometry, narrow down-dip width, and moderate surface slip. Based on a compilation of 25 supershear strike-slip earthquakes ( M w 6.5–8.6; 1979–2025), we find that the rupture length of supershear earthquakes does not fit empirical scaling relationships for strike-slip earthquakes that predict length from magnitude. A length-magnitude scaling relationship based on supershear earthquakes has a best fit of log 10 (surface rupture length) = 0.89 M w – 4.44, indicating that supershear earthquakes tend to be longer than their subshear counterparts for any given magnitude and thus may expose a greater population to shaking.

Mandalay↗
USGS2024

Evidence of nitrate attenuation in intertidal and subtidal groundwater in a subterranean estuary at a Cape Cod embayment, East Falmouth, Massachusetts, 2015–16

Nitrogen dynamics in intertidal and nearshore subtidal groundwater (subterranean estuary) adjacent to the Seacoast Shores peninsula, Falmouth, Massachusetts, were investigated during 2015–16 by the U.S. Geological Survey. The peninsula is a densely populated residential area with septic systems and cesspools that are substantial sources of nitrogen to groundwater. The study area is in the Eel River, an estuarine saltwater embayment connected to the ocean adjacent to the western shore of the peninsula, that was the subject of an earlier study by Colman and others (2018, https://doi.org/10.3133/sir20185095 ) on nitrogen transport and transformations in groundwater between onshore and offshore locations. The previous study documented the distribution of nitrate concentrations and nitrate attenuation reactions in fresh groundwater beneath the peninsula and the estuary. The current study extended those observations with more detailed sampling and analysis of shallow groundwater from wells near discharge sites beneath the estuary. The current field investigation included sampling of existing wells and installation and sampling of clusters of wells and temporary sampling points in the subterranean estuary, including (1) shallow transects 0.3 to 1.2 meters (m) deep extending from 1 to 13.5 m offshore and (2) deeper wells (from 1.83 to 4.88 m deep) extending from 4.3 to 14.3 m offshore. Measurements of hydraulic-head gradients 2–5 m below the sediment/water interface in the intertidal and nearshore subtidal zones indicated that groundwater flow generally was upwards (towards the estuary) under all tide conditions in October 2016. The magnitude of the gradient was greatest during low tide conditions, indicating that groundwater discharge likely decreased during high tides. Measurements of specific conductance in shallow groundwater in the subterranean estuary in three transects perpendicular to shore were consistent with the existence of saltwater flow cells (infiltration of overlying saline water, mixing with fresh groundwater, and discharge to the overlying saline water) in the intertidal and nearshore subtidal regions. The size of these flow cells was variable in space and time and dependent on the elevation of the tide (spring or neap). At this location in the Eel River subterranean estuary, and offshore to at least 13.5 m, offshore flow of fresh groundwater apparently prevented a deeper saltwater wedge from discharging to the surface. Nitrate concentrations in shallow groundwater (30 to 122 centimeters [cm] depth) were variable in space and time, ranging from not detectable to 600 micromoles per liter (μmol/L) (8.4 milligrams per liter as N), and were highest in June 2016 at depths from 61 to 122 cm below the sediment/water interface and from 4 to 9 m offshore. Nitrate generally was not detectable in saline shallow groundwater at 30-cm depth or at any depth from 30 to 122 cm from 10 to 13.5 m offshore. Dissolved oxygen concentrations were suboxic (less than 16 μmol/L) in 60 percent of the sampled subterranean groundwater beneath the intertidal and subtidal zones. In the remaining sites, the range of dissolved oxygen concentrations was from 18 to 272 μmol/L and the median concentration was 43 μmol/L. Evidence for microbial nitrate reduction (denitrification and possibly anammox) was provided by the distribution of the reaction product nitrogen gas (excess N 2 , or N 2MIC ), as determined from analysis of the dissolved nitrogen gas and argon gas (Ar) concentrations in groundwater samples. Excess nitrogen gas provided evidence for nitrate reduction in shallow groundwater below the subtidal and, to a lesser extent, intertidal zones adjacent to the Seacoast Shores peninsula. These zones, where evidence for nitrate reduction was detected, were in fresh and brackish groundwater near subtidal or intertidal saltwater cells where discharging fresh groundwater mixed with infiltrating saline water. Infiltrating seawater may have supplied organic carbon, one of several potential electron donors that are required for denitrification. Other potential electron donors, such as organic carbon, iron, manganese, hydrogen, methane, ammonium, elemental sulfur, or sulfide phases, may have been supplied by the estuarine sediments. Drainage from surface runoff near the shore also may have supplied organic carbon to fresh groundwater near the intertidal saltwater cell. The highest amounts of nitrate converted to excess nitrogen gas were estimated to be in the range of 230 to 430 μmol/L in nearly fresh groundwater near the subtidal saltwater cell at depths of 61 to 122 cm below the sediment/water interface and from 10 to 13.5 m offshore. Evidence of denitrification within 10 m of the shore was sparse (generally limited to less than 50 μmol/L of N 2 -N) despite the presence of high nitrate concentrations. The spatial distribution of estimated nitrate reduction in the intertidal and nearshore subtidal fresh and brackish groundwater may be related to local variability in the distribution of reactive electron donors in those zones. Variations in the amount of nitrate reduction to nitrogen gas were not clearly related to potential aqueous electron donors such as dissolved organic carbon, nor to potential reaction products such as alkalinity, but may have been controlled by combinations of aqueous and solid-phase reactants. The distribution of relatively shallow fresh groundwater containing nitrate could indicate potential nitrate discharge areas in the lower intertidal zone and uncertain locations farther offshore; however, the data did not extend all the way to the sediment/water interface or to the offshore freshwater limit. This study confirmed substantial loss of nitrate from some of the fresh and brackish groundwater in shallow subestuarine sediments prior to discharge but did not quantify how much nitrate eventually discharged to the estuary.

Massachusetts↗
USGS2026

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

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

International Journal of Disaster Risk Reduction↗
USGS2025

False positives in the identification of dynamic earthquake triggering

Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with global M ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial-temporal windows in Southern California tested in DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 , thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023, https://doi.org/10.1029/2023jb026487 ), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor-Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.

JGR Solid Earth↗
USGS2024

Limited preservation of strike-slip surface displacement in the geomorphic record

Offset geomorphic markers are commonly used to interpret slip history of strike-slip faults and have played an important role in forming earthquake recurrence models. These data sets are typically analyzed using cumulative probability methods to interpret average amounts of slip in past earthquakes. However, interpretation of the geomorphic record to infer surface slip history is complicated by slip variability, measurement uncertainty, and modification of offset features in the landscape. To investigate how well geomorphic data record surface slip, we use offset measurements from recent strike-slip surface ruptures ( n = 39), faults with geomorphic evidence of multiple strike-slip earthquakes ( n = 29), and synthetic slip distributions with added noise ( n > 10,000) to examine the constraints of the geomorphic record and the underlying assumptions of the cumulative offset probability distribution analysis method. We find that the geomorphic record is unlikely to resolve more than two paleo-slip distributions, except in specific cases with low slip variability, high slip-per-event, and semiarid climate. In cases where site-specific conditions allow for interpretation of more than two earthquakes, lateral extrapolation along a fault is not straightforward because on-fault displacement and distributed deformation may be spatially variable in each earthquake. We also find that average slip in modern earthquakes is adequately recovered by probability methods, but the reported prevalence of strike-slip faults with characteristic slip history is not supported by geomorphic data. We also propose updated methods to interpret slip history and construct uncertainty bounds for paleo-slip distributions.

Journal of Geophysical Research: Solid Earth↗
USGS2025

Evaluating Laramide orogenesis via flexural basin response in the San Juan basin, New Mexico and Colorado

A challenge in interpreting the location, timing, and magnitude of ancient orogenic events is that ongoing uplift and erosion in the hinterlands often destroys much of the primary record of these events. However, basin-thickness patterns in the sedimentary record can provide complimentary evidence of uplift via flexural effects. Here, we deploy well-log correlation, isochores, basin modeling, flexural modeling, and subcrop mapping to evaluate the Late Cretaceous to Paleogene basin response to Laramide tectonism in the San Juan basin. A wedge of upper Campanian to Maastrichtian sedimentary rock thickens from 200 to 800 meters from southeast to northwest in the basin. This pattern can be successfully simulated via flexural modeling if we infer early Laramide uplift along the northwest basin flank that produced a 0.8 km high topographic load. The Laramide unconformity bounds the top of this Upper Cretaceous sedimentary wedge and truncates progressively older strata to the east, further supporting a westward tilt of the basin. The onset of Campanian Laramide flexure may have also contributed to the profound transgression from the upper Menefee Formation to the Lewis Shale. The Paleocene isochore map displays an approximately symmetrical pattern, with thickening towards the center of the basin. This suggests the possibility of competing flexural loads. The base Eocene structure indicates an asymmetric deep on the northeast flank of the basin, providing flexural evidence of contemporaneous uplift/loading of the Nacimiento uplift and Archuleta arch; this has been modeled as ~2.1 km load height. Both Cretaceous and Paleocene sedimentary wedges are narrow, suggesting low flexural rigidity; modeled effective elastic thicknesses (EET) are 20-30 km, comparable to estimates of modern EET for the region.

Colorado, New Mexico↗