A nonlinear, fractional, surface wave equation with a spatial derivative of second order was developed by Kappler, Shrivastava, Schneider, and Netz [Phys. Rev. Fluids 2, 114804 (2017)] for propagation along an elastic interface coupled to a viscous incompressible liquid. Linear theory for the attenuation and dispersion was developed originally by Lucassen [Trans. Faraday Soc. 64, 2221 (1968)]. Kappler et al. introduced a fractional time derivative to account for the Lucassen wave attenuation and dispersion, and they included quadratic and cubic nonlinearity associated with compression of the elastic interface. Presented here is an integrated form of their time domain equation for progressive waves that is first order in the spatial derivative. Solutions of this evolution equation capture the main features of waveforms predicted by the full model equation of Kappler et al., especially the formation and propagation of shocks, while the evolution equation can be solved numerically with substantially less computational cost. Approximate analytical expressions obtained from the evolution equation for the nonlinear propagation speed and attenuation of a compression pulse reveal that a threshold phenomenon discussed by Kappler et al. is due to competition between quadratic and cubic nonlinearity associated with a lipid monolayer interface.