Geology topics

Shuo Ma

Publications and source records attributed to Shuo Ma.

4 recordsLinked to original sources

USGS2026

CRESCENT earthquake dynamic rupture, earthquake cycle, and tsunami code verification platform

Physics-based simulations are critical for understanding natural hazards. The increasing complexity of numerical codes requires benchmark exercises to verify that different computational methods yield consistent results when solving the same governing equations. Here, we present an open-access web platform designed for the verification of earthquake dynamic rupture, seismic cycle, and tsunami simulations. The platform architecture utilizes a modular, serverless backend on Amazon Web Services (AWS) to provide scalable file processing and visualization. A lightweight static web application provides a secure interface for uploading and managing results, while the browser-based data visualization enables interactive analysis of time series and surface grid data. By using structured JavaScript Object Notation (JSON) text files to define benchmark structures, the system remains fully extensible, allowing the addition of new scenarios without modifying the underlying software logic. The platform hosts the "The Tsunami Problem Versions" (TTPV) 1 & 2, two benchmarks for 3D fully coupled earthquake dynamic rupture and tsunami generation, and provides a framework for earthquake cycle models. This community resource aims to build trust in numerical simulations and facilitate long-term collaborative code verification as modeling software continues to evolve.

Seismica↗

USGS2018

A suite of exercises for verifying dynamic earthquake rupture codes

We describe a set of benchmark exercises that are designed to test if computer codes that simulate dynamic earthquake rupture are working as intended. These types of computer codes are often used to understand how earthquakes operate, and they produce simulation results that include earthquake size, amounts of fault slip, and the patterns of ground shaking and crustal deformation. The benchmark exercises examine a range of features that scientists incorporate in their dynamic earthquake rupture simulations. These include implementations of simple or complex fault geometry, off‐fault rock response to an earthquake, stress conditions, and a variety of formulations for fault friction. Many of the benchmarks were designed to investigate scientific problems at the forefronts of earthquake physics and strong ground motions research. The exercises are freely available on our website for use by the scientific community.

Seismological Research Letters↗

USGS2010

Ground-motion modeling of Hayward fault scenario earthquakes, part II: Simulation of long-period and broadband ground motions

We simulate long-period ( T> 1.0–2.0 s) and broadband ( T> 0.1 s) ground motions for 39 scenario earthquakes ( M w 6.7–7.2) involving the Hayward, Calaveras, and Rodgers Creek faults. For rupture on the Hayward fault, we consider the effects of creep on coseismic slip using two different approaches, both of which reduce the ground motions, compared with neglecting the influence of creep. Nevertheless, the scenario earthquakes generate strong shaking throughout the San Francisco Bay area, with about 50% of the urban area experiencing modified Mercalli intensity VII or greater for the magnitude 7.0 scenario events. Long-period simulations of the 2007 M w 4.18 Oakland earthquake and the 2007 M w 5.45 Alum Rock earthquake show that the U.S. Geological Survey’s Bay Area Velocity Model version 08.3.0 permits simulation of the amplitude and duration of shaking throughout the San Francisco Bay area for Hayward fault earthquakes, with the greatest accuracy in the Santa Clara Valley (San Jose area). The ground motions for the suite of scenarios exhibit a strong sensitivity to the rupture length (or magnitude), hypocenter (or rupture directivity), and slip distribution. The ground motions display a much weaker sensitivity to the rise time and rupture speed. Peak velocities, peak accelerations, and spectral accelerations from the synthetic broadband ground motions are, on average, slightly higher than the Next Generation Attenuation (NGA) ground-motion prediction equations. We attribute much of this difference to the seismic velocity structure in the San Francisco Bay area and how the NGA models account for basin amplification; the NGA relations may underpredict amplification in shallow sedimentary basins. The simulations also suggest that the Spudich and Chiou (2008) directivity corrections to the NGA relations could be improved by increasing the areal extent of rupture directivity with period.

California↗

USGS2010

Inelastic off-fault response and three-dimensional dynamics of earthquake rupture on a strike-slip fault

Large dynamic stress off the fault incurs an inelastic response and energy loss, which contributes to the fracture energy, limiting the rupture and slip velocity. Using an explicit finite element method, we model three-dimensional dynamic ruptures on a vertical strike-slip fault in a homogeneous half-space. The material is subjected to a pressure-dependent Drucker-Prager yield criterion. Initial stresses in the medium increase linearly with depth. Our simulations show that the inelastic response is confined narrowly to the fault at depth. There the inelastic strain is induced by large dynamic stresses associated with the rupture front that overcome the effect of the high confining pressure. The inelastic zone increases in size as it nears the surface. For material with low cohesion (~5 MPa) the inelastic zone broadens dramatically near the surface, forming a "flowerlike" structure. The near-surface inelastic strain occurs in both the extensional and the compressional regimes of the fault, induced by seismic waves ahead of the rupture front under a low confining pressure. When cohesion is large (~10 MPa), the inelastic strain is significantly reduced near the surface and confined mostly to depth. Cohesion, however, affects the inelastic zone at depth less significantly. The induced shear microcracks show diverse orientations near the surface, owing to the low confining pressure, but exhibit mostly horizontal slip at depth. The inferred rupture-induced anisotropy at depth has the fast wave direction along the direction of the maximum compressive stress.

Journal of Geophysical Research B: Solid Earth↗