Global land cover types in 2001 and 2010 were mapped at 250m resolution with multiple year time series Moderate Resolution Imaging Spectrometer (MODIS) data. The map for each single year was produced not only from data of that particular year but also from data acquired in the preceding and subsequent years as temporal context. Slope data and geographical coordinates of pixels were also used. The classification system was derived from the finer resolution observation and monitoring of global land cover (FROM-GLC) project. Samples were based on the 2010 FROM-GLC project and samples for other years were obtained by excluding those changed from 2010. A random forest classifier was used to obtain original class labels and to estimate class probabilities for 2000–2002, and 2009–2011. The overall accuracies estimated from cross validation of samples are 74.93% for 2001 and 75.17% for 2010. The classification results were further improved through post processing. A spatial-temporal consistency model, Maximum a Posteriori Markov Random Fields (MAP-MRF), was first applied to improve land cover classification for each 3 consecutive years. The MRF outputs for 2001 and 2010 were then processed with a rule-based label adjustment method with MOD44B, slope and composited EVI series as auxiliary data. The label adjustment process relabeled the over-classified forests, water bodies and barren lands to alternative classes with maximum probabilities.