Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “New Zealand Journal of Geology and Geophysics”

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

Persistent deep long-period seismicity near the Lassen Volcanic Center

Deep long-period (DLP) earthquakes have been observed at many volcanic settings around the world and linked to the magmatic processes that drive volcanic unrest. At the Lassen Volcanic Center (LVC) of the Cascade arc, limited detection of DLP activity hinders classification of anomalous seismic behavior and its relationship to the LVC magmatic system. This study uses a template matching approach with seismic data from temporary nodal and permanent stations to detect and locate DLP earthquakes near the LVC between 2017 and 2024. Within the DLP catalog of 611 events, a transition occurs from scattered (2017–2020) to oscillatory (2021–2024) occurrence rates. During the oscillatory period, regional earthquakes with large amplitude velocity waveforms observed near the LVC are associated with abrupt changes in DLP occurrence rates. Continued monitoring of DLP activity has the potential to better define the processes that drive volcanic unrest at the LVC in the future.

California↗
USGS2025

A northeast-dipping zone of low frequency earthquakes at the southern edge of Cascadia subduction

Tectonic tremor monitoring occasionally detects events in an anomalous zone in southern Cascadia, 50–100 km west of the main tremor band, near the expected southern edge of the subducting Gorda slab at the Mendocino triple junction. To investigate the geometry and temporal behavior of this tremor, we examine its constituent low-frequency earthquakes (LFEs) by developing 27 stacked LFE waveform templates that we use to detect events from 2018 to 2024. We then relocate LFE sources together with regional seismicity. We find that LFE hypocenters form a northeast-dipping alignment at 22–29 km depth, extending eastward from a zone of micro-earthquakes, ∼15 km south of the southern edge of Gorda slab seismicity. These LFE families exhibit small bursts of activity every few days. Considering the strong world-wide association of tremor and LFEs with high slip-rate, plate-bounding faults, we hypothesize these LFEs may demark the southern edge of Cascadia subduction.

California↗
USGS2026

Dynamic coupling between faulting, rifting and magmatism during 2021-2025 unrest on Reykjanes Peninsula, Iceland

Interactions among faulting, earthquakes, and eruptions are fundamental to plate tectonics and hazard forecasting yet rarely observed along mid-ocean ridges. On Iceland's Reykjanes Peninsula, seismotectonic–volcanic unrest resumed after nearly 800-year hiatus, providing an opportunity to observe these interactions during 2021–2025 activity. By integrating high-resolution seismicity, focal mechanisms, satellite geodesy, surface deformation, and eruption data, we document ∼4 m of total extension accommodated through 14 rifting episodes. The largest, in 2023, involved graben reactivation and diking, with seismic swarms and earthquake faulting that matched the surface ruptures, where strike-slip faulting preceded normal-faulting earthquakes and extension. The accrued extension was released by extension fractures triggered by magma accumulation. Long-term observations show no correlation between erupted magma volume, seismicity, and crustal extension. This highlights dynamic relation between rifting, faulting, and magmatism in transtensional settings and their implications for hazard assessment.

Reykjanes Peninsula↗
USGS2024

Hydroacoustic observations reveal drivers of mixing and salinization of a karst subterranean estuary during intense precipitation

Karst subterranean estuaries within globally ubiquitous carbonate aquifers are coastal groundwater ecosystems that provide an essential water resource for human populations. To understand the drivers of salinization within a coastal aquifer in the Yucatan Peninsula (Mexico), we employed hydroacoustics in flooded caves to observe how oceanic and atmospheric events facilitate mixing between the meteoric lens (fresh-brackish groundwater) and the saline groundwater on tidal and episodic timescales. Precipitation during Tropical Storm Carlotta increased the flow and salinity of the meteoric lens without evidence for vertical mixing across the halocline. We postulate that vertical migration of haloclines in the conduit relative to those within the rock matrix during precipitation creates lateral density gradients that drive mixing, and ultimately creates a brackish layer within the meteoric lens. These results provide a mechanistic explanation for vertical and lateral exchange in a coastal carbonate aquifer, which has implications for groundwater response to future climatic change.

Ox Bel Ha cave system, Yucatan peninsula↗
USGS2025

Strong shaking from past Cascadia Subduction Zone earthquakes encoded in coastal landforms

