Accurate groundwater depth forecasting is particularly important for human life and sustainable groundwater management in arid and semi-arid areas. To improve the groundwater forecasting accuracy, in this paper, a hybrid groundwater depth forecasting model using configurational entropy spectral analyses (CESA) with the optimal input is constructed. An original groundwater depth series is decomposed into subseries of different frequencies using the variational mode decomposition (VMD) method. Cross-correlation analysis and Shannon entropy methods are applied to select the optimal input series for the model. The ultimate forecasted values of the groundwater depth can be obtained from the various forecasted values of the selected series with the CESA model. The applicability of the hybrid model is verified using the groundwater depth data from four monitoring wells in the Xi'an of Northwest China. The forecasting accuracy of the models was evaluated based on the average relative error (RE), root mean square error (RMSE), coefficient of determination (R ) and Nash-Sutcliffe coefficient (NSE). The results indicated that comparing with the CESA and Autoregressive model, the hybrid model has higher prediction performance. This article is protected by copyright. All rights reserved.