USGS Science⌕ Search

SEARCH · USGS Science

Results for “Marine Ornithology: Journal of Seabird Research and Conservation”

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

A glimpse into the future of tectonic tremor monitoring

Tectonic tremor is a weak, long-duration seismic signal often observed in subduction zones and on some other plate-bounding faults. Because of tremor's characteristically low amplitude (and low signal-to-noise) and lack of clear phase arrivals, detecting and locating tremor usually requires techniques distinct from those applied to typical earthquakes. Major advances in detection and understanding of tremor have derived in the past from a powerful combination of new data and new analysis techniques. In a recent study, Sagae et al. (2025, https://doi.org/10.1029/2025jb031348 ) exploit that combination again, developing a new machine-learning based workflow and applying it to the S-net cabled seismic network in the Japan trench offshore northern Honshu. Their approach, although complex, succeeds in detecting several times more tremor activity than earlier studies, resulting in new insights and providing a blueprint for similar approaches that could be applied elsewhere. As real-time earthquake monitoring adopts similar tools, it may present an opportunity to bring tremor monitoring into operational workflows. In turn, this could solidify tremor monitoring as a component of future operational earthquake forecasting.

JGR Solid Earth↗
USGS2025

Characterization and validation of tidally calibrated strains from the Alto Tiberina Near Fault Observatory Strainmeter Array (TABOO-NFO-STAR)

Six horizontal borehole tensor strainmeters (TSM1-6) installed from Fall 2021 to Spring 2022 comprise the Alto Tiberina Near Fault Observatory Strainmeter Array (STAR), providing an unprecedented opportunity to investigate seismic and aseismic deformation from hazardous high- and low-angle normal faults in Italy. Prior to use in tectonic applications, they require in-situ calibration and correction for non-tectonic signals. We tidally calibrate the instruments, characterize the calibration uncertainty, and test the results against environmental and earthquake signals originating from local to teleseismic distances. The STAR sites demonstrably deviate from assumptions common to the standard manufacturer's calibrations, including negative areal coupling at TSM3-6. While the tidally calibrated strains have ~3-56% uncertainty, the calibrated dynamic strains show interstation precision and accuracy to nanostrain levels, and static coseismic offsets in the array footprint are within uncertainty. TSM3 records a complex series of strains that may arise from dynamically triggered near-borehole fracture slip and fluid flow that does not appear to affect its sensitivity to lower strain rate deformation. Future calibration improvement may be afforded with longer stable timeseries, particularly for TSM4. Overall, our analyses demonstrate expanded geodetic capability for detecting deformation in the Alto Tiberina Near Fault Observatory.

Alto Tiberina Near Fault Observatory↗
USGS2025

Framework for assessing intranational risk of critical mineral sectors: Case study of India

Recent geopolitical actions regarding the supply of critical minerals have highlighted both the threats posed by countries dominating a critical mineral supply chain, and the potential mitigations served by alternative supplier countries. However, no systematic quantification exists for the intranational factors that determine such threats or mitigations. This study outlines a methodological framework for quantifying the intranational risk profile of critical mineral sectors and compares them with the sector's international potential for leverage. The four intranational factors of concentration, synchronization, private sector involvement, and regulatory obstacles were analyzed across multiple scales in order to determine an intranational risk score that quantifies the resilience or fragility of the critical mineral sector in question. Network analysis was also separately conducted to identify the most important entities operating across all critical mineral sectors. Finally, this study assessed the critical mineral sector's global share of reserves and production to determine an international risk score that quantifies the country's critical mineral sector's leverage potential. A case study of India showed that two critical mineral sectors with high leverage potential, chromite and barite, were respectively intranationally resilient and fragile. Of the remaining sectors with minimal leverage, five sectors such as aluminum were intranationally resilient and may serve as viable mitigation options for future critical mineral supply. Three sectors such as titanium were intranationally fragile, but recent regulatory developments may both decrease intranational risk and increase international risk, potentially making them major critical mineral sectors in the future.

Resources Policy↗
USGS2025

Structural analysis of brittle-plastic shear zones in the Sangre de Cristo Range, southern Colorado USA: Superposition of Rio Grande rift extension on Laramide contraction

