William J. Stephenson

Seismic reflection images beneath Puget Sound, western Washington State: The Puget Lowland thrust sheet hypothesis

Seismic reflection data show that the densely populated Puget Lowland of western Washington state is underlain by subhorizontal Paleogene and Neogene sedimentary rocks deformed by west and northwest trending faults and folds. From south to north beneath the Lowland, features seen on the seismic data include: the horizontally-stratified, 3.5 km thick Tacoma sedimentary basin; the Seattle uplift with south dipping (∼20°) strata on its south flank and steeply (50° to 90°) north dipping strata and the west-trending Seattle fault on its north flank; the 7.5 km thick, northward-thinning Seattle sedimentary basin; the antiformal Kingston arch; and the northwest trending, transpressional Southern Whidbey Island fault zone (SWIF). Interpreting the uplifts as fault-bend and fault-propagation folds leads to the hypothesis that the Puget Lowland lies on a north directed thrust sheet. The base of the thrust sheet may lie at 14 to 20 km depth within or at the base of a thick block of basaltic Crescent Formation; its edges may be right-lateral strike-slip faults along the base of the Cascade Range on the east and the Olympic Mountains on the west. Our model suggests that the Seattle fault has a long-term slip rate of about 0.25 mm/year and is large enough to generate a M7.6 to 7.7 earthquake.

Blind Shear-Wave Velocity Comparison of ReMi and MASW Results with Boreholes to 200 m in Santa Clara Valley: Implications for Earthquake Ground-Motion Assessment

Multichannel analysis of surface waves (MASW) and refraction micro- tremor (ReMi) are two of the most recently developed surface acquisition techniques for determining shallow shear-wave velocity. We conducted a blind comparison of MASW and ReMi results with four boreholes logged to at least 260 m for shear vel- ocity in Santa Clara Valley, California, to determine how closely these surface meth- ods match the downhole measurements. Average shear-wave velocity estimates to depths of 30, 50, and 100 m demonstrate that the surface methods as implemented in this study can generally match borehole results to within 15% to these depths. At two of the boreholes, the average to 100 m depth was within 3%. Spectral amplifi- cations predicted from the respective borehole velocity profiles similarly compare to within 15% or better from 1 to 10 Hz with both the MASW and ReMi surface-method velocity profiles. Overall, neither surface method was consistently better at matching the borehole velocity profiles or amplifications. Our results suggest MASW and ReMi surface acquisition methods can both be appropriate choices for estimating shear- wave velocity and can be complementary to each other in urban settings for hazards assessment.

Exploration and discovery in Yellowstone Lake: Results from high-resolution sonar imaging, seismic reflection profiling, and submersible studies

Abstract ‘No portion of the American continent is perhaps so rich in wonders as the Yellow Stone’ (F.V. Hayden, September 2, 1874) Discoveries from multi-beam sonar mapping and seismic reflection surveys of the northern, central, and West Thumb basins of Yellowstone Lake provide new insight into the extent of post-collapse volcanism and active hydrothermal processes occurring in a large lake environment above a large magma chamber. Yellowstone Lake has an irregular bottom covered with dozens of features directly related to hydrothermal, tectonic, volcanic, and sedimentary processes. Detailed bathymetric, seismic reflection, and magnetic evidence reveals that rhyolitic lava flows underlie much of Yellowstone Lake and exert fundamental control on lake bathymetry and localization of hydrothermal activity. Many previously unknown features have been identified and include over 250 hydrothermal vents, several very large (>500 m diameter) hydrothermal explosion craters, many small hydrothermal vent craters (∼1–200 m diameter), domed lacustrine sediments related to hydrothermal activity, elongate fissures cutting post-glacial sediments, siliceous hydrothermal spire structures, sublacustrine landslide deposits, submerged former shorelines, and a recently active graben. Sampling and observations with a submersible remotely operated vehicle confirm and extend our understanding of the identified features. Faults, fissures, hydrothermally inflated domal structures, hydrothermal explosion craters, and sublacustrine landslides constitute potentially significant geologic hazards. Toxic elements derived from hydrothermal processes also may significantly affect the Yellowstone ecosystem.

Three-Dimensional Simulations of Ground Motions in the Seattle Region for Earthquakes in the Seattle Fault Zone

