SEARCH · Geology Reports

Results for “Journal of the American Water Resources Association”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,691 records · Page 13Linked to original sources

USGS2025

VIPER site analysis

We needed to evaluate available orbital data of NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) mission area in order to derive a variety of maps to help the science team identify scientifically interesting places for the rover to visit and to provide scientific context for our mission. Some of these maps also fulfilled engineering and mission design needs to enable safe and efficient landing and roving. We incorporated data from the Lunar Reconnaissance Orbiter Camera, the Lunar Orbital Laser Altimeter, the Mini-RF instrument, the Chandrayaan-2 Orbital High Resolution Camera, the Korean Pathfinder Lunar Orbiter’s Shadowcam, the Kaguya Spectral Profiler and Multiband Imager, and the Chandrayaan-1 Moon Mineralogy Mapper. We used a variety of techniques to build these maps, including stereogrammetry, shape-from-shading, ice stability depth and surface temperature calculations, and the horizon method for solar illumination and direct-to-Earth communications maps. Altogether, these maps allowed us to survey for boulders, evaluate features in permanently shadowed regions that VIPER might explore, provide mineralogic context for what VIPER’s instruments may learn, estimate the ages and radar properties of craters in the VIPER mission area, and evaluate the potential for gravity traverses with the rover. These data and techniques provided a rich set of information from which both the VIPER science team and engineering teams were able to draw in order to plan a safe landing and to plan a VIPER surface mission that will be both scientifically valuable and robust from an operational perspective.

The Planetary Science Journal↗

USGS2025

Latitudinal gradients of snow contamination in the Rocky Mountains associated with anthropogenic sources

Seasonal snow is an important source of drinking water and recreation, and for agriculture in the Rocky Mountain region. Monitoring snow-water quality can inform on the effects to the albedo and energy balance of the snowpack, and the sources of natural and anthropogenic aerosol and gases. This study analyzed metals in the seasonal snowpack from water year (WY) 2018 for 49 sites. Calcium, lanthanum, and cerium concentrations support the importance of mineral dust to the southern Rocky Mountains. Mercury (Hg), zinc (Zn), and cadmium (Cd) concentrations showed a similar spatial pattern to mineral dust, whereas antimony (Sb) concentrations were highest in the northern Rocky Mountains. To assess the relative contributions from dust versus anthropogenic contaminant sources, enrichment factors (EF) were calculated, with values above 10 indicating anthropogenic contamination. For Cd, Hg, Sb, and Zn, EF values exceeded 10 at northern sites. These observations were compared to spatial trends of EF values of Hg from WY2009 to WY2018, regional monitoring networks, and back trajectory analyses. The agreement between these datasets revealed temporally consistent contaminant sources and/or transport processes to the northern Rocky Mountains snowpack. Sources include current and historical mining and smelting in the region. Strategies to limit the emissions of these metals to the Northern Rockies could benefit from focusing on remediation of contaminated sites, and continued monitoring and mitigation of active mining and smelting.

Colorado , Idaho, Montana, New Mexico, Utah, Wyomi↗

USGS2024

New U-Pb geochronology and geochemistry of Paleozoic metaigneous rocks from western Yukon and eastern Alaska, cross-border synthesis, and implications for tectonic models

