Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Water Process Engineering”

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

Refining the earthquake history of south-central Alaska through lake records

The Alaska–Aleutian subduction zone (AASZ) is one of the world's most seismically active plate boundaries and the source of the 1964 M w 9.2 Great Alaska earthquake–the second largest instrumentally recorded earthquake in the world. Understanding the nature and frequency of such earthquakes is necessary for seismic and tsunami hazard assessment, but instrumental and historical records that span less than 150 years are too short to allow a statistically reliable analysis of earthquake recurrence times. This calls for studies of evidence of past earthquakes, extending the earthquake catalog further back in time. Subduction-zone paleoseismology in south-central Alaska is predominantly based on coastal evidence of land-level changes and tsunamis generated by megathrust earthquakes and preserved in the geological record. A complementary approach is lacustrine paleoseismology, which is still a relatively young discipline in Alaska. However, globally, lake basins are well-established high-resolution and continuous recorders of paleoseismic activity along subduction zones, relying on the identification of underwater landslide deposits and turbidites generated by seismic shaking. As a result, lake basins not only register ground shaking from megathrust earthquakes, but also from intraslab and crustal earthquakes, which are typically not accompanied by significant land-level changes. In this review paper, we combine coastal and lacustrine paleoseismology approaches to refine the south-central Alaskan earthquake history by comparing the paleoseismic records from two lakes (i.e., Eklutna Lake, located in the Chugach Mountain Range, and Skilak Lake, situated on the Kenai Peninsula) with the coastal and crustal earthquake catalog in Alaska. The resulting age ranges of all known megathrust earthquakes involving the Alaskan megathrust between the Kodiak and Prince William Sound (PWS) sections are more precise and accurate for the last 1.3 kyrs BP than the previously published age ranges from coastal records. As a result, this study supports the following key conclusions: (1) The 1964 CE earthquake was an exceptionally strong and unique event in the last 2000 years, rupturing the PWS, Kenai, Barren Islands, and Kodiak sections simultaneously. (2) The high-resolution and seasonal markings of the varved lake records now disentangle for the first time closely timed earthquakes, which was not possible based on the coastal evidence alone. (3) No persistent megathrust rupture boundaries exist. So, the possibility of a full rupture of the entire eastern AASZ, from PWS to Semidi cannot be excluded. (4) The rupture pattern in the eastern AASZ reveals superimposed cycles of multi-asperity ruptures (1964 earthquake) and clustered complementary partial ruptures, or rupture cascades. (5) The PWS section hosts the largest asperity in the eastern AASZ. (6) The shaking record of megathrust earthquakes indicates a time-dependent (quasiperiodic) behavior for the study area, but the observation of complementary clusters means that the hazards will not drop to zero but instead may even increase for a neighboring section. (7) The time-independent behavior of intraplate earthquakes implies that the intraslab hazard did not decrease following the 2016 and 2018 earthquakes. This study utilizes an integrated approach for subduction zone paleoseismology as a solution for unraveling recurrence and rupture patterns in Alaska, which can be applied worldwide.

Alaska↗
USGS2025

Leveraging wildfire to augment forest management and amplify forest resilience

Successive catastrophic wildfire seasons in western North America have escalated the urgency around reducing fire risk to communities and ecosystems. In historically frequent-fire forests, fuel buildup as a result of fire exclusion is contributing to increased fire severity. The probability of high-severity fire can be reduced by active forest management that reduces fuels, prompting federal and state agencies to commit significant resources to increase the pace and scale of fuel reduction treatments. However, lower severity areas of wildfires also have the potential to act as “treatments,” and even catastrophic fires with large areas of high severity can still have substantial areas of lower severity fire that may be improving forest conditions locally. We quantified active management and wildfire severity across yellow pine and mixed conifer (YPMC) forests in the Sierra Nevada of California over a 22-year period (2001–2022). We did not detect increases in the area treated through time, but the area of beneficial wildfire (low to moderate severity) increased substantially, exceeding active treatment area in 8 of 22 years. Overall, beneficial wildfire treated ~17% more area than all treatments combined, and roughly four times more area than fire-related treatments alone. We then used disturbance history to evaluate resistance to high-severity wildfire and forest loss across the YPMC range. Of the 2.3 million ha YPMC of forests in 2001, 20% lost mature forests due to high-severity fire by 2022, which is nearly half of all YPMC area burned. Most of the landscape (47%) remains at risk of high-severity fire because it had no restorative disturbances, but 33% of the study area has some level of resistance to high-severity wildfire. In these areas, resistance will need to be enhanced and maintained over time via active management or managed wildfire, but these treatment needs will likely outpace capacity even under optimistic implementation scenarios. Given limited resources for implementing active management and the likelihood of a more fiery future, incorporating beneficial wildfire into landscape-level treatment planning has the potential to amplify the impact of active management treatments.

California, Nevada↗
USGS2026

Spatial variation in effectiveness and costs of sagebrush restoration treatments across the western United States

Natural resource managers strive to improve restoration efficiency across the vast network of public lands of the United States, seeking the best return on investment through innovation, geographic prioritization, and adaptive management. This challenge is exemplified by ecosystems dominated by sagebrush ( Artemisia spp.), which occur across the western United States and have been degraded by increasing wildfire frequency driven by invasive plant spread, among other factors. However, both costs and effectiveness of sagebrush restoration can vary spatially due to management practices in addition to biotic and abiotic factors, and characterizing this variation could inform broad-scale planning. We leveraged published models of sagebrush recovery and treatment costs to predict and evaluate cost-effectiveness of aerial and drill seeding Artemisia spp. across 429 718 km 2 that overlapped 12 US states. Compared to natural recovery, effectiveness for both methods was generally low, yet effectiveness was often greater, and costs were always greater, for drill than aerial seeding, resulting in slightly higher cost-effectiveness on average for drill seeding. Cost-effectiveness for both treatments increased near major roads and in areas with more repeated burns. Cost-effectiveness also increased with predisturbance cover of sagebrush but decreased with herbaceous cover. However, we also identified areas where aerial seeding was more cost-effective than drill seeding, particularly with greater slope. Our results provide spatially explicit estimates and potential mechanisms of cost-effectiveness for two common seeding methods for sagebrush, which can help prioritize limited resources, guide land use, and improve restoration efficiency and effectiveness across public rangelands of the western United States.

western United States↗
USGS2026

A temporal look at the influence of topographic amplification on earthquake-triggered landslides in 3D seismic simulations

Earthquakes are a primary trigger for landslides, often leading to catastrophic consequences. While numerous studies have explored the spatial distribution of earthquake-triggered landslides, understanding the interaction between seismic waves and topography remains a critical challenge. Topographic irregularities can cause seismic wave amplification, altering ground shaking, and can trigger landslides that are challenging for predictive models to anticipate. This study investigates the spatial and temporal evolution of topographically amplified landslides, focusing on coseismic landslides triggered by the Mw 7.5 mainshock of the 2018 Papua New Guinea earthquake and post-seismic landslides associated with its four aftershocks, each exceeding Mw 6.0. We employ low-frequency, three-dimensional numerical ground shaking simulations and data-driven multivariate analyses to examine how landslides evolved from the coseismic to post-seismic periods. Our findings reveal a spatial shift in landslide distribution, in which the mainshock triggered slope failures predominantly on steep hillslopes, and the aftershocks triggered landslides on gentler slopes, often near geologic boundaries. We attribute this transition partly to the earthquake legacy effect of the mainshock, where the mainshock caused weakening of these hillslopes, making them more prone to failure when aftershocks occur. Additionally, the concentration of failures along geologic contacts in the post-seismic phase suggests that site amplification, stemming from contrasts in subsurface materials, exerts a key influence on landslide occurrence. Although not explicitly captured in our current numerical simulations, this mechanism warrants further investigation for more accurate hazard modeling.

JGR Solid Earth↗
USGS2025

Seismic moment and local magnitude scales in Ridgecrest, CA from the SCEC/USGS Community Stress Drop Validation Study

