We have implemented and verified a parallel-series Iwan-type nonlinear model in a 3D fourth-order staggered-grid velocity–stress finite-difference method. The Masing unloading and reloading behavior is simulated by tracking an overlay of concentric von Mises yield surfaces. Lamé parameters and failure stresses pertaining to each surface are calibrated to reproduce the stress–strain backbone curve, which is controlled by the reference strain assigned to a given depth level. The implementation is successfully verified against established codes for 1D and 2D SH-wave benchmarks. The capabilities of the method for large-scale nonlinear earthquake modeling are demonstrated for an Mw 7.8 dynamic rupture ShakeOut scenario on the southern San Andreas fault. Although ShakeOut simulations with a single yield surface reduces long-period ground-motion amplitudes by about 25% inside a waveguide in greater Los Angeles, Iwan nonlinearity further reduces the values by a factor of 2. For example, inside the Whittier Narrows corridor spectral accelerations at a period of 3 s are reduced from 1g in the linear case to about 0.8 in the bilinear case and to 0.3–0.4g in the multisurface Iwan nonlinear case, depending on the choice of reference strain. Normalized shear modulus reductions reach values of up to 50% in the waveguide and up to 75% in the San Bernardino basin at the San Andreas fault. We expect the implementation to be a valuable tool for future nonlinear 3D dynamic rupture and ground-motion simulations in models with coupled source, path, and site effects.
more » « less- PAR ID:
- 10526790
- Publisher / Repository:
- Seismological Society of America
- Date Published:
- Journal Name:
- Bulletin of the Seismological Society of America
- Edition / Version:
- 6
- Volume:
- 113
- Issue:
- 6
- ISSN:
- 0037-1106
- Page Range / eLocation ID:
- 2275 to 2291
- Subject(s) / Keyword(s):
- Nonlinear ground motions simulation of ground motion finite difference modeling
- Format(s):
- Medium: X Size: 2MB Other: pdf
- Size(s):
- 2MB
- Sponsoring Org:
- National Science Foundation
More Like this
-
AWP-ODC is a 4th-order finite difference code used by the SCEC community for linear wave propagation, Iwan-type nonlinear dynamic rupture and wave propagation, and Strain Green Tensor simulation. We have ported and verified the CUDA-version of AWP-ODC-SGT, a reciprocal version used in the SCEC CyberShake project, to HIP so that it can also run on AMD GPUs. This code achieved sustained 32.6 Petaflop/s performance and 95.6% parallel efficiency at full scale on Frontier, a Leadership Computing Facility at Oak Ridge National Laboratory. The readiness of this community software on AMD Radeon Instinct GPUs and EPYC CPUs allows SCEC to take advantage of exascale systems to produce more realistic ground motions and accurate seismic hazard products. We have also deployed AWP-ODC to Azure to leverage the array of tools and services that Azure provides for tightly coupled HPC simulation on commercial cloud. We collaborated with Internet 2/Azure Accelerator supporting team, as part of Microsoft Internet2/Azure Accelerator for Research Fall 2022 Program, with Azure credits awarded through Cloudbank, an NSF-funded initiative. We demonstrate the AWP performance with a benchmark of ground motion simulation on various GPU based cloud instances, and a comparison of the cloud solution to on-premises bare-metal systems. AWP-ODC currently achieves excellent speedup and efficiency on CPU and GPU architectures. The Iwan-type dynamic rupture and wave propagation solver faces significant challenges, however, due to the increased computational workload with the number of yield surfaces chosen. Compared to linear solution, the Iwan model adds 10x-30x more computational time plus 5x-13x more memory consumption that require substantial code changes to obtain excellent performance. Supported by NSF’s Characteristic Science Applications (CSA) program for the Leadership-Class Computing Facility (LCCF) at Texas Advanced Computing Center (TACC), we are porting and improving the performance of this nonlinear AWP-ODC software, preparing for the next generation NSF LCCF system called Horizon, to be installed at TACC. During Texascale days on the current TACC’s Frontera, we carried out an Iwan-type nonlinear dynamic rupture and wave propagation simulation of a Mw7.8 scenario earthquake on the southern San Andreas fault. This simulation modeled 83 seconds of rupture with a grid spacing of 25 m to resolve frequencies up to 4 Hz with a minimum shear-wave velocity of 500 m/s.more » « less
-
We have conducted three-dimensional (3D) 0–7.5 Hz physics-based wave propagation simulations to model the seismic response of the Long Valley Dam (LVD), which has formed Lake Crowley in Central California, to estimate peak ground motions and settlement of the dam expected during maximum credible earthquake (MCE) scenarios on the nearby Hilton Creek Fault (HCF). We calibrated the velocity structure, anelastic attenuation model, and the overall elastic properties of the dam via linear simulations of a Mw3.7 event as well as the Mw6.2 Chalfant Valley earthquake of 1986, constrained by observed ground motions on and nearby the LVD. The Statewide California Earthquake Center (SCEC) Community Velocity Model CVM-S4.26.M01 superimposed with a geotechnical layer using [Formula: see text] information tapered from the surface to a 700-m depth was used in the simulations. We found optimal fit of simulated and observed ground motions at the LVD using frequency-independent attenuation of [Formula: see text] ([Formula: see text] in m/s). Using the calibrated model, we simulated 3D nonlinear ground motions at the LVD for Mw6.6 rupture scenarios on the HCF using an Iwan-type, multi-yield-surface technique. We use a two-step method where the computationally expensive nonlinear calculations were carried out in a small domain with the plane wave excitation along the bottom boundary obtained from a full-domain 3D linear finite-fault simulation. Our nonlinear MCE simulation results show that peak ground velocities (PGVs) and peak ground accelerations (PGAs) as high as 72 cm/s and 0.55 g, respectively, can be expected at the crest of the LVD. Compared with linear ground motion simulation results, our results show that Iwan nonlinear damping reduces PGAs on the dam crest by up to a factor of 8 and increasingly depletes the high-frequency content of the waves toward the dam crest. We find horizontal relative displacements of the material inside the dam of up to [Formula: see text] and up to [Formula: see text] of vertical subsidence, equivalent to 1% of the dam height.
-
Abstract We first explore a series of retrospective earthquake interactions in southern California. We find that the four Mw≥7 shocks in the past 150 yr brought the Ridgecrest fault ∼1 bar closer to failure. Examining the 34 hr time span between the Mw 6.4 and Mw 7.1 events, we calculate that the Mw 6.4 event brought the hypocentral region of the Mw 7.1 earthquake 0.7 bars closer to failure, with the Mw 7.1 event relieving most of the surrounding stress that was imparted by the first. We also find that the Mw 6.4 cross-fault aftershocks shut down when they fell under the stress shadow of the Mw 7.1. Together, the Ridgecrest mainshocks brought a 120 km long portion of the Garlock fault from 0.2 to 10 bars closer to failure. These results motivate our introduction of forecasts of future seismicity. Most attempts to forecast aftershocks use statistical decay models or Coulomb stress transfer. Statistical approaches require simplifying assumptions about the spatial distribution of aftershocks and their decay; Coulomb models make simplifying assumptions about the geometry of the surrounding faults, which we seek here to remove. We perform a rate–state implementation of the Coulomb stress change on focal mechanisms to capture fault complexity. After tuning the model through a learning period to improve its forecast ability, we make retrospective forecasts to assess model’s predictive ability. Our forecast for the next 12 months yields a 2.3% chance of an Mw≥7.5 Garlock fault rupture. If such a rupture occurred and reached within 45 km of the San Andreas, we calculate it would raise the probability of a San Andreas rupture on the Mojave section by a factor of 150. We therefore estimate the net chance of large San Andreas earthquake in the next 12 months to be 1.15%, or about three to five times its background probability.more » « less
-
The ShakeOut simulations of a M7.8 earthquake on the southern San Andreas fault (2008) studies predicted unexpectedly large ground motions throughout Southern California due to waveguide e#ects from interconnected sedimentary basins in 3D velocity model. Here, we re-examine the ground motion predictions from the ShakeOut scenario using the most recent updates on the velocity models and the realistic surface topography from the Digital Elevation Model. The exceptional scalability and performance of the AWP-ODC due to the most recent advancements allows for examination of plenty earth models including the irregular surface topography at low computational cost.more » « less
-
Abstract Understanding mechanical conditions that lead to complexity in earthquakes is important to seismic hazard analysis. In this study, we simulate physics‐based multicycle dynamic models of the San Andreas fault (Carrizo through San Bernardino sections) and the San Jacinto fault (Claremont and Clark strands). We focus on a complex fault geometry based on the Southern California Earthquake Center Community Fault Model and its effect over multiple earthquake cycles. Using geodetically derived strain rates, we validate the models against geologic slip rates and recurrence intervals at various paleoseismic sites. We find that the interactions among fault geometry, dynamic rupture and interseismic stress accumulation produce stress heterogeneities, leading to rupture segmentation and variability in earthquake recurrence. Our models produce earthquakes with rupture extents similar to a recent comprehensive paleoseismic catalog. The “earthquake gates” of the Big Bend and the Cajon Pass occasionally impede dynamic ruptures. The angle of compression, which is the subtraction of the maximum shear strain rate direction from the local fault strike, can better determine the likelihood of the impedance of restraining bends to dynamic ruptures. Because the Big Bend has an angle of compression of ∼20°, ruptures that traverse the Big Bend, like the 1857 Fort Tejon earthquake, are more frequent than expected based on empirical relations which predict the ∼40° restraining bend to terminate most ruptures. Our models indicate that large ruptures tend to initiate north of the Big Bend and propagate southwards, similar to the 1857 earthquake, providing critical information for ground shaking assessment in the region.