Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “JGR Earth Surface”

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 37 records · Page 2Linked to original sources
USGS2025

Outwash events inhibit vegetation recovery and prolong coastal vulnerability

Overwash, when high ocean water levels and waves flood a coastline, is a common phenomenon that can lead to washover deposits and barrier rollover. Outwash, by contrast, involves seaward flow, often driven by high back-barrier water levels, and can produce washout channels and nearshore deposition. Our observations show that washout channels were quickly (days to weeks) filled and reshaped into beaches, berms, and washover deposits and ponds often formed at the landward ends. However, there was a significant delay in revegetation of former washout areas compared with washover areas. North Core Banks, North Carolina, was affected by repeat hurricanes in different ways: Hurricane Florence (2018) deposited large washover fans 0.5–1 m thick, and Hurricane Dorian (2019) removed 1–4 m of sediment from washout channels. Aerial surveys captured vegetation recolonization on the Florence washover fans within a year but, after Dorian, surveys showed that although the washout channels and ponds quickly filled with marine sand, the channel throats and new washover platforms remained mostly unvegetated for five years. New vegetation growth was associated with the washout ponds and was characteristic of low-elevation hydrophilic environments. We observed comparable outcomes at washout and washover locations on the coasts of Texas and New York and suggest that outwash interrupts the normal cycle of vegetation and dune growth that is key to rebuilding barrier islands after storms. The lack of vegetation in the former washout channels prolongs vulnerability to overwash, further delaying recovery. Our findings have implications for best-management practices and modeling of coastal geomorphic evolution.

JGR Earth Surface↗
USGS2025

Dynamic feedbacks between river meandering and landsliding in northwestern Washington glacial terraces

Landsliding in river valleys poses unique risks for cascading hazards and can damage infrastructure and cause fatalities. In postglacial valleys, many landslides are posited to occur in relation to lateral river erosion, but the dynamics of fluvial-hillslope interactions are not well understood. Here, we investigate a section of the Nooksack River in western Washington State where the channel is flanked by landslide-prone glacial terraces similar to those that failed in the 2014 State Route 530 “Oso” landslide. We map 216 landslides through time across 17 aerial imagery data sets (1933–2022) and analyze them in relation to river meandering and curvature. We observe dynamic feedbacks between lateral river meandering and valley-adjacent landsliding. Terrace lateral retreat rates of up to 25 m/year owing to combined fluvial erosion and slope failure occur on pinned, outer meander bends immediately downstream from peaks in river curvature (>0.0075 1/m); these locations are predisposed to both shallow and deep-seated landslides. Deep-seated landslides extending 17%–32% of the active valley width into the floodplain can displace the river away from the floodplain margin and change the channel planform. River-displacing landslides relocate meanders up- or downstream, thereby conditioning the location of subsequent landslides. This conceptual model of coupled landslide-driven meander displacement and valley-adjacent landsliding is exemplified across western Washington river systems. The distance between up- and downstream valley-adjacent landsliding scales with valley width, meander wavelength, and terrace height. Our results can advance our understanding of the river-hillslope interface in landscape evolution and can be used to inform hazard management in river corridors.

Washington↗
USGS2025

Numerical simulation of sound-side barrier-island inundation and breaching during Hurricane Dorian (2019)

Hurricane-induced morphological changes and associated community hazards along sandy, barrier-island coastlines have been studied primarily from the perspective of ocean-side attack by storm-driven ocean surge and large waves. Thus, our understanding of long-term barrier island morphological change focuses on beach erosion, overwash, and inlet formation. In contrast, outwash events with inundation from the sound side, such as one that occurred in Cape Lookout National Seashore, North Carolina, USA during Hurricane Dorian (September 2019), are understudied. Studying such events can improve understanding of barrier island response and stability for a broader range of conditions. Here, we model the hydrodynamics and morphological evolution of a barrier island using a coupled wave-current-sediment transport modeling system. Wind-driven surge in Pamlico Sound led to overtopping from the sound side, which eroded outwash channels and transported sediment seaward into the nearshore. Simulations reproduce the channel features observed with aerial imagery and provide information not available from the remote-sensing observations, including channel depths (>2 m) and the fate of the eroded sand. We found that >99% of the eroded sand was deposited in the nearshore, within 1,000 m of the shoreline in depths <10 m, suggesting that the deposited sediment remains available for littoral transport and beach recovery. Simulations with combinations of coarse or fine sediment and vegetated or unvegetated landcover indicate that channel position did not vary with grain size or vegetation, while volume of erosion and channel morphology were more responsive to variations in grain size and less responsive to presence of vegetation.