Strong earthquakes along subduction zones are often devastating events, but sparse records along some tectonic margins limit our understanding of seismic hazards. Constraining shaking intensities is critical, especially in subduction zones with infrequent but large-magnitude earthquakes like the Cascadia Subduction Zone (CSZ), where the lack of recorded ground motions has led to uncertainty in the severity and potential impacts of future earthquakes. Here we fill this observational gap with a novel inventory of quantitative estimates of past shaking intensities from geotechnical modeling of coastal landforms. One hundred fifty-four deep-seated landslides and 65 fragile geologic features constrain minimum and maximum peak ground accelerations, respectively. These estimates are broadly consistent with model predictions of M9 ruptures, suggesting strong shaking of 0.4–0.8 g during past CSZ earthquakes. Local discrepancies between our geologic shaking constraints and earthquake simulations may inform past rupture behavior, leading to better predictions of shaking intensity for future earthquakes.

California, Oregon, Washington↗
USGS2025

Improved prediction of postfire debris flows through rainfall anomaly maps

Predicting where runoff-generated debris flows might occur during rainfall on steep, recently burned terrain is challenging. Studies of mass-movement processes in unburned areas indicate that event locations are well-predicted by rainfall anomaly, R* , in which peak observed rainfall is normalized by local rainfall climatology. Here, we use remote and field methods to map debris flows triggered within the 2020 Dolan Fire burn area in coastal California, demonstrate that a short-duration R* metric predicts debris-flow occurrence more effectively than absolute peak intensity or longer-duration rainfall metrics, and show that incorporating an R* criterion into an existing debris-flow likelihood model can reduce false positive predictions and improve accuracy. We test R * at three other climatically distinct fires in California, demonstrating its utility for mapping likely debris-flow locations in different climates. We also consider how R* can benefit postfire debris-flow prediction given recent increases in climatological variability within individual burn perimeters.

Callifornia↗
USGS1976

Solute transport and modeling of water quality in a small stream

An injection of chloride, sodium, and stable strontium was made at a constant rate for 3 hours into Uvas Creek, Santa Clara County, Calif., to determine the mass transport processes in a small stream. Five observation points were selected within a 610-metre reach of the stream below the injection site. Water samples were collected at the observation points during and immediately after the injection. A mathematical model of the stream was obtained by solving analytically and optimally the one-dimensional mass transport equation of the solutes in the stream. Comparison of field results with a simplified mathematical model indicates the dominance of convection in the behavior of sodium and chloride. The concentration of chloride and sodium can be closely simulated by the model. However, strontium cannot be well represented by the simplified model, which contains a first-order decay-type sink.

California↗
USGS2026

Remote compositional analyses of space-weathered lunar maria

Visible-to-shortwave infrared (VSWIR) reflectance spectroscopy has revolutionized our understanding of planetary surface compositions. However, space-weathering processes on airless bodies complicate quantitative compositional analyses. Here, we present a framework to isolate the signatures of space weathering in VSWIR spectra of lunar maria by leveraging radiative transfer modeling under the assumptions that (i) a space-weathered target can be expressed as a mixture of fresh and fully space-weathered components and (ii) remaining signatures can be modeled by including agglutinates as an end-member component. We first validate this approach against laboratory spectra of space-weathered Apollo mare soils of known mineral compositions using a probabilistic Markov Chain Monte Carlo implementation of the Hapke radiative transfer model. Second, we illustrate how this approach can be applied to orbital Moon Mineralogy Mapper data. The proposed space-weathering correction workflow for lunar maria could be expanded to other lunar lithologies and applied to existing and future data sets.

Planetary Science Journal↗
USGS2025

Ubiquity and causes of soil water preferential flow across 17 ecoregions

Preferential flow (PF) in soil causes the rapid transport of water, nutrients, and contaminants into the subsurface, influencing groundwater recharge and streamflow. Data scarcity has hindered the quantification of PF occurrence and the identification of its drivers across diverse ecoregions. We address this gap by analyzing high-frequency, multi-depth soil moisture data across 17 ecoregions in the USA, using ∼1,500 sensors at 40 sites. We discovered that PF is widespread, with sites experiencing PF in up to 60% of rainfall events ≥2 mm. Multiple approaches consistently show that PF is more likely to occur with increased peak rainfall intensity, finer textured material, low soil moisture variability, humid climate, and higher net primary productivity. This suggests that PF patterns could shift with projected climate changes, increasing uncertainty in predictions of groundwater recharge, water quality, and streamflow generation.

conterminous United States and Puerto Rico↗
USGS2025

Bølling-Allerød productivity in the subarctic Pacific driven by seasonal upwelling