We illustrate the systematic difference between moment magnitude and local magnitude caused by underlying earthquake source physics, using seismic moments submitted to the Statewide California Earthquake Center/United States Geological Survey Community Stress Drop Validation Study 2019 Ridgecrest data set. While the relationship between seismic moment and moment magnitude ( M or M w ) of log 10 ( M 0 ) ~ 1.5* M is uniformly valid for all earthquake sizes by definition (Hanks and Kanamori, 1979), the relationship between local magnitude M L and moment is itself magnitude dependent. For moderate events, ~3< M < ~6, M and M L are coincident; for earthquakes smaller than ~3, M L ~ 1.0 log 10 M 0 (Hanks and Boore, 1984). This is a physical consequence of the corner frequency fc becoming larger than the upper frequency of observation and implies that M L and M differ systematically by a factor of 1.5 for these small events. While this idea is not new, we propose a new, continuous relationship between local magnitude and moment, for magnitudes 2 to 6 which extrapolates to smaller and larger magnitudes, applicable to southern California specific to the Ridgecrest region. We make use of the plethora of seismic moments as submitted by many participants of the Community Stress Drop study, compared to the Southern California Seismic Network (SCSN) catalog magnitudes. Overall, the seismic moments in the Community Study recover moment magnitude well, so we use our new M L - M 0 to convert M L to M , refining the SCSN operational M Lr scale. This systematic difference of 50% in slope between local and moment magnitude at small magnitudes has implications for spectral stress drop estimates, earthquake ground motion modeling, as well as other magnitude scales and earthquake occurrence statistics.

California↗
USGS2026

Computing flow-field distortion coefficients from well-construction and formation properties

Direct measurements of groundwater velocity made with borehole flowmeters in screened wells must be compensated for the effects of flow-field distortion (also known as borehole acceleration). A theoretical equation developed by Drost et al. (1968) and simple inputs describing hydraulic properties of well construction and geologic formation were programmed into an Excel workbook to facilitate computation by groundwater-flowmeter users. Tables describing the physical and hydraulic properties for well constructions and gravel pack media are provided with an example to facilitate use of the workbook. Groundwater flowlines converge or diverge as they pass from a geologic formation, through a gravel pack and well screen. The extent of flowline convergence or divergence and the value of the flow-field distortion coefficient is related to the relative changes in hydraulic conductivity of the well screen, gravel pack, and geologic formation. Convergence or divergence is accompanied by acceleration or deceleration of groundwater. Direct measurements of groundwater velocity at the center of the monitoring well can be adjusted to provide a more accurate estimate of velocity in the formation by applying a correction for flow-field distortion. Variables required to compute the flow-field distortion coefficient include the hydraulic conductivity of the gravel pack, well screen, and the geologic formation surrounding the well screen; the borehole radius, and the inside radius and outside radius of the well screen.

Groundwater↗
USGS2024

Assessment and characterization of ephemeral stream channel stability and mechanisms affecting erosion in Grand Valley, western Colorado, 2018–21

