Geology topics

William J. Stephenson

Publications and source records attributed to William J. Stephenson.

At least 19 recordsLinked to original sources

USGS2026

The Cascadia Region Earthquake Science Center (CRESCENT): Advancing understanding of Cascadia's earthquake hazards

The Cascadia Subduction Zone (CSZ) is capable of producing large megathrust earthquakes and tsunamis, posing significant risks to the Pacific Northwest. This paper examines the seismic characteristics of the CSZ, including its offshore coupling, sparse seismicity, and a complex paleoseismic record. While recent advances in geophysical methods and data collection have improved our understanding of subduction zone processes, important uncertainties remain regarding rupture dynamics, ground shaking, and tsunami impacts. The Cascadia Region Earthquake Science Center (CRESCENT) has been established to address these challenges through coordinated research, education, and community engagement. By integrating scientific studies, fostering collaboration, and leveraging new technologies, CRESCENT seeks to improve hazard assessment and resilience in the region. This paper outlines the primary challenges and opportunities in advancing subduction zone science and enhancing societal preparedness for future CSZ events.

California, Oregon, Washington↗

USGS2026

Cascadia Subduction Zone science: Call for the next generation community seismic velocity model

The Cascadia subduction zone (CSZ) hosts major seismic and tsunami hazards, yet key questions persist about the relationship between margin structure, fluid distribution, episodic tremor and slip, shallow megathrust behavior, shaking and tsunamigenesis, and the resulting hazard estimates. Addressing these problems requires an empirically grounded, three‐dimensional seismic velocity model to illuminate subsurface structure and properties and to provide a basis for geophysical studies such as earthquake simulations and ground‐motion estimation. In May 2024, the National Science Foundation‐funded Cascadia Region Earthquake Science Center (CRESCENT) community velocity model (CVM) working group, with U.S. Geological Survey and regional partners, convened a workshop to identify priorities for such a model. Participants emphasized the features necessary for addressing key science questions, including implementing findability, accessibility, interoperability, and reusability (FAIR) access, capturing along‐strike and along‐dip structural heterogeneity, resolving shallow offshore–onshore structure, constraining elastic properties and quantifying their uncertainties for numerical wave propagation simulations, their validation benchmarks, and supporting associated accurate earthquake ground‐motion simulations and hazard assessments. This article describes the priorities defined in the workshop, and a description of how, guided by these needs, CRESCENT plans to develop multiple generations of a CVM to advance CSZ science and improve seismic and tsunami hazard modeling across the Pacific Northwest. The CVM will span the CSZ from the surface to ∼100 km depth, offshore and east of the Cascades into Idaho (∼132°–110° W) and the southern and northern tectonic regime transitions (∼36°–52° N) to capture the entire tectonic system as well as its surroundings.

Cascadia Subduction Zone↗

USGS2026

Ground-motion characterization for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands

We develop the ground-motion characterization (GMC) for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands (NSHM-PRVI) for earthquakes in active crustal, subduction interface, and subduction intraslab regimes. Using ground-motion models (GMMs) from the Next-Generation Attenuation (NGA)-West2 and NGA-Subduction projects, the GMC is parameterized by scaled-backbone models for median ground motions and by independent logic trees of aleatory variability. We introduce several novel GMC features into the U.S. NSHM: (1) use of regional ground-motion data for modeling median ground motions; (2) development of scaled-backbone models for median ground motions; and (3) development of independent logic trees of aleatory variability from variance components of GMMs and computed from multiple ground-motion datasets, and incorporating regional ground-motion data effects on variability. We compute probabilistic seismic hazard curves and maps to evaluate the contributions from the GMC components and for comparison with the 2003 GMC. Contributions of the GMC to epistemic uncertainty in seismic hazard are evaluated through spatial variations in epistemic uncertainty in hazard maps, comparisons between mean hazard curves and fractiles, as well as investigations of the contributions of logic-tree branches to hazard maps and curves. Comparisons between seismic hazard from the 2025 and 2003 GMCs allow examination of the changes in hazard curves and mapped values with 2% and 10% probabilities of exceedance in 50 years. The 2025 GMC exhibits modest changes in median predictions, relative to the 2003 GMC; overall, values of aleatory variability are higher, except at long periods ( ), resulting in changes to probabilistic ground motions at low probabilities of exceedance (<10% probability of exceedance in 50 years). Changes in hazard at 2% and 10% probabilities of exceedance in 50 years are also relatively modest (within 20%) at most sites, with the impacts from the 2025 GMC exhibiting minor reductions and increases at 0.2 and 1.0 s periods, respectively, relative to the 2003 GMC.

Puerto Rico, U.S. Virgin Islands↗

USGS2026

Quantifying seismic properties of a river channel at Mount Rainier for use in debris flow monitoring and analysis

Theoretical models that relate debris flow properties to their seismic signature suggest seismic methods may be used to remotely characterize properties of these events. However, the complexity of debris flow sources and poorly constrained material properties in near‐surface environments limit our ability to determine important attributes of debris flows, such as volume and particle size distribution, from seismic records alone. In this study, we explore the sensitivity of debris flow seismic signals to subsurface seismic characteristics using an established debris flow seismicity model and subsurface properties derived from active‐source measurements in the Tahoma Creek stream channel at Mount Rainier, United States. Using refraction and multichannel analysis of surface waves, we estimate 1D primary and secondary wave velocity profiles to a depth of 11 m and calculate the frequency‐varying phase velocity ( ⁠⁠ v c ) and Rayleigh‐wave quality factor (Q R ⁠⁠ ) for frequencies between 9 and 50 Hz. We find that the Tahoma Creek stream channel has low v c ⁠ , varying with frequency between 226 and 434 m/s, and is highly attenuating, with Q R estimated below 13.3 at all analyzed frequencies. We model the seismic signal of a hypothetical debris flow in Tahoma Creek using our measured values and compare the results against the same model but varying v c and Q R over a range of frequency‐independent, single values commonly used in the literature when measurements are not available ( ⁠⁠ v c = 250-750 m/s, Q R = 3-33 ⁠ ). We find that the power spectral densities of the modeled debris flows vary by orders of magnitude within the subset of our test values, highlighting the benefits of measuring material properties when using modeled debris flow seismic signals for quantitative monitoring.

Washington↗

USGS2026

Subsurface structure across the Tacoma Basin, Washington State, using trans-dimensional Bayesian inversion of fundamental mode spatial autocorrelation data

Spatial autocorrelation (SPAC), the azimuthal average of the normalized cross-correlation between equidistant station pairs deployed in a 2-D array, is widely used to image the subsurface structure. However, the rigorous estimate of subsurface structure and its uncertainties as a function of depth using SPAC data is challenging due to the nonlinear relation between the SPAC data and Earth structure as well as the trade-off between depth and velocity. Additionally, data noise is strongly correlated due to data processing (e.g. filtering, stacking from multiple time segments and azimuthal averaging). Most studies do not account for the correlated noise and fix the ratio of compressional-wave velocity ( V P ) to shear-wave velocity ( V s ) (i.e. V P / V s ratio) and the number of layers, both of which are typically unknown. To address these challenges, we develop a hierarchical trans-dimensional Bayesian inversion of fundamental mode of SPAC data that properly accounts for the correlated data noise, samples the V P / V s ratio and relaxes the number of layers (i.e. model parametrization) to be unknown in the inversion. We further examine the limitation of using only fundamental modes in the inversion. Our synthetic experiments show that the inversion recovers an incorrect model unless we sample the correlated noise and V P / V s ratio in the inversion. The inversion is then applied to SPAC data acquired at 19 sites across the Tacoma basin in Washington State to characterize the V s and the time-averaged V s over 30-m depth ( V s 30 ). Our results show that the V s 30 varies from ∼200 to 800 m s −1 . The V s 30 within the basin is higher in the middle and lower on the east and west sides. We find that these V s 30 values vary with geologic unit. The uncertainties for V s 30 are within 20 m s −1 in average except for the most eastern site TB28. Additionally, the uncertainties are greater for deeper depths beneath most of the sites as the sensitivity decreases as a function of depth. The Vs structure as a function of depth is also complex beneath some sites, possibly because the SPAC curves are affected by higher order Rayleigh modes that are not considered in the inversion. To better constrain the deeper V s structure, V s 30 and/or other average measures of V s over depth, additional constraints from complementary data, such as ellipticity or geologic data are needed. Moreover, our synthetic experiments show that higher order modes can have significant effect in the inversion results, particularly when there is a low-velocity layer.

