Exploring Long-Term Anomalies in the Vegetation Cover of Peri-Urban Parks Using the Fisher-Shannon Method

The main goal of this study was to evaluate the potential of the Fisher-Shannon statistical method applied to the MODIS satellite time series to search for and explore any small multiyear trends and changes (herein also denoted as inner anomalies) in vegetation cover. For the purpose of our investigation, we focused on the vegetation cover of three peri-urban parks close to Rome and Naples (Italy). For each of these three areas, we analyzed the 2000–2020 time variation of four MODIS-based vegetation indices: evapotranspiration (ET), normalized difference vegetation index (NDVI), leaf area index (LAI), and enhanced vegetation index (EVI). These data sets are available in the Google Earth Engine (GEE) and were selected because they are related to the interactions between soil, water, atmosphere, and plants. To account for the great variability exhibited by the seasonal variations while identifying small multiyear trends and changes, we devised a procedure composed of two steps: (i) application of the Singular Spectrum Analysis (SSA) to each satellite-based time series to detect and remove the annual cycle including the seasonality and then (ii) analysis of the detrended signals using the Fisher-Shannon method, which combines the Shannon entropy and the Fisher Information Measure (FIM). Our results indicate that among all the three pilot test areas, Castel Volturno is characterized by the highest Shannon entropy and the lowest FIM that indicate a low level of order and organization of vegetation time series. This behaviour can be linked to the degradation phenomena induced by the parasite (Toumeyella parvicornis) that has affected dramatically the area in recent years. Our results were nicely confirmed by the comparison with in situ analyzed and independent data sets revealing the existence of subtle, small multiyear trends and changes in MODIS-based vegetation indices.


Introduction
Natural capital (NC) through the ecosystem process provides ecosystem services that are vital and critical to the functioning of the Earth's life-support system, such as (but not only) air, water, fertile soil, pollination, and hazard protection. The concept of NC and associated ecosystem services reflects a recognition that environmental systems are fundamental not only for providing resources and services but also contributing to economic outputs and social well-being. The biophysical evaluation of the ecosystem then leads directly to the ecological and monetary evaluation, since it depends on its state of conservation. Nevertheless, if any form of vegetation cover (cropland, grassland, forest, etc.) provides numerous ecosystem services the estimation of the status and trends of natural capital poses critical challenges due to the diversity of environmental assets, stocks, Entropy 2022, 24, 1784 3 of 17 spiration, and taking precedent (prior month) conditions into account [15][16][17][18]. For the calculation of the PDSI, four inputs are needed: precipitation, temperature, latitude, and the soil's available water capacity (AWC) of the study area, which is a constant also known as the field capacity [19]. These four inputs are used to compute a water balance for the study area, which then serves as the basis for the calculation of the PDSI. For a detailed explanation of the calculation of the PDSI, refer to [20].
The main goal of this study was to evaluate the potential of the Fisher-Shannon statistical method to explore any anomalies happening for the vegetation cover around large urban areas using soil-water-atmosphere-plant-related satellite products available in the Google Earth Engine cloud database (i.e., LAI, NDVI, EVI, and ET from MODIS). Periurban parks were selected for our investigations as particularly significant areas because they play a key role not only in improving environmental quality and life but also in facing climatic change and mitigating climate change effects.

Study Areas and Dataset
For the purpose of this study, the following three study areas in Italy were selected: Appia Park and Castel Porziano in the center, and Castel Volturno in the south ( Figure 1). These areas were selected because they are representative of diverse vegetation covers, as detailed in the following Sections 2. 1-2.3. across different regions is the Palmer Drought Severity Index (PDSI), originally developed by Palmer [14]. Various studies show that PDSI is very effective in determining long-term drought, considering the basic effects of global warming through potential evapotranspiration, and taking precedent (prior month) conditions into account [15][16][17][18]. For the calculation of the PDSI, four inputs are needed: precipitation, temperature, latitude, and the soil's available water capacity (AWC) of the study area, which is a constant also known as the field capacity [19]. These four inputs are used to compute a water balance for the study area, which then serves as the basis for the calculation of the PDSI. For a detailed explanation of the calculation of the PDSI, refer to [20].
The main goal of this study was to evaluate the potential of the Fisher-Shannon statistical method to explore any anomalies happening for the vegetation cover around large urban areas using soil-water-atmosphere-plant-related satellite products available in the Google Earth Engine cloud database (i.e., LAI, NDVI, EVI, and ET from MODIS). Periurban parks were selected for our investigations as particularly significant areas because they play a key role not only in improving environmental quality and life but also in facing climatic change and mitigating climate change effects.

Study Areas and Dataset
For the purpose of this study, the following three study areas in Italy were selected: Appia Park and Castel Porziano in the center, and Castel Volturno in the south ( Figure 1). These areas were selected because they are representative of diverse vegetation covers, as detailed in the following Sections 2.1-2.3.  The characteristics of the study areas (including longitude, latitude, area, annual precipitation, annual mean temperature, vegetation description, and climate system) are presented in Table 1.

Castel Volturno
Castel Volturno is a natural reserve which occupies a total area of 268 hectares; it extends along the sandy coast of the municipality of Castel Volturno (CE), in a strip between the mouth of the Regi Lagni to the north and the mouth of Lago Patria to the South ( Figure 2). The site includes the protected area named ZSC IT8010021 "Pineta di Patria" and the Regional Nature Reserve "Foce Volturno-Costa di Licola" made up of pines.

Castel Porziano
The Presidential Estate of Castelporziano is about 25 km from the center of Rome ( Figure 3) and covers an area of 60 km 2 (6039 hectares), consisting of humid areas behind the dunes and areas with low and high scrub with the typical evergreen and aromatic species. Most of the extension is occupied by lowland hygrophilous wood (lowland wood linked to humid environments), characterized by the presence of evergreen and deciduous oaks and by more purely hygrophilous species, near the wetlands. The

Castel Porziano
The Presidential Estate of Castelporziano is about 25 km from the center of Rome ( Figure 3) and covers an area of 60 km 2 (6039 hectares), consisting of humid areas behind the dunes and areas with low and high scrub with the typical evergreen and aromatic species. Most of the extension is occupied by lowland hygrophilous wood (lowland wood linked to humid environments), characterized by the presence of evergreen and deciduous oaks and by more purely hygrophilous species, near the wetlands. The peculiarity of Castelporziano is above all linked to the interpenetration of the oak grove typical of the Mediterranean climate and the oak grove typical of the continental climate. Among the evergreen oaks, the holm oak, the cork oak, and the crenata oak, a hybrid between turkey oak and cork oak, are widely diffused. Among the deciduous oaks we note the turkey oak, the English oak, and the farnetto, while in the cooler wetlands we can find poplar, ossifillo ash, maple, hornbeam, and Oriental hornbeam typical of Mediterranean coastal environments. The wood (mixed plain), one of the most delicate ecosystems to be protected, extends for about 2300 hectares; the Mediterranean scrub environments, low and high, cover an area of about 500 hectares; the holm oak occupies an area of 261 hectares above all in the back dune area; and the cork oak wood covers an area of about 460 hectares. The woods alternate between clearings and natural grasslands, forming plant associations of great environmental variety. The stone pine forests, created with artificial reforestation, extend for about 750 hectares with the purpose of consolidating the sandy dunes and protecting the rear dunes from sea winds.

The Appia Antica Regional Park
The Appia Antica Regional Park, with its 4580 hectares, is the largest urban protected area in Europe. A green wedge runs from the city center towards the Castelli Romani ( Figure 4). This green wedge, a vast 4580 hectares (following the last extension in October 2018), is characterized by various areas of interest: the Via Appia Antica and its adjacencies; the Caffarella Valley; the archaeological area of the Via Latina and of the Aqueducts; the Tenuta di Tormarancia; the Tenuta Farnesiana; and the areas of Divino Amore, Falcognana and Mugilla. The park is so vast that it affects three municipalities: those of Rome, Ciampino, and Marino. The Park today looks like a mosaic of different environments: large spaces intended for cultivation and extensive grazing are interrupted by uncultivated areas; residual wooded strips, where agricultural exploitation has not arrived or has long since ceased; ditches with the presence of riparian vegetation; and

The Appia Antica Regional Park
The Appia Antica Regional Park, with its 4580 hectares, is the largest urban protected area in Europe. A green wedge runs from the city center towards the Castelli Romani ( Figure 4). This green wedge, a vast 4580 hectares (following the last extension in October 2018), is characterized by various areas of interest: the Via Appia Antica and its adjacencies; the Caffarella Valley; the archaeological area of the Via Latina and of the Aqueducts; the Tenuta di Tormarancia; the Tenuta Farnesiana; and the areas of Divino Amore, Falcognana and Mugilla. The park is so vast that it affects three municipalities: those of Rome, Ciampino, and Marino. The Park today looks like a mosaic of different environments: large spaces intended for cultivation and extensive grazing are interrupted by uncultivated areas; residual wooded strips, where agricultural exploitation has not arrived or has long since ceased; ditches with the presence of riparian vegetation; and some wet areas. These seminatural environments and the agricultural contexts now represent the agroecosystem of the Roman countryside. It is a system of considerable naturalistic and scientific interest, due to the presence of wildlife communities and plant associations, consistent with the ecological potential of the area. The Appia Antica Park is a substantial part of the Ecological Network of the city of Rome and is the most important protected periurban area of the Lazio Region.

