Next Article in Journal
Fuel Drivers of Fire Behaviour in Coastal Mallee Shrublands
Next Article in Special Issue
Assessment of Pinus halepensis Forests’ Vulnerability Using the Temporal Dynamics of Carbon Stocks and Fire Traits in Tunisia
Previous Article in Journal / Special Issue
Site Quality Models and Fuel Load Dynamic Equation Systems Disaggregated by Size Fractions and Vegetative States in Gorse and High Heath Shrublands in Galicia (NW Spain)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Burn Severity and Postfire Salvage Logging Effects on Vegetation and Soil System in a Short-Term Period in Mediterranean Pine Forests

by
Esther Peña-Molina
1,*,
Daniel Moya
1,
Álvaro Fajardo-Cantos
1,
Fuensanta García-Orenes
2,
Jorge Mataix-Solera
2,
Victoria Arcenegui
2,
Manuel Esteban Lucas-Borja
1 and
Jorge de las Heras
1
1
Forest Ecology Research Group (ECOFOR), High Technical School of Agricultural and Forestry Engineering and Biotechnology, University of Castilla-La Mancha, University Campus, s/n, 02071 Albacete, Spain
2
GETECMA-Soil Science and Environmental Technologies Group, Department of Agrochemistry and Environment, Miguel Hernández University of Elche, Avenida de la Universidad, s/n, 03202 Elche, Spain
*
Author to whom correspondence should be addressed.
Fire 2024, 7(4), 127; https://doi.org/10.3390/fire7040127
Submission received: 16 March 2024 / Revised: 5 April 2024 / Accepted: 8 April 2024 / Published: 9 April 2024

Abstract

:
Wildfires are a natural part of the dynamics of Mediterranean forest ecosystems. The fire patterns in the Mediterranean basin have been altered mainly due to changes in land use and climate change. In 2017, a wildfire in Yeste (Spain) burned 3200 hectares of two Mediterranean pine forests. We investigated the effects of burn severity and postfire salvage logging practices on vegetation and soil properties in four experimental areas distributed within the wildfire perimeter. These areas included unburned, low, high, and high burn severity with salvage logging, all located under Pinus halepensis Mill and Pinus pinaster Aiton stands. Salvage logging was applied 18 months after the fire. We established 72 circular plots (nine per treatment and pine species). We collected soil samples to analyze physicochemical and biological soil properties, including pH, electrical conductivity (EC), soil organic matter (SOM) content, carbon from microbial biomass (CBM), basal soil respiration (BSR), metabolic quotient (qCO2), and two enzymatic activities: β-glucosidase (GLU) and phosphatase (PHP). To understand how vegetation changed after fire, we implemented three linear transects per plot to calculate α-diversity indices (richness, Shannon, and Simpson), vegetation coverage (COBV), fraction of bare soil (BSOIL), the number of postfire seedlings (NSeed) and their average height (Hm), and we grouped vegetation into different postfire adaptive strategies: facultative seeder (R+S+), obligate resprouter (R+S−), obligate seeder (R−S+), and non-fire-adapted (R−S−). We ran ANOVA and Tukey’s HSD post hoc tests to evaluate the differences between burn severity and salvage logging practices on the variables examined for each pine stand. We used PCA and correlation analysis to identify plant-soil interactions. Our results suggest that Pinus halepensis stands were more affected by the wildfire than Pinus pinaster stands due to the distinct characteristics of each species (morphology of the leaves, bark thickness, cone structure, etc.) and the significant differences observed in terms of pH, SOM, CBM, qCO2, GLU, PHP, and Nseed. The proportion of obligate resprouter species was higher in Pinus halepensis stands, and the obligate seeder species were higher in Pinus pinaster stands. The study highlighted the importance of monitoring burn severity and postfire management practices to promote forest recovery and reduce wildfire risk. Limiting the negative impact of postfire salvage logging practices can enhance the resilience of vulnerable ecosystems.

Graphical Abstract

1. Introduction

Wildfires are a natural phenomenon in Mediterranean areas [1,2]. However, large wildfires have increased worldwide in recent years due to climate change. There are many examples of bushfires worldwide, like the Canary Islands wildfires in 2019, with more than 12,000 hectares burned; the Amazon Rainforest fires in 2019, which burned more than 800,000 hectares; the California Dixie Fire in 2021, which burned more than 380,000 hectares and more recently, the Australian bushfire season of 2019–2020, where more than 18.6 million hectares were burned. For European countries, the 2022 wildfire season was the second worst on record [3], according to European Forest Fire Information System (EFFIS) data. These high-severity wildfires severely impact the biosphere by causing deaths of species, habitat destruction, and forest loss, among other results [4]. They also decrease air quality and population health and emit large amounts of CO2 [5]. Due to increasingly hot, dry conditions, the number of these perturbations will continue to increase [6,7], and it will be necessary to improve the application of forward-looking fire management policies. These efforts should include prevention, suppression, and restoration strategies, mainly investing in the early extinction stage [8,9,10].
In this context, forest management and restoration activities are critical tools for enhanced ecological resilience in fire-prone areas [11,12,13,14]. Regeneration following a wildfire is directly related to burn severity [15,16,17,18], which is crucial to establishing practical postfire management areas [19] and can be defined as the effect of fire on the environment. It refers to the extent of the immediate fire effects on vegetation (aboveground) and soils (belowground) due to the loss or decomposition of organic matter. The composite burn index (CBI) [20] combines fire severity metrics and ecosystem responses to fire, including substrates (litter, duff, fuel, and soil), herbs, shrubs, trees, and seedling trees. It ranges from 0.0 (unburned) to 3.0 (high-severity) [21]. Burn severity directly impacts both the regeneration of plants and the recovery of the soil following wildfires by limiting natural plant regeneration and reducing soil functionality [22]. According to the review analysis conducted by Agbeshie et al. (2022) [23], soil is the most valuable natural resource. It intervenes in essential processes such as carbon sequestration, nutrient cycling, mineral storage, and plant growth sustainability. Soil organic carbon (SOC), soil reaction (pH), nitrogen (N), phosphorus (P), soil water repellency (SWR), and biological parameters are within the most fire-sensitive soil properties.
There is substantial research on forest management following wildfires and their effects on soil properties and vegetation recovery [16,23,24,25]. Forest management and restoration practices can help to reduce the impacts of burn severity. Still, it is essential to carefully consider the tradeoffs between economic objectives and fire resilience when managing large-scale forest restoration activities [26,27]. Postfire interventions encompass a range of actions that aim to mitigate the impacts of wildfires. These interventions include implementing soil erosion control measures utilizing natural and artificial materials on various scales, such as erosion barriers and check-dams, which can be deployed at hillslope and stream levels. Additionally, postfire activities may involve establishing tree plantations, salvage logging operations to remove burnt trees, or facilitating natural ecosystem regeneration. Of all the postfire restoration activities, salvage logging is one of the more contentious actions in the scientific community due to its potentially harmful effects on soil and vegetation properties [28,29,30,31,32,33]. On the one hand, salvage logging is often used to recover part of the economic value of forests and reduce insect outbreaks [34,35]. On the other hand, salvage logging can impede plant cover regeneration [36] and contribute to soil disturbance and erosion [31]. While the impacts mentioned above are typically contingent on specific contextual factors, further research is necessary to assist managers and policymakers worldwide in determining whether salvage logging will promote restoring values and processes in disturbed forests under their particular local conditions [37].
The primary goal of this study was to assess the ecological response of the ecosystem properties following a wildfire in Yeste (Spain) on 27 July 2017, which resulted in the burning of 3200 hectares with varying degrees of burn severity. The specific aims included:
(i)
Evaluating the effects of burn severity and the immediate impacts of postfire salvage logging interventions on soil system properties in the designated stands: Pinus halepensis Mill and Pinus pinaster Aiton;
(ii)
Examining vegetation population dynamics by considering two pine stands in the fire-affected perimeter. This analysis assessed the dominant tree species alongside the associated thermo-Mediterranean shrub vegetation.
In summary, this study aimed to evaluate whether salvage logging exacerbates the damage caused by wildfires in terms of soil functionality and vegetation recovery or if salvage logging has no more substantial impact than the degradation caused by wildfires. Our initial hypothesis was that burn severity and salvage logging practices would affect vegetation and soil recovery, and the different pine stands would differ in their response to fire due to their inherent characteristics and strategies in response to fire (Pinus halepensis Mill. is more serotinous, and Pinus pinaster Aiton regeneration depends more on soil seed bank) [38].

2. Materials and Methods

2.1. Study Area

