Motivated by the significant topographic influence on land surface processes, this study explored two methods to discretize watersheds into two types of subgrid structures to capture spatial heterogeneity for land surface models. Adopting geomorphologic concepts in watershed discretization yields improved capability in capturing subgrid topographic heterogeneity, which also allowed climatic and land cover variability to be better represented with a nominal increase in computational requirements.