Effects of postfire restoration interventions on oak forest recovery
We found a tendency for an increase in NDVI over time, which indicates a generalized postfire oak forest recovery in the study areas. NDVI maximum annual values increased during the spring meteorological season, which corresponds to the period between the start of the Mediterranean growing season (March) and the maximum of photosynthetic activity (June) (Piedallu et al. 2019). We also observed a shift in annual maximum NDVI values amid the second and the third year (based on partial effects model with values below 0.03 in the first year and over 0.1 after the second year). Such trends agree with other studies that show a gradual recovery of oak forest communities, dominated by the recovery of herbaceous species during the first 2 years, with the recovery of woody species (shrubs and trees) becoming stronger from the third year onwards (Calvo et al. 2003; Capitanio and Carcaillet 2008).
Postfire restoration interventions showed a significant positive impact on oak forest recovery, however with a small model effect. The small effect observed may be the result of the diversity of postfire restoration interventions (e.g., emergency or medium/long-term interventions) and corresponding execution periods included in the selected projects. Indeed, the type of postfire restoration treatment and its moment of implementation is extremely important since distinct interventions and moments of implementation may promote opposite impacts. For example, emergency interventions such as protection of the burned soil to avoid erosion (such as mulching or barriers) should be implemented a few months after the fire to promote an effective restoration and help to prevent fire negative impacts (Vega et al. 2013; Fernández and Vega 2014; Bontrager et al. 2019; Fernández et al. 2019; Fernández 2021; Girona-García et al. 2021), hence facilitating the recovery of burned vegetation (Calvo et al. 2003). However, the execution of interventions outside the optimal timeframe is common in projects with public funding in Portugal (Lopes et al. 2022), and can have a negative impact on vegetation, including inhibition of vegetation recovery (Bautista et al. 1996; Alloza et al. 2016).
Another example is salvage logging, which is a common postfire restoration technique implemented in the medium and long-term in the Mediterranean Basin, in order to remove burnt trees for subsequent reforestation and to generate income (Leverkus et al. 2012). Nevertheless, this intervention creates disturbances in the vegetation, and may lead to changes in the ecosystem structure and functions (Bautista et al. 2010; Morgan et al. 2015; García-Orenes et al. 2017). The magnitude of the effects depends on the type, intensity, frequency and spatial pattern of logging, combined with the characteristics of the fire and post-fire weather conditions (Lindenmayer et al. 2008), with possible negative impacts on vegetation regeneration, loss of organic matter or increased soil erosion and compaction (Bautista et al. 2010). Nevertheless, in some cases there has been no negative impact on vegetation recovery (Fernández and Vega 2016). Particularly in Portugal, 40% of managed oak stands are clear-cut after the fire, of which 24% are cut within the first 3 months and 57% between 3 and 6 months after fire (Sousa 2011). Sowing or planting are also interventions executed in the medium-long term after fire (some months or years after the fire) but may delay vegetation recovery since rates of survival and growth of artificially sowed/planted trees are lower than survival and growth of trees originated from natural regeneration, which starts a few weeks after the fire in the case of broadleaved species (Moreira et al. 2009b; Catry et al. 2010a).
Effects of postfire drought events and fire characteristics on oak forest recovery
The drought index PDSI was the variable with the higher effect in the postfire recovery of NDVI. During the analysed period, PDSI monthly values were mostly negative (reflecting mild to extreme drought levels), with severe drought months observed in the first months after the fire. In fact, the hydrological year 2016–2017 was the 9th driest year since 1931 and in 2017–2018 most of mainland Portugal was under severe to extreme drought. In addition, this drought event was distinct from former drought years since drought severity was aggravated in the autumn of 2017 (IPMA 2020). Our results showed that postfire recovery of NDVI was negatively affected by severe drought levels and positively affected by increased wetness, in accordance with available literature. Indeed, the initial stages of the oak life cycle are especially vulnerable to the reduction of precipitation (Montagnoli et al. 2016; Marañón et al. 2020), and severe droughts will likely affect postfire regeneration capacity of oak forests (Acácio et al. 2017; Marañón et al. 2020). A recent study also showed that precipitation deficits were associated with changes from deciduous oak forests to other land cover types in Portugal (Acácio et al. 2017). In general, post-fire climatic conditions have been pointed out as one of the most important predictive factors for postfire vegetation recovery (Pausas et al. 1999; Torres et al. 2018; Nolan et al. 2021).
Regarding burn severity, our results showed that NDVI responded positively to both low and to high burn severities, while intermediate severities showed a reduced effect on NDVI recovery rate. Although both extreme categories appear to have a similar outcome, the NDVI response can be justified by different reasons. On one hand, low burn-severity fires will, in general, provide beneficial consequences, since trees will not be top-killed, surface litter will be partially consumed, and the soil organic layer will remain largely intact (Keeley 2009; Fernandes et al. 2010). This will result in availability of soil nutrients and habitat rejuvenating, favouring certain plants and life forms (Castro Rego et al. 2021), since plant mortality caused by low-moderate severity fires is highly selective and dependent of plant species and individual size (Marzano et al. 2012).
On the other hand, higher burn severity levels (registered by more than 50% of analysed points) will cause higher impacts on the soil and vegetation when compared to lower severity fires (Viana-Soto et al. 2017; Castro Rego et al. 2021). High-severity fires (often crown fires) originate near total mortality of the above-ground vegetation, including a significant amount of post-fire stem mortality in the deciduous oak trees (Catry et al. 2010b; Catry et al. 2013), and total consumption of the forest floor (Keeley 2009; Tepley et al. 2018), and also lead to changes in ecosystem recovery dynamics, decreasing seed dispersal, viability, and success rate (Marzano et al. 2012). Therefore, high burn severities likely imply that vegetation will require more time to recover (Keeley 2009; Vega et al. 2013; Tepley et al. 2018). As such, NDVI increase under high burn severity is partially explained by the recovery of the understorey vegetation that dominates the burned area (Castro Rego et al. 2021) and also by the greater canopy severity (top kill) that stimulates basal budding of trees (Moreira et al. 2009a; Frelich et al. 2015). Such differences in the recovery of vegetation cannot be detected by NDVI measurements, which show vegetation greenness and do not guarantee that the same type of pre-fire vegetation is being regenerated (Gouveia et al. 2010; Meneses 2021). In agreement with our results, postfire NDVI values in pine forests mixed with oaks and other woody species under moderate/high burn severities increased quickly over time, achieving those of the lightly burned areas two years after the fire (Lee and Chow 2015).
Lastly, number of fires showed a small negative effect on NDVI recovery rate, which decreased with more than 6 fires. Indeed, recurrent fire events may lead to a decrease in oak dominance in favour of other tree species, with the magnitude of this effect being dependent on the number of fires (Burton et al. 2010; Frelich et al. 2015). Furthermore, higher fire frequency may lead to a diminishing capacity of oaks to resprout (Frelich et al. 2015). Until a certain frequency of fires, its effect is increasingly negative, by reducing sapling density, with Pedunculate oak and Pyrenean oak saplings being especially sensitive to fire frequency (Monteiro-Henriques and Fernandes 2018), and by altering the stand structure (Burton et al. 2010). For example, repeated burning can transform an uneven-aged oak stand into an even-aged stand (Knapp et al. 2017). On the contrary, and as our results indicate, when fire frequency is above a certain threshold, it leads to lower fuel loading and lower burn severity, and consequently lower impacts in oak forest recovery (Burton et al. 2010; Steel et al. 2015; Knapp et al. 2017).
Future research needs
Spectral indices are increasingly used in fire ecology (Szpakowski and Jensen 2019) and the majority of the authors considers NDVI as an extremely valuable tool, applicable without additional field validation (Gitas et al. 2012). Spectral indices also show relatively high accuracies of estimation. For example, the index dNBR showed a global accuracy of 81% to estimate burn severity with Sentinel-2 satellite imagery, in comparison with in situ measurements (Sobrino et al. 2019). Nevertheless remote sensing data are approximations of reality (Szpakowski and Jensen 2019). There is a wide range of possible indices available, each with its own strengths and weaknesses which will have an impact on the results (Fernández-Guisuraga et al. 2018; João et al. 2018; Szpakowski and Jensen 2019). In addition, remote sensing data are generally limited by atmospheric effects, saturation, and sensor characteristics, and as such, a careful interpretation of results is recommended (Huang et al. 2021). Nevertheless, newer techniques are approaching maturity (e.g. hyperspectral imagery or Unmanned Aerial Vehicles to obtain ultra-high resolution imagery), which will enable more accurate collection of data (Fernández-Guisuraga et al. 2018; Veraverbeke et al. 2018; van Gerrevink and Veraverbeke 2021), and hence better models.
Our study was also limited by lack of available data. We used data from postfire restoration projects submitted for public funding, which do not specify the interventions implemented or the period of execution (Lopes et al. 2022). In addition, projects data was only available for the most recent years; a longer period of analysis may have led to clearer results. Furthermore, available spatial data on oak land cover for the analysed period aggregates deciduous oak species into a single class, which can also include other broadleaved species, without detailing the dominant oak species or the level of canopy cover. As already discussed, distinct species composition or dominance may alter the overall rate of vegetation recovery (Botequim et al. 2017), and impact the results.
Future changes in fire regimes are expected worldwide as a result of climate change, in particular for southern Europe (e.g., higher frequency and severity) (Dupuy et al. 2020). Although resprouting species such as oaks are likely the most resilient to changing fire regimes, postfire oak recruitment and forest recovery may be hindered by overlapping climate-driven stressors such as pest outbreaks, droughts and heatwaves, which may significantly limit the capacity of vegetation to recover in the future (Nolan et al. 2021). It is therefore necessary to continue studying the impacts of postfire restoration interventions on forest recovery and integrate the possible ecological consequences of restoration into the decision-making process (Robichaud et al. 2009), and accept that, for example, non-interventions can be the most adequate postfire management strategy for the Mediterranean broadleaved forests (Carrari et al. 2022). Ultimately, a better understanding of postfire vegetation dynamics will lead to better decisions in the future regarding postfire forest management.