Assessing Canopy Responses to Thinnings for Sweet Chestnut Coppice with Time-Series Vegetation Indices Derived from Landsat-8 and Sentinel-2 Imagery

Forest management treatments often translate into changes in forest structure. Understanding and assessing how forests react to these changes is key for forest managers to develop and follow sustainable practices. A strategy to remotely monitor the development of the canopy after thinning using satellite imagery time-series data is presented. The aim was to identify optimal remote sensing Vegetation Indices (VIs) to use as time-sensitive indicators of the early response of vegetation after the thinning of sweet chestnut (Castanea Sativa Mill.) coppice. For this, the changes produced at the canopy level by different thinning treatments and their evolution over time (2014&ndash:2019) were extracted from VI values corresponding to two trials involving 33 circular plots (r = 10 m). Plots were subjected to one of the following forest management treatments: Control with no intervention (2800&ndash:3300 stems ha&minus:1), Treatment 1, one thinning leaving a living stock density of 900&ndash:600 stems ha&minus:1 and Treatment 2, a more intensive thinning, leaving 400 stems ha&minus:1. Time series data from Landsat-8 and Sentinel-2 were collected to calculate values for different VIs. Canopy development was computed by comparing the area under curves (AUCs) of different VI time-series annually throughout the study period. Soil-Line VIs were compared to the Normalized Vegetation Index (NDVI) revealing that the Second Modified Chlorophyll Absorption Ratio Index (MCARI2) more clearly demonstrated canopy evolution tendencies over time than the NDVI. MCARI2 data from both L8 and S2 reflected how the influence of treatment on the canopy cover decreases over the years, providing significant differences in the thinning year and the year after. Metrics derived from the MCARI2 time-series also demonstrated the capacity of the canopy to recovery to pretreatment coverage levels. The AUC method generates a specific V-shaped time-signature, the vertex of which coincides with the thinning event and, as such, provides forest managers with another tool to assist decision making in the development of sustainable forest management strategies.