The Bølling-Allerød deglacial event is marked by high diatom productivity and opal deposition throughout the subarctic Pacific. This opal could either constitute a strengthened biological pump and thus carbon sequestration, or a weakened biological pump and release of marine-sequestered CO 2 to the atmosphere. We quantify silicic acid supply at IODP Site U1340 in the Bering Sea using biogenic opal and δ 30 Si of Coscinodiscus , a diatom genus. These records, along with diatom environmental indicators, suggest the Bølling-Allerød had high silicic acid availability related to a shift from stratification to seasonal upwelling dynamics. We thus propose the primary cause of the high productivity event was increased macronutrient supply from vertical exchange that injected old, nutrient-rich, CO 2 -rich waters into the surface. Enhanced CO 2 release from the subarctic Pacific may help explain critical intervals of CO 2 rise that occur at the onsets of the Bølling-Allerød and PreBoreal.

Geophysical Research Letters↗
USGS2025

Mapping a Carrington storm

A map is presented of median 1-min-resolution peak geoelectric-field strength across the United States as would be induced by magnetic storms as intense as the 2 September 1859 Carrington storm. The map is constructed from two data sets: Magnetometer time series from 22 ground-based observatories recording 40 magnetic storms, and surface impedance tensors derived from magnetotelluric measurements acquired at 1616 survey sites across the contiguous United States. Carrington-class storm geoelectric fields are likely to be very strong in the United States East and Midwest; > 5.00 V/km at many places. In Virginia, strengths would likely range from 30.30 V/km, with a 68% confidence interval of [19.44,47.20] V/km, to as low as 0.05 [0.03,0.07] V/km. Comparison of model geopotentials with those measured on 30 long lines, indicates errors of about 18%. A Carrington-class storm would likely induce geoelectric fields with strengths 55% greater than for the 13–14 March 1989 storm.

contiguous United States↗
USGS2025

Apatite geo-thermochronology and geochemistry constrain Oligocene-Miocene growth and geodynamics of the northeastern Tibetan Plateau

Understanding the geodynamics of plateau evolution requires examining the spatial and temporal aspects of mountain building in the northeastern Tibetan Plateau, which are still under debate. Here we integrate apatite geo-thermochronological and geochemical data from the Oligocene-Miocene succession of the Xunhua Basin to elucidate the evolution of the regional topography. The results suggest sediment provenance changes, at ca. 28, 20 and 12 Ma, indicating topographic growth of the West Qinling, Laji Shan, and Jishi Shan source areas, respectively. Our findings reveal Oligocene-Miocene stepwise deformation and middle Miocene stress reorganization within the northeastern Tibetan Plateau. We attribute this to Oligocene lithospheric foundering beneath the south-central Tibetan Plateau and the middle Miocene initiation of sinistral strike-slip faults due to the underthrusting of India and North China. This study highlights the roles of lithospheric removal and continental underthrusting in plateau growth, and the value of multi-proxy apatite analyses in provenance studies.

Northeastern Tibetan Plateau↗
USGS2025

The impact of 3D structure on coseismic coastal land-level change and tsunami generation in the Cascadia Subduction Zone

Estimates of coseismic vertical displacements from past and potential future subduction zone earthquakes provide critical constraints on regional seismic and tsunami hazard. Many studies use elastic homogeneous half-space models to calculate vertical displacements from a specified earthquake slip distribution, neglecting complexities of 3D structure in subduction zone settings. Here, we use 3D ground motion simulations of potential Cascadia Subduction Zone earthquakes to investigate the impact of realistic 3D Earth structure on estimated vertical displacements. We show that offshore uplift increases when including 3D structure due to high fault slip in low-rigidity accretionary wedge sediments and coastal subsidence decreases due to higher rigidity material at depth. These larger offshore uplifts cause increased tsunami maximum wave heights, and coastal subsidence at paleoseismic sites decreases on average by ∼17 cm, or ∼60%, which can have important implications for the amplitude and extent of slip in earthquake scenarios developed based on paleoseismic data constraints.

California, Oregon, Washington↗
USGS2025

Slow rupture, long rise times, and multi-fault geometry: The 2020 M6.4 southwestern Puerto Rico mainshock

The M 6.4 mainshock of the southwestern Puerto Rico seismic sequence on 7 January 2020, was one of the most impactful modern earthquakes in the northeastern Caribbean. Due to its offshore location and complex aftershock distribution, its source kinematics remain poorly constrained. This active sequence illuminated a complex set of previously unrecognized structures that indicate multiple causative faults may have slipped during its rupture. Here, we utilize seismic and geodetic observations to enhance model resolution, estimate the finite slip of the mainshock, and test a multi-segment, geologically realistic fault geometry. Our refined model finds a lower rupture velocity and longer rise times than typical for an event of this magnitude. This indicates a slow-evolving rupture process that resembles characteristics of a tsunami earthquake. Although this normal/strike-slip faulting event was not tsunamigenic, these qualities, if pervasive for this region, have important implications for future seismic monitoring and hazards in southwestern Puerto Rico.