Data Sets
Four different satellite products available in the cloud storage of Google Earth Engine were selected because they are related to the interactions between soil, water, atmosphere, and plants. For the purpose of this study, MODIS products were chosen due to their global coverage and long duration of data acquisition. The supplementary file summarizes the characteristics of the studied datasets. The area-averaged time-series of studied parameters were extracted for the polygons representing the study regions (see  using the GEE JavaScript API for the common period of 2001-2020.

The Singular Spectrum Analysis
There are several techniques for decomposing a time series into a certain number of independent components; among these, Singular Spectrum Analysis (SSA) [21] represents an efficient and well-known decompositional method that is based on phase-lagged copies of the series. The independent components obtained by applying the SSA can be easily recognizable as slowly changing trend, oscillatory components, and structureless noise [22].
Let us consider a time series yi (i = 1, …, N) and a lag M; then the Toeplitz lagged

Data Sets
Four different satellite products available in the cloud storage of Google Earth Engine were selected because they are related to the interactions between soil, water, atmosphere, and plants. For the purpose of this study, MODIS products were chosen due to their global coverage and long duration of data acquisition. The Supplementary File summarizes the characteristics of the studied datasets. The area-averaged time-series of studied parameters were extracted for the polygons representing the study regions (see Figures 1-4) using the GEE JavaScript API for the common period of 2001-2020.

The Singular Spectrum Analysis
There are several techniques for decomposing a time series into a certain number of independent components; among these, Singular Spectrum Analysis (SSA) [21] represents an efficient and well-known decompositional method that is based on phase-lagged copies of the series. The independent components obtained by applying the SSA can be easily recognizable as slowly changing trend, oscillatory components, and structureless noise [22].
Let us consider a time series y i (i = 1, . . . , N) and a lag M; then the Toeplitz lagged correlation matrix can be constructed: Sorting its eigenvalues λ k in decreasing order, the corresponding eigenvectors E k,j where j and k vary from 1 to M are used to calculate the k-th principal component i and the k-th reconstructed component of the time series: Since the eigenvalue λ k represents the fraction of the total variance of the original series explained in the k-th reconstructed component R k , the decreasing order of the eigenvalues also reflects the decreasing order of the reconstructed components by the fraction of the total variance of the series [23]. SSA requires that the lag M is properly selected. Khan and Poskitt [24] The minimum description length (MDL) criterion [25]: is used to separate the series into two parts that we can define as trend and detrended series; λ k are the eigenvalues, p is the number of eigenvalues, identical to M, and N is the length of the original series. The separation occurs at the value of k ∈ {0, 1, 2, . . . , p − 1} for which the MDL is minimized.

The Fisher-Shannon Method
By the Fisher-Shannon method, the informational properties of a time series can be investigated, namely the Fisher Information Measure (FIM) and Shannon entropy (SE), which are used to quantify the local and global smoothness of the distribution of a series. The FIM and SE can be employed to characterize the complexity of non-stationary time series described in terms of order and organization [26]. The FIM measures the order and organization of the series, and the SE its uncertainty or disorder [27]. The FIM and SE are defined by the following formulae: where f (x) is the distribution of the series x. The Shannon entropy power N X is generally used instead of SE: to avoid dealing with negative quantities. FIM and N X are not independent of each other due to the isoperimetric inequality FIM·N X ≥ D [28], where D is the dimension of the space, which, for the time series, is 1. FIM and N X depend on f (x), whose accurate estimation is crucial in order to obtain reliable values of informational quantities. For calculating FIM and N X , we applied the kernel-based approach that Telesca and Lovallo [29] demonstrated to be better than the discrete-based approach. Thus, we apply the kernel density estimator method for f(x) [30,31] as shown in the following formula: where M and b denote the length of the series and the bandwidth, respectively, while K(u) is the kernel that is a continuous, symmetric, and non-negative function satisfying the two following constraints: f (x) is estimated by means of an optimized integrated procedure using the algorithms of Troudi et al. [32] and Raykar and Duraiswami [33], with a Gaussian kernel: Due to the isoperimetric inequality, the Fisher-Shannon information plane (FSIP), which has the N X as the x-axis and FIM as the y-axis, represents a very useful tool to investigate the time dynamics of signals [34]. For scalar signals, the curve FIM·N X = 1 separates the FSIP into two parts, and each signal can be represented by a point located only in the space FIM·N X > 1.

Results
We analyzed the 2000-2020 time variation of four vegetation indices: evapotranspiration (ET), normalized difference vegetation index (NDVI), leaf area index (LAI), and enhanced vegetation index (EVI).
First, the SSA was applied to each time series, and the value of the phase lag M was selected, taking into account the sampling time of the series (8 days for ET and LAI; 16 days for EVI and NDVI). To detect at least the annual cycle, M was set as 45 for the ET and LAI series, and as 24 for the EVI and NDVI series; moreover, these values fit well with Khan and Poskitt's [24] criterion, varying the length of the data from 503 values (EVI and NDVI) to 965 (ET) and 1003 (LAI). Figure 5 shows the application of the SSA to the ET time series of Appia as an example. Before applying the SSA, the original time series was normalized. Figure 5a shows the eigenvalue spectrum of the SSA decomposition; each eigenvalue corresponds to a reconstructed component and represents the fraction of the total variance of the original series explained by that component. Figure 5b shows all the obtained reconstructed components, whose behaviour varies from oscillatory with amplitude modulation to apparently noisy.
Applying the MDL criterion, the signal is separated into a trend and a detrended series; the value of k min corresponding to the minimum MDL represents the number of the first reconstructed components to sum up for obtaining the trend (Table 2). Applying this criterium to the ET time series of Appia, the MDL curve is shown in Figure 5c, and the minimum MDL is at k min = 11; thus, the trend is obtained by summing up the first 11 reconstructed components (Figure 5d) and the detrended series by subtracting the trend from the original normalized series (Figure 5e). Figures 6 and 7 show, similarly to Applying the MDL criterion, the signal is separated into a trend and a detrende series; the value of kmin corresponding to the minimum MDL represents the number of th first reconstructed components to sum up for obtaining the trend ( Table 2). Applying th criterium to the ET time series of Appia, the MDL curve is shown in Figure 5c, and th     The trend is characterized by an oscillatory behaviour that explains the seasonal cy- The trend is characterized by an oscillatory behaviour that explains the seasonal cycles of the series, very likely linked with the meteoclimatic variability. The detrended series, although apparently noisy, would represent the inner time dynamics of the series that might be not influenced by external driving mechanisms. Table 2 shows for all the investigated indices the value of the minimum of the MDL criterion.
Our aim is to characterize the time dynamics of the inner vegetation of the investigated sites by using the Fisher-Shannon method. Thus, for each site we focused on the detrended series, since this represents the inner time variability of vegetation not influenced by external meteoclimatic factors. Figure 8 shows the FSIP of ET (Figure 8a), EVI (Figure 8b), LAI (Figure 8c), and NDVI (Figure 8d). The FSIP indicates that the Castel Volturno site is characterized by the highest Shannon entropy power and the lowest FIM that suggests a low level of order and organization of vegetation indices; Appia Park, except for the ET, is characterized by the lowest Shannon entropy power and the lowest FIM that reveal a relative high level of order and organization of vegetation indices; Castel Porziano is generally characterized by an "intermediate" behaviour, since the vegetation indices, except ET, are located in the FSIP between Castel Volturno and Castel Porziano. series, although apparently noisy, would represent the inner time dynamics of the series that might be not influenced by external driving mechanisms. Table 2 shows for all the investigated indices the value of the minimum of the MDL criterion. Our aim is to characterize the time dynamics of the inner vegetation of the investigated sites by using the Fisher-Shannon method. Thus, for each site we focused on the detrended series, since this represents the inner time variability of vegetation not influenced by external meteoclimatic factors. Figure 8 shows the FSIP of ET (Figure 8a), EVI (Figure 8b), LAI (Figure 8c), and NDVI (Figure 8d). The FSIP indicates that the Castel Volturno site is characterized by the highest Shannon entropy power and the lowest FIM that suggests a low level of order and organization of vegetation indices; Appia Park, except for the ET, is characterized by the lowest Shannon entropy power and the lowest FIM that reveal a relative high level of order and organization of vegetation indices; Castel Porziano is generally characterized by an "intermediate" behaviour, since the vegetation indices, except ET, are located in the FSIP between Castel Volturno and Castel Porziano.

Discussion
The potential of satellite systems for the monitoring of vegetation r recognized and nowadays the most recent ICT technological developm of artificial intelligence and EO, and the growing availability of infor data from the free cloud such as GEE) have opened new frontiers and a The use of EO-based indicators for the monitoring of vegetation i vant, and recently the RS-derived data have been shown to be useful a

