Phase-resolving numerical models are frequently used tools to investigate short and long wave transformation, nonlinear wave interactions, and wave runup. Moreover, nearshore morphodynamics can be explored with the recent advancement of the models and computational resources. Sea surface elevation time series that force phase-resolving models at the offshore boundary are often unavailable. Therefore, time series are usually recreated from wave energy-frequency spectra through the superposition of harmonics. The wave phases of the harmonics are unknown and therefore assumed to be randomly distributed. This implies that an infinite number of time series with different sequencing of waves can be recreated from a single wave-energy spectrum and, for that reason, recreated time series are a source of uncertainty in model predictions. This intrinsic uncertainty has been found to cause variability in wave overtopping of structures (e.g., Pearson et al, 2002; Williams et al., 2014; Romano et al., 2015) and in setup and runup at beaches (McCabe et al., 2011; Torres-Freyermuth et al., 2019). Torres-Freyermuth et al. (2019) investigated the effect of intrinsic uncertainty on runup at planar beaches for different wave conditions and beach slopes and suggested that uncertainty is especially important under dissipative conditions. Yet unknown is the effect of intrinsic uncertainty on bed evolution. Here we assess the effect of intrinsic uncertainty on inner surf and swash zone evolution at three beaches with different beach morphology.Recorded Presentation from the vICCE (YouTube Link): https://youtu.be/w3zi3Yoo170