The study area was in the Albacete province (Spain), in the SE region of Spain (Figure 1). On 27 July 2017, a forest fire started in the area known as “La Parrilla”, located in the Yeste municipality (Spain). After several days, the fire was finally controlled on 9 August at 5:48 p.m., and 3200 hectares had been burned. On the days when fire continued and was out of control, the atmospheric conditions were characterized by stability and dry air, influenced by the Saharan continental air mass and the high-pressure systems prevailing over the Iberian Peninsula. These conditions were only disrupted by a SW advection episode and the localized low-pressure areas that resulted from heating the air close to the land surface. Given these conditions, the predominant west and northwest wind increased in magnitude and was further aided by local winds. The situation led to higher intensity and faster fire spread on the afternoons of 28 and 30 July, estimated at around 30–35 m/min. The Castilla-La Mancha Regional Forest Service provided all this information. The fire perimeter was calculated using remote-sensing techniques, supported by fieldwork and aerial photographs taken from a helicopter [39]. The burn severity was obtained following [40].
Based on the Köppen-Geiger climate classification [41], the study area has a hot summer Mediterranean climate (“Cs” type). Climate data, including the mean annual precipitation from 400–600 mm and the mean temperatures from 10–14 °C, were obtained from the Agricultural Geographic Information System [42]. The soils in the area were classified as Aridisols and Inceptisols according to the USDA Soil Taxonomy [43]. The natural vegetation community in the study area consisted of forests dominated by Pinus halepensis Mill. and Pinus pinaster Aiton pine trees, accompanied by various shrub species, including Quercus coccifera L., Pistacia lentiscus L., Viburnum tinus L., Phyllirea angustifolia L., Arbutus unedo L., Lonicera implexa Ait., Daphne gnidium L., Rubia peregrina L., Juniperus oxycedrus L., among others.
This study followed the postfire activities undertaken by the Castilla-La Mancha Regional Forest Service following the 2017 fire. Eighteen months after the wildfire (Figure 2), the Regional Forest Service implemented several postfire emergency actions in the burned area, including stream stabilization and salvage logging. Check-dams were constructed on slopes along primary streams in publicly and privately owned areas to decrease water velocity and promote the accumulation of sediments. The check dams had a width of 6 m on either side of the barracks’ axis, measured in a straight line. Within these regions, charred plant material repair constructions were built from fallen, shortened, and fragmented tree trunks. Salvage logging practices consisted of removing tree trunks to reduce the amount of valuable timber lost and prevent insect outbreaks. No salvage logging treatment was conducted in low burn severity areas because the trees remained alive, making it unnecessary to justify extensive treatment. This type of postfire management action is typically implemented in areas of high burn severity to mitigate economic losses, minimize visual impact, prevent tree falls, and reduce the risk and severity of future fires.
Based on burn severity, pine forest (dominated by Pinus halepensis or Pinus pinaster), and postfire emergency action (salvage logging treatment or no treatment), we monitored several areas:
(1)
High-severity plots with no intervention;
(2)
High-severity plots with salvage logging treatment;
(3)
Low-severity plots with no intervention;
(4)
Unburned areas with no intervention at all.
In April 2019, the first torrential rainfall (>60 mm h−1) happened after the salvage logging treatment (Figure 2). Our study was conducted two months later, in the summer of 2019.

2.2. Experimental Design

The experimental areas (unburned, low-severity, high-severity, and high-severity with salvage logging treatment) were replicated at two pine forest sites (Pinus halepensis and Pinus pinaster stands). In each pine forest, nine plots were established per experimental area. The selected plots shared similar aspects (SE), slopes (20–40%), and soil types (loam). As a result, 72 circular plots of variable size were established (36 plots designated per pine species). The radius of plots was calculated as the ratio of the crown diameter to its base diameter from the pines selected from unburned areas. The base diameter of the Pinus pinaster plots was multiplied by 6.755 (R2 = 99.87%) to acquire the crown diameter in those areas where salvage logging was implemented, while the Pinus halepensis plots were multiplied by 2.31 (R2 = 99.33%). The center of each plot was georeferenced using a high-precision GPS device (Garmin International, Inc., Olathe, KS, USA), with a randomly selected tree serving as the center point.

2.3. Soil System Properties

According to earlier research conducted by Moya et al. (2019) and Jiménez-Morillo et al. (2020) [45,46], three soil samples were gathered from each plot in June 2019. Each sample was composed of six subsamples randomly collected from different locations in the plot. After removing surface litter, the top 2–3 cm mineral soil was excavated using a 20 × 20 cm box. Afterwards, subsamples were combined inside a plastic bag and refrigerated at 4 °C until they were sent to be analyzed.
The physicochemical characterization was based on a soil analysis conducted in the summer of 2019 to record texture, pH, electrical conductivity (EC, µS cm−1), and soil organic matter (SOM, %). Since soil quality indicators are sensitive to fire [47,48], the parameters selected in this study include the carbon from the soil microbial biomass (CBM, mg C Kg soil−1), basal soil respiration (BSR, mg C-CO2 Kg soil−1 h−1), metabolic quotient (qCO2, mg C-CO2 Kg−1 biomass C h−1) and two enzymatic activities: β-glucosidase (GLU, µmol p-NP g soil−1 h−1), and phosphatase (PHP, µmol p-NP g soil−1 h−1). BSR was measured with an automated impedance meter (BacTrac 4200 Microbiological Analyser, Sylab, Austria), using CO2 emissions by soil microorganisms at 30 °C for 24 h, as detected by indirect impedance measurement. CBM was also measured in the impedance meter as substrate-induced respiration with glucose (3 mg per gram of soil) as carbon substrate according to the Anderson and Domsch (1978) [49] method. The metabolic quotient (qCO2) was calculated as the ratio BSR:CBM.

2.4. Vegetation Index Calculations

To characterize the existing vegetation, three linear transects were set up from the center of each plot to determine α-diversity with three indices: species richness (S) as the total number of species; species abundance (Shannon Index, H) [50]; and species dominance (Simpson, D) [51]. The first transect was oriented along the maximum slope line, and the remaining two transects were spaced 120 degrees apart from the first one. The intersecting length of each species that perpendicularly touched the line along the linear transects was measured following Canfield’s linear transect method [52]. The complete species list is described in Supplementary Materials (Table S1). A visual vegetation cover estimation (COBV, %) and percentage of bare soil (BSOIL, %) were calculated on each plot. COBV was assessed using the improved Braun-Blanquet method [53]. Additionally, three square subplots (1 m2) were placed in the middle of each transect. At the square subplots (1 m2), the number of seedlings (NSeed) and their average height (Hm) were counted and measured, respectively.
Using the plant trait database for Mediterranean Basin species (BROT), we classified species according to their fire-adapted traits as facultative seeder (R+S+), obligate resprouter (R+S−), obligate seeder (R−S+), and non-fire-adapted (R−S−) [54]. The visually estimated species coverages for each trait type were merged and averaged.

2.5. Statistical Analyses

Every analysis was conducted utilizing version 4.3.2 of Rstudio [55] with a 95% level of statistical significance. We tested the normality and homoscedasticity assumptions using the Shapiro-Wilk and Levene tests. Data transformations were performed as required. For each pine forest, a one-way analysis of variance (ANOVA) test was conducted to examine the impact of burn severity and salvage logging on the soil properties and vegetation indices. We excluded the pine forest component from the ANOVA tests due to the variation in pine species composition (P. halepensis vs. P. pinaster), as both species differ in physical, ecological, and biological characteristics and temperaments, which may influence their responses to fire. A Tukey’s HSD post hoc test was run to determine which group’s means exhibited statically significant differences when comparing one another. We used a principal component analysis (PCA) to describe the variance–covariance structure of the studied variables and correlation analyses to identify the plant–soil interactions. Pearson’s correlations were utilized to analyze plant–soil interactions. Additional packages were used to generate the results: ‘dplyr’ [56] for data manipulation and ‘ggplot2’ [57] for making graphics and visualizations.

3. Results

3.1. Soil System Properties Analyses

Regarding the pH levels in the Pinus halepensis stands, there were no significant variations between the unburned and low-severity plots or between the high-severity and high-severity plots with salvage logging treatment. However, between them (unburned and low-severity vs. high-severity and high-severity with salvage logging treatment), a slight increase in pH (up to 0.2 differential points) was noted as the severity level increased. No significant changes were identified in any case for the Pinus pinaster stands.
In EC terms, no significant differences were observed for the Pinus halepensis stands. The only significant differences were observed between high-severity with salvage logging treatment and unburned plots in the Pinus pinaster stands.
No significant differences were noted for SOM between the high-severity and high-severity with salvage logging treatment plots, but important differences were found between these and the unburned and low-severity plots in the Pinus halepensis stands. In addition, a considerable reduction occurred as severity increased. The only significant differences were observed between high-severity and low-severity plots in the Pinus pinaster stands. All the results are summarized in Table 1.
Regarding CBM, no statistically significant differences were noted in the burned (low-severity, high-severity, and high-severity with salvage logging treatment) plots for both Pinus halepensis and Pinus pinaster. Significant differences appeared only between the unburned plots and all the other groups. No significant differences were observed in any stand for BSR.
Significant changes in GLU were noted between the unburned and low-severity plots compared to the high-severity and high-severity with salvage logging treatment plots in the Pinus halepensis stands. No significant variations were observed among the plots in the Pinus pinaster stands.
Regarding PHP, a significant reduction occurred as severity increased in the Pinus halepensis stands. Although this reduction was still noteworthy in the Pinus pinaster stands, it was less pronounced. All the results are summarized in Figure 3.
In Pinus halepensis stands, qCO2 showed a positive trend as burn severity increased (Figure 4). In Pinus pinaster stands, qCO2 increased with burn severity, but only the unburned and high-severity plots showed significant differences.

