Seasonal and interannual changes in the Earth's gravity field are mainly due to mass exchange among the atmosphere, ocean, and continental water sources. The terrestrial water storage changes, detected as gravity changes by the Gravity Recovery and Climate Experiment (GRACE) satellites, are mainly caused by precipitation, evapotranspiration, river transportation and downward infiltration processes. In this study, a land data assimilation system LDAS-G was developed to assimilate the GRACE terrestrial water storage (TWS) data into the Community Land Model (CLM3.5) using the POD-based ensemble four-dimensional variational assimilation method PODEn4DVar, disaggregating the GRACE large-scale terrestrial water storage changes vertically and in time, and placing constraints on the simulation of vertical hydrological variables to improve land surface hydrological simulations. The ideal experiments conducted at a single point and assimilation experiments carried out over China by the LDAS-G data assimilation system showed that the system developed in this study improved the simulation of land surface hydrological variables, indicating the potential of GRACE data assimilation in large-scale land surface hydrological research and applications. data assimilation, land surface model, terrestrial water storage, ensemble four-dimensional variational data assimilation method Citation:Sun Q, Xie Z H, Tian X J, et al. 2015. GRACE terrestrial water storage data assimilation based on the ensemble four-dimensional variational method