Identification of complete drivers for phenology changes is crucial for developing prediction models of plant phenology. In addition to climatic factors, the interaction among phenological events has recently been reported as an important driver for the phenology changes of forests, savannas, and grasslands. However, open questions remain as to whether the phenological interaction exists in agricultural ecosystems, among which winter wheat plays a vital role in feeding human beings. In this study, we investigated the interaction among the phenological events of winter wheat in the North China Plain (NCP) using both field and satellite data. Considering the large discrepancies between the existing satellite estimation and field measurements of winter wheat phenology, we first improved the MODIS-based estimation of green-up date (GUD), heading date (HD), and maturity date (MD) through a re-calibrated relative threshold method (RTM) in the NCP. The GUD, HD, and MD were accurately estimated with the mean absolute errors (MAE) and root mean squared errors (RMSE) lower than 7.5 days, compared with the RMSEs ranging from 12.0 to 36.1 days in previous studies. Then, the relationships among the GUD, HD, and MD were analyzed using the field data collected at agricultural meteorological stations. The GUD (HD) showed a significantly positive correlation with the HD (MD). Quantitatively, a one-day earlier GUD (HD) would result in an earlier HD (MD) of 0.57 days (0.60 days). Furthermore, we applied the partial correlation analysis to the improved MODIS estimation of GUD, HD, and MD to investigate their interactions by considering the simultaneous influences from climatic factors. The results showed that the HD (MD) with 85.2% (94.5%) of all winter wheat pixels presented a significantly positive correlation with the GUD (HD). Meanwhile, the GUD (HD) with 84.2% (33.3%) of the entire winter wheat area presented a significantly negative correlation with pre-season temperature. These results suggest that both the climatic factors and phenological interactions should be included in the future development of winter wheat phenology models to improve the prediction accuracies. on NDVI, a novel snow-free type of VIs, including the empirical normalized difference phenology index (NDPI,[27]) and the semi-analytical normalized difference greenness index (NDGI,[26]), were recently proposed. Based on the time-series of the computed VIs, vegetation phenology can be estimated by extracting the day of year (DOY) corresponding to pre-determined absolute or relative thresholds of the VI values [1,15] or the maximum change rate of the trajectories of VIs [11,28]. For winter wheat, the GUDs were previously estimated based on either a relative threshold of 20% or a maximum change rate of the fitted curve of the time series VIs [1,2,10,11]. These studies mainly focused on changing trends of winter wheat phenology and its responses to climatic factors from the perspective of regarding the wheat fields as important agricultural ecosystems to the global...