North Carolina↗
USGS2025

Modeling the influence of upper and lower shoreface dynamics on barrier island evolution

Barrier island resilience to climate impacts depends on sediment redistribution between the subaqueous shoreface and subaerial barrier during sea-level rise and storms. However, autogenic interactions between the upper and lower shoreface and their influence on the subaerial barrier are poorly characterized. Here, we explore the influences of various shoreface components on barrier morphology using a model of barrier and shoreface evolution under sea-level rise, the Articulated Barrier Shoreface (ABSF) Model. This reduced-complexity model divides the shoreface into upper and lower shoreface panels that respond independently to sea-level rise and deviations from the equilibrium slope. We couple the ABSF with the Lorenzo-Trueba & Ashton, 2014, https://doi.org/10.1002/2013jf002941 model (LTA), a barrier island evolution model driven by overwash and sea-level rise. Through this coupled framework, we examine the influences of upper and lower shoreface slopes, their respective depths, and sensitivity to wave climate on long-term barrier evolution. Results show that the relative depths of the upper and lower shoreface toes influence barrier response to rising seas, alongside overwash flux and closure depth. Notably, the lower shoreface response to sea-level change lags that of the upper shoreface over decades, diminishing the resilience of the barrier over centennial timescales by slowing the overall barrier response. In fact, the ABSF model predicts barriers will drown faster and more than predicted with a linear shoreface. Results highlight the shoreface as an important sediment reservoir for barrier islands and that differences in upper and lower shoreface responses can reduce barrier resilience to sea-level rise due to limited lower shoreface sediment accessibility.

JGR Earth Surface↗
USGS2025

Are equilibrium shoreline models just convolutions?

Yes. Equilibrium shoreline models, which simulate wave-driven cross-shore erosion and accretion, are mathematically equivalent to a discrete convolution (i.e., a weighted, moving average) of a time series of wave-forcing conditions with a parameterized memory-decay kernel function. The direct equivalence between equilibrium shoreline models and convolutions reveals key theoretical aspects of equilibrium behavior. Convolutions (representing quasi-low-pass filter operations) provide an intuitive theoretical description of shoreline erosion and accretion behavior in response to waves: that is, shoreline position often mirrors the weighted moving average of wave time series. Model-convolution equivalence also provides a conceptual basis to interpret, evaluate, and construct data-driven Machine-Learning/Deep-Learning (ML/DL) models that use convolutions to extract features from data and then apply them for prediction (e.g., Convolutional Neural Networks (CNNs)). Finally, our findings provide a methodological pathway (based on Fourier transforms) for future understanding of wave-driven shoreline change, which can be used to interpret the coherence between the frequency spectrum of the processes of waves and shoreline change and construct more computationally efficient and effective shoreline-modeling approaches.

JGR Earth Surface↗
USGS2026

Impacts of transient critical zone evolution on landslide susceptibility and hazard Across the Utuado batholith, Puerto Rico

Landscapes in humid climates often experience rapid chemical and mechanical weathering that can change groundwater fluxes, drainage pathways, and frictional strength via alterations in subsurface mineral composition and matrix porosity. These changes to hillslope critical zone structure and strength should control landslide susceptibility, size, and failure style. Here we test the impact of differential weathering on landslide failure processes by exploiting a natural experiment along the dioritic Utuado batholith in Puerto Rico's Cordillera Central, where a sharp weathering gradient follows an upstream-propagating knickzone into an uplifted shoreline platform that was affected by Hurricane Maria in 2017. Through a combination of landslide mapping, field geophysical and geotechnical surveys, topographic analysis, calibrated stream profile inversion, 3D back-analysis of hillslope strength, and debris-flow inundation assessments, we show that differential weathering induces changes to critical zone structure and consequently landslide susceptibility and hazard. Our mapping reveals that the hurricane-induced landslides preferentially cluster along less weathered hillslopes despite a marked increase in frictional strength. Jointly inverted seismic refraction and electrical resistivity profiles show that the frictionally weak, heavily weathered hillslopes are highly porous and dry, which we posit acts to stabilize slopes through limiting deeper hydrologic connectivity that can produce shallow water tables with destabilizing lateral seepage forces during storms. Debris-flow inundation mapping and analysis reveals that major contributors to debris-flow growth, channel steepness and landslide coalescence, scale linearly with hillslope friction angle. This leads to higher-mobility debris flows along knickzone reaches, which disproportionately affects communities at the base of the plateau.

Puerto Rico↗
USGS2026

Exploring future ecogeomorphic state transitions of an undeveloped barrier island: Multidecadal model projections for North Core Banks, NC, USA

Ecogeomorphic state transitions in coastal barrier systems will be impacted and likely amplified by projected accelerated sea-level rise and enhanced storm intensity over the next several decades. When, where, and how these transitions will take place remain uncertain and complex, yet are of primary interest to coastal managers looking to prepare for future changes. Here, we apply the new Mesoscale Explicit Ecogeomorphic Barrier model to generate probabilistic projections of future ecogeomorphic state transitions up to year 2050 for the undeveloped North Core Banks (NCB), North Carolina, United States. In our full probabilistic projection, which accounts for stochasticity from uncertainty in sea-level rise and high-water event (HWE) projections, loss of barrier interior to passive sea-level rise submergence dominates ecogeomorphic state transitions, suggesting that the biggest threat to many undeveloped barriers in the next several decades will come from the back-barrier rather than ocean-side. Additionally, this projection suggests that dune and unvegetated interior states are most volatile, with the intensity of state transitions varying significantly alongshore. In projections simulating a range of deterministic sea-level rise rates, higher rates enhance: overwash frequency; loss of beach and dune states; and interior submergence. We also find that state transitions are less sensitive to HWE intensity than sea-level rise, and that more frequent, moderate HWEs—rather than the largest events—may be the principal drivers of mesoscale ecogeomorphic state change. Our projections highlight important future mesoscale trends of ecogeomorphic state transitions for NCB that are likely relevant to other stable, sediment-rich barrier systems of U.S. coastlines.

North Carolina↗
USGS2023

On the scale-dependence of fault surface roughness

Defining roughness as the ratio of height to length, the standard approach to characterize amplitudes of single fault, joint and fracture surfaces is to measure average height as a function of profile length. Empirically, this roughness depends strongly on scale. The ratio is approximately 0.01 at a few mm but 10× smaller at a few tens of meters. Surfaces are rougher at small scales. However, these conclusions are metric-dependent. If instead height is averaged over wavelength, roughness is nearly Brown spatial noise, having almost scale-independent apparent surface height to wavelength ratio. The small deviation from scale-independence is of the opposite sense than found using the standard metric; surfaces are slightly rougher at long wavelengths. Some natural surfaces may be Brownian within the measurement uncertainties. These contradictions are curiosities of surfaces that have Hurst exponents between 0.5 and 1, as natural fault surfaces do. The wavelength-based analysis of roughness and how it changes with scale are straight-forward; a normalized Fourier transform approximately preserves amplitude and its scale dependence in the wavelength domain. Among the conclusions from reconsideration of scale dependence are that the scale dependence is weak and much smaller than that of other fault and shear zone properties. Background and aftershock seismicity, jogs and step-overs indicate strong localization (smoothing) with slip and scale. The lack of strong scale dependence to surface roughness suggests it is not the dominant control on brittle shear zone evolution.

JGR Solid Earth↗
USGS2023

Microscale spatial variations in coseismic temperature rise on hematite fault mirrors in the Wasatch fault damage zone

Coseismic temperature rise activates fault dynamic weakening that promotes earthquake rupture propagation. The spatial scales over which peak temperatures vary on slip surfaces are challenging to identify in the rock record. We present microstructural observations and electron backscatter diffraction data from three small-displacement hematite-coated fault mirrors (FMs) in the Wasatch fault damage zone, Utah, to evaluate relations between fault properties, strain localization, temperature rise, and weakening mechanisms during FM development. Millimeter- to cm-thick, matrix-supported, hematite-cemented breccia is cut by ∼25–200 μm-thick, texturally heterogeneous veins that form the hematite FM volume (FMV). Grain morphologies and textures vary with FMV thickness over μm to mm lengthscales. Cataclasite grades to ultracataclasite where FMV thickness is greatest. Thinner FMVs and geometric asperities are characterized by particles with subgrains, serrated grain boundaries, and(or) low-strain polygonal grains that increase in size with proximity to the FM surface. Comparison to prior hematite deformation experiments suggests FM temperatures broadly range from ≥400°C to ≥800–1100°C, compatible with observed coeval brittle and plastic deformation mechanisms, over sub-mm scales on individual slip surfaces during seismic slip. We present a model of FM development by episodic hematite precipitation, fault reactivation, and strain localization, where the thickness of hematite veins controls the width of the deforming zones during subsequent fault slip, facilitating temperature rise and thermally activated weakening. Our data document intrasample coseismic temperatures, resultant deformation and dynamic weakening mechanisms, and the length scales over which these vary on slip surfaces.

Utah↗
USGS2025

Multi-scale geophysical imaging of a hydrothermal system in Yellowstone National Park, USA

Little is known about the local plumbing systems that fuel Yellowstone’s famous hot springs, geysers and mud pots. A multi-method, multi-scale geophysical investigation was carried out in the Obsidian Pool Thermal Area (OPTA) to: (i) delineate the lateral extent of the hydrothermal area and associated surface features; (ii) estimate the dimensions of the upflow zone and identify its main controlling structures; (iii) assess fluids circulation pathways from depth to surface. Ground and airborne geophysical data were acquired to connect local and regional scales, from shallow to large depths. Maps of surface electrical resistivity show a strong correlation with hydrothermal features. At in-termediate depths, electrical resistivity permits delineating the upper limit of the upflow zone, while Poisson’s ratio highlights differences in subsurface fluid content. Combining these results with surface observations and topographic information, we speculate that differential mixing of hydrothermal and fresh water could explain the wide diversity of features observed at OPTA. Low electrical resistivity observed at large depths also suggest that a vast upflow zone, controlled by rhyolite flows and conjugate faults, underlies the OPTA. We speculate that hydrothermal fluids rise along fractures and reach the surface in topographic lows to form hydrothermal features. Our results show that synoptic, multi-scale geophysical measurements provide a roadmap for understanding where and how geologic heterogeneity, topography, fluid-gas separation, and the mixing of thermal and meteoric waters conspire to produce the wide variety of Yellowstone’s renowned hydrothermal features.

Wyoming↗
USGS2022

Power-law viscoelastic flow of the lower accretionary prism in the Makran subduction zone following the 2013 Baluchistan Earthquake