The Sangre de Cristo Range in southern Colorado exposes some of the deepest Cenozoic structural levels in the Rocky Mountain region, including mylonitic shear zones associated with both the Laramide orogeny and Rio Grande rift. We investigated the relation between Laramide contraction and Rio Grande rift extension with detailed geologic mapping, kinematic analysis, and geochronometry in a 50 km 2 area centered on the Independence Mine shear zone (IMSZ). The 15−100-m-thick IMSZ is one of several shallowly to moderately (∼45° ± 20°) W-SW−dipping brittle-plastic shear zones along the western flank of the range. These shear zones display microstructural evidence of initiation as top-NE contractional mylonite zones, consistent with regional Laramide kinematics, which have been pervasively overprinted by shear fabrics indicating top-SW extensional reactivation. Both top-NE and top-SW shear fabrics involve cataclasis and quartz dislocation creep, although top-SW shear is more commonly localized along phyllosilicate-lined shear bands. Shear zones are hosted predominately within Proterozoic gneiss, and contain abundant chlorite and white mica derived from alteration of hornblende and feldspar, which indicates that weakening driven by fluid reactions played an important role in localizing strain. Extensional overprinting appears to be most pervasive along more steeply dipping portions of shear zones and where secondary phyllosilicates form an interconnected weak phase, which suggests that reactivation was primarily controlled by geometry and rheological contrasts inherited from contraction. One top-SW shear zone adjacent to the IMSZ cuts a late Oligocene gabbro stock, and monazite grains synkinematic with top-SW shear in the IMSZ yielded late Oligocene to Early Miocene U-Th-Pb dates that correspond with initiation of the Rio Grande rift. Reactivation of weak reverse faults may represent an important structural control during initial extension in the middle crust, prior to slip along the high-angle Sangre de Cristo normal fault system.

Colorado↗
USGS2026

Evidence for a biological origin of uranium-rich carbon masses within the Ediacaran Salt Range Formation of Pakistan

Thucholites are unique organic structures found in igneous and sedimentary rocks composed of a U-C-rich interior enclosed by an organic outer shell. Their formation and occurrence have perplexed scientists for over 100 years. Typically, thucholites are sparse in sedimentary rocks but where found in abundance, they may be the result of rapid paleoecological disruptions, e.g., volcanic ashfall. Here, we evaluated thucholites from the Ediacaran Salt Range Formation of the Indus Basin in Pakistan using field emission scanning electron microscopy (FESEM) and nanoscale secondary ion mass spectroscopy (NanoSIMS) to propose a two-stage mechanism for thucholite development. NanoSIMS results suggest organic matter in thucholite cores formed by biological fractionation (represented by the presence of lighter 13 C/ 12 C and 34 S/ 32 S) while the outer organic mantle formed via radiolysis-induced polymerization. FESEM elemental analysis confirms compositional differences between the two thucholite components (core and mantle), further implying their contrasting origins. Dimensional comparison of thucholite cores and their U-bearing mineral morphologies to ancient and modern U-biomineralized microbes suggests that the thucholite cores formed from metal-biological interactions. The presence of volcanogenic biotite, alkali feldspar, and clay spherules (interpreted as devitrified volcanic glass) with thucholite suggests that the thucholite cores in the Salt Range Formation are biological responses to rapid paleoenvironmental change from volcanic eruption which preserved the fossilized morphologies of Proterozoic microorganisms.

Salt Range↗
USGS2026

The potential impact of three-dimensional distributed slip models derived from real-time GNSS data on the performance of the ShakeAlert earthquake early warning system for slab interface earthquakes

The ShakeAlert® earthquake early warning (EEW) system is designed to warn users of imminent strong ground motion with sufficient time to take protective actions. ShakeAlert currently uses three algorithms to characterize the earthquake source. One estimates the location and magnitude using the first few seconds of the P wave and, while fast, tends to underestimate magnitude for M w 7.0+ earthquakes. A second estimates the location, orientation, length, and corresponding magnitude of a line source using observed peak ground acceleration and contributes primarily to M w 5.5+ earthquakes. The third infers earthquake magnitude from peak ground displacement measured using Global Navigation Satellite System (GNSS) data and offers nonsaturating magnitudes for M w 7.0+ earthquakes. Other EEW algorithms exist that infer temporally evolving spatially variable slip on a 3D fault surface from real‐time GNSS data, information that might enable more accurate and timely alerts in the event of large‐magnitude subduction interface earthquakes. Here, we evaluate the potential contribution of one such algorithm, BEFORES ( Minson et al. , 2014 ), to improve ShakeAlert performance through a simulated real‐time implementation of Bayesian evidence‐based fault orientation and real‐time earthquake slip (BEFORES) and other ShakeAlert algorithms using data for eight M w 7.6+ earthquakes. The test results demonstrate that BEFORES can produce well‐constrained and accurate magnitude estimates as soon as or sooner than other EEW algorithms, in turn enabling it to increase the amount of warning time users receive in many cases. However, with a modified Mercalli intensity (MMI) threshold of 3.5, which is commonly used for issuing alerts, BEFORES would tend to alert large geographic regions that did not feel strong shaking (MMI 6+). This effect can be mitigated using a higher alert threshold of MMI 4.5 without negative impact on the amount of warning time obtainable with BEFORES.

