Next Article in Journal
Remotely Sensed Changes in Vegetation Cover Distribution and Groundwater along the Lower Gila River
Previous Article in Journal
Spatial Transformation of a New City in 2006–2020: Perspectives on the Spatial Dynamics, Environmental Quality Degradation, and Socio—Economic Sustainability of Local Communities in Makassar City, Indonesia
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Land
Volume 9
Issue 9
10.3390/land9090325
land-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Use of Sentinel-2 Data to Monitor Heterogeneous Effects of Contextual Drought and Heatwaves on Mediterranean Forests

by
Rosa Coluzzi
1,*,
Simonetta Fascetti
2,
Vito Imbrenda
1,
Santain Settimio Pino Italiano
2,
Francesco Ripullone
2 and
Maria Lanfredi
1
1
IMAA—CNR (Institute of Methodologies for Environmental Analysis—Italian National Research Council), c.da Santa Loja snc, 85050 Tito Scalo (PZ), Italy
2
School of Agricultural, Forest, Food and Environmental Sciences, University of Basilicata, Viale dell’Ateneo Lucano 10, I-85100 Potenza, Italy
*
Author to whom correspondence should be addressed.
Land 2020, 9(9), 325; https://doi.org/10.3390/land9090325
Submission received: 21 July 2020 / Revised: 1 September 2020 / Accepted: 11 September 2020 / Published: 14 September 2020
Browse Figures
Versions Notes

Abstract

:
The use of satellite data to detect forest areas impacted by extreme events, such as droughts, heatwaves, or fires is largely documented, however, the use of these data to identify the heterogeneity of the forests’ response to determine fine scale spatially irregular damage is less explored. This paper evaluates the health status of forests in southern Italy affected by adverse climate conditions during the hot and dry summer of 2017, using Sentinel-2 images (10m) and in situ data. Our analysis shows that the post-event—NDVI (Normalized Difference Vegetation Index) decrease, observed in five experimental sites, well accounts for the heterogeneity of the local response to the climate event evaluated in situ through the Mannerucci and the Raunkiaer methods. As a result, Sentinel-2 data can be effectively integrated with biological information from field surveys to introduce continuity in the estimation of climate change impacts even in very heterogeneous areas whose details could not be captured by lower resolution observations. This integration appears to be a successful strategy in the study of the relationships between the climate and forests from a dynamical perspective.

1. Introduction

The health of forests has been under increased scrutiny over the last few decades and climate change is a major contributor [1,2]. In particular, the stochasticity, frequency and magnitude of extreme weather events outside the range of normal values or rare within a statistical range at a particular place and time of year [3] have increased on a global scale [4,5,6], and their damaging impacts are expected to increase in the future [7]. These negative trends are mainly increasing in Europe [8], especially in the Mediterranean regions, where the climate has been warmer and drier [9,10] since the end of the last century. Mediterranean areas, commonly considered as rather resistant to drought event [11], have been recently characterized by more frequent and severe drought episodes which have already significantly impacted forests [12,13], triggering forest decay and mortality episodes [14,15,16,17]. For this reason, the beginning of the new millennium has been identified as a general period of Mediterranean forest decline related to drought [18]. Especially when drought is contextual to a prolonged period of abnormally hot weather (heatwave), the effects on forests can be destructive.
Under limited availability of soil water, the first plant response is to close stomata, which decreases the CO2 inflow into the leaves and determines the consequent reduction of structural carbohydrates [19]. Thus, forest ecosystem responses to drought and their adaptation depend on the drought intensity and duration [20,21]. Moreover, some tree species are more sensitive to drought-triggered growth decline and damage [21,22,23,24] and sometimes the response of forests to drought can be exacerbated by local environmental and site conditions [24,25,26,27,28]. Although, several studies have highlighted a diverse vulnerability of vegetation to drought damage, a full knowledge of the response of Mediterranean forest ecosystems, characterized by environments with high heterogeneity, is still lacking. In short, there are many reasons why the local response of vegetation can be highly heterogeneous, and the assessment of such heterogeneity is a challenging task. In particular, local vegetation cover can be composed of different plant types (trees, shrubs, herbaceous species) which add further complexity to the total forest response. The most accurate evaluations are commonly carried out through data collected in situ, however, the high costs of field surveys do not allow the continuous collection of data in situ, especially for large areas.
Nowadays, remote sensing technologies are an essential tool to study all types of anthropic (agricultural lands, e.g., see [29,30]; seminatural (grassland, e.g., see [31])), and natural covers (shrublands, forests, prairies, e.g., [32,33,34]). Satellite remote sensing signals retrieve indirect indicators of forest health which are useful in the detection or monitoring of the stress phenomena due, for example, to drought, parasite attacks, soil pollution [35,36,37,38] however, they provide no direct measures of forest state variables such as height or diameter of trees.
The use of satellite-derived NDVI (Normalized Difference Vegetation Index) to identify areas impacted by climate events has been widely discussed in literature (e.g.,[39]). A decrease in the NDVI is expected in areas where external stress causes a decline in photosynthetic activity and damages occur, e.g., after severe fires, the NDVI time series offers helps to identify impacted zones [40,41,42,43,44].
In southern Japan, for example, the NDVI evaluation before and after a typhoon and its negative correlation with electric conductivity indicated that the reduction of NDVI values is attributable to increased salinity from flooding, proving the NDVI as an indicator of salinity damage to vegetation [45]. Moreover, a variation in NDVI was found to be the most effective in predicting total floc soil phosphorus values, due to its capacity to capture small variations in chlorophyll associated with soil total phosphorus levels in periphyton, especially in aquatic/non-impacted areas [46].
The integration of in situ and remote sensing data can improve the identification of vegetation processes (i.e., senescence, photosynthesis, evapotranspiration). Temporal, thematic, and spatial aspects for operational implementation of satellite remote-sensing data in forestry are described well in Pause et al. [47].
The main drawback of these technologies is that remote indicators combine information from spatial units (pixels), which generally cover several square meters. Low resolution removes local details and the average information obtained from these observations could not reflect differences in the response of different areas. The response of forests to adverse meteorological events is generally irregular and the damage can be non-uniform even in small areas, particularly due to the heterogeneity of the vegetation cover. Although the link between the NDVI and the amount of photosynthetic active biomass is well established and many studies testify as to its benefits in the detection of areas damaged by external stresses [39], the ability of satellite observations to provide actual quantitative estimates of the damage at the local scale cannot be assumed directly as satellite sensors retrieve spatial averaged responses.
Very high spatial resolution (a few meters), images such as WorldView, RapidEye, etc., which have been widely used for forest research [48,49] are often very expensive and unsuitable for continuous monitoring.
The launch of the new Sentinel-2 satellite, one of the core missions of the Copernicus Earth Observation Programme of the European Union, provided opportunities for new research activities into the assessment and monitoring of forests. Sentinel 2, freely available, is characterized by a wide swath (290 km), a short revisit time (5 days with 2 satellites under cloud-free conditions which results in 2–3 days at mid-latitudes), 13 spectral bands (4 bands at 10 m of spatial resolution; 6 bands at 20 m of spatial resolution, and 3 bands at 60 m of spatial resolution) at a wider scale [50] that guarantee consistent time series minimizing any artifacts introduced by atmospheric variability and allowing for the systematic monitoring of the land surface [20]. The resolution (10 m) of this sensor could be sufficient to retain most of the details of the local response of forests to external stresses accounting both for damaged crowns and understory vegetation.
This explorative work analyzes the response of forest areas to negative and prolonged meteorological conditions, in five case studies investigating the relationships between the decrease in Sentinel-2 NDVI observed in the post-event period and in situ observations concerning both the vitality of single trees and the structural layers of the vegetation species which are present in the forest communities.
Our analysis aims at evaluating the consistency between the drop in NDVI and the average vitality of the plants in different vegetation layers (trees, shrubs, herbaceous plants). Plants are firstly classified according to the Raunkieur [51] biological spectrum (percent of species belonging to each life-form in flora) and their vitality is then classified according to the Mannerucci method [52].
The final goal is to evaluate the use of Sentinel 2 data to follow the evolution over time of Mediterranean forests under the impact of climate change by retaining information on the local heterogeneity.