3.2. Vegetation Index Analyses

The COBV was higher in the unburned and low-severity plots than in the high-severity and high-severity with salvage logging treatment plots (Table 2). However, significant differences were observed only between the unburned and low-severity plots compared to the high-severity with salvage logging treatment plots in the Pinus halepensis stands. The unburned plots in the Pinus pinaster stands exhibited statistically significant differences compared to all the other locations.
At all the sites and for the pine tree species, the unburned and low-severity plots had a lower proportion of BSOIL than the high-severity and high-severity with salvage logging treatment plots.
In S terms, there were no significant variations among sites, regardless of the pine tree species. According to the vegetation indices calculations, no index showed substantial differences, irrespective of the pine tree species (Table 2).
The NSeed was significantly higher in the high-severity with salvage logging treatment plots in Pinus halepensis stands, while it was not markedly higher in Pinus pinaster stands. In both stands, unburned plots did not show pine tree regeneration. Similarly, the Hm increased with burn severity, which was higher in high-severity plots than in high-severity with salvage logging treatment, low-severity, and unburned plots. High-severity plots offered significant differences from the unburned and low-severity plots, but no significant differences were detected between the high-severity and high-severity with salvage logging treatment plots. With regard to diversity indices, no significant differences were found in any case, regardless of pine species or area (Figure 5).
Species were classified based on their specific reactions to fire (Figure 6 and Table 3). After the fire, there was a higher proportion of resprouter species (R+S−) in the Pinus halepensis stands compared to the Pinus pinaster stands. For the Pinus halepensis stands, the only significant differences were noted between the low-severity and high-severity plots, in which the proportion of obligate seeders (R−S+) increased (from 22.56% to 53.29%, respectively) and the ratio of facultative seeders decreased from 32.38% to 10.44%, respectively. In the Pinus pinaster stands, the proportion of non-fire-adapted traits (R−S−) was significantly higher in the low-severity plots (8%) than for the rest of the zones (0.92% for unburned, 0.49% for high-severity, 1.62% for high-severity with salvage logging treatment). In addition, the proportion of obligate seeders (R−S+) was lower in all fire-affected plots compared to the unburned plots, with a statistically significant reduction. No significant differences were found in the proportion of facultative seeders (R+S+).

3.3. Plant–Soil Interactions

Soil and vegetation parameters were subjected to a PCA to decrease dimensionality and retain the maximum amount of their variability, which enhanced data visualization and exploration by minimizing the number of variables involved. The dimensions were reduced to five components (with eigenvalues over 1). However, for the biplot component, we opted for a simplified model (Figure 7). The first component (Dim1) exhibited 27.5% explanatory power for variability, whereas the second component (Dim2) provided 20.6%. The variables that contributed the most to the Dim1 axis were PHP, GLU, EC, CBM, and SOM (0.867, 0.825, 0.824, 0.816, and 0.825, respectively). Conversely, H, D, S, and COBV contributed more to the Dim2 component by showing cosine squares of 0.917, 0.884, 0.882, and 0.463, respectively.
In Pinus halepensis stands, solid and significant correlations (r  >  0.5 ***, where *** denotes a significance level of 0.99) were identified on pH, EC, GLU, PHP, SOM, qCO2, BSOIL, and COBV. The pH was strongly correlated with BSOIL (r = 0.85) and negatively with EC (r = −0.92 ), GLU (r = −0.82 ), PHP (r = −0.72 ), and SOM (r = −0.68 ) (Figure 8A). Concerning the biological soil properties, the CBM was strongly correlated with PHP, GLU, and SOM (r  ≥  0.7) and negatively correlated with qCO2. The enzymatic activities were significantly correlated with BSOIL and qCO2 (r > 0.5). Regarding the vegetation properties, the S, H, and D were correlated to each other (r > 0.7). The obligate seeders (R−S+) were negatively correlated with obligate resprouter species (R+S−), COBV, and EC, and positively correlated with pH. No significant correlation was found between the rest of the physicochemical, biochemical soil properties, or vegetation indices.
The results were less significant in Pinus pinaster stands (Figure 8B). The pH was strongly correlated with PHP (r = 0.69). The EC was negatively correlated with BSOIL (r = −0.72) and positively correlated with COBV (r = 0.69). The CBM was correlated with other biological soil properties (GLU, PHP, and qCO2) but not with BSR. The GLU was also correlated with COBV (r = 0.53). Regarding vegetation indices, S, H, and D were strongly correlated with each other (r > 0.65). The obligate seeders (R−S+) were correlated with COBV (r = 0.57) and negatively correlated with obligate resprouter species (R+S−), as shown in Pinus halepensis stands.

4. Discussion

4.1. Effects of Wildfire and Salvage Logging on Soil Properties

The physicochemical soil properties exhibited minimal variation, slightly intensified by salvage logging, especially in the Pinus halepensis stands. The pH showed a positive relation with the degree of burn. This was attributed to oxides, basic cations, and carbonates in ashes, which abound following wildfires in high-severity areas [23]. Raising pH levels can reduce nutrient accessibility (i.e., phosphorus) for plants. Still, in our case, the difference in pH for Pinus halepensis stands was minimal, and no substantial alterations in pH were detected in the Pinus pinaster stands. Other studies have also revealed that pH was not affected by burning and salvage logging [58].
In contrast to the Pinus pinaster stands, EC showed no significant changes in any case in the Pinus halepensis stands, where the salvage logging effect significantly lowered the EC values (from 174.35 to 71.00 µS cm−1). These results are consistent with Muñoz-Rojas et al. (2016) [59], where EC values significantly decreased with time (5 years after the fire).
As other authors have pointed out [38,47,58], SOM content is one of the most critical soil quality indicators because of its plant growth-related functions (water retention, nutrient exchange, soil structure). In the Pinus halepensis stands, SOM content showed significant differences between the unburned plots and the other burn severities. As burn severity increased, the reduction in SOM content became more noticeable, and salvage logging amplified this effect. The decrease in SOM content was less relevant in the Pinus pinaster stands, possibly due to soil erosion. However, the salvage logging values lowered in the Pinus pinaster stands and showed significant differences with the unburned plots. In this case, salvage logging should not be applied if regeneration can be compromised.
The major effects of fire were observed principally in the microbiological soil properties. CBM was lowered at all the sites regardless of the pine tree stand being affected by the fire. Decreased microbial activity is attributed to the high sensitivity of microorganisms to heat from forest fires. The duration and intensity of wildfires can strongly impact soil microbiology [31,47]. In contrast, BSR (related to microbial activity) remained similar between the unburned and burned plots (low-severity, high-severity, and high-severity with salvage logging treatment). This phenomenon is important because a reduction in CBM while BSR remains constant indicates that microorganisms are experiencing stress and exhibiting higher respiration rates [47]. These findings are supported by Figure 4, which showed that both severity and salvage logging increased qCO2, especially in Pinus halepensis stands.
Of all the enzymatic activities, PHP and GLU showed significant reductions as severity increased, especially in the Pinus halepensis stands. These enzymatic activities are crucial for catalyzing biological processes [47]. According to other authors [22,23], these declines in GLU and PHP may be linked with reduced microbial activity (CBM).

4.2. Effects of Wildfire and Salvage Logging on Vegetation Indices

