One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier Dmul shows Dmul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of Dmul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger Dmul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

Broadband physics-based simulated earthquake ground motions are utilized to characterize the regional-scale seismic risk to modern reinforced concrete (RC) structures. A highly dense dataset of ground motions covering a 100-km (Formula presented.) 40-km domain was generated using kinematic fault rupture models with varying rupture characteristics to represent shallow crustal earthquakes and resolved up to frequencies of 5 Hz. Over 40,000 nonlinear response history simulations of short- and mid-rise RC special moment frame buildings were conducted using simulation models that are capable of representing nonlinear behavior and component deterioration effects. The spatial variability of structural risk within a single earthquake scenario and between different rupture scenarios is examined, and the regions of strongest directivity effects and highest structural demands are identified. The structural demands may vary by factor of up to 8.0 at very short distances from the fault, and the large dispersion in the demands decreases significantly beyond a distance of 15 km. The interstory drift and member rotation demands are substantially impacted by important features of the geological structure and the characteristics of the rupture scenarios, particularly the presence of localized high-slip regions. The frequency characteristics of the structures are found to play an important role in determining the effects of near-fault ground motions on the structural response and expected damage. The results of this study suggest that the simulated ground motions, particularly those generated using the hybrid rupture approach, may offer reasonable structural risk estimates for low-frequency structures and conservative estimates for high-frequency structures.

Ground-motion simulations generated from physics-based wave propagation models are gaining increasing interest in the engineering community for their potential to inform the performance-based design and assessment of infrastructure residing in active seismic areas. A key prerequisite before the ground-motion simulations can be used with confidence for application in engineering domains is their comprehensive and rigorous investigation and validation. This article provides a four-step methodology and acceptance criteria to assess the reliability of simulated ground motions of not historical events, which includes (1) the selection of a population of real records consistent with the simulated scenarios, (2) the comparison of the distribution of Intensity Measures (IMs) from the simulated records, real records, and Ground-Motion Prediction Equations (GMPEs), (3) the comparison of the distribution of simple proxies for building response, and (4) the comparison of the distribution of Engineering Demand Parameters (EDPs) for a realistic model of a structure. Specific focus is laid on near-field ground motions (<10km) from large earthquakes (Mw7), for which the database of real records for potential use in engineering applications is severely limited. The methodology is demonstrated through comparison of (2490) near-field synthetic records with 5 Hz resolution generated from the Pitarka et al (2019) kinematic rupture model with a population of (38) pulse-like near-field real records from multiple events and, when applicable, with NGA-W2 GMPEs. The proposed procedure provides an effective method for informing and advancing the science needed to generate realistic ground-motion simulations, and for building confidence in their use in engineering domains.

We investigated the effects of fault geometry and assumed minimum shear wavespeed (VS min) on 3D ground-motion simulations (0–2.5 Hz) in general, using a moment magnitude (Mw) 6.5 earthquake on the Hayward fault (HF). Simulations of large earthquakes on the northeast-dipping HF using the U.S. Geological Survey (USGS) 3D seismic model have shown intensity asymmetry with stronger shaking for the Great Valley Sequence east of the HF (hanging wall) relative to the Franciscan Complex to the west (footwall). We performed simulations with three fault geometries in both plane-layered (1D) and 3D models. Results show that the nonvertical fault geometries result in larger motions on the hanging wall relative to the vertical fault for the same Earth model with up to 50% amplifications in single-component peak ground velocity (PGV) within 10 km of the rupture. Near-fault motions on the footwall are reduced for the nonvertical faults, but less than they are increased on the hanging wall. Simulations assuming VS min values of 500 and 250 m=s reveal that PGVs are on average 25% higher west of the HF when using the lower VS min, with some locations amplified by a factor of 3. Increasing frequency content from 2.5 to 5 Hz increases PGV values. Spectral ratios of these two VS min cases show average amplifications of 2–4 (0.5–1.5 Hz) for the lower VS min west of the fault. Large differences (up to 2×) in PGV across the HF from previous studies persist even for the case with a vertical fault or VS min of 250 m=s. We conclude that assuming a VS min of 500 m=s underestimates intensities west of the HF for frequencies above 0.5 Hz, and that low upper crustal (depth < 10 km) shear wavespeeds defined in the 3D model contribute most to higher intensities east of the HF.

