No abstract

No abstract

For more information on the USGS-the Federal source for science about the Earth, its natural and living resources, natural hazards, and the environment-visit https://www.usgs.gov or call 1-888-ASK-USGS.For an overview of USGS information products, including maps, imagery, and publications, visit https://store.usgs.gov.Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.Although this information product, for the most part, is in the public domain, it also may contain copyrighted materials as noted in the text. Permission to reproduce copyrighted items must be secured from the copyright owner. Suggested citation:Stephenson, W.J., Reitman, N.G., and Angster, S.J., 2017, P-and S-wave velocity models incorporating the Cascadia subduction zone for 3D earthquake ground motion simulations, version 1. AcknowledgmentsThe V P and V S property volumes of model V1.6 were greatly improved by feedback from end-users including Art Frankel (U.S. Geological Survey [USGS]), Andy Delorey (Los Alamos National Laboratory), and John Vidale (University of Washington). We thank Jack Odum and Robert Williams (USGS) for their technical reviews, which greatly improved this manuscript. Special thanks to Jordan Bretthauer (USGS) for assistance in developing figure 1 of this report. Discussions with numerous collaborators and other interested parties seeking components of the model for their research helped prompt us to formally complete this updated documentation. This research was supported by funding from the USGS Earthquake Hazards Program. IntroductionIn support of earthquake hazards studies and ground motion simulations in the Pacific Northwest, threedimensional (3D) P-and S-wave velocity (V P and V S , respectively) models incorporating the Cascadia subduction zone were previously developed for the region encompassed from about 40.2°N. to 50°N. latitude, and from about 122°W. to 129°W. longitude ( fig. 1). This report describes updates to the Cascadia velocity property volumes of model version 1.3 ([V1.3]; Stephenson, 2007), herein called model version 1.6 (V1.6). As in model V1.3, the updated V1.6 model volume includes depths from 0 kilometers (km) (mean sea level) to 60 km, and it is intended to be a reference for researchers who have used, or are planning to use, this model in their earth science investigations. To this end, it is intended that the V P and V S property volumes of model V1.6 will be considered a template for a community velocity model of the Cascadia region as additional results become available. With the recent and ongoing development of the National Crustal Model (NCM; Boyd and Shah, 2016), we envision any future versions of this model will be directly integrated with that effort. BackgroundThe Cascadia subduction zone stretches for over 1,000 km, from the Mendocino Triple Junction off the northern California coast northward to Vancouver Island, Canada ( fig. 2). The primary reasons for developing these model volumes are (1) for simulat...

Slip distribution, slip rate, and slip per event for strike‐slip faults are commonly determined by correlating offset stream channels—under the assumption that they record seismic slip—but offset channels are formed by the interplay of tectonic and geomorphic processes. To constrain offset channel development under known tectonic and geomorphic conditions, we use numerical landscape evolution simulations along a theoretical strike‐slip fault with uniform and steady uplift, erosion, and diffusion. We investigate the influence of four tectonic parameters (fault zone width, earthquake recurrence interval, variance of the recurrence interval, and total slip relative to channel spacing) on offset channel development through multiple earthquake cycles. Analysis of >3,000 automatically measured offsets from >135 simulations suggests ~30% variability in individual measurements, but modeled displacement is recovered by averaging multiple measurements. However, the average of multiple offset measurements systematically underestimates modeled slip except when the fault zone is less than ~5 m wide, total slip is less than channel spacing, and offsets are measured shortly after an earthquake. In these simulations, postearthquake landscape evolution widens the geomorphic expression of the fault zone and modifies apparent channel offsets. We distinguish this “geomorphic fault zone” from the tectonic fault zone (zone of coseismic distributed deformation). This study highlights the capability of landscape evolution models to explore a range of conditions not easily defined in natural examples and the importance of averaging multiple measurements. Our results verify that paleoseismic studies must consider how geomorphic change has modified offset markers and use caution interpreting slip histories with multiple earthquakes.

Recent earthquakes involving multiple fault ruptures highlight the need to evaluate complex coastal deformation mechanisms, which are important for understanding plate boundary kinematics and seismic and tsunami hazards. We compare ages and uplift of the youngest Holocene marine terraces at Puatai Beach and Pakarae River mouth (∼10  km apart) in the northern Hikurangi subduction margin to examine whether uplift is the result of subduction earthquakes or upper-plate fault earthquakes. From stepped platform-cliff morphology, we infer uplift during 2–3 earthquakes and calculate an average uplift-per-event of 2.9±0.5  m at Puatai Beach and 2.0±0.5  m at Pakarae River mouth. Radiocarbon ages from the youngest beach deposit shells on each terrace and a tephra coverbed on one terrace constrain the timing of earthquakes to 1770–1710, 1100–910, and 420–250 cal. B.P. at Puatai Beach, and 1490–1290 and 660–530 cal. B.P. at Pakarae River mouth. The ages differ at each site indicating uplift is neither the result of subduction earthquakes nor single upper-plate fault earthquakes. A reinterpretation of new and existing bathymetry and seismic reflection data, combined with dislocation modeling, indicates that near-shore fault segmentation is more complex than previously thought and ruptures likely involve multiple upper-plate faults. Future updates of the New Zealand National Seismic Hazard Model should revise the northern Hikurangi subduction seismic sources so that rupture does not uplift Puatai Beach and Pakarae River mouth and include new near-shore upper-plate faults as multifault sources.