Next Article in Journal
Sensitivity Analysis of Influencing Factors of Fire Smoke Transport on Subway Station Platforms
Previous Article in Journal
Forest Fire Object Detection Analysis Based on Knowledge Distillation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Low-Severity Fire on Soil Physico-Chemical Properties in an Andean Páramo in Southern Ecuador

by
Vinicio Carrión-Paladines
1,*,
Andreas Fries
2,
María Belén Hinojosa
3,
Andy Oña
1,
Leticia Jiménez Álvarez
1,
Ángel Benítez
1,
Fausto López Rodríguez
1 and
Roberto García-Ruiz
4
1
Departamento de Ciencias Biológicas, Universidad Técnica Particular de Loja, San Cayetano Alto s/n C.P., Loja 11-01-608, Ecuador
2
Departamento de Geología, Minas e Ingeniería Civil (DGMIC), Universidad Técnica Particular de Loja, San Cayetano Alto s/n, Loja 11-01-608, Ecuador
3
Departamento de Ciencias Ambientales, Universidad de Castilla-La Mancha, Campus Fábricas de Armas, 45071 Toledo, Spain
4
Departamento de Biología Animal, Biología Vegetal y Ecología, Sección de Ecología, Universidad de Jaén, Campus Las Lagunillas, 23071 Jaén, Spain
*
Author to whom correspondence should be addressed.
Submission received: 9 October 2023 / Revised: 12 November 2023 / Accepted: 16 November 2023 / Published: 22 November 2023

Abstract

:
The high Andean páramos (AnP) are unique ecosystems that harbor high biodiversity and provide important ecosystem services, such as water supply and regulation, as well as carbon sequestration. In southern Ecuador, this ecosystem is threatened by anthropogenic burning activities to create pastures and agricultural land. However, knowledge of the effects of fire on soil properties and nutrient availability is still limited. This study conducted an experimental burn with different ignition patterns on an AnP plateau in southern Ecuador. Fire behavior (flame height, propagation speed, temperature reached on the soil), and fire severity were evaluated. In addition, soil samples were collected at 10 cm depth both 24 h and one year after the burns to measure the effects of fire on the main physico-chemical properties. The results indicate that the low severity of the experimental burns did not affect the physico-chemical properties of the soils, and therefore, soil quality was not altered. These results can help decision makers in the design of policies, regulations, and proposals for the conservation and environmental restoration of AnPs affected by wildfires in southern Ecuador.

1. Introduction

Andean páramos (AnPs) are one of the world’s most rapidly evolving biodiversity hotspots [1]. These ecosystems are found at altitudes ranging from 3000 to 3800 m above sea level and are characterized by treeless vegetation and a diversity of lakes and peatlands [2]. AnPs are found in the South American Andean Mountain ranges, including Colombia (14,434 km2), Venezuela (2660 km2), Ecuador (12,600 km2), and northern Peru (4200 km2) [3]. They are also found in small areas of the Cordillera de Talamanca in Costa Rica and Panama (80 km2) [4,5]. Due to their altitudinal and topographic positioning, the AnPs are adapted to a wide range of environmental conditions [3]. AnPs are key to the capture, regulation, and supply of water for nearby areas [6]. This is due to the hydromorphic characteristics of the soil that allow them to retain water, as they develop in water-saturated conditions due to the constant precipitation and relatively low evapotranspiration of this ecosystem throughout the year [7,8]. In addition, these soils are characterized by low bulk density, high porosity, and very friable consistency, which limits their use for agriculture activities and makes them vulnerable to soil compaction by continuous trampling of livestock [9]. According to their physical characteristics, páramo soils have many properties in common with peatlands [2,6]. They can contain up to 44% organic matter, which can reach up to 90% in peatlands [10].
The climate at these high elevations is cold and humid, with sudden fog and drizzle, and rapid changes in temperature, solar radiation, and humidity during a day [2,11]. The cold and humid climate favors the accumulation of organic matter in the soil, which plays an important role in carbon storage.It explains its high porosity and microporosity compared to other soil types [12,13], and it is recognizable by its dark humic soil layer [14]. This accumulation of organic matter and carbon lead to a high water regulation capacity of the AnP soils, which is also facilitated by the topography containing many local concavities and depressions forming peatlands (wetlands) and small lakes [9]. The high water retention capacity of these soils helps to reduce downstream flood risks and to ensure a constant river flow for domestic, industrial, and agricultural uses throughout the year [15,16].
AnPs are also fragile ecosystems that face threats from human activities such as conversion to pasture for livestock production [17,18]. These ecosystems are home to many unique plant species; however, a significant portion of them are lost due to environmental changes brought on by rural development [19]. According to Buytaert et al. [9], AnPs house approximately 5000 plant species, with 60% being endemic, showcasing their exceptional biological diversity [20]. The native species are highly adapted to the prevailing edaphic and climatic conditions, often causing irreversible impacts to the ecosystems [21,22].
In Ecuador, AnPs occupy about 7% of the country’s territory [23]. There are six different types of AnPs in Ecuador. In the northern and central regions, there are the herbaceous, the dry, the espeletias (with thick, hairy leaves that allow them to adapt to extreme conditions), and the cushion-forming plant (plants in the form of compact cushions that protect themselves from the cold and wind and retain water) páramos. In the southern region, there are herbaceous and shrubby páramos [24]. Despite their limited geographical extent, AnPs face significant anthropogenic pressures, including burning practices to create pasture and agricultural land [25]. These burning practices have historical roots [26,27], but the decline in traditional techniques has led to wildfires spreading more frequently to other ecosystems [28]. However, despite research regarding the hydrological capacity of soils [29,30], fire ecology [31], and traditional burning practices [32], knowledge about the effects of fire on soil properties and nutrient availability in AnPs is still limited, and this is because the behavior of fire in this ecosystem is unknown. Moreover, fire intensity and severity depend on vegetation type, fire climatology, soil type, and soil moisture content at the time of burning, which also determine ash production and quality [33,34].
In general, high-, medium-, and low-severity fires can occur when vegetation is burned. Of these, low-severity fires can have positive effects by accumulating organic matter in the soil due to the incomplete combustion of biomass [35]. However, the ashes produced could also have impacts on soil properties, such as affecting bulk density, as they can clog soil pores, reducing its water regulation capacity [33,36].
The objective of this study was to evaluate the effect of low-severity fire on the physico-chemical soil properties in an AnP in southern Ecuador through an experimental burn. For this purpose, a meteorological station (San Lucas_UTPL) was installed and permanent sampling plots were established. In addition, different ignition patterns and techniques were applied to produce different fire severities in the AnP and to analyze their effects on physico-chemical soil properties.

2. Materials and Methods

2.1. Study Area

The experimental burn area is located within the natural AnP of southern Ecuador, where the Saraguro indigenous parish resides (latitude 3°43′38.20″ S; longitude 79°12′30.25″ W; mean elevation 3050 m asl) (Figure 1a). The ecosystem can be classified as a humid montane forest with the presence of herbaceous páramo at higher elevations (Figure 1b) [24]. The soils are identified as Umbric Ferralsol, characterized by a surface layer rich in organic matter, an acidic pH, and a high iron and aluminum content [37,38]. The páramo has a regular topography due to its position within a high mountain plateau, as with other AnPs in Central and South America [39]. This AnP is part of the northwestern buffer zone of the protected area “Corazón de Oro”, which is managed by the local Saraguro community.
The vegetation of this herbaceous páramo is composed of a variety of plants and grasses [6]. The main plant species include páramo straw (Calamagrostis effusa), fleshy-rooted plants such as Eryngium humile, Xenophyllum humile, and Paepalanthus spp., achupallas rosettes (Espeletia sp.), blackberries (Rubus ulmifolius), miconias (Miconia calvescens), puya (Puya hamata L.B.Sm.), páramo straw (Calamagrostis intermedia [J.Presl] Steud.), páramo flower (Oritrophium crocifolium [Kunth] Cuatrec.), naurapo (Myrteola nummularia [Poir.] O. Berg.), sirius (Xyris subulata Ruiz & Pav.), and joyapas (Disterigma rimbachii). These plants serve as food sources for a variety of mammals such as the spectacled bear (Tremarctos ornatus) [40,41].
The herbaceous páramo provides crucial ecosystem services, with the most significant benefit being the water supply for the local population [32]. The main settlement in this region, San Lucas, is located at an altitude of 2800 m above sea level and has a temperate climate with average temperatures around 13.5 °C and an annual rainfall range of 600 to 1000 mm [41]. The primary sources of livelihood for the local population are agriculture, livestock, forestry, and artisan products [42]. Agricultural production mainly consists of crops of corn (Zea mays), beans (Phaseolus vulgaris), broad beans (Vicia faba L.), potatoes (Solanum tuberosum L.), and peas (Pisum sativum L.), grown for self-consumption or the local market, as well as medicinal plants and flowers [41]. The continued expansion of agricultural land and livestock production are the main environmental challenges in the area due to the land use changes through slash-and-burn activities [41,42]. Moreover, the páramo vegetation is frequently burned to regenerate the grass vegetation, which serves as a food source for domestic animals, particularly for cattle (Bos taurus) [32].

2.2. Experimental Design

