Not long after discovery of the Groningen field, gas-production-induced compaction and consequent land subsidence was recognised to be a potential threat to groundwater management in the province of Groningen, in addition to the fact that parts of the province lie below sea level. More recently, NAM's seismological model also pointed to a correlation between reservoir compaction and the observed induced seismicity above the field. In addition to the already existing requirement for accurate subsidence predictions, this demanded a more accurate description of the expected spatial and temporal development of compaction.Since the start of production in 1963, multiple levelling campaigns have gathered a unique set of deformation measurements used to calibrate geomechanical models. In this paper we present a methodology to model compaction and subsidence, combining results from rock mechanics experiments and surface deformation measurements. Besides the optical spirit-levelling data, InSAR data are also used for inversion to compaction and calibration of compaction models. Residual analysis, i.e. analysis of the difference between measurement and model output, provides confidence in the model results used for subsidence forecasting and as input to seismological models.