Puerto Rico↗
USGS2026

Late Miocene Colorado River arrival in the Bidahochi basin supports spillover origin of Grand Canyon

The timing and mechanism of the integration of the Colorado River and incision of the Grand Canyon remain among geology’s enduring controversies. A key question is the configuration of the upper Colorado River watershed between 11 and 6 million years ago. In this study, we present new evidence from zircon uranium-lead geochronology for the arrival of distinctive Colorado–Green River sediment in the Bidahochi basin by 6.6 million years ago derived from the Browns Park Formation. This is coeval with an order-of-magnitude increase in depositional rate, an increase in carbonate strontium isotope ( 87 Sr/ 86 Sr) ratios, the appearance of large fish species characteristic of fast-flowing waters, and other sedimentological changes. This evidence is consistent with the Colorado River supplying water and sediment to the Bidahochi basin before spillover integration of the river through the Grand Canyon.

Arizona, Nevada, Utah↗
USGS2025

Prospectivity modeling of the NASA VIPER landing site at Mons Mouton near the Lunar South Pole

We use a high-resolution digital elevation model and a numerical thermal model to produce a variety of inputs for a water-ice prospectivity model for the Volatiles Investigating Polar Exploration Rover (VIPER) landing site. These input data are maps of topography, surface slope, surface aspect, surface curvature, maximum temperature, depth to ice stability, permanently shadowed regions (PSRs), distance to PSRs, and PSR density. This model predicts where water ice is most likely within the top meter of regolith, assuming plausible relationships between ice concentration and the various inputs. The model is designed to be adjusted in near-real time as data are collected during the VIPER mission. As such, it is a tool for both analyzing data from the mission as well as planning operations. Since the current model, at this point, relies only on orbital remote sensing, the final version will also be a tool to extrapolate the VIPER mission results across the lunar poles.

Planetary Science Journal↗
USGS2024

Evaluating the impact of uncertainty in ground motion forecasts for post-earthquake impact modeling applications

The US Geological Survey’s (USGS) ShakeMap system provides a rapid characterization of strong ground shaking in areas directly affected by an earthquake. This study focuses on studying the aggregate effects of macroseismic shaking estimates from ShakeMap, expressed in terms of modified Mercalli intensity (MMI), when accounting for the uncertainty in forecasted ground motions. We use a Monte Carlo approach to generate numerous spatially correlated realizations of ground motions by utilizing a combination of circulant embedding and kriging techniques for efficiently handling the correlations. We then assessed the aggregate effects of shaking by looking at bin counts across these realizations. We demonstrate that the aggregate shaking regarding the mean macroseismic intensity estimates (from the ShakeMap output) is a biased representation of the aggregate shaking when shaking uncertainty is included. Incorporating shaking uncertainty can help to improve various downstream earthquake impact applications, such as the USGS Prompt Assessment of Global Earthquakes for Response (PAGER) overall earthquake fatality distribution or estimates of shaking-induced ground failure impacts from consequential earthquakes.

Earthquake Spectra Journal↗
USGS2026

Deformed submarine terraces in Puget Sound, Pacific Northwest, indicate only one M >~7.5 earthquake on the Seattle fault zone in the past 11,000 yr

Submerged marine terraces in Puget Sound, deformed across the Seattle fault zone (SFZ), indicate that only one earthquake as large as M~7.5 has occurred in at least the past 11 kyr. Previous paleoseismic studies document a M~7.5 earthquake between 923–4 CE, which uplifted coastal marine terraces by as much as 8 m. We demonstrate that this earthquake was the only such event since ~11 ka by mapping and quantifying deformation of older marine terraces, now submerged in Puget Sound. The submerged terraces, attributed to a late-glacial sea-level lowstand, record both glacial isostatic rebound and tectonic deformation. Vertical offset of the ~11 ka terraces within the SFZ is comparable to that of the marine terraces uplifted in 923 CE, implying no additional large (M>~7.5) earthquake on the SFZ since ~11 ka. This result implies a longer recurrence interval than current hazard estimates, which assumes recurrence of M>7.1 events every 5 kyr. Our mapping of SFZ deformation since ~11 ka also supports fault segmentation and contiguous block uplift between the Seattle and Tacoma fault zones.

Washington↗