USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Science: Processes and Impacts”

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,662 records · Page 67Linked to original sources
USGS2024

Static and dynamic strain in the 1886 Charleston, South Carolina, earthquake

During the 1886 Mw 7.3 Charleston, South Carolina, earthquake, three railroads emanating from the city were exposed to severe shaking. Expansion joints in segmented railroad tracks are designed to allow railroad infrastructure to withstand a few parts in 10,000 of thermoelastic strain. We show that, in 1886, transient contractions exceeding this limiting value buckled rails, and transient extensions pulled rails apart. Calculated values for dynamic strain in the meizoseismal region are in reasonable agreement with those anticipated from the relation between strain and moment magnitude proposed by Barbour et al. (2021) and exceed estimated tectonic strain released by the earthquake by an order of magnitude. Almost all of the documented disturbances of railroad lines, including evidence for shortening of the rails, can thus be ascribed to the effects of dynamic strain changes, not static strain. Little or no damage to railroads was reported outside the estimated 10 −4 dynamic strain contour. The correspondence between 10 −3 and 2×10 −4 contours of dynamic strain and Mercalli intensity 9 and 8, anticipated from the dependence of each quantity on peak ground velocity, suggests it may be possible to use railroad damage to quantitatively estimate shaking intensity. At one location, near Rantowles, ≈20 km west of Charleston, a photograph of buckled track taken one day after the earthquake has been cited as evidence for shallow dextral slip and has long focused a search for a causal fault in this region. Photogrammetric analysis reveals that the buckle was caused by transient contraction of <10 cm with no dextral offset. Our results further weaken the evidence for faulting in the swamps and forests south of the Ashley River in 1886, hitherto motivated by the photograph and limited macroseismic evidence for high‐intensity shaking.

South Carolina↗
USGS2024

Improving crop-specific groundwater use estimation in the Mississippi Alluvial Plain: Implications for integrated remote sensing and machine learning approaches in data-scarce regions

Study region The Mississippi Alluvial Plain (MAP) in the United States (US). Study focus Understanding local-scale groundwater use, a critical component of the water budget, is necessary for implementing sustainable water management practices. The MAP is one of the most productive agricultural regions in the US and extracts more than 11 km 3 /year for irrigation activities. Consequently, groundwater-level declines in the MAP region pose a substantial challenge to water sustainability, and hence, we need reliable groundwater pumping monitoring solutions to manage this resource appropriately. New hydrological insights for the region We incorporate remote sensing datasets and machine learning to improve an existing lookup table-based model of groundwater use previously developed by the U.S. Geological Survey (USGS). Here, we employ Distributed Random Forests, an ensemble machine learning algorithm to predict annual and monthly groundwater use (2014–2020) throughout this region at 1-km resolution, using pumping data from existing flowmeters in the Mississippi Delta. Our model compares favorably with the existing USGS model, with higher R 2 (0.51 compared to 0.42 in the previous model), and lower root mean square error (RMSE) and mean absolute error (MAE)— 0.14 m and 0.09 m, respectively in our model, compared to 0.15 m and 0.1 m in the previous model. Therefore, this work advances our ability to predict groundwater use in regions with scarce or limited in-situ groundwater withdrawal data availability.

Journal of Hydrology Regional Studies↗
USGS1990

Effects of storm-water runoff on local ground-water quality, Clarksville, Tennessee

