AbstractOne critical challenge of exploring flora phenology is on characterizing ecosystem phenological diversity (EPD), and thus how EPD’s performance is influenced by climate changes has also been an open macro-ecological question. To fill these two gaps, we proposed an innovative method for reflecting EPD, by taking the advantage of the often-classified inverse factor of spatial resolution discrepancy between the used remote sensing datasets of vegetation phenological dates (green-up and brown-up) and snow cover phenological dates (SPDs) (onset and end) around the Arctic, and further, we examined the cross response/feedbacks of the two kinds of EPDs to the two categories of SPDs. We found that the circumpolar green-up and brown-up EPDs both were shrinking, driven more by the delaying of the onset SPDs than the advancing of the end SPDs; North America and North Eurasia performed with inconsistent EPD response/feedbacks to the related SPD anomalies; and further, the EPD-SPD response/feedbacks in some locations exhibited the time-lag effect, e.g., the green-up EPDs made the strongest response to the onset SPDs of two years earlier. Overall, the validated method and the new findings are of implications for improving the phenology modules in Earth system models, and the contributions of the present study have enlightening significance for kicking off the new EPD branch in macrosystem phenological ecology.