The Grand Valley in western Colorado is in the semiarid Southwest United States. The north side of the Grand Valley has many ungaged ephemeral streams, which are of particular interest because (1) the underlying bedrock geology, Late Cretaceous Mancos Shale, is a sedimentary rock deposit identified as a major salinity contributor to the Colorado River and (2) despite infrequent streamflows of short duration, monsoon-derived floods in these ephemeral streams can carry substantial amounts of sediment downstream, affecting upstream and downstream banks and channel cross sections. The study area is of interest, because salinity, or the total dissolved solids concentration, in the Colorado River causes an estimated $300 million to $400 million per year in economic damages in the United States, and it is estimated 62 percent of the Upper Colorado River Basin’s total dissolved solid loads originate from geologic sources. In an effort to minimize salt contributions to the Colorado River from public lands administered by the Bureau of Land Management, a comprehensive salinity control approach is typically used to reduce nonpoint sources of salinity through land management techniques and practices. In 2018, the U.S. Geological Survey, in cooperation with the Bureau of Land Management, began an assessment of ephemeral streams located on the north side of the Grand Valley, western Colorado, to characterize stream channel stability and identify mechanisms affecting erosion. The U.S. Geological Survey developed a method for automatically extracting channel cross-section geometry from existing remotely sensed terrain models. Based on estimated flood stage and surrogate streamflows, hydraulic characteristics were calculated. Furthermore, the channel geometries and hydraulic characteristics were used to estimate channel stability using a statistical model. Cross-section stabilities were determined from a stream channel stability assessment for a subset of 1,406 visited (field observed) locations out of 13,415 cross sections, which were delineated from remotely sensed terrain models. The application of Manning’s resistance equation in combination with multiple logistic regression models demonstrated channel stability can be estimated with a 0.845 goodness of fit for a validation dataset when using a combination of drainage area, width-to-depth ratio, sinuosity, and shear stress as the explanatory variables. Stream channel stability was extrapolated for 13,415 unvisited (not field observed) cross sections using the multiple logistic regression model and defined explanatory variables. Mapping of the ephemeral streams and their associated stabilities may be used by the Bureau of Land Management to prioritize areas for remediation or changes in management strategies to reduce sediment and salinity loading to the Colorado River. The study found channel stability within the ephemeral streams to be spatially variable, longitudinally discontinuous, and dictated by changes in channel bed slope. The stable ephemeral streams were relatively wide and shallow and often had smaller drainage areas with less potential for producing shear stresses capable of overcoming channel adhesion. A change in channel bed slope can provide the means necessary to generate shear stresses appropriate to initiate erosion and a subsequent stability transition to incising channels. Channel widening happens when either or both banks of an incising channel reach a critical height for mass wasting, or when channel curvature causes higher sidewall stress. Regardless, widening channels can promote increases in sinuosity and subsequently reduce steep channel bed slopes. Consequently, stable and widening channels can have comparable bed slopes, making channel bed slope a poor explanatory variable to predict channel stability overall, despite its function to initiate channel instability. The results were based on a surrogate 0.10 annual exceedance probability (AEP; return period equal to the 10-year flood) interval streamflow, although it was recognized fluctuations in streamflow would also affect channel stability. Past and current changes within the study area affect streamflow; therefore, mechanisms affecting erosion include land use disturbances, soil compaction, loss of vegetation cover, drought, less frequent and more extreme precipitation, and fires—which all intensify the potential runoff and erosion within the study area.

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

Benthic habitat map of Olowalu Reef, Maui, Hawaii—Geomorphological structure, biological cover, and geologic zonation determined with spectral, lidar, and acoustic data

The fringing coral reef off Olowalu, Maui, Hawaii, has been identified as a local conservation priority site. In 2007, the National Oceanic and Atmospheric Administration (NOAA) produced a benthic habitat map of the Hawaiian Islands that was used as a foundation for this study. To support place-based management of the reef in the future, the U.S. Geological Survey (USGS) mapped the geologic zone, major and dominant geomorphological structure, biological cover type, and percent of biological cover for 11 square kilometers (km 2 ) of Olowalu Reef at a minimum mapping unit (MMU) of 100 square meters (m 2 ) to create a benthic habitat map. Heads-up digitization was employed on 0.50-meter (m) natural color satellite orthoimagery with ancillary 1-m acoustic backscatter imagery from single-scan sonar (sound navigation and ranging). A 1-m, 4-m, and 8-m digital bathymetric model (DBM) was interpolated from bathymetric lidar (light detection and ranging), and various geomorphometric layers derived from the DBMs were used for habitat interpretation. Still-frame imagery of the seafloor extracted from vessel-towed underwater video transects on Olowalu Reef served as ground validation points ( n =870) during active mapping and accuracy assessment points ( n =216) for thematic accuracy assessment. Thematic accuracy was cross-validated by the Hawai‘i Department of Land and Natural Resources Division of Aquatic Resources. Final thematic accuracy was 88.8 percent for major structure, 85.6 percent for dominant structure, 86.0 percent for major biological cover, and 78.6 percent for type and percent of major biological cover. Reef and hardbottom constituted 52 percent of the total mapped habitat, comprising mostly aggregate reef (31 percent) and pavement (11 percent), with large swaths of spur-and-groove (9 percent). Of this hardbottom, 17 percent was covered with moderate (10 to <50 percent) coral and 27 percent with high coral cover (50 to <90 percent). High (50 to <90 percent) macroalgae cover dominated the continuous sand sheets in offshore bank/shelf zones. The map created in this study supplements the NOAA 2007 map and expands on the observations made by USGS sampling of the reef. The NOAA 2007 map and our map differed in total areal extent by a negligible 6 m 2 and were in general thematic agreement. Our map is intended to serve as a baseline for public access, general research, local-level management, and reef change for future studies.