The tectonic evolution of and relation between the Yukon-Tanana terrane and the Lake George assemblage, as well as other associated tectonic assemblages in western Yukon and eastern Alaska, have been debated for decades. The Yukon-Tanana terrane is widely considered to be an allochthonous rifted fragment derived from the Laurentian continental margin, whereas the Lake George assemblage and associated assemblages are currently interpreted to be part of the parautochthonous continental margin of western North America (Laurentia). To address these topics, we present 40 new U-Pb zircon ages and 20 new whole-rock geochemical analyses. We incorporate these data into a new compilation of available geological mapping for a large area that straddles the Alaska-Yukon border, together with 34 previously published U-Pb age determinations and an extensive geochemical database of metaigneous rocks from Late Devonian to Early Mississippian and middle to late Permian assemblages in this area. Magmatism in the Lake George assemblage and related assemblages occurred in two pulses from about 371 to 360 and from about 358 to 347 million years ago (Ma); geochemical discrimination diagrams indicate a large crustal component, possibly indicative of arc magmatism, for felsic metaigneous rocks and a range of tectonic environments for mafic rocks. Magmatism in the Fortymile River and related assemblages, and parts of the Nasina assemblage—all parts of the Yukon-Tanana terrane—are mainly Early Mississippian and span a crystallization age range from about 361 to 343 Ma; geochemical discrimination diagrams for these rocks indicate primarily arc geochemical signatures for both mafic and felsic rocks. Middle to late Permian crystallization ages (about 261–253 Ma) are indicated for felsic metaigneous rocks in the Klondike assemblage and some of the felsic metaigneous rocks in the Nasina assemblage. Based on our mapping, we propose the existence of a possible unconformity between the Mississippian and Permian felsic metavolcanic rocks within the Nasina assemblage that is marked by sporadic occurrences of stretched-pebble conglomerate. Our combined database supports the well-established model of a magmatic arc comprising the Fortymile River and Finlayson assemblages of the rifted Yukon-Tanana terrane continental fragment on which a middle to late Permian arc (Klondike assemblage) was later built. The assemblages of the Yukon-Tanana terrane were subsequently intruded by Late Triassic to Early Jurassic granitoids, presumably during reaccretion of the Yukon-Tanana terrane to the continental margin. Permian and Late Triassic to Early Jurassic intrusions have not been mapped in the now structurally lower plate Lake George assemblage; their absence is one of the lines of evidence that have been used to support the parautochthonous, rather than allochthonous, origin of the Lake George assemblage and related assemblages. Our new data, together with previously published ranges of igneous crystallization ages and geochemical tectonic signatures of the Late Devonian to Early Mississippian magmatic rocks in the Lake George assemblage and associated assemblages and in the Fortymile River, Nasina, and correlated assemblages of the Yukon-Tanana terrane, indicate that the currently accepted interpretation of the Lake George assemblage and associated rocks being part of parauthochthonous North America is not the only possible interpretation of this tectonic entity. Approximately half of the dated intrusive rocks in the Lake George assemblage are contemporaneous with the metaigneous rocks of the Yukon-Tanana terrane arc (<361 Ma). We speculate that our approximately 361 Ma U-Pb age for quartz syenite in part of the North American continental margin in south-central Yukon defines the beginning of rifting of the Laurentian margin. Although the currently favored model of prolonged middle Paleozoic subduction and extension in both the Yukon-Tanana terrane and parautochthonous North America allows for simultaneous middle Paleozoic magmatism on both sides of the Slide Mountain Ocean, we now propose an alternative hypothesis in which the Lake George assemblage represents a deeper part of the rifted Yukon-Tanana terrane arc. If this is the case, the absence of Permian and Late Triassic to Early Jurassic arc rocks in the Lake George assemblage could be explained either by the arcs of these ages not being wide enough to have affected the Lake George assemblage or by tectonic displacement of these arc rocks away from the Lake George assemblage. Our approximately 259 Ma U-Pb zircon age and geochemical analyses of metarhyolite in the Seventymile terrane in Alaska, which comprises remnants of the back-arc basin that separated the Yukon-Tanana terrane from the Laurentian continental margin, confirm the presence of a late middle Permian volcanic arc component to the terrane. Our approximately 319 Ma U-Pb zircon age from the Chicken assemblage (as redefined in this study) in eastern Alaska, combined with previously reported fossil ages and a U-Pb zircon age from this assemblage, indicate that it is a Late Mississippian to Early Pennsylvanian arc assemblage. We propose several other relatively young, locally developed arc assemblages outboard of the ancient continental margin of Laurentia that may correlate with the Chicken assemblage, but we consider its origin to remain an enigma.

Alaska↗

USGS2024

Detrital zircons and the magmatic history of Viti Levu, Fiji

We integrate the existing detrital zircon data from multiple modern river sediment samples on Viti Levu, Fiji, with the most current available geological and topographic mapping of the respective river drainage basins to compare detrital populations with potential bedrock sources. The temporal and spatial variations in zircon geochemistry supplement what is known from igneous rocks and confirm the petrological differences between plutonic and volcanic rocks from the Eocene to early Oligocene (Yavuna age, >30 Ma), middle Oligocene to middle Miocene (Wainimala age, 30–12.5 Ma), late Miocene (Colo age, 12.5–6.5 Ma) and latest Miocene (Namosi age, 6.5–5 Ma). The >30 Ma Yavuna-age zircons are restricted to areas that drain the previously mapped Yavuna Group. The 30–12.5 Ma zircons are found across central Viti Levu from west to east, and the 30–15 Ma zircons have distinctively low U/Yb and high Dy/Yb ratios. They are the best radiometric evidence of widespread early to middle Miocene arc magmatism in Fiji that was relatively U-poor. Peak deconvolution of the Colo age zircons from individual basins suggests the following ages for undated or poorly dated plutons from central Viti Levu. The large Mavuvu pluton is probably composed of multiple intrusions in the 12–10 Ma range, the Waiqa pluton is probably ca 10 Ma, and the Noikoro pluton is probably ca 9 Ma. There are zircons from unknown plutonic or volcanic sources between 8 and 7 Ma in western Viti Levu that have distinct Eu/Eu* ratios. We attribute the highest U/Yb ratios in some Colo age zircons to crustal anatexis. Namosi-age zircons are abundant in the Medrausucu Group and can be found scattered across Viti Levu.

