The land subsidence due to groundwater exploitation has an obvious hysteretic nature with respect to the decrease of the under groundwater level,and the uneven settlement often causes ground fissures.To study these important features,a visco-elastic plastic constitutive relationship with consideration of the coupling of seepage and soil deformation is proposed,and a finite element model with variable coefficients based on the Biot’s consolidation theory is built.With the groundwater exploitation and the land subsidence control in Cangzhou City,Hebei Province as an example,the variations of the under groundwater level and the development of the land subsidence due to the groundwater exploitation are simulated and ground fissures are predicted by the horizontaldisplacement calculation.The results show that the lag time between the land subsidence and the under groundwater level descent is about a month,and the simulated results of fissures agree well with the observed data.The model can well reveal the characterization of the interaction between the land subsidence and the groundwater exploitation.