Washington↗

USGS2025

Neotectonic origins for the Meadow Bank scarp, Wabash Valley seismic zone USA

The Meadow Bank scarp (MBS) in southeastern Illinois is a linear geomorphic expression, ∼10 km long and ∼8 m high above a relatively flat landscape. It parallels an underlying northeast‐oriented Late‐Precambrian–Early‐Cambrian structural fabric, called the Wabash Valley fault zone, and is within an area of modern, historic, and paleo seismicity, called the Wabash Valley seismic zone. In addition, the MBS acts as a boundary of the Wabash River floodplain, as well as Pleistocene glacial outwash channels, which show evidence of frequent outburst flood events. To better understand the MBS’s equivocal origin in this complex geologic environment, we acquired a 917‐m‐long seismic‐reflection survey across its axis to assess the subsurface geologic configuration. The resultant image indicates a complex set of faults that antiformally fold and displace the top of Paleozoic bedrock by ∼12 m across the survey. Moreover, fault and/or fold deformation extends into the shallowest imaged Quaternary strata at ∼6 m below the ground surface. This suggests the MBS origin is related to underlying Quaternary reactivated, positively inverted faults rather than exclusively to glacial outburst flood erosion. These results provide rare paleoearthquake spatial constraints for central U.S. regional seismic hazard consideration.

Illinois↗

USGS2025

Three-dimensional seismic velocity model for the Cascadia Subduction Zone with shallow soils and topography, version 1.7

The U.S. Geological Survey’s seismic velocity model for the Cascadia Subduction Zone provides P- and S-wave velocity ( V P and V S , respectively) information from 40.2° to 50.0° N. latitude and −129.0° to −121.0° W. longitude, and is used to support a variety of research topics, including three-dimensional (3D) earthquake simulations and seismic hazard assessment in the Pacific Northwest. This report describes an update to the previous version (v) 1.6 of the 3D seismic velocity model for the Cascadia Subduction Zone. This new model (herein referred to as v1.7) contains more detailed near-surface structure for improved earthquake ground motion modeling. Updated features include the addition of a new shallow soil velocity model in the top few hundred meters and the option of adding user-specified topography. Although v1.6 of the Cascadia seismic velocity model has a minimum V S of 600 meters per second (m/s), the new model (v1.7) has a minimum V S of approximately 40 m/s. Overall, this update will allow for more accurate ground motion estimates from 3D simulations of scenario earthquakes in the Cascadia Subduction Zone region.

British Columbia, California, Oregon, Washington↗

USGS2024

Shallow faulting and folding beneath south‐central Seattle, Washington State, from land‐based high‐resolution seismic‐reflection imaging

The geologic framework of the Seattle fault zone (SFZ) has been extensively studied, but the structure and fault strand locations in the central portion of the fault zone through the city of Seattle have remained controversial. Much of what is known about the SFZ has come from light detection and ranging (lidar)‐topographic surveys and paleoseismic investigations of fault scarps primarily west of Puget Sound, regional gravity and aeromagnetic modeling, and multiscale marine seismic imaging in waters both west and east of Seattle. We analyze ∼24 km of land‐based P ‐wave seismic‐reflection data that fill in a critical gap in our understanding of the SFZ beneath the urban areas of West Seattle, south‐central Seattle, and Mercer Island. These data image deformed strata in the upper 1 km, including upwarped Tertiary rock and younger sediments. Collectively, these data provide evidence for multiple Quaternary‐active thrust faults, back thrusts, and sub‐basins within the SFZ beneath the city of Seattle. The results indicate that multiple and potentially active back thrusts in the upper ∼500 m extend across the length of the SFZ and the entire urban corridor that may be analogous to those on Bainbridge Island west of Puget Sound.