The experimental burning site is located on a high mountain plateau with slopes ranging from 0 to 20%. We established nine experimental plots in an area of the herbaceous páramo facing southeast with a slope of 15% (Figure 2). The selection of this site was based on the knowledge of the Saraguros indigenous people, who have vast experience in fire management in this area, which also allowed us to determine the oldest area in which a fire occurred and therefore the area of greatest fuel accumulation [32].
Individual plots were 4 m wide and 20 m long, with an area of 80 m2. The individual plots were spaced 3 m apart, or 8 m, according to the scheme shown in Figure 2, to prevent fire spread. The plant material that was collected during plot delimitation was used as additional fuel, which was evenly distributed in the plots intended for controlled burning. This addition was justified by the fact that this technique reflects the traditional burning practices used in the region. In these practices, the Saraguro indigenous people build mounds with dried grass, typical of the herbaceous vegetation of the páramo, which are later incinerated [32]. Therefore, no additional material was incorporated in the control plots, as has been reported in previous studies [43]. The plots were protected with wooden posts and barbed wire to prevent access by cattle (Bos taurus) that occasionally graze in the area and native species such as the spectacled bear (Tremarctos ornatus) [41] (Figure 1). The main purpose was to protect natural regeneration, avoiding browsing and trampling in the experimental plots (Figure 2a).
The first plots of each block were burned from the bottom to the top (uphill) to create a frontal fire behavior (higher fire intensity and speed), while the second plots of each block were burned from the top to the bottom (downhill) to produce a backing fire behavior (lower fire intensity and speed) [44,45,46]. In addition, the third plot was kept unburned as a control plot (Figure 2a,b). A strip burn technique was applied for ignition (Ig), with the starting line set 20 cm inside the plot to prevent fire spreading into the cut areas [47]. The ignition was induced using a mixture of diesel (3/4) and gasoline (1/4) in a 5 L fuel container. All plots were ignited using a single starting line (Ig, Figure 2b).

2.3. Evaluation of Fire Behavior and Severity during Experimental Burns

The experimental burning was conducted in 2021 during the Veranillo del Niño phenomenon (VdN), a period of dry and sunny weather, which typically occurs for approximately 15 days between October and November [32,48]. Meteorological conditions were monitored using a Vantage Pro Plus automatic weather station from Davis Instruments USA (San Lucas_UTPL), located close to the experimental site [49]. The monitored meteorological variables included temperature (°C), relative humidity (RH%), solar radiation (W/m2), wind speed (m/s), and wind direction (°), which were recorded at 15 min intervals and transmitted in real time to the Technical University of Loja (Universidad Técnica Particular de Loja (UTPL)). The experimental burning was executed in all the plots to be burned on the same day between 12 p.m. and 2 p.m., as recommended by Geron and Hays [50], which is characterized as the time when the highest values of solar radiation (W/m2) and temperature (°C) are typically observed, as well as the lowest relative humidity (%).
Before the experimental burning, the fine dead fuel moisture was calculated by applying the procedures outlined in the Interagency Fire Use Module Field Guide [51], which is widely used internationally. Fuel load was calculated by random sampling in the AnP area near the experimental site but outside the plots, using 10 wood squares of 1 m2 each. The squares were randomly distributed, and all herbaceous vegetation was cut with a sickle at ground level. The samples were preserved in paper bags and transported to the laboratory of the UTPL, where they were dried in an oven for 48 h at 60 °C. After dehydrating, the samples were weighed using a Rice Lake TC balance. The fuel load per m2 was expressed in grams of dry matter and the fuel moisture percentage was calculated based on the difference between wet weight and dry weight [52].
Flame height during the experimental fire was measured by placing 4 metal equidistant bars in each plot, with a 5 m interval between each bar. The metal bars were painted with fire-resistant paint and were marked with 5 cm increments to accurately determine the height of the flames in each section. Additionally, photographs were taken during the experiment to confirm the flame height at each bar and calculate the average frame height. The images were carefully coded and analyzed using the UTHSCSA-Imagetool software (The University of Texas Health Science Center, San Antonio, TX, USA) [53].
To measure the soil temperature during the experimental burning, 4 thermocouples were placed per plot at a depth of 5 cm (using EasyLog data logger EL-USB-TC Thermocouples and K-TYPE PROBE 1M5 probe) [54]. The thermocouples were positioned along the centerline of each plot at the same level as the metal bars. The fire propagation speed was determined by using a stopwatch to measure the time it took for the flame front to travel between 2 bars (5 m or 20 m2) and was calculated using the formula provided by Aguirre [55].
Finally, once natural flame extinction occurred, 3 ash samples were collected randomly within each PSP. The samples were carefully stored in labeled plastic Petri dishes. At the laboratory, the ash samples were ground to powder to estimate fire severity using the Munsell chart color method, which has been widely adopted for such purposes [56,57,58].

2.4. Soil Physico-Chemical Properties Analyses

To evaluate the impact of fire on soil physico-chemical properties, samples were collected from burned and unburned plots. One set of samples was taken 24 h after the experimental burn (S1), and the other set was collected 365 days later (S2), following the recommendation of recent studies [59]. This choice was because the most significant changes in soil properties occur in the short and medium term, especially due to the climatic characteristics of the páramo that presents constant precipitation throughout the year and erosive processes can be enhanced in this ecosystem [32,34,48].
Sampling was performed by a single person using the same sampling cylinder, to ensure consistency in the measurements and avoid possible variations that could influence the results. In addition, during samplings S1 and S2, before collecting soil samples, the ash layer was removed according to Santín et al. [60]. In S1, the ash had an average thickness of 3 cm, while in S2, the layer had decreased considerably, presenting a few millimeters of thickness when it was present. For this purpose, the surface of the soil that had been in contact with the ash layer was carefully scraped with a razor to remove any carbonized residue, and then soil samples were collected.
Sampling consisted of collecting 3 soil samples (at the beginning, middle, and end of each PSP), resulting in nine soil subsamples per block, to determine soil bulk density (Bd) and soil water content. We also collected an additional set of nine soil samples at each PSP during collection times S1 and S2, for soil chemical characterization. These nine subsamples were mixed to obtain a composite sample for T1, T2, and the control plots, respectively, for each block. All soil samples were collected at a depth of 10 cm using a standard 283 cm3 metal cylinder (6 cm diameter and 10 cm height) [18]. Each of the samples obtained were placed in plastic bags and labeled appropriately for transport to the laboratory.
In the laboratory, Bd was immediately estimated [61]. Soil water content was determined by weight difference between wet and dry samples (105 °C for 48 h; [61]). Subsequently, an aliquot of the other set of samples intended for the determination of soil texture and chemical properties was dried at room temperature, sieved through a 2 mm mesh, and all visible root debris removed [62]. Soil texture was determined using the hydrometer method [63], while soil pH was measured with a pH meter applying the standard method [61]. Porosity was determined based on the assumption that the density of stones was approximately 2.65 g cm−3 [64]. Soil organic carbon (SOC) and soil organic matter (SOM) were determined by the method of Walkley and Black [65,66], for which the sample was placed in an oven at 125 °C for 45 min, after oxidation in a solution of K2Cr2O7/H2SO4. Soil ammonia nitrogen concentration (SAN; mg/kg) was measured by the colorimetric method [67], while soil available phosphorus (P; mg/kg), potassium (K; meq/100 g), calcium (Ca; meq/100 g), and magnesium (Mg; meq/100 g) were determined by the modified Olsen method [68].
Finally, the equation proposed by Walteros et al. [69] was used to determine soil organic carbon (SOC) content.
S O C = O C B d h
where SOC is the stock of soil organic carbon (tC ha−1 top 10 cm of soil), OC is the total organic carbon concentration (%), Bd is the bulk density (g cm−3), and h is the depth at which the sample was taken (cm).

2.5. Statistical Analysis

Statistical analyses were performed to evaluate fire behavior and its effect on soil physico-chemical properties. Differences in fire behavior descriptors (flame height, soil temperature, and propagation velocity) were evaluated using a one-way ANOVA (n = 3), with a significance level set at p < 0.05 and Tukey’s HSD post hoc test (p < 0.05). For this purpose, the mean value along different lengths (5 m, 10 m, 15 m, and 20 m) in the experimental plots was considered.
On the other hand, to determine the effect of fire on soil physico-chemical properties at different sampling times (S1, 24 h after burning; S2, one year after burning), a two-way analysis of variance with repeated measures (p < 0.05) was used. In addition, to verify the differences between treatments within each sampling time, one-factor ANOVA and Tukey’s HSD post hoc test (p < 0.05) were performed. Normality and homogeneity of variance assumptions of the data were corroborated by Shapiro–-Wilk and Levene tests (Shapiro–Wilk), respectively, before applying parametric tests.
Finally, to evaluate the integrated effects of fire treatments in soil physical properties, such as bulk density (Bd), porosity, and soil moisture, together with changes in chemical properties (SOM, SOC, pH, SAN, P, K, Ca, and Mg) between the two sample collection times (S1 and S2), a principal component analysis (PCA) was performed. The Pearson correlation coefficient between the soil physico-chemical properties and the scores of PC1 and PC2 was also estimated. All statistical analyses were performed with the statistical program PAST version 3 [70].

3. Results

3.1. Fire Behavior and Severity of the Burns

The San Lucas_UTPL automatic weather station recorded an average relative humidity of 72.6%, an average air temperature of 15.6 °C, a maximum solar radiation of 1135.1 (W/m2), and an average wind speed of 8.9 m/s (32.0 km/h). It is essential to mention that during the experiment, the wind direction concerning the longitudinal axis of each plot was always in favor of the downhill treatment, with directions mostly from the east to southeast. The absence of precipitation for several days before the experiment indicated good conditions for burning the PSPs with herbaceous páramo vegetation (Table 1).
According to the danger index of the Interagency Fire Use Module Field Guide [51], the atmospheric conditions during the experimental burn were classified as alerting, with low ignition probability (40%) and a relatively high wind speed of 8.9 m/s (Table 2). Moreover, the combustible moisture content of the vegetation (23.6%) was relatively high, with a wet biomass of 2213.7 g m−2 (22.1 t/ha) and a dry biomass of 1790.2 g m−2 (17.9 t/ha). It should be noted that soils in the AnPs are generally saturated [71,72], which leads to a classification of moderate fuel moisture content for the herbaceous páramo ecosystem.
There were no statistically significant differences between the flame heights of the plots burned uphill (T1) and downhill (T2) (p-value 0.08). The mean flame height for T1 was 59 cm and for T2 it was 71.6 cm, with the greatest differences observed in the center of the plots after 10 m (Figure 3a), where T1 reached a flame height of 51.7 cm and T2 reached a flame height of 75.5. Regarding soil temperature, no significant statistical difference was found between T1 and T2 (p-value 0.31) (Figure 3b). T1 reached a mean soil temperature of 26.4 °C, while T2 reached a mean soil temperature of 25.8 °C. The greatest differences were registered in the center of the burned plots after 10 m, where T1 had a soil temperature of 28.1 °C and T2 had a soil temperature of 25.5 °C. Soil temperature in the unburned control remained stable at around 13.3 °C. The mean propagation speed for T1 was 0.7 m/minute and for T2 it was 1.0 m/minute, with no significant statistical difference between them (p-value 0.06) (Figure 3c,d). However, the propagation speed within the first 5 m was slightly higher in the plots burned downhill (T2: 1.3 m/minute), compared to the T1 plots (T1: 0.7 m/minute), due to the predominate wind direction being from the southeast.
The ash samples from all burned plots had a very dark brown color according to the Munsell color code (=5 yr 2.5/2), which corresponds to low-severity fires [56,57].

3.2. Effects of the Experimental Fire on Soil Physico-Chemical Properties

Table 3 and Figure 4 summarize the effects of low-severity fires on the soil physico-chemical properties at a depth of 10 cm in the two sampling times: S1 (4 h after burning) and S2 (one post-fire year). Shortly after burning, Bd was not significantly different between the burned and unburned plots, including burn directions (uphill slope S1: 0.31 g/cm3; downhill slope S1: 0.30 g/cm3; and control S1: 0.33 g/cm3: p-value 0.90). One year after burning, Bd values increased in all plots, indicating a slight soil compaction, but with no significant differences between them (uphill S2: 0.50 g/cm3; downhill S2: 0.46 g/cm3; and control S2: 0.43 g/cm3: p-value 0.17) (Table 3; Figure 4a). This increase in Bd over time (S1 vs. S2) was statistically significant (p-value 0.02), which was also confirmed by porosity. In addition, porosity showed no significant differences between the burned and control plots in S1 (T1: 88.3%; T2: 88.8%; control: 87.7%: p-value 0.75), as well as in S2 (T1: 81.1%; T2: 82.5%; control: 83.87%: p-value 0.36). However, a significant decrease over time was observed (p-value 0.04; Table 3; Figure 4b). As with Bd and porosity, no significant differences were found for soil moisture in S1 between the burned and control plots (T1: 170.9%; T2: 169.7%; control: 141.2%: p-value 0.59), as well as in S2 (T1: 132.8%; T2: 154.0%; control: 126.7%: p-value 0.61). However, a sligtdecrease in soil moisture was observed over the different sampling dates (S1 vs. S2; p-value 0.05) (Table 3; Figure 4c). In general, as the statistical analyses indicated, fire had no significant effect on the physical soil properties, regardless of the treatment. However, the time factor showed significant changes (S1 vs. S2), in which an increase in Bd and a decrease in porosity, as well as in soil moisture, were observed.
The fire had no significant effects on SOM or SOC, and this was true for S1 and S2 (Table 3; Figure 4). Both tended to increase, although not significantly, during one year as the values in S2 tended to be higher than those of S1 (Figure 4d,e, respectively).
Although soil pH showed no significant differences between plots in S1 (p-value 0.38) and in S2 (p-value 0.20), the values were slightly higher in the burned plots. In addition, the pH values decreased with time in all plots, with the control plots always having the lowest values (Table 3, Figure 4f).
SAN showed no significant differences between burning treatments in S1 and in S2 (p-value 0.34 and 0.61, respectively) (Table 3). One year after burning (S2), SAN values were significantly higher than those of S1 (Figure 4g), and this was true independently of the fire. Values were as high as 144 mg/kg. No significant differences were found for soil magnesium among the plots in S1 and S2, but it increased significantly over time in all plots (p-value 0.02), regardless of the effect of fire. Mg values in S2 were almost four times higher than those measured in S1 in all plots (Table 3; Figure 4h).
No significant differences were found for P, either regarding the treatments (p-value 0.20) or the time factor (p-value 0.47). Only slight increases were observed in all plots one year after burning (Figure 4i). Soil exchangeable potassium (K, Figure 4j) showed an initial increase in the burned plots in S1, followed by a decrease one year later (S2). However, no significant differences in K were observed regarding the treatment (p-value 0.13) or the time factor (p-value 0.11), but a significant effect was found for the treatment–time interaction (p-value 0.01), as the values in the unburned plots did not change over time. Ca (Table 3; Figure 4k) showed an initial slight increase (S1) and a decrease one year later (S2) in the plots burned downhill (T2), while in the T1 and control plots, a slight increase in S2 compared to S1 was observed. However, the differences were not statistically significant either regarding the treatment (p-value 0.24) or the time factor (p-value 0.83) or the interaction between them.

3.3. Principal Component Analysis (PCA) in the Burned and Unburned Plots

Figure 5 shows the results of the principal component analysis (PCA) of the uphill burning, downhill burning, and control plots after one day (S1) and one year (S1) of the low-severity fire. The first two PCA components, PC1 and PC2, explain 57.3% and 23.4% of the total variance, respectively.
Bd, SOM, SOC, SAN, and Mg were all significantly and positively correlated with PC1 scores, and porosity and pH were negatively correlated with PC1 scores (Table 4). On the other hand, soil water content, K, and Ca were significantly correlated with PC2 scores (Table 4). PC1 differentiates soils according to their sampling time (Figure 5).
On the other hand, PC2 differentiates shortly after the fires between the burned and unburned soil, without differences between uphill and downhill burning.

4. Discussion

4.1. Fire Behavior and Severity during Experimental Burning

Determining the optimal time to carry out the burning of herbaceous vegetation in the páramo presents a challenge due to the frequent rainfall that characterizes this region [73]. Burning requires a rain-free period of at least one week [32], but unfortunately, accurate meteorological information is limited in high mountain areas [62]. Nevertheless, during the execution of the experimental burn, the atmospheric conditions were favorable (VdN phenomenon was recorded), which facilitated the successful ignition of herbaceous vegetation in the páramo (see Table 1) [32]. In this context, the Interagency Fire Use Module Field Guide Hazard Index [51] classified the atmospheric conditions on the day of the burn as alert, with moderate probability of ignition (40%) and a relatively high wind speed (see Table 2). It would have been advisable to have other meteorological conditions to carry out the experimental burns that were within a range of 15–19% relative humidity, 21–31 °C, and a wind speed of 5–11 m/s, based on which the Field Guide danger index would have classified conditions as alarm. Therefore, it is relevant to highlight that this guide was designed for U.S. ecosystems and its applicability may be limited in the unique conditions of the AnP ecosystems. Thus, it is essential to consider its use with caution in this specific region.
Compared to other studies in contrasting ecosystems, flame height in the herbaceous páramo vegetation was lower, mainly due to a lower vegetation height (40 cm) and lower fuel load (22.1 t/ha). For example, Rodríguez-Trejo et al. [74] reported higher flame heights in a 1 m artificial pasture in Chiapas (Mexico) with a fuel load of 6214 t/ha, and under drier climatic conditions. Anderson et al. [75] also found a positive correlation between fuel load and flame height in tropical savannas, indicating that higher fuel loads produce higher flames. Therefore, the flame height and propagation speed observed during the experimental burn are specific to the herbaceous páramo vegetation of the study area and may not be directly comparable to other ecosystems with different vegetation cover and fuel loads.
The severity of the experimental fire could be classified as low, as indicated by the very dark brown color of the ashes. According to Bodí et al. [57], this color indicates incomplete combustion of vegetation, including uncharred organic particles, caused by dehydration and oxidation of iron components. Fires of moderate to high severity are infrequent in AnPs due to the cool and humid weather conditions, which is why soils are generally saturated [76]. However, even low-severity fires could inhibit the regeneration of the AnP ecosystem when heavy rains occur after burning, which enhance erosive processes, causing the loss of the topsoil layer where organic matter and nutrients are concentrated. Therefore, careful management and monitoring of burning activities in AnPs is crucial to minimize the risk of soil degradation and promote ecosystem health. Nevertheless, low-severity fires can also provide benefits to ecosystems, as fire is a natural perturbation of most ecosystems on Earth, and many plant communities have adapted to it [77].

4.2. Low-Severity Fire Did Not Affect the Physico-Chemical Soil Properties

The low-severity of the burns did not manifest adverse effects on soil physical and chemical properties in the San Lucas AnP. These findings are consistent with previous research, such as that of Mehdi et al. [78], who reported that low-severity experimental burns did not cause significant alterations in soil properties compared to control sites. These results also align with studies conducted in various global ecosystems, such as those carried out by Agbeshie et al. [79], which corroborate that low-severity burns do not induce deteriorations in soil physical and chemical properties due to incomplete biomass combustion [80].
Low Bd and high porosity values in AnP soils are due to the high SOM concentration, which explains its high water storage capacity [81]. The increased Bd and decreased porosity one year after burning could be due to pore clogging by ash and clay particles, a phenomenon that was previously documented by Woods and Balfour [82]. The presence of these particles in the soil may reduce the water- and air-holding capacity of the pores, resulting in a decrease in porosity over time. This effect is corroborated by Larsen et al. [83], who showed that ash can clog soil macropores, leading to an increase in Bd, reducing soil infiltration rates, and thus increasing surface flow. This effect is especially evident in soils with high porosity, such as AnP soils, in the first 10 months after burning [82]. Later, however, soil porosity recovers, especially in regions with high amounts of previous precipitation [57]. Likewise, Bd values did not reach the level of soil compaction, which could inhibit germination [84], as values higher than 1.45 g cm−3 indicate this condition.
Consequently, low-severity burns in AnP do not negatively affect soil physical properties, as the destruction of soil aggregates by low-severity fires is not to be expected [85]. Furthermore, low-severity fires do not compromise the hydrological properties of AnP soils, even though fire usually reduces soil moisture by evaporation [86]. Although soil moisture decreased over time in the study area, this cannot be related to the experimental fire, as it depends on the fire temperature and soil heating, which were minimal in this AnP study, and due to climatic conditions, which show high annual and interannual variations in the Andes of southern Ecuador [87].
In addition, the low-severity experimental burns did not cause negative effects on the chemical properties of the studied AnP soils, as reported in recent research [79]. For example, the observed slight increase in SOM in all plots one year after burning (S2: upslope 33.7%; downslope 43.4%; and control 35.2%) could be due to the fact that the fencing of the study area led to increased litter deposition and facilitated natural succession, as the burned area was protected from grazing by domestic and wild animals (browsing of vegetation was prevented), which also contributed to SOM accumulation over time [14]. However, these values are relatively lower than those reported by Patiño et al. [88], who, analyzing 33 scientific papers in South American AnPs, determined an average of 43% SOM in the soil. Therefore, to better understand these results, it is recommended to continue monitoring soil organic matter content over time and to consider other factors such as soil erosion that could be influencing the observed changes.
On the other hand, the cold and humid climate of the AnP also favors SOM accumulation [2], and this process is enhanced by the formation of organometallic complexes, known as an association between organic matter and soil mineral particles (Al3+; Fe3+) [89]. This increase in SOM could also be favored by low fire severity, the slow decomposition of burned biomass, and ash addition [79]. However, it should be noted that a high SOM content has benefits for the edaphic ecology of this páramo since it is known that SOM allows the storage of a large amount of water, improving infiltration rates and hydraulic conductivity [90], which improves soil structure and leads to greater granularity and increased root development [91]. These edaphic conditions can especially favor the growth of native páramo species that are adapted to the adverse climatic conditions of these high-altitude areas, such as puya (Puya aequatorialis), páramo straw (Calamagrostis intermedia), páramo flower (Oritrophium crocifolium), naurapo (Myrteola nummularia), and joyapas (Disterigma rimbachii) [40,41].
Moreover, the climate in the studied AnP and the high water storage capacity of the soils lead to slow rates of the decomposition and accumulation of organic matter, which makes this ecosystem an important carbon sink [6]. This is corroborated by Medina and Mena [92], who reported that AnPs store six times more carbon than tropical forests, which is illustrated by the high SOC values calculated in this study at a soil depth of 10 cm (average S2: 93.7 tC ha−1; Table 3). These values are similar to the SOC values published by Santín and Vidal [60] for different soil types in a shrubby páramo area in Podocarpus National Park (PNP) in southern Ecuador (91.52 tC ha−1 for inceptisol soils and 68.37 tC ha−1 for entisol soils). However, the type of management and geographic location of the AnP also affect SOC. Quiroz Dahik et al. [93] only reported values of around 1.5 tC ha−1 in pastoral páramos (with extensive and intensive grazing) near the equator, while Cargua et al. [94] reported values of up to 277.81 tC ha−1 in páramos located on volcanoes due to frequent ash expulsion.
Another parameter that was affected by low-severity burns was pH. During the two sample collection times (S1 and S2), pH values in all PSPs ranged between 4.1 and 4.5, indicating suboptimal acidity for agricultural purposes, as the optimal soil pH for agriculture is between 5.0 and 8.0 [95]. According to Neina [96], low pH can negatively affect soil biogeochemistry and limit the availability of essential nutrients for plant growth. However, native, and endemic plant species in AnPs, which include native vascular plants, grasses, bryophytes, and lichens, may be adapted to these acidic soils. To elucidate this question, some researchers used acid phosphatase enzyme activity to determine soil quality/health in the tropical Andes [18]. For example, Turner [97] found that in tropical soils, acid phosphatase activity increases at a pH optimum of around 4.0, which is close to the values observed in this study. This may further explain the relatively high values of available phosphorus (P) in páramo soil (Figure 4i) and highlights the importance of soil pH in páramo ecology. Therefore, further investigation of the role of pH in this unique ecosystem is needed.
Similarly, low-severity burning did not have a negative effect on soil nutrients, as previous studies have shown [35]. In addition, it is important to note that low-severity burning did not result in nitrogen losses, as would sometimes be expected, but rather transformed it into a more plant-accessible form, thanks to suitable soil temperature and pH conditions [98].
As with SAN, magnesium (Mg; Figure 4h) also showed significant increases one year after burning (S2), even though it was expected that values would increase immediately after burning (S1) due to the combustion of organic materials, and decrease subsequently due to ash erosion, leaching, and plant uptake [99]. However, the slow decomposition process at these high altitudes and the increase in SOM led to an accumulation of Mg in the soils, which was not only observed in the burned plots, but also in the unburned control plots, probably due to the distribution of Mg over wider areas by surface and subsurface flows and by the protection of the study area with wooden posts and barbed wire.
To a lesser extent, available phosphorus (P; Figure 4i), which is the second most limiting plant nutrient after nitrogen, showed a slightincrease during S2. P becomes available through the transformation of organic P into organic matter during fire (i.e., mineralization), which increases P in the upper soil horizons [100]. The high P concentration observed during S2 (40.0 mg kg−1) could be attributed to the low soil pH and the presence of mycorrhizal fungi that release organic acids, which promote weathering and the decomposition of organic matter, leading to increased soil P [101].
The soil exchangeable potassium (K; Figure 4j) showed the expected pattern of an initial increase immediately after burning (S1), followed by a decrease one year later (S2) due to leaching and plant uptake [102]. Interestingly, calcium (Ca; Figure 4k) showed an increase in the burned plots in S1, but the decrease in S2 was only observed in the plots burned downhill (T2), while the values in the plots burned uphill (T1) stayed stable and the values in the control plot increased. This could confirms the nutrient transport and distribution by surface and subsurface flows in this precipitous terrain (see Figure 2).
In summary, this study demonstrates that low-severity burns do not have negative effects on soil physico-chemical properties in the herbaceous páramo ecosystems of southern Ecuador. In addition, one year after burning (S2), an increase in SOM, SOC, SAN, and Mg contents was observed. Likewise, the available phosphorus (P) content increased due to the low soil pH (see Figure 5, Table 4), which increases acid phosphatase activity [97]. These results are consistent with recent studies, which indicate that low-severity fires increase nutrient availability and favor post-fire natural regeneration as well as plant community growth [35,100].
However, our study cannot verify any increase or decrease in nutrient availability that can be confidently attributed to fire. Therefore, further research on burning practices in AnP, including fire severity and frequency, is needed to better understand their impact on post-fire soil and plant development [103].
Looking to the future, it is important to control and reduce anthropogenic pressures on AnPs due to their capacities as carbon sinks and for water storage. Conversely, an intensive conversion to agricultural land and cattle ranching should be avoided, as well as avoiding the establishment of pine plantations (Pinus patula and Pinus radiata) and large-scale mining as occurs in other AnPs in Ecuador [6]. According to Avellaneda-Torres et al. [104], all of these activities have negative impacts on water supply and nutrient availability, which subsequently alter natural succession and intensify soil erosion processes, which could lead to the loss of the biological structure of the AnP when the topsoil is completely lost [22]. It is therefore crucial to adopt sustainable management practices that consider the fragility of the ecosystem and the services it provides to humans and wildlife (water regulation, carbon sink), in addition to biodiversity conservation. Such practices should give priority to the conservation and restoration of AnPs, including alternative income for the livelihoods of local people, such as sustainable tourism in this unique ecosystem.

5. Conclusions

This study analyzed the impact of low-severity fire on soil physico-chemical properties in a herbaceous páramo in southern Ecuador. The results indicate that the optimal period for burning is during the Veranillo del Niño (VdN) phenomenon, which is characterized by favorable climatic conditions involving a decrease in relative humidity, increased temperature, and solar radiation, and, especially, the absence of precipitation. The study concluded that low-severity burns have no significant effects on soil physico-chemical properties, in terms of soil organic matter and nutrient availability, even one year after burning. However, long-term monitoring is necessary to analyze soil erosion processes and succession dynamics, as well as nutrient availability due to leaching over time. Conserving this natural ecosystem and protecting it from anthropogenic activities, such as intensive agriculture and cattle ranching, is essential to ensure ecosystem services such as water supply for the local population and the reduction in food risks in the valleys. In addition, AnPs are important carbon sinks that absorb significant amounts of the anthropogenic greenhouse gas emissions [105], which is essential for buffering the effects of ongoing climate change. Finally, it is crucial to implement comprehensive fire management policies to ensure sustainable management of AnPs.

Author Contributions

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

Funding

This research was carried out with the technical support of the trilateral cooperation project “Strengthening of Technical and Institutional Capacities in Ecuador for Integrated Fire Management in Conservation Areas” implemented by the German Cooperation for Sustainable Development through the Deutsche Gesellschaft für Internationale Zusammenarbeit (GIZ) GmbH, the Amazonía Sin Fuego Ecuador Program of the Ministry of Environment and Water of Ecuador, by the Brazilian Cooperation Agency (ABC), the National Center for Prevention and Combat of WildFires—Prevfogo of the Brazilian Institute of Environment and Renewable Natural Resources (Ibama of Brazil), and the National System of Conservation Areas (SINAC) of Costa Rica. The funding code is PROY_CONS_CCBIO_2020_2751.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

Our thanks go to the Universidad Técnica Particular de Loja and GIZ_Ecuador for funding this research (PROY_CONS_CCBIO_2020_2751).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Madriñán, S.; Cortés, A.J.; Richardson, J.E. Páramo is the world’s fastest evolving and coolest biodiversity hotspot. Front. Genet. 2013, 4, 192. [Google Scholar] [CrossRef] [PubMed]
  2. Castañeda-Martín, A.E.; Montes-Pulido, C.R. Carbono almacenado en páramo andino. Entramado 2017, 13, 210–221. [Google Scholar] [CrossRef]
  3. Jiménez-Rivillas, C.; García, J.J.; Quijano-Abril, M.A.; Daza, J.M.; Morrone, J.J. A new biogeographical regionalisation of the Páramo biogeographic province. Aust. Syst. Bot. 2018, 31, 296–310. [Google Scholar] [CrossRef]
  4. Hofstede, R. Los Páramos en el Mundo: Su Diversidad y sus Habitantes; Hofstede, R., Segarra, P., Vásconez, P.M., Eds.; Global Peatland Initiative/NC-IUCN/EcoCiencia: Quito, Ecuador, 2003. [Google Scholar]
  5. Serrano-Giné, D.; Galárraga-Sánchez, R. El páramo andino: Características territoriales y estado ambiental. Aportes interdisciplinarios para su conocimiento. Estud. Geogr. 2015, 76, 369–393. [Google Scholar] [CrossRef]
  6. Chuncho, C.; Chuncho, G. Páramos del Ecuador, importancia y afectaciones: Una revisión. Bosques Latid. Cero. 2019, 9, 71–83. [Google Scholar]
  7. Acevedo Amaya, O.J. Identificación de las Principales Variables de los Suelos Asociadas a Condiciones de Humedad: Discriminación Taxonómica de Suelos Hidromórficos; Instituto de Investigación de Recursos Biológicos Alexander von Humboldt: Bogotá, Colombia, 2014. [Google Scholar]
  8. Daza Torres, M.C.; Hernández Florez, F.; Alba Triana, F. Efecto del Uso del Suelo en la Capacidad de Almacenamiento Hídrico en el Páramo de Sumapaz-Colombia. Rev. Fac. Nal. Agrom. Medellín 2014, 67, 7189–7200. [Google Scholar] [CrossRef]
  9. Buytaert, W.; Célleri, R.; De Bièvre, B.; Cisneros, F.; Wyseure, G.; Deckers, J.; Hofstede, R. Human impact on the hydrology of the Andean páramos. Earth-Sci. Rev. 2006, 79, 53–72. [Google Scholar] [CrossRef]
  10. Llambí, L.D.; Soto, A.; Celleri, R.; De Bievre, B.; Ochoa, B.; Borja, P. Páramos Andinos: Ecología, Hidrología y Suelos de Páramos. Los Suelos del Páramo. 2012. Available online: https://biblio.flacsoandes.edu.ec/libros/digital/56475.pdf (accessed on 15 May 2022).
  11. Van der Hammen, T.; Otero-García, J.; Morales, M.; Otero, J. Los páramos: Archipiélagos terrestres en el norte de los Andes. In Atlas de Paramos de Colombia; Instituto de Investigación de Recursos Biológicos Alexander von Humboldt: Bogotá, Colombia, 2007. [Google Scholar]
  12. Iñiguez, V.; Morales, O.; Cisneros, F.; Bauwens, W.; Wyseure, G. Analysis of the drought recovery of Andosols on southern Ecuadorian Andean páramos. Hydrol. Earth Syst. Sci. 2016, 20, 2421–2435. [Google Scholar] [CrossRef]
  13. Pinos-Morocho, D.; Morales-Matute, O.; Durán-López, M.E. Suelos de páramo: Análisis de percepciones de los servicios ecosistémicos y valoración económica del contenido de carbono en la sierra sureste del Ecuador. Rev. Cienc. Ambient. 2021, 55, 151–173. [Google Scholar] [CrossRef]
  14. Shoji, S.; Nanzyo, M.; Dahlgren, R.A. Volcanic Ash Soils: Genesis, Properties, and Utilization; Elsevier: Amsterdam, The Netherlands, 1994. [Google Scholar]
  15. Buytaert, W. The Properties of the Soils of the South Ecuadorian Páramo and the Impact of Land Use Changes on Their Hydrology. Ph.D. Thesis, Katholieke Universiteit Leuven, Leuven, Belgium, 2004. [Google Scholar]
  16. Lazo, P.X.; Mosquera, G.M.; McDonnell, J.J.; Crespo, P. The role of vegetation, soils, and precipitation on water storage and hydrological services in Andean Páramo catchments. J. Hydrol. 2019, 572, 805–819. [Google Scholar] [CrossRef]
  17. Cortés-Duque, J.; Sarmiento Pinzón, C.E. Visión Socioecosistémica de los Páramos y la alta Montaña Colombiana. Memorias del Proceso de Definición de Criterios para la Delimitación de Páramos; Instituto de Investigación de Recursos Biológicos Alexander von Humboldt: Bogotá, Colombia, 2013. [Google Scholar]
  18. Carrión-Paladines, V.; Benítez, Á.; García-Ruíz, R. Conversion of Andean montane forest to exotic forest plantation modifies soil physicochemical properties in the buffer zone of Ecuador’s Podocarpus National Park. For. Ecosyst. 2022, 9, 100076. [Google Scholar] [CrossRef]
  19. Diaz-Acevedo, C.J.; Romero-Alarcon, L.V.; Miranda-Esquivel, D.R. Páramos neotropicales como unidades biogeográficas. Rev. Biol. Trop. 2020, 68, 508–516. [Google Scholar] [CrossRef]
  20. DeMay, S.M.; Roon, D.A.; Rachlow, J.L.; Cisneros, R. Selective foraging on bromeliads by Andean bears in the Ecuadorian páramo. Ursus 2014, 25, 139–147. [Google Scholar] [CrossRef]
  21. Otero, J.D.; Figueroa, A.; Muñoz, F.A.; Peña, M.R. Loss of soil and nutrients by surface runoff in two agro-ecosystems within an Andean paramo area. Ecol. Eng. 2011, 37, 2035–2043. [Google Scholar] [CrossRef]
  22. Thompson, J.B.; Zurita-Arthos, L.; Müller, F.; Chimbolema, S.; Suárez, E. Land use change in the Ecuadorian páramo: The impact of expanding agriculture on soil carbon storage. Arct. Antarct. Alp. Res. 2021, 53, 48–59. [Google Scholar] [CrossRef]
  23. Carrillo-Rojas, G.; Silva, B.; Rollenbeck, R.; Célleri, R.; Bendix, J. The breathing of the Andean highlands: Net ecosystem exchange and evapotranspiration over the páramo of southern Ecuador. Agric. For. Meteorol. 2019, 265, 30–47. [Google Scholar] [CrossRef]
  24. Sierra, R.; Cerón, C.; Palacios, W.; Valencia, R. Criterios para la Clasificación de la Vegetación del Ecuador; Sierra, R., Ed.; Propuesta preliminar de un sistema de clasificación de vegetación para el Ecuador continental; Proyecto INEFAN/GEF-BIRF y EcoCiencia: Quito, Ecuador, 1999. [Google Scholar]
  25. Hofstede, R. El páramo como espacio para la fijación de carbono atmosférico. In El Páramo Como Espacio de Mitigación de Carbono Atmosférico; Serie Páramo 1; GTP/Abya-Yala: Quito, Ecuador, 1999. [Google Scholar]
  26. Messina, J.P.; Cochrane, M.A. The forests are bleeding: How land use change is creating a new fire regime in the Ecuadorian Amazon. J. Latin Am. Geogr. 2007, 6, 85–100. [Google Scholar] [CrossRef]
  27. Schritt, H.; Beusch, C.; Ríos Guayasamín, P.; Kaupenjohann, M. Transformation of traditional shifting cultivation into permanent cropping systems: A case study in Sarayaku, Ecuador. Ecol. Soc. 2020, 25, 10. [Google Scholar] [CrossRef]
  28. Cochrane, M.A.; Horn, S.P.; Kappelle, M. Fire in the páramo ecosystems of Central and South America. In Tropical Fire Ecology: Climate Change, Land Use, and Ecosystem Dynamics; Springer: Cham, Switzerland, 2009; pp. 505–539. [Google Scholar]
  29. Farley, K.A.; Kelly, E.F. Effects of afforestation of a páramo grassland on soil nutrient status. For. Ecol. Manag. 2004, 195, 281–290. [Google Scholar] [CrossRef]
  30. Buytaert, W.; Iniguez, V.; De Bievre, B. The effects of afforestation and cultivation on water yield in the Andean páramo. For. Ecol. Manag. 2007, 251, 22–30. [Google Scholar] [CrossRef]
  31. Armenteras, D.; González, T.M.; Vargas Ríos, O.; Meza Elizalde, M.C.; Oliveras, I. Incendios en ecosistemas del norte de Suramérica: Avances en la ecología del fuego tropical en Colombia, Ecuador y Perú. Caldasia 2020, 42, 1–16. [Google Scholar] [CrossRef]
  32. Díaz, S.C.; Quezada, L.C.; Álvarez, L.J.; Loján-Córdova, J.; Carrión-Paladines, V. Indigenous use of fire in the paramo ecosystem of southern Ecuador: A case study using remote sensing methods and ancestral knowledge of the Kichwa Saraguro people. Fire Ecol. 2023, 19, 5. [Google Scholar] [CrossRef] [PubMed]
  33. Minervini, M.G.; Morrás, H.J.; Taboada, M.A. Efectos del fuego en la matriz del suelo: Consecuencias sobre las propiedades físicas y mineralógicas. Ecol. Aust. 2018, 28, 12–27. [Google Scholar] [CrossRef]
  34. Carrión-Paladines, V.; Hinojosa, M.B.; Jiménez Álvarez, L.; Reyes-Bueno, F.; Correa Quezada, L.; García-Ruiz, R. Effects of the Severity of Wildfires on Some Physical-Chemical Soil Properties in a Humid Montane Scrublands Ecosystem in Southern Ecuador. Fire 2022, 5, 66. [Google Scholar] [CrossRef]
  35. Chandra, K.K.; Bhardwaj, A.K. Incidence of forest fire in India and its effect on terrestrial ecosystem dynamics, nutrient, and microbial status of soil. Int. J. Res. Agric. For. 2015, 5, 69–78. [Google Scholar]
  36. Mataix-Solera, J.; Cerdà, A.; Arcenegui, V.; Jordán, A.; Zavala, L. Fire effects on soils aggregation: A review. Earth-Sci. Rev. 2011, 109, 44–60. [Google Scholar] [CrossRef]
  37. Podwojewski, P.; Poulenard, J.; Toulkeridis, T.; Gräfe, M. Polygenic soils in the southern central Ecuadorian highlands as the result of long-lasting pedogenesis, geodynamic processes and climate change. J. S. Am. Earth Sci. 2022, 120, 104096. [Google Scholar] [CrossRef]
  38. Calispa, M.; Vásconez, F.; Santamaría, S.; Samaniego, P. Los Suelos de los Páramos del Ecuador En: Los Páramos del Ecuador: Pasado, Presente y Futuro; Hofstede, R., Suárez, E., Mena-Vásconez, P., Eds.; Universidad San Francisco: Quito, Ecuador, 2023. [Google Scholar] [CrossRef]
  39. Kappelle, M.; Horn, S.P. The Páramo ecosystem of Costa Rica’s highlands. In Costa Rican Ecosystems; The University of Chicago Press: Chicago, IL, USA, 2016; pp. 492–523. [Google Scholar]
  40. Galeas, R.; Guevara, J.E.; Medina-Torres, B.; Chinchero, M.Á.; Herrera, X. Sistema de Clasificación de los Ecosistemas del Ecuador Continental; Ministerio del Ambiente del Ecuador: Quito, Ecuador, 2013. [Google Scholar]
  41. San Lucas GAD Parroquial; LOJA, Cantón. Actualización del Plan de Desarrollo y Ordenamiento Territorial 2015–2019; Plan de Desarrollo y Ordenamiento Territorial: Loja, Ecuador, 2015. [Google Scholar]
  42. Coronel-Alulima, T.N. Los sistemas agroecológicos de la parroquia San Lucas (Loja). Prácticas resilientes ante el cambio climático. Let. Verdes Rev. Latinoam. Estud. Socioambientales 2019, 26, 191–212. [Google Scholar] [CrossRef]
  43. Hueso-González, P.; Martínez-Murillo, J.F.; Ruiz-Sinoga, J.D. Prescribed fire impacts on soil properties, overland flow and sediment transport in a Mediterranean forest: A 5 year study. Sci. Total Environ. 2018, 636, 1480–1489. [Google Scholar] [CrossRef]
  44. Nelson, R.M., Jr. An effective wind speed for models of fire spread. Int. J. Wildland Fire. 2002, 11, 153–161. [Google Scholar] [CrossRef]
  45. Gillon, D.; Houssard, C.; Valette, J.C.; Rigolot, E. Nitrogen and phosphorus cycling following prescribed burning in natural and managed Aleppo pine forests. Can. J. For. Res. 1999, 29, 1237–1247. [Google Scholar] [CrossRef]
  46. Calderón, R.D.J.B.; Flores, R.E.E. Análisis de la problemática de los incendios forestales y sus implicaciones en la pérdida de cobertura forestal en Cintalapa, Chiapas. Rev. Geogr. Agríc. 2020, 65, 63–83. [Google Scholar]
  47. Waldrop, T.A.; Goodrick, S.L. Introduction to Prescribed Fires in Southern Ecosystems. Science Update SRS-054; US Department of Agriculture Forest Service, Southern Research Station: Asheville, NC, USA, 2012; Volume 80, p. 54.
  48. Bendix, J.; Rollenbeck, R.; Fabian, P.; Emck, P.; Richter, M.; Beck, E. Climate variability. In Gradients in a Tropical Mountain Ecosystem of Ecuador; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  49. Esquivel-Hernández, G.; Sánchez-Murillo, R.; Vargas-Salazar, E. Chirripó hydrological research site: Advancing stable isotope hydrology in the Central American Páramo. Hydrol. Process. 2021, 35, e14181. [Google Scholar] [CrossRef]
  50. Geron, C.; Hays, M. Air emissions from organic soil burning on the coastal plain of North Carolina. Atmos. Environ. 2013, 64, 192–199. [Google Scholar] [CrossRef]
  51. National Wildfire Coordinating Group. Interagency Fire Use Module Field Guide. NWCG Guide for Wildland Fire Modules, PMS 431, Boise, Idaho, USA. 2014. Available online: https://training.nwcg.gov/pre-courses/s490/s490pcw/Interagency%20Fire%20Use%20Module.pdf (accessed on 22 October 2022).
  52. Manals-Cutiño, E.M.; Penedo-Medina, M.; Salas-Tort, D. Caracterización del bagazo de caña como biomasa vegetal. Tecnol. Quím. 2015, 35, 244–255. [Google Scholar]
  53. Cardoso, A.Á.; Pereira, F.J.; Pereira, M.P.; Corrêa, F.F.; de Castro, E.M.; Santos, B.R. Anatomy of stems, leaves, roots, and the embryo of Garcinia brasiliensis Mart.-Clusiaceae. Rev. Ciênc. Agrár. Amaz. J. Agric. Environ. Sci. 2013, 56, 23–29. [Google Scholar] [CrossRef]
  54. Rocha, D.R.; Melquiades, F.L.; Thomaz, E.L. Modeling the soil burn effect for temperature prediction by energy dispersive X ray Fluorescence in an haplic cambisol soil. Appl. Radiat. Isot. 2019, 150, 26–30. [Google Scholar] [CrossRef] [PubMed]
  55. Aguirre Briones, F. Manual de Prevención de Riesgos en Incendios Forestales; TRAGSA: Madrid, Spain, 1996. [Google Scholar]
  56. Pereira, P.; Jordán, A.; Cerdà, A.; Martin, D. Editorial: The role of ash in fire-affected ecosystems. Catena 2014, 135, 337–339. [Google Scholar] [CrossRef]
  57. Bodí, M.; Martin, D.A.; Santín, C.; Balfour, V.; Doerr, S.H.; Pereira, P.; Cerdà, A.; Mataix Solera, J. Wildland fire ash: Production, composition, and eco-hydro-geomorphic effects. Earth Sci. Rev. 2014, 130, 103–127. [Google Scholar] [CrossRef]
  58. Alexakis, D.E. Contaminated land by wildfire effect on ultramafic soil and associated human health and ecological risk. Land 2020, 9, 409. [Google Scholar] [CrossRef]
  59. Valdes, L. Efectos del fuego en algunas características de suelos de pinares, Macurije, Pinar del Río, Cuba. Rev. Latinoam. Recur. Nat. 2016, 12, 60–65. [Google Scholar]
  60. Santín, A.; Vidal, E. Generación de una Línea Base de los Reservorios de Carbono de los Páramos del PNP y Evaluación de su Aplicación como Mecanismos de Mitigación al Cambio Climático. Ph.D. Thesis, Tesis de Ingeniería en Manejo y Conservación del Medio Ambiente, Loja, Ecuador, 2012. [Google Scholar]
  61. Sandoval, M.; Fernández, J.; Seguel, O.; Becerra, J.; Salazar, D. Métodos de Análisis Físicos de Suelos. Sociedad Chilena de la Ciencia del Suelo. Ph.D. Thesis, Universidad de Concepción, Concepción, Chile, 2011; pp. 1–75. [Google Scholar]
  62. Fries, A.; Silva, K.; Pucha-Cofrep, F.; Oñate-Valdivieso, F.; Ochoa-Cueva, P. Water balance and soil moisture deficit of different vegetation units under semiarid conditions in the andes of southern Ecuador. Climate 2020, 8, 30. [Google Scholar] [CrossRef]
  63. Black, C.A.; Evans, D.D.; White, J.L.; Ensminger, L.E.; Clark, F.E. Methods of Soil Analysis. Part 2. Chemical and Microbiological Properties; John Wiley & Sons: Hoboken, NJ, USA, 1965. [Google Scholar]
  64. Mingming, G.U.O.; Wenlong, W.; Hongliang, K.; Bo, Y. Changes in soil properties and erodibility of gully heads induced by vegetation restoration on the Loess Plateau, China. J. Arid Land 2018, 10, 712–725. [Google Scholar] [CrossRef]
  65. Page, A.L.; Miller, R.H.; Keeney, D.R. Methods of Soil Analysis, Part II; American Society of Agronomy: Madison, WI, USA, 1982. [Google Scholar]
  66. Vela Correa, G.; López Blanco, J.; Rodríguez Gamiño, M.D.L. Niveles de carbono orgánico total en el Suelo de Conservación del Distrito Federal, centro de México. Investig. Geogr. 2012, 77, 18–30. [Google Scholar]
  67. Keeney, D.R.; Nelson, D.W. Nitrogen-Inorganic Forms. In Methods of Soil Analysis, Agronomy Monograph 9, Part 2, 2nd ed.; Page, A.L., Ed.; ASA, SSSA: Madison, WI, USA, 1982; pp. 643–698. [Google Scholar] [CrossRef]
  68. Chowdhury, S.; Manjón-Cabeza, J.; Ibáñez, M.; Mestre, C.; Broncano, M.J.; Mosquera-Losada, M.R.; Plaixats, J.; Sebastià, M.T. Responses in Soil Carbon and Nitrogen Fractionation after Prescribed Burning in the Montseny Biosphere Reserve (NE Iberian Peninsula). Sustainability 2022, 14, 4232. [Google Scholar] [CrossRef]
  69. Walteros-Torres, I.; Palacios-Pacheco, S.; Cely, G.E.; Serrano, P.A.; Moreno-Pérez, D.F. Influencia del cambio de uso del suelo sobre las reservas de carbono orgánico en el Parque Natural Regional Cortadera, Boyacá (Colombia). Rev. UDCA Actual. Divulg. Cient. 2022, 25, 1–9. [Google Scholar]
  70. Hammer, Ø.; Harper, D.A.; Ryan, P.D. PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 2001, 4, 9. [Google Scholar]
  71. Yuste, J.C.; Hereş, A.M.; Ojeda, G.; Paz, A.; Pizano, C.; García-Angulo, D.; Lasso, E. Soil heterotrophic CO2 emissions from tropical high-elevation ecosystems (Páramos) and their sensitivity to temperature and moisture fluctuations. Soil Biol. Biochem. 2017, 110, 8–11. [Google Scholar] [CrossRef]
  72. Oleszczuk, R.; Truba, M. The analysis of some physical properties of drained peat-moorsh soil layers. Annals of Warsaw University of Life Sciences-SGGW. Land Reclam. 2013, 45, 41–48. [Google Scholar]
  73. Luteyn, J.L.; Churchill, S.P.; Griffin, D.; Gradstein, S.R.; Sipman, H.J.; Mauricio, R.; Gavilanes, A. Páramos: A Checklist of Plant Diversity, Geographical Distribution, and Botanical Literature; New York Botanical Garden Press: New York, NY, USA, 1999; Volume 278. [Google Scholar]
  74. Rodríguez-Trejo, D.A.; Martínez-Muñoz, P.; Pulido-Luna, J.A.; Martínez-Lara, P.J.; Cruz-López, J.D. Combustibles, comportamiento del fuego y emisiones en un pastizal y una sabana artificiales en Chiapas, México. Rev. Biol. Trop. 2020, 68, 641–654. [Google Scholar] [CrossRef]
  75. Anderson, A.N.; Cook, G.D.; Williams, R.J. Fire in Tropical Savannas; Springer: New York, NY, USA, 2003. [Google Scholar]
  76. García, V.J.; Márquez, C.O.; Rodríguez, M.V.; Orozco, J.J.; Aguilar, C.D.; Ríos, A.C. Páramo Ecosystems in Ecuador’s Southern Region: Conservation State and Restoration. Agronomy 2020, 10, 1922. [Google Scholar] [CrossRef]
  77. Richter, C.; Rejmánek, M.; Miller, J.E.; Welch, K.R.; Weeks, J.; Safford, H. The species diversity× fire severity relationship is hump-shaped in semiarid yellow pine and mixed conifer forests. Ecosphere 2019, 10, e02882. [Google Scholar] [CrossRef]
  78. Mehdi, H.; Ali, S.; Ali, M.; Mostafa, A. Effects of different fire severity levels on soil chemical and physical properties in Zagros forests of western Iran. Folia For. Pol. 2012, 54, 241–250. [Google Scholar]
  79. 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]
  80. 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 afect the biological and chemical soil properties of Pinus halepensis Mill. stands in the short and mid-terms after wildfre. J. Environ. Manag. 2019, 235, 250–256. [Google Scholar] [CrossRef] [PubMed]
  81. Minasny, B.; McBratney, A.B. Limited effect of organic matter on soil available water capacity. Eur. J. Soil Sci. 2018, 69, 39–47. [Google Scholar] [CrossRef]
  82. Woods, S.W.; Balfour, V.N. The effects of soil texture and ash thickness on the post-fire hydrological response from ash-covered soils. J. Hydrol. 2010, 393, 274–286. [Google Scholar] [CrossRef]
  83. Larsen, I.; MacDonald, L.H.; Brown, E.; Rough, D.; Welsh, M.J.; Pietraszek, J.H.; Libohova, Z.; Benavides-Solorio, J.D.; Schaffrath, K. Causes of post-fire runoff and erosion: Water repellency, cover, or soil sealing? Soil Sci. Soc. Am. J. 2009, 73, 1393–1407. [Google Scholar] [CrossRef]
  84. Laureda, D.A.; Botta, G.F.; Tolón Becerra, A.; Rosatto, H.G. Compactación del suelo inducida por la maquinaria en campos de polo en Argentina. Rev. Fac. Cienc. Agrar. 2016, 48, 79–99. [Google Scholar]
  85. Lucas-Borja, M.E.; Ortega, R.; Miralles, I.; Plaza-Álvarez, P.A.; González-Romero, J.; Peña-Molina, J.; Moya, D.; Zema, D.A.; 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]
  86. Erkovan, Ş.; Koc, A.; Güllap, M.K.; Erkovan, H.I.; Bilen, S. The effect of fire on the vegetation and soil properties of ungrazed shortgrasssteppe rangeland of the Eastern Anatolia region of Turkey. Turk. J. Agric. 2016, 40, 290–299. [Google Scholar] [CrossRef]
  87. Rollenbeck, R.; Bendix, J. Rainfall distribution in the Andes of southern Ecuador derived from blending weather radar data and meteorological field observations. Atmos. Res. 2011, 99, 277–289. [Google Scholar] [CrossRef]
  88. Patiño, S.; Hernández, Y.; Plata, C.; Domínguez, I.; Daza, M.; Oviedo-Ocaña, R.; Buytaer, W.; Ochoa-Tocachi, B.F. Influence of land use on hydro-physical soil properties of Andean páramos and its effect on streamflow buffering. Catena 2021, 202, 105227. [Google Scholar] [CrossRef]
  89. Wagai, R.; Kajiura, M.; Asano, M.; Hiradate, S. Nature of soil organo-mineral assemblage examined by sequential density fractionation with and without sonication: Is allophanic soil different? Geoderma 2015, 241, 295–305. [Google Scholar] [CrossRef]
  90. Harden, C.P. Human impacts on headwater fluvial systems in the northern and central Andes. Geomorphology 2006, 79, 249–263. [Google Scholar] [CrossRef]
  91. Osorio, A.; Bahamon, C. Dinámica de la Humedad del Suelo en Bosques Altoandinos en el Páramo de Guerrero, Cundinamarca-Colombia. Ph.D. Thesis, Universidad Nacional de Colombia sede Medellín, Medellín, Colombia, 2008. [Google Scholar]
  92. Medina, G.; Mena, P.V. Páramo como Espacio de Mitigación de Carbono Atmosférico; GTP: Quito, Ecuador, 1999. [Google Scholar]
  93. Quiroz Dahik, C.; Crespo, P.; Stimm, B.; Mosandl, R.; Cueva, J.; Hildebrandt, P.; Weber, M. Impacts of pine plantations on carbon stocks of páramo sites in southern Ecuador. Carbon Balance Manag. 2021, 16, 1–15. [Google Scholar] [CrossRef] [PubMed]
  94. Cargua, F.E.; Rodríguez, M.V.; Recalde, C.G.; Vinueza, L.M. Cuantificación del contenido de carbono en una plantación de pino insigne (Pinus radiata) y en estrato de páramo de ozogoche bajo, Parque Nacional Sangay, Ecuador. Inf. Tecnol. 2014, 25, 83–92. [Google Scholar] [CrossRef]
  95. McCauley, A.; Jones, C.; Jacobsen, J. Soil pH and organic matter. Nut. Manag. Modul. 2009, 8, 1–12. [Google Scholar]
  96. Neina, D. The role of soil pH in plant nutrition and soil remediation. Cienc. Suelo Apl. Ambient. 2019, 2019, 1–9. [Google Scholar] [CrossRef]
  97. Turner, B.L. Variation in pH optima of hydrolytic enzyme activities in tropical rain forest soils. Appl. Environ. Microbiol. 2010, 76, 6485–6493. [Google Scholar] [CrossRef]
  98. Rodriguez-Cardona, B.M.; Coble, A.A.; Wymore, A.S.; Kolosov, R.; Podgorski, D.C.; Zito, P.; Spencer, R.G.M.; Prokushkin, A.S.; McDowell, W.H. Wildfires lead to decreased carbon and increased nitrogen concentrations in upland arctic streams. Sci. Rep. 2020, 10, 8722. [Google Scholar] [CrossRef] [PubMed]
  99. Lucas-Borja, M.E.; Plaza-Alvarez, P.A.; Xu, X.; Carra, B.G.; Zema, D.A. Exploring the factors influencing the hydrological response of soil after low and high-severity fires with post-fire mulching in Mediterranean forests. Int. Soil Water Conserv. Res. 2023, 11, 169–182. [Google Scholar] [CrossRef]
  100. Liu, J.; Qiu, L.; Wang, X.; Wei, X.; Gao, H.; Zhang, Y.; Cheng, J. Effects of wildfire and topography on soil nutrients in a semiarid restored grassland. Plant Soil 2018, 428, 123–136. [Google Scholar] [CrossRef]
  101. Mendoza, L.G.; Vera, V.; Giler, J.M.; Simbaña, K. Características fisicoquímicas de suelos de uso agrícola y forestal. Caso: San Pablo de Tarugo, Chone–Ecuador. Rev. Cient. Cienc. Nat. Ambient. 2022, 16, 334–341. [Google Scholar] [CrossRef]
  102. Akburak, S.; Son, Y.; Makineci, E.; Çakir, M. Impacts of low-intensity prescribed fire on microbial and chemical soil properties in a Quercus frainetto forest. J. For. Res. 2018, 29, 687–696. [Google Scholar] [CrossRef]
  103. Bowman, D.M.; Perry, G.L.; Higgins, S.L.; Johnson, C.N.; Fuhlendorfand, S.D.; Murphy, B.P. Pyrodiversity is the coupling of biodiversity and fire regimes in food webs. Philos. Trans. R. Soc. B Biol. Sci. 2016, 371, 20150169. [Google Scholar] [CrossRef] [PubMed]
  104. Avellaneda-Torres, L.M.; Sicard, T.E.L.; Rojas, E.T. Impact of potato cultivation and cattle farming on physicochemical parameters and enzymatic activities of Neotropical high Andean Páramo ecosystem soils. Sci. Total Environ. 2018, 631, 1600–1610. [Google Scholar] [CrossRef]
  105. IPCC. Summary for Policymakers. In Climate Change 2022: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Pörtner, H.-O., Roberts, D.C., Tignor, M., Poloczanska, E.S., Mintenbeck, K., Alegría, A., Craig, M., Langsdorf, S., Löschke, S., Möller, V., et al., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2022; pp. 3–33. [Google Scholar] [CrossRef]
Figure 1. Digital elevation model (a,b) and principal land uses in the San Lucas parish in southern Ecuador. The blue rhombus indicates the experimental burning site.
Figure 1. Digital elevation model (a,b) and principal land uses in the San Lucas parish in southern Ecuador. The blue rhombus indicates the experimental burning site.
Fire 06 00447 g001
Figure 2. Experimental burning design in the study area: (a) photograph of the experimental site after burning, and (b) the burning pattern design, including the burning direction (upslope or downslope) and ignition starting line (Ig).
Figure 2. Experimental burning design in the study area: (a) photograph of the experimental site after burning, and (b) the burning pattern design, including the burning direction (upslope or downslope) and ignition starting line (Ig).
Fire 06 00447 g002
Figure 3. Evaluation of fire behavior during experimental burning through an ANOVA (n = 3), with a significance level set at p < 0.05. (a) Average flame height (cm); (b) average soil temperature reached at a depth of 5 cm; (c) average fire propagation speed (m/minute). Different letters (a, b) indicate significant difference among burning treatments (p ≤ 0.05%, HSD Tukey).
Figure 3. Evaluation of fire behavior during experimental burning through an ANOVA (n = 3), with a significance level set at p < 0.05. (a) Average flame height (cm); (b) average soil temperature reached at a depth of 5 cm; (c) average fire propagation speed (m/minute). Different letters (a, b) indicate significant difference among burning treatments (p ≤ 0.05%, HSD Tukey).
Fire 06 00447 g003
Figure 4. Effects of fire treatments at different sampling times on bulk density (Bd) (a), porosity (b), soil moisture (c), soil organic matter (SOM) (d), soil organic carbon (SOC) (e), soil pH (f), soil ammonium nitrogen SAN (g), Mg (h), P (i), K (j), and Ca (k). Mean values and standard deviation are shown. Different letters indicate significant differences among burning treatments for different sampling times (p < 0.05, Tukey HSD).
Figure 4. Effects of fire treatments at different sampling times on bulk density (Bd) (a), porosity (b), soil moisture (c), soil organic matter (SOM) (d), soil organic carbon (SOC) (e), soil pH (f), soil ammonium nitrogen SAN (g), Mg (h), P (i), K (j), and Ca (k). Mean values and standard deviation are shown. Different letters indicate significant differences among burning treatments for different sampling times (p < 0.05, Tukey HSD).
Fire 06 00447 g004
Figure 5. Ordination of the different burning treatments for the two sampling times in the space defined by the PC1 and PC2 axis of the PCA analysis carried out with physico-chemical soil properties. Coordinates are the mean of three replicates and bars represent the standard error of the mean.
Figure 5. Ordination of the different burning treatments for the two sampling times in the space defined by the PC1 and PC2 axis of the PCA analysis carried out with physico-chemical soil properties. Coordinates are the mean of three replicates and bars represent the standard error of the mean.
Fire 06 00447 g005
Table 1. Atmospheric conditions during the experimental burning.
Table 1. Atmospheric conditions during the experimental burning.
BlocksPlotTreatment Temperature (°C)Relative Humidity (%)Wind Speed (m/s)Wind Direction (°)Precipitation (mm)Solar Radiation (W/m2)
11T1 = uphill14.47712.19001209
2T2 = downhill14.3798900492
3Control plot14.5756.39001364
24T1 = uphill15.9728.99001376
5T2 = downhill15.7708.9157.501208
6Control plot16.1709.413501243
37T1 = uphill16.4708.913501155
8T2 = downhill16.571813501164
9Control plot16.2699.413501005
Average15.672.68.9117.501135.1
Table 2. Fine dead fuel moisture according to the Interagency Fire Use Module Field Guide and calculation of fuel load and fuel moisture content.
Table 2. Fine dead fuel moisture according to the Interagency Fire Use Module Field Guide and calculation of fuel load and fuel moisture content.
Fine Dead Fuel Moisture (Interagency Fire Use Module Field Guide)Fuel Loading and Humidity
Variables°CRelative humidity (%)Shading (%)Ignition probability (%)Wind speed (m/s)Danger indexWet biomass (kg m−2) Dry biomass (kg m−2)Fuel moisture content (%)
15.6 ± 0.272.6 ± 0.90–10408.9 ± 0.2Alert2.21 ± 0.3317.90 ± 0.3223.6 ± 1.5
Note: ± is the standard error.
Table 3. Effects of low-severity fire on soil physico-chemical properties of herbaceous páramo soil in southern Ecuador. Significant differences were evaluated with repeated measures. Significance level of α = 0.05.
Table 3. Effects of low-severity fire on soil physico-chemical properties of herbaceous páramo soil in southern Ecuador. Significant differences were evaluated with repeated measures. Significance level of α = 0.05.
VariableRepeated Measures ANOVA S1 and S2 (p-Value)
Burning (B)Sampling Time (T)B × T
Bd (g cm−3)0.28440.0263; S2 > S10.5487
Porosity (%)0.6550.0402; S2 < S10.3673
Soil water content (%)0.63760.0472; S2 > S10.2113
SOM (%)0.47220.0701; S2 = S10.8689
SOC (tC ha−1)0.95770.0687; S2 = S10.6164
pH0.01740.0428; S2 > S10.1354
SAN (mg/kg)0.3920.0384; S2 > S10.3496
Mg (Meq/100 g)0.82670.015; S2 > S10.7094
P (mg/kg)0.19990.4701; S2 = S10.7608
K (Meq/100 g)0.13520.1141; S2 = S10.006
Ca (Meq/100 g)0.23790.8305; S2 = S10.2037
Table 4. Pearson product-moment correlation coefficients between soil physico-chemical properties and scores of the PC1 and PC2, which explained, respectively, the 57.3% and 23.4% of the variance in the PCA analysis. Significant coefficients are noted with *.
Table 4. Pearson product-moment correlation coefficients between soil physico-chemical properties and scores of the PC1 and PC2, which explained, respectively, the 57.3% and 23.4% of the variance in the PCA analysis. Significant coefficients are noted with *.
VariablePC1PC2
Bd (g cm−3)0.348 *−0.154
Porosity (%)−0.349 *0.232
Soil water content (%)−0.1010.498 *
SOM (%)0.367 *0.089
SOC (tC ha−1)0.386 *−0.065
pH−0.319 *−0.009
SAN (mg/kg)0.374 *0.165
Mg (Meq/100 g)0.367 *0.050
P (mg/kg)0.2620.298 *
K (Meq/100 g)−0.0800.501 *
Ca (Meq/100 g)0.1050.541 *
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

Carrión-Paladines, V.; Fries, A.; Hinojosa, M.B.; Oña, A.; Álvarez, L.J.; Benítez, Á.; Rodríguez, F.L.; García-Ruiz, R. Effects of Low-Severity Fire on Soil Physico-Chemical Properties in an Andean Páramo in Southern Ecuador. Fire 2023, 6, 447. https://0-doi-org.brum.beds.ac.uk/10.3390/fire6120447

AMA Style

Carrión-Paladines V, Fries A, Hinojosa MB, Oña A, Álvarez LJ, Benítez Á, Rodríguez FL, García-Ruiz R. Effects of Low-Severity Fire on Soil Physico-Chemical Properties in an Andean Páramo in Southern Ecuador. Fire. 2023; 6(12):447. https://0-doi-org.brum.beds.ac.uk/10.3390/fire6120447

Chicago/Turabian Style

Carrión-Paladines, Vinicio, Andreas Fries, María Belén Hinojosa, Andy Oña, Leticia Jiménez Álvarez, Ángel Benítez, Fausto López Rodríguez, and Roberto García-Ruiz. 2023. "Effects of Low-Severity Fire on Soil Physico-Chemical Properties in an Andean Páramo in Southern Ecuador" Fire 6, no. 12: 447. https://0-doi-org.brum.beds.ac.uk/10.3390/fire6120447

Article Metrics

Back to TopTop