Next Article in Journal
Long-Term Changes and Variability of Ecologically-Based Climate Indices along an Altitudinal Gradient on the Qinghai-Tibetan Plateau
Next Article in Special Issue
Global Surface Temperature: A New Insight
Previous Article in Journal
Local Institutions and Climate Change Adaptation: Appraising Dysfunctional and Functional Roles of Local Institutions from the Bilate Basin Agropastoral Livelihood Zone of Sidama, Southern Ethiopia
Previous Article in Special Issue
Development of a Parametric Regional Multivariate Statistical Weather Generator for Risk Assessment Studies in Areas with Limited Data Availability
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Soil Response to Rainstorms after Wildfire and Prescribed Fire in Mediterranean Forests

1
Escuela Técnica Superior Ingenieros Agrónomos y Montes, Universidad de Castilla-La Mancha, Campus Universitario, E-02071 Albacete, Spain
2
Department AGRARIA, Mediterranean University of Reggio Calabria, Loc. Feo di Vito, I-89122 Reggio Calabria, Italy
*
Author to whom correspondence should be addressed.
Submission received: 17 November 2020 / Revised: 14 December 2020 / Accepted: 15 December 2020 / Published: 17 December 2020
(This article belongs to the Special Issue Application of Climatic Data in Hydrologic Models)

Abstract

:
The use of the Soil Conservation Service-curve number (SCS-CN) model for runoff predictions after rainstorms in fire-affected forests in the Mediterranean climate is quite scarce and limited to the watershed scale. To validate the applicability of this model in this environment, this study has evaluated the runoff prediction capacity of the SCS-CN model after storms at the plot scale in two pine forests of Central-Eastern Spain, affected by wildfire (with or without straw mulching) or prescribed fire and in unburned soils. The model performance has been compared to the predictions of linear regression equations between rainfall depth and runoff volume. The runoff volume was simulated with reliability by the linear regression only for the unburned soil (coefficient of Nash and Sutcliffe E = 0.73–0.89). Conversely, the SCS-CN model was more accurate for burned soils (E = 0.81–0.97), also when mulching was applied (E = 0.96). The performance of this model was very satisfactory in predicting the maximum runoff. Very low values of CNs and initial abstraction were required to predict the particular hydrology of the experimental areas. Moreover, the post-fire hydrological “window-of-disturbance” could be reproduced only by increasing the CN for the storms immediately after the wildfire. This study indicates that, in Mediterranean forests subject to the fire risk, the simple linear equations are feasible to predict runoff after low-intensity storms, while the SCS-CN model is advisable when runoff predictions are needed to control the flooding risk.

1. Introduction

The importance of forests in the climatic context at the planetary scale is well known, since forests produce oxygen and store carbon, regulate water and energy fluxes, support biodiversity and provide other fundamental ecosystem services [1,2,3]. However, the fundamental role of forests is threatened by some natural and anthropogenic agents, such as the extreme weather events and fire, with a long history of influence on forest ecosystems [4]. Extreme weather events (e.g., heavy storms and severe drought) are more and more intensified by climate change trends and occur in all regions of the world [5,6]. The fire effects extend to several components of forests, such as soil, vegetation, air and surface water [7]. With regard to the effects on surface water and soil, fire strongly alters the hydrological response of recently burnt areas, increasing by many folds the soil’s aptitude to generate runoff and erosion compared the unburned forest areas [8,9]. High runoff and erosion rates in forests increase the flooding risk, debris flow occurrence and water quality alteration in downstream areas, with possible loss of human lives and heavy damage to infrastructures and environment [10,11]. For instance, with regard to the Mediterranean climate, exceptionally high erosion rates (up to 100 tons per ha) have been reported immediately after wildfire by Menendez-Duarte et al. [12] in the Iberian peninsula, while Lopez-Batalla et al. [13] measured increases in flood runoff by 30% and in peak discharges by 120% in the same environments. These studies together with a large body of literature (see several examples in the milestone review of Shakesby [14]) clearly demonstrate the need to control and mitigate the hydrological response of forest soils after the wildfire adopting prediction tools and post-fire management actions.
The hydrological processes in forests are influenced by several factors, among which fire severity is important [15]. In other words, the more severe the fire, the more susceptible the soil to increases in runoff and erosion and worsening in water quality changes in the downstream ecosystems [7]. For instance, Lucas-Borja et al. [7] found increases in runoff and erosion in soil burned by wildfire by about 20% and even 200%, respectively, compared to unburned soils in Central-Eastern Spain.
The hydrology of burned forests is very complex, since it is the product of several factors, such as climate and edapho-climatic conditions, fire severity, soil, vegetation, morphology and land management after fire [15,16,17,18,19]. The needs to predict the impacts of fire on runoff and erosion and to implement the most effective actions for rehabilitation of fire-affected areas have increased the demand for hydrological models [16,20]. This demand is particularly important for forest managers working in the Mediterranean Basin, which is characterized by dry and hot summers followed by frequent and high-intensity rains immediately after the wildfire season [15,21]. In Mediterranean areas, increases in wildfire frequency and burned areas are expected under the forecasted climate scenarios [22,23].
In Mediterranean forests, the literature reports applications and verifications of hydrological models with different nature and complexity: from the simplest empirical models (such as the Soil Conservation Service (SCS)-curve number model to predict runoff, the universal soil loss equation, USLE, to simulate soil erosion) through the semiempirical models (e.g., the Morgan–Morgan–Finney model, MMF) until the most complex physically-based models (for instance, the Water Erosion Prediction Project, WEPP) or even the artificial neural networks [24,25]. Nonetheless, the empirical models are sometime more commonly used compared to the more complex models, mainly in data-poor environments (that is, in those situations with limited availability of parameter inputs) and for quick identification of sources of water, sediments and pollutants in forests [24,25,26].
Accurate predictions of surface runoff are fundamental to achieve reliable estimations of erosion rates and water quality parameters using hydrological models [27], particularly in forests subject to climate change and fire, since the latter factors play a large influence of the soil’s hydrological response [28,29]. The Soil Conservation Service (SCS)-curve number (CN) model (hereinafter “SCS-CN model”) is one of the most common methods for estimating the runoff volume generated by a rainstorm [30,31]. The popularity of this method is due to its simplicity, ease of use, widespread acceptance and large availability of input data [32]. Moreover, the SCS-CN model takes into account most of the factors that influence runoff generation, such as soil type, land use, hydrologic condition and antecedent moisture of the soil [33]. The model has also been incorporated as a hydrological submodel in several distributed rainfall–runoff models at the watershed scale (e.g., AnnAGNPS—annual agricultural non-point source pollution model, SWAT—soil and water assessment tool model and HEC-HMS—hydrologic engineering center-hydrologic modeling system model), supporting its robustness and popularity [33,34]). To date, there is no any alternative model that offers as many advantages as the SCS-CN model, which therefore is still commonly used in the large majority of environments and climatic conditions [35].
However, various studies conducted throughout the world on the applicability of the SCS-CN method have suggested a need for further improvement or overhauling of the model [32,36], since in some environments the method provides unsatisfactory predictions, particularly when the soil’s hydrological response does not follow the runoff generation mechanism by saturation-excess. Moreover, in spite of the large number of requisites, the model has been surprisingly little used for hydrological predictions in fire-affected areas, and the CN values are not completely known in burned areas [37]. The hydrological research has been mainly carried out at the watershed scale, where post-fire runoff has been predicted using the SCS-CN method incorporated in watershed-scale models (e.g., WILDCAT4 Flow Model [38] and FIRE HYDRO [39]). For instance, the SCS-CN model was used by Candela et al. [40], who analyzed the flood frequency curves for pre- and post-fire conditions, showing an increase in the average curve numbers and a decrease in the catchment time lag. Increases by 25 units in post-fire CNs were estimated by Soulis [37] in a small Greek watershed using pre-fire and post-fire rainfall–runoff datasets. A daily-constant CN in the SWAT model was used by Nunes et al. [41] to simulate the effects of soil water repellency on runoff from burnt hillslopes in a Mediterranean forest throughout three years after fire. It is therefore evident that the modeling experiences are scarce at the plot scale. At this scale, modeling of soil hydrology is less complex compared to the watershed scale, where the hydrological response to Mediterranean storms is further complicated by a combination and overlaying of several hydrological processes (e.g., water routing in the channel network, ponding and uneven soil properties) other than surface runoff generation. Furthermore, the studies about the hydrological effects of post-fire management on runoff in forests using the SCS-CN model are scarce.
Therefore, there is a need of further studies that must evaluate the runoff prediction capacity of the SCS-CN model in forest hillslopes or plots affected by fire of different severity—a fire parameter referring to the effects of wildfire on plant communities—in comparison with simpler models that estimate runoff directly from precipitation, such as the linear regression equations. In other words, is the SCS-CN model accurate to predict surface runoff in Mediterranean burned forests? Is it able to simulate post-fire hydrology with or without rehabilitation management actions (such as straw mulching) after a wildfire? When is its use convenient compared to a simpler linear regression between rainfall and runoff?
This study aims to reply to these questions, evaluating the hydrological performance of the SCS-CN model in two pine forests of Central-Eastern Spain affected by a wildfire and a prescribed fire, respectively, having different fire severity. More specifically, observations of rainfall–runoff patterns collected throughout one year in undisturbed soils (assumed as control) and in plots subject to prescribed fire/wildfire (the latter with or without a mulching treatment) are compared with the corresponding predictions of the SCS-CN model and linear rainfall–runoff regressions. The outcomes of this study help land managers to adopt strategies to control the hydrological effects of fire in Mediterranean forests.

