Ecosystem state can be characterized by a set of attributes that are related to the ecosystem functionality, which is a relevant issue in understanding the quality and quantity of ecosystem services and goods, adaptive capacity and resilience to perturbations. This study proposes a major identification of Ecosystem Functional Types (EFTs) in Spain to characterize the patterns of ecosystem functional diversity and status, from several functional attributes as the Normalized Difference Vegetation Index (NDVI), Land Surface Temperature (LST) and Albedo. For this purpose, several metrics, related to the spatial variability in seasonal and annual patterns (e.g., relative range), have been derived from remote sensing time series of 1 km MODIS over the period 2000-2009. Moreover, precipitation maps from data provided by the AEMet (Agencia Estatal de Meteorología) and the corresponding aridity and humidity indices were also included in the analysis. To create the EFTs, the potential of the joint use of Kohonen's Self-Organizing Map (SOM) and the k-means clustering algorithm was tested. The EFTs were analyzed using different remote sensing (i.e., Gross Primary Production) and climatic variables. The relationship of the EFTs with existing land cover datasets and climatic data were analyzed through a correspondence analysis (CA