The mechanical models for earthquake fault dislocations and volcano magma chamber inflation or deflation and dike injection can be used to give some idea about subsidence due to water or other fluid extraction or injection. They will not be completely accurate, though, because the volume change due to fluid poroelastic effects is generally over a volume, not a thin fault plane or dike surface. Paul Segall, the book author, and his student used arrays of Mogi volume change sources to model fluid extraction ( Mossop, A., and P. Segall (1997), Subsidence at The Geysers geothermal field, N. California from a comparison of GPS and leveling surveys, Geophysical Research Letters , 24 (14), 1839–1842. ). I used a similar approach for modeling volume change at shallow depths for Bam postseismic deformation:
Fielding, E. J., P. R. Lundgren, R. Bürgmann, and J. Funning Gareth (2009), Shallow fault-zone dilatancy recovery after the 2003 Bam earthquake in Iran, Nature , 458 (7234), 64-68, doi:10.1038/nature07817.