2. Materials and Methods

2.1. Study Areas

Two experimental areas were selected in pine forests of the Province of Albacete, Castilla—La Mancha Region, Central Eastern Spain. The first area (Sierra de las Quebradas, municipality of Liétor) was affected by a wildfire in July 2016. The second forest (municipality of Lezuza) was subjected to a prescribed fire in March 2016, to reduce fuel loading and thus the potential risk and severity of subsequent fires (Figure 1). Prior to wildfire, the soil cover of the forest was mainly composed of plants, litter and stones with variable composition.
Both study areas have a semiarid Mediterranean climate, BSk according to the Köppen–Geiger classification [42]. The average annual rainfall and medium annual temperature are 282 (Liétor) and 450 (Lezuza) mm and 13.5 (Liétor) and 16 (Lezuza) °C, respectively (Spanish National Meteorological Agency, 1950–2016). According to historical data (1990–2014) provided by the Spanish Meteorological Agency, the maximum precipitation is concentrated in October (44.5 mm) and the minimum in May (39.6 mm); from June to September a hot and dry period (air relative humidity below 50%) occurs. The mean minimum temperature of the coldest month is −0.9 °C and the mean maximum temperature of the hottest month is nearly 32 °C.

2.1.1. Wildfire-Affected Forest (Liétor)

The experimental area of Liétor is located at an elevation between 520 and 770 m a.s.l. with W-SW and N aspect and mean slope of 15–20% (Table 1). Soils are classified as Inceptisols and Aridisols with sandy-loam texture [43].
The wildfire burned about 830 ha of the forestland (mainly Pinus halepensis Mill). The mean tree density of this forest was between 500 and 650 trees per hectare with height between 7 and 14 m. Rosmarinus officinalis L., Brachypodium retusum (Pers.) Beauv., Cistus clusii Dunal, Lavandula latifolia Medik., Thymus vulgaris L., Helichrysum stoechas (L.) Moench, Macrochloa tenacissima (L.) Kunth, Quercus coccifera L. and Plantago albicans L. are the shrub or herbaceous species of the forest. The wildfire, classified as high-severity fire by the local forest managers according to the methodology proposed by Vega et al. [44], determined a tree mortality of 100% (Table 1). Forest floor was about 3–5 cm deep. This forest floor, as happened also in Lezuza, was blanketed with decaying Pinus halepensis M. needles and twigs and other wood debris such as cones or branches coming from trees. In both sites, the forest floor was blanketed with decaying Pinus halepensis M. needles and twigs and other wood debris such as cones or branches coming from trees in both sites, Lezuza and Liétor. More information related to this suggestion is provided in Table 1.
In September 2016 (two months after the wildfire), a mulching treatment was carried out in some areas of the burned site. Barley straw was spread manually on soil at a depth of 3 cm and a rate of 0.2 kg m−2 (dry weight), following the dose proposed by different authors for forests of Northern Spain, to achieve a burned soil cover over 80% [45].

2.1.2. Forest Subjected to Prescribed Fire (Lezuza)

The forest area was selected in a relatively hilly area at an elevation from about 1010 to 1040 m a.s.l., with a 50-year old mixed plantation of Pinus halepensis and Pinus pinaster. The mean slope is around of 15% and the aspect is N-NE. Soils are classified as Alfisols with Xeralf Rhodoxeralf horizon with clay texture [43] (Table 1).
The tree density of this forest area was about 500 trees per hectare with a mean height of 6.40 m. The understory was dominated by Quercus faginea Lam. L., Quercus ilex subsp. ballota, Quercus coccifera L., Juniperus oxycedrus, Brachypodium retusum P. and Thymus sp. (Table 1). Forest floor depth was about 5–7 cm.

2.2. Description of Experimental Plots and Measurement of Runoff Volume

In the two selected forests, experimental plots were installed in unburned (control) and burned areas. Plots were randomly distributed in the experimental site in areas with the same morphological and ecological characteristics to ensure comparability.
More specifically, in the wildfire-affected forest of Liétor, eighteen rectangular plots (20-m long × 10-m wide) were installed at a distance between 200 and 500 m. Six plots were located in the forest outside of the burned site and assumed as a control. Twelve plots were instead located in the burned area, of which six were not treated, while six plots were mulched. In the forest subjected to the prescribed fire in Lezuza, twelve plots (4-m long and 2-m wide) were isolated, of which six were located in the unburned area and the other six in the burned site. The prescribed fire was carried out under controlled air conditions in the forests (wind speed of 14 km/h, air temperature of 14 °C and relative humidity of 63%), which are reference values for applying the prescribed fire as a forest protection measure. The upper and side borders of all plots in both areas were hydraulically isolated by geotextile fabric pounded into the ground, to prevent external water inputs. At the plot bottom, runoff was collected using a metal fence conveying the water into a pipe, which discharged to a plastic tank of 25 (Liétor) or 50 (Lezuza) liters. In these plots, immediately the runoff volumes were measured after each rainfall event throughout an observation period of about one year (Table 2).
A weather station (WatchDog 2000 Series model) with a tipping bucket rain gauge was placed in each study area to measure total daily precipitation, storm duration, air temperature and rain intensity during the study period. In the hourly rainfall series of the experimental database, two consecutive events were considered separate, if no rainfall was recorded for 6 h or more [46,47]. The mean rainfall intensity was the total rainfall divided by the storm duration. The main variables characterizing the observed events are reported in Table 2.

2.3. Outlines on the SCS-CN Model

This model, proposed by the Soil Conservation Service of the United States Department of Agriculture in 1972, hypothesizes that:
V P n = W S
being:
  • V = runoff volume (mm);
  • Pn = net precipitation (mm);
  • W = water volume stored into the soil (mm);
  • S = maximum water storage capacity of soil (mm).
Pn is the difference between the total precipitation (P) and the initial losses (Ia, as the water storage in the soil dips, interception, infiltration and evapotranspiration) prior to surface runoff. Ia is assumed to be proportional to S through a coefficient λ:
I a = λ S
S is given by:
S = 25.4 · ( 1000 C N 10 )
when the parameter CN is the so called “curve number”. The CN can be considered as the soil’s aptitude to produce runoff and is a function of the hydrological properties and conditions of soil, and land use. The CN varies between 0 and 100 (0 means that the soil does not produces runoff, 100 means that all the precipitation turns into surface runoff and then the hydrological losses are zero).
According to this model, the runoff volume V is:
V = ( P λ S ) 2 P + ( 1 λ ) S
To estimate CN in agroforest areas, the soil hydrological class, vegetation cover, hydrological condition (good, medium and poor) and cultivation practice and the antecedent moisture condition (AMC) of the soil must be determined.
The soil hydrological class (A to D) is related to the soil’s capability to produce runoff, on its turn due to the soil infiltration capacity. A low runoff production capability corresponds to the A soil hydrological class, while the highest runoff capability is typical of less permeable soils D.
The actual AMC of the soil subject to a rainfall/runoff event was estimated as a function of the total height of precipitation in the five days before the event in the two different conditions of crop dormancy or the growing season. On this regard, three AMCs are identified:
  • AMCI: dry condition and minimum surface runoff;
  • AMCII: average condition and surface runoff;
  • AMCIII: wet condition and maximum surface runoff.
