The potential for building damage and personal injury due to induced earthquakes in the Groningen gas field is being modeled in order to inform risk management decisions. To facilitate the quantitative estimation of the induced seismic hazard and risk, a ground motion prediction model has been developed for response spectral accelerations and duration due to these earthquakes that originate within the reservoir at 3 km depth. The model is consistent with the motions recorded from small-magnitude events and captures the epistemic uncertainty associated with extrapolation to larger magnitudes. In order to reflect the conditions in the field, the model first predicts accelerations at a rock horizon some 800 m below the surface and then convolves these motions with frequency-dependent nonlinear amplification factors assigned to zones across the study area. The variability of the ground motions is modeled in all of its constituent parts at the rock and surface levels.