Bulletin of the Seismological Society of America↗
USGS2026

Latest Pleistocene to 19th-century earthquakes on bending-moment reverse faults of the Seattle fault zone, Washington

Fault-related folds and their associated secondary faults play a critical yet often underrecognized role in accommodating strain and generating earthquakes in active fold-and-thrust belts. In the Seattle fault zone (SFZ), Washington, USA, we present new paleoseismic, geomorphic, and geophysical evidence for late Pleistocene and Holocene earthquakes on shallow, south-dipping secondary faults—the Lytle Beach and Vasa Park faults—that lie within the hanging wall of the greater SFZ and are on trend with the primary, blind Blakely Harbor fault. Our data show that these structures have ruptured independently, producing localized uplift and deformation at the surface, with the most recent event (RH2) likely occurring in the early nineteenth century. While a temporal overlap between the late Pleistocene RH1 and VP1 earthquakes raises the possibility of a ≥35 km rupture along the Blakely Harbor fault, structural and temporal evidence instead supports independent rupture on individual faults related to folding. We interpret these faults as bending-moment reverse faults that formed within a synclinal hinge zone of the main fault, reflecting mechanical and kinematic influences of the broader fault system. Combined with prior studies, our findings indicate that faulting related to folding dominates the mode of strain release within the SFZ since the late Pleistocene with more frequent earthquake recurrence (∼350 yr) over the past ∼2500 yr.

Washington↗
USGS2025

Near-surface material and topography generate anomalous high-frequency ground motion amplification in Chugiak, Alaska

An ∼3 km long nodal array oriented approximately east–west was deployed in Chugiak, Alaska, by the U.S. Geological Survey during 2021. The array intersects with the permanent NetQuakes station NP.ARTY, where peak ground acceleration (PGA) value of 1.98 g was recorded during the 2018 M w 7.1 Anchorage, Alaska, earthquake, in sharp contrast to the PGA of ∼0.3 g at a site just 4 km to the west. Seismic data for M w 1.8–4.3 aftershocks from the M w 7.1 event recorded by the nodal array confirm the anomalously large ground motions obtained at NP.ARTY as well as similar amplifications at nodes within ∼1 km to the east. Here, we performed 0–10 Hz 3D finite‐difference simulations, including high‐resolution surface topography, to explore the cause of the unexpectedly large amplification. As expected, the simulations computed with a regional 3D tomography velocity model severely underpredict the 0–10 Hz acceleration records at almost all sites. Adding a near‐surface low‐velocity taper to 300 m depth amplifies the accelerations by up to a factor of 5 and enables a reasonable match between the nodal data and simulations at sites to the west of NP.ARTY. However, this model still underpredicts the spectral energy in the area covered by glacial sediments by up to an order of magnitude. The addition of a till layer using a depth‐dependent shear‐wave velocity ( ⁠⁠ V s ) profile along with a homogeneous, 8 m thick low‐velocity layer with V s = 250 m/s representing the kame terraces improves the fit to data to within a factor of 2 at nodes located on top of the glacial sediments. Our study shows that the anomalously large high‐frequency amplification recorded at and near NP.ARTY can be explained by a combination of topographic effects and near‐surface low‐velocity material with amplification effects on the high‐frequency ground motion by up to about 40% and an order of magnitude, respectively.

Alaska↗
USGS2026

Revisiting the utility of regional-scale, high-quality geophysical data in mineral exploration - A case study featuring the Mammoth Magnetic Anomaly, Pinal County, Arizona

Regional aeromagnetic surveys passively measure the total magnetic intensity (TMI) and are a foundational tool used in mineral exploration (Airo, 2015). With the increased global demand and the number of critical mineral resources required for manufacturing high-tech devices, developing high-quality, regional-scale geophysical surveys could aid critical mineral exploration efforts and geologic mapping. In 2019, the U. S. Geological Survey launched the Earth Mapping Resources Initiative (Earth MRI) to modernize the geologic and geophysical mapping of regions that have the potential to contain critical mineral resources within the United States. In support of planning Earth MRI geophysical surveys, Drenth and Grauch (2019) defined five aeromagnetic data quality rankings (rank 1 through rank 5) applying them to the airborne geophysical survey inventory of the United States (Johnson et al., 2021). Rank 1 aeromagnetic surveys are of the highest quality, meeting modern standards and allowing best practices for qualitative and quantitative interpretation; whereas rank 5 aeromagnetic surveys are of the lowest quality, being useful only for qualitative interpretation of broad features. Through the Earth MRI effort, 48 high-quality, regional-scale rank 1 and 2 airborne magnetic and radiometric geophysical surveys have been planned, collected, or publicly release through May 2025 (U. S. Geological Survey, 2025). Here, a portion of a rank 1 Earth MRI aeromagnetic survey in southeast Arizona is presented and compared to a legacy rank 5 aeromagnetic survey over the Mammoth Magnetic Anomaly (MMA), demonstrating how modern, high-quality aeromagnetic data improves our view of crustal geology, aiding mineral exploration.

Arizona↗
USGS2025

Scoping decision-maker needs and science availability to support regional natural capital accounting in the U.S. Colorado River Basin

Natural capital accounting has the potential to yield important policy insights at multiple scales, but there remains a disconnect between regional-scale natural capital accounts and their use for informing policy. In this paper, we propose a roadmap that could lead to the creation of policy-relevant regional accounts, with steps split across an initial scoping phase and a subsequent development phase. We demonstrate the scoping steps in action with an application to the Colorado River Basin (“Basin”), a large watershed in the southwestern United States (U.S.) that has faced aridification and substantial high-profile tradeoffs around the use of its water and other natural resources. Drawing on prior U.S. Geological Survey science co-production efforts, we conducted a series of eight discussion sessions with 41 scientists and science representatives whose work is relevant to Basin water, riparian and riverine ecosystems, upland ecosystems and energy and minerals. We summarise participants' thoughts on key topics and economic linkages, their insights and questions of interest and their recommendations on existing scientific data sources and gaps. We evaluate the suitability of the available data for construction of System of Environmental-Economic Accounting (SEEA) Central Framework and SEEA Ecosystem Accounting accounts, including those for land, water, forests, energy and minerals and ecosystems (covering extent, condition and ecosystem services). We present a series of lessons learned during the scoping phase, as well as lessons that could be relevant for future practitioners engaging in the development phase. The information can help guide the development of timely and relevant regional-scale environmental-economic accounts in the U.S. and beyond.

Arizona, California, Colorado, Nevada, New Mexico,↗
USGS2024

Federal lands greenhouse gas emissions and sequestration in the United States: Estimates for 2005–22

