Parameter Files and Simulation Inputs for "When Earthquakes Cross the Gap' Dynamic Rupture Models"

Description This public repository contains all parameter files, input configurations, and meshes used to explore the dynamic rupture scenarios investigated in our study on rupture jumps across step-overs in normal-faulting earthquakes. These datasets reproduce precisely the conditions used in the simulations presented in our manuscript. This repository contains two main folders: parameter_files: Folder containing all input parameter files and meshes. figures: Folder containing the Python scripts used to reproduce the manuscript figures. Below is a concise description of the contents and naming conventions. parameter_files folder Main Control File parameters.par Control file defining: paths to all YAML parameter files initializing spatially dependent fields (wave propagation properties, initial stress, physical properties). It includes paths to mesh, coordinates of synthetic receivers. additional input variables governing simulation control and output generation (paths, formats). this control file can be modified in order to prescribed a desired parametric configuration. Read the following sections and naming code to better understand how to modified this parameter file to recreate a given result. Parametric Input Files These files contain the variables that change between scenarios. fault_material_off*_gap*.yamlMaterial properties on the fault planes, varying according to the relative geometry and location between faults. fault_stress_off*_gap*.yamlInitial stress conditions on the faults, varying with each configuration. Common Fixed Files Files used identically in all simulations: material.yamlBackground medium properties (shared across all scenarios). 4. Mesh Files laga_off*_gap*.xdmf laga_off*_gap* (associated mesh files to the XDMF) These provide the spatial discretization for each geometric configuration of the step-over system. Naming Convention Files whose names vary with the pattern off**_gap** correspond to different parametric configuration depending on the geometric settings of the two-fault system. The naming convention is: off(m/p)(0–5)_gap(m/p)(0–5).EXT Where: off = along-strike offset of the first (source) fault relative to the second fault gap = gap (perpendicular to the strike direction) distance between the two faults m / p = negative or positive sign (minus or plus) 0–5 = length of the distance in kilometers .EXT = file extension (yaml, xdmf, etc.) Interpretation example offp5p0_gapm2p0 → the emitting (first) fault is shifted +5.0 km along strike and –2.0 km across strike relative to the second fault (which is alwasy fixed at its position). These codes appear consistently in all parametric, stress, material, and mesh files. Here some other examples Offset | Gap | Code | -2.5 -1.0 offm2p5_gapm1p0 -5.0 +5.0 offm5p0_gapp5p0 0.0 -2.0 offp0p0_gapm2p0 +2.5 +3.0 offp2p5_gapp3p0 +5.0 -5.0 offp5p0_gapm5p0The remaining files (*.dat) describe the fault and off-fault receivers where the particle velocity and stress fields are saved during the simulations. figures folder Python scripts to reproduce the figures The figures/ folder contains Python scripts used to generate the figures of the manuscript from the simulation outputs. Main figures Fig2_MC_initial_state.py: Computes the initial stress field assuming Andersonian stress conditions and plots the Mohr–Coulomb stress state together with topographic lithostatic load over the two faults. Fig3_Static_DCFF.py: Computes and plots the static Coulomb stress change (ΔCFF) for two different cases (HW and FW sides) using SeisSol simulation output files. Fig4_summary_results.py: Produces a summary plot of rupture scenarios showing the dependence of earthquake magnitude and rupture jumps on fault offset, gap distance, initial stress levels and SHmax orientation. Profiles and cross-sections Fig5 Depth_profiles.py: Extracts and plots depth cross-sections of the stress variation along selected vertical profiles from the simulations at different time instants. Profile1.py: Generates the mapviews and vertical cross-sections showing stress and displacement evolution across the fault and the propagation medium for the HW and FW cases described in the manuscript. This figures are associated to the first vertical profile to the south. Profile2.py: Generates the same mapviews and the second vertical cross-sections for comparison with Profile 1. Supplementary figures FigS3-4_2D_sections.py: Produces 2D horizontal representations of the resulting outcomes of the dynamic rupture simulations: rupture jumps, arrested ruptures, magnitudes. FigS5_cumulative.py: Computes and plots the cumulative number of rupture jumps for all tested scenarios and SHmax directions. Model geometry Fig1_model.* Schematic representation of the fault geometry and model configuration used in the simulations. videos folder (Videos of rupture scenarios) The videos/ folder contains animations showing the time evolution of dynamic stress transfer between the emitter and receiver faults for two representative rupture scenarios. Video_Profile1.mp4: Snapshot sequence of the rupture evolution, comparing the footwall (FW) and hanging wall (HW) configurations. The animation includes map views, vertical cross-sections perpendicular to strike, and fault-parallel views of the receiver fault. Colors show shear stress along dip, grayscale shading indicates vertical particle velocity, and Coulomb stress change is displayed on the receiver fault. A short preliminary rupture on the receiver fault is visible in the HW case. Video_Profile2.mp4: Snapshot sequences, showing the evolution of the same rupture scenarios but the vertical profiles are located furhter north. The animation highlights the development of stress transfer patterns and the nucleation of the second rupture in the footwall (FW) case, visible in map view, cross-sections, and fault-parallel views. Both scenarios correspond to simulations with a ±5 km gap perpendicular to strike and a 2.5 km along-strike offset between faults (25% overlap of fault surface area). The videos illustrate the spatiotemporal evolution of dynamic stress interaction discussed in the Results section of the manuscript.