Viti Levu↗

USGS2025

Identification of representative earthquakes for probabilistic tsunami hazard analysis (PTHA) using earthquake rupture forecasts and machine learning

As probabilistic tsunami hazard analysis (PTHA) focuses more on assessments for localized, populous regions, techniques are needed to identify a subsample of representative earthquake ruptures to make the computational requirements for producing high-resolution hazard maps tractable. Moreover, the greatest epistemic uncertainty in seismic PTHA is related to source characterization, which is often poorly defined and subjective. We address these two salient issues by applying streamlined earthquake rupture forecasts (ERFs), based on combinatorial optimization methods, to an unsupervised machine learning workflow for identifying representative ruptures. ERFs determine the optimal distribution of a millennia-scale sample of earthquakes by inverting the observed slip rate on major faults. We use two previously developed combinatorial optimization ERFs, integer programming and greedy sequential, to produce the optimal location of ruptures with seismic moments sampled from a regional Gutenberg–Richter magnitude–frequency distribution. These ruptures in turn are used to calculate peak nearshore tsunami amplitude, using computationally efficient tsunami Green's functions. An unsupervised machine learning workflow is then used to identify a small subsample of the earthquakes input to ERFs for onshore PTHA analysis. We eliminate epistemic uncertainty related to source distribution under traditional PTHA analysis; in its place, a quantifiable, less subjective and generally smaller uncertainty related to the input to ERFs is included. The Nankai subduction zone is used as a test case, where previous ERFs have been conducted. Results indicate that the locations of representative earthquakes are sensitive to choice of magnitude–area relation and to whether a minimum cumulative stress objective is imposed on the fault. In general, incorporating ERFs into PTHA provide a physically self-consistent method to incorporate fault slip information in determining representative earthquakes for onshore PTHA, eliminating a major source of epistemic uncertainty.

Nankai subduction zone↗

USGS2026

Rapid seismic and infrasound assessment of large landslides: A case study from Denali National Park and Preserve (Alaska)

Large, rapid landslides are a global hazard that can occur in remote, mountainous areas. Eyewitness reports of landslides and satellite imagery can often be limited or delayed, particularly during inclement weather. However, landslide-generated seismic and infrasound (low-frequency atmospheric sound) waves can be remotely detected in near real-time. This information can significantly expedite characterization and possible landslide response activities. Here, we highlight these capabilities using a > 4 million m³ ice–rock avalanche in Denali National Park and Preserve (Alaska). This event was detected via a landslide-specific seismic location and volume estimation algorithm deployed in Alaska, and — notably — by standard earthquake monitoring systems. Following rapid detection of this event, we combined its seismic and infrasound dataset with optical, synthetic aperture radar, and oblique aerial imagery, multitemporal digital elevation models, and a numerical flow model to reconstruct its failure timeline and dynamics. We apply array processing to infrasound signals traveling > 250 km and find that two precursory events occurred minutes prior to the main failure. We use long-period seismic signals to infer the force exerted by the landslide on the Earth and constrain the rheological parameters of our numerical flow simulation with this result and deposit morphology. The main failure produced a steeply-dipping impulsive initial downward force and reached speeds exceeding 60 m/s. This impulsive force generated relatively strong seismic body waves, which contributed to the earthquake system detection. This large, remote Alaska landslide underscores the key value of seismic and infrasound analysis for rapid landslide assessment and motivates efforts to further operationalize these approaches.

Alaska↗

USGS2026

Upper Triassic igneous rocks of the southern Kenai Peninsula, Alaska – Prelude to Early Jurassic subduction along the western Wrangellia composite terrane margin