Regeneration is favored in both pine forests by high burn severity. According to the Shannon (H) and Simpson (D) indices, there were no significant differences, and salvage logging did not increase the negative impacts of wildfire. Salvage logging implied more remarkable regeneration in the Pinus halepensis plots, although salvage logging did not significantly influence the regeneration in the Pinus pinaster plots. According to the results, the high-severity plots favored pine regeneration, possibly due to poor competition and the fact that nutrients, water, and light were more available, as reported by Erdozain et al. (2023) [60]. These results favored seedling development and growth. The number of regenerated seedlings was also more prominent in the Pinus halepensis stands than those dominated by the Pinus pinaster stands in all the studied cases. This could be because Pinus halepensis is more heliophilous than Pinus pinaster and regenerates better when exposed to more sunlight [61]. In addition, the aerial seed bank, protected by serotinous pine cones, is usually more abundant in Pinus halepensis stands [16,17,62]. For both habitats (Pinus halepensis and Pinus pinaster stands), the results showed that high-severity implied lower COBV and a higher percentage of BSOIL but no differences in richness. These results are similar to those of other studies [63], which report no significant differences.
Regarding fire-adapted traits, our results showed that the proportion of resprouting species was higher in Pinus halepensis stands, while the proportion of obligate seeder species was higher in Pinus pinaster stands. Burn severity favored seedling dispersion in Pinus halepensis stands, while obligate seeder species in Pinus pinaster stands were reduced. Resprouting from remaining vegetative structures is the primary method of surviving following a wildfire [64]. As severity increased, Pinus halepensis stands showed no significant changes in (R−S−) proportions. However, there were differences in Pinus pinaster, with low-severity plots exhibiting the most significant proportions. R+S+ was minimized in high-severity plots at Pinus halepensis stands. In contrast, in Pinus pinaster stands, the (R+S+) proportion increased in tandem with the severity. Salvage logging in Pinus pinaster stands exacerbated this effect. According to Moya et al. (2020) [63], obligate seeders effectively deal with intense fires in these Mediterranean environments. In addition, resprouter species occurrence in all conditions suggests an effective technique beyond fire adaptation. Our results showed that pine forests respond positively to severe fires even in salvage logging areas, leading to more remarkable regeneration in areas with high burn severity. This was particularly evident in the case of Pinus halepensis stands. Therefore, in addition to the fire intensity, the morphological and temperamental differences between the two pine species may also impact the recovery capacity after a fire [38]. These different responses in regeneration patterns could be related to the pine serotiny level of each species. According to Cruz et al. (2019) [65], the researchers identified a significant correlation between the serotiny level and the tree’s age. If crown fires occur at a greater frequency, species resistant to fire and adapted for success in fire-prone environments could be at risk unless they adapt correspondingly [66].
In addition, salvage logging increases the possibility of mortality due to increased exposure to sunlight and direct seedling damage from mechanical treatments, which leaves more soil exposed. Furthermore, dead canopies can raise soil moisture and enhance seedling water availability while lowering soil temperature and solar radiation [67]. More studies would be necessary to analyze whether areas with recurrent fires (like the wildfire of 1994) occurring in intervals less than 25 years can affect these values.

4.3. Effects of Wildfire and Salvage Logging on Plant–Soil Interactions

The PCA results showed that the high-variance components captured significant data patterns and could extract the most relevant information. Among all the variables analyzed, the most significant ones always include PHP, GLU, EC, CBM, SOM, and COBV. These values were concentrated in non-burned areas (Dim1), while species richness and diversity indices (S, H, and D) were focused on areas with high-severity and salvage logging (Dim2). These results are consistent with those of Agbeshie et al. (2022) [23], which emphasizes that wildfires negatively affect physicochemical and biological soil properties.
The correlation matrix results provided crucial information on how soil properties relate to each other and the biodiversity present in the ecosystem. These results included pH, qCO2, BSOIL, and the variables identified before in the PCA (EC, GLU, PHP, SOM, and COBV). In the Pinus halepensis stands, COBV is related to pH, EC, SOM, GLU, PHP, and S. As mentioned before, SOM is one of the soil properties most sensitive to fire [47], and it is related to CBM and enzyme properties (PHP and GLU). In turn, CBM is related to qCO2, a ratio that can be used to identify the stress level of soil microorganisms [68]. In contrast, in the Pinus pinaster stands, SOM showed no significant relationship with any other soil property or vegetation index. We also observed important relationships between COBV, EC, CBM, GLU, and PHP. It was demonstrated that burn severity negatively impacted enzymatic activities and CBM [69].
Based on these results, analyzing how enzymatic activities (PHP and GLU) relate to recalcitrance SOM compounds would be valuable, as it could reduce soil functionality and compromise vegetation recovery. Understanding the long-term interactions between soil properties, plant biodiversity, and vegetation indices can improve the development of more effective and sustainable management strategies, the preservation of biological diversity, and the promotion of ecosystem restoration.

5. Conclusions

This study provides insights into how burn severity and salvage logging postfire management affect soil physicochemical and microbiological responses and vegetation recovery after wildfires in the short term. The high burn severity negatively impacted the ecosystem functions of soil and plant coverage. This study also demonstrated that some ecological properties recovered without being affected by salvage logging (species richness, abundance, or dominance). As burn severity increased, the decrease in SOM content became more prominent, and salvage logging exacerbated this impact. Furthermore, organisms remaining after the fire have shown elevated stress levels, particularly in areas with high burn severity and salvage logging practices. Highlighting the importance of studying the responses of different pine species to wildfires is vital for enhancing our understanding of the diverse effects observed on soil and vegetation recovery.
Nevertheless, it is essential to acknowledge the limitations of the current research before concluding. The study is constrained by the singular fire event in Yeste (Spain), which restricts the generalizability of the findings to regions with different environmental conditions and pine species. Moreover, the analysis examines the immediate impacts (18 months following the fire and salvage logging), which may not provide a comprehensive understanding of the long-term restoration of pine forest ecosystems. Specific effects, such as changes in species composition and vegetation diversity, may only become apparent later. Another challenge is that factors like fauna, complex interactions between organisms, and genetic effects may significantly influence postfire recovery beyond that of soil and vegetation characteristics. These challenges highlight the necessity of adopting a holistic and long-term approach when examining soil functionality and vegetation restoration in Mediterranean ecosystems impacted by forest fires.
In conclusion, it is imperative to continue researching the long-term effects of wildfires on Mediterranean forest ecosystems to enhance our understanding and conservation efforts.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/fire7040127/s1, Table S1: Inventory of species within the study area.

Author Contributions

Conceptualization, E.P.-M. and D.M.; methodology, E.P.-M., D.M. and Á.F.-C.; validation, E.P.-M. and Á.F.-C.; formal analysis, E.P.-M.; investigation, E.P.-M. and Á.F.-C.; resources, J.d.l.H.; data curation, E.P.-M.; writing—original draft preparation, E.P.-M.; writing—review and editing, E.P.-M., D.M., J.M.-S., F.G.-O., V.A. and M.E.L.-B.; visualization, E.P.-M. and Á.F.-C.; supervision, D.M., J.M.-S., F.G.-O., V.A., M.E.L.-B. and J.d.l.H.; project administration, J.d.l.H.; funding acquisition, J.d.l.H. All authors have read and agreed to the published version of the manuscript.

Funding

Esther Peña is supported by a doctoral fellowship from the Internal Research Plan of the University of Castilla-La Mancha, co-financed by the European Social Fund (2020-PREDUCLM-16032). This study was also funded by the Spanish Institute for Agricultural Research and Technology (INIA) with National Research Project “VIS4FIRE”: Integrated vulnerability of forest systems to wildfire: implications for forest management tools (RTA2017-0042-C05-05) and regional funds from the “Junta Comunidades Castilla-La Mancha (PRESFIRE: SBPLY/19/180501/000130/1)”, MCIN/AEI/10.13039/501100011033 “ERDF a way of making Europe” (ENFIRES: PID2020-116494RR-C43).

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Acknowledgments

