Soil Vulnerability Indicators to Degradation by Wildfires in Torres del Paine National Park (Patagonia, Chile)

Wildfires are a common phenomenon globally, nonetheless some ecosystems are more adapted than others. In this work, we show results of a study conducted in Torres del Paine National Park (Chile) that suffered a big forest fire in 2011, affecting 17,666 ha. The Park is located in the transitional forest-steppe zone whose annual rainfall varies between 300 and 1,500 mm. Plant communities goes from Patagonian steppe, pre-Andean scrub to Magellanic forest. Based on vegetation community and trying to cover all representative sectors affected by the fire, 5 areas of the park were sampled in 2019 following the transects where a vegetation recovery study has been monitored since the fire occurred. Soil parameters used to evaluate were: texture, soil water repellency (WR), organic matter (OM), and aggregation, including total content of macroaggregates (TCA) and their stability (AS). The results showed high values of OM, however three of the 5 areas showed statistically lower values of OM in burned samples. WR was present (from slight to severe) in the 75% of the samples, and without differences between burned and unburned. The correlations analyses indicated that WR is more related with the OM type than with quantity, since better correlations were obtained when only samples from the same area were included in the analyses. The results of aggregation indicated that soils have a poor structural development. The TCA varies from 16 to 50%, and the AS is not very high (average of 66%), being the highest in the area with lowest TCA and more affected by the fire and erosion processes. This suggests that the higher values of AS are a consequence of the destruction and loss of the less resistant fraction after the fire. The combination of the high sand content (low specific surface area) and high OM make them very susceptible to develop WR. Since these soils have a scarce development with a poor structure, the combination of the WR and the poor soil structure make them very vulnerable to erosion processes after the fire. Measures to protect the soils or accelerate the recovery are recommended in these areas when new human caused wildfires could occur.


