Search USGS⌕ Search

Geology topics

William J. Stephenson

Publications and source records attributed to William J. Stephenson.

At least 55 records · Page 3Linked to original sources

Preliminary assessment of a previously unknown fault zone beneath the Daytona Beach sand blow cluster near Marianna, Arkansas

We collected new high‐resolution P ‐wave seismic‐reflection data to explore for possible faults beneath a roughly linear cluster of early to mid‐Holocene earthquake‐induced sand blows to the south of Marianna, Arkansas. The Daytona Beach sand blow deposits are located in east‐central Arkansas about 75 km southwest of Memphis, Tennessee, and about 80 km south of the southwestern end of the New Madrid seismic zone (NMSZ). Previous studies of these sand blows indicate that they were produced between 10,500 and 5350 yr B.P. (before A.D. 1950). The sand blows are large and similar in size to those in the heart of the NMSZ produced by the 1811–1812 earthquakes. The seismic‐reflection profiles reveal a previously unknown zone of near‐vertical faults imaged in the 100–1100‐m depth range that are approximately coincident with a cluster of earthquake‐induced sand blows and a near‐linear surface lineament composed of air photo tonal anomalies. These interpreted faults are expressed as vertical discontinuities with the largest displacement fault showing about 40 m of west‐side‐up displacement at the top of the Paleozoic section at about 1100 m depth. There are about 20 m of folding on reflections within the Eocene strata at 400 m depth. Increasing fault displacement with depth suggests long‐term recurrent faulting. The imaged faults within the vicinity of the numerous sand blow features could be a causative earthquake source, although it does not rule out the possibility of other seismic sources nearby. These newly located faults add to a growing list of potentially active Pleistocene–Holocene faults discovered over the last two decades that are within the Mississippi embayment region but outside of the historical NMSZ.

Arkansas↗

Comparisons between vs30 and spectral response for 30 sites in Newcastle, Australia from collocated seismic cone penetrometer, active- and passive-source vs data

Although the time‐averaged shear‐wave velocity down to 30 m depth ( V S 30 ) can be a proxy for estimating earthquake ground‐motion amplification, significant controversy exists about its limitations when used as a single parameter for the prediction of amplification. To examine this question in absence of relevant strong‐motion records, we use a range of different methods to measure the shear‐wave velocity profiles and the resulting theoretical site amplification factors (AFs) for 30 sites in the Newcastle area, Australia, in a series of blind comparison studies. The multimethod approach used here combines past seismic cone penetrometer and spectral analysis of surface‐wave data, with newly acquired horizontal‐to‐vertical spectral ratio, passive‐source surface‐wave spatial autocorrelation (SPAC), refraction microtremor (ReMi), and multichannel analysis of surface‐wave data. The various measurement techniques predicted a range of different AFs. The SPAC and ReMi techniques have the smallest overall deviation from the median AF for the majority of sites. We show that V S 30 can be related to spectral response above a period T of 0.5 s but not necessarily with the maximum amplification according to the modeling done based on the measured shear‐wave velocity profiles. Both V S 30 and AF values are influenced by the velocity ratio between bedrock and overlying sediments and the presence of surficial thin low‐velocity layers (<2 m thick and <150 m/s), but the velocity ratio is what mostly affects the AF. At 0.2< T <0.4 s, the AFs are largely controlled by the surficial geology of a particular site. AF maxima are the highest in the hard classes, which is the inverse of the findings used in the Australian Building Code. Only for T >0.5 s do the amplification curves consistently show higher values for soft site classes and lower for hard classes.

Bulletin of the Seismological Society of America↗

Seismic imaging beneath an InSAR anomaly in eastern Washington State: Shallow faulting associated with an earthquake swarm in a low-hazard area

In 2001, a rare swarm of small, shallow earthquakes beneath the city of Spokane, Washington, caused ground shaking as well as audible booms over a five‐month period. Subsequent Interferometric Synthetic Aperture Radar (InSAR) data analysis revealed an area of surface uplift in the vicinity of the earthquake swarm. To investigate the potential faults that may have caused both the earthquakes and the topographic uplift, we collected ∼3 km of high‐resolution seismic‐reflection profiles to image the upper‐source region of the swarm. The two profiles reveal a complex deformational pattern within Quaternary alluvial, fluvial, and flood deposits, underlain by Tertiary basalts and basin sediments. At least 100 m of arching on a basalt surface in the upper 500 m is interpreted from both the seismic profiles and magnetic modeling. Two west‐dipping faults deform Quaternary sediments and project to the surface near the location of the Spokane fault defined from modeling of the InSAR data.

Washington↗

Compilation of V S 30 Data for the United States

V S 30 , the time-averaged shear-wave velocity ( V S ) to a depth of 30 meters, is a key index adopted by the earthquake engineering community to account for seismic site conditions. V S 30 is typically based on geophysical measurements of V S derived from invasive and noninvasive techniques at sites of interest. Owing to cost considerations, as well as logistical and environmental concerns, V S 30 data are sparse or not readily available for most areas. Where data are available, V S 30 values are often assembled in assorted formats that are accessible from disparate and (or) impermanent Web sites. To help remedy this situation, we compiled V S 30 measurements obtained by studies funded by the U.S. Geological Survey (USGS) and other governmental agencies. Thus far, we have compiled V S 30 values for 2,997 sites in the United States, along with metadata for each measurement from government-sponsored reports, Web sites, and scientific and engineering journals. Most of the data in our V S 30 compilation originated from publications directly reporting the work of field investigators. A small subset (less than 20 percent) of V S 30 values was previously compiled by the USGS and other research institutions. Whenever possible, V S 30 originating from these earlier compilations were crosschecked against published reports. Both downhole and surface-based V S 30 estimates are represented in our V S 30 compilation. Most of the V S 30 data are for sites in the western contiguous United States (2,141 sites), whereas 786 V S 30 values are for sites in the Central and Eastern United States; 70 values are for sites in other parts of the United States, including Alaska (15 sites), Hawaii (30 sites), and Puerto Rico (25 sites). An interactive map is hosted on the primary USGS Web site for accessing V S 30 data ( https://earthquake.usgs.gov/data/vs30/us/ ).

Data Series↗

Kinematics of shallow backthrusts in the Seattle fault zone, Washington State

Near-surface thrust fault splays and antithetic backthrusts at the tips of major thrust fault systems can distribute slip across multiple shallow fault strands, complicating earthquake hazard analyses based on studies of surface faulting. The shallow expression of the fault strands forming the Seattle fault zone of Washington State shows the structural relationships and interactions between such fault strands. Paleoseismic studies document an ∼7000 yr history of earthquakes on multiple faults within the Seattle fault zone, with some backthrusts inferred to rupture in small (M ∼5.5–6.0) earthquakes at times other than during earthquakes on the main thrust faults. We interpret seismic-reflection profiles to show three main thrust faults, one of which is a blind thrust fault directly beneath downtown Seattle, and four small backthrusts within the Seattle fault zone. We then model fault slip, constrained by shallow deformation, to show that the Seattle fault forms a fault propagation fold rather than the alternatively proposed roof thrust system. Fault slip modeling shows that back-thrust ruptures driven by moderate (M ∼6.5–6.7) earthquakes on the main thrust faults are consistent with the paleoseismic data. The results indicate that paleoseismic data from the back-thrust ruptures reveal the times of moderate earthquakes on the main fault system, rather than indicating smaller (M ∼5.5–6.0) earthquakes involving only the backthrusts. Estimates of cumulative shortening during known Seattle fault zone earthquakes support the inference that the Seattle fault has been the major seismic hazard in the northern Cascadia forearc in the late Holocene.

Washington↗

The use of wavenumber normalization in computing spatially averaged coherencies (KRSPAC) of microtremor data from asymmetric arrays

The SPAC method of processing microtremor noise observations for estimation of Vs profiles has a limitation that the array has circular or triangular symmetry in order to allow spatial (azimuthal) averaging of inter-station coherencies over a constant station separation. Common processing methods allow for station separations to vary by typically &plusmn;10% in the azimuthal averaging before degradation of the SPAC spectrum is excessive. A limitation on use of high-wavenumbers in inversions of SPAC spectra to Vs profiles has been the requirement for exact array symmetry to avoid loss of information in the azimuthal averaging step. In this paper we develop a new wavenumber-normalised SPAC method (KRSPAC) where instead of performing averaging of sets of coherency versus frequency spectra and then fitting to a model SPAC spectrum, we interpolate each spectrum to coherency versus k.r, where k and r are wavenumber and station separation respectively, and r may be different for each pair of stations. For fundamental mode Rayleigh-wave energy the model SPAC spectrum to be fitted reduces to Jo(kr). The normalization process changes with each iteration since k is a function of frequency and phase velocity and hence is updated each iteration. The method proves robust and is demonstrated on data acquired in the Santa Clara Valley, CA, (Site STGA) where an asymmetric array having station separations varying by a factor of 2 is compared with a conventional triangular array; a 300-mdeep borehole with a downhole Vs log provides nearby ground truth. The method is also demonstrated on data from the Pleasanton array, CA, where station spacings are irregular and vary from 400 to 1200 m. The KRSPAC method allows inversion of data using kr (unitless) values routinely up to 30, and occasionally up to 60. Thus despite the large and irregular station spacings, this array permits resolution of Vs as fine as 15 m for the near-surface sediments, and down to a maximum depth of 2.5 km.

Conference Paper↗

Modeling the effects of source and path heterogeneity on ground motions of great earthquakes on the Cascadia Subduction Zone Using 3D simulations

We ran finite‐difference earthquake simulations for great subduction zone earthquakes in Cascadia to model the effects of source and path heterogeneity for the purpose of improving strong‐motion predictions. We developed a rupture model for large subduction zone earthquakes based on a k −2 slip spectrum and scale‐dependent rise times by representing the slip distribution as the sum of normal modes of a vibrating membrane. Finite source and path effects were important in determining the distribution of strong motions through the locations of the hypocenter, subevents, and crustal structures like sedimentary basins. Some regions in Cascadia appear to be at greater risk than others during an event due to the geometry of the Cascadia fault zone relative to the coast and populated regions. The southern Oregon coast appears to have increased risk because it is closer to the locked zone of the Cascadia fault than other coastal areas and is also in the path of directivity amplification from any rupture propagating north to south in that part of the subduction zone, and the basins in the Puget Sound area are efficiently amplified by both north and south propagating ruptures off the coast of western Washington. We find that the median spectral accelerations at 5 s period from the simulations are similar to that of the Zhao et al. (2006) ground‐motion prediction equation, although our simulations predict higher amplitudes near the region of greatest slip and in the sedimentary basins, such as the Seattle basin.

Bulletin of the Seismological Society of America↗

Ground-motion site effects from multimethod shear-wave velocity characterization at 16 seismograph stations deployed for aftershocks of the August 2011 Mineral, Virginia earthquake

We characterize shear-wave velocity versus depth (Vs profile) at 16 portable seismograph sites through the epicentral region of the 2011 M w 5.8 Mineral (Virginia, USA) earthquake to investigate ground-motion site effects in the area. We used a multimethod acquisition and analysis approach, where active-source horizontal shear (SH) wave reflection and refraction as well as active-source multichannel analysis of surface waves (MASW) and passive-source refraction microtremor (ReMi) Rayleigh wave dispersion were interpreted separately. The time-averaged shear-wave velocity to a depth of 30 m (Vs30), interpreted bedrock depth, and site resonant frequency were estimated from the best-fit Vs profile of each method at each location for analysis. Using the median Vs30 value (270&ndash;715 m/s) as representative of a given site, we estimate that all 16 sites are National Earthquake Hazards Reduction Program (NEHRP) site class C or D. Based on a comparison of simplified mapped surface geology to median Vs30 at our sites, we do not see clear evidence for using surface geologic units as a proxy for Vs30 in the epicentral region, although this may primarily be because the units are similar in age (Paleozoic) and may have similar bulk seismic properties. We compare resonant frequencies calculated from ambient noise horizontal:vertical spectral ratios (HVSR) at available sites to predicted site frequencies (generally between 1.9 and 7.6 Hz) derived from the median bedrock depth and average Vs to bedrock. Robust linear regression of HVSR to both site frequency and Vs30 demonstrate moderate correlation to each, and thus both appear to be generally representative of site response in this region. Based on Kendall tau rank correlation testing, we find that Vs30 and the site frequency calculated from average Vs to median interpreted bedrock depth can both be considered reliable predictors of weak-motion site effects in the epicentral region.

Virginia↗

Site response in the eastern United States: A comparison of Vs30 measurements with estimates from horizontal:vertical spectral ratios

Earthquake damage is often increased due to local ground-motion amplification caused by soft soils, thick basin sediments, topographic effects, and liquefaction. A critical factor contributing to the assessment of seismic hazard is detailed information on local site response. In order to address and quantify the site response at seismograph stations in the eastern United States, we investigate the regional spatial variation of horizontal:vertical spectral ratios (HVSR) using ambient noise recorded at permanent regional and national network stations as well as temporary seismic stations deployed in order to record aftershocks of the 2011 Mineral, Virginia, earthquake. We compare the HVSR peak frequency to surface measurements of the shear-wave seismic velocity to 30 m depth (Vs30) at 21 seismograph stations in the eastern United States and find that HVSR peak frequency increases with increasing Vs30. We use this relationship to estimate the National Earthquake Hazards Reduction Program soil class at 218 ANSS (Advanced National Seismic System), GSN (Global Seismographic Network), and RSN (Regional Seismograph Networks) locations in the eastern United States, and suggest that this seismic station–based HVSR proxy could potentially be used to calibrate other site response characterization methods commonly used to estimate shaking hazard.

GSA Special Papers↗

Quaternary extensional growth folding beneath Reno, Nevada, imaged by urban seismic profiling

We characterize shallow subsurface faulting and basin structure along a transect through heavily urbanized Reno, Nevada, with high‐resolution seismic reflection imaging. The 6.8 km of P‐wave data image the subsurface to approximately 800 m depth and delineate two subbasins and basin uplift that are consistent with structure previously inferred from gravity modeling in this region of the northern Walker Lane. We interpret two primary faults that bound the uplift and deform Quaternary deposits. The dip of Quaternary and Tertiary strata in the western subbasin increases with greater depth to the east, suggesting recurrent fault motion across the westernmost of these faults. Deformation in the Quaternary section of the western subbasin is likely evidence of extensional growth folding at the edge of the Truckee River through Reno. This deformation is north of, and on trend with, previously mapped Quaternary fault strands of the Mt. Rose fault zone. In addition to corroborating the existence of previously inferred intrabasin structure, these data provide evidence for an active extensional Quaternary fault at a previously unknown location within the Truckee Meadows basin that furthers our understanding of both the seismotectonic framework and earthquake hazards in this urbanized region.

Nevada↗

Vs30 and spectral response from collocated shallow, active- and passive-source Vs data at 27 sites in Puerto Rico

Shear‐wave velocity (V S ) and time‐averaged shear‐wave velocity to 30 m depth (V S30 ) are the key parameters used in seismic site response modeling and earthquake engineering design. Where V S data are limited, available data are often used to develop and refine map‐based proxy models of V S30 for predicting ground‐motion intensities. In this paper, we present shallow V S data from 27 sites in Puerto Rico. These data were acquired using a multimethod acquisition approach consisting of noninvasive, collocated, active‐source body‐wave (refraction/reflection), active‐source surface wave at nine sites, and passive‐source surface‐wave refraction microtremor (ReMi) techniques. V S ‐versus‐depth models are constructed and used to calculate spectral response plots for each site. Factors affecting method reliability are analyzed with respect to site‐specific differences in bedrock V S and spectral response. At many but not all sites, body‐ and surface‐wave methods generally determine similar depths to bedrock, and it is the difference in bedrock V S that influences site amplification. The predicted resonant frequencies for the majority of the sites are observed to be within a relatively narrow bandwidth of 1–3.5 Hz. For a first‐order comparison of peak frequency position, predictive spectral response plots from eight sites are plotted along with seismograph instrument spectra derived from the time series of the 16 May 2010 Puerto Rico earthquake. We show how a multimethod acquisition approach using collocated arrays compliments and corroborates V S results, thus adding confidence that reliable site characterization information has been obtained.

Bulletin of the Seismological Society of America↗

Origin of the Blytheville Arch, and long-term displacement on the New Madrid seismic zone, central United States

The southern arm of the New Madrid seismic zone of the central United States coincides with the buried, ~110 km by ~20 km Blytheville Arch antiform within the Cambrian–Ordovician Reelfoot rift graben. The Blytheville Arch has been interpreted at various times as a compressive structure, an igneous intrusion, or a sediment diapir. Reprocessed industry seismic-reflection profiles presented here show a strong similarity between the Blytheville Arch and pop-up structures, or flower structures, within strike-slip fault systems. The Blytheville Arch formed in the Paleozoic, but post–Mid-Cretaceous to Quaternary strata show displacement or folding indicative of faulting. Faults within the graben structure but outside of the Blytheville Arch also appear to displace Upper Cretaceous and perhaps younger strata, indicating that past faulting was not restricted to the Blytheville Arch and New Madrid seismic zone. As much as 10–12.5 km of strike slip can be estimated from apparent shearing of the Reelfoot arm of the New Madrid seismic zone. There also appears to be ~5–5.5 km of shearing of the Reelfoot topographic scarp at the north end of the southern arm of the New Madrid seismic zone and of the southern portion of Crowley's Ridge, which is a north-trending topographic ridge just south of the seismic zone. These observations suggest that there has been substantial strike-slip displacement along the Blytheville Arch and southern arm of the New Madrid seismic zone, that strike-slip extended north and south of the modern seismic zone, and that post–Mid-Cretaceous (post-Eocene?) faulting was not restricted to the Blytheville Arch or to currently active faults within the New Madrid seismic zone.

New Madrid seismic zone↗

Characterization of intrabasin faulting and deformation for earthquake hazards in southern Utah Valley, Utah, from high-resolution seismic imaging

We conducted active and passive seismic imaging investigations along a 5.6-km-long, east&ndash;west transect ending at the mapped trace of the Wasatch fault in southern Utah Valley. Using two-dimensional (2D) P-wave seismic reflection data, we imaged basin deformation and faulting to a depth of 1.4 km and developed a detailed interval velocity model for prestack depth migration and 2D ground-motion simulations. Passive-source microtremor data acquired at two sites along the seismic reflection transect resolve S-wave velocities of approximately 200 m/s at the surface to about 900 m/s at 160 m depth and confirm a substantial thickening of low-velocity material westward into the valley. From the P-wave reflection profile, we interpret shallow (100&ndash;600 m) bedrock deformation extending from the surface trace of the Wasatch fault to roughly 1.5 km west into the valley. The bedrock deformation is caused by multiple interpreted fault splays displacing fault blocks downward to the west of the range front. Further west in the valley, the P-wave data reveal subhorizontal horizons from approximately 90 to 900 m depth that vary in thickness and whose dip increases with depth eastward toward the Wasatch fault. Another inferred fault about 4 km west of the mapped Wasatch fault displaces horizons within the valley to as shallow as 100 m depth. The overall deformational pattern imaged in our data is consistent with the Wasatch fault migrating eastward through time and with the abandonment of earlier synextensional faults, as part of the evolution of an inferred 20-km-wide half-graben structure within Utah Valley. Finite-difference 2D modeling suggests the imaged subsurface basin geometry can cause fourfold variation in peak ground velocity over distances of 300 m.

Utah↗

Recent faulting in western Nevada revealed by multi-scale seismic reflection

The main goal of this study is to compare different reflection methods used to image subsurface structure within different physical environments in western Nevada. With all the methods employed, the primary goal is fault imaging for structural information toward geothermal exploration and seismic hazard estimation. We use seismic CHIRP (a swept-frequency marine acquisition system), weight drop (an accelerated hammer source), and two different vibroseis systems to characterize fault structure. We focused our efforts in the Reno metropolitan area and the area within and surrounding Pyramid Lake in northern Nevada. These different methods have provided valuable constraints on the fault geometry and activity, as well as associated fluid movement. These are critical in evaluating the potential for large earthquakes in these areas, and geothermal exploration possibilities near these structures.

Nevada↗

High-resolution seismic reflection imaging of growth folding and shallow faults beneath the Southern Puget Lowland, Washington State

Marine seismic reflection data from southern Puget Sound, Washington, were collected to investigate the nature of shallow structures associated with the Tacoma fault zone and the Olympia structure. Growth folding and probable Holocene surface deformation were imaged within the Tacoma fault zone beneath Case and Carr Inlets. Shallow faults near potential field anomalies associated with the Olympia structure were imaged beneath Budd and Eld Inlets. Beneath Case Inlet, the Tacoma fault zone includes an &sim;350-m wide section of south-dipping strata forming the upper part of a fold (kink band) coincident with the southern edge of an uplifted shoreline terrace. An &sim;2 m change in the depth of the water bottom, onlapping postglacial sediments, and increasing stratal dips with increasing depth are consistent with late Pleistocene to Holocene postglacial growth folding above a blind fault. Geologic data across a topographic lineament on nearby land indicate recent uplift of late Holocene age. Profiles acquired in Carr Inlet 10 km to the east of Case Inlet showed late Pleistocene or Holocene faulting at one location with &sim;3 to 4 m of vertical displacement, south side up. North of this fault the data show several other disruptions and reflector terminations that could mark faults within the broad Tacoma fault zone. Seismic reflection profiles across part of the Olympia structure beneath southern Puget Sound show two apparent faults about 160 m apart having 1 to 2 m of displacement of subhorizontal bedding. Directly beneath one of these faults, a dipping reflector that may mark the base of a glacial channel shows the opposite sense of throw, suggesting strike-slip motion. Deeper seismic reflection profiles show disrupted strata beneath these faults but little apparent vertical offset, consistent with strike-slip faulting. These faults and folds indicate that the Tacoma fault and Olympia structure include active structures with probable postglacial motion.

Washington↗

The Quaternary Silver Creek Fault Beneath the Santa Clara Valley, California

The northwest-trending Silver Creek Fault is a 40-km-long strike-slip fault in the eastern Santa Clara Valley, California, that has exhibited different behaviors within a changing San Andreas Fault system over the past 10-15 Ma. Quaternary alluvium several hundred meters thick that buries the northern half of the Silver Creek Fault, and that has been sampled by drilling and imaged in a detailed seismic reflection profile, provides a record of the Quaternary history of the fault. We assemble evidence from areal geology, stratigraphy, paleomagnetics, ground-water hydrology, potential-field geophysics, and reflection and earthquake seismology to determine the long history of the fault in order to evaluate its current behavior. The fault formed in the Miocene more than 100 km to the southeast, as the southwestern fault in a 5-km-wide right step to the Hayward Fault, within which the 40-km-long Evergreen pull-apart basin formed. Later, this basin was obliquely cut by the newly recognized Mt. Misery Fault to form a more direct connection to the Hayward Fault, although continued growth of the basin was sufficient to accommodate at least some late Pliocene alluvium. Large offset along the San Andreas-Calaveras-Mt Misery-Hayward Faults carried the basin northwestward almost to its present position when, about 2 Ma, the fault system was reorganized. This led to near abandonment of the faults bounding the pull-apart basin in favor of right slip extending the Calaveras Fault farther north before stepping west to the Hayward Fault, as it does today. Despite these changes, the Silver Creek Fault experienced a further 200 m of dip slip in the early Quaternary, from which we infer an associated 1.6 km or so of right slip, based on the ratio of the 40-km length of the strike-slip fault to a 5-km depth of the Evergreen Basin. This dip slip ends at a mid-Quaternary unconformity, above which the upper 300 m of alluvial cover exhibits a structural sag at the fault that we interpret as a negative flower structure. This structure implies some continuing strike slip on the Silver Creek Fault in the late Quaternary as well, with a transtensional component but no dip slip. Our only basis for estimating the rate of this later Quaternary strike slip on the Silver Creek Fault is to assume continuation of the inferred early Quaternary rate of less than 2 mm/yr. Faulting evident in a detailed seismic reflection profile across the Silver Creek Fault extends up to the limit of data at a depth of 50 m and age of about 140 ka, and the course of Coyote Creek suggests Holocene capture in a structural depression along the fault. No surface trace is evident on the alluvial plain, however, and convincing evidence of Holocene offset is lacking. Few instrumentally recorded earthquakes are located near the fault, and those that are near its southern end represent cross-fault shortening, not strike slip. The fault might have been responsible, however, for two poorly located moderate earthquakes that occurred in the area in 1903. Its southeastern end does mark an abrupt change in the pattern of abundant instrumentally recorded earthquakes along the Calaveras Fault-in both its strike and in the depth distribution of hypocenters-that could indicate continuing influence by the Silver Creek Fault. In the absence of convincing evidence to the contrary, and as a conservative estimate, we presume that the Silver Creek Fault has continued its strike-slip movement through the Holocene, but at a very slow rate. Such a slow rate would, at most, yield very infrequent damaging earthquakes. If the 1903 earthquakes did, in fact, occur on the Silver Creek Fault, they would have greatly reduced the short-term future potential for large earthquakes on the fault.

Open-File Report↗

Near-Surface Shear Wave Velocity Versus Depth Profiles, VS30, and NEHRP Classifications for 27 Sites in Puerto Rico

In 2004 and 2005 the Puerto Rico Seismic Network (PRSN), Puerto Rico Strong Motion Program (PRSMP) and the Geology Department at the University of Puerto Rico-Mayaguez (UPRM) collaborated with the U.S. Geological Survey to study near-surface shear-wave (Vs) and compressional-wave (Vp) velocities in and around major urban areas of Puerto Rico. Using noninvasive seismic refraction-reflection profiling techniques, we acquired velocities at 27 locations. Surveyed sites were predominantly selected on the premise that they were generally representative of near-surface materials associated with the primary geologic units located within the urbanized areas of Puerto Rico. Geologic units surveyed included Cretaceous intrusive and volcaniclastic bedrock, Tertiary sedimentary and volcanic units, and Quaternary unconsolidated eolian, fluvial, beach, and lagoon deposits. From the data we developed Vs and Vp depth versus velocity columns, calculated average Vs to 30-m depth (VS30), and derived NEHRP (National Earthquake Hazards Reduction Program) site classifications for all sites except one where results did not reach 30-m depth. The distribution of estimated NEHRP classes is as follows: three class 'E' (VS30 below 180 m/s), nine class 'D' (VS30 between 180 and 360 m/s), ten class 'C' (VS30 between 360 and 760 m/s), and four class 'B' (VS30 greater than 760 m/s). Results are being used to calibrate site response at seismograph stations and in the development of regional and local shakemap models for Puerto Rico.

Open-File Report↗

Seismic Hazard Maps for Seattle, Washington, Incorporating 3D Sedimentary Basin Effects, Nonlinear Site Response, and Rupture Directivity

This report presents probabilistic seismic hazard maps for Seattle, Washington, based on over 500 3D simulations of ground motions from scenario earthquakes. These maps include 3D sedimentary basin effects and rupture directivity. Nonlinear site response for soft-soil sites of fill and alluvium was also applied in the maps. The report describes the methodology for incorporating source and site dependent amplification factors into a probabilistic seismic hazard calculation. 3D simulations were conducted for the various earthquake sources that can affect Seattle: Seattle fault zone, Cascadia subduction zone, South Whidbey Island fault, and background shallow and deep earthquakes. The maps presented in this document used essentially the same set of faults and distributed-earthquake sources as in the 2002 national seismic hazard maps. The 3D velocity model utilized in the simulations was validated by modeling the amplitudes and waveforms of observed seismograms from five earthquakes in the region, including the 2001 M6.8 Nisqually earthquake. The probabilistic seismic hazard maps presented here depict 1 Hz response spectral accelerations with 10%, 5%, and 2% probabilities of exceedance in 50 years. The maps are based on determinations of seismic hazard for 7236 sites with a spacing of 280 m. The maps show that the most hazardous locations for this frequency band (around 1 Hz) are soft-soil sites (fill and alluvium) within the Seattle basin and along the inferred trace of the frontal fault of the Seattle fault zone. The next highest hazard is typically found for soft-soil sites in the Duwamish Valley south of the Seattle basin. In general, stiff-soil sites in the Seattle basin exhibit higher hazard than stiff-soil sites outside the basin. Sites with shallow bedrock outside the Seattle basin have the lowest estimated hazard for this frequency band.

Open-File Report↗