Washington↗

USGS2024

Empirical ground-motion basin response in the California Great Valley, Reno, Nevada, and Portland, Oregon

We assess how well the Next-Generation Attenuation-West 2 (NGA-West2) ground-motion models (GMMs), which are used in the US Geological Survey’s (USGS) National Seismic Hazard Model (NSHM) for crustal faults in the western United States, predict the observed basin response in the Great Valley of California, the Reno basin in Nevada, and Portland and Tualatin basins in Oregon. These GMMs rely on site parameters such as the time-averaged shear-wave velocity ( V S ) in the upper 30 m of Earth’s crust ( V S30 ) and depths to 1.0 and 2.5 km/s shear-wave isosurfaces ( Z 1.0 and Z 2.5 ) to capture basin effects and were developed using observations and simulations primarily from the Los Angeles region in southern California. Using ground-motion records from mostly small-to-moderate earthquakes and mixed-effects regression analysis, we find that the GMMs perform well with our local basin-depth models for the California Great Valley. With our local basin-depth models for Reno, the GMMs do not perform as well for this relatively shallow basin and exhibit little sensitivity to the basin parameters used in the NGA-West2 GMMs. We also find good performance for the local Z 1.0 model across the Portland region, whereas the local Z 2.5 model provides little predictive power except at sites in the deepest part of the Tualatin basin. Additional work could improve the performance of the site and basin terms in the NGA-West2 GMMs for regions with geologic structure different than the deep basins in southern California and the Great Valley. In addition, we find significant discrepancies among the GMMs in how the uncertainty in the ground motion varies with basin depth and pseudospectral period. Our results can help guide seismic hazard analyses on whether to include these local basin-depth models.

Earthquake Spectra↗

USGS2024

The 2023 US National Seismic Hazard Model: Ground-motion characterization for the conterminous United States

We update the ground-motion characterization for the 2023 National Seismic Hazard Model (NSHM) for the conterminous United States. The update includes the use of new ground-motion models (GMMs) in the Cascadia subduction zone; an adjustment to the central and eastern United States (CEUS) GMMs to reduce misfits with observed data; an updated boundary for the application of GMMs for shallow, crustal earthquakes in active tectonic regions (i.e. western United States (WUS)) and stable continental regions (i.e. CEUS); and the use of improved models for the site response of deep sedimentary basins in the WUS and CEUS. Site response updates include basin models for the California Great Valley and for the Portland and Tualatin basins, Oregon, as well as long-period basin effects from three-dimensional simulations in the Greater Los Angeles region and in the Seattle basin; in the CEUS, we introduce a broadband (0.01- to 10-s period) amplification model for the effects of the passive-margin basins of the Atlantic and Gulf Coastal Plains. In addition, we summarize progress on implementing rupture directivity models into seismic hazard models, although they are not incorporated in the 2023 NSHM. We implement the ground-motion characterization for the 2023 NSHM in the US Geological Survey’s code for probabilistic seismic hazard analysis, nshmp-haz-v2 , and present the sensitivity of hazard to these changes. Hazard calculations indicate widespread effects from adjustments to the CEUS GMMs, from the incorporation of Coastal Plain amplification effects, and from the treatment of shallow-basin and out-of-basin sites in the San Francisco Bay Area and Los Angeles region, as well as locally important changes from subduction-zone GMMs, and from updated and new WUS basins.

conterminous United States↗

USGS2023

The 2023 US 50-State National Seismic Hazard Model: Overview and implications