Accurate measurements of the time-dependent deformations of a building during earthquake excitation are essential for interpretation of the dynamic response of the as-built system and for quantifying the seismic demands. Traditional approaches for monitoring building systems are based on strong motion accelerometers mounted at selected elevations. However, accelerometer-based systems do not directly measure the deformations of the structure, and can have significant limitations that make it challenging to correctly measure deformations, particularly permanent deformations from inelastic response. In the study described herein, computational simulations and experiments were combined to evaluate the potential of a new optically based sensor to directly measure time-dependent deformations of a building, including inelastic deformations. The sensor methodology includes corrections for localized structural member rotations and can provide estimates of the absolute accelerations at each floor. A laser-based system utilizing a recently developed discrete diode position sensor (DDPS) is evaluated, and the ability of such a system to measure earthquake induced transient deformations characterized by building interstory drift is demonstrated.

Interstory drift (ID) is a key response parameter for buildings subjected to lateral loads and is used to define performance-based limit states, allowable deformations and damage states in a number of seismic design codes and standards. An ability to rapidly and accurately measure both transient and residual ID during an earthquake would provide important observables for understanding the seismic demands and post-earthquake condition of a building. Accurate retrieval of ID from accelerometer-based instrumentation systems can be very challenging, if not impossible, as a result of instrumentation limitations and the post-processing associated with strong motion accelerometer data. This is particularly true for the case in which residual drifts occur during inelastic building response. In the study presented herein, a newly developed optical sensor system, designed specifically for directly measuring both transient and residual ID, was experimentally evaluated through shake table testing and computational simulations. The ability of the sensor to accurately measure ID is demonstrated and key operational considerations for sensor system deployment are examined through a model-based investigation.

In soil-structure interaction modeling of systems subjected to earthquake motions, it is classically assumed that the incoming wave field, produced by an earthquake, is unidimensional and vertically propagating. This work explores the validity of this assumption by performing earthquake soil-structure interaction modeling, including explicit modeling of sources, seismic wave propagation, site, and structure. The domain reduction method is used to couple seismic (near-field) simulations with local soil-structure interaction response. The response of a generic nuclear power plant model computed using full earthquake soil-structure interaction simulations is compared with the current state-of-the-art method of deconvolving in depth the (simulated) free-field motions, recorded at the site of interest, and assuming that the earthquake wave field is spatially unidimensional. Results show that the 1-D wave-field assumption does not hold in general. It is shown that the way in which full 3-D analysis results differ from those which assume a 1-D wave field is dependent on fault-to-site geometry and motion frequency content. It is argued that this is especially important for certain classes of soil-structure systems of which nuclear power plants subjected to near-field earthquakes are an example.

We performed fully deterministic broadband (0–4 Hz) high-performance computing ground motion simulations of a magnitude 7.0 scenario earthquake on the Hayward Fault (HF) in the San Francisco Bay Area of Northern California. Simulations consider average one-dimensional (1-D) and three-dimensional (3-D) anelastic structure with flat and topographic free surfaces. Ground motion intensity measures (GMIMs) for the 3-D model display dramatic differences across the HF due to geologic heterogeneity, with low wave speeds east of the HF amplifying motions. The median GMIMs agree well with Ground Motion Prediction Equations (GMPEs); however, the 3-D model generates more scatter than the 1-D model. Ratios of 3-D/1-D GMIMs from the same source allow isolation of path and site effects for the 3-D model. These ratios show remarkably similar trends as site-specific factors for the GMPE predictions, suggesting that wave propagation effects in our 3-D simulations are on average consistent with empirical data.

Accurate understanding and quantification of the risk to critical infrastructure posed by future large earthquakes continues to be a very challenging problem. Earthquake phenomena are quite complex and traditional approaches to predicting ground motions for future earthquake events have historically been empirically based whereby measured ground motion data from historical earthquakes are homogenized into a common data set and the ground motions for future postulated earthquakes are probabilistically derived based on the historical observations. This procedure has recognized significant limitations, principally due to the fact that earthquake ground motions tend to be dictated by the particular earthquake fault rupture and geologic conditions at a given site and are thus very site-specific. Historical earthquakes recorded at different locations are often only marginally representative. There has been strong and increasing interest in utilizing large-scale, physics-based regional simulations to advance the ability to accurately predict ground motions and associated infrastructure response. However, the computational requirements for simulations at frequencies of engineering interest have proven a major barrier to employing regional scale simulations. In a U.S. Department of Energy Exascale Computing Initiative project, the EQSIM application development is underway to create a framework for fault-to-structure simulations. This framework is being prepared to exploit emerging exascale platforms in order to overcome computational limitations. This article presents the essential methodology and computational workflow employed in EQSIM to couple regional-scale geophysics models with local soil-structure models to achieve a fully integrated, complete fault-to-structure simulation framework. The computational workflow, accuracy and performance of the coupling methodology are illustrated through example fault-to-structure simulations.