Subduction zone accretionary prisms are commonly modeled as elastic structures where permanent deformation is accommodated by faulting and folding of otherwise elastic materials, yet accretionary prisms may exhibit other deformation styles over relatively short time scales. In this study, we use 6.5-year (2014–2021) Sentinel-1 interferometric synthetic aperture radar (InSAR) time-series of post-seismic deformation in the Makran accretionary prism of southeast Pakistan to characterize non-linear viscoelastic deformation within an active accretionary prism on short timescales (months to years). We constructed a series of 3-D finite-element models of the Makran subduction zone, including an accretionary prism, and constrained the elastic thickness of the upper wedge and the flow-law parameters (power-law exponent, activation enthalpy, and pre-exponential constant) of the lower wedge through forward model fits to the InSAR time-series. Our results show that the prism is elastically thin (8–12 km) and the non-linear viscoelastic relaxation of the deep portions of the prism alone can sufficiently explain the post-seismic surface deformation. Our best fitting flow-law parameters ( n = 3.76 ± 0.39, Q = 82.2 ± 37.73 kJ mol −1 , and A = 10 −3.36±4.69 ) are consistent with triggering of low temperature dislocation creep within fluid-saturated siliciclastic rocks. We believe that the fluids necessary for this weakening originate from sedimentary underplating and/or the presence the hydrocarbons. The presence of power-law rheology within the lower wedge impacts the estimated plate coupling and the stress state in the subduction system, with respect to the conventional elastic wedge model, and hence should to be considered in future earthquake cycle models.

Makran accretionary prism↗
USGS2020

Seismic reflection imaging of the low-angle Panamint normal fault system, eastern California

Shallowly dipping (<30°) low‐angle normal faults (LANFs) have been documented globally; however, examples of active LANFs in continental settings are limited. The western margin of the Panamint Range in eastern California is defined by a LANF that dips west beneath Panamint Valley and has evidence of Quaternary motion. In addition, high‐angle dextral‐oblique normal faults displace middle to late Quaternary alluvial fans near the range front. To image shallow (<1 km depth), crosscutting relationships between the low‐ and high‐angle faults along the range front, we acquired two high‐resolution P wave seismic reflection profiles. The northern, 4.6‐km‐long profile crosses the 2‐km‐wide Wildrose graben and the southern, 0.8‐km‐long profile extends onto the Panamint Valley playa, ~7.5 km S of Ballarat, CA. The profile across the Wildrose graben reveals a robust, low‐angle reflector interpreted to represent the LANF separating Plio‐Pleistocene alluvial fanglomerate and Proterozoic metasedimentary deposits. High‐angle faults interpreted in the seismic profile correspond to fault scarps on Quaternary alluvial fan surfaces. Interpretation of the reflection data suggests that the high‐angle faults vertically displace the LANF up to 80 m within the Wildrose graben. Similarly, the profile south of Ballarat reveals a low‐angle reflector, which appears both rotated and displaced up to 260 m by high‐angle faults. These results suggest that near the Panamint range front, the high‐angle faults are the dominant active structures. We conclude that at least at shallow (<1 km) depths, the LANF we imaged is not active today.

California↗
USGS2026

Impacts of cross-shore winds, waves, and sand ridges on inner shelf circulation

Understanding of circulation on inner continental shelves needs to be advanced to accurately predict the transport and connection between the shelves and nearshore. The circulation over the storm-dominated Fire Island Shelf, New York and a series of in situ shoreface connected sand ridges was investigated through high-resolution field observations, numerical simulations of hydrodynamics, and principal component analyses. The along-shelf subtidal flows were more energetic than the cross-shelf subtidal flows. Alongshore winds were driving the along-shelf subtidal flow and the variability of water levels. The cross-shelf subtidal flows were dominantly onshore near the seabed and offshore near the surface: the processes significant in this two-layer exchange were cross-shore winds which were commonly overlooked in literature, alongshore winds through Coriolis, and bottom friction. The ridges were significant in changing the magnitude and even the direction of the pressure gradients and associated flows. Wave impacts including water level changes, Stokes drift, vortex force, and breaking became significant in the inner shelf and the nearshore during wave-energetic events. Depending on the directions of waves and winds, waves caused to accelerate or decelerate and even reverse the wind-induced flow. In case of more shore-normal waves, the cross-shore Stokes drift became more significant; its interaction with Earth's rotation and flow vorticity generated significant Stokes-Coriolis impacts and vortex force, respectively, in the alongshore. With results applicable to other shelf settings, this study advanced the understanding of spatial and temporal variations of flows along and across inner shelves under varying forcing conditions.