In 2016, the Secretary of the U.S. Department of the Interior requested that the U.S. Geological Survey (USGS) produce a publicly available and annually updated database of estimated greenhouse gas emissions associated with the extraction and use of fossil fuels from Federal lands. The first report in this series included emissions estimates from 2005 to 2014 and were reported for 29 States and two offshore areas. Native American and Tribal lands were not included in that analysis. This report recalculates those previous years (2005–14) with updated data and methods and extends the estimates to 2022. Nationwide emissions from fossil fuels produced on Federal lands in 2022 were 1,081.2 million metric tons of carbon dioxide equivalent (MMT CO 2 Eq.) for CO 2 , 33.4 MMT CO 2 Eq. for methane (CH 4 ), and 4.3 MMT CO 2 Eq. for nitrous oxide (N 2 O). Compared to 2005, the 2022 totals represent decreases in emissions for all three greenhouse gases (by 17 percent for CO 2 , 37 percent for CH 4 , and 30 percent for N 2 O). Emissions from fossil fuels produced on Federal lands represent, on average, 21.8 percent of U.S. emissions for CO 2 , 6.1 percent for CH 4 , and 1.3 percent for N 2 O over the 18 years included in this estimate. The trends and relative magnitudes of the greenhouse gas emissions estimated are roughly parallel to the Federal lands fossil fuel production volumes. In 2021, Federal lands of the conterminous United States stored 70,424 MMT CO 2 Eq. in terrestrial ecosystems. Soils stored most of the terrestrial ecosystem carbon (66 percent), followed by live vegetation (25 percent), deadwood (5 percent), and litter (4 percent). Carbon sequestration on Federal lands was highly variable over time, owing primarily to interannual variability in climate and weather, and variability in land use and land cover (LULC) change and disturbances, among these are wildfires and logging. Between 2005 and 2021, Federal lands sequestered an average of 83 MMT CO 2 Eq./yr. By subtracting the cumulative effects of LULC and disturbance-related CO 2 losses to the atmosphere from the total, we estimate that ecosystems at the national level sequestered CO 2 at an annual mean rate of 17 MMT CO 2 Eq./yr in a term called the net ecosystem exchange (NEE). This annual NEE sequestration value represents about 1.4 percent of average fossil fuel emissions over the same period. The USGS estimates presented in this report represent an accounting for the emissions resulting from fossil fuel extraction on Federal lands and the end-use combustion of those fuels, as well as for the sequestration of carbon in terrestrial ecosystems on Federal lands. A combined net CO 2 emissions estimate, which is the difference between the emitted and sequestered CO 2 from both the fossil fuel and ecosystems estimates, provides context for evaluating the greenhouse gas contributions of activities on these lands. The estimates included in this report can provide context for future energy decisions, as well as a basis to track change in the future.

Scientific Investigations Report↗
USGS2026

Climate-adaptive urban planning: Quantitative assessment of drought impacts and practical strategies for climate-resilient urban green spaces

Urban green spaces (UGSs) are vital for enhancing a city’s resilience and livability; however, their functionality is increasingly jeopardized by drought, particularly in water-scarce regions. This study evaluates drought impact on UGSs in Metropolitan Adelaide, Australia, a representative semi-arid urban system, using satellite-derived Normalized Difference Vegetation Index (NDVI) time-series data spanning 2000–2020. Vegetation dynamics were analyzed through Seasonal-Trend decomposition using Loess (STL), standardized anomaly assessment, lagged Pearson correlation, Ordinary Least Squares (OLS) regression, and Mann–Kendall trend analysis. To isolate climatically sensitive signals, 29 urban lawn patches were examined separately from mixed urban canopy, given their shallow root systems and direct dependence on surface moisture. NDVI declined by approximately 0.09 units during the Millennium Drought (2001–2009), with summer greenness deficits reaching 24% below the 20-year benchmark. Temperature was the dominant driver of lawn NDVI variability (r = −0.863, R 2 = 74.5%), substantially exceeding the effect of rainfall (r = 0.156, R 2 = 2.4%). El Niño–Southern Oscillation (ENSO) cycles modulated vegetation responses, with La Niña years supporting recovery and El Niño years amplifying decline. Post-drought recovery remained incomplete, with NDVI deficits of 8–20% persisting through 2020; full recovery was observed only in 2017, coinciding with the highest recorded summer rainfall. No significant directional trend was detected over the full study period (Mann–Kendall τ = 0.005, p = 0.908). These findings demonstrate that heat, rather than water limitation alone, is the primary driver of vegetation stress in urban systems, highlighting the benefits of integrated management strategies that address both warming and moisture deficits to sustain urban green infrastructure under future climate conditions. We introduce the concept of “urban greenery drought,” referring to a form of vegetation stress in managed urban landscapes where greenness is reduced primarily by elevated temperature and atmospheric demand despite water availability.

Adelaide↗
USGS2025

Grand Canyon landslide-dam and paleolake triggered by the Meteor Crater impact at 56 ka