New U–Pb zircon geochronology identifies a latest Triassic (ca 214–201 Ma) igneous suite of tuff, hypabyssal dikes, and a pluton on the southern Kenai Peninsula, Alaska. The igneous suite was emplaced within Upper Triassic sedimentary rocks along the southern margin of Western Wrangellia, the western-most fragment of the Wrangellia composite terrane. The igneous rocks range from mafic (50.6% SiO 2 ) to felsic (78.3% SiO 2 ), characteristically have less than 1.55% K 2 O, and generally have low trace element abundances. The tonalitic and trondhjemitic magmas were largely sourced in mafic-rich lower crust and incompletely assimilated quartz and other mineral xenocrysts are common. Fractionation involving plagioclase and amphibole is indicated for some magmas and composite intrusions and igneous xenoliths indicate magma mixing was possible. Paleozoic and Precambrian inherited zircons and initial 87 Sr/ 86 Sr (0.704103–0.705609) and 143 Nd/ 144 Nd (0.512396–0.512777) ratios indicate that the Western Wrangellia crustal sources are heterogeneous and contain sialic components. The latest Triassic magmatism reflects processes that preceded Early Jurassic subduction along the Wrangellia composite terrane and Pacific Ocean plate boundary. These processes involved heating and melting of mantle lithosphere and lower crust as mantle instabilities accompanied the breaking of the plate boundary linkages. The Late Triassic transition to subduction along the Wrangellia composite terrane margin coincided with the transition to subduction cessation in the Late Triassic arcs of the western Intermontane terranes of Canada. The shift to subduction along the outboard Wrangellia composite terrane margin marks the beginning of the Pacific Ocean–Cordillera plate interactions that came to dominate the tectonic evolution of the northern Cordillera from the Early Jurassic to today.

Alaska↗

USGS2025

Refinements to the attenuated propagation of local earthquake shaking (APPLES) ground-motion-based earthquake early warning algorithm

We refined the Attenuated ProPagation of Local Earthquake Shaking (APPLES) ground-motion-based earthquake early warning (EEW) approach, and directly compare APPLES performance with that of the source-characterization-based U.S. ShakeAlert EEW system for a suite of historical earthquakes in the U.S. West Coast and Japan. APPLES is an extension of the Propagation of Local Undamped Motion (PLUM) algorithm in which observed shaking intensity at seismic stations is used to forward-predict intensity distributions to surrounding areas using an attenuation model derived from an intensity prediction equation. We test new configuration options within APPLES, such as using the second highest estimated ground motion rather than the maximum, to better match median ground-motion observations and reduce alerts for small magnitude earthquakes, both of which are key alerting priorities within ShakeAlert. We evaluate these configurations alongside ShakeAlert by comparing the ground-motion estimation accuracy and available warning times relative to station observations and ShakeMap distributions. Our preferred APPLES configuration produces accurate ground-motion estimates and corresponds better with median observations compared to ShakeAlert’s estimates. This preferred configuration substantially reduces alert issuance for M < 5.0 earthquakes compared to the previous APPLES configuration, and alert-release criteria can further restrict alerts to primarily M ≥ 5.5 earthquakes without requiring magnitude estimation. Prioritizing matching median-observed ground motions may reduce APPLES warning times compared to configurations that were tuned to avoid missed alerts (such as those that use the maximum estimated ground motions), which can lead to shorter warning times compared to ShakeAlert for the same alert threshold. However, station-based warning time assessments demonstrate that APPLES can outperform ShakeAlert for high target thresholds. APPLES is a simple, independent EEW approach that may improve the robustness of EEW for the West Coast of the U.S.

California, Oregon, Washington↗

USGS2026

Conventional hydrothermal power-producing systems of the Great Basin, USA

As part of the update to the electric-grade conventional hydrothermal assessment of the Great Basin, USA, Monte Carlo analyses of identified resources within explored regions will be performed to make estimates of discovered resources and associated uncertainty. Analyses use conditional statistics where estimates are conditioned upon a hydrothermal favorability map, allowing for the likelihood that more resources exist in regions of higher hydrothermal favorability. For these analyses, a dataset of identified hydrothermal systems is compiled, and the new compilation is described herein. Recognizing that a single hydrothermal system may be developed with multiple power plants, and that the hydrothermal upflow zone may be several kilometers across with many measurements characterizing a single hydrothermal system, a procedure was developed and employed to create clusters of points (power plants, measurements, etc.) that are associated with a single system, and a new central point was defined as the best estimator of the center of the hydrothermal system. Hydrothermal systems were uniquely identified by grouping electric-grade hydrothermal measurements and operating power plants within a distance of 10 km. Groups that are >10 km apart are assumed to be different electric-grade hydrothermal systems. While 10 km was used as the threshold, most systems were significantly further apart, and most points within groups were typically within 5 km of each other. A well measurement was considered an electric-grade measurement of a hydrothermal system if it had two properties: a measured temperature of >85 °C and evidence of hydrothermal convection. Other points that were added to the dataset are locations of operating powerplants or locations that have been classified as an electric-grade hydrothermal resource by either the U.S Geological Survey (USGS) or the Great Basin Center for Geothermal Energy. After all points are assigned to systems, new points were computed with the goal of identifying the center of the throat of the hydrothermal upflow zone. If operating powerplants exist for a system, then the arithmetic average of all power plant locations is used. Otherwise, if USGS made an estimate, that location is used. In the absence of both powerplants or USGS estimates, the arithmetic average of all electric-grade measurement locations is used. An example is shown of how these newly compiled locations might be ranked for uncertainty analyses, where higher confidence is assumed if measured temperature is higher and there are many supporting measurements indicating an electric-grade resource. In summary, 28 systems have operating power plants, an additional 78 systems are known identified electric-grade hydrothermal resources, and 99 new systems were identified as probable electric-grade systems with varying levels of confidence. These 205 locations are shown as a function of a recent hydrothermal favorability map, conceptually illustrating the conditional statistics that can be used to make estimates of the undiscovered resources of the Great Basin. An accompanying data release provides summaries of developed capacity by system and USGS estimates of likely total capacity and associated uncertainty.

Arizona, California, Idaho, Nevada, Oregon, Utah↗

USGS2026

Discrete element investigation of the influence of shallow soil density on the manifestations of strike-slip surface fault rupture

This study investigates the influence of soil relative density on strike-slip surface fault rupture manifestations using three-dimensional numerical simulations performed with the discrete element method (DEM). The simulations capture the formation of distinctive fault strands within complex flower structures using tens of millions of grains. The tendency for dense soils to localize shear manifests as multiple localized shear bands within wide zones of deformation in strike-slip fault rupture, whereas diffuse shear deformation develops within narrow zones in loose soils. The spatial extents of soil deformation are consistent between simulations having similar relative densities but with different quantities of grains. However, individual shears are more distinguishable in assemblages of finer grains than in assemblages of coarser grains. The simulations show the progressive development of new shears within the bounds of previously developed shears. Shear activity transitions inward as fault activity diminishes along the outermost shears and continues along newly developed shears until a vertically dipping throughgoing shear structure develops that accommodates most of the fault displacement thereon. The throughgoing fault develops at smaller fault displacements in looser soils because the first shear rupture propagates closer to the vertical direction and does not undergo as much inward translation of shear activity, as is observed in denser soils. In all simulations, ground surface uplift develops between nonintersecting active shears, and ground surface subsidence tends to develop where new shears intersect previous shears. The surface traces in these simulations are shown to be consistent with analog models and case histories of surface fault rupture occurring in different shallow subsurface materials. Although computationally costly, these modeling results are valuable for providing a strong numerical supplement to traditional analog models used to represent the mechanics of strike-slip zones in soil, and they provide quantifiable stresses and large-strain deformations throughout the model domain.

Journal of Geotechnical and Geoenvironmental Engin↗

USGS2024

Following the tug of the audience from complex to simplified hazards maps at Cascade Range volcanoes

Volcano-hazard maps are broadly recognized as important tools for forecasting and managing volcanic crises and for disseminating spatial information to authorities and people at risk. As scientists, we might presume that hazards maps can be developed at the time and with the methods of our discretion, yet the co-production of maps with stakeholder groups, who have programmatic needs of their own, can sway the timing, usability, and acceptance of map products. We examine two volcano hazard map-making efforts by staff at the U.S. Geological Survey. During the 1990s and early 2000s scientists developed a series of hazard assessments and maps with detailed zonations for volcanoes in Washington and Oregon. In 2009, the National Park Service expressed the need for simplified versions of the existing hazard maps for a high-profile visitor center exhibit. This request created an opportunity for scientists to rethink the objectives, scope, content, and map representations of hazards. The primary focus of this article is a discussion of processes used by scientists to distill the most critical information within the official parent maps into a series of simplified maps using criteria specified. We contextualize this project with information about development of the parent maps, public response to the simplified hazard maps, the value of user engagement in mapmaking, and with reference to the abundance of guidance available to the next generation of hazard-mapmakers. We argue that simplified versions of maps should be developed in tandem with any hazard maps that contain technical complexities, not as a replacement, but as a mechanism to broaden awareness of hazards. We found that when scientists endeavor to design vivid and easy-to-understand maps, people in many professions find uses for them within their organization’s information products, resulting in extensive distribution.