New York↗
USGS2025

Magmatic volatile flux drives non-eruptive volcano-tectonic seismicity at Mount St. Helens, USA from 2008-2024

Seismicity during non-eruptive periods is useful for observing stress changes related to magmatic transport and volatile exsolution within active volcanoes. Mount St. Helens in Washington, USA, is the most active volcano in the continental United States and has been in quiescence since 2008. To explore the processes driving seismicity at Mount St. Helens, we create a high-resolution seismicity catalog of non-eruptive seismicity from 2008 through 2023, consisting of 31,133 events. We find persistent shallow seismicity (−2.2 to 2 km below sea level (BSL), 0–4.2 km below the surface) throughout the entire study period that concentrates beneath the dacite dome complex from the 1980–1986 and 2004–2008 eruptions. Additionally, there is frequent deeper seismicity (2–8 km BSL, 4.2–12.2 km below the surface) beginning in 2016. We examine a selection of deep earthquake swarms and find complex along-depth seismicity patterns. Within swarms, increases in shallow seismicity rates can precede or are concurrent with increases of deep seismicity rates. Lastly, we discover a series of semi-periodic, shallow, burst-like swarms, consisting of low-amplitude, repetitive similar earthquakes, indicating periodic valve-like release of fluid pressure from the conduit. Increased seismic activity beginning in 2016 indicates ongoing repressurization within the magmatic system driven by recharge or crystallization-induced second boiling within the upper-crustal reservoir after 2008. The data indicate that non-eruptive seismicity at Mount St. Helens is controlled by fluid pressure changes from gas flux sourced from the magma reservoir that migrates through crack networks.

Washington↗
USGS2025

Three-dimensional seismic velocity models for the San Francisco Bay region, California from joint body-wave and surface-wave tomography validated by waveform simulation

A high-resolution seismic velocity model for both P and S waves is essential for the San Francisco Bay (SFB) region to accurately simulate earthquake ground motions and assess seismic hazards in this high-risk region. We present new three-dimensional P-wave and S-wave velocity (Vp and Vs) models of the SFB region developed using joint tomographic inversion with a rich set of earthquake and active-source body-wave data and ambient-noise surface-wave data. We designed a dynamic balance strategy to effectively balance the contributions of body-wave and surface-wave data during each iteration of the inversion. Owing to the complementary resolution of body-wave and surface-wave data, our new models are significantly improved compared to previous tomographic models. Major geologic features in the SFB region, including low-velocity basins, high-velocity Franciscan Complex, granite, and gabbro rocks, and across-fault velocity contrasts, are imaged in our Vp and Vs models. Compared to the geology-based U.S. Geological Survey velocity model, our new model has better defined structures for some important areas that are of great concern for seismic hazard assessment, such as the Sacramento-San Joaquin Delta, and the 2014 magnitude 6.0 South Napa earthquake area. We validated our model by simulating waveforms for a set of moderate-magnitude earthquakes. The good performance of our model in waveform simulation indicates its potential for improving ground motion simulation and seismic hazard assessment in the SFB region.

California↗
USGS2018

Euler-vector clustering of GPS velocities defines microplate geometry in southwest Japan

