Vegetation plays an important role in atmospheric, hydrologic and biochemical cycles and is an important indicator of the impact of climate and human factors on the environment. In this paper, a method, which combines the empirical orthogonal function (EOF) and temporal unmixing analysis (TUA) methods, is applied to monitor the phenological characteristcs and spatial distribution of vegetation phenology in the middle part of the Huai River region. Based on the variance and EOF curves, the EOF provides the number of phenology modes, information which is the basis for an accurate temporal unmixing model. The TUA describes the temporal vegetation phenological details and spatial distribution. Importantly, this approach does not require assumptions, prior information or pre-defined thresholds. The vegetation phenology curves derived from the MODIS EVI data using the combined EOF and TUA methods display much more detail than the curves from Landsat TM using spectral mixture analysis (SMA). Additionally, the vegetation phenology spatial distribution from MODIS EVI is consistent with the field survey data. The combination method of EOF and TUA can be used to monitor vegetation phenology spatiotemporal change in a large area from time series of MODIS EVI data.