jats:titleAbstract</jats:title>jats:pThe Groningen gas reservoir, situated in the northeast of the Netherlands is western Europe's largest gas reservoir. Due to gas production measurable subsidence and seismicity has been detected across this region, attributed to the deformations induced by reservoir pore pressure depletion. We investigate the surface displacement history using a principal component analysis‐based inversion method to combine a diverse set of optical leveling, interferometric synthetic aperture radar, and Global Positioning System data to better constrain reservoir compaction and subsidence history. The generated compaction model is then used in combination with prior pressure depletion models to determine a reservoir uniaxial compressibility. The best fitting model of uniaxial compressibility is time independent but spatially variable. The absence of evidence for any significant time delay between changes in depletion and compaction rates supports an instantaneous poroelastic reservoir response. The absence of nonlinear yielding at the largest compaction strains suggests that anelastic deformations are a minor part of reservoir compaction.</jats:p>