The classification of plant associations and their mapping play a key role in defining habitat biodiversity management, monitoring, and conservation strategies. In this work we present a methodological framework to map Mediterranean forest plant associations and habitats that relies on the application of the Functional Principal Component Analysis (FPCA) to the remotely sensed Normalized Difference Vegetation Index (NDVI) time series. FPCA, considering the chronological order of the data, reduced the NDVI time series data complexity and provided (as FPCA scores) the main seasonal NDVI phenological variations of the forests. We performed a supervised classification of the FPCA scores combined with topographic and lithological features of the study area to map the forest plant associations. The supervised mapping achieved an overall accuracy of 87.5%. The FPCA scores contributed to the global accuracy of the map much more than the topographic and lithological features. The results showed that (i) the main seasonal phenological variations (FPCA scores) are effective spatial predictors to obtain accurate plant associations and habitat maps; (ii) the FPCA is a suitable solution to simultaneously express the relationships between remotely sensed and ecological field data, since it allows us to integrate these two different perspectives about plant associations in a single graph. The proposed approach based on the FPCA is useful for forest habitat monitoring, as it can contribute to produce periodically detailed vegetation-based habitat maps that reflect the “current” status of vegetation and habitats, also supporting the study of plant associations.

The main purpose of the 92/43/EEC Habitats Directive is to contribute to the conservation of biodiversity, understood as habitat types and species of the flora and fauna of the European Union. To achieve this goal, natural and semi-natural biodiversity as a whole must be recognized and included in its annexes. As for the conservation of biotopes, named habitat types, Italy is unfortunately lacking as it the Annex I does not include important ecosystems that are typical of its territory, rare for biogeographical reasons or threatened. Therefore, the opportunity to identify a first list of significant habitats for central Italy is discussed here. For each of the new proposed types (new habitats or new subtypes) a sheet has been prepared to highlight their salient characteristics. The new proposals concern seven habitat types and one subtype: sedge and reeds formations (Freshwater large sedge and reed beds), willow shrublands (Shrubby willow formations of river banks and fens), Apennine garrigues (Apennine hilly and montane garrigues), a new subtype of Annex I Habitat 6130 (Communities of herbaceous and dwarf shrub-suffrutescent plants of Italian ultramafic substrates), ancient olive groves ("Centuries-old olive groves" with evergreen Quercus spp. and arborescent matorral), secondary meadows (Italian submontane and montane pastured meadows dominated by Cynosurus cristatus), badlands (Pioneer halophilous and sub-halophilous communities of “calanchi” and “biancane” badlands) and hop-hornbeam woods (Italian-Balkan hop-hornbeam woods).

The aim of the present study is to demonstrate the value of some nitrophilous plants as bioindicators related to global warming. As a case study, the space–time dynamics of the populations of some nitrophilous and ruderal species have been investigated along the southern-central Adriatic side of the Italian Peninsula. These have been examined according to the range of their distribution areas, to their ecological needs, and to the availability of past data, which have then been related to the data on global warming across the same territories. The choice for this investigation was for nitrophilous species with a Mediterranean distribution – of the Stenomediterranean type – and for some recent entries into the local flora as alien species. The spread of these species occurs in areas with intense human activity, where they have exploited the conditions of the warmer niche, such that their presence and their spatial spread observed over time are clearly linked to global warming. For all of the species in question, rapid increases in the population numbers have been observed, along with a northward shift of their distribution areas. These changes correspond to the increase in average annual temperature as revealed by the thermometric measurements.