We used the 3D finite-difference method to model observed seismograms of two earthquakes ( M L 4.9 and 3.5) in the Seattle region and to simulate ground motions for hypothetical M 6.5 and M 5.0 earthquakes on the Seattle fault, for periods greater than 2 sec. A 3D velocity model of the Seattle Basin was constructed from studies that analyzed seismic-reflection surveys, borehole logs, and gravity and aeromagnetic data. The observations and the simulations highlight the importance of the Seattle Basin on long-period ground motions. For earthquakes occurring just south of the basin, the edge of the basin and the variation of the thickness of the Quaternary deposits in the basin produce much larger surface waves than expected from flat-layered models. The data consist of seismograms recorded by instruments deployed in Seattle by the USGS and the University of Washington (UW). The 3D simulation reproduces the peak amplitude and duration of most of the seismograms of the June 1997 Bremerton event ( M L 4.9) recorded in Seattle. We found the focal mechanism for this event that best fits the observed seismograms in Seattle by combining Green9s functions determined from the 3D simulations for the six fundamental moment couples. The February 1997 event ( M L 3.5) to the south of the Seattle Basin exhibits a large surface-wave arrival at UW whose amplitude is matched by the synthetics in our 3D velocity model, for a source depth of 9 km. The M 6.5 simulations incorporated a fractal slip distribution on the fault plane. These simulations produced the largest ground motions in an area that includes downtown Seattle. This is mainly caused by rupture directed up dip toward downtown, radiation pattern of the source, and the turning of S waves by the velocity gradient in the Seattle basin. Another area of high ground motion is located about 13 km north of the fault and is caused by an increase in the amplitude of higher-mode Rayleigh waves caused by the thinning of the Quaternary deposits.

Simulation of Broadband Ground Motion Including Nonlinear Soil Effects for a Magnitude 6.5 Earthquake on the Seattle Fault, Seattle, Washington

The Seattle fault poses a significant seismic hazard to the city of Seattle, Washington. A hybrid, low-frequency, high-frequency method is used to calculate broadband (0–20 Hz) ground-motion time histories for a M 6.5 earthquake on the Seattle fault. Low frequencies ( 1 Hz) are calculated by a stochastic method that uses a fractal subevent size distribution to give an ω -2 displacement spectrum. Time histories are calculated for a grid of stations and then corrected for the local site response using a classification scheme based on the surficial geology. Average shear-wave velocity profiles are developed for six surficial geologic units: artificial fill, modified land, Esperance sand, Lawton clay, till, and Tertiary sandstone. These profiles together with other soil parameters are used to compare linear, equivalent-linear, and nonlinear predictions of ground motion in the frequency band 0–15 Hz. Linear site-response corrections are found to yield unreasonably large ground motions. Equivalent-linear and nonlinear calculations give peak values similar to the 1994 Northridge, California, earthquake and those predicted by regression relationships. Ground-motion variance is estimated for (1) randomization of the velocity profiles, (2) variation in source parameters, and (3) choice of nonlinear model. Within the limits of the models tested, the results are found to be most sensitive to the nonlinear model and soil parameters, notably the overconsolidation ratio.

Multiscale seismic imaging of active fault zones for hazard assessment; a case study of the Santa Monica fault zone, Los Angeles, California

High‐resolution seismic reflection profiles at two different scales were acquired across the transpressional Santa Monica Fault of north Los Angeles as part of an integrated hazard assessment of the fault. The seismic data confirm the location of the fault and related shallow faulting seen in a trench to deeper structures known from regional studies. The trench shows a series of near‐vertical strike‐slip faults beneath a topographic scarp inferred to be caused by thrusting on the Santa Monica fault. Analysis of the disruption of soil horizons in the trench indicates multiple earthquakes have occurred on these strike‐slip faults within the past 50 000 years, with the latest being 1000 to 3000 years ago. A 3.8-km-long, high‐resolution seismic reflection profile shows reflector truncations that constrain the shallow portion of the Santa Monica Fault (upper 300 m) to dip northward between 30° and 55°, most likely 30° to 35°, in contrast to the 60° to 70° dip interpreted for the deeper portion of the fault. Pr...

Near-surface structural model for deformation associated with the February 7, 1812, New Madrid, Missouri, earthquake

The February 7, 1812, New Madrid, Missouri, earthquake (M [moment magnitude] 8) was the third and final large-magnitude event to rock the northern Mississippi Embayment during the winter of 1811–1812. Although ground shaking was so strong that it rang church bells, stopped clocks, buckled pavement, and rocked buildings up and down the eastern seaboard, little coseismic surface deformation exists today in the New Madrid area. The fault(s) that ruptured during this event have remained enigmatic. We have integrated geomorphic data documenting differential surficial deformation (supplemented by historical accounts of surficial deformation and earthquake-induced Mississippi River waterfalls and rapids) with the interpretation of existing and recently acquired seismic reflection data, to develop a tectonic model of the near-surface structures in the New Madrid, Missouri, area. This model consists of two primary components: a north-northwest–trending thrust fault and a series of northeast-trending, strike-slip, tear faults. We conclude that the Reelfoot fault is a thrust fault that is at least 30 km long. We also infer that tear faults in the near surface partitioned the hanging wall into subparallel blocks that have undergone differential displacement during episodes of faulting. The northeast-trending tear faults bound an area documented to have been uplifted at least 0.5 m during the February 7, 1812, earthquake. These faults also appear to bound changes in the surface density of epicenters that are within the modern seismicity, which is occurring in the stepover zone of the left-stepping right-lateral strike-slip fault system of the modern New Madrid seismic zone.