Washington↗

USGS2026

The petrogenesis of Þingmúli volcano, East Fjords, Iceland

In this work we revisit Þingmúli volcano (Þ = Th), a classic locality known as an example of a complete tholeiitic differentiation. Þingmúli is a ~ 9.5 Ma extinct central volcano located in the East Fjords of Iceland, in which the whole compositional spectrum from basalt to rhyolites have erupted. These volcanic products have been previously considered as petrogenetically related by an ideal fractionation trend, regardless any temporal relationship or volumetric considerations. Here we report new whole-rock geochemistry, mineral chemistry, isotope analyses, estimation of residence times of the different eruptive deposits, and an update of the original petrogenetic model. Our results highlight that an enriched source, likely spinel lherzolites, generated transitional-alkaline basaltic melts after 15–20% of partial melting at depths of 40–45 km. Many of these basaltic melts erupted at various stages of the volcano's history, while others remained longer in the volcanic plumbing system. These evolved by fractional crystallisation into basaltic andesite magmas with a residence time of ~5 years based on the crystal size distribution of the plagioclase population. Isotopic differences between the basalts/basaltic andesites ( 87 Sr/ 86 Sr ~ 0.7034; 143 Nd/ 144 Nd ~ 0.51315) and the erupted rhyolites ( 87 Sr/ 86 Sr ~ 0.7037; 143 Nd/ 144 Nd ~ 0.51304) indicate that the latter are not petrogenetically related to the former. Therefore, instead of a fractional crystallisation mechanism to generate the rhyolites, we propose the partial melting of ignimbrite layers located beneath the volcano. The broad range of trace element concentrations in andesites and dacites and their different isotopic values compared to the basalts strongly suggest that these magmas have been generated by magma mixing between basaltic and rhyolitic melts, similar to modern day Icelandic volcanoes such as Hekla. These results highlight the need to revisit previously studied Icelandic classic localities and reassess their traditionally proposed petrogenetic models.

Þingmúli volcano, East Fjords↗

USGS2024

Framework for mapping liquefaction hazard–Targeted design ground motions

Liquefaction-induced ground failure poses substantial challenges to geotechnical earthquake engineering design. Current approaches for designing against liquefaction hazards, as specified in most seismic provisions, focus on estimating a liquefaction factor of safety ( 𝐹⁢𝑆𝐿 ) and typically characterize earthquake loading using design parameters based on probabilistic or deterministic ground motion levels. Because 𝐹⁢𝑆𝐿 is estimated deterministically, this basis of design neglects considerable uncertainties for estimating liquefaction triggering and its consequences and results in a lack of liquefaction-specific design criteria, particularly as structural design has advanced toward risk-targeted performance objectives. This study presents a framework for developing liquefaction-targeted design criteria based on a minimum acceptable return period of liquefaction, informed by probabilistic liquefaction hazard analysis (PLHA). PLHA quantifies annualized rates of liquefaction by considering contributions from (1) the full ground-motion probability space, and (2) uncertainties in liquefaction triggering using probabilistic models. PLHA is used in this study to characterize the current, effective return periods of 𝐹⁢𝑆𝐿 ( 𝑇𝑅,𝐹⁢𝑆 ) obtained from conventional liquefaction hazard analysis (CLHA) using uniform-hazard ground motions. 𝑇𝑅,𝐹⁢𝑆 is evaluated in a parametric study of nearly 100 sites throughout the conterminous United States. The results indicate large geographic variations in acceptable liquefaction hazard levels, with implied 𝑇𝑅,𝐹⁢𝑆 ranging between approximately 1,000 to 3,000 years. To address these inconsistencies without the computational demands of full PLHA, a framework is proposed for developing a liquefaction-targeted design peak ground acceleration, 𝑃⁢𝐺⁢𝐴𝐿 , for use in liquefaction models that result in consistent liquefaction design levels across all geographic locations. The mapped 𝑃⁢𝐺⁢𝐴𝐿 is shown to be somewhat sensitive to site-specific properties, and adjustment factors are developed and presented. The proposed 𝑃⁢𝐺⁢𝐴𝐿 mapping procedure produces 𝐹⁢𝑆𝐿 estimates that are consistent with those obtained from full PLHA at a target 𝑇𝑅,𝐹⁢𝑆 , providing a promising roadmap to incorporating PLHA concepts into current liquefaction design methods.

Journal of Geotechnical and Geoenvironmental Engin↗

USGS2026

Permafrost thaw controls iron flux from wetlands and sulfide-bearing rocks to Arctic rivers and streams

Recent warming has caused widespread iron mobilization into Arctic waterbodies that degrades ecosystems and threatens natural resources. Yet, understanding where and when iron flux occurs remains limited. Here, we investigate iron loading across regional to local scales in Arctic Alaska using climate, water chemistry, and borehole data together with mapped geology and permafrost presence. We show that both anoxic microbial iron reduction and acid rock drainage from iron-sulfide oxidation mobilize iron. Iron influx is strongly associated with lowland wetlands, sulfide-rich upland bedrock, and near-surface permafrost. Acid rock drainage chemistry correlates very strongly with the depth of seasonal thaw above permafrost from the previous year, indicating a one-year lag. These findings clarify the spatial and temporal dynamics of Arctic river rusting, provide a mechanistic understanding of the phenomenon, and may allow anticipation of its occurrence and assessment of its implications for aquatic ecosystem health and subsistence resources under ongoing climate change.

Alaska↗

USGS2026

Nonstationary demographic state-space models using unreplicated counts for species undergoing environmental stressors

A fundamental task in ecological statistics is to estimate abundance and growth rate distributions from wildlife monitoring data to inform conservation management. Modeling time series of wildlife populations presents a number of challenges from both statistical and ecological perspectives, including discreteness; lack of replication; nonstationarity; and observation, demographic, and other phenomenological processes. Nonstationary dynamics are often exhibited by populations undergoing environmental stressors. Models must account for these characteristics to produce reliable estimates of abundance and trends, yet estimation can be challenging with unreplicated data. We propose nonstationary demographic state-space models using unreplicated counts for populations undergoing environmental stressors. A reduced growth rate model matches the complexity of the unreplicated count data, and a fecundity bound on growth rate distributions allows the separation of processes affecting growth rates like environmental stressors from those affecting abundance external to growth rates like migration. NDSSMs allow for the embedding of nonstationary model components, and we explore the use of changepoints, volatility clustering, and migration processes. We apply the proposed nonstationary models in case studies of herons affected by predator/competitor reestablishment and three bat species affected by a fungal pathogen causing white-nose syndrome. Nonstationary models outperform stationary models and generalized linear mixed effects models according to model scoring and visual inspection of predictions, and provide estimates more consistent with published values. Incorporating migration improves model fit universally, even with approximate one-way immigration, most likely because populations are extirpated, recolonized, and increase multiple-fold over the upper bound set by species fecundity. In addition, estimates of the timing and severity of the environmental stressor differed for models with migration. Including nonstationary and demographic components in a fecundity-bounded growth rate model improves inference and benefits interpretability of hyperparameters. In turn, this adjusts uncertainties in predictions of abundance and growth rates over time, providing the ingredients needed for informed conservation analysis and for directing future monitoring of at-risk species.

Journal of Agricultural, Biological and Environmen↗

USGS2026

Shallow differentiation of primitive arc magmas at the Jurassic Emigrant Gap mafic complex, Sierra Nevada, California

The Emigrant Gap composite pluton exposes ultramafic to silicic intrusive rocks that preserve the chemical evolution of primitive mafic arc magmas and their open-system interactions in the upper crust during mid-Jurassic growth of the Sierra Nevada batholith (California). We present field and petrographic observations and mineral and whole-rock chemistry of the ~35-km 2 ultramafic to dioritic Emigrant Gap mafic complex and an adjacent penecontemporaneous ~90-km 2 granodiorite that together make up the composite pluton. In the Emigrant Gap mafic complex, four roughly central masses of dunite, wehrlite, and olivine clinopyroxenite are surrounded by weakly layered gabbronorite and non-layered diorite. The ultramafic rocks are cumulates formed from near-liquidus minerals of primitive arc magmas that accumulated in steep feeder zones with substantial modification by melt–mush reaction as primitive liquids repeatedly transited the mush-filled conduits. The dominant gabbronoritic rocks are the variably accumulative products of more advanced crystallization–differentiation of arc tholeiitic basalts and basaltic andesites. The adjacent granodiorite intrusion originated separately and preserves field and geochemical evidence for assimilation of metasedimentary rocks. Open-system hybridization between the gabbronoritic mushes and the granodioritic magma produced an intervening body of two-pyroxene diorite. We infer that the ultramafic rocks and gabbronorite of the Emigrant Gap mafic complex crystallized from near-primitive arc basaltic to basaltic andesitic magmas at ~0.15–0.3 GPa, with estimated f O 2 of ≥FMQ +1 and dissolved H 2 O concentrations of only ~0.5–2 wt %. Notably, the Emigrant Gap composite pluton is distinct from other Mesozoic plutons in the Sierra Nevada batholith because of (1) its abundance of mafic and ultramafic rocks that crystallized from relatively primitive mafic melts and (2) the low inferred H 2 O concentrations of its parental magmas, indicated by a near absence of igneous amphibole and by the intermediate rather than calcic compositions of plagioclase. A Jurassic regional extension event probably accounts for the formation of relatively dry primitive arc magmas, as well as for their ascending to the upper crust.