The authors thank Helen Warburton for revising the English language and express their gratitude to the anonymous reviewers, whose helpful suggestions and comments improved the quality of the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Pausas, J.G.; Llovet, J.; Rodrigo, A.; Vallejo, R.; Pausas, J.G.; Llovet, J.; Rodrigo, A.; Vallejo, R. Are Wildfires a Disaster in the Mediterranean Basin?—A Review. Int. J. Wildland Fire 2009, 17, 713–723. [Google Scholar] [CrossRef]
  2. Bowman, D.M.J.S.; Kolden, C.A.; Abatzoglou, J.T.; Johnston, F.H.; van der Werf, G.R.; Flannigan, M. Vegetation Fires in the Anthropocene. Nat. Rev. Earth Environ. 2020, 1, 500–515. [Google Scholar] [CrossRef]
  3. San-Miguel-Ayanz, J.; Durrant, T.; Boca, R.; Maianti, P.; Liberta’, G.; Oom, D.; Branco, A.; De, R.D.; Ferrari, D.; Roglia, E.; et al. Advance Report on Forest Fires in Europe, Middle East and North Africa 2022. Available online: https://publications.jrc.ec.europa.eu/repository/handle/JRC133215 (accessed on 21 June 2023).
  4. Geary, W.L.; Buchan, A.; Allen, T.; Attard, D.; Bruce, M.J.; Collins, L.; Ecker, T.E.; Fairman, T.A.; Hollings, T.; Loeffler, E.; et al. Responding to the Biodiversity Impacts of a Megafire: A Case Study from South-Eastern Australia’s Black Summer. Divers. Distrib. 2022, 28, 463–478. [Google Scholar] [CrossRef]
  5. Xu, R.; Yu, P.; Abramson, M.J.; Johnston, F.H.; Samet, J.M.; Bell, M.L.; Haines, A.; Ebi, K.L.; Li, S.; Guo, Y. Wildfires, Global Climate Change, and Human Health. N. Engl. J. Med. 2020, 383, 2173–2181. [Google Scholar] [CrossRef] [PubMed]
  6. Ruffault, J.; Curt, T.; Moron, V.; Trigo, R.M.; Mouillot, F.; Koutsias, N.; Pimont, F.; Martin-StPaul, N.; Barbero, R.; Dupuy, J.-L.; et al. Increased Likelihood of Heat-Induced Large Wildfires in the Mediterranean Basin. Sci. Rep. 2020, 10, 13790. [Google Scholar] [CrossRef] [PubMed]
  7. Ellis, T.M.; Bowman, D.M.J.S.; Jain, P.; Flannigan, M.D.; Williamson, G.J. Global Increase in Wildfire Risk Due to Climate-Driven Declines in Fuel Moisture. Glob. Chang. Biol. 2022, 28, 1544–1559. [Google Scholar] [CrossRef]
  8. Alcasena, F.J.; Ager, A.A.; Bailey, J.D.; Pineda, N.; Vega-García, C. Towards a Comprehensive Wildfire Management Strategy for Mediterranean Areas: Framework Development and Implementation in Catalonia, Spain. J. Environ. Manag. 2019, 231, 303–320. [Google Scholar] [CrossRef]
  9. Duane, A.; Castellnou, M.; Brotons, L. Towards a Comprehensive Look at Global Drivers of Novel Extreme Wildfire Events. Clim. Chang. 2021, 165, 43. [Google Scholar] [CrossRef]
  10. Turco, M.; Llasat, M.-C.; von Hardenberg, J.; Provenzale, A. Climate Change Impacts on Wildfires in a Mediterranean Environment. Clim. Chang. 2014, 125, 369–380. [Google Scholar] [CrossRef]
  11. Vallejo, V.R.; Arianoutsou, M.; Moreira, F. Fire Ecology and Post-Fire Restoration Approaches in Southern European Forest Types. In Post-Fire Management and Restoration of Southern European Forests; Moreira, F., Arianoutsou, M., Corona, P., de las Heras, J., Eds.; Managing Forest Ecosystems; Springer: Dordrecht, The Netherlands, 2012; pp. 93–119. ISBN 978-94-007-2208-8. [Google Scholar]
  12. Newton, A.C.; Cantarello, E. Restoration of Forest Resilience: An Achievable Goal? New For. 2015, 46, 645–668. [Google Scholar] [CrossRef]
  13. O’Donnell, F.C.; Flatley, W.T.; Springer, A.E.; Fulé, P.Z. Forest Restoration as a Strategy to Mitigate Climate Impacts on Wildfire, Vegetation, and Water in Semiarid Forests. Ecol. Appl. 2018, 28, 1459–1472. [Google Scholar] [CrossRef] [PubMed]
  14. Alayan, R.; Rotich, B.; Lakner, Z. A Comprehensive Framework for Forest Restoration after Forest Fires in Theory and Practice: A Systematic Review. Forests 2022, 13, 1354. [Google Scholar] [CrossRef]
  15. Arianoutsou, M.; Ne’eman, G. Post-Fire Regeneration of Natural Pinus halepensis Forests in the East Mediterranean Basin. In Ecology, Biogeography and Management of Pinus halepensis and P. brutia Forest Ecosystems in the Mediterranean Basin; Backhuys Publishers: Leiden, The Netherlands, 2000; pp. 269–290. [Google Scholar]
  16. De las Heras, J.; Moya, D.; Vega, J.A.; Daskalakou, E.; Vallejo, V.R.; Grigoriadis, N.; Tsitsoni, T.; Baeza, J.; Valdecantos, A.; Fernández, C.; et al. Post-Fire Management of Serotinous Pine Forests. In Post-Fire Management and Restoration of Southern European Forests; Moreira, F., Arianoutsou, M., Corona, P., de las Heras, J., Eds.; Managing Forest Ecosystems; Springer: Dordrecht, The Netherlands, 2012; pp. 121–150. ISBN 978-94-007-2208-8. [Google Scholar]
  17. Fernández-García, V.; Fulé, P.Z.; Marcos, E.; Calvo, L. The Role of Fire Frequency and Severity on the Regeneration of Mediterranean Serotinous Pines under Different Environmental Conditions. For. Ecol. Manag. 2019, 444, 59–68. [Google Scholar] [CrossRef]
  18. Rodríguez-García, E.; Santana, V.M.; Alloza, J.A.; Ramón Vallejo, V. Predicting Natural Hyperdense Regeneration after Wildfires in Pinus halepensis (Mill.) Forests Using Prefire Site Factors, Forest Structure and Fire Severity. For. Ecol. Manag. 2022, 512, 120164. [Google Scholar] [CrossRef]
  19. Keeley, J.E. Fire Intensity, Fire Severity and Burn Severity: A Brief Review and Suggested Usage. Int. J. Wildland Fire 2009, 18, 11. [Google Scholar] [CrossRef]
  20. Key, C.H.; Benson, N.C. Landscape Assessment: Ground measure of severity, the Composite Burn Index; and Remote sensing of severity, the Normalized Burn Ratio. In FIREMON: Fire Effects Monitoring and Inventory System; USDA Forest Service, Rocky Mountain Research Station: Ogden, UT, USA, 2006. [Google Scholar]
  21. Lutes, D.C.; Keane, R.E.; Caratti, J.F.; Key, C.H.; Benson, N.C.; Sutherland, S.; Gangi, L.J. FIREMON: Fire Effects Monitoring and Inventory System. In General Technical Report RMRS-GTR-164-CD; U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station. 1 CD: Fort Collins, CO, USA, 2006; p. 164. [Google Scholar] [CrossRef]
  22. Moya, D.; González-De Vega, S.; Lozano, E.; García-Orenes, F.; Mataix-Solera, J.; Lucas-Borja, M.E.; de las Heras, J. The Burn Severity and Plant Recovery Relationship Affect the Biological and Chemical Soil Properties of Pinus halepensis Mill. Stands in the Short and Mid-Terms after Wildfire. J. Environ. Manag. 2019, 235, 250–256. [Google Scholar] [CrossRef] [PubMed]
  23. Agbeshie, A.A.; Abugre, S.; Atta-Darkwa, T.; Awuah, R. A Review of the Effects of Forest Fire on Soil Properties. J. For. Res. 2022, 33, 1419–1441. [Google Scholar] [CrossRef]
  24. Girona-García, A.; Vieira, D.C.S.; Silva, J.; Fernández, C.; Robichaud, P.R.; Keizer, J.J. Effectiveness of Post-Fire Soil Erosion Mitigation Treatments: A Systematic Review and Meta-Analysis. Earth-Sci. Rev. 2021, 217, 103611. [Google Scholar] [CrossRef]
  25. López-Vicente, M.; Cerdà, A.; Kramer, H.; Keesstra, S. Post-Fire Practices Benefits on Vegetation Recovery and Soil Conservation in a Mediterranean Area. Land. Use Policy 2021, 111, 105776. [Google Scholar] [CrossRef]
  26. Gómez-Sánchez, E.; de las Heras, J.; Lucas-Borja, M.; Moya, D. Assessing fire severity in semi-arid environments: Application in Donceles 2012 wildfire (SE Spain). Rev. Teledetección 2017, 103–113. [Google Scholar] [CrossRef]
  27. Gómez Sánchez, M.E.; Lucas Borja, M.E.; Plaza Álvarez, P.A.; González Romero, J.; Sagra Cozar, J.; Moya Navarro, D.; de las Heras Ibáñez, J. Efecto de los trabajos de restauración forestal post-incendio en ladera sobre la recuperación de la funcionalidad del suelo. Cuad. Soc. Española Cienc. For. 2019, 45, 35–44. [Google Scholar] [CrossRef]
  28. Donato, D.C.; Fontaine, J.B.; Campbell, J.L.; Robinson, W.D.; Kauffman, J.B.; Law, B.E. Post-Wildfire Logging Hinders Regeneration and Increases Fire Risk. Science 2006, 311, 352. [Google Scholar] [CrossRef]
  29. Heil, L.J.; Burkle, L.A. The Effects of Post-Wildfire Salvage Logging on Plant Reproductive Success and Pollination in Symphoricarpos Albus, a Fire-Tolerant Shrub. For. Ecol. Manag. 2019, 432, 157–163. [Google Scholar] [CrossRef]
  30. Povak, N.A.; Churchill, D.J.; Cansler, C.A.; Hessburg, P.F.; Kane, V.R.; Kane, J.T.; Lutz, J.A.; Larson, A.J. Wildfire Severity and Postfire Salvage Harvest Effects on Long-Term Forest Regeneration. Ecosphere 2020, 11, e03199. [Google Scholar] [CrossRef]
  31. García-Carmona, M.; García-Orenes, F.; Mataix-Solera, J.; Roldán, A.; Pereg, L.; Caravaca, F. Salvage Logging Alters Microbial Community Structure and Functioning after a Wildfire in a Mediterranean Forest. Appl. Soil. Ecol. 2021, 168, 104130. [Google Scholar] [CrossRef]
  32. Leverkus, A.B.; Buma, B.; Wagenbrenner, J.; Burton, P.J.; Lingua, E.; Marzano, R.; Thorn, S. Tamm Review: Does Salvage Logging Mitigate Subsequent Forest Disturbances? For. Ecol. Manag. 2021, 481, 118721. [Google Scholar] [CrossRef]
  33. Mantero, G.; Morresi, D.; Negri, S.; Anselmetto, N.; Lingua, E.; Bonifacio, E.; Garbarino, M.; Marzano, R. Short-Term Drivers of Post-Fire Forest Regeneration in the Western Alps. Fire Ecol. 2023, 19, 23. [Google Scholar] [CrossRef]
  34. Müller, J.; Noss, R.F.; Thorn, S.; Bässler, C.; Leverkus, A.B.; Lindenmayer, D. Increasing Disturbance Demands New Policies to Conserve Intact Forest. Conserv. Lett. 2019, 12, e12449. [Google Scholar] [CrossRef]
  35. Leverkus, A.B.; Gustafsson, L.; Lindenmayer, D.B.; Castro, J.; Rey Benayas, J.M.; Ranius, T.; Thorn, S. Salvage Logging Effects on Regulating Ecosystem Services and Fuel Loads. Front. Ecol. Environ. 2020, 18, 391–400. [Google Scholar] [CrossRef]
  36. Leverkus, A.B.; Lorite, J.; Navarro, F.B.; Sánchez-Cañete, E.P.; Castro, J. Post-Fire Salvage Logging Alters Species Composition and Reduces Cover, Richness, and Diversity in Mediterranean Plant Communities. J. Environ. Manag. 2014, 133, 323–331. [Google Scholar] [CrossRef]
  37. Leverkus, A.B.; Gustafsson, L.; Rey Benayas, J.M.; Castro, J. Does Post-Disturbance Salvage Logging Affect the Provision of Ecosystem Services? A Systematic Review Protocol. Environ. Evid. 2015, 4, 16. [Google Scholar] [CrossRef]
  38. Lucas-Borja, M.E.; Plaza-Álvarez, P.A.; González-Romero, J.; Miralles, I.; Sagra, J.; Molina-Peña, E.; Moya, D.; de las Heras, J.; Fernández, C. Post-Wildfire Straw Mulching and Salvage Logging Affects Initial Pine Seedling Density and Growth in Two Mediterranean Contrasting Climatic Areas in Spain. For. Ecol. Manag. 2020, 474, 118363. [Google Scholar] [CrossRef]
  39. De Santis, A.; Chuvieco, E. GeoCBI: A Modified Version of the Composite Burn Index for the Initial Assessment of the Short-Term Burn Severity from Remotely Sensed Data. Remote Sens. Environ. 2009, 113, 554–562. [Google Scholar] [CrossRef]
  40. Viedma, O.; Chico, F.; Fernández, J.J.; Madrigal, C.; Safford, H.D.; Moreno, J.M. Disentangling the Role of Prefire Vegetation vs. Burning Conditions on Fire Severity in a Large Forest Fire in SE Spain. Remote Sens. Environ. 2020, 247, 111891. [Google Scholar] [CrossRef]
  41. Kottek, M.; Grieser, J.; Beck, C.; Rudolf, B.; Rubel, F. World Map of the Köppen-Geiger Climate Classification Updated. Meteorol. Z. 2006, 15, 259–263. [Google Scholar] [CrossRef] [PubMed]
  42. SIGA. Available online: https://sig.mapama.gob.es/siga/ (accessed on 30 November 2022).
  43. Soil Survey Staff. Keys to Soil Taxonomy, 13th ed.; USDA-Natural Resources Conservation Service: Washington, DC, USA, 2022. [Google Scholar]
  44. AEMET Weather Data from Yeste Station (Albacete-Spain) 2016–2021 State Meteorological Agency of Spanish Government (AEMET), Madrid, Spain. Available online: https://datosclima.es/Aemet2013/Precipitastad2013.php (accessed on 20 February 2024).
  45. Moya, D.; Certini, G.; Fulé, P.Z.; Moya, D.; Certini, G.; Fulé, P.Z. Fire Regime and Ecosystem Responses: Adaptive Forest Management in a Changing World (Part 2). Int. J. Wildland Fire 2019, 28, 471–472. [Google Scholar] [CrossRef]
  46. Jiménez-Morillo, N.T.; Almendros, G.; de la Rosa, J.M.; Jordán, A.; Zavala, L.M.; Granged, A.J.P.; González-Pérez, J.A. Effect of a Wildfire and of Post-Fire Restoration Actions in the Organic Matter Structure in Soil Fractions. Sci. Total Environ. 2020, 728, 138715. [Google Scholar] [CrossRef]
  47. Mataix-Solera, J.; Guerrero, C.; García-Orenes, F.; Bárcenas, G.M.; Torres, M.P. Forest Fire Effects on Soil Microbiology. In Fire Effects on Soils and Restoration Strategies; CRC Press: Boca Raton, FL, USA, 2009; pp. 149–192. [Google Scholar]
  48. Moya, D.; González-De Vega, S.; García-Orenes, F.; Morugán-Coronado, A.; Arcenegui, V.; Mataix-Solera, J.; Lucas-Borja, M.E.; de las Heras, J. Temporal Characterisation of Soil-Plant Natural Recovery Related to Fire Severity in Burned Pinus halepensis Mill. Forests. Sci. Total Environ. 2018, 640–641, 42–51. [Google Scholar] [CrossRef]
  49. Anderson, J.P.E.; Domsch, K.H. A Physiological Method for the Quantitative Measurement of Microbial Biomass in Soils. Soil. Biol. Biochem. 1978, 10, 215–221. [Google Scholar] [CrossRef]
  50. Shannon, C.E. A Mathematical Theory of Communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef]
  51. Simpson, E.H. Measurement of Diversity. Nature 1949, 163, 688. [Google Scholar] [CrossRef]
  52. Canfield, R. Application of the Line Interception Method in Sampling Range Vegetation. J. For. 1941, 39, 388–394. [Google Scholar]
  53. van der Maarel, E. Transformation of Cover-Abundance Values for Appropriate Numerical Treatment-Alternatives to the Proposals by Podani. J. Veg. Sci. 2007, 18, 767–770. [Google Scholar] [CrossRef]
  54. Paula, S.; Arianoutsou, M.; Kazanis, D.; Tavsanoglu, Ç.; Lloret, F.; Buhk, C.; Ojeda, F.; Luna, B.; Moreno, J.M.; Rodrigo, A.; et al. Fire-Related Traits for Plant Species of the Mediterranean Basin. Ecology 2009, 90, 1420. [Google Scholar] [CrossRef]
  55. Posit Team. RStudio: Integrated Development Environment for R. Posit Software, PBC, Boston, MA. Available online: http://www.posit.co/2024 (accessed on 31 January 2024).
  56. Wickham, H.; François, R.; Henry, L.; Müller, K.; Vaughan, D. Dplyr: A Grammar of Data Manipulation, R package version 1.1.4; Posit Software, PBC: Boston, MA, USA, 2023. [Google Scholar]
  57. Wickham, H. Ggplot2; Use R! Springer International Publishing: Cham, Switzerland, 2016; ISBN 978-3-319-24275-0. [Google Scholar]
  58. Lucas-Borja, M.E.; Ortega, R.; Miralles, I.; Plaza-Álvarez, P.A.; González-Romero, J.; Peña-Molina, E.; Moya, D.; Zema, D.A.; Wagenbrenner, J.W.; de las Heras, J. Effects of Wildfire and Logging on Soil Functionality in the Short-Term in Pinus halepensis M. Forests. Eur. J. For. Res. 2020, 139, 935–945. [Google Scholar] [CrossRef]
  59. Muñoz-Rojas, M.; Erickson, T.E.; Martini, D.; Dixon, K.W.; Merritt, D.J. Soil Physicochemical and Microbiological Indicators of Short, Medium and Long Term Post-Fire Recovery in Semi-Arid Ecosystems. Ecol. Indic. 2016, 63, 14–22. [Google Scholar] [CrossRef]
  60. Erdozain, M.; Bonet, J.A.; Martínez de Aragón, J.; de-Miguel, S. Forest Thinning and Climate Interactions Driving Early-Stage Regeneration Dynamics of Maritime Pine in Mediterranean Areas. For. Ecol. Manag. 2023, 539, 121036. [Google Scholar] [CrossRef]
  61. Ribeiro, S.; Cerveira, A.; Soares, P.; Fonseca, T. Natural Regeneration of Maritime Pine: A Review of the Influencing Factors and Proposals for Management. Forests 2022, 13, 386. [Google Scholar] [CrossRef]
  62. Verkaik, I.; Espelta, J.M. Post-Fire Regeneration Thinning, Cone Production, Serotiny and Regeneration Age in Pinus halepensis. For. Ecol. Manag. 2006, 231, 155–163. [Google Scholar] [CrossRef]
  63. Moya, D.; Sagra, J.; Lucas-Borja, M.E.; Plaza-Álvarez, P.A.; González-Romero, J.; de Las Heras, J.; Ferrandis, P. Post-Fire Recovery of Vegetation and Diversity Patterns in Semiarid Pinus halepensis Mill. Habitats after Salvage Logging. Forests 2020, 11, 1345. [Google Scholar] [CrossRef]
  64. Keeley, J.E.; Pausas, J.G. Evolutionary Ecology of Fire. Annu. Rev. Ecol. Evol. Syst. 2022, 53, 203–225. [Google Scholar] [CrossRef]
  65. Cruz, O.; García-Duro, J.; Casal, M.; Reyes, O. Role of Serotiny on Pinus Pinaster Aiton Germination and Its Relation to Mother Plant Age and Fire Severity. iForest-Biogeosci. For. 2019, 12, 491. [Google Scholar] [CrossRef]
  66. Romero, B.; Ganteaume, A. Does Recent Fire Activity Impact Fire-Related Traits of Pinus halepensis Mill. and Pinus sylvestris L. in the French Mediterranean Area? Ann. For. Sci. 2020, 77, 106. [Google Scholar] [CrossRef]
  67. Castro, J.; Allen, C.D.; Molina-Morales, M.; Marañón-Jiménez, S.; Sánchez-Miranda, Á.; Zamora, R. Salvage Logging Versus the Use of Burnt Wood as a Nurse Object to Promote Post-Fire Tree Seedling Establishment. Restor. Ecol. 2011, 19, 537–544. [Google Scholar] [CrossRef]
  68. Panico, S.C.; Memoli, V.; Santorufo, L.; Aiello, S.; Barile, R.; de Marco, A.; Maisto, G. Soil Biological Responses under Different Vegetation Types in Mediterranean Area. Int. J. Environ. Res. Public Health 2022, 19, 903. [Google Scholar] [CrossRef]
  69. Fernández-García, V.; Miesel, J.; Baeza, M.J.; Marcos, E.; Calvo, L. Wildfire Effects on Soil Properties in Fire-Prone Pine Ecosystems: Indicators of Burn Severity Legacy over the Medium Term after Fire. Appl. Soil. Ecol. 2019, 135, 147–156. [Google Scholar] [CrossRef]
Figure 1. (A) Location of the Albacete Province in Spain; (B) the study area with the Yeste wildfire (black) and burn severity; (C) Pinus pinaster Aiton stands; (D) Pinus halepensis Mill. stands. Plots on each stand are classified as unburned, low-burn severity, high-burn severity, high-burn severity with salvage logging treatment, and background with the Google satellite. The reference coordinate system was ETRS89 (EPSG:25830).
Figure 1. (A) Location of the Albacete Province in Spain; (B) the study area with the Yeste wildfire (black) and burn severity; (C) Pinus pinaster Aiton stands; (D) Pinus halepensis Mill. stands. Plots on each stand are classified as unburned, low-burn severity, high-burn severity, high-burn severity with salvage logging treatment, and background with the Google satellite. The reference coordinate system was ETRS89 (EPSG:25830).
Fire 07 00127 g001
Figure 2. Historical intensity of precipitation (mm h−1) for the research area (Yeste, Spain) from 2016 to 2021 [44]. The year of the fire is marked in red; the salvage logging treatment is colored in green, and the field sampling is in yellow.
Figure 2. Historical intensity of precipitation (mm h−1) for the research area (Yeste, Spain) from 2016 to 2021 [44]. The year of the fire is marked in red; the salvage logging treatment is colored in green, and the field sampling is in yellow.
Fire 07 00127 g002
Figure 3. Soil biological properties for the unburned, low-severity, high-severity, and high-severity with salvage logging treatment sites. The significant differences between treatments (HSD, p < 0.05) are indicated by letters above the group label. Having the same letter in the group label suggests a similar response of means. The hanging bars in the graph correspond to standard deviations.
Figure 3. Soil biological properties for the unburned, low-severity, high-severity, and high-severity with salvage logging treatment sites. The significant differences between treatments (HSD, p < 0.05) are indicated by letters above the group label. Having the same letter in the group label suggests a similar response of means. The hanging bars in the graph correspond to standard deviations.
Fire 07 00127 g003
Figure 4. Metabolic quotient (qCO2, mg C-CO2 Kg−1 biomass C h−1) for the unburned, low-severity, high-severity, and high-severity with salvage logging treatment sites. The significant differences between treatments (HSD, p < 0.05) are indicated by letters above the group label. Having the same letter in the group label suggests a similar response of means. The hanging bars in the graph correspond to standard deviations.
Figure 4. Metabolic quotient (qCO2, mg C-CO2 Kg−1 biomass C h−1) for the unburned, low-severity, high-severity, and high-severity with salvage logging treatment sites. The significant differences between treatments (HSD, p < 0.05) are indicated by letters above the group label. Having the same letter in the group label suggests a similar response of means. The hanging bars in the graph correspond to standard deviations.
Fire 07 00127 g004
Figure 5. Vegetation indices for the unburned, low-severity, high-severity, and high-severity with salvage logging sites. The significant differences between treatments (HSD, p < 0.05) are indicated by letters above the group label. Having the same letter in the group label suggests a similar response of means. The hanging bars in the graph correspond to standard deviations.
Figure 5. Vegetation indices for the unburned, low-severity, high-severity, and high-severity with salvage logging sites. The significant differences between treatments (HSD, p < 0.05) are indicated by letters above the group label. Having the same letter in the group label suggests a similar response of means. The hanging bars in the graph correspond to standard deviations.
Fire 07 00127 g005
Figure 6. Average coverages according to the postfire trait adaptive strategy: facultative seeder (R+S+), obligate resprouter (R+S−), obligate seeder (R−S+), and non-fire-adapted (R−S−) in the study area: unburnt plots, low-severity plots, high-severity plots, and high-severity plots with salvage logging treatment.
Figure 6. Average coverages according to the postfire trait adaptive strategy: facultative seeder (R+S+), obligate resprouter (R+S−), obligate seeder (R−S+), and non-fire-adapted (R−S−) in the study area: unburnt plots, low-severity plots, high-severity plots, and high-severity plots with salvage logging treatment.
Fire 07 00127 g006
Figure 7. The principal component analysis (PCA) biplot diagram in the Pinus halepensis and P. pinaster stands. S: species richness; H: Shannon Index; D: Simpson Index; (R−S−): non-fire-adapted trait; (R+S−): obligate seeder; (R−S+): obligate resprouter; (R+S+): facultative seeder; BSOIL: bare soil; COBV: vegetation coverage; CBM: carbon from soil microbial biomass; BSR: basal soil respiration; SOM: soil organic matter; GLU: β-glucosidase enzymatic activity; PHP: phosphatase enzymatic activity.
Figure 7. The principal component analysis (PCA) biplot diagram in the Pinus halepensis and P. pinaster stands. S: species richness; H: Shannon Index; D: Simpson Index; (R−S−): non-fire-adapted trait; (R+S−): obligate seeder; (R−S+): obligate resprouter; (R+S+): facultative seeder; BSOIL: bare soil; COBV: vegetation coverage; CBM: carbon from soil microbial biomass; BSR: basal soil respiration; SOM: soil organic matter; GLU: β-glucosidase enzymatic activity; PHP: phosphatase enzymatic activity.
Fire 07 00127 g007
Figure 8. Correlation matrix among physicochemical, biological, and vegetation properties for the unburned, low-severity, high-severity, and high-severity with salvage logging plots in Pinus halepensis Mill. stands (A) and Pinus pinaster Aiton stands (B) (n = 36). Significant correlations are represented with asterisks. Abbreviations: electrical conductivity (EC); organic matter (SOM, %); carbon from the soil microbial biomass (CBM, mg C Kg soil−1); basal soil respiration (BSR, mg C-CO2 Kg soil−1 h−1), metabolic quotient (qCO2, mg C-CO2 Kg−1 biomass C h−1); β-glucosidase (GLU, µmol p-NP g soil−1 h−1); phosphatase (PHP, µmol p-NP g soil−1 h−1); species richness (S); Shannon Index (H); Simpson index (D); vegetation cover estimation (COBV, %); percentage of bare soil (BSOIL, %); facultative seeder (R+S+); obligate resprouter (R+S−); obligate seeder (R−S+); and non-fire-adapted (R−S−). The symbols *, **, and *** refer to the levels of statistical significance of 0.90, 0.95, and 0.99, respectively.
Figure 8. Correlation matrix among physicochemical, biological, and vegetation properties for the unburned, low-severity, high-severity, and high-severity with salvage logging plots in Pinus halepensis Mill. stands (A) and Pinus pinaster Aiton stands (B) (n = 36). Significant correlations are represented with asterisks. Abbreviations: electrical conductivity (EC); organic matter (SOM, %); carbon from the soil microbial biomass (CBM, mg C Kg soil−1); basal soil respiration (BSR, mg C-CO2 Kg soil−1 h−1), metabolic quotient (qCO2, mg C-CO2 Kg−1 biomass C h−1); β-glucosidase (GLU, µmol p-NP g soil−1 h−1); phosphatase (PHP, µmol p-NP g soil−1 h−1); species richness (S); Shannon Index (H); Simpson index (D); vegetation cover estimation (COBV, %); percentage of bare soil (BSOIL, %); facultative seeder (R+S+); obligate resprouter (R+S−); obligate seeder (R−S+); and non-fire-adapted (R−S−). The symbols *, **, and *** refer to the levels of statistical significance of 0.90, 0.95, and 0.99, respectively.
Fire 07 00127 g008
Table 1. Physicochemical soil properties (VAR); pH, electrical conductivity (EC, µS cm−1), and soil organic matter (SOM, %); concerning tree species (SP): Pinus halepensis Mill. and Pinus pinaster Aiton at different levels of burn severity: unburned plots, low-severity plots, high-severity plots, and high-severity plots with salvage logging. Each box includes the mean values ± standard deviations. Groups that do not share a letter differed significantly (p < 0.05).
Table 1. Physicochemical soil properties (VAR); pH, electrical conductivity (EC, µS cm−1), and soil organic matter (SOM, %); concerning tree species (SP): Pinus halepensis Mill. and Pinus pinaster Aiton at different levels of burn severity: unburned plots, low-severity plots, high-severity plots, and high-severity plots with salvage logging. Each box includes the mean values ± standard deviations. Groups that do not share a letter differed significantly (p < 0.05).
SPVARUnburnedLow-SeverityHigh-SeverityHigh-Severity with Salvage Logging TreatmentFPr (>F)
Pinus halepensispH8.07 ± 0.07 b8.04 ± 0.05 b8.2 ± 0.01 a8.28 ± 0.01 a14.5400.013
EC195.85 ± 7.28 a246.95 ± 82.05 a148.90 ± 4.94 a154.5 ± 34.64 a2.0550.249
SOM10.09 ± 2.37 c7.24 ± 1.82 b4.17 ± 0.73 a5.24 ± 1.00 a28.5200.000
Pinus pinasterpH8.12 ± 0.15 a7.99 ± 0.08 a7.24 ± 0.94 a7.72 ± 0.31 a14.5400.013
EC174.35 ± 9.82 a118.05 ± 7.42 ab137.30 ± 34.64 ab71.00 ± 5.37 b2.0550.249
SOM6.73 ± 1.37 ab4.79 ± 1.50 b7.26 ± 2.19 a5.50 ± 2.32 ab3.1890.037
Table 2. Vegetation variables (VAR): plant coverage (COBV, %), bare soil (BSOIL, %), vegetation richness (S), number of pine seedlings (NSeed), and average height (Hm, m) concerning tree species (SP): Pinus halepensis Mill. and Pinus pinaster Aiton at different levels of burn severity: unburned plots, low-severity plots, high-severity plots, and high-severity plots with salvage logging. Each box includes the mean values ± standard deviations. Groups that do not share a letter differ significantly (p < 0.05).
Table 2. Vegetation variables (VAR): plant coverage (COBV, %), bare soil (BSOIL, %), vegetation richness (S), number of pine seedlings (NSeed), and average height (Hm, m) concerning tree species (SP): Pinus halepensis Mill. and Pinus pinaster Aiton at different levels of burn severity: unburned plots, low-severity plots, high-severity plots, and high-severity plots with salvage logging. Each box includes the mean values ± standard deviations. Groups that do not share a letter differ significantly (p < 0.05).
SPVARUnburnedLow-SeverityHigh-SeverityHigh-Severity with Salvage Logging TreatmentFPr (>F)
Pinus halepensisCOBV137.02 ± 41.73 a134.64 ± 50.41 a100.06 ± 34.49 ab88.13 ± 33.31 b4.9510.006
BSOIL0 ± 0 b3.60 ± 5.59 b35.42 ± 14.65 a39.55 ± 19.30 a35.7300.000
S5.56 ± 1.19 a5.08 ± 1.17 a5.81 ± 1.80 a6.33 ± 2.61 a0.2760.842
NSeed0 ± 0 b1 ± 1 b1 ± 1 b3 ± 2 a6.5900.0009
Hm0 ± 0 b5.69 ± 3.23 b15.20 ± 8.69 a13.60 ± 5.84 a15.530.0000
Pinus pinasterCOBV214.2 ± 32.55 a77.78 ± 40.71 b104.02 ± 49.93 b73.09 ± 24.55 b11.4800.000
BSOIL0.79 ± 2.05 b0 ± 0 b13.59 ± 17.13 b44.01 ± 14.75 a35.7600.000
S7.42 ± 1.25 a6.44 ± 2.57 a5.72 ± 2.03 a5.72 ± 1.42 a1.3670.271
NSeed0 ± 0 a1 ± 1 a1 ± 1 a1 ± 1 a1.1400.3430
Hm0 ± 0 c5.67 ± 7.16 bc16.20 ± 14.00 a15.10 ± 6.18 ab6.8980.0007
Table 3. Average coverages according to the postfire trait adaptive strategy: facultative seeder (R+S+), obligate resprouter (R+S−), obligate seeder (R−S+), and non-fire-adapted (R−S−) in the study area: unburned plots, low-severity plots, high-severity plots, and high-severity plots with salvage logging treatment. Each box includes the mean values ± standard deviations. Groups that do not share a letter differ significantly (p < 0.05).
Table 3. Average coverages according to the postfire trait adaptive strategy: facultative seeder (R+S+), obligate resprouter (R+S−), obligate seeder (R−S+), and non-fire-adapted (R−S−) in the study area: unburned plots, low-severity plots, high-severity plots, and high-severity plots with salvage logging treatment. Each box includes the mean values ± standard deviations. Groups that do not share a letter differ significantly (p < 0.05).
SPVARUnburnedLow-SeverityHigh-SeverityHigh-Severity with Salvage Logging TreatmentFPr (>F)
Pinus halepensis(R−S−)1.38 ± 2.56 a1.36 ± 1.70 a1.48 ± 4.54 a2.81 ± 4.58 a0.1420.934
(R−S+)34.60 ± 16.99 ab22.56 ± 13.78 b53.29 ± 21.48 a41.02 ± 17.86 ab3.9550.0166
(R+S−)43.88 ± 18.70 a43.70 ± 17.40 a34.79 ± 23.18 a32.30 ± 16.21 a0.8230.491
(R+S+)20.14 ± 16.46 ab32.38 ± 11.66 a10.44 ± 7.85 b23.87 ± 17.65 ab4.3790.0108
Pinus pinaster(R−S−)0.92 ± 1.43 b8.69 ± 8.00 a0.49 ± 1.10 b1.62 ± 3.88 b5.4240.0039
(R−S+)75.34 ± 9.99 a47.18 ± 19.52 b42.71 ± 21.11 b53.24 ± 19.25 b5.530.0036
(R+S−)7.79 ± 5.66 a18.29 ± 17.89 a24.21 ± 24.35 a12.99 ± 8.66 a2.290.0971
(R+S+)15.95 ± 10.49 a25.84 ± 13.34 a32.59 ± 20.45 a32.15 ± 19.07 a0.8620.471
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Peña-Molina, E.; Moya, D.; Fajardo-Cantos, Á.; García-Orenes, F.; Mataix-Solera, J.; Arcenegui, V.; Lucas-Borja, M.E.; de las Heras, J. Burn Severity and Postfire Salvage Logging Effects on Vegetation and Soil System in a Short-Term Period in Mediterranean Pine Forests. Fire 2024, 7, 127. https://doi.org/10.3390/fire7040127

