A comprehensive estimation of water budget components, particularly groundwater storage (GWS) and fluxes, is crucial. In this study, we evaluate the terrestrial water budget of the Donga basin (Benin, West Africa), as simulated by three land surface models (LSMs) used in the African Monsoon Multidisciplinary Analysis Land Surface Model Intercomparison Project, phase 2 (ALMIP2): CLM4, Catchment LSM (CLSM), and Minimal Advanced Treatments of Surface Interaction and Runoff (MATSIRO). All three models include an unconfined groundwater component and are driven by the same ALMIP2 atmospheric forcing from 2005 to 2008. Results show that all three models simulate substantially shallower water table depth (WTD) with smaller seasonal variations, approximately 1-1.5 m compared to the observed values that range between 4 and 9.6 m, while the seasonal variations of GWS are overestimated by all the models. These seemingly contradictory simulation results can be explained by the overly high specific yield prescribed in all models. All models achieve similar GWS simulations but with different fractions of precipitation partitioning into surface runoff, base flow, and evapotranspiration (ET), suggesting high uncertainty and errors in the terrestrial and groundwater budgets among models. The poor performances of models can be attributed to bias in the hydrological partitioning (base flow vs surface runoff) and sparse subsurface data. This analysis confirms the importance of subsurface hydrological processes in the current generation of LSMs and calls for substantial improvement in both surface water budget (which controls groundwater recharge) and the groundwater system (hydrodynamic parameters, vertical geometry).