The US National Seismic Hazard Model (NSHM) was updated in 2023 for all 50 states using new science on seismicity, fault ruptures, ground motions, and probabilistic techniques to produce a standard of practice for public policy and other engineering applications (defined for return periods greater than ∼475 or less than ∼10,000 years). Changes in 2023 time-independent seismic hazard (both increases and decreases compared to previous NSHMs) are substantial because the new model considers more data and updated earthquake rupture forecasts and ground-motion components. In developing the 2023 model, we tried to apply best available or applicable science based on advice of co-authors, more than 50 reviewers, and hundreds of hazard scientists and end-users, who attended public workshops and provided technical inputs. The hazard assessment incorporates new catalogs, declustering algorithms, gridded seismicity models, magnitude-scaling equations, fault-based structural and deformation models, multi-fault earthquake rupture forecast models, semi-empirical and simulation-based ground-motion models, and site amplification models conditioned on shear-wave velocities of the upper 30 m of soil and deeper sedimentary basin structures. Seismic hazard calculations yield hazard curves at hundreds of thousands of sites, ground-motion maps, uniform-hazard response spectra, and disaggregations developed for pseudo-spectral accelerations at 21 oscillator periods and two peak parameters, Modified Mercalli Intensity, and 8 site classes required by building codes and other public policy applications. Tests show the new model is consistent with past ShakeMap intensity observations. Sensitivity and uncertainty assessments ensure resulting ground motions are compatible with known hazard information and highlight the range and causes of variability in ground motions. We produce several impact products including building seismic design criteria, intensity maps, planning scenarios, and engineering risk assessments showing the potential physical and social impacts. These applications provide a basis for assessing, planning, and mitigating the effects of future earthquakes.

Earthquake Spectra↗

USGS2022

Introduction to the special issue of the Consortium of Organizations for Strong Motion Observation Systems (COSMOS) international guidelines for applying noninvasive geophysical techniques to characterize seismic site conditions

Knowledge about local seismic site conditions provides critical information to account for site effects that are commonly observed in strong motion recordings. Certainly, other wave propagation effects can influence these observations, which are attributable to variations in material properties of the paths traveled by the waves, as well as the characteristics of the seismic source. However, local geologic conditions, particularly, when under shear-wave excitation, are known to have a strong influence on the behavior of ground shaking in the frequency range that is expected to directly affect the built environment. Thus, shear waves traveling in the shallow subsurface—defined here as tens to hundreds of meters beneath the ground surface—are the main foci for application and research in the earthquake engineering community.

Journal of Seismology↗

USGS2022

Flexible multimethod approach for seismic site characterization

We describe the flexible multimethod seismic site characterization technique for obtaining shear-wave velocity ( V S ) profiles and derivative information, such as the time-averaged V S of the upper 30 m ( V S30 ). Simply stated, the multimethod approach relies on the application of multiple independent noninvasive site characterization acquisition and analysis techniques utilized in a flexible field-based approach through a decision process based on primary factors such as the investigator’s available equipment, field logistics (e.g., available array deployment space, urban versus rural), geologic setting, and funding among other primary factors. In a multimethod approach, sites ideally should be characterized using both active and passive noninvasive (i.e., no drilling and only minimal disturbance to the ground surface) methods. Almost without exception, we recommend the use of at least one active-source technique for analyzing surface waves, which in the current state of practice would include one or more of the following: spectral analysis of surface waves (SASW; commonly Rayleigh waves) and multichannel analysis of surface waves (MASW; Rayleigh and/or Love waves). In addition, passive-source microtremor array methods (MAMs) using linear (one-dimensional) and two-dimensional arrays may be appropriate or even required for characterizing a given site. Similarly, the microtremor horizontal-to-vertical spectral ratio (mHVSR) method can be valuable for identifying shallow rock sites. Finally, P- and SH-wave seismic refraction methods warrant consideration, especially at rock and shallow-rock sites.

Journal of Seismology↗

USGS2022

Microtremor array method using spatial autocorrelation analysis of Rayleigh‑wave data

Microtremor array measurements (MAM) and passive surface wave methods in general, have been increasingly used to non-invasively estimate shear-wave velocity structures (Vs) for various purposes. The methods estimate dispersion curves and invert them for retrieving S-wave velocity profiles. This paper summarizes principles, limitations, data collection and processing methods. It intends to enable students and practitioners to understand the principles needed to plan a microtremor array investigation, record and process the data, and evaluate the quality of investigation result. The paper focuses on the spatial autocorrelation (SPAC) processing method among microtremor array processing methods because of its relatively simple calculation and stable applicability.

