Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Planets”

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.

At least 955 records · Page 53Linked to original sources
USGS2024

Vortex trapping of suspended sand grains over ripples

Coastal hydrodynamics and morphodynamics integrate the effects of small-scale fluid-sediment interactions; yet, these small-scale processes are not well understood. To investigate sediment trapping by turbulent coherent structures or vortices, the transport of coarse sand over ripples was analyzed in a small-oscillatory flow tunnel with phase-separated Particle Image and Tracking Velocimetry. Results from one of the first direct measurements of vortex-trapped sand grains under oscillatory flows are presented. The vortices mobilized sand grains along the ripple slopes just prior to flow reversal and transported the suspended sediment grains. During several flow cycles, some sand grains were temporarily trapped in the vortex, prescribing semi-circular trajectories off-center from the vortex core in quadrants of the vortex that were closest to the ripple slope, as illustrated by Nielsen (1992, https://doi.org/10.1142/1269 ). Comparisons of the horizontal sediment grain velocity with the horizontal fluid velocity yielded a linear relationship with a slope of 0.87. The vertical grain velocities also varied linearly with the vertical fluid velocity with a slope of approximately 1 and an offset of −0.08 m s −1 . The offset is close to the still water settling velocity for coarse sand grains, as hypothesized during vortex trapping. Additionally, estimates of the off-center distance, between the centers of the semi-circular sediment paths and vortex cores, compared well with the ratio of the settling velocity to the radian frequency of the vortex yielding a linear regression slope of 0.99. Improved understanding of vortex trapping effects on sediment dynamics may decrease uncertainty in model predictions of large-scale coastal hydrodynamics and sediment transport.

Journal of Geophysical Research: Earth Surface↗
USGS2023

Comparing sediment microbial communities of Arctic beaver ponds to tundra lakes and streams

In recent decades the habitat of North American beaver ( Castor canadensis ) has expanded from boreal forests into Arctic tundra ecosystems. Beaver ponds in Arctic watersheds are known to alter stream biogeochemistry, which is likely coupled with changes in the activity and composition of microbial communities inhabiting beaver pond sediments. We investigated bacterial, archaeal, and fungal communities in beaver pond sediments along tundra streams in northwestern Alaska (AK), USA and compared them to those of tundra lakes and streams in north-central Alaska that are unimpacted by beavers. β -glucosidase activity assays indicated higher cellulose degradation potential in beaver ponds than in unimpacted streams and lakes within a watershed absent of beavers. Beta diversity analyses showed that dominant lineages of bacteria and archaea in beaver ponds differed from those in tundra lakes and streams, but dominant fungal lineages did not differ between these sample types. Beaver pond sediments displayed lower relative abundances of Crenarchaeota and Euryarchaeota archaea and of bacteria from typically anaerobic taxonomic groups, suggesting differences in rates of fermentative organic matter (OM) breakdown, syntrophy, and methane generation. Beaver ponds also displayed low relative abundances of Chytridiomycota (putative non-symbiotic) fungi and high relative abundances of ectomycorrhizal (plant symbionts) Basidiomycota fungi, suggesting differences in the occurrence of plant and fungi mutualistic interactions. Beaver ponds also featured microbes with taxonomic identities typically associated with the cycling of nitrogen and sulfur compounds in higher relative abundances than tundra lakes and streams. These findings help clarify the microbiological implications of beavers expanding into high latitude regions.

Alaska↗
USGS2023

Aquatic carbon export and dynamics in mountain headwater streams of the western U.S.

Mountain headwater streams actively cycle carbon, receiving it from terrestrial landscapes and exporting it through downstream transport and gas exchange with the atmosphere. Although their importance is now widely recognized, aquatic carbon fluxes in headwater streams remain poorly characterized. In this study, aquatic carbon fluxes were measured in 15 mountain headwater streams and were used in a geostatistical analysis to determine how landscape characteristics influence aquatic carbon fluxes. In-stream sensors were used to measure fluorescent dissolved organic matter (fDOM) (a surrogate for dissolved organic carbon (DOC)) at a subset of sites to characterize dynamic responses to hydroclimatic events. Wetlands have a positive influence on aquatic carbon fluxes, whereas perennial snow/ice has the opposite effect, reflecting differences in soil organic matter content in these landscapes. Mean annual temperature (MAT) has a complex influence on DOC, with peak DOC exports in basins with MAT of 0–2°C. Precipitation has a strong positive influence on aquatic carbon fluxes, and declining snowpacks in the western United States may reduce future aquatic carbon exports. fDOM (and by implication DOC) and

Journal of Geophysical Research: Biogeosciences.↗
USGS2024

Spatial, seasonal, and diel controls of nitrogen-carbon-oxygen cycling during lake-water infiltration to an aquifer

Many freshwater lakes are groundwater flow-through systems. Although lakes commonly are considered to be sinks for nitrogen inputs, relatively little is known about carbon and nitrogen export from lakes to groundwater. The current study focused on lake-bottom biogeochemical processes accompanying the transport of nitrogen, dissolved oxygen (O 2 ), and dissolved organic carbon (DOC) during lake-water recharge from a groundwater flow-through lake. Lake-water and porewater (15–100 cm below lakebed) samples were collected along transects within the lake downwelling zone. Infiltrating porewater O 2 and DOC concentrations decreased with depth while nitrate (NO 3 − ) concentrations increased, indicating nitrification of organic matter within the profiles. The depth of NO 3 − production and transport was seasonally dependent. In winter, NO 3 − and O 2 were exported beyond 100-cm depth; whereas in summer, shallow nitrification zones were underlain by deeper NO 3 − reduction zones, and diel patterns of O 2 and NO 3 − penetration depths were observed. Microbial community compositions and stable isotope profiles (δ 15 N[NO 3 − ], δ 18 O[NO 3 − ], δ 18 O[O 2 ]) were consistent with apparent C–N–O reaction stoichiometries indicating O 2 reduction and nitrification in shallower porewater, followed by varying NO 3 − reduction at depth. Maximum porewater NO 3 − concentrations (∼10–20 μM) were limited by infiltrating O 2 concentrations and C/N ratios of reacting organic matter. Lake-water level variations caused changes in shoreline position and porewater velocities, while variations in lake-water temperature, DOC, and O 2 contributed to changes in reaction rates and depth of O 2 and NO 3 − penetration into the lakebed. The quality of groundwater recharged by lake water reflected temporally and spatially varying physical and biogeochemical processes in the sediment porewater.

Journal of Geophysical Research: Biogeosciences↗
USGS2024

Controls on stable methane isotope signatures in northern peatlands and potential shifts in signatures under permafrost thaw scenarios

Northern peatlands are a globally significant source of methane (CH 4 ), and emissions are projected to increase due to warming and permafrost loss. Understanding the microbial mechanisms behind patterns in CH 4 production in these systems will be key to predicting annual emissions changes, with stable carbon isotopes (δ 13 C-CH 4 ) being a powerful tool for characterizing these drivers. Given that δ 13 C signatures of CH 4 are used in top-down atmospheric inversion models to partition sources, our ability to model CH 4 production pathways and associated δ 13 C-CH 4 signatures in peatland types impacted by a changing climate is critical. We sought to characterize the role of environmental conditions, including both hydrologic and vegetation patterns associated with permafrost thaw, on δ 13 C-CH 4 signatures from a diverse set of high-latitude peatlands. We measured porewater and emitted CH 4 stable isotopes, pH, and vegetation composition from five boreal-Arctic peatlands. Porewater δ 13 C-CH 4 was strongly associated with peatland type, with δ 13 C enriched values obtained from more minerotrophic fens (-61.2 ± 9.1‰) compared to permafrost-free bogs (-74.1 ± 9.4‰) and raised permafrost bogs (-81.6 ± 11.5‰). Variation in porewater δ 13 C-CH 4 was best explained by sedge cover, CH 4 concentration, and the interactive effect of peatland type and pH ( r 2 = 0.50, p < 0.001). Emitted δ 13 C-CH 4 varied greatly but was positively correlated with porewater δ 13 C-CH 4 , suggesting that porewater data can be used to predict changing emissions signatures from these systems. We calculated a weighted mean mixed atmospheric CH 4 signature for northern peatlands of -65.3 ± 7‰ and show that this signature is more sensitive to landscape drying (4 to 10 % depletion in δ 13 C) than wetting (1.5 to 5% enrichment in δ 13 C) under permafrost thaw scenarios. Our results suggest northern peatland δ 13 C-CH 4 signatures are likely to shift in the future which has important implications for source partitioning in atmospheric inversion models.

Journal of Geophysical Research: Biogeosciences↗
USGS2025

Quantifying the effect of petrogenic carbon on SOC turnover for two Rocky Mountain soils: When are petrogenic carbon corrections required?

Petrogenic organic carbon (OC petro ), derived from sedimentary rocks, is an often overlooked and poorly quantified source of soil organic carbon (SOC), which may influence measured or modeled SOC composition, age, and stability. In this study, we exploited differences in thermochemical stability between OC petro and biogenic SOC (OC bio ) using stepped elemental analysis to quantify the fractional contribution of OC petro to total SOC (f petro ), and we conducted a sensitivity analysis to estimate the effects of OC petro on modeled SOC transit times and system ages. Specifically, we compared the effects of accounting for OC petro inputs in SOC turnover modeling (using SoilR) for two montane meadow soils that are underlain by Cretaceous Mancos Shale. At these sites, we estimate that OC petro comprises 7%–9% of the total SOC stock (f petro = 0.07–0.09). However, accounting for OC petro as a mixture of inert and passive C or as completely inert C had negligible effects on SOC transit times and system ages, suggesting that there is a threshold of OC petro content under which there is minimal effect on calculated SOC turnover. Based on our sensitivity analysis, we estimate this threshold to be f petro = 0.125, further supporting that the accurate calculation of OC petro remains an important factor in estimating SOC turnover.

Colorado↗
USGS2024

Responses of marginal and intrinsic water-use efficiency to changing aridity using FLUXNET observations

According to classic stomatal optimization theory, plant stomata are regulated to maximize carbon assimilation for a given water loss. A key component of stomatal optimization models is marginal water-use efficiency (mWUE), the ratio of the change of transpiration to the change in carbon assimilation. Although the mWUE is often assumed to be constant, variability of mWUE under changing hydrologic conditions has been reported. However, there has yet to be a consensus on the patterns of mWUE variabilities and their relations with atmospheric aridity. We investigate the dynamics of mWUE in response to vapor pressure deficit (VPD) and aridity index using carbon and water fluxes from 115 eddy covariance towers available from the global database FLUXNET. We demonstrate a non-linear mWUE-VPD relationship at a sub-daily scale in general; mWUE varies substantially at both low and high VPD levels. However, mWUE remains relatively constant within the mid-range of VPD. Despite the highly non-linear relationship between mWUE and VPD, the relationship can be informed by the strong linear relationship between ecosystem-level inherent water-use efficiency (IWUE) and mWUE using the slope, m *. We further identify site-specific m * and its variability with changing site-level aridity across six vegetation types. We suggest accurately representing the relationship between IWUE and VPD using Michaelis–Menten or quadratic functions to ensure precise estimation of mWUE variability for individual sites.

Journal of Geophysical Research Biogeosciences↗
USGS2024

A semi-mechanistic model for partitioning evapotranspiration reveals transpiration dominates the water flux in drylands

Popular evapotranspiration (ET) partitioning methods make assumptions that might not be well-suited to dryland ecosystems, such as high sensitivity of plant water-use efficiency (WUE) to vapor pressure deficit (VPD). Our objectives were to (a) create an ET partitioning model that can produce fine-scale estimates of transpiration (T) in drylands, and (b) use this approach to evaluate how climate controls T and WUE across ecosystem types and timescales along a dryland aridity gradient. We developed a novel, semi-mechanistic ET partitioning method using a Bayesian approach that constrains abiotic evaporation using process-based models, and loosely constrains time-varying WUE within an autoregressive framework. We used this method to estimate daily T and weekly WUE across seven dryland ecosystem types and found that T dominates ET across the aridity gradient. Then, we applied cross-wavelet coherence analysis to evaluate the temporal coherence between focal response variables (WUE and T/ET) and environmental variables. At yearly scales, we found that WUE at less arid, higher elevation sites was primarily limited by atmospheric moisture demand, and WUE at more arid, lower elevation sites was primarily limited by moisture supply. At sub-yearly timescales, WUE and VPD were sporadically correlated. Hence, ecosystem-scale dryland WUE is not always sensitive to changes in VPD at short timescales, despite this being a common assumption in many ET partitioning models. This new ET partitioning method can be used in dryland ecosystems to better understand how climate influences physically and biologically driven water fluxes.

New Mexico↗
USGS2024

Challenging ring-current models of the Carrington storm

A detailed analysis is made of horizontal-component geomagnetic-disturbance data acquired at the Colaba observatory in India recording the Carrington magnetic storm of September 1859. Prior to attaining its maximum absolute value, disturbance at Colaba increased with an e -folding timescale of 0.46 hr (28 min). Following its maximum, absolute disturbance at Colaba decreased as a trend having an e -folding timescale of 0.31 hr (19 min). Both of these timescales are much shorter than those characterizing the drift period of ring-current ions. Furthermore, over one 28-min interval when absolute disturbance was increasing, the data indicate an absolute rate of change of ≥2,436 nT/hr. If this is representative of disturbance generated by a symmetric magnetospheric ring current, then, assuming a standard and widely used parameterization, an interplanetary electric field of ≥451 mV/m is indicated. An idealized and extreme solar-wind dynamic pressure could, conceivably, reduce this bound on the interplanetary electric field to ≥202 mV/m. If the parameterization for electric-field extrapolation is accurate, but the field strengths obtained are deemed implausible, then it can be concluded that the Colaba disturbance data were significantly affected by partial-ring, field-aligned, or ionospheric currents. The same conclusion is supported by the shortness of the e -folding timescales characterizing the Colaba data. Several prominent studies of the Carrington event need to be reconsidered.

Journal of Geophysical Research Space Physics↗
USGS2024

Limited preservation of strike-slip surface displacement in the geomorphic record

Offset geomorphic markers are commonly used to interpret slip history of strike-slip faults and have played an important role in forming earthquake recurrence models. These data sets are typically analyzed using cumulative probability methods to interpret average amounts of slip in past earthquakes. However, interpretation of the geomorphic record to infer surface slip history is complicated by slip variability, measurement uncertainty, and modification of offset features in the landscape. To investigate how well geomorphic data record surface slip, we use offset measurements from recent strike-slip surface ruptures ( n = 39), faults with geomorphic evidence of multiple strike-slip earthquakes ( n = 29), and synthetic slip distributions with added noise ( n > 10,000) to examine the constraints of the geomorphic record and the underlying assumptions of the cumulative offset probability distribution analysis method. We find that the geomorphic record is unlikely to resolve more than two paleo-slip distributions, except in specific cases with low slip variability, high slip-per-event, and semiarid climate. In cases where site-specific conditions allow for interpretation of more than two earthquakes, lateral extrapolation along a fault is not straightforward because on-fault displacement and distributed deformation may be spatially variable in each earthquake. We also find that average slip in modern earthquakes is adequately recovered by probability methods, but the reported prevalence of strike-slip faults with characteristic slip history is not supported by geomorphic data. We also propose updated methods to interpret slip history and construct uncertainty bounds for paleo-slip distributions.

Journal of Geophysical Research: Solid Earth↗
USGS2024

Controls on the stratigraphic architecture of the US Atlantic margin: Processes forming the accommodation space

Accommodation space governs the spatial and temporal distributions of sediments in continental margins. Mapping the sedimentation patterns, therefore, offers insights into the solid-Earth processes that shape accommodation space. We assembled an unprecedented amount of seismic and borehole data along the Eastern North American Margin and used it to divide the margin's sedimentary package into eight chronostratigraphic intervals, identifying temporal shifts in depocenters under the continental shelf, slope, and rise. The Jurassic depocenters follow the syn-rift structure and its thermal subsidence loci. The Long Island Platform is the only margin segment where the early post-rift sediment thickness matches subsidence predictions from uniform-stretching models, whereas in Georges Bank Basin (GBB) and Baltimore Canyon Trough (BCT), sediment thickness is 1.5–3 times higher than predicted, pointing to other factors at play. A margin-wide Jurassic transient shoulder uplift is inferred from the occurrence of stratigraphic onlaps above thinned crust. Unlike the Jurassic, the Cretaceous and Cenozoic depocenters disregard the inherited subsidence pattern. The accommodation space over the shelf and coastal plain during the Cretaceous was affected by regional isostatic compensation of the sedimentary loads accumulated on the shelf and rise. Accommodation space development in the GBB was interrupted during the Cretaceous after the margin crossed the Great Meteor Hotspot track, resulting in a widespread permanent uplift, erosion, and sediment redistribution. The distribution of anomalous Neogene subsidence in the BCT challenges previous suggestions of mantle dynamic control on the accommodation space and favors flexural downwarping of the shelf by sediment accumulation on the rise.

Journal of Geophysical Research: Solid Earth↗
USGS2024

Earthquake cycle mechanics during caldera collapse: Simulating the 2018 Kīlauea eruption

In multiple observed caldera-forming eruptions, the rock overlying a draining magma reservoir dropped downward along ring faults in sequences of discrete collapse earthquakes. These sequences are analogous to tectonic earthquake cycles and provide opportunities to examine fault mechanics and collapse eruption dynamics over multiple events. Collapse earthquake cycles have been studied with zero-dimensional slider-block models, but these do not account for the complicated interplay between fluid and elastic dynamics or for factors such as the heterogeneous fault properties and non-vertical ring fault geometries often inferred at volcanoes. We present two-dimensional axisymmetric mafic piston-like collapse earthquake cycle models that include rate-and-state friction, fully-dynamic elasticity, and compressible viscous fluid magma flow. We demonstrate that collapse earthquake intervals and magnitudes are highly sensitive to inertial effects, evolving stress fields, fault geometry, and depth-varying fault friction. Given the consistent earthquake cycles observed in most eruptions, this suggests that ring faults can quickly stabilize and often become nearly vertical at depth. We use the well-monitored 2018 collapse sequence at Kı̄lauea as a case study. Our model can produce many features of Kı̄lauea seismic and geodetic observations, except for a significant amount of interseismic slip, which cannot be readily explained with simple rate-and-state friction parameterizations.

Hawaii↗
USGS2024

Afterslip and creep in the rate-dependent framework: Joint inversion of borehole strain and GNSS displacements for the Mw 7.1 Ridgecrest earthquake

The elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub-nanostrain precision and sub-second sampling rate of BSMs bridges a gap between conventional seismologic and geodetic methods, exemplified by atypical postseismic shear strain reversals observed at nearfield (<2 km) station B921 that remain unexplained. We jointly invert GNSS displacements and BSM strains for coseismic and postseismic slip spanning hours to months over 7 independent periods. Cosiesmically, our model resolves the largest slip magnitudes of up to 6.6 m on the mainshock rupture plane, with similar patterns to other inferred slip distributions. The foreshock fault appears to slip coincidently with mainshock, revealing potential asperities activated during the preceding Mw 6.4 event. Postseismically, the best-fitting models adhere to mechanical rate-and-state expectations of logarithmically decaying slip adjacent to the coseismic rupture terminus, and where deep rheologic conditions favor creep. Most spatial variation occurs in the early postseismic timeframe (<1–2 weeks), with evidence for regional rheologic control and static stress dependence. Triggered creep on the neighboring Garlock Fault unexpectedly persists for >178 days—further highlighting the importance of fault networks in postseismic stress redistribution, critical to assessing future hazard.

Journal of Geophysics Research↗
USGS2024

Subduction zone geometry modulates the megathrust earthquake cycle: Magnitude, recurrence, and variability

Megathrust geometric properties exhibit some of the strongest correlations with maximum earthquake magnitude in global surveys of large subduction zone earthquakes, but the mechanisms through which fault geometry influences subduction earthquake cycle dynamics remain unresolved. Here, we develop 39 models of sequences of earthquakes and aseismic slip (SEAS) on variably-dipping planar and variably-curved nonplanar megathrusts using the volumetric, high-order accurate code tandem to account for fault curvature. We vary the dip, downdip curvature and width of the seismogenic zone to examine how slab geometry mechanically influences megathrust seismic cycles, including the size, variability, and interevent timing of earthquakes. Dip and curvature control characteristic slip styles primarily through their influence on seismogenic zone width: wider seismogenic zones allow shallowly-dipping megathrusts to host larger earthquakes than steeply-dipping ones. Under elevated pore pressure and less strongly velocity-weakening friction, all modeled fault geometries host uniform periodic ruptures. In contrast, shallowly-dipping and sharply-curved megathrusts host multi-period supercycles of slow-to-fast, small-to-large slip events under higher effective stresses and more strongly velocity-weakening friction. We discuss how subduction zones' maximum earthquake magnitudes may be primarily controlled by the dip and dimensions of the seismogenic zone, while second-order effects from structurally-derived mechanical heterogeneity modulate the recurrence frequency and timing of these events. Our results suggest that enhanced co- and interseismic strength and stress variability along the megathrust, such as induced near areas of high or heterogeneous fault curvature, limits how frequently large ruptures occur and may explain curved faults' tendency to host more frequent, smaller earthquakes than flat faults.

Journal of Geophysical Research: Solid Earth↗
USGS2024

Brittle regime slip partitioned damage and deformation mechanisms along the eastern Denali fault zone in southwestern, Yukon

Rare bedrock exposures of the eastern Denali fault zone in southwestern Yukon allow for the measurement, sampling, and analyses of brittle regime fault slip data and deformation mechanisms to explore relations to far field, oblique plate motions. Host rock lithologies and associated slip surfaces show episodic damage zone‐related deformation and calcite ± hematite ± chlorite related hydrothermal fluid flow. This regional scale network of asymmetric fault damage is spatially and kinematically linked to a discrete and narrow fault core. Fault network observations, orientations, slip data, and strain inversions document a slip partitioned strike‐slip fault system with locally and mutually overprinting strike‐, oblique‐, and dip‐slip components. Microstructural analyses reveal crystal plastic and co‐seismic brittle deformation mechanisms active in a narrow range of upper crustal temperature, pressure, fluid, and chemical conditions. The net damage related slip is not exclusively formed by a single kinematic system, but rather a fully partitioned, time integrated system likely operative for much of the fault's brittle regime evolution temporally constrained by previously published thermochronometric data. Although the fault slip data was collected from outcrop‐scale exposures at sites tens of kilometers apart, results show remarkable correlation between fault kinematics and plate motions along the ∼580 km long eastern Denali fault segment. End member, subhorizontal, northeast directed reverse and north directed dextral strike slip fault strain axes closely reflect relative plate motion interactions over at least the last 30 m.y. and act as a proxy for far‐field stresses compatible with the kinematics of the damage zone network.

Alaska↗
USGS2025

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington↗
USGS2025

Timescales of surface faulting preservation in low-strain intraplate regions from landscape evolution modeling and the geomorphic and historical record

Large surface-rupturing intraplate earthquakes in stable continental regions (SCRs) are uncommon globally and have recurrence intervals of thousands to hundreds of thousands of years based on the paleoseismic and geomorphic record, challenging accurate active fault identification in these regions. To constrain the timescales of preservation for scarps created by surface ruptures from dip-slip earthquakes, we use a two-dimensional scarp diffusion model for typical intraplate settings and explore which parameters influence fault scarp preservation. These parameters include the coseismic vertical surface offset, the recurrence interval of similar magnitude earthquakes, diffusivity (as a proxy for mean annual precipitation rate), and the erodibility of the surficial material. We constrain parameter ranges from a compilation of historical surface ruptures in intraplate settings in a variety of climates, including the Central and Eastern United States, Australia, Europe, Central Asia (Mongolia, China), India, and West Africa. The timescales of scarp preservation from landscape evolution modeling agree well with observations of scarp preservation in low-strain SCR and intraplate tectonic settings, with some notable exceptions for Australian scarps. We find that the erodibility of the surficial material and earthquake recurrence interval have a stronger effect on the timescales of scarp preservation than diffusivity or coseismic vertical surface offset. Our model results may aid in identifying and characterizing subtle, slow-moving active faults in low-strain SCR and intraplate tectonic settings for different tectonic, geomorphic, and climatic characteristics. Accurate fault locations and characterization from the landscape record has implications for both probabilistic seismic and fault displacement hazard analyses.

Journal of Geophysical Research Solid Earth↗
USGS2025

Crustal structure along and surrounding the Hawaiian Islands: Volcanic island construction across scarred oceanic crust

The Hawaiian Ridge is a classic example of an intraplate volcanic island chain emplaced on oceanic lithosphere. We seek to constrain both the deformation from island loading around the Hawaiian Ridge and the influence of the oceanic lithosphere, including the Molokaʻi fracture zone (MFZ), on off-axis volcanic emplacement. To examine these processes, we conducted a marine geophysical experiment in 2018 that included the acquisition of eight multichannel seismic reflection lines and used a 6,600 in 3 tuned air gun array and an ultra-long hydrophone streamer cable towed behind R/V Marcus G. Langseth across and around the southern Hawaiian Ridge. We image both the top of igneous oceanic crust and the Moho, and we observe significant variations in regional crustal structure and thickness variations, Moho characteristics, the locations of the buried MFZ, and the flexure of the Pacific oceanic lithosphere under the Hawaiian Ridge. We observe up to ∼4.5 km of deflection of the igneous oceanic crust in response to the volcanic load with sediment thickness increasing to ∼3–3.4 km near the ridge but not correlating with the deflection. A systematic difference in igneous oceanic crustal thickness is found north (average ∼5.2 km) and south (average ∼5.9 km), of the MFZ that also correspond to changes in Moho characteristics. The MFZ itself is associated with the largest crustal thickness variations (∼3.7–7.6 km). Magma ascent through these different crustal regions may account for some variations of magmatic flux to the surface along the Hawaiian Ridge.

Hawaii↗