This Zenodo repository contains supplementary data and information for the paper "Sub- and supershear ruptures during the 2023 Mw 7.8 and Mw 7.6 earthquake doublet in SE Türkiye" by Melgar et al. in Seismica (https://doi.org/10.26443/seismica.v2i3.387). If you use any of the data or models please cite it accordingly For any questions or issues, email Diego Melgar at dmelgarm@uoregon.edu. Files in this repository: An aftershock catalog from hypoDD relocations - aftershocks.xlsx - the file has appropriate column headers, which should be self-explanatory. GNSS and strong motion waveforms used in the slip inversion in folder - Waveforms - the waveforms are in miniSEED format and are in physical units. Acceleration data is in m/s/s, and GNSS data is in m. The coordinates for the stations are in the stations.txt file. Strong motion data is in north, east, up coordinates with channel codes HNN, HNE, and HNZ, respectively. HR-GNSS is in the same reference system but has channel codes LYN, LYE, and LYZ. Best fitting inverse models are in folder - Models - this folder includes two sub-folders, one for each rupture. Each subfolder contains: A .inv file, which has the full result of the multi-time window inversion A .inv.total file, which contains the total slip at each subfault A .log file that has inversion details A .mshout file that has the definition of the 3D mesh used for inversion An stf.txt file with the source time function Animation_S1: shown is the rupture propagation for both events. Plotted is the slip rate, which then fades into the final slip. Inset are the source time functions for both events. Figure S1: Station code names (A) Strong motion and (B) HR-GNSS. Fault traces are the same as in Figure 1. Figure S2: Waveform fits for the Mw7.8 event when considering rupture on the EAFZ triggering at the time of the arrival of S-waves from the hypocenter. Compare to Figure 4 for reference to the delayed triggering of the EAFZ. Notice how stations close to the hypocenter (HR-GNSS stations ANTE, ONIY, and strong motion stations 2712, 2718, locations in Figure S1) have particularly poor fits. This is most obvious in the east component of ONIY highlighted in the dashed blue line. Here the initial phase of ground motion in the synthetic (red) is much larger than in the observed (black), hinting at the need for delaying slip until rupture from the Nurdagi-Pazarcik fault reaches the EAFZ, as in Animation S1. Figure S3: Same as Figure 4 but considering a maximum rupture speed of 3.8 km/s instead of 3.2 km/s. This would be equivalent to a supershear rupture.