SUMMARY We present a kinematic slip model and a simulation of the ensuing tsunami for the 2020 Mw 7.0 Néon Karlovásion (Samos, Eastern Aegean Sea) earthquake, generated from a joint inversion of high-rate GNSS, strong ground motion and InSAR data. From the inversion, we find that the source time function has a total duration of ∼20 s with three peaks at ∼4, 7.5 and 15 s corresponding to the development of three asperities. Most of the slip occurs at the west of the hypocentre and ends at the northwest downdip edge. The peak slip is ∼3.3 m, and the inverted rake angles indicate predominantly normal faulting motion. Compared with previous studies, these slip patterns have essentially similar asperity location, rupture dimension and anticorrelation with aftershocks. Consistent with our study, most published papers show the source duration of ∼20 s with three episodes of increased moment releases. For the ensuing tsunami, the eight available gauge records indicate that the tsunami waves last ∼18–30 hr depending on location, and the response period of tsunami is ∼10–35 min. The initial waves in the observed records and synthetic simulations show good agreement, which indirectly validates the performance of the inverted slip model. However, the synthetic waveforms struggle to generate long-duration tsunami behaviour in simulations. Our tests suggest that the resolution of the bathymetry may be a potential factor affecting the simulated tsunami duration and amplitude. It should be noted that the maximum wave height in the records may occur after the decay of synthetic wave amplitudes. This implies that the inability to model long-duration tsunamis could result in underestimation in future tsunami hazard assessments.