Storm-related water-quality data were collected at a drainage-well site and at a spring site in Clarksville, Tennessee, to define the effects of storm-water runoff on the quality of ground water in the area. A dye-trace test verified the direct hydraulic connection between the drainage well and Mobley Spring. Samples of storm run off and spring flow were collected at these sites for nine storms during the period February to October 1988. Water samples were collected also from Mobley Spring and two other springs and two observation wells in the area during dry-weather conditions to assess the general quality of ground water in an urban karst terrain. Evaluation of the effect of storm-water runoff on the quality of local ground water is complicated by the presence of other sources of contaminants in the area Concentrations and load for most major constituents were much smaller in storm-water runoff at the drainage well than in the discharge of Mobley Spring, indicating that much of the chemical constituent load discharged from the spring comes from sources other than the drainage well. However, for some of the minor constituents associated with roadway runoff (arsenic, copper, lead, organic carbon, and oil and grease), the drainage well contributed relatively large amounts of these constituents to local ground water during storms. The close correlation between concentrations of total organic carbon and concentrations of most trace metals at the drainage-well and Mobley Spring sites indicates that these constituents are transported together. Many trace metals were flushed early during each runoff event. Mean storm loads for copper, lead, zinc, and four nutrient species (total nitrogen, ammonia nitrogen, total phosphorus, and orthophosphorus) in storm-water runoff at the drainage-well site were lower than mean storm load predicted from an existing regression model. The overprediction by the model may be a result of the small size of the drainage area relative to the range of drainage areas used in the development of the models, or to the below-normal amounts of rainfall during the period of sampling for this investigation. Loads& in storm-water runoff for 22 constituents were extrapolated from sampled storms to total loads for the period February to October 1988. Calculated loads for trace metals for the period ranged from 0.030pound.s for cadmium to 12pound.s for strontium. Loads of the primary nutrients ranged from 0.97pounds for nitrite as nitrogen to 34pounds of organic nitrogen. Storm-water quality at the drainage-well and Mobley Spring sites was compared to background water quality of the local aquifer; as characterized by dry-weather samples from three springs and two observation wells in the Clarksville area. Concentrations of total-recoverable cadmium, chromium, copper, lead, and nickel were higher in many stormwater samples from both the drainage-well and Mobley Spring sites than in samples from any other site. In addition, concentrations of total organic carbon, methylene blue active substances, and total-recoverable oil and grease were generally higher in storm-water samples from the drainage-well site than in any ground-water sample. Densities of fecal coliform and fecal streptococcus bacteria and concentrations of total recoverable iron, manganese, and methylene blue active substances in storm samples from the drainage-well site exceeded the maximum contaminant levels listed in Tennessee’s drinking-water standards (1988) by as much as 2,500 and 5,500 colonies per 100 milliliters, and 2.7, 0.29, and 0.05 milligrams per liter, respectively. Densities of fecal coliform and fecal streptococcus bacteria and concentrations of total-recoverable iron, manganese, and lead in storm samples from Mobley Spring exceeded the maximum contaminant levels by as much as 500 and 4,500 colonies per 100 milliliters, and 18.7,0.65, and 0.02 milligrams per liter, respectively. For iron, manganese, and bacteria, these undesirable levels are not necessarily attributable to storm-water recharge, because concentrations of these constituents also exceeded drinking-water standards in one or more of the dry-weather samples from selected springs and observation wells in the area.

Tennessee↗
USGS2025

Local, regional, and distal recordings of seismic unrest at Tau Island volcano, American Samoa

A seismic swarm near Taʻū Island, a volcanic island in eastern American Samoa, occurred from July to October 2022. The earliest unrest was noted as felt shaking reports in late July, and instrumentation varied in the beginning of the sequence as the U.S. Geological Survey (USGS) Hawaiian Volcano Observatory responded by installing temporary and then permanent seismometers to monitor the activity. This network variability made it difficult to characterize the earliest seismicity and contextualize the entire sequence to discriminate between an underlying tectonic or volcanic source. Here, we present results analyzing hydroacoustic detections from an International Monitoring System hydrophone array near Wake Island, 4500 km northwest of Taʻū Island volcano. Using least-squares beamforming analysis, we create a catalog of T-wave detections from the direction of Taʻū Island to track the earthquakes, some of which were located by the USGS National EarthquakeInformation Center. Both the rate and hydroacoustic pressures, which we interpret as a proxy for earthquake size, gradually increased from late July to August, peaking on August 19 (rate) and August 24 (size), before decreasing to background in late September. Minutes-long bursts of tremor were also contemporaneously recorded as local network data became avail-able on August 20. Tremor activity continued throughout the rest of August, peaking on August 25, before ending in earlySeptember. These tremor bursts were band-limited to ~ 1–5 Hz and recorded as S waves at a regional station on the island of Upolu in Samoa, 250 km to the west of Taʻū Island. Our results do not constrain the tremor locations, but comparisons of earthquake and tremor reduced displacements recorded locally and regionally suggest a deeper tremor source. We interpret the increase in earthquake size and rate, together with the occurrence, characteristics, and relative depth of the tremor to be the result of magmatic activity beneath Taʻū Island volcano.

