Analysis of Olive Grove Destruction by Xylella fastidiosa Bacterium on the Land Surface Temperature in Salento Detected Using Satellite Images

Agricultural activity replaces natural vegetation with cultivated land and it is a major cause of local and global climate change. Highly specialized agricultural production leads to extensive monoculture farming with a low biodiversity that may cause low landscape resilience. This is the case on the Salento peninsula, in the Apulia Region of Italy, where the Xylella fastidiosa bacterium has caused the mass destruction of olive trees, many of them in monumental groves. The historical land cover that characterized the landscape is currently in a transition phase and can strongly affect climate conditions. This study aims to analyze how the destruction of olive groves by X. fastidiosa affects local climate change. Land surface temperature (LST) data detected by Landsat 8 and MODIS satellites are used as a proxies for microclimate mitigation ecosystem services linked to the evolution of the land cover. Moreover, recurrence quantification analysis was applied to the study of LST evolution. The results showed that olive groves are the least capable forest type for mitigating LST, but they are more capable than farmland, above all in the summer when the air temperature is the highest. The differences in the average LST from 2014 to 2020 between olive groves and farmland ranges from 2.8 ◦C to 0.8 ◦C. Furthermore, the recurrence analysis showed that X. fastidiosa was rapidly changing the LST of the olive groves into values to those of farmland, with a difference in LST reduced to less than a third from the time when the bacterium was identified in Apulia six years ago. The change generated by X. fastidiosa started in 2009 and showed more or less constant behavior after 2010 without substantial variation; therefore, this can serve as the index of a static situation, which can indicate non-recovery or non-transformation of the dying olive groves. Failure to restore the initial environmental conditions can be connected with the slow progress of the uprooting and replacing infected plants, probably due to attempts to save the historic aspect of the landscape by looking for solutions that avoid uprooting the diseased plants. This suggests that social-ecological systems have to be more responsive to phytosanitary epidemics and adapt to ecological processes, which cannot always be easily controlled, to produce more resilient landscapes and avoid unwanted


Introduction
Mediterranean landscapes are strongly influenced by human activities that shape the structure and functions of the landscape in consideration of human needs. Currently, the olive grove if it is not appropriately managed or replaced, would over time be characterized as that of herbaceous vegetation. Therefore, a Recurrence Quantification Analysis (RQA) was performed to detect LST evolution using an ex-ante and ex-post approach to the spread of X. fastidiosa by applying a window of 18 years. The RQA is non-linear and useful for highlighting and detecting vegetation indices that show the resilience of the system at to disasters [14]. The recurrences analysis can help interpret land surface variation on a provincial scale to show the cumulative effect of X. fastidiosa on the olive groves, highlight its potential start, and show the human response to the epidemic.

Study Area
An LST analysis of the agro-ecosystem that characterizes the Salento Peninsula in the Apulia region of southern Italy was performed [5] (Figure 1). The agro-ecosystem is characterized mainly by olive groves, farmland and vineyards. Natural vegetation is characterized by isolated forested nuclei, grass and small wetlands [26]. At present, the whole province of Lecce is declared to be an X. fastidiosa infected area (Commission Implementing Regulation (EU) 2020/1201). The agro-ecosystem is characterized mainly by olive groves, farmland and vineyards. Natural vegetation is characterized by isolated forested nuclei, grass and small wetlands [26]. At present, the whole province of Lecce is declared to be an X. fastidiosa infected area (Commission Implementing Regulation (EU) 2020/1201).

Satellite Data
The study was carried out using Landsat 8 data and Moderate Resolution Imaging Spectroradiometer (MODIS) images provided by the National Aeronautics and Space Administration (NASA) and the U.S. Geological Survey (USGS) [27,28]. Data were analysed using the QGIS, Octave and Matlab software. Band 4, Band 5 (with spatial resolution of 30 m) and Band 10 (spatial resolution of 100 m) from Landsat 8 were used to analyze LST variation in the spring and summer from 2014 to 2020 [29,30].
The MOD11A2 product was used to build the LST temporal profile from 2003 to 2020 with a 1 km spatial resolution and an 8-day time step [31].
Land cover from the Apulia Territorial Information System (SIT Puglia) for 2011 was used to build a Region of Interests (RoI) and extrapolate LST data for olive groves, farmland, vineyards and forest [32]. Furthermore, to evaluate the effect of climate conditions on the satellite data, we used air temperature and precipitation data provided by Civil Protection of the Apulia region [33].