INTRODUCTION
Wildfires are a common and natural phenomenon globally and have played a vital role in the history and evolution of ecosystems (Pausas and Keeley, 2009). However, some ecosystems are more adapted to this perturbation than others due to the fire regime relationship with climate conditions. Chile is a country with large forest areas, where anthropogenic forest fires are widespread; they not only occur in the Mediterranean climate region in central Chile but more recently also in more humid and cold climatic zones like the southern Patagonian Region (Salas et al., 2016;Úbeda and Sarricolea, 2016;Mcwethy et al., 2018). It is known that there were fires during the Pleistocene-Holocene in Patagonia-Tierra de Fuego (Moreno et al., 2018). Huber et al. (2004) suggest that the causes of those wildfires were both natural and anthropogenic. There is evidence of human occupation in Tierra del Fuego-Patagonia dated to ca 11 14 C ka BP (Borrero and McEwan, 1997). The most likely source of these early fires was the indigenous population (Holz and Veblen, 2011). These authors documented an increase in the number of fires after the colonization and permanent settlement of the region by Euro Chileans. However, there are evidences that the climate conditions of the area would have inhibited the spread and dimension of most of those fires (Mundo et al., 2017). Fire does not form part of the ecosystem dynamics in southern Patagonia, meaning a lack of resilience to wildfire, and patches of the Patagonian steppe are exposed and more susceptible to be colonized by invasive species (Vidal, 2012;Kitzberger et al., 2016). The ecotonal zones between eastern semi-arid steppes and Nothofagus spp. forests in western Patagonia results from broad landscape and ecosystem changes, intensified in the last 140 years, with anthropogenic forest fires, logging, exotic pasture establishment, and mining as the leading causes (Zegers et al., 2020).
Forest fires can affect several soil properties depending on fire intensity, severity, and soil type (Certini, 2005). Vegetation and soil recovery are directly affected by heat and/or indirectly as a consequence of the temporal elimination of soil biota, vegetation coverage, and the effects of ash incorporation on soil conditions (Bodí et al., 2014). Among the soil properties affected, some of the most relevant are soil organic matter (OM), water repellency (WR), and soil aggregation (Doerr et al., 2000;Mataix-Solera et al., 2011). It is known that some soils are much more vulnerable to fire and post-fire erosion processes than others due to their nature and properties, and therefore with a very different response to fire perturbations (García-Corona et al., 2004). The OM content can decrease or increase due to fire depending on other factors such as type of fire, its intensity, and severity (Almendros and González-Vila, 2012). Soil OM reduction can be produced in high intensity fires as a consequence of the combustion processes and the remaining OM suffers quality changes in its composition (González-Pérez et al., 2004). Low intensity fires and fires that spread through the crowns tend to increase the OM content by the incorporation of partially burned vegetation residuals (Mataix-Solera et al., 2002).
Following fires, runoff and soil erosion can increase, reducing the OM content in the topsoil (Larsen et al., 2009;Bonilla et al., 2014). Soil WR is a property that reduces the affinity of soils to water, affecting the infiltration-runoff balance and therefore the erosion rates . It can be naturally caused by some soil organic compounds or induced by the effect of fires (DeBano, 2000;Doerr et al., 2000). Previous studies conducted in Mediterranean soils demonstrated that their susceptibility to developing WR depends not only on soil temperatures produced by fires but also on soil texture, mineralogy of clay fraction and the quantity and quality of the soil organic matter . The measurement of soil structure and aggregate stability described soil vulnerability to erosion processes after wildfires. Previous studies have found different patterns of aggregation in fire-affected soils (Mataix-Solera et al., 2011). The decreases of aggregate stability are related to high fire severity and loss of OM content (Badía and Martí, 2003), and the increases could be explained by changes in WR at certain temperatures (García-Corona et al., 2004;Terefe et al., 2008), or some specific mineral transformations with high fire severity (Giovannini et al., 1988;Jiménez-Pinilla et al., 2016). Consequently, there is no general pattern of soil behaviour to wildfire since many factors are involved.
Wildfires are relatively rare in cold temperate rainforest, but they can strongly change forest landscape (Zegers et al., 2020). In Chile, there are reported evidence about the mid-and short-term effect of fire mostly on vegetation in the Mediterranean region (Úbeda and Sarricolea, 2016). Under new climate change scenarios, the frequency of fires will increase in Patagonia, with anomalous dry conditions (Holz and Veblen, 2012). Understanding the mid-term effects of fire on soil and vegetation is necessary to implement effective restoring practices. Torres del Paine National Park, a world biosphere reserve, is one of the most emblematic parks in the Chilean Patagonia. It is located in the southwestern part of the region and preserves a sensitive ecosystem in a landscape dominated by evergreen shrubs, Nothofagus forests and grasses that resist harsh wind and weather. Between 1985 and2011, 29 fires caused by humans were recorded in the Park, affecting more than 36,000 ha, three of them were the most significant events produced in recent years (1985, 2005, and 2011). Following the last fire, Bonilla et al. (2014) reported no difference in soil organic content when comparing burned with unburned sites. However, the burned sites showed a high soil erosion potential. In this work, we analysed the results of soil properties from samples collected in 2019 in the Torres del Paine National Park (Chilean Patagonia) from different areas that suffered a large forest fire during the end of 2011 (Medina, 2013).
The objectives of this research are: 1) to evaluate the medium-term impacts of forest fire in some soil properties in a cold temperate rainy climate, like Torres del Paine, and 2) to compare the vulnerability of these soils to this perturbation with other environments also affected by wildfires but much more studied and resilient such as the Mediterranean.

Study Area, Fieldwork and Soil Sampling
The Torres del Paine National Park has a total surface close to 181,000 ha and it is in a transitional forest-steppe zone in the Region of Magallanes and Chilean Antarctic (50°58′59″S, 72°57′59″W) ( Figure 1). The climate is a temperate cold rainy zone without dry season, the temperatures range below 15°C, being this the maximum temperature for the warmest month (January), and a fairly rigorous winter with an average minimum temperature for the coldest month (July) never lower than −2.5°C. Precipitation from westerly winds shows a marked gradient from west to east, and annual rainfall varies between 300 mm in the east part of the park to 1,500 mm in the west side (Pisano, 1974). Within the park, it is possible to find four plant communities: Patagonian steppe, pre-Andean scrub (Xerophilous and Mesophilic), deciduous Magellanic Forest, and Andean desert. The Park was completely covered by ice during the last glacial maximum and therefore the soils are young. According to Marden and Clapperton (1995) the steppe area of Torres del Paine was not colonized by peat until after by ca 12 ka BP. The soils of the study areas, mostly developed over recent glaciofluvial sediments and with scarce development, are typical of cold regions. According to World Reference Base for Soil Resources system of classification the dominant soil types in the study areas are Phaeozems, Regosols and Leptosols (IUSS Working Group WRB, 2015).
In December 2011, a big forest fire occurred in the National Park, affecting 17,666 ha. Based on geographical location and vegetation coverage, in 2019, 5 areas of the Park were sampled following the transects where a vegetation recovery study had been previously monitored, to understand the status of the ecosystem and how fire and post-fire conditions affected the area ( Figure 1) (Armesto et al., 2012). The record was concentrated on sites affected by the 2011 fire considering: Cover of native and exotic herbaceous plants and regeneration of seedlings, and regrowth of shrubs and trees. Landsat satellite images (FORMOSAT-2) from 2011 and maps made by the National Park administration from areas previously affected by wildfires were used to select the vegetation and soil monitoring sites. These areas were chosen to cover all the range of plant communities, considering the accessibility and the presence of unburned control areas on each of the selected locations. The main difference between the study areas was the vegetation, being a transition between Patagonian steppe and xerophilous pre-Andean scrub in study area 1, the driest sector, area being more dominated by scrub (xerophilous and mesophilic), study areas three and four were predominantly of mesophilic scrub and forest, corresponding with the western sector of the Park, as indicated more humid, and area 5 low-density of Magellanic forest ( Figure 2). The main species present in each of the study areas are included in Table 1. After the fire in 2011, a field assessment of soil fire severity was done using the field guide for mapping post-fire developed by the USDA-Forest Service. The burned sites were classified into three levels; low severity, moderate severity and, high severity of above-ground vegetation and below-ground (Robichaud et al., 2000;Pearsons et al., 2010). The burned soil severity evaluation also considered three categories (Pearsons et al., 2010). 1) Low soil burn severity: Surface organic layers are not entirely consumed and are still recognizable. Structural aggregate stability is not changed from its unburned condition. Roots are generally unchanged because the heat pulse below the soil surface was not significant enough to consume or char any underlying organics, 2) Moderate soil burn severity: Up to 80% of the pre-fire ground cover (litter and ground fuels) may be consumed but generally not all of it. Fine roots (∼0.1 inches or 0.25 cm diameter) may be scorched but rarely consumed over most of the areas, and 3) High soil burn severity: All or nearly all of the pre-fire ground cover and surface organic matter (litter, duff, and fine roots) is generally consumed, and charring may be visible on more prominent roots. According to this and with the results of the assessments in all the sampling sites, for the five study areas, fire severity was classified as follows: low severity for area 1, low to moderate in area 2, moderate to high in area 3, and moderate severity in areas 4 and 5. Each study area includes a control (unburned) site as a reference (numbers 1, 13, 28, 32, and 34 of Figure 1). The others were affected by the fire of 2011, and only three of the points from area 2 (numbers 14, 15, and 16) were also affected by a previous fire in 1985. There is no overlapping with the fire of 2005 that affected an area of the park located in the Park's northern section. Vegetation cover was estimated in a 2 m radius (12.56 m 2 ) circular plots randomly located in each point (n 4 plots/ sampling point × 34 soil sampling points 136). We estimated the percentages of soil occupied by litter, standing trees, tree stumps, and/or burned shrubs, trunks, herbaceous cover and bare soil in each one of the plots. All individuals with phanerophyte habits were classified into five categories: burned with no leaves present, burned with burned leaves present, burned with green leaves, burned with sprouts, and not burned. This manuscript use the percentage of soil cover (or the percentage of bare soil) in 2019 to compare with the initial data collected in 2012, just after the fire occurred.
At each of the 34 sampling points (Figure 1), composite soil samples from the first 5 cm depth of mineral soil were taken in November 2019. Two or three composite samples per sampling point, by pooling subsamples from different sites, were collected for laboratory analyses, resulting in a total of 74 soil samples, including burned and unburned soils in the 2011 fire.

Analytical Methods
In the laboratory, the following soil properties were analysed: pH, electrical conductivity (EC), soil texture, soil water repellency (WR), organic matter content (OM), and aggregation, including total content of macroaggregates (TCA; % of sample that are forming macroaggregates) and their stability (AS; % of macroaggregates that resist the energy of a rainfall simulation of known energy). Soil pH was measured in deionised water (1: 2.5 w/v). EC was also measured in deionised water (1:5 w/v). Texture determined by the Bouyoucos method (Gee and Bauder, 1986). OM content was determined by potassium dichromate oxidation (Nelson and Sommers, 1982). The water drop penetration time (WDPT) test (Wessel, 1988) was used to measure the persistence of WR. Prior to WR assessment, approximately 15 g of each soil sample was placed on separate 50-mm diameter plastic dishes and exposed to a controlled laboratory atmosphere (20°C, ∼50% relative humidity) for 1 week in order to eliminate potential effects of any variations in preceding atmospheric humidity on soil WR and in accordance with Doerr et al. (2005). This involved placing three drops of distilled water (∼0.05 ml) on the sample surface and recording the time required for their complete penetration. The average time for triplicate drops is taken as the WDPT value of a sample, and then classified according to Bisdom et al. (1993) and Doerr et al. (1998), as follows: wettable (<5 s), slightly water repellent (5-60 s), strongly water repellent (60-600 s), severely water repellent (600-3,600 s), and extremely water repellent (>3,600 s). Soil aggregate stability (AS) and total content of macroaggregates (TCA) were studied in the macroaggregate fraction of 4-0.25 mm and determined with the rainfall simulator method according to Roldán et al. (1994) and based on the method of Benito and Díaz-Fierros (1989). This method examines the proportion of macroaggregates that remained intact after a soil sample was subjected to artificial rainfall of a specific energy. Four grams of soil material (4-0.25 mm) were placed into a 0.25-mm sieve and exposed to an energy of 270 J m −2 applied by an artificial rainfall. After treatment, the material remaining within the sieve was dried and weighed. The material on the sieve, consisting of remaining macroaggregates, mineral particles and organic debris, was then washed until only mineral particles and organic debris remained. The dry weight of the remaining mineral particles and organic debris material allowed the calculation of the difference in aggregate weight (%) within a sample before and after artificial rainfall (AS). The method also allows one to calculate the proportions of sample that forms macroaggregates by differences of weights (TCA). Results are based on the mean value of three replicate experiments per soil sample.

Statistical Procedures
The fitting of the data to a normal distribution was checked with the Kolmogorov-Smirnov test. When necessary, analytical data was transformed using logarithms to assure normal distribution. The t-test for independent samples was used to compare the means of parameters between the burned and control soil groups per study area. Pearson correlation coefficients (r) were calculated to quantify the linear relationship between parameters. A Principal Component Analysis (PCA) was performed with results of all soil samples in order to study the relationship between the soil parameters measured regarding the study areas using the "FactoMineR" package implemented in R (Le et al., 2008). The vegetation cover percentage comparison (2012 vs 2019) per study area was evaluated with non-parametric Kruskal-Wallis test. All statistical analyses were performed using the RStudio v.3.6.2 (RStudio Team, 2020).

RESULTS
The results of soils characterization (Table 2) showed that two of the areas have a loam soil texture and sandy loam in the other three, pH ranging from 5.4 to 6.7, and low values of EC. Area 1 was observed as the most vulnerable to erosion processes being the area where higher percentage of bare soil (43% in average) remained 8 years after the fire ( Figure 3A). It is also the only area where under unburned conditions soil bare (18%) was found. The other study areas were recovered by vegetation and after 8 years of the fire, the percentage of bare soil is low, ranging from 4 to 14% ( Figure 3A), however there is still significant differences also in areas 3 and 4. Table 3 show the comparison of vegetation cover in burned soils between 2012 and 2019 per study area. In all study areas except area 1 it was observed statistical differences with a significant increase in vegetation cover. The area 1 however, despite the low fire severity registered and starting in 2012 with a soil coverage close to 50%, the values did not reach 60% in 2019 without significant differences between years.
The OM content of the soil samples is in general high, with an average of 10.5% and with some soil areas reaching more than 20% ( Figure 3B). Three of the five areas showed statistically lower values of OM in burned samples (areas 1, 3 and 4), and in one of them, area 5, levels were higher after burning. No differences were found for area 2. Soil WR was present in the 75% of the samples with WDPT classes ranging from slight to severe water repellent, and without differences between burned (222 ± 392 s) and unburned samples (298 ± 318 s; p > 0.05). All study areas showed WR (79, 67, 76, 75, and 100% of water repellent samples for areas 1 to 5, respectively; Figure 4). The correlations analyses (Table 4) indicated that WR is positively correlated with OM content ( Figure 5A), but this correlation is higher when only samples from areas where the vegetation structure and composition is more homogeneous are included. As an example, it is showed the correlation between both parameters in area 1 (r 0.852; p < 0.01, Figure 5B). The results of aggregation indicated that these soils have a poor structural development. The TCA varies from 16 to 50% ( Figure 6A), and the AS measured were not very high (average of 67%; Figure 6B), being the highest in the area 1 which is the area with lowest TCA. In the area 1 the TCA was statistically lower for burned soils and AS higher. The OM content was positively correlated with TCA but not with AS and a negative correlation was found between TCA and AS. Not correlations were found between AS and OM and between AS and WR ( Table 4). The PCA performed explained a total of 76.8% of variability (Figure 7) showing that some areas like 5, and specially 1, have more homogeneous data and they clearly differentiate from the other areas. Area 1 is clearly separated for the soil aggregation parameters studied (TCA and AS). Area 2 seems to collect the most variability of all study areas. Area 5 is related with higher values of WR and AS, while areas 3 and 4 are related with high values of TCA and OM.
Regarding the effect of fire recurrency on soil properties, for this study, it is only limited to a section of area 2. The compartive of the results for those sampling points (burned in 1985 and 2011), with those of the same area only affected by the fire in 2011 showed no statistical differences for the parameters selected for this study.

DISCUSSION
Our soil sampling results in various vegetation areas burned and not burned in 2011 show that the effect of fires on soil properties depends on local site conditions and vegetation recovery.

Soil Organic Matter
The OM content of the studied soils is high but revealed decreases in three of the five study areas as consequence of the fire and/or post-fire erosion processes. Fire can affect directly to the OM content depending on fire intensity and severity, in fact previous studies have observed both increases or decreases depending on the fire behaviour (Mataix-Solera et al., 2002). The area 1 showed statistically lower OM content in burned soils despite this is the area where the fire had low severity, therefore this decrease should be more probable due as consequence of the post-fire erosion processes. In the cases of areas 3 and 4 where fire severity was higher, the decrease observed could be caused for both, higher combustion and post-fire erosion. Bonilla et al. (2014) studied burned areas in Torres del Paine, some of them affected in the same fire of this research but just 1 year after the fire. They reported no difference on soil organic content comparing burned with control (unburned) sites, one reason could be because the depth of the soil samples collected (0-20 and 20-40 cm). The fire effects on soils tend to affect only the first mm or cm of mineral soil and these effects are masked if soil samples that are taken with more depth (Neary et al., 2005). In this case the samples are from 0 to 5 cm, thus more sensitive to reflect the changes. In any case, taking into account the fire severity reached, it is considered that the detected losses of OM in the areas 1, 3, and 4, eight years after the fire are more probably due to post-fire erosion processes than from a direct effect of combustion during the fire. This is also supported from the results of the erosion rates estimated in the same study area by Bonilla et al. (2014) using data from fire affected sites in 2005 with the Water Erosion Prediction Project (WEPP Model). Their results showed that water erosion increased after the fire, mainly due to the effect of vegetation removal but also due to changes in soil properties and erodibility. These authors reported that the sediment yield after fire increased from 0.02 to 9.25 t ha −1 year −1 and produced an estimated increase in soil organic carbon losses from 0 to 0.54 t ha −1 year −1 . Other studies have reported reduction in OM content by post-fire erosion when vulnerable soil conditions and low vegetation cover persist in time (García-Orenes et al., 2017). Area 1 was observed as the most vulnerable to erosion processes being the area where higher percentage of bare soil (43% in average) remained 8 years after the fire (Figures 2A,B, 3A). It is also the only area where under unburned conditions some of soil bare (18%) is found. This also supports the affirmation that especially in the area 1 the OM loss has been produced by topsoil erosion. The other study areas were recovered by vegetation and 8 years after fire the percentage of bare soil observed were low (Figure 2A). On the other hand, the increased OM content in area 5 and the results of area 2 without statistical differences are indicative of how variable the effect of fire can be in the organic fraction of the soil.

Soil Water Repellency
A high percentage of soil samples (75%) showed WR, without differences between burned and unburned soils. This reveals that WR is a natural property in these soils. Many studies have found natural soil WR under different vegetation types, mainly associated with evergreen tree types and particularly, trees with considerable amounts of resins, waxes or aromatic oils. However, soil WR has also been detected under deciduous trees, shrubs and grass species (Doerr et al., 2000). The combination of the high proportion of sand, high OM content and acid pH make the soils of Torres del Paine very susceptible to develop WR (Doerr et al., 2000), although the variability of WR persistence observed is elevated (Figure 4). Mataix-Solera et al. (2014), found that small variations in soil properties such as texture, mineralogy of clay fraction, and organic matter quantity and quality can control the development and persistence of WR, being the sand content, one of the most relevant factors. These soils have a high sand content with textures that varies from loam to sandy loam. Sandy soils are more prone to develop WR because of the lower specific surface area to be covered by hydrophobic substances (González-Peñaloza et al., 2013). The pH of the soils of this study varies from 5.7 to 6.7 (Table 2), and this can also contribute to the presence of WR. Many studies showed that acidic soils are more prone to develop WR than alkaline (e.g., Dekker and Jungerius, 1990;Doerr et al., 1998;Varela et al., 2005). There are different reasons that can contribute to this; alkaline conditions enhancing the solubility of the humic acid fraction (implied in WR  TCA, total content of aggregates; AS, aggregate stability; OM, organic matter; WR, soil water repellency. **, *: significant level for p < 0.01, p < 0.05, respectively; ns, not significant at p > 0.05. Spanish Journal of Soil Science | Published by Frontiers July 2021 | Volume 11 | Article 10008 development), the dependency of the pH of water contact angles of fulvic and humic acids and their negative relationship with pH (Lin et al., 2006), and the enhanced activity of fungi releasing hydrophobic compounds in soils with lower pH (Jex et al., 1985;Paul and Clark, 1996;Weyman-Kaczmarkowa and Pędziwilk, 2000). Regarding the relationship between WR and OM, numerous researchers have also attempted to establish general relationships, but with some inconsistent results. Some authors found positive correlations (e.g., Wallis et al., 1990;McKissock et al., 1998). However, Teramura (1980) found a negative relationship and DeBano (1991) and Wallis et al. (1993) no relation between the two variables. This inconsistency has been explained because only some compounds of soil organic matter induce WR, and therefore not always is correlated with the amount, particularly when soils with different organic matter quality are compared. In this study case the correlations analyses ( Table 4) support also the idea that WR is more related with the type of OM than with the quantity, since better correlations were obtained when only samples from areas with similar vegetation were included in the analyses, and a worse correlation was obtained for areas with more heterogeneity in forest (herbs, scrub, trees) or when all the samples from different sites were pooled together. An example is showed in Figure 5A, with all samples from the different areas pooled together, (r 0.356; p < 0.01), and how the correlation improves when only samples from area 1 (steppe vegetation, Figure 5B) are included (r 0.852; p < 0.01).
Fire effects on WR can be very variable since they can induce it in wettable soils, enhancing or destroying in original water repellent ones, depending on many factors, being the temperatures reached in the soil one of the most determinant (DeBano, 2000), but in this case, 8 years after the fire, no significant differences were found in WR between burned and unburned soils. The natural presence of WR in these soils and the time passed since the fire could be masking any fire effect on soil WR.
Different studies have suggested that the natural presence of soil WR caused by plants is thought to be a competitive strategy of species, improving water conservation by channelling water to greater depth (e.g., Scott, 1992;Moore and Blackwell, 1998), and reducing evaporation of water in the surface layer (Doerr et al., 2000). This property, however, becomes a problem during the temporal window from the fire event (bare soil) to the point when the recovery of vegetation occurs, since erosion processes are accelerated due to increased runoff rates (Shakesby and Doerr, 2006).

Soil Aggregation
These soils showed a low structural development. The proportion of soil that are forming macroaggregates (TCA) varies from very low (16% of average in area 1) to medium (50% in area 3), and the stability of macroaggregates (AS) showed results with high variability among study areas. There are different methods to measure AS and between them the data are not always comparable, since they are relative to the energy applied and the way to do it (drop impacts, water immersion, ultrasounds, etc.). In this case a laboratory rainfall simulator method was used, and to compare the results previous studies that adhere to this methodology are used. Chrenková et al. (2014) using this methodology, studied aggregation in different types of Mediterranean soils (including different range of soil textures) and with different land use (forest and agricultural). They found that the TCA in Mediterranean forest soils varied from 55 to 95%, much higher that the results obtained here for soils of Torres del Paine. Regarding AS, results from different studies in Mediterranean forest soils showed a range from 45 to 95% (Zornoza et al., 2007;Chrenková et al., 2014), or from 46 to 76% in A horizons of reforested soils affected bt gully erosion in Slovakia (Dlapa et al., 2012). Fire can affect in a different way to the soil parameters related to aggregation (Mataix-Solera et al., 2011). Many authors have found that high severity fires can destroy directly FIGURE 6 | Box-plots of (A) total content of macroaggregates (TCA %) and (B) aggregate stability (AS, %) per study area comparing burned soils vs controls. **, *: significant level for p < 0.01 and p < 0.05, respectively; ns, not significant at p > 0.05.

Spanish Journal of Soil Science | Published by Frontiers
July 2021 | Volume 11 | Article 10008 macroaggregates or reduce their stability making them less resistant to erosion (Zavala et al., 2010;Jiménez-Pinilla et al., 2016). However, others didn't found differences (Llovet et al., 2008;Zavala et al., 2010) and even in some cases AS increases have been reported (Josa et al., 1994;García-Corona et al., 2004;Mataix-Solera and Doerr, 2004). The different patterns of aggregation after the fire have been explained depending on the type of soil, the temperatures reached, and how this can affect to other soil properties related to aggregation such as fungi content, organic matter, water repellency or mineralogy (Mataix-Solera et al., 2011). In this study, 8 years after the fire, differences were found for both aggregation parameters only in area 1, with a decrease in TCA and an increase in AS, being the higher values precisely in this area with lower TCA. Taking into account the results of both parameters and the reduced observed in OM content in burned soils of this area, this suggests that the higher values of AS found are consequence of the destruction and loss of the less resistant fraction after the fire and erosion processes. This is also supported with the negative correlation obtained between TCA and AS (Table 4), the PCA (Figure 7), and with the results of percentage of bare soil ( Figure 3A), higher also in area 1, thus more soil directly exposed to rainfall and wind events. Despite many previous works has found positive correlations between AS and OM (Soto et al., 1991;Roldán et al., 1994;Badía and Martí, 1999;Dlapa et al., 2012), and between AS and WR (Chenu et al., 2000;Benito et al., 2003;García-Corona et al., 2004), this is not the case of this study (Table 4), and OM was only positively correlated with TCA and with WR. These relationships used to be strong in undisturbed and well aggregated soils (Martí et al., 2001), but becoming inconsistent and more complex in burned and degraded soils because all the factors involved, showing contrasted responses in some cases (Mataix-Solera et al., 2011). The lower proportion of soil forming macroaggregates, and the erosion processes, breaking some of them and increasing the relative proportion of the stronger ones, could be explaining the lack of correlation between OM and AS. For this reason, the evaluation of the aggregation only with AS data can lead us to erroneous interpretations. It is necessary to include measurements such as the TCA. This becomes of particular interest and necessary when studying soils with low degree of structure development and subjected to disturbances such as fires and subsequent erosive processes.

General Implications
The role that vegetation is playing in these soils, as a protector against erosive processes and through the dense network of roots, is fundamental given the low percentage of soil that is forming macroaggregates. There are many areas of the Park where vegetation was able to regrowth, however changes in soil properties might affect the regrowth of key species that FIGURE 7 | Scores and loadings for PCA performed with results of all soil samples per study area. AS, aggregate stability; OM, organic matter content; TCA, total content of macroaggregates; WR, water repellency.
Spanish Journal of Soil Science | Published by Frontiers July 2021 | Volume 11 | Article 10008 depend on seed production. One of the main concerns in the Park is that some tree species like Nothofagus pumilio and Nothofagus betuloides do not regenerate naturally and could be replaced by shrubs that resist fire and resprout vigorously (Mermoz et al., 2005). In addition to that, this study indicates that the soil should also be an object of concern given the vulnerability of degradation shown by some indicators. Evaluations with analysis of NDVI series for the 2011 burned areas showed that for the year 2014, the herbaceous cover of Patagonian steeps had an important recovery of coverage, not so the pre-Andean scrubs areas (Peña and Ulloa, 2017). For these reasons, measures to protect the soils such as mulching or accelerating the recovery are recommended in these areas when new wildfires caused by humans occur (Díaz Vidal et al., 1959;Keeley, 2009;CONAF, 2015;Díaz-Raviña et al., 2018;Girona-García et al., 2021). This will be of particular interest in the areas of the national park where annual precipitation is lower and both the soil and the ash from the fire are highly exposed to erosive processes, not only water but also wind erosion, since blizzard episodes are frequent in the area. This work promotes the achievement of Sustainable Development Goal 15, by addressing how sustainably manage forests and reverse land degradation.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
JM-S funding acquisition to conduct the study, conceived and designed research, conducted the fieldwork, collected data, supervised analyses and wrote the manuscript. EA and PV also participated in the funding acquisition, design of the research and fieldwork. EA participated in the fieldwork and collected data. LO, JG, VA, MG-C, and NG-F did soil laboratory and statistical analysis. All authors contributed to discussions of the results, revised and approved the final version of the manuscript.

FUNDING
This work was supported by funding by the POSTFIRE_CARE project of the Spanish Research Agency (AIE) and the European Union through European Funding for Regional Development (FEDER) (Ref.: CGL2016-75178-C2-1-R), UMH grant for international mobility and ANID PIA/BASAL FB0002 (Chile).