We have applied the empirical Green’s function (EGF) method to 53 pairs of earthquakes, with magnitudes ranging from M = 0.4 to M = 3.4, induced by gas production from the Groningen field in the Netherlands. For a subset of the events processed, we find that the relative source time functions obtained by the EGF deconvolution show clear indications of a horizontal component of rupture propagation. The earthquake monitoring network used has dense azimuthal coverage for nearly all events such that wavelet duration times can be picked as a function of source-station azimuth and inverted using the usual Doppler broadening model to estimate rupture propagation strike, distance, and velocity. Average slip velocities have also been estimated and found to be in agreement with typical published values. We have used synthetic data, from both a simple convolutional model of the seismogram and more sophisticated finite difference rupture simulations, to validate our data processing workflow and develop kinematic models which can explain the observed characteristics of the field data. Using a measure based on the L1-norm to discriminate results of differing quality, we find that the highest quality results show very good alignment of the rupture propagation with directions of the detailed fault map, obtained from the full-field 3D seismic data. The dip direction rupture extents were estimated from the horizontal rupture propagation distances and catalogue magnitudes showing that, for all but the largest magnitude event (the M = 3.4 event of 8th January 2018), the dip-direction extent is sufficiently small to be contained wholly within the reservoir.