The present study classified global Ecosystem Functional Types (EFTs) derived from seasonal vegetation dynamics of the GIMMS3g NDVI time-series. Rotated Principal Component Analysis (PCA) was run on the derived phenological and productivity variables, which selected the Standing Biomass (approximation of Net Primary Productivity), the Cyclic Fraction (seasonal vegetation productivity), the Permanent Fraction (permanent surface vegetation), the Maximum Day (day of maximum vegetation development) and the Season Length (length of vegetation growing season) variables, describing 98% of the variation in global ecosystems. EFTs were created based on Isodata classification of the spatial patterns of the Principal Components and were interpreted via gradient analysis using the selected remote sensing variables and climatic constraints (radiation, temperature, and water) of vegetation growth. The association of the EFTs with existing climate and land cover classifications was demonstrated via Detrended Correspondence Analysis (DCA). The ordination indicated good description of the global environmental gradient by the EFTs, supporting the understanding of phenological and productivity dynamics of global ecosystems. Climatic constraints of vegetation growth explained 50% of variation in the phenological data along the EFTs showing that part of the variation in the global phenological gradient is not climate related but is unique to the Earth Observation derived variables. DCA demonstrated good correspondence of the EFTs to global climate and also to land use classification. The results show the great potential of Earth Observation derived parameters for the quantification of ecosystem functional dynamics and for providing reference status information for future assessments of ecosystem changes.

Observing trends in global ecosystem dynamics is an important first step, but attributing these trends to climate variability represents a further step in understanding Earth system changes. In the present study, we classified global Ecosystem Response Types (ERTs) based on common spatio-temporal patterns in time-series of Standardized Precipitation Evapotranspiration Index (SPEI) and FPAR3g anomalies (1982–2011) by using an extended Principal Component Analysis. The ERTs represent region specific spatio-temporal patterns of ecosystems responding to drought or ecosystems with decreasing severity in drought events as well as ecosystems where drought was not a dominant factor in a 30-year period. Highest explanatory values in the SPEI12-FPAR3g anomalies and strongest SPEI12-FPAR3g correlations were seen in the ERTs of Australia and South America whereas lowest explanatory value and lowest correlations were observed in Asia and North America. These ERTs complement traditional pixel based methods by enabling the combined assessment of the location, timing, duration, frequency and severity of climatic and vegetation anomalies with the joint assessment of wetting and drying climatic conditions. The ERTs produced here thus have potential in supporting global change studies by mapping reference conditions of long term ecosystem changes.