This paper hypothesizes that the Meteor Crater impact in Arizona, USA, 56,000 years ago triggered landslides in Grand Canyon that dammed the Colorado River and formed Nankoweap paleolake. This is compatible with shock and earthquake physics for the impact that infer a M5.4 seismic event, attenuated to an effective magnitude of M3.5 at Grand Canyon. Results that support the hypothesis include radiocarbon dating of driftwood and luminescence dating of associated slack-water lake sediments that are preserved in caves up to 60 m above the modern Colorado River. Radiocarbon ages from two locations, including Stanton’s Cave, date the driftwood as 55.25 ± 2.44 ka (n = 4). Sediments associated with the driftwood gave a luminescence age of 56.00 ± 6.39 ka (n = 2). These six Grand Canyon dates, and three published ages for the Meteor Crater impact, show statistically indistinguishable results that support the hypothesis for a geologically instantaneous series of events with a mean age of 55.60 ± 1.30 ka. This work highlights the value of radiocarbon dating near the limits of the technique, integration of multiple dating methods, and seismic and landslide hazards associated with meteorite impacts in regions of extreme topography like Grand Canyon.

Arizona↗
USGS2025

Formation of the Mount Weld rare earth element deposit, Western Australia: A carbonatite-derived laterite

Carbonatite-hosted rare earth element (REE) deposits are the primary source of the world’s light REEs. The Mount Weld REE deposit in Western Australia is hosted in a lateritic sequence that reflects supergene enrichment of the underlying carbonatite. Water-rock interaction is a key to the formation of this world-class deposit. REE enrichment in the laterite is controlled by the breakdown of primary minerals, the release and transport of REEs, and the formation of secondary minerals. Secondary REE-bearing phosphate minerals are the primary REE-host phases in the laterite ore with monazite as the dominant phase; other REE-bearing phases include rhabdophane, cerianite, churchite, florencite, and crandallite subgroup minerals. Profiles through the laterite show that in the REE-rich zone, apatite and primary calcite and dolomite have broken down such that the loss of Ca and Mg, as well as Si and K, leads to a relative increase in the REEs. Sequestering of REEs in secondary mineral phases formed by groundwater further enhances the REE concentration.

Mount Weld Mine↗
USGS2026

From collision to subduction: Thermal-kinematic inversions constrain plate boundary structure and dip-slip activity in southwestern New Zealand

Constraints on three-dimensional fault system geometry and kinematics are required to understand both long-term tectonics and contemporary seismic hazard. In many cases, these parameters must be inferred from a combination of different data or analyses. Here, we invert low-temperature thermochronological data to quantitatively parameterize both the crustal-scale geometry and late Cenozoic dip-slip activity of the plate boundary system in southwest New Zealand. We find that the modern transpressive plate boundary system has accommodated the bulk of Australian-Pacific Plate convergence since establishment in the middle Miocene-early Pliocene. Shortening has occurred via (a) dextral-reverse slip on the moderate angle, listric central Alpine Fault décollement, and (b) reverse-slip on either the sharply curved, anti-listric northern Puysegur subduction interface or low-angle offshore thrust fault system that splays from the high-angle southern Alpine Fault at depth. We resolve late Cenozoic shortening rates on these structures that are indistinguishable, and interpret these to reflect the transfer of convergent deformation outboard of the predominantly strike-slip southern Alpine Fault. Our work helps elucidate the transition from strain localized oblique continental collision to strain partitioned subduction. We also provide an estimate of the across-strike width of the central-southern Alpine Fault section boundary at seismogenic depths. This is the first quantitative constraint on the geometry of a structural feature inferred to have arrested several major prehistoric earthquakes. This study furthers understanding of fault structure and dip-slip kinematics along the greatest sources of seismic hazard in and offshore of southwest New Zealand.

South Island↗
USGS2026

Rupture into slow-slip fault regime during the 2018 Mw 6.9 Island of Hawaiʻi earthquake is followed by modest postseismic slip

On 4 May 2018, a M w 6.9 earthquake occurred on the south flank of Kīlauea, in the midst of an historic event that included a voluminous eruption from Kīlauea’s lower East Rift zone and caldera collapse at its summit. The earthquake was a consequence of both short‐ and long‐term stress buildup due to magmatic activity associated with the eruption and steady flank motion, respectively, and it revealed features of Kīlauea’s décollement fault that can inform understanding of future earthquake activity. We used geodetic data to determine the distributions of slip during the coseismic and postseismic periods and compared these with areas of known fault slip during past earthquakes and slow‐slip events (SSEs). The 2018 earthquake ruptured into an area of the décollement fault that was active during quasi‐regular SSEs that occurred in the two decades prior to 2018 but that have not been observed since. The coseismic slip model indicates that the amount of motion on the décollement fault was several times greater than what typically occurred during SSEs, suggesting that it may take decades for the fault to rebuild stress to the point at which SSEs will occur again. Postseismic afterslip also occurred in an area of the fault known to experience slow slip; however, unlike at other creeping faults, postseismic afterslip was rapid, being largely over within 2–3 days. The rapid nature and small magnitude of the postseismic afterslip may be due to the lack of a viscoelastic relaxation component, which is possibly a result of the shallow dip of the décollement fault not transferring stress efficiently into the lower crust.

