Search USGS⌕ Search

Geology topics

Nasser A. Marafi

Publications and source records attributed to Nasser A. Marafi.

10 recordsLinked to original sources
USGS2021

A generic soil velocity model that accounts for near-surface conditions and deeper geologic structure

Near-surface soil conditions can significantly alter the amplitude and frequency content of incoming ground motions – often with profound consequences for the built environment – and are thus important inputs to any ground-motion prediction. Previous soil-velocity models (SVM) have predicted shear-wave velocity profiles based on the time-averaged shear-wave velocity in the upper 30 m ( V S30 ). This article presents a generic soil-velocity model that accounts both for near-surface conditions ( V S30 ) and deeper geologic structure, as represented to the depth at which the profile reaches a velocity of 1.0 km/s ( Z 1.0 ). To demonstrate the advantages of our new SVM, we apply it to the Cascadia Region of North America, where numerous geologic basins and glaciated landscapes give rise to a wide range of V S30 and Z 1.0 combinations. This soil velocity model yields good estimates of site response across all site conditions, and significantly improves upon a model calibrated using only V S30 data. In conjunction with existing models that describe the deep velocity structure of the region (e.g., (Stephenson et al., 2017) [27]; the proposed model is particularly suited for use in regional-scale predictions of site response, liquefaction, landslides, infrastructure damage, and loss. The proposed methodology is broadly applicable to the development of SVMs elsewhere, and with improved understanding of near-surface and deep velocity structures, can facilitate more accurate ground-motion predictions globally.

Washington↗
USGS2020

Ensemble ShakeMaps for magnitude 9 earthquakes on the Cascadia Subduction Zone

We develop ensemble ShakeMaps for various magnitude 9 ( ⁠ M "> M M 9) earthquakes on the Cascadia megathrust. Ground‐shaking estimates are based on 30 M "> M M 9 Cascadia earthquake scenarios, which were selected using a logic‐tree approach that varied the hypocenter location, down‐dip rupture limit, slip distribution, and location of strong‐motion‐generating subevents. In a previous work, Frankel et al. (2018) used a hybrid approach (i.e., 3D deterministic simulations for frequencies &lt; 1 &#x2009;&#x2009; Hz "> < 1 Hz <1 Hz and stochastic synthetics for frequencies &gt; 1 &#x2009;&#x2009; Hz "> > 1 Hz >1 Hz ⁠ ) and uniform site amplification factors to create broadband seismograms from this set of 30 earthquake scenarios. Here, we expand on this work by computing site‐specific amplification factors for the Pacific Northwest and applying these factors to the ground‐motion estimates derived from Frankel et al. (2018) . In addition, we use empirical ground‐motion models (GMMs) to expand the ground‐shaking estimates beyond the original model extent of Frankel et al. (2018) to cover all of Washington State, Oregon, northern California, and southern British Columbia to facilitate the use of these ensemble ShakeMaps in region‐wide risk assessments and scenario planning exercises. Using this updated set of 30 M "> M M 9 Cascadia earthquake scenarios, we present ensemble ShakeMaps for the median, 2nd, 16th, 84th, and 98th percentile ground‐motion intensity measures. Whereas traditional scenario ShakeMaps are based on a single hypothetical earthquake rupture, our ensemble ShakeMaps take advantage of a logic‐tree approach to estimating ground motions from multiple earthquake rupture scenarios. In addition, 3D earthquake simulations capture important features such as strong ground‐motion amplification in the Pacific Northwest’s sedimentary basins, which are not well represented in the empirical GMMs that compose traditional scenario ShakeMaps. Overall, our results highlight the importance of strong‐motion‐generating subevents for coastal sites, as well as the amplification of long‐period ground shaking in deep sedimentary basins, compared with previous scenario ShakeMaps for Cascadia.

California, Oregon, Washington↗
USGS2019

Impacts of simulated M9 Cascadia Subduction Zone motions on idealized systems

Ground motions have been simulated for a magnitude 9 (M9) Cascadia Subduction Zone earthquake, which will affect the Puget Lowland region, including cities underlain by the Seattle, Everett, and Tacoma sedimentary basins. The current national seismic maps do not account for the effects of these basins on the risk-targeted Maximum Considered Earthquake (MCE R ). The simulated motions for Seattle had large spectral accelerations (at a period of 2 s, 43% of simulated M9 motions exceeded the MCE R ), damaging spectral shapes (particularly at periods near 1 s), and long durations (5%–95% significant durations near 110 s). For periods of 1 s or longer, the resulting deformation demands and collapse likelihood for four sets of single-degree-of-freedom systems exceeded the corresponding values for motions consistent with the conditional mean spectra at the MCE R intensity (MCE R ). The regional variation of damage was estimated by combining probabilistic characterizations of the seismic resistance of structures and of the effective spectral acceleration, S a,eff , which accounts for the effects of spectral acceleration, spectral shape, and ground-motion duration. For high-strength, low-ductility systems located above deep basins ( Z 2.5 > 6 km), the likelihood of collapse during an M9 earthquake averaged 13% and 18% at 1.0 s and 2.0 s periods, respectively. For low-strength, high-ductility systems, the corresponding likelihoods of collapse averaged 18% and 7%.

Washington↗
USGS2018

Broadband synthetic seismograms for magnitude 9 earthquakes on the Cascadia Megathrust based on 3D simulations and stochastic synthetics, Part 2: Rupture parameters and variability

We used a combination of 3D finite‐difference simulations ( &lt; 1 &#x2009;&#x2009; Hz "> < 1 Hz ) and 1D stochastic synthetics ( &gt; 1 &#x2009;&#x2009; Hz "> > 1 Hz ) to generate broadband (0–10 Hz) synthetic seismograms for numerous M w "> Mw 9 earthquake rupture scenarios on the Cascadia megathrust. Slip consists of multiple high‐stress‐drop subevents ( M w "> Mw 8) with short rise times on the deeper portions of the fault, superimposed on a background slip distribution with longer rise times. We identify key rupture parameters that control the intensity of ground shaking and resulting seismic hazard; these include the hypocenter location, the down‐dip limit of slip, the average rupture velocity, and the character (i.e., location, magnitude, and stress drop) of subevents. Extending the down‐dip limit of rupture to the top of the nonvolcanic tremor zone results in localized regions with a factor of 5–10 increase in spectral acceleration (SA) for periods &lt; 5 &#x2009;&#x2009; s "> <5 s , compared to a rupture that is completely offshore. This is primarily due to the closer proximity of high‐stress‐drop subevents to inland locations when the rupture is allowed to extend deeper. Similarly, we find that the hypocenter location can result in a variation in the intensity of ground motions of a factor of &gt; 10 "> >10 , due to the effects of rupture directivity (i.e., SA at periods &gt; 1 &#x2009;&#x2009; s "> >1 s ). We also observe a coupling between rupture directivity and basin amplification. The intensity of ground motions is also strongly affected by the magnitude, stress drop, and location of high‐stress‐drop subevents, which are poorly constrained. Overall, our results quantify the effect of kinematic rupture parameters on ground motions for an M w "> Mw 9 earthquake in Cascadia and emphasize the need for further constraints on these parameters to improve seismic hazard estimates in the Pacific Northwest.

Bulletin of the Seismological Society of America↗
USGS2018

Broadband synthetic seismograms for magnitude 9 earthquakes on the Cascadia megathrust based on 3D simulations and stochastic synthetics, part 1: Methodology and overall results

We have produced a large set of broadband (0–10 Hz) synthetic seismograms for Mw 9.0 earthquakes on the Cascadia megathrust by combining synthetic seismograms derived from 3D finite‐difference simulations ( ⁠≤1 Hz⁠ ) with finite‐source, stochastic synthetics ( ⁠≥1 Hz⁠ ). We used a compound rupture model consisting of high stress drop M w 8 subevents superimposed on large, shallower slip with long‐slip duration, informed by observations of the M w 9.0 Tohoku, Japan, and M w 8.8 Maule, Chile, earthquakes. Thirty 3D simulations were run, considering a variety of rupture parameters, to determine the range of expected ground motions. For sites not in sedimentary basins, the spectral accelerations of the synthetics are similar to the BC Hydro ground‐motion prediction equations (GMPEs) for periods of 0.1–6 s, but exceed them at periods greater than 6 s. Response spectra from the synthetics at sites in the Seattle and Tacoma sedimentary basins show large amplifications of factors of 2–5 at periods of 1–10 s. This basin amplification is substantially larger than that found for crustal earthquakes in the Next Generation Attenuation‐West2 (NGA‐West2) GMPEs. Basin amplification is caused by basin‐edge generated surface waves and by amplification and focusing of S waves and surface waves by the 3D basin structure. The synthetic seismograms show effective average durations of strong motions of about 70 s for coastal sites, increasing to about 120 s at 200 km distance. We find that the interevent and intraevent standard deviations of the spectral amplitudes of the synthetics are larger for sites closer to the rupture, because they are more sensitive to the location of subevents and rupture directivity.

California, Oregon, Washington↗
USGS2018

Effects of simulated magnitude 9 earthquake motions on structures in the Pacific Northwest

The Cascadia Subduction Zone (CSZ) produces long-duration, large-magnitude earthquakes that could severely affect structures in the Pacific Northwest (PNW). The impact of synthetic M9.0 CSZ earthquakes on buildings in the Pacific Northwest is studied using eight reinforced concrete wall archetypes that range from 4 to 40 stories. These archetypes were subjected to an ensemble of simulated ground-motions from 30 M9 earthquakes for a location in Seattle that overlies a ~8km deep sedimentary basin and an equivalent location outside the basin. Long-period (1-7s) ground motions are strongly amplified in the CSZ synthetics within the Seattle Basin, leading to a 6-fold increase in deformation demand (inter-story drift) compared to an equivalent out of basin site. The variability in demand was also found to be much larger inside the basin than outside the basin.

Conference Paper↗
USGS2017

Effects of deep basins on structural collapse during large subduction earthquakes

Deep sedimentary basins are known to increase the intensity of ground motions, but this effect is implicitly considered in seismic hazard maps used in U.S. building codes. The basin amplification of ground motions from subduction earthquakes is particularly important in the Pacific Northwest, where the hazard at long periods is dominated by such earthquakes. This paper evaluates the effects of basins on spectral accelerations, ground-motion duration, spectral shape, and structural collapse using subduction earthquake recordings from basins in Japan that have similar depths as the Puget Lowland basin. For three of the Japanese basins and the Puget Lowland basin, the spectral accelerations were amplified by a factor of 2 to 4 for periods above 2.0 s. The long-duration subduction earthquakes and the effects of basins on spectral shape combined, lower the spectral accelerations at collapse for a set of building archetypes relative to other ground motions. For the hypothetical case in which these motions represent the entire hazard, the archetypes would need to increase up to 3.3 times its strength to compensate for these effects.

Earthquake Spectra↗
USGS2017

3-D simulations of M9 earthquakes on the Cascadia Megathrust: Key parameters and uncertainty

Geologic and historical records indicate that the Cascadia subduction zone is capable of generating large, megathrust earthquakes up to magnitude 9. The last great Cascadia earthquake occurred in 1700, and thus there is no direct measure on the intensity of ground shaking or specific rupture parameters from seismic recordings. We use 3-D numerical simulations to generate broadband (0-10 Hz) synthetic seismograms for 50 M9 rupture scenarios on the Cascadia megathrust. Slip consists of multiple high-stress drop subevents (~M8) with short rise times on the deeper portion of the fault, superimposed on a background slip distribution with longer rise times. We find a >4x variation in the intensity of ground shaking depending upon several key parameters, including the down-dip limit of rupture, the slip distribution and location of strong-motion-generating subevents, and the hypocenter location. We find that extending the down-dip limit of rupture to the top of the non-volcanic tremor zone results in a ~2-3x increase in peak ground acceleration for the inland city of Seattle, Washington, compared to a completely offshore rupture. However, our simulations show that allowing the rupture to extend to the up-dip limit of tremor (i.e., the deepest rupture extent in the National Seismic Hazard Maps), even when tapering the slip to zero at the down-dip edge, results in multiple areas of coseismic coastal uplift. This is inconsistent with coastal geologic evidence (e.g., buried soils, submerged forests), which suggests predominantly coastal subsidence for the 1700 earthquake and previous events. Defining the down-dip limit of rupture as the 1 cm/yr locking contour (i.e., mostly offshore) results in primarily coseismic subsidence at coastal sites. We also find that the presence of deep subevents can produce along-strike variations in subsidence and ground shaking along the coast. Our results demonstrate the wide range of possible ground motions from an M9 megathrust earthquake in Cascadia, and the potential to further constrain key rupture parameters using geologic and geophysical observations, ultimately improving our estimation of seismic hazard associated with the Cascadia megathrust.

Conference Paper↗