Journal of Seismology↗

USGS2022

Seismic site characterization with shear wave (SH) reflection and refraction methods

Reflection and critically refracted seismic methods use traveltime measurements of body waves propagating between a source and a series of receivers on the ground surface to calculate subsurface velocities. Body wave energy is refracted or reflected at boundaries where there is a change in seismic impedance, defined as the product of material density and seismic velocity. This article provides practical guidance on the use of horizontally propagating shear wave (SH-wave) refraction and reflection methods to determine shear wave velocity as a function of depth for near-surface seismic site characterizations. Method principles and the current state of engineering practice are reviewed, along with discussions of limitations and uncertainty assessments. Typical data collection procedures are described using basic survey equipment, along with information on more advanced applications and emerging technologies. Eight case studies provide examples of the techniques in real-world seismic site characterizations performed in a variety of geological settings.

Journal of Seismology↗

USGS2022

An assessment of uncertainties in VS profiles obtained from microtremor observations in the phased 2018 COSMOS blind trials

Site response is a critical consideration when assessing earthquake hazards. Site characterization is key to understanding site effects as influenced by seismic site conditions of the local geology. Thus, a number of geophysical site characterization methods were developed to meet the demand for accurate and cost-effective results. As a consequence, a number of studies have been administered periodically as blind trials to evaluate the state-of-practice on-site characterization. We present results from the Consortium of Organizations for Strong Motion Observation Systems (COSMOS) blind trials, which used data recorded from surface-based microtremor array methods (MAM) at four sites where geomorphic conditions vary from deep alluvial basins to an alpine valley. Thirty-four invited analysts participated. Data were incrementally released to 17 available analysts who participated in all four phases: (1) two-station arrays, (2) sparse triangular arrays, (3) complex nested triangular or circular arrays, and (4) all available geological control site information including drill hole data. Another set of 17 analysts provided results from two sites and two phases only. Although data from one site consisted of recordings from three-component sensors, the other three sites consisted of data recorded only by vertical-component sensors. The sites cover a range of noise source distributions, ranging from one site with a highly directional microtremor wave field to others with omni-directional (azimuthally distributed) wave fields. We review results from different processing techniques (e.g., beam-forming, spatial autocorrelation, cross-correlation, or seismic interferometry) applied by the analysts and compare the effectiveness between the differing wave field distributions. We define the M index as a quality index based on estimates of the time-averaged shear-wave velocity of the upper 10 ( V S10 ), 30 ( V S30 ), 100 ( V S100 ), and 300 ( V S300 ) meters and show its usefulness in quantitative comparisons of V S profiles from multiple analysts. Our findings are expected to aid in building an evidence-based consensus on preferred cost-effective arrays and processing methodology for future studies of seismic site effects.

Seismological Research Letters↗

USGS2022

Implementation of basin models and sediment depth terms in the 2023 update of the U.S. National Seismic Hazard Model: Example from Reno, Nevada

We present a framework to evaluate the inclusion of candidate basin depth models in the U.S. Geological Survey National Seismic Hazard Model. We compute intensity measures (peak and spectral amplitudes) from uniformly processed earthquake ground motions in and around the basin of interest and compare these to ground-motion model (GMM) estimates over a range of oscillator periods. The GMMs use depth to specific shear-wave velocity isosurfaces (Z x ) as a proxy for basin depth. We quantify whether the GMM estimates using Z x from the candidate basin depth model outperform the default estimates based on V S30 (the time-averaged shear wave velocity in the upper 30 m). We partition GMM residuals into event and site terms and compare site terms for stations within the basin to non-basin sites. We apply this framework to the greater Reno, Nevada, region, which has shallow basin depths (less than 450 m) and empirical amplifications of up to a factor of 2.2 at 1.0 s. There are no strong trends in site terms with basin depth, and the use of Z x values from the new candidate basin model shows marginal improvement of GMM total residuals compared to the default V S30 -based model (the condition of no relative basin amplification).

Nevada↗

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↗