Surface Seismic Measurements of Near‐Surface P‐ and S‐Wave Seismic Velocities at Earthquake Recording Stations, Seattle, Washington

We measured P- and S-wave seismic velocities to about 40-m depth using seismic-refraction/reflection data on the ground surface at 13 sites in the Seattle, Washington, urban area, where portable digital seismographs recently recorded earthquakes. Sites with the lowest measured V s correlate with highest ground motion amplification. These sites, such as at Harbor Island and in the Duwamish River industrial area (DRIA) south of the Kingdome, have an average V s in the upper 30 m (V¯ s30 ) of 150 to 170 m/s. These values of V¯ s30 place these sites in soil profile type E (V¯ s30 < 180 m/s). A “rock” site, located at Seward Park on Tertiary sedimentary deposits, has a V¯ s30 of 433 m/s, which is soil type C (V¯ s30 : 360 to 760 m/s). The Seward Park site V¯ s30 is about equal to, or up to 200 m/s slower than sites that were located on till or glacial outwash. High-amplitude P- and S-wave seismic reflections at several locations appear to correspond to strong resonances observed in earthquake spectra. An S-wave reflector at the Kingdome at about 17 to 22 m depth probably causes strong 2-Hz resonance that is observed in the earthquake data near the Kingdome.

An example of neotectonism in a continental interior - Thebes Gap, Midcontinent, United States

Abstract Some of the most intense neotectonic activity known in the continental interior of North America has been recently discovered on a fault zone in the Thebes Gap area, Missouri and Illinois. This faulting almost assuredly was accompanied by large earthquakes. The zone is located approximately 30 km north of the New Madrid seismic zone and consists of complex north-northeast- to northeast-striking, steeply dipping faults that have had a long-lived history of reactivation throughout most of the Phanerozoic. Geophysical studies by others suggest that the faults are rooted in the deeply buried Late Proterozoic and Early Cambrian Reelfoot rift system. Quaternary deposits are cut by at least four episodes of faulting, two of which occurred during the Holocene. The overall style of neotectonic deformation is interpreted as right-lateral strike-slip faulting. At many locations, however, near-surface displacements have stepped from one fault strand to another and produced normal and oblique-slip faults in areas of transtension and high-angle reverse faults, thrust faults, and folds in areas of transpression. There is evidence of reactivation of some near-surface fault segments during the great 1811–1812 New Madrid earthquakes. Quaternary faulting at Thebes Gap demonstrates that there are additional seismic-source zones in the Midcontinent, U.S., other than New Madrid, and that even in the absence of plate-margin orogenesis, intense neotectonic activity does occur over long time periods along crustal weaknesses in continental interiors.

Correlation of 1- to 10-Hz earthquake resonances with surface measurements of S-wave reflections and refractions in the upper 50 m

Resonances observed in earthquake seismograms recorded in Seattle, Washington, the central United States and Sherman Oaks, California, are correlated with each sites respective near-surface seismic velocity profile and reflectivity determined from shallow seismic-reflection/refraction surveys. In all of these cases the resonance accounts for the highest amplitude shaking at the site above 1 Hz. These results show that imaging near-surface reflections from the ground surface can locate impedance structures that are important contributors to earthquake ground shaking. A high-amplitude S -wave reflection, recorded 250-m northeast and 300-m east of the Seattle Kingdome earthquake-recording station, with a two-way travel time of about 0.23 to 0.27 sec (about 18- to 22-m depth) marks the boundary between overlying alluvium ( V S < 180 m/sec) and a higher velocity material ( V S about 400 m/sec). This reflector probably causes a strong 2-Hz resonance that is observed in the earthquake data for the site near the Kingdome. In the central United States, S -wave reflections from a high-impedance boundary (an S -wave velocity increase from about 200 m/sec to 2000 m/sec) at about 40-m depth corresponds to a strong fundamental resonance at about 1.5 Hz. In Sherman Oaks, strong resonances at about 1.0 and 4 Hz are consistently observed on earthquake seismograms. A strong S -wave reflector at about 40-m depth may cause the 1.0 Hz resonance. The 4.0-Hz resonance is possibly explained by constructive interference between the first overtone of the 1.0-Hz resonance and a 3.25- to 3.9-Hz resonance calculated from an areally consistent impedance boundary at about 10-m depth as determined by S -wave refraction data.