I have used Euler-vector clustering to assign 469 GEONET stations in southwest Japan to k clusters ( k = 2, 3,..., 9) so that, for any k , the velocities of stations within each cluster are most consistent with rigid-block motion on a sphere. That is, I attempt to explain the raw (i.e., uncorrected for strain accumulation), 1996–2006 velocities of those 469 Global Positioning System stations by rigid motion of k clusters on the surface of a spherical Earth. Because block geometry is maintained as strain accumulates, Euler-vector clustering may better approximate the block geometry than the values of the associated Euler vectors. The microplate solution for each k is constructed by merging contiguous clusters that have closely similar Euler vectors. The best solution consists of three microplates arranged along the Nankaido Trough-Ryukyu Trench between the Amurian and Philippine Sea Plates. One of these microplates, the South Kyushu Microplate (an extension of the Ryukyu forearc into the southeast corner of Kyushu), had previously been identified from paleomagnetic rotations. Relative to ITRF2000 the three microplates rotate at different rates about neighboring poles located close to the northwest corner of Shikoku. The microplate model is identical to that proposed in the block model of Wallace et al. (2009, https://doi.org/10.1130/G2522A.1) except in southernmost Kyushu. On Shikoku and Honshu, but not Kyushu, the microplate model is consistent with that proposed in the block models of Nishimura and Hashimoto (2006, https://doi.org/10.1016/j.tecto.2006.04.017) and Loveless and Meade (2010, https://doi.org/10.1029/2008JB006248) without the low-slip-rate boundaries proposed in the latter.

southwest Japan↗
USGS2026

Fully contained laboratory earthquakes: The effect of asperity aspect ratio and free surfaces

Corner frequency ( f c ) and seismic moment ( M 0 ) are key parameters derived from seismic signals that are used to characterize earthquake stress drop, rupture area, and slip. These parameters are also affected by fault geometry and boundary conditions. However, the systematic study of these effects in laboratory settings has been challenging. This study presents laboratory earthquake experiments that examine how rupture dynamics are influenced by (a) the aspect ratio of rectangular PMMA velocity-weakening (VW) asperities surrounded by the Teflon velocity-strengthening (VS) patches, and (b) whether the sides of a VW asperity are confined with VS patches or are free surfaces. We found that increasing confinement by reducing free surfaces or increasing the VW asperity aspect ratio stabilizes fault slip, so that higher normal stress is required to transition from aseismic to seismic slip. Increased confinement and high aspect ratios also reduced M 0 and increased f c , both of which were determined from the radiated seismic waves. M 0 and f c were primarily controlled by the shorter dimension of the VW asperity. Analysis of high-frequency acoustic emission signals revealed that ruptures on high-aspect-ratio VW asperities propagated more unidirectionally, whereas ruptures on square VW asperities were more complex. Further, the high-aspect-ratio asperities were more likely to be eroded by surrounding VS regions while low-aspect-ratio asperities were more likely to rupture into the VS surroundings. These results demonstrate that both the confinement from surrounding stable areas and the geometry of the seismogenic patch can affect rupture nucleation, propagation, and seismic source characteristics.

JGR Solid Earth↗
USGS2021

Geophysical constraints on the crustal architecture of the transtensional Warm Springs Valley fault zone, northern Walker Lane, western Nevada, USA

The Walker Lane is a zone of distributed transtension where normal faults are overprinted by strike-slip motion. We use two newly-acquired high-resolution seismic reflection profiles and a reprocessed Consortium for Continental Reflection Profiling (COCORP) deep crustal reflection profile to assess the subsurface geometry of the Holocene-active, transtensional Warm Springs Valley fault zone (WSVFZ) near Reno, Nevada, USA. Our multi-scale observations extend to 12 km depth and suggest that the WSVFZ is more complex in the subsurface than implied by late Pleistocene surface fault traces. Two ~4-km-long high-resolution profiles image to a depth of ~2 km and reveal moderately dipping reflections and truncations, some of which project to mapped scarps formed in late Pleistocene surfaces. The shallow lines are co-located with COCORP profile NV 08 along ~40° N latitude. Re-analysis of the COCORP data reveals previously unidentified coherent reflections to a depth of ~12 km and a previously mapped ~30 west-dipping fault at 8-12 km. From these seismic profiles, the WSVFZ is not a simple, sub-vertical fault zone extending through the entire seismogenic crust. Instead, the reflections are consistent with a zone of steep- and moderately-dipping faults that simplify and steepen with depth before intersecting a mid-crustal, low angle (~25-30°) fault. The complex fault geometry of the WSVFZ implies that crustal shear is accommodated by a mix of dipping and subvertical faults in the transtensional northern Walker Lane. If so, transtensional fault zones may present challenges to paleoseismic and geodetic studies and require careful treatment when included in seismic hazard analyses.

Nevada↗