Ground Fire Legacy Effects on Water-Dynamics of Volcanic Tropical Soils

The forest ﬂ oor is a critical component in maintaining the life cycles of forest ecosystems. It normally includes organic soil horizons, known as duff and litter, which are prone to be rapidly consumed after ﬂ aming and smoldering ﬁ res. This work aims to understand the legacy effects of surface and ground ﬁ res on the in ﬁ ltration capacity of a volcanic forest soil. We studied ﬁ ve sites with ﬁ res recorded in the last 20 years. All of them are located in pine-oak forests of the volcanic mountain region in central Mexico with a temperate climate and Andic soil properties. Tension-in ﬁ ltration tests were carried out to determine hydraulic conductivity and the number of active macropores. After each test, cores were taken to evaluate in a laboratory setting, where soil water repellency at different moisture concentrations and the integrative dynamic repellency index were determined. Field-saturated hydraulic conductivity was moderately high in all sites, with mean values of 13 and 42 mm · h − 1 for burned and control plots, respectively. A non-linear relationship was found between recurrence and type of ﬁ res with the concentration of active pores and several dynamic water repellency parameters. This work con ﬁ rmed the presence of latent combustion in these temperate neotropical forests. The changes in soil water repellency and hydraulic conductivity detected do not necessarily imply an exceeded soil in ﬁ ltration capacity. However, many of the ﬁ res in this region are associated with increasing agricultural activities, so further studies are needed to determine if higher ﬁ re frequencies could exceed the resilience capacity of the soils triggering land degradation.


INTRODUCTION
Tropical mountain ecosystems are characterized by significant biological and lithological diversity, which plays a critical role in maintaining ecosystem services and supporting a variety of life forms (Bruijnzeel et al., 2011;Poorter et al., 2015).Among the various components of these ecosystems, the forest floor plays a critical role in maintaining the life cycles of forest ecosystems (Descroix et al., 2001;Keith et al., 2010;Neris et al., 2013).This typically includes organic soil horizons, known as duff (O e and O a horizons) and litter (O i horizons), and in some cases of abundant moisture, even peat (H a , H e , and H i horizons).Most meso-and microfauna and ectomycorrhizal fungi concentrate their life cycles in these first centimeters of forest soils (Prescott and Vesterdal, 2021;Peng et al., 2022).
Surface fires, which spread at a high speed, rapidly consume the understory without significantly changing the soil properties (Agee, 1993;Keeley, 2009).However, when the forest floor includes thick layers of peat, duff, or litter, slow smoldering combustion can take place after flame combustion, which can generate profound soil changes (Ogle and Schumacher, 1998;Rein, 2015;Santoso et al., 2022).These effects have been studied in tropical, boreal, and sub-boreal peatland ecosystems; however, there is a gap in this knowledge in neotropical mountain ecosystems (Santoso et al., 2019).
Different time intervals between fires could affect tree-seed sprouting (Busby et al., 2020) and, thus, the forest's community composition and ecological services (Miller et al., 2011).However, the effects on soil properties are still an emerging field where nonlinear relationships have been observed (Moghli et al., 2022).
Moreover, while it is well established that low soil water content in fire-affected water-repellent soils can exacerbate the risk of erosion (Doerr et al., 2000;Shakesby and Doerr, 2006;Bento-Gonçalves et al., 2012), the emergence and severity of water-repellency do not occur in a linear progression as it dries.Rather, it occurs as a complex phenomenon where slightly moist soils can be even more repellent than completely dry soils (Wallis and Horne, 1992;Majid et al., 2023).Parameters such as IRDI (Integrative Repellency Dynamic Index), gravimetric water content at maximum water repellency levels (w.max), or gravimetric water content where the water repellency starts to be negligible (w.min) are alternatives to capture this dynamism (Regalado and Ritter, 2006;Regalado and Ritter, 2008).
In this sense, this work arises as a response to understanding the legacy effects of surface and ground fires on the infiltration capacity of mountain forest soils and to evaluate the degree of resistance of the Andosols to up to two recurring fires.We will also validate the presence of smoldering fire in two plots with duff in their forest floor.The central hypothesis of this research is that there will be a negative effect of the fire recurrence on soil hydraulic properties with a positive correlation of the number of fires with soil water repellency, and a negative relationship with hydraulic conductivity.

Study Sites and Soil Description
We studied five sites, one of them an unburned control, and four affected by different fires between 2000 and 2018; showing blackening of the trunks up to a height of 1.5 m (Rodríguez Rodríguez, 2014;Olivares-Martínez, 2020).Their land tenure was collective (ejido) and private (Table 1 and Figure 1), so all landowners of all sites were contacted to obtain their informed consent before fieldwork.
Due to the serious insecurity caused by drug cartels in the region, ensuring a minimum security condition necessary to carry out the fieldwork was also an important criterion.This limited the final number of sites and samples studied, in order to ensure the safety and physical integrity of the staff who worked in the field.
The vegetation is composed by pine-oak forest over volcanic soils with slopes between 15 °and 35 ° (INEGI, 2017;Olivares-Martinez et al., 2023).The forest has dominance of Pinus pseudostrobus and co-dominance of Quercus rugosa, Quercus laurina, and Quercus crassipes (Olivares-Martínez, 2020).Tree densities are between 170 and 620 trees per hectare, canopies of more than 80%, and in these sites there are some forestry activities like resinating pine trees or collecting firewood without logging.
The burned sites were affected by small fires of less than 5 ha and the burned trees had flame scars up to 2 m.In two plots, in addition to the evidence of flames, there were charcoal residues between the deeper duff and down to 10 cm in the mineral soil, which suggested the presence of latent combustion (Olivares-Martínez, 2020).They were assigned codes according to their fire history status: capital letters to express if they were Surface (S), Ground (G) or Control (C) plots, and number for the number of fires in the 18 years period.
The burned soils have loamy textures with clay percentages between 9% and 20% and very high organic matter contents (Olivares-Martínez, 2020); the first 5 cm of the mineral soil have soil organic matter levels between 7% and 19.5%.The sites affected by ground fires had duff forest floors up to 12 cm in depth, with assumed soil organic carbon levels of more than 20%, by definition (IUSS Working Group WRB, 2022).On the other hand, the reference control plot has a silty loam texture, with 6% of clay, 7% of organic matter, thick organic layers of duff, and litter of 14 cm in depth (Rodríguez Rodríguez, 2014).All the soils have Andic properties and were classified as Silandic Andosols (IUSS Working Group WRB, 2022).At each site, the type of mulch was evaluated according to its thickness and mineralization/humification processes (Siebe et al., 1996), a soil pit was dug, and the genetic horizons were described (FAO, 2009).

Infiltration Tests
INDI tension infiltrometers were used to perform 15 to 20 infiltration tests per site to ensure adequate representativeness according to its log normal statistical distribution [Gómez-Tagle (Jr.)Ch. et al., 2011].The tests were conducted on the forest floor, which comprised mainly mineral A horizons, and in some cases organic O a horizons after removing the litter.The INDI is a low-cost self-made device to characterize the saturated & unsaturated field hydraulic conductivity and estimate the effective soil porosity (Eqs A1-A3) (Watson and Luxmoore, 1986;Gómez-Tagle Chávez et al., 2014).The INDI tension infiltrometer is inspired by designs from Špongrová et al. (2009) and Perroux and White (1988).Further details can be reviewed by the reader in Gómez-Tagle et al. (2011).We applied −9.0, −3.0, and −0.5 cm tensions.We considered geometric corrections and the non-linear numerical approximation from Wooding's and Gardner's equations (Eqs A1, A2) to know the saturated and unsaturated hydraulic conductivity of the field (Logsdon and Jaynes, 1993).The number of active pores and the volumetric percentage were estimated from the unsaturated infiltration curves and the Jurin Law (Eqs A3, A4) (Watson and Luxmoore, 1986;Gómez-Tagle Chávez et al., 2014).We assumed the dynamic water viscosity as 1.002E-03 kg•m −1 •s −1 (standard at 20 °C), the water density as 998.2 kg•m 3 (standard at 20 °C), and the gravity acceleration as 9.81 m•s −2 , and a macropore size of more than 0.5 mm (according to Luxmoore, 1981).
To ensure the complete contact of the base of the infiltrometer and the soil surface, we used a fine sandy marble contact bed (90% sand, 8% silt, and 2% clay; "Marmolina comercial Tipo I, Fina") (Gómez-Tagle Chávez et al., 2014;Rodríguez Rodríguez, 2014).In the case of steep slopes (more than 35%), we made small terraces to level the infiltrometer, ensuring the measurement occurred within the target soil depth.This is a reliable method with no risk of losses by lateral flow when there are no abrupt textural changes (Miyazaki, 1988).

Soil Sampling and Lab Analysis
A core of about 1 cm in thickness and 5 cm in diameter was taken inside the area where each infiltration test was performed.The cores were carefully put in Petri dishes to measure water repellency at different moisture levels, according to the Water Drop Penetration Time test and the Molarity of an Ethanol Droplet (Woudt, 1959;Letey et al., 2000).In each case, the integrative dynamic repellency index (IRDI), maximum water repellency angle (a.max or α max ), water repellency angle after drying the samples at 105 °C (a.105 or α 105 °C), the trapezoid integral area below the α(θ g ) curve (S), and gravimetric water content where the water repellency starts to be neglectable (w.min or θ g.min ), were calculated (Regalado and Ritter, 2005).
Equation Solving and Statistical Analysis R v4.2.2 statistical software was used to solve the Wooding & Gardner non-linear equations and to perform statistical tests (R Core Team, 2022).One-way ANOVAs were done to assess the effects of type of fire and recurrence onto the hydrophysical soil variables comparing with control plot.Two-way ANOVAs were carried out, excluding the control plot to explore the effect of the interaction between the factors of type (surface, and ground fire), and recurrence (one, or two fires inside the period).Log-normal adjustments to compare maximum likelihood estimators were performed to the four hydraulic conductivity-derived variables, the water repellency angle after drying the samples at 105 °C, and the maximum water repellency angle.A post hoc honestlysignificant-difference test by means of Tukey (HSD) was done in the significant models.

Morphology of the Burned Volcanic Soils
The four burned sites showed fire scars with charcoal residues in the bark of their trees and soil water repellency.The volumetric field water contents of the sites, when sampled, were from 50% to 102%.The main burned trees were pines, but in the G(1) plot, scars were also found on some oaks.
The S(1) and S(2) plots presented a Mull humus (duff).On the other hand, the G(1) and G(2) plots had a Mor humus (duff), and there was evidence of smoldering combustion in both.There were visible charcoal residues found to a depth of more than 10 cm below the surface.A system of galleries and macropores with the rest of the root charcoal down to more than 50 cm was observed on the G(2) plot (Figure 2).

Hydro-Physical Effects From Ground and Surface Fire
A total of 85 water infiltration tests were performed.Nine tests failed due to high water repellency levels on the field, so only 76 tests were completed.The field-saturated hydraulic conductivity was from moderately high to high in all sites, the control plot had 42 mm per hour, much higher compared to the burned plots, with a mean hydraulic conductivity of 12.88 mm per hour (Table 2).
The water flux via macropores decreased after fires, like the field-saturated hydraulic conductivity.Burned soils could have a macropore flux 20 times less than the reference soil.Nevertheless, the number of active pores (N M•m-² ) had contrasting scenarios across the different types and fire recurrencies (Tables 2, 3).
A total of 63 water repellency tests were performed throughout different levels of soil moisture content.Almost all the samples showed water repellency to some extent; only tree samples in G(2) plot were completely hydrophilic.The mean a. max throughout the samples was 100.7 °with no significant differences when type and recurrence of the fire were considered (p < 0.05, d. f. = 1, F = 5.53), significant differences could be found only considering a one-way ANOVA model (p > 0.001, d. f. = 4, F = 6.55).The rest of the water-repellency parameters showed statistically significant differences (Figure 3).
The control site showed a higher a.105 dry-soil water repellency than the burned plots: 100.6 °against a range of 90.7 °-96.5 °, with only one significant difference between C(0) and S(1).The a.max, and the S parameters showed similar patterns.The w.min showed significant differences after the first fire, according to the type of fire but not when compared to the control plot (Figure 3).On the other hand, the IRDI showed significant differences in the control plot, specifically after a second surface fire.However, a second ground fire reestablished the IRDI to reference values.
The hydraulic soil properties and the active porosity showed a consistent and significant relationship to the type and the cumulative number of fires (Table 3).On the other hand, the soil water repellency dynamic parameters showed differences mainly according to the type of fire, and just the integrated index IRDI showed a significant relationship with fire type and recurrence (p > 0.01, d.f.= 1, and F = 11.06).
After two fires, the water repellency could have contrasting dynamics depending on the type of fire and the type of forest floor.For example, the duff-covered G(2) soil showed a similar total water repellency to the control plot with the IRDI index but with higher maximums at a. max and even a lower dry-soil water repellency at a.105 values; on the other hand, the S(2) showed a decrease in the IRDI, similar values to the control plot on w. min and a. max, but a lower a.105 (Figure 3).    1.

Fire Legacies on Soil Memory
The presence of residues of organic duff layers (O a and O e genetic horizons) were found only on ground-fire-affected plots.When the soil has high organic matter (nearly 20%), it is highly likely to develop smoldering combustion after a surface fire (Rein, 2015).It is possible that the increased number of galleries and macropores on the G(2) plot after a first fire may have enhanced the oxygen availability so that a second fire could burn deeper layers and leave deep charcoal residues such as those found in the field.Of all the components of a fire regime, the rate of spread of fire (and therefore the time of residence of the heat) is the one that can have the greatest effect on the soil (Agee, 1993;Keeley, 2009).In this sense, the two sites with clear evidence of smoldering combustion presented higher hydrophobicity at deeper layers (Olivares-Martínez, 2020).These plots had scars on pines and even on oaks (Quercus spp.), whose bark is less prone to be burned unless there are lower fire spread rates.The longer the time of residence, the deeper and greater the fire severity in the soil (Neary et al., 1999;Varner et al., 2005;Rein, 2015).
Smoldering fires can increase soil temperatures by more than 60 °C down to 20 cm on soil mineral horizons (Varner et al., 2005) and more than 50 cm on peat layers (Wösten et al., 2006).So, it is logical that a ground fire would show lower standard errors than a surface fire in its water-repellency-related properties since the former has more homogeneous fires.The G(2) plot could have presented deep charcoal residues and big galleries due to the smoldering of whole trees' root systems adding to the organic soil smoldering combustion.
Organic layers can start to smolder even with gravimetric water contents up to 120% (Neary et al., 1999).When sampled, the field water contents of the ground-fire-burned sites were from 79% to 102%, so they could have burned at that time even when they seemed wet.Soil humidity is a factor that regulates both ease of burning and soil water repellency and water infiltration rates (Jordán et al., 2013;Mataix-Solera et al., 2014;Marín-Castro et al., 2016).It suggests that the humus type of the forest floor drives the water thresholds on smoldering to fire dynamics and the charcoal imprinted legacies.
To understand better why the control plot and G(2) plot presented the higher IRDI values is necessary to understand how the index is structured.The index is constructed with the S and the w.min variables to describe the dynamic of water repellency according to soil humidity, so its changes can be due to the size of the S parameter, which tells us about the total area under the water repellency curve, and in an inversely proportional way, the index changes with the w.min parameter (Regalado and Ritter, 2005).
Moreover, other indirect effects on the IRDI could be explained by the size of the differences between a.max and a.105 (difference known as α err ) (Regalado and Ritter, 2005).Although a. max did not show great differences between sites except for the Sup(1) plot, a.105 was sensitive to the number and type of fire.It is possible that the presence of resins and other tree exudates could increase the dry water repellency at the control plot even more than the burned plots (Rodríguez Rodríguez, 2014).

Dynamic Water Repellency Remarks
The IRDI index was proposed as an attempt to understand better the non-linear dynamic of water repellency, and it should be interpreted as a shape index more than a total repellency index (e.g., S parameter), which tells where there are major differences between the hydrophobicity of dry and moist soil in narrow moisture gradients (Regalado and Ritter, 2005).A high IRDI value does not necessarily imply a stronger water repellency, it tells us about the relationship between the total water repellency and how easy it disappears with humidity.
Although no significant differences were seen in S and w.min, when combined in the IRDI it was possible to find significant differences generated by both recurrence and type of fire (Figure 3).The soils have imprinted different shapes in the repellency curve according to their fire history, so we can understand the IRDI as a "repellency signature."Moreover, although the plot two ground fires had increases in the persistence of soil water repellency, it reestablishes the form of this relationship, so the IRDI is very similar to the control plot.It would be interesting to further expand on this pattern by comparing with more sites that have a similar fire history.
It could be confusing when trying to consider the IRDI in hydrologic terms; the most prone sites to have problems with water infiltration and erosion will be those with the greatest water repellency area (S), a combination between high maximum water repellencies (a.max) and high minimum humidities (w.min) with low standard deviations.Water repellency is a patchy-distributed spatial variable (Woods et al., 2007), so greater standard errors such as those found at surface fires (bigger error bars at S and w.min on Figure 3) suggest a strong relationship with Hortonian water fluxes to a more hydrophilic matrix.Surface fires tend to have higher spread rates, so heat transmissions to soil are more scattered, and, consequently, the soil changes (Agee, 1993;Mataix-Solera et al., 2007).
Although without noticeable effects on the number of active macropores on these volcanic soils, the amount of water flowing through them was affected by the type and number of fires.These reductions may be linked to the persistence of water repellency without significant changes in organic matter content (Olivares-Martínez, 2020).The w.min, as a measure of the persistence of hydrophobicity in humid conditions, is important to show that even after the first rains, soil can still have water repellency.That could be important when modeling surface water runoff and soil erosion processes (Cerdà and Robichaud, 2009).

Resilience and Degradation
Although only the C(0) and the G(2) plots were above the global hydraulic conductivity average according to their textural class (Gupta et al., 2020), it could be common that Andosols have low hydraulic conductivities, especially when water repellency is present (Farahnak et al., 2019).Neris et al. (2013) have shown that Andosols beneath pine vegetation cover tend to have lower infiltration capacity due to the water repellency associated with the type of duff of its forest floor.This makes them more susceptible to erosion by runoff, despite the high resistance to degradation that these soils may have.
Average rainfall in the region can be up to 9 mm daily during the rainy season (CNA and SEMARNAT, 2020).If we assume that convective precipitation events could concentrate the day's rainfall in 1 h, all the rainwater can be infiltrated with the cumulative effect of one or two surface fires.The assumed runoff exceeded the water infiltration capacity only in the G(1) plot, which has a mean hydraulic conductivity of 6.5 mm per hour, the rest of the burnt plots can buffer this mean rainfall (Table 2).So, this could indicate that even with a reduced infiltration and enhanced water repellency, one or even two fires in 18 years have no serious repercussions for surface runoff and water erosion.In any case, further studies are necessary to confirm this result.
Despite the infiltration increase after a second ground fire, it should be considered that a higher standard deviation (patchiness), coupled with a low percentage of macropore flow and a higher repellency, could imply preferential finger flows (Ruthrof et al., 2019).This type of water movement may lead to higher risks of nutrient and contaminant leaching at recursive soil points (Ritsema and Dekker, 1994;Debano, 2000) but higher water infiltration potential at the landscape level (Doerr et al., 2000).Therefore, there may be an increase in the ecological function of water supply but not in water purification.
In addition to preferential flow pathways, humidity's nonlinear effect on the soil hydraulic conductivity becomes negligible only when soil water content surpasses a certain threshold (Regalado and Ritter, 2005).Nevertheless, it is essential to note that a high soil water content is not always a safeguard against water repellency since certain soil types, including Andosols, can retain their hydrophobic properties even when moistened (Wallis and Horne, 1992;Regalado and Ritter, 2008;Jordán et al., 2013).

Landscape Concerns
As shown in the previous sections, sites with sporadic fires in these temperate neotropical forests of the avocado belt will not face major degradation problems, especially on surface fires, where scattered water repellency does not imply infiltration reductions at landscape levels (Doerr et al., 2000).Moreover, water repellency could enhance water storage in the soil by diminishing direct soil evaporation rates (Ruthrof et al., 2019).So, we suggest that inherited water repellency in the soil can be part of the dynamics of a healthy forest ecosystem when fire regimes do not exceed resilience thresholds on the soil.
Water repellency on soils is a common phenomenon that could be found even on reference sites (Weninger et al., 2019).
The reference forest of this work could have water repellency even during the rainy seasons (Rodríguez Rodríguez, 2014).Activities like resining pine trees could increase the water repellency of the soils, decreasing their water infiltration potential; so recurrent low-severity fires could ameliorate this situation when they increase the patchiness of the soil water repellency distribution as long as they do not exceed land degradation thresholds.
However, sporadic fires tend to be an exception more than a rule, many of the fires in the forests of this region are associated with agricultural activities that are rapidly increasing their surface area (Barsimantov and Navia Antezana, 2012;Tejera Hernández et al., 2013;Olivares-Martinez et al., 2023).In this sense, it is pertinent to study whether areas with a higher fire frequency, enhanced by activities such as the international avocado trade, could exceed the soil's resilience, and trigger degrading processes.

CONCLUSION
In this work, we ascertained the presence of latent combustion in two temperate neotropical forests.The infiltration capacity of the Andosols showed a high decrease after ground or surface lowseverity fires with high resistance to change in periods of up to 20 years.The IRDI soil water repellency parameter was sensible for fire histories with up to two recurrent fires.Parameters such as the percentage of charcoal, pore space structure, and water repellency were useful indicators for understanding the fire history of these soils.
Parameters such as moisture at minimum hydrophobicity w. min and area S under the hydrophobicity curve were found to be edaphic memory traits related to fire history in the Andosols of subtropical regions.Although reduced, the infiltration capacity of S(1), S(2), and G(2) plots is not diminished to such an extent that surface runoff is generated in an average rainfall.
funding was provided by the projects: CONACyT Project CONACyT/PDCPN2016/Proyecto3053 "Impactos y consecuencias del desarrollo de la franja aguacatera sobre aspectos hidrológicos, funcionales, genéticos y de biodiversidad en ecosistemas templados de México" (Impacts and consequences of avocado growing-belt development on hydrological, functional, genetic and biodiversity aspects in temperate ecosystems of Mexico), ICTI-Project PFCTI/ ICTI/2019/A/315 "Aplicación del conocimiento ecológico para favorecer la sustentabilidad del cultivo de aguacate en el estado de Michoacán: Aspectos de suelo, hidrológicos e interacciones bióticas" (Application of ecological knowledge to promote the sustainability of avocado crops in the state of Michoacán: Soil, hydrological aspects and biological interactions).
Copyright © 2023 Olivares-Martinez, Gomez-Tagle and Mataix-Solera.This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY).The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice.No use, distribution or reproduction is permitted which does not comply with these terms.

FIGURE 1 |
FIGURE 1 | Study sites inside the volcanic region of Central Mexico.

FIGURE 2 |
FIGURE 2 | Morphology comparison between plots with strong (B) and insignificant (A) presence of smoldering combustion.Note the galleries with root charcoal until 80 cm.

FIGURE 3 |
FIGURE 3 | Water repellency across different fire histories.IRDI as a dimensionless index, w.min as volumetric water content, a.105 and a.max as the estimated contact angles.The bars represent standard errors, and the letters indicate the HSD statistical groups.Codes in the horizontal axis according to Table1.

TABLE 1 |
Summary of environmental features of the study sites.

TABLE 2 |
Hydraulic conductivity and active porosity summary.
Mean values ± standard deviations, where n is the number of water infiltration essays.Ks is the hydraulic conductivity with Gardner's α parameter.Bolded letters indicate the HSD groups for each variable (log-normal distribution in the four cases), same letter means no significant differences.

TABLE 3 |
Significative ANOVA models between fire histories.