American Samoa, Taʻū Island↗
USGS2025

REDPy: A Python tool for automated repeating earthquake detection and visualization

Detecting and cataloging seismic events are among the most fundamental tasks in seismology. Many standardized tools for these tasks exist, including the open‐source package repeating earthquake detector in Python (REDPy). REDPy generates an organized catalog of seismic events from continuous waveform data, in which events are automatically separated into groups (“families”) by their waveform similarity through cross‐correlation. REDPy also automatically generates various outputs that allow a user to visualize important trends in the catalog, which may be used in real time or in retrospective analyses to allow rapid identification of interesting features. The code was designed for near‐real‐time volcano monitoring but is applicable across a broad range of use cases in seismology and seismoacoustics. In this article, the utility and performance of REDPy are demonstrated on two highly seismogenic volcanic eruption sequences: the onset of the dome‐building eruption of Mount St. Helens, Washington, from 2004 to 2005, and the entirety of the summit caldera collapse sequence of Kīlauea, Hawai‘i, in 2018. This article is meant to be a companion to the documentation of the code; in addition to detailing the basic required inputs, script functionality, and resulting outputs, the reasonings behind several important design decisions are also discussed.

Seismological Research Letters↗
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

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↗
USGS1992

Public water-supply systems and water use in Tennessee, 1988

This report summarizes the results of a study conducted by the U.S. Geological Survey, in cooperation with the Tennessee Department of Environment and Conservation (TDEC), Division of Water Supply in 1988. Data gathered during an inventory by the TDEC were collated to determine water use, supply sources, population served, and design and storage capacities of the systems. The inventory was limited to systems that were active on June 30, 1988. Results of a survey of the systems conducted by the Tennessee Department of Health and Environment during 1988 were a primary source of data for this report. Data from computer and manual files maintained by the Tennessee Department of Health and Environment and the U.S. Geological Survey also were used. The Division of Water Supply, TDEC, surveyed 541 public water-supply systems. These systems served 81 percent of the population of the State, or 3.95 million people. The gross per capita use statewide for public-supplied water was 179 gallons per day. Total water withdrawals for public supply increased about 39 percent from 510 million gallons per day (Mgal/d) in 1980, to 708 Mgalld in 1988. During the same period, the population increased about 7 percent. Surface-water withdrawals accounted for 63 percent (446 Mgal/d) of the total water withdrawn in the State. All of these withdrawals occurred in the Tennessee (56 percent or 249 Mgal/d) and the Ohio (44 percent or 197 Mgalld) hydrologic regions. Ground water supplied 262 Mgal/d or 37 percent of the total water withdrawn by public-supply systems statewide. Of that amount, 79 percent, or 208 Mgalld, was used in western Tennessee.

Tennessee↗
USGS1995

Methods for estimating the magnitude and frequency of peak discharges of rural, unregulated streams in Virginia

Methods are presented for estimating the peak discharges of rural, unregulated streams in Virginia. A Pearson Type III distribution is fitted to the logarithms of the unregulated annual peak-discharge records from 363 stream-gaging stations in Virginia to estimate the peak discharge at these stations for recurrence intervals of 2 to 500 years. Peak-discharge characteristics for 284 unregulated stations are divided into eight regions based on physiographic province, and regressed on basin characteristics, including drainage area, main channel length, main channel slope, mean basin elevation, percentage of forest cover, mean annual precipitation, and maximum rainfall intensity. Regression equations for each region are computed by use of the generalized least-squares method, which accounts for spatial and temporal correlation between nearby gaging stations. This regression technique weights the significance of each station to the regional equation based on the length of records collected at each cation, the correlation between annual peak discharges among the stations, and the standard deviation of the annual peak discharge for each station. Drainage area proved to be the only significant explanatory variable in four regions, while other regions have as many as three significant variables. Standard errors of the regression equations range from 30 to 80 percent. Alternate equations using drainage area only are provided for the five regions with more than one significant explanatory variable. Methods and sample computations are provided to estimate peak discharges at gaged and engaged sites in Virginia for recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years, and to adjust the regression estimates for sites on gaged streams where nearby gaging-station records are available.

Virginia↗
USGS2024

Simulated effects of projected 2014–40 withdrawals on groundwater flow and water levels in the New Jersey Coastal Plain

Abstract Groundwater flow between 2014 through 2040 was simulated in the New Jersey Coastal Plain based on three withdrawal scenarios. Two of the scenarios were based on projected population trends and the assumption of water conservation; the nominal water-loss scenario projected a status quo in the efficiency of water loss in the delivery systems whereas the optimal water-loss scenario projected a better water-loss efficiency resulting in less withdrawals. The third scenario assumes that all wells will withdraw water at their full allocation level which is generally much more than reported withdrawals in 2013 or projected under the other two scenarios. Maps and summaries of heads and drawdowns are presented for nine confined aquifers. All the aquifers have areas with heads below sea level by 2040. Of the three scenarios, the drawdowns are most extreme in the full allocation scenarios; there are large areas of head decline greater than 20 feet in 5 of the 9 confined aquifers. The exceptions are the Vincentown aquifer, despite some areas of large drawdown in the vicinity of wells, and the three Potomac-Raritan-Magothy (PRM) aquifers where withdrawals are regulated by Critical Area restrictions. The nominal and optimal water-loss scenarios have some areas of head declines; most are less than 15 feet. The simulation of these scenarios shows some extensive areas of head recovery as well—especially in the aquifers that are regulated by the Critical Area restrictions. Budgets of inflow and outflow components were calculated for 44 hydrologic budget areas (HBAs). The budget analysis shows that the water movement is complex and varies based on the aquifer geometry and location of pumping wells. Flow components between the unconfined and confined parts of the system were summarized by HUC11 (hydrologic unit code 11) basins.

New Jersey↗
USGS2024

Status and understanding of groundwater quality in the Mojave Basin Domestic-Supply Aquifer study unit, 2018—California GAMA Priority Basin Project

Groundwater quality in the western part of the Mojave Desert in San Bernardino County, California, was investigated in 2018 as part of the California State Water Resources Control Board Groundwater Ambient Monitoring and Assessment Program Priority Basin Project. The Mojave Basin Domestic-Supply Aquifer study unit (MOBS) region was divided into two study areas—floodplain and regional—to assess differences between the two major aquifers used for drinking water supply in the area. This assessment characterized the quality of ambient groundwater and not the quality of treated drinking water. The study included three components: (1) a status assessment, which characterized the quality of groundwater resources used for domestic drinking-water supply in the floodplain and regional study areas; (2) a brief understanding assessment, which evaluated factors that could potentially affect the quality of groundwater used by domestic wells in the region; and (3) a comparative assessment between the groundwater resources used by domestic wells and public-supply wells in the two study areas. The domestic-well assessment was based on data collected by the U.S. Geological Survey from 48 domestic wells in January–May 2018. The public-supply assessment was based on data for samples from 322 public-supply wells in 2008–18, either collected by the U.S. Geological Survey or compiled from the California State Water Resources Control Boards Division of Drinking Water publicly available database. Concentrations of water-quality constituents in ambient groundwater were compared to regulatory and non-regulatory benchmarks typically used by the State of California and Federal agencies as health-based or aesthetic standards for public drinking water. Relative concentrations, defined as the measured concentration divided by the benchmark concentration, were classified as high (greater than 1.0), moderate (greater than 0.5 for inorganic constituents or 0.1 for organic and special-interest constituents, and not high), or low (concentrations lower than moderate). The floodplain and regional study areas were divided into 15 and 35 grid cells, respectively, and grid-based methods were used to compute the areal proportions of the two study areas with high, moderate, or low relative concentrations of individual constituents and classes of constituents. For the domestic-supply assessment, one or more inorganic constituents with health-based benchmarks were detected at high relative concentrations in 58 percent of the regional study area and 13 percent of the floodplain study area. The inorganic constituents with health-based benchmarks detected at high relative concentrations in the regional study area were arsenic, chromium and hexavalent chromium, fluoride, adjusted gross alpha particle activity, uranium, molybdenum, strontium, and nitrate; only arsenic was detected at high relative concentrations in the floodplain study area. One or more inorganic constituents with secondary maximum contaminant level benchmarks were detected at high concentrations in 15 and 6.7 percent of the regional and floodplain study areas, respectively. The constituents detected at high relative concentrations in the regional study area were total dissolved solids, chloride, sulfate, and iron; only total dissolved solids and sulfate were detected at high relative concentrations in the floodplain study area. Organic constituents were not detected at moderate or high relative concentrations in either the regional or floodplain study areas. Volatile organic compounds were detected at low relative concentrations in 21 and 27 percent of the regional and floodplain study areas, respectively, and pesticides were detected at low relative concentrations in 9.1 and 20 percent of the regional and floodplain study areas, respectively. The only individual organic constituent detected in more than 10 percent of either study area was the trihalomethane trichloromethane. Total coliform bacteria were detected in 15 and 27 percent of the grid wells in the regional and floodplain study areas, respectively. The greater prevalence of high relative concentrations of many inorganic constituents in the regional study area compared to the floodplain area likely indicates the greater diversity of geologic material at depth in aquifer material and generally finer-grained alluvium compared to the floodplain study area combined with generally older groundwater that has had more contact time with aquifer materials. In general, trace element concentrations (1) increased with increasing groundwater age, (2) increased with distance from recharge sources in the mountains, and (3) increased with closer proximity to some types of geological units. In general, groundwater from domestic wells in the floodplain study area is young, with most samples containing a component of modern groundwater based on tritium and unadjusted carbon-14 activities, whereas groundwater from domestic wells in the regional study area generally is old, with most samples having unadjusted carbon-14 ages of 5,000–40,000 years. Public-supply wells in MOBS generally were deeper than domestic wells and presumably are in contact with older, more weathered alluvium that may have more mobile trace elements, such as arsenic or uranium. However, only 26 percent of the public-supply regional study area had high relative concentrations of inorganic constituents, compared to 58 percent for the domestic regional study area. The percentages of the public-supply and domestic floodplain study areas with high relative concentrations of inorganic constituents were 11 and 13 percent, respectively. The ages of groundwater used by public-supply and domestic wells in each study area were similar, which was not expected given the greater depth of the public-supply wells. Three potential factors may contribute to these results: (1) greater spatial footprint of domestic well network, which may result in domestic wells pumping groundwater from fractured bedrock or mineralized areas not used by public-supply wells; (2) greater pumping rates in public-supply wells, resulting in more water being withdrawn from coarse-grained, heterogeneous alluvium than finer-grained layers, which may have higher concentrations of (or more mobile) inorganic constituents; and (3) a greater degree of well management with public-supply wells, which may include pausing use of or decommissioning wells if treating or blending water is not feasible to lower constituent concentrations.

California↗
USGS2025

Lessons in business recovery following the 2023 Kahramanmaraş earthquake sequence, Türkiye informed by women entrepreneurs

On 6 February 2023, Southern Türkiye was hit by devastating earthquakes, directly affecting over 14 million people in 11 cities, causing more than 50,000 deaths and the destruction of more than 800,000 buildings. This article goes beyond the physical damage imposed by the catastrophe to discuss the effects of the earthquakes on the operations of women-owned businesses. The mixed-method study with entrepreneurs belonging to a women’s business association operating in a moderately disrupted part of the region explores their struggles and recovery expectations. Thirty-five questionnaires were analyzed to identify the reasons for business closure, challenges, and needs faced in the post-disaster period and their recovery strategies. In addition, 23 entrepreneurs participated in roundtable discussions to provide a broader context to their responses to survey topics as well as lessons learned. Across both the survey and roundtables, while many respondents reported minor physical damage to their building, they also experienced financial and personal challenges from disruption to equipment, infrastructure, services, supply chains, institutional decisions, employee well-being, and customer base. Many used their business resources and personal savings to assist employees and others in the community. The women entrepreneurs often felt their recovery needs were ignored by government and private relief organizations and encountered barriers to receiving assistance from public and private institutions. Organizing together as women in business, even informally, provided mutual support during the crisis and recovery periods and catalyzed their role in support of their communities. The results illuminate functional community recovery as a balance of recovery of built infrastructure functionality and recovery of the broader social and economic fabric of the community.

southern Turkey↗
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

Decadal-scale characteristics of natural and anthropogenic dune morphology along North Carolina Barrier Islands (SE Atlantic Coast)

Dunes offer protection to coastal communities from storms, yet quantitative research characterizing human influences on dune systems across decades and 10–100 km spatial scales constitutes a known data gap. This study analyzed a multidecadal archive of dune crest elevations and positions to map mesoscale spatiotemporal trends in foredune morphology across the barrier island coastline of North Carolina, SE Atlantic coast. Results are summarized within four categories of coastline representing a range of coastal management intensities, integrating beach width and dune lateral accommodation space to assess the comparative influence of differing intensities of anthropogenic intervention on dune morphology. These categories account for the presence of structures, roads, and alongshore connectivity to anthropogenic sediment sources (beach nourishment). These categories show distinct dune and beach system characteristics. Dunes are largest in the areas with the longest legacy of coastal management. As categories decrease in development intensity, dune and beach system characteristics generally depart from the large, static dunes and narrow beaches consistent with hold-the-line management. Where small or transgressive dunes exist on heavily developed coasts, they coincide with relatively strong hydrodynamic forcing or sediment-limiting geology. There is evidence of indirect coastal management influence on undeveloped barrier islands with alongshore connectivity to developed barriers. While this analysis was conducted at the mesoscale, the results provide researchers with the framework to make insights at smaller scales than those demonstrated in this study.

North Carolina↗
USGS2026

Stress states on the eve of past earthquakes inform earthquake rupture through fault complexity along the San Andreas and San Jacinto faults

Estimating the evolving state of stress along active fault systems can provide insight into the conditions that generated past ground‐rupturing earthquakes and influenced their ability to propagate through areas of geometric complexity, such as fault branches and stepovers. We use quasi‐static forward numerical models that incorporate the 3D complex configuration of active faults in southern California to estimate shear tractions on the geometrically complex southern San Andreas and San Jacinto faults from 1000 to 1900 C.E. These tractions include interseismic accumulation of traction due to tectonic loading, viscoelastic relaxation of shear stress within the upper crust between earthquakes, and effects of other earthquakes on the fault network. We simulate ground‐rupturing earthquakes based on the along‐strike earthquake extents modeled by Scharer and Yule (2020) , assuming that stress drop is complete in each earthquake. We use Monte Carlo simulations to estimate uncertainty in evolving shear tractions due to uncertainties in earthquake timing and in upper‐crustal viscosity. Pre‐earthquake shear tractions typically do not exceed ∼2 MPa. Although ruptures with length <200 km have pre‐earthquake shear tractions that range from near zero to ∼1.75 MPa, these tractions are not less than ∼0.4 MPa for earthquakes with rupture length >200 km. Earthquakes with long (>200 km) ruptures occur only in the single‐stranded part of the system, whereas those with short (<125 km) rupture length and high pre‐earthquake shear traction occur near fault stepovers and branches. This suggests that high accumulated shear traction encourages longer rupture propagation, but may not be sufficient to overcome geometric complexities. This modeling approach informs our understanding of rupture propagation and provides estimates of fault shear tractions that are unavailable from direct measurements.