Hawaii↗
USGS2025

Alteration mapping in granitic gneiss using handheld geophysical and geochemical instruments: Implications for iron oxide-apatite and rare earth elements exploration

The Adirondack Mountains of New York, U.S.A. contain iron oxide-apatite (IOA) mineral deposits with variable concentrations of rare earth elements (REE). The IOA mineral deposits are typically hosted in the Lyon Mountain Granite Gneiss and are spatially correlated with extensive Na metasomatism (albitization) of the surrounding country rocks, although some mineral deposits also occur in metagabbro, paragneiss, and anorthosite hosts. The location of albitization zones is key to finding new IOA mineral deposits and better understanding the mineralization processes associated with their genesis. However, the Na alteration zones are generally not visibly identifiable in outcrop or hand specimens because the color and textural changes are subtle and are thus difficult to map. Here, we discuss the results from testing two spectroscopic techniques (handheld, field gamma-ray spectroscopy, and portable X-ray fluorescence (pXRF)) to characterize albitization zones surrounding IOA mineral deposits. The gamma-ray and pXRF spectrometers successfully distinguished the relative intensity of alteration in the Lyon Mountain Granite Gneiss based on the K abundance. The measured K content decreases towards the IOA mineral deposits, and the magnitude and width of the gradient are similar along the strike of each mineral deposit. Elevated Th and U values are present in host rocks adjacent to REE-bearing IOA mineral deposits. The pXRF and gamma-ray spectrometer K values are remarkably consistent with laboratory-based, whole-rock XRF compositional data and, therefore, useful for semi-quantitative analysis. Notably, albitization aureoles are consistent around REE-rich and REE-poor mineral deposits, suggesting that saline fluids are essential for the petrogenesis of IOA mineral deposits but may not be critical to REE mineralization. Ultimately, this study demonstrates the utility of handheld gamma-ray and pXRF spectrometry for identifying otherwise cryptic albitization gradients associated with IOA mineral deposits in granitic gneiss.

New York↗
USGS2026

Efficient physics‐informed ground‐motion simulations with reduced‐order models: CyberShake implications and high‐resolution site terms for southern San Andreas fault earthquakes

Recent advances in Probabilistic Seismic Hazard Analysis (PSHA) leverage physics‐based ground‐motion simulations to estimate seismic hazard, such as the CyberShake project. However, computational costs quickly escalate when performing PSHA for numerous faults or sites and can become prohibitively expensive. To reduce computational demands, CyberShake uses reciprocity and interpolates physics‐informed corrections from simulations conducted at fewer locations, but the accuracy of these interpolations remains poorly quantified. To quantify the interpolation accuracy, we derive high‐resolution, frequency‐dependent site terms for southern California and compare them with interpolated site terms using the CyberShake approach. We accomplish this by performing a set of earthquake point‐source simulations distributed along the nonplanar fault geometry for the southern San Andreas fault (SSAF) extending from Bombay Beach to Lake Hughes. Using SeisSol, we simulate three minutes of viscoelastic seismic wave propagation for these sources and store the horizontal‐component Green’s functions for 480,000 sites. We then use a scientific machine learning approach based on interpolated proper orthogonal decomposition to construct an accurate reduced‐order model of the Green’s functions to efficiently predict effective amplitude spectra (EAS) for finite‐source rupture models of SSAF earthquakes. Using minimum curvature interpolation with tension, as used in CyberShake, we compare the interpolated site terms against our high‐resolution site terms. We identify local discrepancies with EAS differing by up to a factor of approximately three. Furthermore, we identify locations where unexpectedly high or low ground motions are missed when using the interpolated dataset for these earthquakes. We estimate that our approach may be used within CyberShake to reduce the time‐to‐solution by a factor of 336 for the entire earthquake rupture forecast. Our analysis of physics‐based site terms provides more insight into the seismic hazard due to SSAF ruptures and guides future developments by combining high‐performance computing and reduced‐order modeling techniques for PSHA.

California↗