Hawaii↗
USGS2024

Shallow storage of the explosive Earthquake Flat Pyroclastics magma body, Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand: Evidence from phase-equilibria experiments

Rhyolitic tuffs range widely in their crystal contents from nearly aphyric to crystal-rich, and their crystal cargoes inform concepts of upper crustal magma reservoirs. The Earthquake Flat pyroclastics (Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand) are 10 km 3 of rhyolitic tuffs with abundant (~ 40 vol.%) plagioclase and quartz, minor biotite, hornblende, and orthopyroxene, and accessory Fe-Ti oxides, apatite, and zircon, set in high-silica rhyolitic glass. Major minerals form large, euhedral phenocrysts and abundant glomerocrysts with few disequilibrium textures excepting some faintly resorbed quartz. Plagioclase phenocrysts have thick rims of nearly constant composition near An 30 , and hornblende is weakly zoned or unzoned. The abundant and texturally complex mineral assemblage contrasts with the nearby (~ 25 km), nearly synchronous, but more voluminous and crystal-moderate rhyolite tuffs from Rotoiti caldera. New H 2 O-saturated phase-equilibria results on the erupted Earthquake Flat melt (glass) determine its co-saturation with the partial phenocryst assemblage of plagioclase, quartz, biotite, and Fe-Ti oxides at: 140 MPa, 755 ºC. These closely approximate the conditions of the pre-eruptive magma body assuming it was saturated with nearly pure H 2 O and at an f O 2 of ~ Ni–NiO. Absence of hornblende and orthopyroxene from the synthesized assemblages may result from those minerals being in a peritectic reaction relation with melt to produce biotite, so they would not grow from the liquid used as starting material. Experimental results on Rotoiti rhyolite (Nicholls et al. 1992) show that the two bodies resided at similar pressures, temperatures, and f O 2 s. Lower crystal abundance of the Rotoiti tuffs may result from slight compositional differences. We interpret that the Earthquake Flat pyroclastics were sourced from the crystal-rich periphery of a mushy reservoir system with the Rotoiti occupying a more melt-rich central location. Uncertain is whether this was a single intrusion zoned continuously in crystallinity, or discrete adjacent intrusions, but our results illustrate and quantify complexities of magma storage across relatively short distances.

North Island, Taupo Volcanic Zone↗
USGS1976

Uranium and trace elements in stream sediments as an exploration tool

More than 45 trace elements have been reported in anomalous amounts in the uranium ores. The specific suite of elements associated with any one uranium deposit varies according to deposit type and geologic province. The primary geochemical halo of uranium and associated trace elements in the host rock, together with secondary dispersion halo in soils and alluvium, offers a potential geochemical exploration target. Sediment from streams near low-grade uranium occurrences in arkosic sandstone of the Denver basin, Colorado, and sediment from streams in the igneous and metamorphic terrace near the Midnite mine of eastern Washington have been analyzed for acid-extractable U, Cu,Pb, Zn, Ni, Co, Fe, Mn, Mo, Cd, Cr, Ag, V, and Se, and for total Fe, Mn, Hg, As, and organic C. The -80 mesh fraction of the alluvium downstream from a small uraniferous limonite occurrence in the Denver basin shows anomalous concentrations of U, Pb, As, Ni, Mo, V, Zn, and Se. Anomalous concentrations of U persist farther downstream (as much as 1,000 m) than any of the trace elements. In another stream system southeast of Denver, above-background levels of U, Pb, As, Zn, V, Fe, and Mn occur near a uranium anomaly (defined by other methods) which has no surface mineralization. In the Midnite mine area, sediments of streams draining prospect areas showed anomalous concentrations of uranium (as much as 10 times background), but no accompanying anomalous trace-element concentrations, although analyses of Midnite mine ores show that a broad suite of trace elements are associated with uranium.

