This paper presents a methodology for studying ground vibration in which the railway track is modelled in the space-time domain using the finite element method (FEM) and, for faster computation, discretisation of the ground using either FEM or the boundary element method (BEM) is avoided by modelling it in the wavenumber-frequency domain. The railway track is coupled to the ground through a series of rectangular strips located at the surface of the ground; their vertical interaction is described by a frequency-dependent dynamic stiffness matrix whose elements are represented by discrete lumped parameter models. The effectiveness of this approach is assessed firstly through frequency domain analysis using as excitation a stationary harmonic load applied on the rail. The interaction forces at the ballast/ground interface are calculated using the FE track through-ground coupling in the lumped parameter model, which has been found to be necessary for both track dynamics and ground vibration predictions.