Increases in global mean temperature mainly caused by emissions of greenhouse gases are widely acknowledged by the scientific community. Further increases are expected in the future, with higher frequency and severity of extreme drought events and heat waves1 and related impacts and damages. Extreme event is generally defined as either taking maximum values or exceedance above pre-existing high thresholds,2 thus forecasting and evaluating the consequences of these events on forest ecosystems are emerging issues for both ecologists and forest managers.3 Drought events limit the development and growth rate of plants, their strength, and then biomass accumulation and carbon sink capacity at ecosystem level. Moreover, it induces tree mortality and forest degradation.4
In the recent decades, the frequency, duration, and severity of drought have increased substantially at a global scale,5 especially in arid and semiarid regions.1 Under limited availability of soil water, the first plant response is to close stomata,6 which decreases the inflow rate into the leaves.7 Thus, forest ecosystem responses to drought and their internal adjusting mechanisms vary according to drought intensity and duration, which add uncertainty of forest responses to climatic events.8 This assumes relevance also for Mediterranean forest species, commonly considered resistant to drought events.9
Extreme drought effects on forest ecosystems have been assessed at regional,10 national,11,12 and continental scales. Distinctively, Ciais et al.13 have demonstrated a reduction of about 30% in gross primary productivity, with consequences on net source of carbon dioxide to the atmosphere and reversed the effect of four years of net ecosystem carbon sequestration. Experimental studies carried out to assess the effects of such climatic events on tree growth, resilience, and mortality14–16 pointed out the need to investigate forest vulnerability and risk of dieback, especially in the Mediterranean basin, in the frame of global warming.
The assessment of forest health and tree growth is commonly carried through data collection by forest inventories and monitoring programs like the International Co-operative Program on assessment and monitoring of air pollution effects on forests.17 The main drawback of these methods, based on field plot sampling approach, is the low spatial information. On the other hand, high or very high-resolution images (e.g., WorldView-2, RapidEye satellite imagery) have been widely used for research issues, but these datasets are often expensive and not freely available. Medium-resolution imagery (particularly Landsat) was easier accessible and often cost-free for a broad majority of users, leading to many research projects. With the launch of the Sentinel-2 (S2) series, in June 2015, possibilities for research came into existence with great potential for forest monitoring because of its spectral, geometrical, and temporal resolutions. S2 combines a large swath, frequent revisit, and systematic acquisition of all land surfaces at high-spatial resolution and, with 13 spectral bands that guarantee consistent time series, showing variability in land surface conditions and minimizing any artifacts introduced by atmospheric variability. The increased swath width, along with the short revisit time, allows rapid changes to be monitored. In fact, S2 may systematically map different classes of cover, such as forest, crops, grassland, water surfaces, and artificial covers, like roads and buildings, without any need to plan a network of sample plots or point of interests investigated. Furthermore, S2 has the potential to highlight the health condition, growth, and productivity of terrestrial ecosystems at a wider scale.18
In this perspective, the objective of this study was to evaluate the capability of S2 data to detect the crown dieback in Mediterranean forest ecosystems in Tuscany (central Italy) after an extreme drought event occurred in summer 2017. Multitemporal (both intra-annual and interannual) S2 images and field data have been collected within three different forest stands, two deciduous broadleaved species, namely beech (Fagus sylvatica L.) and downy oak (Quercus pubescens Willd.), and an evergreen species, holm oak (Quercus ilex L.), in order to investigate on the response strategies of the species to the climatic event and test the capability for multitemporal monitoring climate-related processes of forest conditions at regional scale.
Material and Methods
Observational samples (118 in total) have been carried out within three different forest stands representative of the Tuscany region (Fig. 1): beech (Fagus sylvatica L., FS, covering about of the forest regional area), downy oak (Quercus pubescens Willd., QP, covering about ), and holm oak (Quercus ilex L., QI, covering about ), distributed along an ecological gradient, from Mediterranean (coastal) to continental (hilly and mountain zones), from west to east part of the region. For each forest type, stands with different soil bed rock and depth, altitude, slope, and aspect have been considered to represent the different conditions in which the species can grow. Information on topographic, structural, and community compositional data have also been collected within all the observational samples. Using the information on drought-induced leaves and branches desiccation and crown defoliation, each observational sample plot has been classified as “healthy” or “unhealthy.”
Considering the direct relation between air temperature with aspect and elevation, and between slope and drought, three additional topography layers (namely elevation, slope, and aspect) have also been considered. A digital terrain model at 20-m resolution was then used to create these layers.
S2 mission consists of two satellites, S2A launched in 2015 and S2B launched in 2017, which operate at an altitude of 705 km in the same orbit, phased at 180 deg to each other and with the orbit inclination of 98.5 deg. The satellites are equipped with modern multispectral high-resolution scanners, 13 spectral channels, and it has resolutions of 10, 20, and 60 m and swath width of 290 km. The revisit time is 10 days for one satellite and 5 days for the two satellites. In this study, 30 S2 tiles have been processed and used for the analysis (Table 1). According to Hadjimitsis et al.19, Sen2Cor processor has been used to derive surface reflectance data (bottom of atmosphere) from S2 level 1C data, provided in the top of atmosphere reflectance. The processing was performed by sen2R package.20
S2 tiles used in this work.
As a proxy of the photosynthetic vegetation activity, the three more common VIs that can be generated using the four channels with 10-m spatial resolution (centered at 490, 560, 665, and 842 nm wavelengths; blue, green, red, and NIR bands, respectively) have been calculated: normalized difference vegetation index (NDVI), red edge NDVI (RENDVI), and simple ratio index (SRI).
The analysis has been carried out focusing on the following different situations:
(i) each VI has been calculated for all summer months in 2017 (June, July, and August) to evaluate the drought effect on the forest tree health (crown) condition;
(ii) for each month, a comparison between observational samples classified as “healthy” and “unhealthy” has been carried out to highlight differences over calculated VIs;
(iii) FS and QI stands have been further analyzed by means of a multitemporal analysis considering the year before (2016) and after (2018) the drought event, to evaluate the presence of damages occurred before or after the extreme drought event.
The statistical significance of differences among VI values groups has been analyzed through -values obtained from -test.
The applied VIs have been able to analyze the dynamics of vegetation conditions in summer 2017 within the pixels referred to the observational samples (Fig. 2). Higher mean of VI values has been registered in June and lower in August, and the analysis of S2 data highlights a significant reduction () of VI values in August 2017 (for both NDVI, RENDVI, and SRI), mainly for deciduous forests (FS and QP). As expected, the unhealthy stands showed a reduction of the photosynthetic activity in August 2017 (Fig. 2). VI values were statistically similar () for both healthy and unhealthy stands in June and July, whereas they showed difference in August 2017 ().
In a further step, the photochemical properties of the studied stands in June, July, and August 2016 and, for the same period in 2018 (i.e., one year before and one after the extreme drought event considered, respectively), have been analyzed (Fig. 3). For FS forests, NDVI values registered for healthy and unhealthy stands in both July 2016, August 2016, and from April to August 2018 were similar (). For this species, similar mean NDVI values have been measured in healthy and unhealthy stands in July (for healthy stands: in 2016, in 2017, and in 2018; for unhealthy stands: in 2016, in 2017, and in 2018) and August (for healthy stands: in 2016, in 2017, and in 2018; for unhealthy stands: in 2016, in 2017, and in 2018).
Conversely, QI forests have shown similar values of mean NDVI in healthy and unhealthy stands till August 2016. Despite QP stands are dominated by broadleaved species, such as FS stands, their response to the drought event was more comparable to the one had by QI (evergreen) stands.
Discussion and Conclusions
This case study addressed the sensitivity of S2 imagery to detect extreme drought effects on tree crown health status in central Italy, comparing the response of three different forest types, characterized by two broadleaved and one evergreen tree species, typical of Mediterranean forest environments. Detecting of such effects has been based on seasonal variations of three VIs (namely NDVI, RENDVI, and SRI) in three years: 2016, 2017, and 2018. The three VIs have similar performance (Fig. 3) and can be alternatively used for operative purposes, e.g., tools for automatic detection from S2 imageries. These findings, obtained exploiting yearly VIs trajectories, have been validated by visual inspection of all 118 observational samples performed after drought event in 2018 and are consistent with phenological behavior of the examined species. For example, beech is considered a relatively drought-susceptible species, able to optimize its photosynthetic assimilation rates under mild drought conditions. This behavior became disadvantageous under severe drought events, as the one occurred in 2017, leading to desiccation damages in leaves and branches. Leaves of this species are stress-sensitive particularly in spring, during leaf emergence.21,22 However, Polle et al.23 found that mature leaves in late summer will become stress-sensitive if an extended period of drought and elevated air temperatures occurs. The anticipated reduction of tree vitality in beech mature leaves was detected by S2 images as a strong reduction of VI values in August 2017. Despite this, as pointed out by recent studies,24 an important aspect of drought tolerance in trees is their ability to recover after drought release and quickly resume full physiological activity. The analysis of S2 data confirmed, without field surveys, that beech can recover the photosynthetic activity due to a leaf mass, developed after the drought events (Fig. 3).
In parallel, as typical Mediterranean sclerophyllous species, holm oak is tolerant to short and occasional extreme events, as confirmed by VI trajectories (Fig. 3, bottom part). These recursive stresses generate morphological, genetic, and biochemical modifications at tree level that can be positive for trees, making them more resistant to future exposure to the same stress factor.25 However, the increased frequency, length, and intensity of drought events experienced in the Mediterranean basin in last two decades, with high rate of evapotranspiration and low precipitation, negatively affect the vegetation growth of holm oaks.26
Compared with other similar satellites (e.g., mainly Landsat) and considering its recent launch, the main limitation of S2 on environmental monitoring relies on historical time series analysis.27 On the other hand, the finer geometric, temporal, and spectral resolutions make S2 multitemporal data suitable to monitor drought-induced canopy damages in forest ecosystems has been demonstrated, confirming the potential of S2 satellites for applications in forestry and vegetation analysis.18
The case study presented shows that S2 can be applied for multitemporal monitoring climate-related processes providing a synthetic overview of forest conditions at regional scale, and, by means of VIs, to point out different response strategies of tree species to extreme climatic events. S2-10-m spatial resolution seems to be adequate to investigate punctual phenomena, like drought-induced damages on forest stands, providing information about potential declines in vegetation health very quickly (early warning) without fieldwork.
If properly established in a web GIS application (e.g., Google Earth Engine), this kind of analysis can be readily used for forest management to detect and map the occurrence of drought-induced damages at forest stand level over wide areas.
We would like to thank Leonardo Tonveronachi (CREA, Research Centre for Forestry and Wood) and Giovanni Iacopetti (University of Florence) for their support in field data collection.