The accurate and reliable numerical approximation of the hyperbolic wave equation is of fundamental importance to the simulation of acoustic, electromagnetic and elastic wave propagation phenomena. Here, we present families of variational space–time discretisation methods of higher order for the acoustic wave equation as a protoype model for the three-dimensional elastic wave equation. In tempo...