Phytosociology is a reference method to classify vegetation that relies on field data. Its classification in hierarchical vegetation units, from plant associations to class level, hierarchically reflects the floristic similarity between different sites on different spatial scales. The development of remotely sensed multispectral platforms as satellites enormously contributes to the detection and mapping of vegetation on all scales. However, the integration between phytosociology and remotely sensed data is rather difficult and little practiced despite being a goal for the modern science of vegetation. In this study, we demonstrate how normalized difference vegetation index (NDVI) time series with functional principal component analysis (FPCA) could support the analyses of phytosociologists. The approach supports the recognition and characterization of forest plant communities identified on the ground by the phytosociological approach by using NDVI time series that encode phenological behaviors. The methodology was evaluated in two study areas of central Italy, and it could characterize and discriminate six different forest plant associations that have similar dominant tree species but distinct specific composition: three dominated by black hornbeam (Ostrya carpinifolia) and three dominated by holm oak (Quercus ilex). The methodology was also able to optimize the ground data collection of unexplored areas (from a phytosociological point of view) by using a phenoclustering approach. The obtained results confirmed that by using remote sensing, it is possible to separate and distinguish plant communities in an objective/instrumental way, thus overcoming the subjectivity intrinsic to the phytosociological method. In particular, FPCA functional components (NDVI seasonalities) were significantly correlated with vegetation abundance data variation (Mantel r = 0.76, p < 0.001).

The mapping and monitoring of natural and semi-natural habitats are crucial activities and are regulated by European policies and regulations, such as the 92/43/EEC. In the Mediterranean area, which is characterized by high vegetational and environmental diversity, the mapping and monitoring of habitats are particularly difficult and often exclusively based on in situ observations. In this scenario, it is necessary to automate the generation of updated maps to support the decisions of policy makers. At present, the availability of high spatiotemporal resolution data provides new possibilities for improving the mapping and monitoring of habitats. In this work, we present a methodology that, starting from remotely sensed time-series data, generates habitat maps using supervised classification supported by Functional Data Analysis. We constructed the methodology using Sentinel-2 data in the Mediterranean Special Area of Conservation “Gola di Frasassi” (Code: IT5320003). In particular, the training set uses 308 field plots with 11 target classes (five forests, two shrubs, one grassland, one mosaic, one extensive crop, and one urban land). Starting from vegetation index time-series data, Functional Principal Component Analysis was applied to derive FPCA scores and components. In particular, in the classification stage, the FPCA scores are considered as features. The obtained results out-performed a previous map derived from photo-interpretation by domain experts. We obtained an overall accuracy of 85.58% using vegetation index time-series, topography, and lithology data. The main advantages of the proposed approach are the capability to efficiently compress high dimensional data (dense remote-sensing time series) providing results in a compact way (e.g., FPCA scores and mean seasonal time profiles) that: (i) facilitate the link between remote sensing with habitat mapping and monitoring and their ecological interpretation and (ii) could be complementary to species-based approaches in plant community ecology and phytosociology.

The identification, classification and mapping of different plant communities and habitats is of fundamental importance for defining biodiversity monitoring and conservation strategies. Today, the availability of high temporal, spatial and spectral data from remote sensing platforms provides dense time series over different spectral bands. In the case of supervised mapping, time series based on classical vegetation indices (e.g., NDVI, GNDVI, …) are usually input characteristics, but the selection of the best index or set of indices (which guarantees the best performance) is still based on human experience and is also influenced by the study area. In this work, several different time series, based on Sentinel-2 images, were created exploring new combinations of bands that extend the classic basic formulas as the normalized difference index. Multivariate Functional Principal Component Analysis (MFPCA) was used to contemporarily decompose the multiple time series. The principal multivariate seasonal spectral variations identified (MFPCA scores) were classified by using a Random Forest (RF) model. The MFPCA and RF classifications were nested into a forward selection strategy to identify the proper and minimum set of indices’ (dense) time series that produced the most accurate supervised classification of plant communities and habitat. The results we obtained can be summarized as follows: (i) the selection of the best set of time series is specific to the study area and the habitats involved; (ii) well-known and widely used indices such as the NDVI are not selected as the indices with the best performance; instead, time series based on original indices (in terms of formula or combination of bands) or underused indices (such as those derivable with the visible bands) are selected; (iii) MFPCA efficiently reduces the dimensionality of the data (multiple dense time series) providing ecologically interpretable results representing an important tool for habitat modelling outperforming conventional approaches that consider only discrete time series.