Understanding swash zone dynamics is of crucial importance for coastal management as the swash motion, consisting of the uprush of the wave on the beach face and the subsequent downrush, is responsible for driving changes the beach morphology trough sediment exchanges between the sub-aerial and sub-aqueous beach. Improved understanding of the probabilistic characteristics of these motions has the potential to allow coastal engineers to develop improved sediment transport models which, in turn, can be further developed into coastal management tools. In this paper, novel descriptors of swash motions are obtained by combining field data and statistical modelling. Our results indicate that the probability distribution function (PDF) of shoreline height (p(ζ)) and trough-to-peak swash heights (p(ρ)) measured at a high energy, sandy beach were both inherently multimodal. Based on the observed multimodality of these PDFs, Gaussian Mixtures are shown to be the best method to statistically model them. Further, our results show that both offshore and surf zone dynamics are responsible for driving swash zone dynamics, which indicates unsaturated swash. The novel methods and results developed in this paper, both data collection and analysis, could aid coastal managers to develop improved swash zone models in the future.