According to the phenological data from Chinese Phenology Observation Network of Chinese Academy of Sciences, we compiled the natural calendar for 1985-2012 with 21 plant species and 99 phenophases in Harbin, Heilongjiang Province. Comparing the calendar with the original one for 1963-1984, the research revealed the variation characteristics of 99 phenophases. Correlation analysis and regression analysis were used to examine the relationships between changes of phenophases and climatic drivers. Since 1985, timing of phenological spring (represented by the timing of bud expansion of Ulmus pumila), summer (represented by the timing of 50% of full flowering of Syringa reticulate) and autumn (represented by the timing of fruit maturity of Lonicera maackii) have been advanced by 7 days, 6 days and 19 days respectively, while timing of phenological winter (represented by the timing of end of leaf fall of Juglans mandshurica) has been delayed by 2 days. Meanwhile, compared with the original calendar, the average dates of phenophases have been advanced by 3 to 11 days in spring, summer and autumn, but delayed by 3 days in winter. The earliest date of phenology showed advances mainly in all seasons, while the latest dates of phenology were delayed in summer and winter. The order of some phenophases in phenological seasons changed with different degrees. Temperature changes before the majority of phenophases is probably the main reason for the changes of phenological season in the last 30 years. Different sensibility of different species and phenophases may result in the change of phenophases order in seasons. The result that the phenological events in spring have been advanced in accordance with many other studies at home and abroad, showing the responses of ecosystems to a warming atmosphere.