2. Study Area

The study areas are located in Basilicata, southern Italy, 46.8% of its territory is mountainous, 45.2% is hilly and only 8% is flat (Figure 1). In general, natural forested areas are found at high elevation, whereas man managed areas are mostly found in low-level areas. The five sites fall within areas negatively impacted by 2017 summer heatwave and represent the major forest vegetation types covering the inner part of Basilicata (mixed broadleaves stands with dominant Quercus and Carpinus genus species) within the sub-montane belt.
Gole del Marmo Platano and Costa dell’Angelo are located in the Marmo Platano-Melandro basin (municipalities of Vietri di Potenza and Savoia di Lucania) an enchanting landscape of forests and mountains of about 87,700 hectares located in the north-western part of the region.
Beech (Fagus sylvatica L.) inhabits at high altitudes, whilst oak species (Quercus cerris L. and Quercus pubescens Willd.), Carpinus orientalis Mill., and Fraxinus ornus L. characterize the sub-montane belt. Abies alba Mill. and other minor conifer species are also present in Costa dell’Angelo due to reforestation activities in the 1970s and 1980s.
The remaining sites are located in the Gallipoli Cognato-Piccole Dolomiti Lucane Park which covers over 27,000 hectares in the heart of Basilicata and includes the municipalities of Accettura, Calciano, and Oliveto Lucano in the province of Matera, and Castelmezzano and Pietrapertosa in the province of Potenza, all areas with a typically Mediterranean climate [53]. The park includes several tree species such as downy oak (Q. pubescens, Quercus ilex L.), oriental hornbeam (C. orientalis, black hornbeam (Ostrya carpinifolia Scop.), manna ash (F. ornus, maples (Acer campestre L., Acer opalus Mill., Acer monspessulanum L. ), hollies (Ilex aquifolium L.), and Turkey oak (Q. scerris ) which are the most widespread species.
According to the bioclimatic classification [54,55], the study areas belong to the Mediterranean pluviseasonal oceanic-weak semicontinental macro-bioclimate, lower subhumid ombrotype and upper meso-mediterranean termotype belt. However, local factors such as altitude and aspect determine the transitions towards the drier and warmer types.
The ecoregional classification includes the study areas into the southern Tyrrhenian section, Lucanian subsection [56].
Additional information on the field survey areas is reported in Table 1.

3. Materials and Methods

The approach we followed to evaluate the coherence between in-situ and Sentinel-2 satellite data is based on the use of NDVI data as proxy for the average photosynthetic activity at the pixel level (10 m for Sentinel-2 data). Three datasets (meteorological, remote, and in situ) were analyzed.

3.1. Meteorological Data

Meteorological data (total monthly precipitation; maximum, minimum, and average monthly temperatures) for 2016–2018 from 4 stations were provided by the Basilicata Italian Department for Civil Protection and were selected on the basis of their proximity to the study areas:
  • Balvano station is closest to the Gole del Marmo Platano site;
  • Sant’Angelo Le Fratte - Fiume Melandro station is closest to the Costa dell’Angelo site;
  • Albano di Lucania station is closest to the Castelmezzano and Pietrapertosa sites;
  • San Mauro Forte station is closest to the Palazzo site.

3.2. Satellite Data

The satellite data selected were based on the phenological phases of the forest stands analyzed: (i) early stage (April/May), (ii) middle stage (July), and (iii) late stage (September/October) of vegetation activity. We processed images acquired in 2017 and before (2016) and after (2018) the extreme climatic event.
Nine Sentinel-2 images with tile 33TWE were retrieved from the European Space Agency’s Scientific Data Hub (https://scihub.copernicus.eu/) as standard Level-1C (Top of Atmosphere reflectance) (Table 2). Data were selected on the basis of the low percentage of cloud cover in the investigated areas (completely free of clouds in the test-sites). Diffuse cloudiness in the study areas limited the availability of images in the corresponding dates for all the three years.
Satellite image processing was fundamentally based on the analysis of the NDVI, which is a proxy for plant photosynthetic activity [57,58] and a valuable tool at different spatial scales for the Mediterranean ecosystems investigated [53,59,60,61]. Several studies have established the efficiency of the NDVI as an indicator for anomalies associated with climatic changes [39,62,63].
The responsiveness of the NDVI to precipitation changes is well established at various spatial and temporal scales [27].
Drought condition effects are described by a deviation of the NDVI from their typical values, usually estimated over extended periods. Based on the concept that healthy vegetation reflects more near-infrared (NIR) light and absorbs more RED light compared to other wavelengths, the NDVI is calculated using the following formula:
N D V I = ρ N I R − ρ R E D ρ N I R + ρ R E D
where ρ N I R and ρ R E D correspond to the reflectances at the NIR and RED wavelengths, respectively [57]. For each acquired image, we calculated the NDVI using bands 8 (near-infrared channel) and 4 (red channel) at 10 m spatial resolution.
Successively, we calculated the average NDVI for each test-area by considering a buffer of 2 × 2 pixels centered on the GPS coordinates (see Table 1).
Finally, in order to capture potential differences in the phenological behavior of the vegetation we compared 2017 NDVI values with corresponding 2016 (pre-event) and 2018 (post-event) values.
Sentinel-2 data were analyzed using the Sentinel Application Platform (SNAP) v7.0.0 (available at http://step.esa.int/main/download/snap-download/) and complementary elaborations were implemented with the desktop software QGIS 3.10.4 [64]. The geographic information was projected to WGS84 UTM zone 33.

3.3. In Situ Measurements

In situ measurements were collected between September and November 2017 (post-event period) in the study areas described in Section 2 in order to evaluate the impacts on forests caused by heat and drought in the summer of 2017.
The sampling sites were randomly chosen as a result of on-site observations in different forest typologies in the Lucanian Apennine [65] which showed heat/drought stress.
The major problems were found in deciduous forests located in an altitude range of 500–800 m a.s.l. They were identified by the decline symptoms of the tree and shrub layers highlighted mainly by shoot dieback, withering and early leaf loss compared to the normal seasonal pattern. In the Central Basilicata territory, the dormancy of deciduous forests occurs from November with a drop in temperature.
In September 2017 the deciduous forests already fully manifested this condition, determined by the simultaneous occurrence of prolonged drought and particularly high temperatures during the spring-summer period.

3.3.1. The Mannerucci Method

For each site, a stand area of about 200 m2 was surveyed, identified by two geographic coordinates (center of the area). For each site, a field-form was completed with the description of: (a) number of trees; (b) tree species; (c) diameter of trees measured at 1.30 m height using a tree calliper; (d) height of the trees measured using a vertex ipsometer; (e) dendrotype to identify the origin of the trees (agamic or gamic); (f) tree vitality.
These measurements are also useful to perform the calculation of the basal area and volume, in order to evaluate the biomass present and to monitor, with subsequent measurements, variations in growth and therefore wood biomass.
The evaluation of the vitality of the vegetation was performed using the Mannerucci method [52] based on visual and qualitative analysis, assigning a vitality class, i.e., a numeric value in the range 1 (healthy plant) to 6 (dead plants), to each observed plant (Figure 2).

3.3.2. Floristic-Vegetational Sampling and Plant Functional Type Classification

The analysis of the biological spectra of the floristic-vegetational surveys and of the vulnerability classes of the species present provides a summary of the consequences of the climatic event on the structure and functional recovery capacity of the forest communities examined.
Available data for vegetational analysis (see Supplementary Materials, Table S1) were collected according to the Braun-Blanquet method [66]. Frequency and vulnerability classes were detected for each species according to the Mannerucci method (see Section 3.3.1). The average of the vulnerability class of the species observed on the site was reported to quantify the structural and functional damage to the forest community.
In order to highlight and evaluate the synecological damage and vulnerability, the structural composition of the plants communities is indicated and the functional type based on the life forms frequency is taken into consideration for each species according to the Raunkiaer method [51]. The plant functional types are considered a possible framework for predicting ecosystem response to human-induced changes such as global warming [67].
The life forms method was based on the observation of the different protection of the vegetative buds as a strategy of adaptation to environmental stress (drought or cold) and allows for the evaluation of each plant and in the plant community of the potential resistance to periods of adverse climatic conditions [68,69].
Plant life-form usually displays an obvious relationship between the plants and the environmental factors [70,71,72].
In Raunkiaer’s system [51], the more pronounced the unfavorable season the more protected the renewing buds. The effectiveness of the protection depends on the position of the perennial buds with respect to the surface of the soil during the adverse seasons (e.g., perules in the aerial buds of woody plants such as fanerophytes and chamaephyta; underground stems such as rhizomes, tubers and bulbs in geophytes; seed in therophytes).
For each site, the analysis of the functional types was performed by calculating the biological spectrum both on the percentage of the life forms and on their average frequency.
Life form for each species is derived from Pignatti [73]. Vascular plants nomenclature followed Bartolucci et al. [74].

4. Results and Discussion

4.1. Analysis of Meteorological Data

Figure 3 shows the monthly mean temperature values recorded during the period April–October in 2016, 2017, and 2018. The positive temperature anomaly observed between June and August 2017 is evident (see Supplementary Materials, Figures S1 and S2 for more details).
The similarity between thermal curves measured in different study sites reflects the non local character of this large scale climatic event that can be classified as a heatwave. This warm anomaly reached its peak in August 2017, with monthly maximum values of about 40 °C in all sites examined, more than 5 °C higher than the values measured in other years. If we consider that the elevation of the sites is higher than 500 m, this temperature appears to be really high [60]. The positive anomaly occurs between early spring and early autumn which, on the contrary, show the lowest values observed in these periods.
As for rainfalls, 2017 was drier than 2016 and 2018 for most of the year and in particular during the heatwave (Figure 4).
Severe drops in precipitation from 2016 to 2017 were observed (90% in the San Mauro Forte meteorological station). Altogether, the contextual occurrence of temperature and rainfall anomalies in summer, which is already the warmest and driest period of the year, defines the conditions observed between June and August 2017 as a climate event potentially harmful to forest health.

4.2. Analysis of Sentinel-2 NDVI

The sequence of the NDVI values estimated for the surveyed areas (Figure 5) summarizes the phenological period April to October observed by Sentinel-2 during 2016, 2017, and 2018. This sequence provides a proxy for the photosynthetic activity before, during, and after the climate event.
The most striking feature which emerges from the comparison of the plots in Figure 5 is the low level of photosynthetic activity in October 2017, after the event, which is the lowest estimated in this month over the three years. According to the NDVI values, at the end of the meteorological critical period a decrease in photosynthetic activity is detected in all sites.
The 2016 spring values are particularly high, as largely testified in literature (e.g., [39]) while the values in August are similar for all sites, except for Castelmezzano where the 2017 spring values were already low and the decrease in vegetation activity in October leads to an NDVI value that is typical of less densely vegetated areas.
Nevertheless, the magnitude of the decrease is site-specific, in other words each site experienced different negative impacts, according to satellite data.
Finally, we can hypothesize vegetation recovery as Figure 5 shows similar NDVI trends for the 2016 and 2018 summer-autumn periods for all sites and the difference between NDVI values for the 2018 autumn/summer period is around zero (Table 3).

4.3. Analysis of In Situ Data

The analysis of in situ data aims to provide a quantitative indicator of the damage to vegetation produced by drought and a heatwave for comparison with the decrease of NDVI in October.

4.3.1. Gole del Marmo-Platano Site

Thermophilous deciduous forest vegetation type, with 55 plants belonging to 4 different tree species were identified: A. monspessulanum (13%), Crategus monogina L. (2%), Q. pubescens (11%), and C.orientalis (74%). According to the Mannerucci method, they were classified as 2 (18%), 3 (29%), and 4 (53%). Each species in this area is present in class 4 (Figure 6).
The Pcaesp life form with C. orientalis dominant (cover about 70–80%) (see Supplementary Materials, Table S1- Relevè number 1) shows the greatest vulnerability confirming the results of the quantitative analysis of the damage. In the undergrowth, the shrubs (Pn, Plian) show less suffering, while the perennial herbaceous species with hypogean vegetative buds (Grhiz, Gbulb) do not report any damage at all (Figure 7).

4.3.2. Costa dell’Angelo Site

This site is characterized by 7 different tree species: A.alba (20%), C. orientalis (45%); Cedrus atlantica (Endl.) G. Manetti ex Carriére (10%); Q. ilex (2%); F. ornus (12%); P. halepensis (1%); Pinus pinaster Aiton (10%).
Trees form the structure of the forest with an evergreen canopy tree layer (height about 15 m; coverage 70%) with lower deciduous broadleaf trees and shrub layers, which showed the greatest damage (Figure 7). Among these plants, the species in the most critical conditions were C. orientalis (94% of this species belonging to class 4 of the Mannerucci method) and especially A. alba, whose young trees and seedlings were dead. The remaining plants were class 1 (40%), 2 (3%), 3 (5%), and 6 (9%) (Figure 6).
No damage was observed on the woody plants in the undergrowth (Pn, Plian) or on the herbaceous plants (Geophytes and Emicriptophytes; see Supplementary Materials, Table S1—Relevè number 2).

4.3.3. Palazzo Site

This site is currently used as pasture and the soil showedmarked erosion. The vegetation cover was discontinuous: the few trees (about 50% cover) and the scattered bushes of Q. pubescens, Spartium junceum L. and C. monogyna suffered significant damaged from the climatic event (Vulnerability class > 3) (Figure 6 and Figure 7). Two different tree species are present in this area: 27 are Q. pubescens specimen, only 1 F. ornus specimen. These species are classified as 3 (50%) and 4 (50%) according to the Mannerucci method (Figure 6).
Stress levels were intense between July and August 2017 (Figure 3 and Figure 4) and the shrubs and herbaceous species with epigean vegetative buds (Hemicriptophytes) were mainly affected (Figure 7 and Supplementary Materials, Table S1—Relevè number 3).

4.3.4. Castelmezzano Site

This site is located on a steep south-facing slope. It was an olive grove taken over by the S. junceum, a pioneering shrub characteristic of the dynamic phases of reconstitution of the oak and broad-leaved forests. The edaphic and topographical conditions increased the impact of stress for all life forms except for plants with underground buds (Geophytes; Figure 7 and Supplementary Materials, Table S1—Relevè number 4). In this site 35 trees of 5 different species were identified: A.r monspessulanum (3%); Cornus mas L. (14%); F. ornus (63%); Q. pubescens (9%); Pistacia terebinthus L. (11%). According to the Mannarucci method they are classified as 1 (20%), 2 (40%), 3 (9%), and 4 (31%). The P.terebinthus, class 1, was not affected by stress (Figure 6).

4.3.5. Pietrapertosa Site

This site is near site 4 (see Section 2, Figure 1), located on the opposite side of the valley, in an impluvium with a north-facing slope. The tree layers are sparse (covers 40%) and 63 plants of 5 different species were classified: O. carpinifolia (2%); C. orientalis (71%); F. ornus (3%); Q. pubescens (22%); P. terebinthus (2%). More than 70% is C. orientalis and according to the Mannerucci method is classified as 4 (Figure 6). High-shrub species were prevalent (Pcaesp = cover 80% and presences 58.8%, see Supplementary Materials, Table S1- Relevè number 5) and were affected by significant damage (average Vulnerability Class = 49.4%; Figure 7).

4.4. Comparative Analysis

The comparative analysis between satellite and in situ data was carried out considering the percentage of damaged plants (from class 2/3 to 6 evaluated for each vegetation layer).
Figure 8 shows the distribution by class of the tree layers independently of the species. Costa dell’Angelo is the site with less average damage, whereas the Palazzo site was the most vulnerable site with all specimens belonging to class 3 and 4.
The percentage of damaged trees (classes 2–6) appears to be proportional to the decrease in the NDVI between July and October 2017 (R2 = 0.98; Figure 9). As such, the decrease in NDVI accounts for the average impact of the climate event on the forest, Costa dell’Angelo is paradigmatic. Despite the presence of many damaged plants and some dead trees (class 6), the high percentage of unaffected trees (40%) limits the decrease of the NDVI in autumn 2017. Similarly, in the Castelmezzano site, the presence of about 20% of unaffected plants produces a decrease in the NDVI which is about twice that observed in Costa dell’Angelo. Finally, for the other three sites, where vulnerability levels are mainly distributed between classes 2, 3, and 4, the decrease of the NDVI is doubled.
The analysis of the remaining vegetation layers is consistent with Figure 9 and indeed completes the explanation of the results. In fact, in Costa dell’Angelo, the healthy condition of the dominant evergreen trees and of the undergrowth and herbaceous plants explains the modest decrease of the NDVI in autumn 2017.
The NDVI recorded very low values in Castelmezzano (Figure 5) as the soil was almost devoid of vegetation due to the widespread death of the leaf systems due to the severe and prolonged drought during September and October 2017. Moreover, in 2018, the spring rainfalls and the short period of summer drought favored the vegetative growth of less damaged plants such as Olea europea L., P. terebinthus, and F. Ornus, as highlighted by the NDVI values thereby favoring the full 2018 recovery. The decrease in the NDVI in Gole del Marmo-Platano can be mainly ascribed to tree and shrub layers, whereas the floristic richness of herbaceous layers has allowed the community to renew itself in a short time and to return in 2018 to slightly higher NDVI values than those in 2016 (Figure 5; Table 3).
The NDVI in Palazzo showed a remarkable recovery in 2018 due to the vegetative growth of the little damaged woody species (e.g., F.ornus,) and to the heliophilous herbaceous plants with hypogean vegetative buds (Geophytes; e.g., Asphodelus microcarpus Salzm. & Viv. ) not eaten by cattle (Figure 5; Table 3).
Finally, the climatic trend over the years in Pietrapertosa is similar to Castelmezzano (Figure 3 and Figure 4), however, the microclimate conditions and the remarkable floristic richness of the undergrowth facilitated efficient and effective recovery of the forest vegetation. This was confirmed by the NDVI values which had already returned to the level close to that of 2016 by autumn 2018, before the climatic stress (Figure 5).

5. Conclusions

This study investigated the consistency between in-situ data and Sentinel-2 imagery to assess the health status of forests affected by adverse climate conditions (rainfall deficit and high temperatures) in 2017, classified as one of the driest and warmest years since 1800 (https://www.ncdc.noaa.gov/sotc/national/201713;http://www.isac.cnr.it/climstor/climate_news.html).
The use of Sentinel-2 for a continuous monitoring which accounts for local damage heterogeneity was evaluated, beyond the mere identification of affected areas which is usually investigated in the current literature. Such a continuous monitoring is fundamental to follow the health status of Mediterranean forests, which are among the most vulnerable to climate risk, not only because of the expected increase of extreme events but also for possible changes in the rainfall regimes [75] that can alter the ecosystem phenological cycles. Our analyses focused on in situ and remote data separately to make independent assessments of the forest status. Overall, field surveys showed different local scale response of forests in the different sites studied. The worst case was Gole del Marmo Platano, where a high percentage of trees and shrubs was damaged. On the other hand, Costa dell’Angelo was characterized by good post-event conditions of the dominant evergreen trees and of the undergrowth and herbaceous plants. This diversity is ascribable both to different vegetation and to different magnitude of the extreme meteorological event. The 2017 spring-summer period in Gole del Marmo Platano was the driest in the area, whereas the warmest conditions were observed in the Palazzo site where, in fact, all the tree crowns were damaged.
Although drought and heatwaves are large-scale events, local differences due, e.g., to elevation and aspect can mitigate or exacerbate their impact on vegetation. In spite of these differences, the overall response of the Mediterranean forests analyzed in this work seems to be rather efficient as only a low percentage of trees, only in one site, died and the recovery was rather fast in all the area investigated.
Our results showed a good consistency between the different assessments including a separate evaluation of the health status of each vegetation layer (trees, shrubs, and herbaceous plants). The role of tree cover in determining the drop in the NDVI was striking as it is proportional to the number of damaged trees estimated using the Mannerucci method. The consistency of the contribution of the other vegetation layers was also impressive and helps to interpret both the response of the forest ecosystems and their recovery favored by the understory vegetation.
These findings suggest that the Sentinel-2 NDVI is a reliable index for the assessment of damage to vegetation at the pixel level (10 m) and shows its suitability for reporting on different levels of damage; this is an encouraging result in view of the climatic changes witnessed in recent years.
Moreover, the availability of a free and robust Sentinel-2 database and the frequent revisit time of its constellation makes these data especially suited to the continuous monitoring of forests after extreme events as they are able to provide information on large forest areas while preserving local heterogeneity. This information is useful to study the complex forest climate dynamics and to support decision makers in the management of Mediterranean forests which are particularly vulnerable to climate change.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-445X/9/9/325/s1, Figure S1: Monthly rainfalls values estimated in the ground meteorological stations from January 2016 to October 2018, Figure S2: Monthly mean temperature values estimated in the ground meteorological stations from January 2016 to October 2018, Table S1: Floristic-vegetational surveys with vulnerability class for each species.

Author Contributions

Conceptualization and methodology, F.R., M.L., and R.C.; analysis of in situ data F.R., S.F., and S.S.P.I.; analysis of meteorological and satellite data, V.I., M.L., R.C., and S.S.P.I.; writing-original draft preparation, M.L. and R.C.; writing-review and editing, F.R., V.I., M.L., R.C., S.F. and S.S.P.I.; funding acquisition, M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the MIUR (Italian Ministry of Education, University and Research), OT4CLIMA PON project.

Acknowledgments

We would like to thank the Italian Civil Protection Department of the Basilicata Region for providing meteorological data and Paula Preston, Grad.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lindner, M.; Fitzgerald, J.B.; Zimmermann, N.E.; Reyer, C.; Delzon, S.; van der Maaten, E.; Schelhaas, M.-J.; Lasch, P.; Eggers, J.; van der Maaten-Theunissen, M.; et al. Climate change and European forests: What do we know, what are the uncertainties, and what are the implications for forest management? J. Environ. Manag. 2014, 146, 69–83. [Google Scholar] [CrossRef] [Green Version]
  2. Millar, C.I.; Stephenson, N.L. Temperate forest health in an era of emerging megadisturbance. Science 2015, 349, 823–826. [Google Scholar] [CrossRef]
  3. The Intergovernmental Panel on Climate Change (IPCC). TAR Climate Change 2001: Impacts, Adaptation, and Vulnerability; Cambridge University Press: New York, NY, USA, 2001. [Google Scholar]
  4. Coumou, D.; Rahmstorf, S. A decade of weather extremes. Nat. Clim. Chang. 2012, 2, 491–496. [Google Scholar] [CrossRef]
  5. Easterling, D.R.; Meehl, G.A.; Parmesan, C.; Changnon, S.A.; Karl, T.R.; Mearns, L.O. Climate extremes: Observations, modeling, and impacts. Science 2000, 289, 2068–2074. [Google Scholar] [CrossRef] [Green Version]
  6. Ummenhofer, C.C.; Meehl, G.A. Extreme weather and climate events with ecological relevance: A review. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2017, 372. [Google Scholar] [CrossRef]
  7. Field, C.B. Climate Change 2014—Impacts, Adaptation and Vulnerability: Regional Aspects; Cambridge University Press: New York, NY, USA, 2014; ISBN 978-1-107-05816-3. [Google Scholar]
  8. Lindner, M.; Maroschek, M.; Netherer, S.; Kremer, A.; Barbati, A.; Garcia-Gonzalo, J.; Seidl, R.; Delzon, S.; Corona, P.; Kolström, M.; et al. Climate change impacts, adaptive capacity, and vulnerability of European forest ecosystems. For. Ecol. Manag. 2010, 259, 698–709. [Google Scholar] [CrossRef]
  9. Mariotti, A.; Pan, Y.; Zeng, N.; Alessandri, A. Long-term climate change in the Mediterranean region in the midst of decadal variability. Clim. Dyn. 2015, 44, 1437–1456. [Google Scholar] [CrossRef]
  10. Polade, S.D.; Gershunov, A.; Cayan, D.R.; Dettinger, M.D.; Pierce, D.W. Precipitation in a warming world: Assessing projected hydro-climate changes in California and other Mediterranean climate regions. Sci. Rep. 2017, 7, 1–10. [Google Scholar] [CrossRef]
  11. Bastida, F.; López-Mondéjar, R.; Baldrian, P.; Andrés-Abellán, M.; Jehmlich, N.; Torres, I.F.; García, C.; López-Serrano, F.R. When drought meets forest management: Effects on the soil microbial community of a Holm oak forest ecosystem. Sci. Total Environ. 2019, 662, 276–286. [Google Scholar] [CrossRef]
  12. Vicente-Serrano, S.M.; Lopez-Moreno, J.-I.; Beguería, S.; Lorenzo-Lacruz, J.; Sanchez-Lorenzo, A.; García-Ruiz, J.M.; Azorin-Molina, C.; Morán-Tejeda, E.; Revuelto, J.; Trigo, R.; et al. Evidence of increasing drought severity caused by temperature rise in southern Europe. Environ. Res. Lett. 2014, 9, 044001. [Google Scholar] [CrossRef]
  13. Dai, A. Increasing drought under global warming in observations and models. Nat. Clim. Chang. 2013, 3, 52–58. [Google Scholar] [CrossRef]
  14. Allen, C.D.; Macalady, A.K.; Chenchouni, H.; Bachelet, D.; McDowell, N.; Vennetier, M.; Kitzberger, T.; Rigling, A.; Breshears, D.D.; Hogg, E.H.T.; et al. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. For. Ecol. Manag. 2010, 259, 660–684. [Google Scholar] [CrossRef] [Green Version]
  15. Zhang, Q.; Shao, M.; Jia, X.; Wei, X. Relationship of Climatic and Forest Factors to Drought—And Heat-Induced Tree Mortality. PLoS ONE 2017, 12, e0169770. [Google Scholar] [CrossRef] [PubMed]
  16. Bréda, N.; Huc, R.; Granier, A.; Dreyer, E. Temperate forest trees and stands under severe drought: A review of ecophysiological responses, adaptation processes and long-term consequences. Ann. For. Sci. 2006, 63, 625–644. [Google Scholar] [CrossRef] [Green Version]
  17. Camarero, J.J.; Manzanedo, R.D.; Sanchez-Salguero, R.; Navarro-Cerrillo, R.M. Growth response to climate and drought change along an aridity gradient in the southernmost Pinus nigra relict forests. Ann. For. Sci. 2013, 70, 769–780. [Google Scholar] [CrossRef]
  18. Gentilesca, T.; Camele, I.; Colangelo, M.; Lauteri, M.; Lapolla, A.; Ripullone, F. Il declino dei soprassuoli di querce nel sud Italia: Il caso di studio del bosco di Gorgoglione. In Proceedings of the Atti del II Congresso Internazionale di Selvicoltura, Progettare il Futuro per il Settore Forestale, Florence, Italy, 26–29 November 2014; Volume 1, pp. 123–129. [Google Scholar]
  19. Ripullone, F.; Borghetti, M.; Raddi, S.; Vicinelli, E.; Baraldi, R.; Guerrieri, M.R.; Nolè, A.; Magnani, F. Physiological and structural changes in response to altered precipitation regimes in a Mediterranean macchia ecosystem. Trees 2009, 23, 823. [Google Scholar] [CrossRef]
  20. Puletti, N.; Mattioli, W.; Bussotti, F.; Pollastrini, M. Monitoring the effects of extreme drought events on forest health by Sentinel-2 imagery. JARS 2019, 13, 020501. [Google Scholar] [CrossRef]
  21. Adams, H.D.; Zeppel, M.J.B.; Anderegg, W.R.L.; Hartmann, H.; Landhäusser, S.M.; Tissue, D.T.; Huxman, T.E.; Hudson, P.J.; Franz, T.E.; Allen, C.D.; et al. A multi-species synthesis of physiological mechanisms in drought-induced tree mortality. Nat. Ecol. Evol. 2017, 1, 1285–1291. [Google Scholar] [CrossRef]
  22. Gazol, A.; Camarero, J.J.; Anderegg, W.R.L.; Vicente-Serrano, S.M. Impacts of droughts on the growth resilience of Northern Hemisphere forests. Glob. Ecol. Biogeogr. 2017, 26, 166–176. [Google Scholar] [CrossRef]
  23. Sánchez-Salguero, R.; Navarro-Cerrillo, R.; Camarero, J.; Fernández-Cancio, Á. Selective drought-induced decline of pine species in southeastern Spain. Clim. Chang. 2012, 113, 767–785. [Google Scholar] [CrossRef]
  24. Colangelo, M.; Camarero, J.J.; Borghetti, M.; Gentilesca, T.; Oliva, J.; Redondo, M.-A.; Ripullone, F. Drought and Phytophthora Are Associated With the Decline of Oak Species in Southern Italy. Front. Plant Sci. 2018, 9. [Google Scholar] [CrossRef] [PubMed]
  25. Arzac, A.; García-Cervigón, A.I.; Vicente-Serrano, S.M.; Loidi, J.; Olano, J.M. Phenological shifts in climatic response of secondary growth allow Juniperus sabina L. to cope with altitudinal and temporal climate variability. Agric. For. Meteorol. 2016, 217, 35–45. [Google Scholar] [CrossRef]
  26. Colangelo, M.; Camarero, J.J.; Borghetti, M.; Gazol, A.; Gentilesca, T.; Ripullone, F. Size Matters a Lot: Drought-Affected Italian Oaks Are Smaller and Show Lower Growth Prior to Tree Death. Front. Plant Sci. 2017, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Rita, A.; Camarero, J.J.; Nolè, A.; Borghetti, M.; Brunetti, M.; Pergola, N.; Serio, C.; Vicente-Serrano, S.M.; Tramutoli, V.; Ripullone, F. The impact of drought spells on forests depends on site conditions: The case of 2017 summer heat wave in southern Europe. Glob. Chang. Biol. 2020, 26, 851–863. [Google Scholar] [CrossRef] [PubMed]
  28. Gentilesca, T.; Camarero, J.J.; Colangelo, M.; Nolè, A.; Ripullone, F. Drought-induced oak decline in the western Mediterranean region: An overview on current evidences, mechanisms and management options to improve forest resilience. iForest Biogeosci. For. 2017, 10, 796. [Google Scholar] [CrossRef] [Green Version]
  29. Azzari, G.; Grassini, P.; Edreira, J.I.R.; Conley, S.; Mourtzinis, S.; Lobell, D.B. Satellite mapping of tillage practices in the North Central US region from 2005 to 2016. Remote Sens. Environ. 2019, 221, 417–429. [Google Scholar] [CrossRef]
  30. Sharma, S.; Dhakal, K.; Wagle, P.; Kilic, A. Retrospective tillage differentiation using the Landsat-5 TM archive with discriminant analysis. Agrosyst. Geosci. Environ. 2020, 3, e20000. [Google Scholar] [CrossRef] [Green Version]
  31. Chang, S.; Wu, B.; Yan, N.; Davdai, B.; Nasanbat, E. Suitability Assessment of Satellite-Derived Drought Indices for Mongolian Grassland. Remote Sens. 2017, 9, 650. [Google Scholar] [CrossRef] [Green Version]
  32. Moffet, C.A.; Taylor, J.B.; Booth, D.T. Very Large Scale Aerial (VLSA) Imagery for assessing postfire bitterbrush recovery. In Proceedings of the-Shrublands under fire: Disturbance and Recovery in a Changing World, Cedar City, UT, USA, 6–8 June 2006; Kitchen, S.G., Pendleton, R.L., Monaco, T.A., Vernon, J., Eds.; U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2008; Volume 52, pp. 161–168. [Google Scholar]
  33. Byer, S.; Jin, Y. Detecting Drought-Induced Tree Mortality in Sierra Nevada Forests with Time Series of Satellite Data. Remote Sens. 2017, 9, 929. [Google Scholar] [CrossRef] [Green Version]
  34. Sharma, S.; Ochsner, T.E.; Twidwell, D.; Carlson, J.D.; Krueger, E.S.; Engle, D.M.; Fuhlendorf, S.D. Nondestructive estimation of standing crop and fuel moisture content in tallgrass prairie. Rangel. Ecol. Manag. 2018, 71, 356–362. [Google Scholar] [CrossRef]
  35. Zhao, H.; Yang, C.; Guo, W.; Zhang, L.; Zhang, D. Automatic Estimation of Crop Disease Severity Levels Based on Vegetation Index Normalization. Remote Sens. 2020, 12, 1930. [Google Scholar] [CrossRef]
  36. Dotzler, S.; Hill, J.; Buddenbaum, H.; Stoffels, J. The Potential of EnMAP and Sentinel-2 Data for Detecting Drought Stress Phenomena in Deciduous Forest Communities. Remote Sens. 2015, 7, 14227–14258. [Google Scholar] [CrossRef] [Green Version]
  37. Bajocco, S.; De Angelis, A.; Salvati, L. A satellite-based green index as a proxy for vegetation cover quality in a Mediterranean region. Ecol. Indic. 2012, 23, 578–587. [Google Scholar] [CrossRef]
  38. D’Emilio, M.; Coluzzi, R.; Macchiato, M.; Imbrenda, V.; Ragosta, M.; Sabia, S.; Simoniello, T. Satellite data and soil magnetic susceptibility measurements for heavy metals monitoring: Findings from Agri Valley (Southern Italy). Environ. Earth Sci. 2018, 77, 63. [Google Scholar] [CrossRef]
  39. Greco, S.; Infusino, M.; De Donato, C.; Coluzzi, R.; Imbrenda, V.; Lanfredi, M.; Simoniello, T.; Scalercio, S. Late Spring Frost in Mediterranean Beech Forests: Extended Crown Dieback and Short-Term Effects on Moth Communities. Forests 2018, 9, 388. [Google Scholar] [CrossRef] [Green Version]
  40. Simoniello, T.; Coluzzi, R.; Imbrenda, V.; Lanfredi, M. Land cover changes and forest landscape evolution (1985–2009) in a typical Mediterranean agroforestry system (high Agri Valley). Nat. Hazards Earth Syst. Sci. 2015, 15, 1201–1214. [Google Scholar] [CrossRef] [Green Version]
  41. Simoniello, T.; Lanfredi, M.; Coppola, R.; Imbrenda, V.; Macchiato, M. Correlation of Vegetation and Air Temperature Seasonal Profiles—Spatial Arrangement and Temporal Variability. In Phenology Climate Change; IntechOpen: Rijeka, Croatia, 2012. [Google Scholar] [CrossRef]
  42. Imbrenda, V.; D’Emilio, M.; Lanfredi, M.; Simoniello, T.; Ragosta, M.; Macchiato, M. Integrated Indicators for the Estimation of Vulnerability to Land Degradation. In Soil Processes and Current Trends in Quality Assessment; IntechOpen: Rijeka, Croatia, 2013. [Google Scholar] [CrossRef] [Green Version]
  43. Bonfiglio, A.; Cuomo, V.; Lanfredi, M.; Macchiato, M. Interfacing NOAA/ANHRR NDVI and soil truth maps for monitoring vegetation phenology at a local scale in a heterogeneous landscape of Southern Italy. Int. J. Remote Sens. 2002, 23, 4181–4195. [Google Scholar] [CrossRef]
  44. Mancino, G.; Nolè, A.; Ripullone, F.; Ferrara, A. Landsat TM imagery and NDVI differencing to detect vegetation change: Assessing natural forest expansion in Basilicata, southern Italy. iForest Biogeosci. For. 2014, 7, 75. [Google Scholar] [CrossRef] [Green Version]
  45. Goto, K.; Goto, T.; Nmor, J.C.; Minematsu, K.; Gotoh, K. Evaluating salinity damage to crops through satellite data analysis: Application to typhoon affected areas of southern Japan. Nat. Hazards 2015, 75, 2815–2828. [Google Scholar] [CrossRef] [Green Version]
  46. Rivero, R.G.; Grunwald, S.; Binford, M.W.; Osborne, T.Z. Integrating spectral indices into prediction models of soil phosphorus in a subtropical wetland. Remote Sens. Environ. 2009, 113, 2389–2402. [Google Scholar] [CrossRef]
  47. Pause, M.; Schweitzer, C.; Rosenthal, M.; Keuck, V.; Bumberger, J.; Dietrich, P.; Heurich, M.; Jung, A.; Lausch, A. In Situ/Remote Sensing Integration to Assess Forest Health—A Review. Remote Sens. 2016, 8, 471. [Google Scholar] [CrossRef] [Green Version]
  48. McCarthy, M.J.; Jessen, B.; Barry, M.J.; Figueroa, M.; McIntosh, J.; Murray, T.; Schmid, J.; Muller-Karger, F.E. Automated High-Resolution Time Series Mapping of Mangrove Forests Damaged by Hurricane Irma in Southwest Florida. Remote Sens. 2020, 12, 1740. [Google Scholar] [CrossRef]
  49. Meiforth, J.J.; Buddenbaum, H.; Hill, J.; Shepherd, J.D.; Dymond, J.R. Stress Detection in New Zealand Kauri Canopies with WorldView-2 Satellite and LiDAR Data. Remote Sens. 2020, 12, 1906. [Google Scholar] [CrossRef]
  50. Coluzzi, R.; Imbrenda, V.; Lanfredi, M.; Simoniello, T. A first assessment of the Sentinel-2 Level 1-C cloud mask product to support informed surface analyses. Remote Sens. Environ. 2018, 217, 426–443. [Google Scholar] [CrossRef]
  51. Raunkiaer, C. The Life Forms of Plants and Statistical Plant Geography; Oxford University Press: London, UK, 1934. [Google Scholar]
  52. Mannerucci, F.; Sicoli, G. Oak decline in Apulia southern Italy: An ecological indicator of landscape evolution? In Patterns and Processes in Forest Landscapes. Consequences of Human Management; Lafortezza, R., Sanesi, G., Eds.; Locorotondo: Bari, Italy, 2006. [Google Scholar]
  53. Coluzzi, R.; D’Emilio, M.; Imbrenda, V.; Giorgio, G.A.; Lanfredi, M.; Macchiato, M.; Ragosta, M.; Simoniello, T.; Telesca, V. Investigating climate variability and long-term vegetation activity across heterogeneous Basilicata agroecosystems. Geomat. Nat. Hazards Risk 2019, 10, 168–180. [Google Scholar] [CrossRef] [Green Version]
  54. Rivas-Martinez, S.; Penas, A.; Diaz, T.E. Biogeographic map of Europe; University of Leon: Leon, Spain, 2004. [Google Scholar]
  55. Pesaresi, S.; Biondi, E.; Casavecchia, S. Bioclimates of Italy. J. Maps 2017, 13, 955–960. [Google Scholar] [CrossRef]
  56. Blasi, C.; Capotorti, G.; Copiz, R.; Guida, D.; Mollo, B.; Smiraglia, D.; Zavattero, L. Classification and mapping of the ecoregions of Italy. Plant Biosyst. Int. J. Deal. Asp. Plant Biol. 2014, 148, 1255–1345. [Google Scholar] [CrossRef]
  57. Rouse, J.W.; Deering, D.W.; Haas, R.H.; Schell, J.A. Monitoring Vegetation Systems in the Great Plains with ERTS; NASA: Washington, DC, USA, 1974; Volume 1. [Google Scholar]
  58. Tucker, C.J. Red and photographic infrared linear combinations for monitoring vegetation. Remote Sens. Environ. 1979, 8, 127–150. [Google Scholar] [CrossRef] [Green Version]
  59. Simoniello, T.; Lanfredi, M.; Liberti, M.; Coppola, R.; Macchiato, M. Estimation of vegetation cover resilience from satellite time series. Hydrol. Earth Syst. Sci. 2008, 12, 1053–1064. [Google Scholar] [CrossRef] [Green Version]
  60. Lanfredi, M.; Coppola, R.; Simoniello, T.; Coluzzi, R.; D’Emilio, M.; Imbrenda, V.; Macchiato, M. Early Identification of Land Degradation Hotspots in Complex Bio-Geographic Regions. Remote Sens. 2015, 7, 8154–8179. [Google Scholar] [CrossRef] [Green Version]
  61. Quaranta, G.; Salvia, R.; Salvati, L.; Paola, V.D.; Coluzzi, R.; Imbrenda, V.; Simoniello, T. Long-term impacts of grazing management on land degradation in a rural community of Southern Italy: Depopulation matters. Land Degrad. Dev. 2020, 2020, 1–16. [Google Scholar] [CrossRef]
  62. Piao, S.; Friedlingstein, P.; Ciais, P.; Zhou, L.; Chen, A. Effect of climate and CO2 changes on the greening of the Northern Hemisphere over the past two decades. Geophys. Res. Lett. 2006, 33. [Google Scholar] [CrossRef] [Green Version]
  63. Wardlow, B.D.; Tadesse, T.; Brown, J.F.; Callahan, K.; Swain, S.; Hunt, E. Vegetation Drought Response Index: An Integration of Satellite, Climate, and Biophysical Data; Drought Mitigation Center Faculty Publications: Boca Raton, FL, USA, 2012. [Google Scholar]
  64. QGIS Development Team QGIS Geographic Information System. Open Source Geospatial Foundation Project (2019). Available online: https://www.qgis.org/en/site/ (accessed on 28 August 2020).
  65. Costantini, G.; Bellotti, A.; Mancino, G.; Borghetti, M.; Ferrara, A. Forest map of Regione Basilicata—Atlas; INEA: Potenza, Italy, 2006; ISBN 88-8145-062-3. [Google Scholar]
  66. Braun-Blanquet, J. Pflanzensoziologie: Grundzüge der Vegetationskunde; Springer: Vienna, Austria; New York, NY, USA, 1964; ISBN 978-3-7091-4078-9. [Google Scholar]
  67. Duckworth, J.C.; Kent, M.; Ramsay, P.M. Plant functional types: An alternative to taxonomic plant community description in biogeography? Prog. Phys. Geogr. Earth Environ. 2000, 24, 515–542. [Google Scholar] [CrossRef]
  68. Godron, M.; Daget, P.; Emberger, L.; Long, G.; Le Floch, E.; Poissonet, J.; Sauvage, C.; Wacquant, J.P. Vade-mecum pour le Relevé Méthodique de la Végétation et du Milieu—CNRS Editions; Éditions du Centre National de la Recherche Scientifique: Paris, France, 1969. [Google Scholar]
  69. Mueller-Dombois, D.; Ellenberg, H. Aims and Methods of Vegetation Ecology; Wiley: New York, NY, USA, 1974; ISBN 978-0-471-62290-1. [Google Scholar]
  70. Cain, S.A. Life-Forms and Phytoclimate. Bot. Rev. 1950, 16, 1–32. [Google Scholar] [CrossRef]
  71. Fitchett, J.M.; Grab, S.W.; Thompson, D.I. Plant phenology and climate change: Progress in methodological approaches and application. Prog. Phys. Geogr. Earth Environ. 2015, 39, 460–482. [Google Scholar] [CrossRef]
  72. Parmesan, C.; Hanley, M.E. Plants and climate change: Complexities and surprises. Ann. Bot. 2015, 116, 849–864. [Google Scholar] [CrossRef]
  73. Pignatti, S. Flora d’Italia; Edagricole: Bologna, Italy, 1982. [Google Scholar]
  74. Bartolucci, F.; Peruzzi, L.; Galasso, G.; Albano, A.; Alessandrini, A.; Ardenghi, N.M.G.; Astuti, G.; Bacchetta, G.; Ballelli, S.; Banfi, E.; et al. An updated checklist of the vascular flora native to Italy. Plant Biosyst. Int. J. Deal. Asp. Plant Biol. 2018, 152, 179–303. [Google Scholar] [CrossRef]
  75. Lanfredi, M.; Coluzzi, R.; Imbrenda, V.; Macchiato, M.; Simoniello, T. Analyzing Space–Time Coherence in Precipitation Seasonality across Different European Climates. Remote Sens. 2020, 12, 171. [Google Scholar] [CrossRef] [Green Version]
Land 09 00325 g001 550
Figure 1. Study areas. Red dots represent the forest areas studied.
Figure 1. Study areas. Red dots represent the forest areas studied.
Land 09 00325 g001
Land 09 00325 g002 550
Figure 2. Tree vitality classes according to the Mannerucci method: the effects due to stress condition start to be evident from the apical shoots to the base of the branches.
Figure 2. Tree vitality classes according to the Mannerucci method: the effects due to stress condition start to be evident from the apical shoots to the base of the branches.
Land 09 00325 g002
Land 09 00325 g003 550
Figure 3. Monthly mean temperature values estimated in the ground meteorological stations selected for the analysis.
Figure 3. Monthly mean temperature values estimated in the ground meteorological stations selected for the analysis.
Land 09 00325 g003
Land 09 00325 g004 550
Figure 4. Monthly rainfalls estimated in the ground meteorological stations selected for the analysis.
Figure 4. Monthly rainfalls estimated in the ground meteorological stations selected for the analysis.
Land 09 00325 g004
Land 09 00325 g005 550
Figure 5. Sentinel-2 Normalized Difference Vegetation Index (NDVI) values estimated for spring, summer, and autumn 2016, 2017, and 2018 in the surveyed areas. On the y-axis NDVI values.
Figure 5. Sentinel-2 Normalized Difference Vegetation Index (NDVI) values estimated for spring, summer, and autumn 2016, 2017, and 2018 in the surveyed areas. On the y-axis NDVI values.
Land 09 00325 g005
Land 09 00325 g006 550
Figure 6. Tree species for each study site and their classification according to the Mannerucci method.
Figure 6. Tree species for each study site and their classification according to the Mannerucci method.
Land 09 00325 g006
Land 09 00325 g007 550
Figure 7. The biological spectrum in the forest communities was developed through the life forms of all the plants present. For each biological form the percentage of presence and the frequency are shown. The vulnerability is reported for each life form (percentage of the average of the vulnerability class according to the Mannerucci method; Supplementary Materials, Table S1).
Figure 7. The biological spectrum in the forest communities was developed through the life forms of all the plants present. For each biological form the percentage of presence and the frequency are shown. The vulnerability is reported for each life form (percentage of the average of the vulnerability class according to the Mannerucci method; Supplementary Materials, Table S1).
Land 09 00325 g007
Land 09 00325 g008 550
Figure 8. Distribution of trees by class (independently of the species) according to the Mannerucci method.
Figure 8. Distribution of trees by class (independently of the species) according to the Mannerucci method.
Land 09 00325 g008
Land 09 00325 g009 550
Figure 9. Percentage of damaged trees according to the Mannerucci method (classes 2–6) and the NDVI values between July and October 2017 (Table 3).
Figure 9. Percentage of damaged trees according to the Mannerucci method (classes 2–6) and the NDVI values between July and October 2017 (Table 3).
Land 09 00325 g009
Table
Table 1. Main features of the study areas.
Table 1. Main features of the study areas.
Study AreaCoordinatesAltitude
(m, a.s.l.)		AspectSlope (%)
Gole del Marmo Platano40°37.039′ N 500North/West30
15°27.600′ E
Costa dell’Angelo40°34.211′ N 770North/Northwest30
15°33.466′ E
Palazzo40°30.948′ N 790West15
16°8.878′ E
Castelmezzano40°32.040′ N682South/Southeast50
16°3.284′ E
Pietrapertosa40°32.015′ N 625North/Northwest35
16°3.476′E
Table
Table 2. List of the processed Sentinel-2 images and corresponding cloud cover percentage.
Table 2. List of the processed Sentinel-2 images and corresponding cloud cover percentage.
Sensing DateCloud Percentage (%)
23 May 20160
22 July 20160
30 September 20162.5
8 April 20170
12 July 20170
5 October 20170
13 April 20180.12
2 July 20180
30 September 20180.47
Table
Table 3. Difference between NDVI values for the period autumn/summer 2016, 2017, and 2018.
Table 3. Difference between NDVI values for the period autumn/summer 2016, 2017, and 2018.
NDVI Autumn Values—NDVI Summer Values
201620172018
Costa dell’Angelo0.01−0.060.05
Gole del Marmo Platano−0.03−0.190.02
Pietrapertosa−0.07−0.18−0.03
Castelmezzano−0.04−0.11−0.03
Palazzo−0.04−0.180.00

Share and Cite

MDPI and ACS Style

Coluzzi, R.; Fascetti, S.; Imbrenda, V.; Italiano, S.S.P.; Ripullone, F.; Lanfredi, M. Exploring the Use of Sentinel-2 Data to Monitor Heterogeneous Effects of Contextual Drought and Heatwaves on Mediterranean Forests. Land 2020, 9, 325. https://doi.org/10.3390/land9090325

AMA Style

Coluzzi R, Fascetti S, Imbrenda V, Italiano SSP, Ripullone F, Lanfredi M. Exploring the Use of Sentinel-2 Data to Monitor Heterogeneous Effects of Contextual Drought and Heatwaves on Mediterranean Forests. Land. 2020; 9(9):325. https://doi.org/10.3390/land9090325

Chicago/Turabian Style

Coluzzi, Rosa, Simonetta Fascetti, Vito Imbrenda, Santain Settimio Pino Italiano, Francesco Ripullone, and Maria Lanfredi. 2020. "Exploring the Use of Sentinel-2 Data to Monitor Heterogeneous Effects of Contextual Drought and Heatwaves on Mediterranean Forests" Land 9, no. 9: 325. https://doi.org/10.3390/land9090325

APA Style

Coluzzi, R., Fascetti, S., Imbrenda, V., Italiano, S. S. P., Ripullone, F., & Lanfredi, M. (2020). Exploring the Use of Sentinel-2 Data to Monitor Heterogeneous Effects of Contextual Drought and Heatwaves on Mediterranean Forests. Land, 9(9), 325. https://doi.org/10.3390/land9090325

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy