The over-exploitation of groundwater has led to a significant drop in groundwater levels, which may lead to a series of geological disasters and ecological environmental problems such as ground subsidence and ground cracks. Therefore, through studying the dynamic change characteristics of groundwater, we can grasp the dynamic changes in groundwater level over time and invert the hydrogeological parameters, which provides an important basis for the management of groundwater resources. In this study, the confined aquifer III groundwater between 2005 and 2014 in Yancheng City was selected as the research object, and the Back Propagation (BP) neural network, Spatial-temporal Auto Regressive and Moving Average (STARMA) model, and BP-STARMA model were used to predict the spatial and temporal evolution trends of groundwater. In order to compare the prediction effectiveness of the BP-STARMA model, the fitting and prediction accuracies of the three models were measured from the perspectives of time and space. The results of the Relative Squared Error (RSE), Normal Mean Squared Error (NMSE), Root-Mean-Squared Error (RMSE), and Mean Absolute Error (MAE) were used to assess the robustness of the BP-STARMA model. The results showed that the fitting of the RMSE of BP-STARMA model was reduced by 39.92%, 38.35%, 30.25%, 31.55%, and 13.57% compared with the STARMA model, and by 22.2%, 8.7%, 15.9%, 28.5%, and 4.42% compared with the BP neural network model, respectively. Collectively, this shows that the BP-STARMA model has a better spatiotemporal prediction of groundwater level than the STARMA and BP neural network models, is more applicable to spatially continuous time-discrete spatiotemporal sequences, and is more applicable to spatiotemporal sequences that respond to natural geographic phenomena.