AMA Style

Peña-Molina E, Moya D, Fajardo-Cantos Á, García-Orenes F, Mataix-Solera J, Arcenegui V, Lucas-Borja ME, de las Heras J. Burn Severity and Postfire Salvage Logging Effects on Vegetation and Soil System in a Short-Term Period in Mediterranean Pine Forests. Fire. 2024; 7(4):127. https://doi.org/10.3390/fire7040127

Chicago/Turabian Style

Peña-Molina, Esther, Daniel Moya, Álvaro Fajardo-Cantos, Fuensanta García-Orenes, Jorge Mataix-Solera, Victoria Arcenegui, Manuel Esteban Lucas-Borja, and Jorge de las Heras. 2024. "Burn Severity and Postfire Salvage Logging Effects on Vegetation and Soil System in a Short-Term Period in Mediterranean Pine Forests" Fire 7, no. 4: 127. https://doi.org/10.3390/fire7040127

APA Style

Peña-Molina, E., Moya, D., Fajardo-Cantos, Á., García-Orenes, F., Mataix-Solera, J., Arcenegui, V., Lucas-Borja, M. E., & de las Heras, J. (2024). Burn Severity and Postfire Salvage Logging Effects on Vegetation and Soil System in a Short-Term Period in Mediterranean Pine Forests. Fire, 7(4), 127. https://doi.org/10.3390/fire7040127

Article Metrics

Back to TopTop