Soil moisture plays an important role in water and energy balance in land-atmospheric interaction, but is impacted directly by the groundwater table. Dynamic variation of the groundwater table can be described mathematically by a moving boundary problem. In this paper, the moving boundary problem is reduced to a fixed boundary problem through a coordinate transformation. A new model of groundwater table simulation is developed using the mass-lumped finite element method and is coupled with the land surface model of Variable Infiltration Capacity (VIC). The simulation results show that the new model not only can simulate the groundwater table dynamically, but also can evade the choice of water table depth scale in computation with a low computation cost.