Callifornia↗
USGS2026

Ductile and brittle Rio Grande Rift deformation in Oligocene granite records a two-stage rift history in southern Colorado

The timing and nature of early deformation in the Rio Grande Rift remains poorly constrained. We present evidence for the earliest structural signature of rift extension in the Sangre de Cristo Range, southern Colorado, based on new geologic mapping, structural analysis, rock magnetic data, and thermochronology. These analyses focus on the ~30.0 Ma granite of Chokecherry Canyon, which hosts discrete low-angle mylonitic shear zones and a distributed, gently SW-dipping protomylonitic fabric. Incremental stretching axes, stretching lineations, and Kmax magnetic lineations plunge gently WSW. Quartz microstructures and crystallographic orientations indicate dominantly coaxial strain in the protomylonite and general shear in the discrete shear zones. Quartz c-axis opening-angle thermometry suggests deformation at ~420–540°C. Thermal modeling of ⁴⁰Ar/³⁹Ar K-feldspar data indicates rapid post magmatic cooling below the brittle–plastic transition, supporting shear-zone formation immediately after emplacement. Slow cooling from ~20–13 Ma was followed by renewed rapid cooling at ~13 Ma, interpreted as the onset of extensional exhumation along the Sangre de Cristo Fault System. These results show that extension in the northern Rio Grande Rift was active by ~30 Ma, earlier than previously recognized. We propose a two-stage model for northern Rio Grande Rift evolution: Stage I (30–23 Ma) records ENE–WSW extension localized in low-angle mylonitic shear zones associated with mid-crustal intrusions; Stage II (≤18 Ma) reflects brittle high-angle normal faulting, focused exhumation, and rift narrowing. Stage I magmatism and deformation along the western range front likely established crustal weaknesses that guided later fault development.

Colorado↗
USGS2025

Geochemical monitoring of volcanic fluids in the twenty-first century

This chapter reviews the state-of-the-art of geochemical volcano monitoring techniques. We cover in-situ monitoring technologies that rely on sampling fluids (direct sampling) and on instrumental analysis of the composition of such fluids in real-time (remote sensing of volcanic fluids is covered elsewhere in this book). We first review key concepts and principles in the field, and then review the results of some selected case studies and applications. We cover the large variety of fluid categories emitted by volcanoes, in both the near-(crater fumaroles and lakes, and plumes) and far-(degassing soils, groundwaters) fields. Our aim is to demonstrate the utility of measuring the chemistry of fluids released by volcanoes, and how these can help characterize volcano unrest, and eventually the increased likelihood of eruption. We conclude with a brief discussion of current challenges and knowledge gaps, and on future directions in geochemical monitoring.

Book chapter↗
USGS2024

Onset of aftershocks: Constraints on the Rate-and-State model

Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R (t)=K(c+t) −p ⁠ , with parameters K , c , and p . The rate‐and‐state (RS) model, which is currently the most popular physics‐based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c ‐value that depends on the size of the coseismic stress change. Because the mainshock‐induced stresses strongly vary in space, the c ‐value should vary accordingly. Short‐time aftershock incompleteness (STAI) in earthquake catalogs has prevented a detailed test of this prediction so far, but the newly developed a ‐positive method for reconstructing the true earthquake rate now allows its testing. Using previously published slip models, we calculate the coseismic stress changes for the six largest mainshocks in Southern California in recent decades and estimate the maximum shear as a scalar proxy of the coseismic stress tensor. Aftershock rates reconstructed for events in different stress ranges show that the rates follow a power law with p = 1 independent of stress with no clear sign of a c ‐value. The onset of the power‐law decay is abrupt and more delayed in areas with smaller stress changes. The observations do not necessarily contradict the RS model, as STAI limits the resolution for early aftershocks, and the RS model can reproduce the observations for specific Aσ values. However, the observations lead to strong constraints, namely Aσ <10 kPa and a power‐law decay of the background rate with distance to the fault, with exponent 2.7.

Seismological Research Letters↗