Discussion
The potential of satellite systems for the monitoring of vegetation resources is widely recognized and nowadays the most recent ICT technological developments, the joint use of artificial intelligence and EO, and the growing availability of information (as well as data from the free cloud such as GEE) have opened new frontiers and application fields.
The use of EO-based indicators for the monitoring of vegetation is particularly relevant, and recently the RS-derived data have been shown to be useful across many fields; however, undoubtedly, large earth observation data such as satellite time series pose several challenges to face in order to transform data in useful and reliable information.
For the purpose of our investigation, the four vegetation indices (ET, NDVI, LAI, and EVI) are various satellite products available in the cloud storage of Google Earth Engine, selected because they are related to the interactions between soil, water, atmosphere, and plants, while MODIS products were chosen due to their global coverage and long duration of data acquisition.
The analyses of satellite time series are generally quite complex and time-consuming due to the amount of data, but they are expected to be as suitable for the identification of both slow and fast changes as, for example, parasites or salinization, deforestation or wildfires, which adversely have been affecting NC during the last decades. In reality, the ability and effectiveness of change detection approaches and methods depend on the ability to account for the great variability exhibited by the seasonal variations (at seasonal and/or intra-annual scales) while identifying small multi-year trends and changes at diverse interannual time scales.
The methodological approach consisted of the following steps: decomposition of each satellite vegetation index through the SSA, detection of annual and seasonal cycles, separation between trend and detrended series, and application of the Fisher-Shannon method to the detrended series. In particular, this approach enabled us to perform the deseasonality and, therefore, to split the stronger seasonal dynamics from the subtle inner time variability of the investigated signals. Identifying and extracting information related to the potential presence of small but significant trends or variations in vegetation is an important issue, and the effectiveness of change detection approaches depends on their ability to account for both the great variability exhibited by the seasonal variations and the small multiyear changes that might be completely veiled by the seasonal dynamics.
Our study highlighted that the trend is characterized by an oscillatory behaviour that explains the seasonal cycles of the series, very likely linked with the meteoclimatic variability. Furthermore, the detrended series, whose variability our study has focused on, although apparently noisy, would represent the inner time dynamics of the series that might be not influenced by external driving mechanisms.
For Castel Volturno, a low level of order and organization of the MODIS time series was observed during the whole investigated period. This behaviour denoted an anomalous vegetational dynamic that can be explained and attributable to the effect of attack by the parasite Toumeyella parvicornis, which in the recent years adversely impacted the Pinus trees of the area, dramatically damaging them.
The reliability of the analytical results obtained from the Fisher-Shannon approach was assessed by comparisons with field surveys and independent data analyses. In fact, the results obtained from the statistical analysis herein conducted fit well with the results obtained from the processing of Sentinel 2 data jointly carried out by CNR and Carabinieri [35] and shown in Figure 9; the dark grey pixels (in Figure 9d) indicate the areas affected by a decreasing trend (site degradation) as a resulting effect of the parasite attacking the pinus trees; the white pixels are related to areas involved in increasing trend, mainly linked to agricultural activities. Finally, Figure 9f,g, acquired during the field survey, clearly provide evidence of the macroscopic effect of the Toumeyella parvicornis on the pinus trees. An example of this behaviour, i.e., grey and white pixels, related to decreasing and increasing trends, is shown in Figure 9e, where the blue and red lines depicted the maximum NDVI The inner dynamic of the vegetation of Appia Antica Park seems quite stable, and this was confirmed by in situ analysis. This site is mainly involved in and characterised by agricultural activities that were conducted systematically maintaining the same cultivation types for the whole period of our analysis [35]. The behaviour of the inner vegetation appears without anomalous dynamics, because the area was not involved in significant changes of vegetation status as well as of the land use and land cover as it can be seen from the Corine land cover updates (freely available online in the framework of the Copernicus initiative; see, for example, CORINE Land Cover-Copernicus Land Monitoring Service).
Castel Porziano, instead, presents FIM and Shannon Entropy values in the middle between those of Castel Volturno and Appia, except for the evapotranspiration. Comparison with independent data sets [35] can confirm that from 2000 to 2020 the area was quite The inner dynamic of the vegetation of Appia Antica Park seems quite stable, and this was confirmed by in situ analysis. This site is mainly involved in and characterised by agricultural activities that were conducted systematically maintaining the same cultivation types for the whole period of our analysis [35]. The behaviour of the inner vegetation appears without anomalous dynamics, because the area was not involved in significant changes of vegetation status as well as of the land use and land cover as it can be seen from the Corine land cover updates (freely available online in the framework of the Copernicus initiative; see, for example, CORINE Land Cover-Copernicus Land Monitoring Service).
Castel Porziano, instead, presents FIM and Shannon Entropy values in the middle between those of Castel Volturno and Appia, except for the evapotranspiration. Compari- son with independent data sets [35] can confirm that from 2000 to 2020 the area was quite stable, as it can be seen from Figure 10, where the Google Earth satellite pictures at higher resolution do not show particular changes in land cover. stable, as it can be seen from Figure 10, where the Google Earth satellite pictures at higher resolution do not show particular changes in land cover. Figure 10. Castel Porziano from GE high-resolution images, which show that from 2000 to 2020 no significant changes occurred (Google Earth Pro courtesy).