California↗

USGS2025

Preliminary depth to basement modeling at Salton Sea, California

The San Andreas Fault – Imperial Fault (SAF-IF) transtensional step-over zone along the southern margin of the Salton Sea hosts substantial geothermal production and lithium brine resources. Recent volcanism at the Salton Buttes and active seismicity along the SAFIF fault system highlight active tectonic and magmatic processes that pose natural hazards and may impact energy and mineral production. Characterizing the subsurface architecture and extent of concealed alteration associated with this tectono-magmatic system enhances understanding of these active processes, associated hazards, and resources. We have compiled a gravity database, consisting of new and re-processed existing data, from which we have constructed a new isostatic residual gravity anomaly map of the Salton trough. We have used this new gravity dataset together with a compilation of publicly available borehole data to develop new depth to basement inversion models for the region. These depth to basement models help to constrain basin geometries, inform alteration mapping, and reveal variations in basement rocks. Due to the concealed nature of the complex tectonic framework at the Salton trough, it is necessary to utilize geophysical methods for subsurface characterization. These new depth to basement models are a first step toward constructing 2D and 3D geophysical and geologic models of the Imperial Valley and Salton Sea geothermal area. This analysis complements other geophysical initiatives, including magnetotelluric (MT) modeling (Tokmakoff et al., 2024), magnetic mapping (Glen and Earney, 2023, 2024) and potential field modeling, and seismic studies focused on hazard and resource investigations in the Imperial Valley.

California↗

USGS2026

Sedimentological and geochemical characterization of lacustrine deposits of the Babouri-Figuil basin, northern Cameroon: Implications for source rocks distribution and petroleum exploration

The West and Central African Rift System (WCARS) refers to the series of Cretaceous rift basins where commercial hydrocarbon accumulations have been discovered. Some of the WCARS frontier basins are currently being investigated to increase our understanding of these basins in light of new commercial discoveries. The present study was performed in the Babouri-Figuil Basin (BFB), which is genetically related to the WCARS and constitutes an area of interest in terms of petroleum prospecting, where the distribution of petroleum source rocks and potential targets for petroleum exploration across the entire basin is poorly understood. For the current study, an integrated facies analysis along with organic and inorganic geochemical techniques were applied to the basin's Cretaceous deposits with the aim of reconstructing the paleodepositional environment, assessing factors that triggered the input of organic matter, and providing a spatial overview of the organic matter accumulation in the basin based on outcrop samples. An alluvial fan-lacustrine-braided river system is inferred from the facies analysis of the stratigraphic sequence consisting of conglomerate, sandstone, siltstone, limestone, marlstone, and claystone. Bulk analysis of organic matter reveals that black shale and massive claystone are the main prospective petroleum source rocks in the basin. Inorganic geochemical analyses reveal the influence of anoxic conditions, moderate to high primary productivity, and low terrigenous inputs in organic matter enrichment. The formations rich in organic matter are predominantly concentrated in the western and eastern parts of the basin which may represent areas with depressions, characterized by high accommodation space. In terms of the regional context of the WCARS rift basins, typical hydrocarbon exploration in the BFB may target basal-conglomerate, sandstone beds situated directly above or/and interbedded with the Lower Cretaceous source rocks, and the Upper Cretaceous sandstone beds. Basement rocks (granite, granodiorite, and gneisses) and oil shale deposits may represent potential unconventional hydrocarbon exploration. The current integrated study provides an insight that should guide future hydrocarbon exploration campaigns in the basin.

Babouri-Figuil basin↗