The SCS guidelines make the CN values available in tables for a soil of given hydrological class and condition, vegetal cover, cultivation practice and average AMC (AMCII). The values of CNs related to AMCI (CNI) or AMCIII (CNIII) can be calculated through the following equations:
C N I = 4.2 C N I I 10 0.058 C N I I
C N I I I = 23 C N I I 10 + 0.13 C N I I
The parameter AMC takes into account the influence of the soil water content on the hydrological response of the soil to the rainstorm and distinguishes “dry” (AMCI), “average” (AMCII) and “wet” (AMCIII) conditions depending on the total rainfall height of the five days before the event.

2.4. Model Implementation

2.4.1. Linear Regression between Rainfall and Runoff

A linear regression model was established between the surface runoff volume (dependent variable) and the rainfall height (independent variable) for each event, as follows:
V = a P
where:
  • V = runoff volume (mm);
  • P = total precipitation (mm);
  • a = slope (-).
The intercept of this linear equation was forced to zero, in order to avoid runoff without any precipitation.

2.4.2. SCS-CN Model

The SCS-CN model was first applied considering the “default” input parameters, that is, the values of λ and CN derived from the SCS guidelines for woods (control plots) or pasture (burned plots) for the soil hydrologic group A of the experimental soils and AMC “I” for all the modeled events (since no or vey low precipitation was recorded in the antecedent five days). However, the runoff prediction capacity was totally unsatisfactory using default CNs, since very large errors between predictions and observations were achieved. Therefore, the SCS-CN model was adjusted by manual trials tuning both λ and CN parameters until the maximum coefficient of efficiency E (see below) was achieved using optimal λ and CN (Table 3).

2.5. Evaluation of Model Prediction Accuracy

The runoff simulations of the linear regression equation and SCS-CN model were analyzed for “goodness-of-fit” with the corresponding observations. More specifically, the observed and simulated runoff volumes were visually compared in scatterplots. Then, the main statistics and the indexes of goodness-of-fit commonly used in literature (e.g., [27,48,49,50]) were adopted (Table 4): (i) the maximum, minimum, mean and standard deviation of both the observed and simulated values; (ii) the coefficient of determination (r2); (iii) the coefficient of efficiency of Nash and Sutcliffe [51] (E); (iv) the root mean square error (RMSE) and (v) the coefficient of residual mass (CRM, also known as “percent bias”, PBIAS). Table 4 reports the equations and the range of variability of these indexes [52,53,54,55]. Generally speaking, these indexes are based on the analysis of the errors between simulations and predictions of the modeled hydrological variables.

3. Results and Discussion

3.1. Hydrological Characterization

3.1.1. Wildfire-Affected Forest (Liétor)

During the observation period, only nine events (total rainfall of 413 mm) produced surface runoff. For these events, precipitation height and mean intensity were in the range 11.6–93.8 mm and 0.98–28.0 mm/h, respectively. Expectedly, in the burned plots the runoff (on average 0.60 mm with a maximum value of 2.20 mm) was higher compared to the unburned soils (average of 0.03 mm and maximum of 0.08 mm). This may be due to the reduced infiltration and some combination of sealing, soil water repellency, loss of surface cover and decrease in soil aggregate stability, for the loss of organic matter [61].
In the mulched soil the mean and maximum runoff was 0.53 and 1.65 mm, respectively (Figure 2). These volumes were lower compared to the runoff generated in the burned plots. This shows the effectiveness of mulching as post-fire management technique to reduce the runoff generation capacity of the burned soils. These results confirm several other studies about the efficacy of mulch application, in order to control the hydrological response of soil after wildfires (e.g., [8,62,63,64,65,66]).
A sudden increase in the runoff generation capacity was evident in the first event after the fire (21 October 2016), presumably due to the ash release (that sealed the soil) and changes in the physicochemical properties (as the depletion in the organic matter content, which reduces the aggregate stability of the soil) [28,67]. A temporal reduction in runoff generation was found in burned soils (treated or not). This indicates a decrease in the hydrological response over time since the fire, also noticed by several authors in the early storms immediately after wildfire (e.g., [68,69,70]). The higher runoff is due to both the changes in soil hydrological properties and to the reduction of the vegetal cover after fire. As a matter of fact, the development of a water-repellent layer (also due to the ash released by fire) over the soil surface and the destruction of soil aggregates reduce water infiltration and thus increase runoff [71,72]. Over time, the shrub and herb vegetation quickly recovery, which decreases the runoff generation on the soil left bare by wildfire [73].

3.1.2. Forest Subjected to Prescribed Fire (Lezuza)

Sixteen storms (totaling 368 mm) produced runoff. The mean runoff from these rainfall events was 0.39 mm and the maximum was 0.69 mm. In the burned plots, the mean and maximum runoff volumes recorded were 0.40 and 0.75 mm, respectively (Figure 3). For few events, runoff from burned soils was lower compared to the control plots, while, for the majority of the monitored storms, the soil subjected to prescribed fire produced noticeably more runoff compared to the unburned plots. This waiving soil response to storms confirms the low impacts of low-intensity fires on the hydrological response of soils, already observed by several authors (e.g., [28,74]). This means that the prescribed fire has a limited potential to change the soil properties that drives the hydrological behavior, such as the repellency and infiltrability. However, as observed for the wildfire, attention should be paid to the first rainfall events occurring immediately after fire, when the removal of almost all the vegetal cover (including the litter) may leave the soil bare and thus exposed to rainfall erosivity.

3.2. Hydrological Modeling

3.2.1. Linear Regression

The simulation of runoff volumes gave satisfactory results for the unburned soil both in Lezuza and Liétor, as shown by the values of E (close or over 0.75) and r2 (>0.62) indexes and the closeness between predictions and observations (mean error of 10-15%). RMSE values (0.01, Lezuza, and 0.07, Liétor) were under the limit of acceptance of Table 4 (half std. dev.) only for the runoff measured in Lezuza (0.165), but not in Liétor (0.005). In general, the linear models tended to overestimate the runoff volumes in unburned and burned and mulched plots (CRM < 0), while underestimating the runoff in burned plots (CRM > 0; Figure 4). In the latter condition, the maximum runoff values were noticeably underestimated (difference between 25 and 40%; Table 5).
The performance of the linear regression Equation (7) was only acceptable but not satisfactory in burned soils (with or without treatment), because E (between 0.52, soil burned by wildfire, and 0.62, soil burned by wildfire and then mulched) was lower than the suggested limit of Table 4 and the differences between the maximum values of observations and predictions were over 20%; only the mean values were close each others (error <12%; Table 5 and Figure 5a,b); the values of RMSE, which was acceptable in Lezuza (0.12 vs. a limit of 0.15), were 0.62 (burned and untreated plots) and 0.46 (burned and mulched soils) and therefore were over the acceptance limit (Table 4). This limited performance is mainly due to inaccurate prediction of the most intense rainfall–runoff event (21 October 2016), immediately following the fire. Moreover, for the burned soils, the RMSE values were higher than 50% of the standard deviation of observed runoff and thus not satisfactory; for soils burned by wildfire, also the coefficient of determination was poor (r2 < 0.39). This unsatisfactory model performance is also visually shown by the large scattering of the simulations around the regression line (Figure 5a,b), which highlights a particular prediction inaccuracy for the first rainfall event (21 October 2016 recorded in Liétor). This inaccuracy is due to the soil changes induced by the wild fire (e.g., water repellency, decreases in infiltration and interception), which alters its hydrological response to rainstorms, but disappear some weeks or a few months after the fire [7,8,9,26]. Moreover, the linear equation was not able to simulate the variability of the hydrological processes with the precipitation, since the same runoff was observed for the same precipitation. This means that linear regressions are not able to simulate with reliability the surface runoff produced in burned plots, although these models may give an indication at least of the magnitude of the hydrological response of soils under different precipitation input and conditions.

