We developed a new global ice sheet model for simple climate model the Global Resolved Energy Balance (GREB) in the project for paleoclimate simulation. In the thesis, we described how the ISM was built and how climate and ice sheet components coupled to have globally coupled Earth system model - the GREB-ISM. After that, a series of benchmark experiments are presented to evaluate the performance of the GREB-ISM in long-term mean state and transition simulations. Last but not least, we applied the GREB-ISM to understand the importance and relevant physics of climate-ice sheet feedback by using a series of sensitivity experiments.