Conclusions
The vegetation of three study areas from the Central (Appia Ancient Park and Castel Porziano) to the Southern (Castel Volturno) parts of Italy were analyzed. The study areas were periurban and specifically selected as key in improving environmental quality: in fact, they are rich in biodiversity and allow urban areas to be more sustainable, helping to combat climate change and make cities more comfortable, as recently strongly emphasized by the COVID-19 pandemic emergency.
Thus, for each site we focused on the detrended series, since this represents the inner time variability of the vegetation not influenced by external meteoclimatic factors.
Results of our analyses highlighted that the (i) trend is characterized by an oscillatory behaviour that explains the seasonal cycles of the series, very likely linked with the meteoclimatic variability, (ii) detrended series, although apparently noisy, would represent the inner time dynamics of the series that might be not influenced by external driving mechanisms.
Among the sites investigated, Castel Volturno was characterized by the highest Shannon entropy power and the lowest FIM that indicate a low level of order and organization of the vegetation indices for this site. Independent analyses and field survey highlighted that Castel Volturno is strongly affected by a parasite, the Toumeyella parvicornis, which has been provoking dramatic damage to the Pinus trees in recent years.

Conclusions
The vegetation of three study areas from the Central (Appia Ancient Park and Castel Porziano) to the Southern (Castel Volturno) parts of Italy were analyzed. The study areas were periurban and specifically selected as key in improving environmental quality: in fact, they are rich in biodiversity and allow urban areas to be more sustainable, helping to combat climate change and make cities more comfortable, as recently strongly emphasized by the COVID-19 pandemic emergency.
Thus, for each site we focused on the detrended series, since this represents the inner time variability of the vegetation not influenced by external meteoclimatic factors.
Results of our analyses highlighted that the (i) trend is characterized by an oscillatory behaviour that explains the seasonal cycles of the series, very likely linked with the meteo-climatic variability, (ii) detrended series, although apparently noisy, would represent the inner time dynamics of the series that might be not influenced by external driving mechanisms.
Among the sites investigated, Castel Volturno was characterized by the highest Shannon entropy power and the lowest FIM that indicate a low level of order and organization of the vegetation indices for this site. Independent analyses and field survey highlighted that Castel Volturno is strongly affected by a parasite, the Toumeyella parvicornis, which has been provoking dramatic damage to the Pinus trees in recent years.
Our results could contribute to the definition of methods suitable for an early diagnosis of deterioration trends, and create operational tools for multiscale, multisensor, multitemporal monitoring of biophysical parameters relating to the state of vegetation.
Our future work will be in the application of robust statistical analyses to satellite time series to define, for example, indicators devised to assess and monitor land degradation that can be applied at different spatial and temporal scales using different satellite time series (MODIS along with Sentinel 2 data sets). In fact, the use of EO-based indicators is particularly relevant because RS-derived data have been shown to be useful across many fields, and at the local to global levels using data freely available from NASA (MODIS, TM) and ESA (Sentinel 1-5). Funding: This work was supported by CNR in the framework of the SMART Forest project, funded by Carabinieri CUFAA, and COELUM, funded by CNR.

Conflicts of Interest:
The authors declare no conflict of interest.