3.2.2. SCS-CN Model

The predictions of runoff volume became more accurate for the majority of fire (in terms of severity) and soil conditions (Figure 6a,b). Compared to the linear regressions, the runoff predictions improved in the unburned plots of Lezuza (E = 0.87 and r2 = 0.92), but slightly worsened in all the plots of Liétor (E = 0.88 and r2 = 0.95; Table 6).
Conversely, the runoff was predicted with greater accuracy in burned soils (with or without treatment), as shown by E > 0.80—with peaks of 0.96-0.97 after wildfire—and r2 > 0.94, these indicators being noticeably over the acceptance limit (Table 4). The RMSE values were always lower than 50% of the observed standard deviation (Table 6). In general, the SCS-CN model always showed a runoff overestimation (see CRM < 0). It should be noticed that the mean runoff values were predicted with less accuracy compared to the linear regression equations (differences between predictions and observations over 20% in some cases), but the SCS-CN model was much more reliable in predicting the maximum runoff volumes (differences lower than 7% for soils burned by wildfire and 15% for plots subjected to prescribed fire; Table 6). This indicates that, when accurate predictions of runoff are required to control the flooding risk (linked to the highest runoff volumes), the SCS-CN should be preferred to the simpler linear regressions.
Some additional considerations about SCS-CN model application in the experimental conditions should be made.
First, very low values of CN and λ were provided in this study as input to the model, in order to predict with accuracy runoff after the two fire-severity conditions. This means that the water losses during and immediately after the rainfall (reflected by Ia and S, such as water storage in the soil dips, interception, infiltration and evapotranspiration) are very high and the storms produce very small runoff volumes. In more detail, the small CN simulates a large water storage capacity (S) of soil through the infiltration process and λ must be decreased even by three orders of magnitude to simulate the very low initial water losses, due to interception and evapo-transpiration.
Second, both for wildfire and prescribed fire, unrealistic input parameters are required to simulate such a minimal runoff generation capacity of these soils. As a matter of fact, values of 15–16 (after wildfire) or even 0.25–3 (after prescribed fire) for CN and 0.0001 for λ against common values over 30 for CN and 0.2 for λ are needed to fit the runoff predictions to the corresponding observations. This should be taken into account when the SCS-CN model must be implemented in soils having a small hydrological response.
Third, a unique CN value as input for the SCS-CN model is not able to reproduce the increase in runoff immediately after wildfire. The worsening of the hydrological response of the burned soil both after wildfire and prescribed fire has been shown by a number of studies (e.g., [9]) and particularly in Mediterranean forests (e.g., Keizer et al. [8], in a Portuguese eucalypt forest; Lucas-Borja et al. [7,62], in pine forests of Central Spain). This increase is mainly due to soil water repellency and vegetation cover removal due to fires, but these effects disappear some months after fire. In order to simulate the hydrological effects of a repellent and almost bare soil, it is necessary to increase the CNs in this so-called “window-of-disturbance” [15,75] up to values that may be noticeably high. Soil mulching “smoothes” the increase in the soil hydrological response and this requires a lower increase in CN values.

4. Conclusions

This study evaluated the runoff prediction capacity of the SCS-CN model in comparison with linear regression equations after storms in two pine forests of Central-Eastern Spain affected by wildfire (with or without a rehabilitation treatment using straw mulching) or prescribed fire.
The simulation of runoff volumes by the linear regression gave satisfactory results only for the unburned soils. Conversely, for the burned plots, the linear regressions failed in simulating the runoff with reliability. The SCS-CN model was instead accurate to predict the runoff volume particularly in burned soils, also when mulching was applied. Although the mean runoff was predicted with less accuracy compared to the linear equations, the model performance was very satisfactory in predicting the maximum volumes. Moreover, all the soil conditions (unburned, burned and burned and mulched) were simulated with reliability. To reproduce the peculiar hydrology of the experimental areas, very low values of CNs and initial abstraction were required, which may appear unrealistic; moreover, the post-fire hydrological window-of-disturbance could be reproduced only by increasing the CN for the storms occurring few months after wildfire.
The performances of the two tested models indicate that, in Mediterranean forests subject to the fire risk, the use of simple linear equations is suggested for predicting runoff generated by relatively low storms, while the SCS-CN model is more reliable and therefore advisable when runoff predictions are needed to control the flooding risk.
Overall, the study has confirmed the viability of the SCS-CN method to reproduce the complex hydrological response of unburned and burned (and treated or not) soils of Mediterranean forests. Although this assumption is limited to the experimental conditions, the results are encouraging towards larger applications of this method in other climatic and geomorphologic conditions. However, further modeling studies are needed, in order to explore the runoff prediction capacity of the model in fire-affected forests with different ecological and soil characteristics. These studies should also be enlarged from the plot to the watershed scale, using more complex hydrological models based on the SCS-CN method. Once validated in a wide range of environmental contexts, the use of these models may support the land managers to control runoff and erosion in mountain forests that are prone to both the wildfire and hydrogeological risks.

Author Contributions

Conceptualization, D.A.Z.; methodology, A.L., D.D., P.D., P.A.P.-A. and B.G.C.; formal analysis, G.B. and M.E.L.-B.; investigation, A.L., D.D., P.D., P.A.P.-A. and B.G.C.; resources, M.E.L.-B.; data curation, A.L., D.D., P.D., P.A.P.-A. and B.G.C.; writing—Original draft preparation, D.A.Z. and M.E.L.-B.; writing—review and editing, D.A.Z., G.B. and M.E.L.-B.; supervision, D.A.Z., G.B. and M.E.L.-B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