Open-File Report↗
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↗
USGS2024

Preliminary ground and airborne-based geophysical mapping and modelling of an active hydrothermal system at Mammoth Lakes, California

Mammoth Lakes, California hosts a productive hydrothermal system within the seismically active south moat of Long Valley Caldera. Surficial evidence of the shallow hydrothermal system includes discrete zones of tree-kill dispersed between Shady Rest Park and the Casa Diablo Geothermal Power Plant (40 MW), as well as east of the power plant. The tree-kill areas are associated with elevated diffuse CO₂ emissions, heated ground, hydrothermal alteration, diffuse soil H₂S emissions, and gas vents. Previous mapping delineates prominent north and northwest trending structures within the south moat along the southwestern edge of the resurgent dome that may accommodate gas and fluid flow at the Shady Rest Park and Basalt Canyon Tree Kill Areas (SRTKA and BCTKA, respectively). Both tree-kill areas are also located along contacts between resurgent rhyolite, mafic lavas, and surficial deposits which may provide additional pathways for gas and fluid migration in the shallow subsurface. Characterizing structure and lithology using geophysical anomalies is critical to determining primary structural controls on the hydrothermal system and the extent of subsurface alteration at these sites. We conducted ground and airborne-based potential field geophysical surveys to map gravity and magnetic anomalies. These anomalies are then used to model subsurface geology, structure, and hydrothermal alteration. Here we present our preliminary geophysical mapping and modelling results at both tree-kill locations. Gravity and magnetic data suggest complex structural intersections are coincident with heated ground and gas emissions at the SRTKA and BCTKA. Hydrothermal systems are often observed or interpreted to exploit fault intersections which can serve as highly permeable pathways for hydrothermal fluid and gas discharge, enabling economic geothermal energy production. Geophysical mapping and modelling are an effective means of investigating such structural complexity at Mammoth Lakes due to the presence of unidentified and concealed structures.

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

Core microbiomes as a potential fingerprinting method of Western USA dust sources

Introduction: Changing frequency and intensity of dust emissions impacts ecosystems and human health. Dust carries microbes, nutrients, heavy metals, and other materials that may change environmental biogeochemistry at deposition sites. Identifying dust sources provides key information on where and when mitigation strategies should be employed. However, commonly used geochemical or isotopic tracers are often not capable of distinguishing between geographic regions. Methods: We explored whether soil bacterial communities may provide distinct fingerprints of dust sources in the western United States. We identified bacterial core communities of dust from ten locations monitored by the National Wind Erosion Research Network (NWERN) with varied land use (cropland, rangeland, and playa), and compared communities to location, soil, and regional characteristics. Samples were collected monthly from Modified Wilson and Cooke (MWAC) samplers, composited by season (spring, summer, and fall), and analyzed using 16S rRNA sequencing. Results: We found distinct bacterial core communities that reflected dust source characteristics. In order of importance, precipitation levels ( p = 0.0001), location ( p = 0.0001), soil texture ( p = 0.0001), seasonality ( p = 0.0001), and elevation (p = 0.0002) were correlated with bacterial community composition. Discussion: Distinct bacterial core communities were associated with site characteristics such as biocrusts, playas, and military base proximity. Our results suggest that the use of core microbiomes may offer a fingerprinting method to identify dust source regions.

Colorado, Nevada, New Mexico, North Dakota, Oklaho↗
USGS2026

Rapid characterization of the 2025 Mw 8.8 Kamchatka, Russia earthquake

The 29 July 2025 M w 8.8 Kamchatka, Russia, earthquake was the sixth largest instrumentally recorded earthquake. This event was seismically well observed at regional and teleseismic distances, but publicly available near‐source data were sparse at the time of the event, presenting unique challenges for rapid source and impact characterization. The U.S. Geological Survey (USGS) National Earthquake Information Center provides global real‐time monitoring for earthquakes, including rapid response information products that estimate source characteristics, shaking, and the resulting impacts. We describe the USGS rapid response earthquake information products following the Kamchatka event and discuss their implications for ongoing hazards in the region. We describe potential improvements to our response workflows motivated by this event, including more rapid constraints on source geometries and the automated selection of fault geometries for finite‐fault inversions. The rapid response products together support the interpretation of a unilateral southwestward rupture with significant slip on the southwestern end of the rupture extent. The M w 8.8–9.0 event in 1952, which ruptured a comparable extent of the Kuril–Kamchatka subduction interface, has many similarities to the 2025 rupture. This illustrates that slip deficits may remain following great earthquakes and highlights the usefulness of comparative studies between historic and modern events.

Kuril-Kamchatka subduction zone↗
USGS2024

Mafic alkaline magmatism and rare earth element mineralization in the Mojave Desert, California: The Bobcat Hills connection to Mountain Pass

Occurrences of alkaline and carbonatite rocks with high concentrations of rare earth elements (REE) are a defining feature of Precambrian geology in the Mojave Desert of southeastern California. The most economically important occurrence is the carbonatite stock at Mountain Pass, which constitutes the largest REE deposit in the United States. A central scientific goal is to understand the genesis of the carbonatite ore body in the context of widespread REE-rich igneous activity. A swarm of mafic alkaline (shonkinite) dikes has been mapped and sampled at Bobcat Hills, 65 km southeast of the Mountain Pass mine. Whole-rock geochemistry and zircon geochronology demonstrate a clear affinity to the ca. 1.4 Ga Mountain Pass intrusive system. Bobcat Hills dikes have comparably high REE concentrations (La ∼1,000× chondritic) and an error-weighted mean 207 Pb/ 206 Pb zircon crystallization age of 1,426 ± 2 Ma (2 σ ). Unlike the alkaline intrusions at Mountain Pass, which have abundant inherited zircon from Paleoproterozoic basement rocks and crustally influenced oxygen isotope compositions (δ 18 O zircon = 6.5–7.5‰), the Bobcat Hills dikes lack any evidence of crustal assimilation and have oxygen isotope values that overlap a mantle range (Bobcat Hills average δ 18 O zircon = 5.6 ± 0.3‰). The dikes were a high-temperature, early center of mafic alkaline magmatism in the Mojave Desert that serve as a snapshot of melt generation from a spatially extensive, metasomatized mantle source. We propose that modification of the crust over many tens of Myr at Mountain Pass created an environment that favored crustal assimilation and enabled ascent of late-stage, REE-rich carbonatite magmas.

California↗
USGS2025

Landsliding follows signatures of wildfire history and vegetative regrowth in a steep coastal shrubland

Five years after the deadly and destructive 9 January 2018 Montecito debris flows (Santa Barbara County, California, USA), an atmospheric river storm on 9 January 2023 triggered widespread landsliding that affected many of the same drainages in the Santa Ynez Mountains. Using high-resolution aerial imagery, we identified >10,000 landslides over an ∼160 km 2 area. Most of the landslides were shallow (<1 m in depth) translational debris slides that initiated on steep (∼40°), south-facing hillslopes, with the highest incidence of landsliding in a sandstone-dominated bedrock unit. The landslides mobilized into debris flows and delivered substantial quantities of sediment downstream, which contributed to costly infrastructure impairments. We detected order-of-magnitude differences in landslide density across the study area that could not be attributed to variations in geomorphology (topographic aspect and slope), geology (bedrock type), or hydrology (seasonal antecedent rainfall, peak hourly storm rainfall intensity, total storm rainfall), which are usually considered relevant factors for shallow landsliding. Rather, we found that vegetation regrowth following wildfire was likely a relevant factor associated with the highly variable landslide densities. Hillslopes with less and different types of vegetation regrowth after fire appear to have been more susceptible to shallow landslides. We identify a possible vegetation control on postfire landsliding, which highlights an opportunity for hypothesis testing using more advanced techniques to track the evolution of vegetation cover and vegetation type in steep shrubland environments following wildfire.

California↗