Analysis of LST Landsat 8 Data for Land-Cover Types
The LST for olive groves, farmland, vineyards and forest was calculated for the spring and summer from 2014 to 2020, one year after official pathogen detection (2013). Landsat 8 images from 1 March to 10 May (for the springtime) and from 21 July to 11 August (for the summertime) were used. For each time window, only images with high quality and low cloud cover were selected.
For this analysis, Band 10 (thermal infrared) was converted in Brightness Temperature, while band 4 (red) and band 5 (infrared) were useful for determining the Normalized Difference Vegetation Index (NDVI) [29,30,34,35] as follows: The index was used to estimate the land cover emissivity (Table 1). For an NDVI in the range of 0.2-0.5 the emissivity was calculated using the following equation [29,30].
where ε v and ε s are the land cover emissivity and C is the surface roughness (considered with an uniform value of 0.005 [30,36]. The LST was calculated using the following equation [29,30,37]: where BT is the brightness temperature.

Recurrence Analysis of LST Time-Series Using MODIS Data
To analyze the climate effect of X. fastidiosa and the potential evolution of regional plans and policies, LST time series from 2013 to 2020 were built for farmland and olive groves using the MOD11A2. Here, we studied their relationship by comparing their their recurrence structure. Each land cover is characterized by the periodic recurrence of vegetative states over a season influenced by climate conditions, such as periodic recurrences caused by the seasonal effect in the vegetation time series [20]. This recurrence state can change in time as a consequence of disturbing events and can be investigated by the recurrence analysis [38].
A recurrence structure change can be highlighted using recurrence plots (RPs), which indicate all pairs of time points (i, j) where the system recurred to a preceding state x i [39,40]: The recurrence matrix R consists of a characteristic pattern and can be considered as a fingerprint of the dynamics represented by the time series. The Cross Recurrence Plot (CRP) highlights the co-occurrence of similar states of two time series, whereas the Joint Recurrence Plot (JRP) highlights the simultaneous occurrence of recurrences between two time series [14,39,40]. The CRP between the variables x i and y i is determined by: and the JRP is determined by the RPs multiplication of each system x i and y i : where x i , y i ∈ R represents the farmland and olive grove LST time series, respectively; I, j = 1, . . . .N; N is x i , y i ; time series length ε is a threshold distance; θ(•) is the Heaviside function and • is a norm [14,39].
In the CRP, the main measure of similarity is represented by the Line of Synchronization (LoS). For two equal signals, the LoS corresponds to a continuous line in the main diagonal in the plot, whereas deviations between two signals correspond to interruptions, or bendings, of the LoS (Figure 2A-C). Instead, the LoS completely vanishes for completely different signals. The LOS is useful for studying synchronization or time distortions [41]. Each land cover is characterized by the periodic recurrence of vegetative states over a season influenced by climate conditions, such as periodic recurrences caused by the seasonal effect in the vegetation time series [20]. This recurrence state can change in time as a consequence of disturbing events and can be investigated by the recurrence analysis [38].
A recurrence structure change can be highlighted using recurrence plots (RPs), which indicate all pairs of time points (i, j) where the system recurred to a preceding state xi [39,40]: The recurrence matrix R consists of a characteristic pattern and can be considered as a fingerprint of the dynamics represented by the time series. The Cross Recurrence Plot (CRP) highlights the co-occurrence of similar states of two time series, whereas the Joint Recurrence Plot (JRP) highlights the simultaneous occurrence of recurrences between two time series [14,39,40]. The CRP between the variables and is determined by: and the JRP is determined by the RPs multiplication of each system and : where x i , y i ∈ R represents the farmland and olive grove LST time series, respectively; I, j =1,….N; N is , ; time series length ε is a threshold distance; θ(•) is the Heaviside function and ‖•‖is a norm [14,39].
In the CRP, the main measure of similarity is represented by the Line of Synchronization (LoS). For two equal signals, the LoS corresponds to a continuous line in the main diagonal in the plot, whereas deviations between two signals correspond to interruptions, or bendings, of the LoS (Figure 2A-C). Instead, the LoS completely vanishes for completely different signals. The LOS is useful for studying synchronization or time distortions [41].  Therefore, for the synchronization analysis, the F-recurrence measure defined for the CRP was used. It is given by 1 N ∑ i C i,i , which quantifies the density of points along the main diagonal. The F-recurrence measure is 1 for total synchronization and tends to very small numbers for highly desynchronized systems [14].
The JRP is always a main diagonal because it is based on individual RP signals. Here, the analysis is carried out by calculating the recurrence points loss in the JRP with respect to the individual RPs. In two signals with the same recurrence pattern, the JRP corresponds to individual RPs because each same recurrence point between two RPs is preserved in the JRP, whereas specific deviations between the two signals in individual RPs can produce a corresponding loss of recurrence point in the JRP. This allows checking for generalized synchronization (Figure 2A,B,D) [14].
The diagonal lines in JRP indicate that the two variables have states that repeat at the same time (each variable recurred at the same time) and therefore the variable presents a similar evolution. Determinism Det = is a measure that quantifies the fraction of recurrence points that form the diagonal lines in the JRP with P(l) the histogram of lengths of such diagonal lines [40][41][42].
The CRP and JRP analysis was combined with a sliding window approach to detect the temporal effect of X. fastidiosa in LST time series evolution and synchronization. For CRP, overlapping was applied, and a sliding temporal windows with a "one-year length" to the LST time series was used to calculate the Det and F-recurrence parameters ex-ante and ex-post the X. fastidiosa spread in consideration of the assumption that it may have been introduced earlier than 2013 [14,16]. For the recurrence threshold, an empirically defined threshold was used to ensure a fixed number of recurrence points (measured by recurrence rate, The selected value did not qualitatively affect the used recurrence measures because of its low qualitative effect in this and similar studies [14,40,41]. The RQA was performed with "CRP Toolbox" in MATLAB software (http://tocsy.pik-potsdam.de/CRPtoolbox/ accessed on 15 June 2021).

Analysis of Climate Data
Since the LST can be affected by meteorological conditions, data for each Landsat image came from 8 meteorological stations managed by the Apulian Civil Protection [33] ( Figure 3). Specifically, the air temperature for six days before image acquisition was considered and the average value was calculated. Moreover, for the same period, the cumulate precipitation was also acquired.

LST Landsat 8 Analysis
The analysis of average LST values in summer for different land cover types showed that forests had a high capacity for reducing LST while farmland had the lowest ( Figure  4). The olive groves and vineyards showed average values of the LST; from 2014 to 2016, the LST of the olive groves was lower than vineyards, while from 2017 to 2020, it was higher.

LST Landsat 8 Analysis
The analysis of average LST values in summer for different land cover types showed that forests had a high capacity for reducing LST while farmland had the lowest ( Figure 4). The olive groves and vineyards showed average values of the LST; from 2014 to 2016, the LST of the olive groves was lower than vineyards, while from 2017 to 2020, it was higher. As reported, the LST value is influenced by climate conditions: the low value of LST for 2014 and 2018 corresponds to a high value of cumulative precipitation and low air temperature, while the highest value for the year 2017 corresponds to the absence of pre-  As reported, the LST value is influenced by climate conditions: the low value of LST for 2014 and 2018 corresponds to a high value of cumulative precipitation and low air temperature, while the highest value for the year 2017 corresponds to the absence of precipitation and the highest air temperature ( Figure 5). The Fisher test showed that the average LST for different land cover types for different years belonged to distributions having a different variance (Table 2), while the average values of the Student's t test were statistically different. Indeed, in the tests the values were higher than the references tabulated to determine the significance. The analysis of average values of LST in the spring for different land cover types shows that the forest land cover always had the highest ability to reduce LST, while the vineyards had the lowest ( Figure 6). The analysis of average values of LST in the spring for different land cover types shows that the forest land cover always had the highest ability to reduce LST, while the vineyards had the lowest ( Figure 6).
The olive groves and farmland had close average LST values. In general, the differences were less pronounced than in the summer. The Fisher test showed that the average LST for different land cover types for the different years belong to distributions with different variance (Table 3), while the average values for the Student's t test were statistically different. Indeed, in the tests the values were higher than the references tabulated to determine significance. Table 2. Fisher test and Student test to analyse the statistical differences among the land covers reported in the Figure 4.  The olive groves and farmland had close average LST values. In general, the differences were less pronounced than in the summer. The Fisher test showed that the average LST for different land cover types for the different years belong to distributions with different variance (Table 3), while the average values for the Student's t test were statistically different. Indeed, in the tests the values were higher than the references tabulated to determine significance. Table 3. Fisher test and Student test to analyse the statistical differences among the land covers reported in Figure 6. In this case, the LST values seem to be mainly influenced by air temperature as high  Table 3. Fisher test and Student test to analyse the statistical differences among the land covers reported in Figure 6. In this case, the LST values seem to be mainly influenced by air temperature as high values of air temperature correspond to high values of LST and vice-versa (Figure 7).

Recurrence Quantification Analysis of LST Time Series
To better analyze the impact of X. fastidiosa on the LST evolution of olive groves, the JRP between olive grove and farmland LSTs was carried out from 2003 to 2020, using a four-year sliding temporal windows. For each time window, the Det parameter was estimated.
The JRP-based measure identified a shift in the generalized synchronization around 2009-2010. This shift could be interpreted as the evolution of olive grove land surface temperature to one of farmland ( Figure 8).

Recurrence Quantification Analysis of LST Time Series
To better analyze the impact of X. fastidiosa on the LST evolution of olive groves, the JRP between olive grove and farmland LSTs was carried out from 2003 to 2020, using a fouryear sliding temporal windows. For each time window, the Det parameter was estimated.
The JRP-based measure identified a shift in the generalized synchronization around 2009-2010. This shift could be interpreted as the evolution of olive grove land surface temperature to one of farmland ( Figure 8).
The synchronization between farmland and olive grove LSTs was evaluated using a one-year sliding temporal windows on CRP before and after the spread of X. fastidiosa. The F-recurrence was calculated in each time step for the two LSTs.
The analysis revealed a pretty small value for the F-recurrences, i.e., the farmland and olive groves LSTs were only weakly related. Moreover, this weak interrelationship did not change over time, as shown by very similar values for the F-recurrences. According to the F-test, there is no statistically significant difference in the distribution of the F-recurrence between the ex-ante and ex-post 2010 time. The p-values of the F-test was 5% higher than the threshold (Table 4), so we could not reject the null hypothesis that the mean values for the F-recurrence distribution before and after the trigger event were referenced at distribution values with a similar variance. This suggests that X. fastidiosa did not have effect on the pattern synchronization of the variables. The synchronization between farmland and olive grove LSTs was evaluated using a one-year sliding temporal windows on CRP before and after the spread of X. fastidiosa. The F-recurrence was calculated in each time step for the two LSTs.
The analysis revealed a pretty small value for the F-recurrences, i.e., the farmland and olive groves LSTs were only weakly related. Moreover, this weak interrelationship did not change over time, as shown by very similar values for the F-recurrences. According to the F-test, there is no statistically significant difference in the distribution of the Frecurrence between the ex-ante and ex-post 2010 time. The p-values of the F-test was 5% higher than the threshold (Table 4), so we could not reject the null hypothesis that the mean values for the F-recurrence distribution before and after the trigger event were referenced at distribution values with a similar variance. This suggests that X. fastidiosa did not have effect on the pattern synchronization of the variables.

Discussion
The LST Landsat images showed that the olive groves did not have the same ability as the forest to mitigate the rise in temperature; indeed, in both times the olive grove LST was higher than the forest LST. This difference was accentuated in the summer and could have been influenced by the differences in their respective canopy structures. Therefore, even if many groves were characterized by monumental plants with trunk diameters greater than 80 cm, the canopy structure was less dense than that of forests, which could explain the low ability to regulate the microclimate of the ecosystem (Figure 9).

Discussion
The LST Landsat images showed that the olive groves did not have the same ability as the forest to mitigate the rise in temperature; indeed, in both times the olive grove LST was higher than the forest LST. This difference was accentuated in the summer and could have been influenced by the differences in their respective canopy structures. Therefore, even if many groves were characterized by monumental plants with trunk diameters greater than 80 cm, the canopy structure was less dense than that of forests, which could explain the low ability to regulate the microclimate of the ecosystem (Figure 9). Furthermore, the LST of the olive groves was influenced by background vegetation, which are grassland or bare soil considering the agriculture practices in spring and summer ( Figure 10). Indeed, we observed some similarity between olive groves and farmland in the spring and with vineyards in the summer. This was confirmed by the F-recurrence parameter between olive grove and farmland LSTs before and after 2010 that showed that these land cover types had the same synchronization independent of the spread of X. fastidiosa. Furthermore, the LST of the olive groves was influenced by background vegetation, which are grassland or bare soil considering the agriculture practices in spring and summer ( Figure 10). Indeed, we observed some similarity between olive groves and farmland in the spring and with vineyards in the summer. This was confirmed by the F-recurrence parameter between olive grove and farmland LSTs before and after 2010 that showed that these land cover types had the same synchronization independent of the spread of X. fastidiosa. Furthermore, the LST of the olive groves was influenced by background vegetation, which are grassland or bare soil considering the agriculture practices in spring and summer ( Figure 10). Indeed, we observed some similarity between olive groves and farmland in the spring and with vineyards in the summer. This was confirmed by the F-recurrence parameter between olive grove and farmland LSTs before and after 2010 that showed that these land cover types had the same synchronization independent of the spread of X. fastidiosa. The Landsat analysis did not show any direct effect of X. fastidiosa on LST even though a drastic increase in olive grove LST from 2014 to 2020 was expected as a consequence of the death of olive trees mainly of the monumental trees. This is because the bacterium diffusion presumably started some years before the 2013 detection but achieved maximum diffusion and effect later, forcing the Apulia region to declare the study area compromised. The meteorological conditions during image acquisition had a strong effect on LSTs. In any case, this was also influenced by different factors. First, infected and dead olive trees remained in the field for a long time because cutting operations only began in recent years. Second, dying plants still developed basal suckers with leaves, which shaded the soil affecting LST values ( Figure 11). The Landsat analysis did not show any direct effect of X. fastidiosa on LST even though a drastic increase in olive grove LST from 2014 to 2020 was expected as a consequence of the death of olive trees mainly of the monumental trees. This is because the bacterium diffusion presumably started some years before the 2013 detection but achieved maximum diffusion and effect later, forcing the Apulia region to declare the study area compromised. The meteorological conditions during image acquisition had a strong effect on LSTs. In any case, this was also influenced by different factors. First, infected and dead olive trees remained in the field for a long time because cutting operations only began in recent years. Second, dying plants still developed basal suckers with leaves, which shaded the soil affecting LST values ( Figure 11).
However, in the summer of 2019 and 2020, the LST of the olive groves was near the value recorded for farmland, whereas it was lower than for vineyards from 2014 to 2016, and the highest values from 2017 to 2020. This could be connected to the low canopy structures of the olive groves which had a major effect in the evolution to a farmland LST.
The evolution of the olive grove LST into a farmland LST because of X. fastidiosa diffusion was confirmed by the Det profile of JRP, which showed an increase after 2010. This seemed to support the argument that X. fastidiosa was introduced before 2013 [16]. The substantial variation of Det in behavior before 2010 can be linked to cultivation practices like olive tree pruning. On the other hand, the Det analysis after 2010 showed a more or less constant behavior without substantial variation; therefore, this can be the index of a static situation in the evolution of this system. It does not indicate a recovery or transformation for the dying olive groves.
The olive groves can have effects in provide microclimate mitigation ecosystem services in the agroecosystem landscape, which are important in the summertime where the air temperatures have the highest values. For example, Table 5 shows how the average LST from 2014 to 2020 between olive groves and farmland ranged from 2.8 • C to 0.8 • C. These differences can affect the rural landscape where there is a strong interaction between urban sprawl and olive groves ( Figure 12).
This study was carried out at a provincial scale; therefore, it can describe general situations but specific analysis at a smaller scale is needed for a specific context, such as the interface between urban areas and olive groves or where the isolated urban settlements extend into olive groves [43].
Indeed, a limitation of the study was the spatial resolution of satellite images considering the fragmentation of the landscape. Therefore, it was not possible to apply the Recurrence Quantification Analysis to forests and vineyards due to the small size of the patch in these land cover types and the poor spatial resolution (1 km) of the MODIS data for LST. The results reported here are in accordance with the literature, indicating that RQA analysis applied at remote sensing vegetation indices and other parameters is an effi-cient approach for analyzing the connection between climate systems and environmental variables over extensive areas [14,[44][45][46][47]. Moreover, these findings show the usefulness of analyzing the response of a landscape system to a specific disaster or to highlight the evolution of a system into a new configuration [14,[44][45][46][47]. However, in the summer of 2019 and 2020, the LST of the olive groves was near the value recorded for farmland, whereas it was lower than for vineyards from 2014 to 2016, and the highest values from 2017 to 2020. This could be connected to the low canopy structures of the olive groves which had a major effect in the evolution to a farmland LST.
The evolution of the olive grove LST into a farmland LST because of X. fastidiosa diffusion was confirmed by the Det profile of JRP, which showed an increase after 2010. This seemed to support the argument that X. fastidiosa was introduced before 2013 [16]. The substantial variation of Det in behavior before 2010 can be linked to cultivation practices like olive tree pruning. On the other hand, the Det analysis after 2010 showed a more or less constant behavior without substantial variation; therefore, this can be the index of a static situation in the evolution of this system. It does not indicate a recovery or transformation for the dying olive groves.
The olive groves can have effects in provide microclimate mitigation ecosystem services in the agroecosystem landscape, which are important in the summertime where the air temperatures have the highest values. For example, Table 5 shows how the average LST from 2014 to 2020 between olive groves and farmland ranged from 2.8 °C to 0.8 °C. These differences can affect the rural landscape where there is a strong interaction between urban sprawl and olive groves ( Figure 12). This study was carried out at a provincial scale; therefore, it can describe general situations but specific analysis at a smaller scale is needed for a specific context, such as the interface between urban areas and olive groves or where the isolated urban settlements extend into olive groves [43].
Indeed, a limitation of the study was the spatial resolution of satellite images considering the fragmentation of the landscape. Therefore, it was not possible to apply the Recurrence Quantification Analysis to forests and vineyards due to the small size of the patch in these land cover types and the poor spatial resolution (1 km) of the MODIS data for LST. The results reported here are in accordance with the literature, indicating that RQA analysis applied at remote sensing vegetation indices and other parameters is an efficient approach for analyzing the connection between climate systems and environmental variables over extensive areas [14,[44][45][46][47]. Moreover, these findings show the usefulness of analyzing the response of a landscape system to a specific disaster or to highlight the evolution of a system into a new configuration [14,[44][45][46][47].

The Vision of Recurrence Analysis in Panarchy Approach
X. fastidiosa produced a fast change in the landscape, destroying a historical agroecosystem that had been developed over centuries. Moreover, the analysis showed how ecological processes act independently and unpredictably. At the same time, the social The Vision of Recurrence Analysis in Panarchy Approach X. fastidiosa produced a fast change in the landscape, destroying a historical agroecosystem that had been developed over centuries. Moreover, the analysis showed how ecological processes act independently and unpredictably. At the same time, the social and economic components of the social-ecological system showed low resilience and a slow adaptative ability to respond adequately to a disaster.
The JRP analysis showed a static evolution of affected olive groves areas and if the institutional components had developed policies to replant olive trees or reproduce a similar canopy structure using different woody plants. Therefore, this can be considered the example of a panarchy effect where the staticity can be corresponding at the slow capacity of social and institutional components to develop adaptative strategies at such fast perturbations actions generated at the landscape scale. For this reason, decisionmakers must act with greater incisiveness and take more effective actions to tackle the environmental problem created by X. fastidiosa ( Figure 13).
Resistance to uprooting the infected monumental trees at the outset of the X. fastidiosa diffusion was linked to historical landscape conservation and may have slowed down landscape regeneration. We cannot say whether a more timely uprooting of the infected plants (as well as insect vector control) could have reduced the spread of the disease, but replacing them with new plants (resistant to the pathogen or non-host plants) would have helped develop new olive vegetation by accelerating the resilience mitigation of the local microclimate and others not studied here, but which are directly connected with the olive culture and the production of extra virgin olive oil.
institutional components had developed policies to replant olive trees or reproduce a similar canopy structure using different woody plants. Therefore, this can be considered the example of a panarchy effect where the staticity can be corresponding at the slow capacity of social and institutional components to develop adaptative strategies at such fast perturbations actions generated at the landscape scale. For this reason, decision-makers must act with greater incisiveness and take more effective actions to tackle the environmental problem created by X. fastidiosa ( Figure 13). Figure 13. Interpretation of the Det evolution of recurrence analysis considering regional and local processes.
Resistance to uprooting the infected monumental trees at the outset of the X. fastidiosa diffusion was linked to historical landscape conservation and may have slowed down Figure 13. Interpretation of the Det evolution of recurrence analysis considering regional and local processes.

Conclusions
In this study, remote sensing data was useful for analyze the Xylella fastidiosa effect in a retrospective way and in the province of Lecce (South Italy), where in situ field variable measurement is not available. This type of analysis is useful for predicting the landscape evolution and ecological processes on a local scale. Furthermore, the study highlighted how nature responded faster than humans to the bacterial invasion. For the aim of the study, the JRP was more relevant than CRP for analyzing the differences caused by X. fastidiosa, possibly because the CRP was influenced by the breakpoint choice. Therefore, the JRP showed a strong ability to study the evolution of the system in a non-linear way to monitor the X. fastidiosa pathogen in a way that was less influenced by the setting of the analysis or without knowing the precise time when it occurred.
The destruction of the monumental olive groves in the province of Lecce was not only an agricultural and economic disaster that had a negative effect on the aesthetics of the landscape, but it also affected Land Surface Temperature (LST) mitigation. Regional landscape regenerations policies should promote new land cover patterns and functions with a high capacity to mitigate microclimates that have an economic and cultural aspect. Furthermore, the recurrence analysis showed a change in the albedo; therefore, this study could be used to focus on the relationship between land cover change caused by X. fastidiosa and biogeochemical processes to develop regional land rehabilitation policies [48].
The easiest landscape regeneration would be to replace the dead olive trees with new ones with a high resistance to X. fastidiosa, as promoted by the regional policy. However, the result may not be as expected in a social-ecological system characterized by non-linear and non-predictable behavior, because social and economic considerations can have a strong impact on land cover pattern evolution. The tree deaths could induce growers to change the use of the land, abandon it or lose interest in olive cultivation, as evidenced by the fact that many lands in the area are for sale [5]. Indeed, the elimination of the dead trees would present an economic cost that some landowners could be not able to sustain; therefore, land evolution is uncertain. Moreover, to return to optimal land surface temperatures, the choice of a new configuration and plant structure for new olive groves would also be important.
It will be interesting to analyze whether the effect of LST variation has any effect on local biogeophysical processes and climate and how the LST variations recorded would be amplified by a global temperature increase of 1-2 degrees if adequate strategies to restore the tree cover are not adopted.