Pedro Antonio Plaza Álvarez was supported by a predoctoral fellowship from the Spanish Ministry of Education, Culture and Sport (FPU16/03296).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hou, Y.; Zhang, M.; Liu, S.; Sun, P.; Yin, L.; Yang, T.; Yide, L.; Qiang, L.; Wei, X. The hydrological impact of extreme weather-induced forest disturbances in a tropical experimental watershed in south China. Forests 2018, 9, 734. [Google Scholar] [CrossRef] [Green Version]
  2. Aryal, Y.; Zhu, J. Effect of watershed disturbance on seasonal hydrological drought: An improved double mass curve (IDMC) technique. J. Hydrol. 2020, 585, 124746. [Google Scholar] [CrossRef]
  3. Lindner, M.; Maroschek, M.; Netherer, S.; Kremer, A.; Barbati, A.; Garcia-Gonzalo, J.; Lexer, M.J. Climate change impacts, adaptive capacity, and vulnerability of European forest ecosystems. For. Ecol. Manag. 2010, 259, 698–709. [Google Scholar] [CrossRef]
  4. Flanagan, S.A.; Hurtt, G.C.; Fisk, J.P.; Sahajpal, R.; Hansen, M.C.; Dolan, K.A.; Sullivan, J.H.; Zhao, M. Potential Vegetation and Carbon Redistribution in Northern North America from Climate Change. Climate 2016, 4, 2. [Google Scholar] [CrossRef] [Green Version]
  5. Nunes, L.J.; Meireles, C.I.; Pinto Gomes, C.J.; Almeida Ribeiro, N.M. Forest Contribution to Climate Change Mitigation: Management Oriented to Carbon Capture and Storage. Climate 2020, 8, 21. [Google Scholar] [CrossRef] [Green Version]
  6. Niemeyer, R.J.; Bladon, K.D.; Woodsmith, R.D. Long term hydrologic recovery after wildfire and post-fire forest management in the interior Pacific Northwest. Hydrol. Process. 2020, 34, 1182–1197. [Google Scholar] [CrossRef]
  7. Lucas-Borja, M.E.; González-Romero, J.; Plaza-Álvarez, P.A.; Sagra, J.; Gómez, M.E.; Moya, D.; de Las Heras, J. The impact of straw mulching and salvage logging on post-fire runoff and soil erosion generation under Mediterranean climate conditions. Sci. Total Environ. 2019, 654, 441–451. [Google Scholar] [CrossRef]
  8. Keizer, J.J.; Silva, F.C.; Vieira, D.C.; González-Pelayo, O.; Campos, I.M.A.N.; Vieira, A.M.D.; Prats, S.A. The effectiveness of two contrasting mulch application rates to reduce post-fire erosion in a Portuguese eucalypt plantation. Catena 2018, 169, 21–30. [Google Scholar] [CrossRef]
  9. Wilson, C.; Kampf, S.K.; Wagenbrenner, J.W.; MacDonald, L.H. Rainfall thresholds for post-fire runoff and sediment delivery from plot to watershed scales. For. Ecol. Manag. 2018, 430, 346–356. [Google Scholar] [CrossRef]
  10. Kinoshita, A.M.; Chin, A.; Simon, G.L.; Briles, C.; Hogue, T.S.; O’Dowd, A.P.; Albornoz, A.U. Wildfire, water, and society: Toward integrative research in the “Anthropocene”. Anthropocene 2016, 16, 16–27. [Google Scholar] [CrossRef] [Green Version]
  11. Smith, H.G.; Sheridan, G.J.; Lane, P.N.; Nyman, P.; Haydon, S. Wildfire effects on water quality in forest catchments: A review with implications for water supply. J. Hydrol. 2011, 396, 170–192. [Google Scholar] [CrossRef]
  12. Menéndez-Duarte, R.; Fernández, S.; Soto, J. The application of 137Cs to post-fire erosion in north-west Spain. Geoderma 2009, 150, 54–63. [Google Scholar] [CrossRef]
  13. Lopéz, R.; Batalla, R.J. Análisis del comportamiento hidrológico de la cuenca mediterránea de Arbúbices antes y después de un incendio forestal. III Congreso Forestal, Español, Granada, Spain. 2001. Available online: ile:///C:/Users/MDPI/AppData/Local/Temp/3CFE01-084.pdf (accessed on 16 November 2020).
  14. Shakesby, R.A. Post-wildfire soil erosion in the Mediterranean: Review and future research directions. Earth Sci. Rev. 2011, 105, 71–100. [Google Scholar] [CrossRef]
  15. Zavala, L.M.M.; de Celis Silvia, R.; López, A.J. How wildfires affect soil properties. A brief review. Cuad. Investig. geográfica/Geogr. Res. Lett. 2014, 40, 311–331. [Google Scholar] [CrossRef] [Green Version]
  16. Moody, J.A.; Shakesby, R.A.; Robichaud, P.R.; Cannon, S.H.; Martin, D.A. Current research issues related to post-wildfire runoff and erosion processes. Earth Sci. Rev. 2013, 122, 10–37. [Google Scholar] [CrossRef]
  17. Nunes, J.P.; BernardJannin, L.; Rodriguez Blanco, M.L.; Santos, J.M.; Coelho, C.D.O.A.; Keizer, J.J. Hydrological and erosion processes in terraced fields: Observations from a humid Mediterranean region in northern Portugal. Land Degrad. Dev. 2018, 29, 596–606. [Google Scholar] [CrossRef]
  18. Nunes, J.P.; Naranjo Quintanilla, P.; Santos, J.M.; Serpa, D.; Carvalho-Santos, C.; Rocha, J.; Keesstra, S.D. Afforestation, subsequent forest fires and provision of hydrological services: A model-based analysis for a Mediterranean mountainous catchment. Land Degrad. Dev. 2018, 29, 776–788. [Google Scholar] [CrossRef]
  19. Nunes, J.P.; Doerr, S.H.; Sheridan, G.; Neris, J.; Santín, C.; Emelko, M.B.; Silins, U.; Robichaud, P.R.; Elliot, W.J.; Keizer, J. Assessing water contamination risk from vegetation fires: Challenges, opportunities and a framework for progress. Hydrol. Process. 2018, 32, 687–1684. [Google Scholar] [CrossRef] [Green Version]
  20. Filianoti, P.; Gurnari, L.; Zema, D.A.; Bombino, G.; Sinagra, M.; Tucciarelli, T. An evaluation matrix to compare computer hydrological models for flood predictions. Hydrology 2020, 7, 42. [Google Scholar] [CrossRef]
  21. Lucas-Borja, M.E.; Zema, D.A.; Carrà, B.G.; Cerdà, A.; Plaza-Alvarez, P.A.; Cózar, J.S.; de las Heras, J. Short-term changes in infiltration between straw mulched and non-mulched soils after wildfire in Mediterranean forest ecosystems. Ecol. Eng. 2018, 122, 27–31. [Google Scholar] [CrossRef] [Green Version]
  22. IPCC. IPCC’s 5th Assessment Report for Europe. 2013. Available online: http://www.ipcc.ch/pdf/assessment-report/ar5/wg2/WGIIAR5-Chap23_FINAL.pdf (accessed on 31 August 2018).
  23. Bedia, J.; Herrera, S.; Camia, A.; Moreno, J.M.; Gutiérrez, J.M. Forest fire danger projections in the Mediterranean using ENSEMBLES regional climate change scenarios. Clim. Chang. 2014, 122, 185. [Google Scholar] [CrossRef] [Green Version]
  24. Zema, D.A.; Lucas-Borja, M.E.; Fotia, L.; Rosaci, D.; Sarnè, G.M.; Zimbone, S.M. Predicting the hydrological response of a forest after wildfire and soil treatments using an Artificial Neural Network. Comput. Electron. Agric. 2020, 170, 105280. [Google Scholar] [CrossRef]
  25. Zema, D.A.; Nunes, J.P.; Lucas-Borja, M.E. Improvement of seasonal runoff and soil loss predictions by the MMF (Morgan-Morgan-Finney) model after wildfire and soil treatment in Mediterranean forest ecosystems. Catena 2020, 188, 104415. [Google Scholar] [CrossRef]
  26. Merritt, W.S.; Letcher, R.A.; Jakeman, A.J. A review of erosion and sediment transport models. Environ. Model. Softw. 2003, 18, 761–799. [Google Scholar] [CrossRef]
  27. Zema, D.A.; Lucas-Borja, M.E.; Carrà, B.G.; Denisi, P.; Rodrigues, V.A.; Ranzini, M.; Soriano, F.C.; de Cicco, V.; Zimbone, S.M. Simulating the hydrological response of a small tropical forest watershed (Mata Atlantica, Brazil) by the AnnAGNPS model. Sci. Total Environ. 2018, 636, 737–750. [Google Scholar] [CrossRef] [Green Version]
  28. Pereira, P.; Francos, M.; Brevik, E.C.; Ubeda, X.; Bogunovic, I. Post-fire soil management. Curr. Opin. Environ. Sci. Health 2018, 5, 26–32. [Google Scholar] [CrossRef]
  29. Lucas-Borja, M.E.; Carrà, B.G.; Nunes, J.P.; Bernard-Jannin, L.; Zema, D.A.; Zimbone, S.M. Impacts of land-use and climate changes on surface runoff in a tropical forest watershed (Brazil). Hydrol. Sci. J. 2020, 65, 1956–1973. [Google Scholar] [CrossRef]
  30. Mishra, S.K.; Singh, V.P. Another look at SCS-CN method. J. Hydrol. Eng. 1999, 4, 257–264. [Google Scholar] [CrossRef]
  31. Mishra, S.K.; Singh, V.P. Soil Conservation Service Curve NUMBER (SCS-CN) methodology (Vol. 42); Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013; p. 277. [Google Scholar]
  32. Suresh Babu, P.; Mishra, S.K. Improved SCS-CN–inspired model. J. Hydrol. Eng. 2012, 17, 1164–1172. [Google Scholar] [CrossRef]
  33. Michel, C.; Andreassian, V.; Perrin, C. Soil conservation service curve number method: How to mend a wrong soil moisture accounting procedure. Water Resour. Res. 2005, 41, W02011. [Google Scholar] [CrossRef]
  34. Jain, M.K.; Mishra, S.K.; Suresh Babu, P.; Venugopal, K. On the Ia–S relation of the SCS-CN method. Hydrol. Res. 2006, 37, 261–275. [Google Scholar] [CrossRef]
  35. Garen, D.; Moore, D.S. Curve number hydrology in water quality modeling: Uses, abuses, and future directions. J. Am. Water Resour. Assoc. 2005, 41, 377–388. [Google Scholar] [CrossRef]
  36. Ponce, V.M.; Hawkins, R.H. Runoff curve number: Has it reached maturity? J. Hydrol. Eng. 1996, 1, 11–19. [Google Scholar] [CrossRef]
  37. Soulis, K.X. Estimation of SCS Curve Number variation following forest fires. Hydrol. Sci. J. 2018, 63, 1332–1346. [Google Scholar] [CrossRef]
  38. Hawkins, R.H.; Greenberg, R.J. WILDCAT4 Flow Model; University of Arizona, School of Renewable Natural Resources: Tucson, AZ, USA, 1990. [Google Scholar]
  39. Cerrelli, G.A. Fire hydro, a simplified method for predicting peak discharges to assist in the design of flood protection measures for western wildfires. In Proceedings of the 2005 Watershed Management Conference—Managing Watersheds for Human and Natural Impacts: Engineering, Ecological, and Economic Challenges, Williamsburg, VA, USA, 19–22 July 2005; Moglen, G.E., Ed.; American Society of Civil Engineers: Alexandria, VA, USA, 2005; pp. 935–941. [Google Scholar]
  40. Candela, A.; Aronica, G.; Santoro, M. Effects of Forest Fires on Flood Frequency Curves in a Mediterranean Catchment/Effets d’incendies de forêt sur les courbes de fréquence de crue dans un bassin versant Méditerranéen. Hydrol. Sci. J. 2005, 50, 1–206. [Google Scholar] [CrossRef]
  41. Nunes, J.P.; Catarina Simões Vieira, D.; Keizer, J.J. Comparing simple and complex approaches to simulate the impacts of soil water repellency on runoff and erosion in burnt Mediterranean forest slopes. In Proceedings of the European Geosciences Union, General Assembly, EGUGA, Vienna, Austria, 23 April 2017; p. 10815. [Google Scholar]
  42. Kottek, M.; Grieser, J.; Beck, C.; Rudolf, B.; Rubel, F. World map of the Köppen-Geiger climate classification updated. Meteorol. Z. 2006, 15, 259–263. [Google Scholar] [CrossRef]
  43. Soil Survey Staff. Soil Taxonomy. Soil Taxonomy: A Basic System of Soil Classification for Making and Interpreting Soil Surveys. Agricultural Handbook 436, 2nd ed.; Natural Resources Conservation Service, USDA: Washington, DC, USA, 1999; p. 869. [Google Scholar]
  44. Vega, J.A.; Fontúrbel, T.; Merino, A.; Fernández, C.; Ferreiro, A.; Jiménez, E. Testing the ability of visual indicators of soil burn severity to reflect changes in soil chemical and microbial properties in pine forests and shrubland. Plant Soil. 2013, 369, 73–91. [Google Scholar] [CrossRef]
  45. Vega, J.A.; Fernández, C.; Fonturbel, T.; González-Prieto, S.; Jiménez, E. Testing the effects of straw mulching and herb seeding on soil erosion after fire in a gorse shrubland. Geoderma 2014, 223, 79–87. [Google Scholar] [CrossRef]
  46. Zema, D.A.; Labate, A.; Martino, D.; Zimbone, S.M. Comparing different infiltration methods of the HEC-HMS model: The case study of the Mésima Torrent (Southern Italy). Land Degrad. Dev. 2017, 28, 294–308. [Google Scholar] [CrossRef]
  47. Wischmeier, W.H.; Smith, D.D. Prediction Rainfall Erosion Losses; Handbook, No. 537; USDA: Washington, DC, USA, 1978. [Google Scholar]
  48. Willmott, C.J. Some comments on the evaluation of model performance. Bull. Am. Meteorol. Soc. 1982, 63, 1309–1313. [Google Scholar] [CrossRef] [Green Version]
  49. Legates, D.R.; McCabe, G.J. Evaluating the use of “goodness of fit” measures in hydrologic and hydroclimatic model validation. Water Resour. Res. 1999, 35, 233–241. [Google Scholar] [CrossRef]
  50. Loague, K.; Green, R.E. Statistical and graphical methods for evaluating solute transport models: Overview and application. J. Contam. Hydrol. 1991, 7, 51–73. [Google Scholar] [CrossRef]
  51. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models: Part I. A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  52. Zema, D.A.; Bingner, R.L.; Govers, G.; Licciardello, F.; Denisi, P.; Zimbone, S.M. Evaluation of runoff, peak flow and sediment yield for events simulated by the AnnAGNPS model in a Belgian agricultural watershed. Land Degrad. Dev. 2012, 23, 205–215. [Google Scholar] [CrossRef]
  53. Krause, P.; Boyle, D.P.; Base, F. Comparison of different efficiency criteria for hydrological model assessment. Adv. Geosci. 2005, 5, 89–97. [Google Scholar] [CrossRef] [Green Version]
  54. Moriasi, D.N.; Arnold, J.G.; Van Liew, M.W.; Bingner, R.L.; Harmel, R.D.; Veith, T.L. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans. ASABE 2007, 50, 885–900. [Google Scholar] [CrossRef]
  55. Van Liew, M.W.; Garbrecht, J. Hydrologic simulation of the Little Washita River experimental watershed using SWAT. J. Am. Water Resour. Assoc. 2003, 39, 413–426. [Google Scholar] [CrossRef]
  56. Santhi, C.; Arnold, J.G.; Williams, J.R.; Dugas, W.A.; Srinivasan, R.; Hauck, L.M. Validation of the SWAT model on a large river basin with point and nonpoint sources. J. Am. Water Resour. Assoc. 2001, 37, 1169–1188. [Google Scholar] [CrossRef]
  57. Van Liew, M.W.; Arnold, J.G.; Garbrecht, J.D. Hydrologic simulation on agri- cultural watersheds: Choosing between two models. Trans. ASABE 2003, 46, 1539–1551. [Google Scholar] [CrossRef]
  58. Vieira, D.C.S.; Prats, S.A.; Nunes, J.P.; Shakesby, R.A.; Coelho, C.O.A.; Keizer, J.J. Modelling runoff and erosion, and their mitigation, in burned Portuguese forest using the revised Morgan-Morgan-Finney model. For. Ecol. Manag. 2014, 314, 150–165. [Google Scholar] [CrossRef]
  59. Singh, J.; Knapp, H.V.; Demissie, M. Hydrologic modeling of the Iroquois River watershed using HSPF and SWAT. ISWS CR 2004-08. Champaign, Ill.: Illinois State Water Survey. 2004. Available online: http://www.sws.uiuc.edu/pubdoc/CR/ISWSCR2004-08.pdf (accessed on 14 February 2018).
  60. Gupta, H.V.; Sorooshian, S.; Yapo, P.O. Status of automatic calibration for hydrologic models: Comparison with multilevel expert calibration. J. Hydrol. Eng. 1999, 4, 135–143. [Google Scholar] [CrossRef]
  61. Neary, D.G.; Ryan, K.C.; DeBano, L.F. Wildland fire in ecosystems: Effects of fire on soils and water. Gen. Tech. Rep. RMRS-GTR-42-vol. 4. Ogden, UT: US Department of Agriculture, Forest Service, Rocky Mountain Research Station. 250 p. 42. 2005. Available online: https://www.fs.fed.us/rm/pubs/rmrs_gtr042_4.pdf (accessed on 16 November 2020).
  62. Lucas-Borja, M.E.; Plaza-Álvarez, P.A.; Gonzalez-Romero, J.; Sagra, J.; Alfaro-Sánchez, R.; Zema, D.A.; de Las Heras, J. Short-term effects of prescribed burning in Mediterranean pine plantations on surface runoff, soil erosion and water quality of runoff. Sci. Total Environ. 2019, 674, 615–622. [Google Scholar] [CrossRef]
  63. Prats, S.A.; González-Pelayo, Ó.; Silva, F.C.; Bokhorst, K.J.; Baartman, J.E.; Keizer, J.J. Post-fire soil erosion mitigation at the scale of swales using forest logging residues at a reduced application rate. Earth Surf. Process. Landf. 2019, 44, 2837–2848. [Google Scholar] [CrossRef]
  64. Prats, S.A.; Malvar, M.C.; Coelho, C.O.A.; Wagenbrenner, J.W. Hydrologic and erosion responses to compaction and added surface cover in post-fire logged areas: Isolating splash, interrill and rill erosion. J. Hydrol. 2019, 575, 408–419. [Google Scholar] [CrossRef]
  65. Lopes, A.R.; Prats, S.A.; Silva, F.C.; Keizer, J.J. Effects of ploughing and mulching on soil and organic matter losses after a wildfire in Central Portugal. Cuad. Investig. Geográfica 2020, 46, 303–318. [Google Scholar] [CrossRef] [Green Version]
  66. Robichaud, P.R.; Lewis, S.A.; Brown, R.E.; Bone, E.D.; Brooks, E.S. Evaluating post-wildfire logging-slash cover treatment to reduce hillslope erosion after salvage logging using ground measurements and remote sensing. Hydrol. Process. 2020, 34, 4431–4445. [Google Scholar] [CrossRef]
  67. Keesstra, S.; Wittenberg, L.; Maroulis, J.; Sambalino, F.; Malkinson, D.; Cerdà, A.; Pereira, P. The influence of fire history, plant species and post-fire management on soil water repellency in a Mediterranean catchment: The Mount Carmel range, Israel. Catena 2017, 149, 857–866. [Google Scholar] [CrossRef]
  68. de Dios Benavides-Solorio, J.; MacDonald, L.H. Measurement and prediction of post-fire erosion at the hillslope scale, Colorado Front Range. Int. J. Wildland Fire 2005, 14, 457–474. [Google Scholar] [CrossRef] [Green Version]
  69. DeBano, L.F.; Debano, L.F.; Neary, D.G.; Ffolliott, P.F. Fire Effects on Ecosystems; John Wiley & Sons: Hoboken, NJ, USA, 1998. [Google Scholar]
  70. MacDonald, L.H.; Sampson, R.; Brady, D.; Juarros, L.; Martin, D. Chapter 4. Predicting post-fire erosion and sedimentation risk on a landscape scale: A case study from Colorado. J. Sustain. For. 2000, 11, 57–87. [Google Scholar] [CrossRef]
  71. DeBano, L.F.; Mann, L.D.; Hamilton, D.A. Translocation of hydrophobic substances into soil by burning organic litter 1. Soil. Sci. Soc. Am. J. 1970, 34, 130–133. [Google Scholar] [CrossRef]
  72. Shakesby, R.A.; Doerr, S.H.; Walsh, R.P.D. The erosional impact of soil hydrophobicity: Current problems and future research directions. J. Hydrol. 2000, 231, 178–191. [Google Scholar] [CrossRef]
  73. Díaz-Delgado, R.; Pons, X. Spatial patterns of forest fires in Catalonia (NE of Spain) along the period 1975–1995: Analysis of vegetation recovery after fire. For. Ecol. Manag. 2001, 147, 67–74. [Google Scholar] [CrossRef]
  74. Keeley, J.E. Fire intensity, fire severity and burn severity: A brief review and suggested usage. Int. J. Wildland Fire 2009, 18, 116–126. [Google Scholar] [CrossRef]
  75. Prosser, I.P.; Williams, L. The effect of wildfire on runoff and erosion in native Eucalyptus forest. Hydrol. Process. 1998, 12, 251–265. [Google Scholar] [CrossRef]
Figure 1. Location and layout of forest plots subject to prescribed fire and wildfire and monitored for hydrological observations (Lezuza and Liétor, Castilla La Mancha, Spain). Geographic coordinates and map source: Lezuza X: 557588E, Y: 4306475N; Liétor: X: 600081 E, Y: 4262798 N (unburned area); X: 598358 E, Y: 4264032 N (burned area); Google Earth, last access on 6/15/2019).
Figure 1. Location and layout of forest plots subject to prescribed fire and wildfire and monitored for hydrological observations (Lezuza and Liétor, Castilla La Mancha, Spain). Geographic coordinates and map source: Lezuza X: 557588E, Y: 4306475N; Liétor: X: 600081 E, Y: 4262798 N (unburned area); X: 598358 E, Y: 4264032 N (burned area); Google Earth, last access on 6/15/2019).
Climate 08 00150 g001
Figure 2. Rainfall and runoff volumes observed in forest plots subject to prescribed fire (Lezuza, Castilla La Mancha, Spain). Vertical bars are the standard deviations.
Figure 2. Rainfall and runoff volumes observed in forest plots subject to prescribed fire (Lezuza, Castilla La Mancha, Spain). Vertical bars are the standard deviations.
Climate 08 00150 g002
Figure 3. Rainfall and runoff volumes observed in forest plots subject to wildfire (Liétor, Castilla La Mancha, Spain). Vertical bars are the standard deviations.
Figure 3. Rainfall and runoff volumes observed in forest plots subject to wildfire (Liétor, Castilla La Mancha, Spain). Vertical bars are the standard deviations.
Climate 08 00150 g003
Figure 4. Linear regressions between observed rainfall and runoff in plots subject to prescribed fire and wildfire (Lezuza, Left, and Liétor, Right; Castilla La Mancha, Spain). V is the runoff volume and P is the rainfall height, while r2 and p are the coefficient of determination and the significance level, respectively.
Figure 4. Linear regressions between observed rainfall and runoff in plots subject to prescribed fire and wildfire (Lezuza, Left, and Liétor, Right; Castilla La Mancha, Spain). V is the runoff volume and P is the rainfall height, while r2 and p are the coefficient of determination and the significance level, respectively.
Climate 08 00150 g004
Figure 5. Scatterplots of runoff observations vs. predictions using linear regressions in plots subject to prescribed fire (Lezuza, a) and wildfire (Liétor, control soils, b, and burned soils c) (Castilla La Mancha, Spain).
Figure 5. Scatterplots of runoff observations vs. predictions using linear regressions in plots subject to prescribed fire (Lezuza, a) and wildfire (Liétor, control soils, b, and burned soils c) (Castilla La Mancha, Spain).
Climate 08 00150 g005
Figure 6. Scatterplots of runoff observations vs. predictions using the curve number model in plots subject to prescribed fire (Lezuza, a), and wildfire (Liétor, control soils, b, and burned soils, c) (Castilla La Mancha, Spain).
Figure 6. Scatterplots of runoff observations vs. predictions using the curve number model in plots subject to prescribed fire (Lezuza, a), and wildfire (Liétor, control soils, b, and burned soils, c) (Castilla La Mancha, Spain).
Climate 08 00150 g006
Table 1. Main characteristics of forests and plots subject to prescribed fire and wildfire (Lezuza and Liétor, Castilla La Mancha, Spain).
Table 1. Main characteristics of forests and plots subject to prescribed fire and wildfire (Lezuza and Liétor, Castilla La Mancha, Spain).
CharacteristicsPrescribed Fire (Lezuza)Wildfire (Liétor)
Soil Condition
ControlBurnedControlBurnedBurned and Mulched
Number of plots66666
Plot area (m2)8200
Elevation (m)1010–1040520–770
Slope (%)15 (±4.4)14.5 (±2.6)15–20
AspectNN-NEW-SW and N
Tree density (n ha−1)477 (±33)529 (±60)500–650
Tree height (m)8–137–14
Tree diameter (cm)27 (±7)32 (±6)25–35
Tree canopy cover (%)75–8560–70
Tree species compositionP. pinaster Ait., P. halepensis M.P. pinaster Ait., P. halepensis M., Q. ilex L.Pinus halepensis Mill
Shrub/herb speciesQuercus coccifera L., Brachypodium retusum P., Thymus vulgaris L., Dactylis glomerata L.Quercus coccifera L., Brachypodium retusum P., Thymus vulgaris L., Sanguisorba minor S.Rosmarinus officinalis L., Brachypodium retusum (Pers.) Beauv., Cistus clusii Dunal, Lavandula latifolia Medik., Thymus vulgaris L., Helichrysum stoechas (L.) Moench, Macrochloa tenacissima (L.) Kunth, Quercus coccifera L. and Plantago albicans L.
Canopy consumed by fire (%)-18 (±5)-90 (±15)
Shrub/herb cover (%)Pre-fire56 (±9)51 (±7)65 (±9)
Post-fire59 (±5)10 (±3)0 (±0)
Litter (%)Pre-fire40 (±8)44 (±7)10 (±5)
Post-fire38 (±9)79 (±6)5 (±2)
Bare soil (%)Pre-fire4 (±1)5 (±1)25 (±8)
Post-fire3 (±1)11 (±3)90 (±12)
Table 2. Main characteristics of precipitation events in plots subject to prescribed fire and wildfire (Lezuza and Liétor, Castilla La Mancha, Spain).
Table 2. Main characteristics of precipitation events in plots subject to prescribed fire and wildfire (Lezuza and Liétor, Castilla La Mancha, Spain).
EventDateDays after FireRainfall Height (mm)Maximum Intensity (mm h−1)
Prescribed fire (Lezuza)
14 Apr 2016520.08.8
26 May 20163720.18.4
318 May 20164910.24.3
412 Oct 201619621.18.8
519 Oct 201620327.45.3
68 Nov 201622317.05.6
72 Dec 201624752.44.2
823 Dec 201626859.611.6
911 Feb 201731838.26.3
104 Apr 201737720.25.7
1128 Apr 201739428.26.8
Wildfire (Liétor)
121 Oct 20169840.03.99
224 Nov 201612941.01.48
38 Dec 201614659.00.98
421 Dec 201615993.82.1
530 Jan 201719928.00.84
622 Feb 201722216.81.14
78 Mar 201723611.61.78
820 Mar 2017248102.616.2
912 May 201730120.73.77
Table 3. Optimal values of the Soil Conservation Service-curve number (SCS-CN) model parameters used for runoff predictions in plots subject to prescribed fire and wildfire (Lezuza and Liétor, Castilla La Mancha, Spain).
Table 3. Optimal values of the Soil Conservation Service-curve number (SCS-CN) model parameters used for runoff predictions in plots subject to prescribed fire and wildfire (Lezuza and Liétor, Castilla La Mancha, Spain).
Input ParameterSoil Condition
Prescribed Fire (Lezuza)Wildfire (Liétor)
UnburnedBurnedUnburnedBurnedBurned and Mulched
Soil hydrologic classA
λ0.0001
CN15160.253 (27) *3 (22) *
AMCI
Note: * indicates the CN value of the first modeled event.
Table 4. Indexes and related equations and range of variability to evaluate the runoff prediction capacity of the linear regression and curve number models in forest plots subject to prescribed fire and wildfire (Lezuza and Liétor, Castilla La Mancha, Spain).
Table 4. Indexes and related equations and range of variability to evaluate the runoff prediction capacity of the linear regression and curve number models in forest plots subject to prescribed fire and wildfire (Lezuza and Liétor, Castilla La Mancha, Spain).
IndexEquationRange of VariabilityAcceptance Limit and Notes
Coefficient of determination (r2) r 2 = [ i = 1 n ( O i O ¯ ) ( P i P ¯ ) i = 1 n ( O i O ¯ ) 2 i = 1 n ( P i P ¯ ) 2 ] 2 0 to 1>0.5 [56,57,58]
Coefficient of efficiency (E, Nash and Sutcliffe [51]) E = 1 i = 1 n ( O i P i ) 2 i = 1 n ( O i O ¯ ) 2 −∞ to 1“Good” model accuracy if E ≥ 0.75, “satisfactory” if 0.36 ≤ E ≤ 0.75 and “unsatisfactory” if E ≤ 0.36 [55]
Root mean square error (RMSE) RMSE = i = 1 n ( P i O i ) 2 n 0 to ∞<0.5 of observed standard deviation [59]
Coefficient of residual mass
(CRM or PBIAS, Loague and Green [50])
CRM = i = 1 n O i i = 1 n P i i = 1 n O i −∞ to ∞<0.25 [54]
CRM < 0 indicates model underestimation
CRM > 0 indicates model overestimation
[60]
Notes: n = number of observations; Oi, Pi = observed and predicted values at the time step i; O ¯ = mean of observed values.
Table 5. Statistics and indexes to evaluate the runoff prediction capacity of linear regression models in forest plots subject to prescribed fire and wildfire (Lezuza and Liétor, Castilla La Mancha, Spain).
Table 5. Statistics and indexes to evaluate the runoff prediction capacity of linear regression models in forest plots subject to prescribed fire and wildfire (Lezuza and Liétor, Castilla La Mancha, Spain).
Runoff VolumeMeanStandardMinimumMaximum r2ECRMRMSE
Prescribed fire (Lezuza)
Control
Observed0.390.250.140.69----
Simulated0.420.310.080.610.620.73-0.080.07
Burned
Observed0.400.160.190.75----
Simulated0.400.300.080.580.750.600.010.12
Wildfire (Liétor)
Control
Observed0.030.000.030.08----
Simulated0.030.010.020.070.900.89-0.070.01
Burned
Observed0.600.040.722.20----
Simulated0.590.150.431.320.220.520.020.62
Burned and mulched
Observed0.490.010.621.66----
Simulated0.550.140.401.230.390.62-0.110.46
Notes: r2 = coefficient of determination; E = coefficient of efficiency; CRM = coefficient of residual mass; RMSE = root mean square error.
Table 6. Statistics and indexes to evaluate the runoff prediction capacity of the curve number model in forest plots subject to prescribed fire and wildfire (Lezuza and Liétor, Castilla La Mancha, Spain).
Table 6. Statistics and indexes to evaluate the runoff prediction capacity of the curve number model in forest plots subject to prescribed fire and wildfire (Lezuza and Liétor, Castilla La Mancha, Spain).
Runoff Volume MeanStandardMinimumMaximumr2ECRMRMSE
Prescribed fire (Lezuza)
Control
Observed0.390.140.250.69----
Simulated0.370.150.200.730.920.870.060.05
Burned
Observed0.390.140.250.69----
Simulated0.410.180.210.790.950.81-0.030.06
Wildfire (Liétor)
Control
Observed0.030.030.0030.08----
Simulated0.020.030.0000.080.950.880.260.01
Burned
Observed0.600.720.042.20----
Simulated0.470.700.002.060.980.970.220.16
Burned and mulched
Observed0.490.620.011.66----
Simulated0.530.610.021.650.940.96-0.070.15
Notes: r2 = coefficient of determination; E = coefficient of efficiency; CRM = coefficient of residual mass; RMSE = root mean square error.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lucas-Borja, M.E.; Bombino, G.; Carrà, B.G.; D’Agostino, D.; Denisi, P.; Labate, A.; Plaza-Alvarez, P.A.; Zema, D.A. Modeling the Soil Response to Rainstorms after Wildfire and Prescribed Fire in Mediterranean Forests. Climate 2020, 8, 150. https://0-doi-org.brum.beds.ac.uk/10.3390/cli8120150

AMA Style

Lucas-Borja ME, Bombino G, Carrà BG, D’Agostino D, Denisi P, Labate A, Plaza-Alvarez PA, Zema DA. Modeling the Soil Response to Rainstorms after Wildfire and Prescribed Fire in Mediterranean Forests. Climate. 2020; 8(12):150. https://0-doi-org.brum.beds.ac.uk/10.3390/cli8120150

Chicago/Turabian Style

Lucas-Borja, Manuel Esteban, Giuseppe Bombino, Bruno Gianmarco Carrà, Daniela D’Agostino, Pietro Denisi, Antonino Labate, Pedro Antonio Plaza-Alvarez, and Demetrio Antonio Zema. 2020. "Modeling the Soil Response to Rainstorms after Wildfire and Prescribed Fire in Mediterranean Forests" Climate 8, no. 12: 150. https://0-doi-org.brum.beds.ac.uk/10.3390/cli8120150

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop