Search for: All records

Abstract With the rise of data volume and computing power, seismological research requires more advanced skills in data processing, numerical methods, and parallel computing. We present the experience of conducting training workshops in various forms of delivery to support the adoption of large-scale high-performance computing (HPC) and cloud computing, advancing seismological research. The seismological foci were on earthquake source parameter estimation in catalogs, forward and adjoint wavefield simulations in 2D and 3D at local, regional, and global scales, earthquake dynamics, ambient noise seismology, and machine learning. This contribution describes the series of workshops delivered as part of research projects, the learning outcomes for participants, and lessons learned by the instructors. Our curriculum was grounded on open and reproducible science, large-scale scientific computing and data mining, and computing infrastructure (access and usage) for HPC and the cloud. We also describe the types of teaching materials that have proven beneficial to the instruction and the sustainability of the program. We propose guidelines to deliver future workshops on these topics. 

Abstract Seismic faults are surrounded by damaged rocks with reduced rigidity and enhanced attenuation. These damaged fault zone structures can amplify seismic waves and affect earthquake dynamics, yet they are typically omitted in physics‐based regional ground motion simulations. We report on the significant effects of a shallow, flower‐shaped fault zone in foreshock‐mainshock 3D dynamic rupture models of the 2019 Ridgecrest earthquake sequence. We find that the fault zone structure both amplifies and reduces ground motions not only locally but at distances exceeding 100 km. This impact on ground motions is frequency‐ and magnitude‐dependent, particularly affecting higher frequency ground motions from the foreshock because its corner frequency is closer to the fault zone's fundamental eigenfrequency. Within the fault zone, the shallow transition to a velocity‐strengthening frictional regime leads to a depth‐dependent peak slip rate increase of up to 70% and confines fault zone‐induced supershear transitions mostly to the fault zone's velocity‐weakening roots. However, the interplay of fault zone waves, free surface reflections, and rupture directivity can generate localized supershear rupture, even in narrow velocity‐strengthening regions, which are typically thought to inhibit supershear rupture. This study demonstrates that shallow fault zone structures may significantly affect intermediate‐ and far‐field ground motions and cause localized supershear rupture penetrating into velocity‐strengthening regions, with important implications for seismic hazard assessment. 

Abstract We develop an earthquake simulator to study the partitioning of seismic/aseismic slip and dynamics of Earthquakes on a Heterogeneous strike‐slip Fault (HFQsim) using a generalized model of a discrete fault governed by static/dynamic friction and creep in an elastic half‐space. Previous versions of the simulator were shown to produce various realistic seismicity patterns (e.g., frequency‐magnitude event statistics, hypocenter and slip distributions, temporal occurrence) using friction levels and creep properties that vary in space but are fixed in time. The new simulator incorporates frictional heat generation by earthquake slip leading to temperature rises, subsequent diffusion cooling into the half space, and time‐dependent creep on the fault. The model assumes a power law dependence of creep velocity on the local shear stress, with temperature‐dependent coefficients based on the Arrhenius equation. Temperature rises due to seismic slip produce increased aseismic slip, which can lead to further stress concentrations, aftershocks, and heat generation in a feedback loop. The partitioning of seismic/aseismic slip and space‐time evolution of seismicity are strongly affected by the temperature changes on the fault. The results are also affected significantly by the difference between the static and kinetic friction levels. The model produces realistic spatio‐temporal distribution of seismicity, transient aseismic slip patterns, mainshock‐aftershock sequences, and a bimodal distribution of earthquakes with background and clustered events similar to observations. The HFQsim may be used to clarify relations between fault properties and different features of seismicity and aseismic slip, and to improve the understanding of failure patterns preceding large earthquakes. 

Abstract The nonlinear mechanical responses of rocks and soils to seismic waves play an important role in earthquake physics, influencing ground motion from source to site. Continuous geophysical monitoring, such as ambient noise interferometry, has revealed co‐seismic wave speed reductions extending tens of kilometers from earthquake sources. However, the mechanisms governing these changes remain challenging to model, especially at regional scales. Using a nonlinear damage model constrained by laboratory experiments, we develop and apply an open‐source 3D discontinuous Galerkin method to simulate regional co‐seismic wave speed changes during the 2015 M_w7.8 Gorkha earthquake. We find pronounced spatial variations of co‐seismic wave speed reduction, ranging from <0.01% to >50%, particularly close to the source and within the Kathmandu Basin, while disagreement with observations remains. The most significant reduction occurs within the sedimentary basin and varies with basin depths, whereas wave speed reductions correlate with the fault slip distribution near the source. By comparing ground motions from simulations with elastic, viscoelastic, elastoplastic, and nonlinear damage rheologies, we demonstrate that the nonlinear damage model effectively captures low‐frequency ground motion amplification due to strain‐dependent wave speed reductions in soft sediments. We verify the accuracy of our approach through comparisons with analytical solutions and assess its scalability on high‐performance computing systems. The model shows near‐linear strong and weak scaling up to 2,048 nodes, enabling efficient large‐scale simulations. Our findings provide a physics‐based framework to quantify nonlinear earthquake effects and emphasize the importance of damage‐induced wave speed variations for seismic hazard assessment and ground motion predictions. 

ABSTRACT Fault zones exhibit geometrical complexity and are often surrounded by multiscale fracture networks within their damage zones, potentially influencing rupture dynamics and near-field ground motions. In this study, we investigate the ground-motion characteristics of cascading ruptures across damage zone fracture networks of moderate-size earthquakes (Mw 5.5–6.0) using high-resolution 3D dynamic rupture simulations. Our models feature a listric normal fault surrounded by more than 800 fractures, emulating a major fault and its associated damage zone. We analyze three cases: a cascading rupture propagating within the fracture network (Mw 5.5), a non-cascading main-fault rupture with off-fault fracture slip (Mw 6.0), and a main-fault rupture without a fracture network (Mw 6.0). Cascading ruptures within the fracture network produce distinct ground-motion signatures with enriched high-frequency content, arising from simultaneous slip of multiple fractures and parts of the main fault, resembling source coda-wave-like signatures. This case shows elevated near-field characteristic frequency (fc) and stress drop, approximately an order of magnitude higher than the estimation directly on the fault of the dynamic rupture simulation. The inferred fc of the modeled vertical ground-motion components reflects the complexity of the radiation pattern and rupture directivity of fracture-network cascading earthquakes. We show that this is consistent with observations of strong azimuthal dependence of corner frequency in the 2009–2016 central Apennines, Italy, earthquake, sequence. Simulated ground motions from fracture-network cascading ruptures also show pronounced azimuthal variations in peak ground acceleration (PGA), peak ground velocity, and pseudospectral acceleration, with average PGA nearly double that of the non-cascading cases. Cascading ruptures radiate high-frequency seismic energy, yield nontypical ground-motion characteristics including coda-wave-like signatures, and may result in a significantly higher seismologically inferred stress drop and PGA. Such outcomes emphasize the critical role of fault-zone complexity in affecting rupture dynamics and seismic radiation and have important implications for physics-based seismic hazard assessment. 

Abstract The 1 January 2024, moment magnitude 7.5 Noto Peninsula earthquake ruptured in complex ways, challenging analysis of its tsunami generation. We present tsunami models informed by a 6‐subevent centroid moment tensor (CMT) model obtained through Bayesian inversion of teleseismic and strong motion data. We identify two distinct bilateral rupture episodes. Initial, onshore rupture toward the southwest is followed by delayed re‐nucleation at the hypocenter, likely aided by fault weakening, causing significant seafloor uplift to the northeast. We construct a complex multi‐fault uplift model, validated against geodetic observations, that aligns with known fault system geometries and is critical in modeling the observed tsunami. The simulations can explain tsunami wave amplitude, timing, and polarity of the leading wave, which are crucial for tsunami early warning. Upon comparison with alternative source models and analysis of 2000 multi‐CMT ensemble solutions, we highlight the importance of incorporating complex source effects for realistic tsunami simulations. 

Abstract Despite routine detection of coseismic acoustic‐gravity waves (AGWs) in Global Navigation Satellite System (GNSS) total electron content (TEC) observations, models of the earthquake‐atmosphere‐ionosphere dynamics, essential for validating data‐driven studies, remain limited. We present the results of three‐dimensional numerical simulations encompassing the entire coupling from Earth's interior to the ionosphere during the 7.8 2016 Kaikoura earthquake. Incorporating the impact of data/model uncertainties in estimating the ionospheric state, the results show a good agreement between observed and simulated slant TEC (sTEC) signals, assessed through a set of metrics. The signals exhibit intricate waveforms, resulting from the integrated nature of TEC and phase cancellation effects, emphasizing the significance of direct signal comparisons along realistic line‐of‐sight paths. By comparing simulation results initialized with kinematic and dynamic source models, the study demonstrates the quantifiable sensitivity of sTEC to AGW source specifications, pointing to their utility in the analysis of coupled dynamics. 

Abstract Previous geodetic and teleseismic observations of the 2021 7.4 Maduo earthquake imply surprising but difficult‐to‐constrain complexity, including rupture across multiple fault segments and supershear rupture. Here, we present an integrated analysis of multi‐fault 3D dynamic rupture models, high‐resolution optical correlation analysis, and joint optical‐InSAR slip inversion. Our preferred model, validated by the teleseismic multi‐peak moment rate release, includes unilateral eastward double‐onset supershear speeds and cascading rupture dynamically triggering two adjacent fault branches. We propose that pronounced along‐strike variation in fracture energy, complex fault geometries, and multi‐scale variable prestress drives this event's complex rupture dynamics. We illustrate how supershear transition has signatures in modeled and observed off‐fault deformation. Our study opens new avenues to combine observations and models to better understand complex earthquake dynamics, including local and potentially repeating supershear episodes across immature faults or under heterogeneous stress and strength conditions, which are potentially not unusual. 

Abstract Understanding the dynamics of microearthquakes is a timely challenge with the potential to address current paradoxes in earthquake mechanics, and to better understand earthquake ruptures induced by fluid injection. We perform fully 3D dynamic rupture simulations caused by fluid injection on a target fault for Fault Activation and Earthquake Ruptures experiments generatingM_w ≤ 1 earthquakes. We investigate the dynamics of rupture propagation with spatially variable stress drop caused by pore pressure changes and assuming different slip‐weakening constitutive parameters. We show that the spontaneous arrest of propagating ruptures is possible by assuming a high fault strength parameter S, that is, a high ratio between strength excess and dynamic stress drop. In faults with high S values (low rupturing potential), even minor variations inD_c(from 0.45 to 0.6 mm) have a substantial effect on the rupture propagation and the ultimate earthquake size. Modest spatial variations of dynamic stress drop determine the rupture mode, distinguishing self‐arresting from run‐away ruptures. Our results suggest that several characteristics inferred for accelerating dynamic ruptures differ from those observed during rupture deceleration of a self‐arresting earthquake. During deceleration, a decrease of peak slip velocity is associated with a nearly constant cohesive zone size. Moreover, the residual slip velocity value (asymptotic value for a crack‐like rupture) decreases to nearly zero. This means that an initially crack‐like rupture becomes a pulse‐like rupture during spontaneous arrest. These findings highlight the complex dynamics of small induced earthquakes, which differ from solutions obtained from conventional crack‐like models of earthquake rupture. 

Abstract The Mediterranean Hellenic Arc subduction zone (HASZ) has generated several 8 earthquakes and tsunamis. Seismic‐probabilistic tsunami hazard assessment typically utilizes uniform or stochastic earthquake models, which may not represent dynamic rupture and tsunami generation complexity. We present an ensemble of ten 3D dynamic rupture earthquake scenarios for the HASZ, utilizing a realistic slab geometry. Our simplest models use uniform along‐arc pre‐stresses or a single circular initial stress asperity. We then introduce progressively more complex models varying initial shear stress along‐arc, multiple asperities based on scale‐dependent critical slip weakening distance, and a most complex model blending all aforementioned heterogeneities. Thereby, regional initial conditions are constrained without relying on detailed geodetic locking models. Varying epicentral locations in the simplest, homogeneous model leads to different rupture speeds and moment magnitudes. We observe dynamic fault slip penetrating the shallow slip‐strengthening region and affecting seafloor uplift. Off‐fault plastic deformation can double vertical seafloor uplift. A single‐asperity model generates a 8 scenario resembling the 1303 Crete earthquake. Using along‐strike varying initial stresses results in 8.0–8.5 dynamic rupture scenarios with diverse slip rates and uplift patterns. The model with the most heterogeneous initial conditions yields a 7.5 scenario. Dynamic rupture complexity in prestress and fracture energy tends to lower earthquake magnitude but enhances tsunamigenic displacements. Our results offer insights into the dynamics of potential large Hellenic Arc megathrust earthquakes and may inform future physics‐based joint seismic and tsunami hazard assessments.