Meta‐analysis of field‐saturated hydraulic conductivity recovery following wildland fire: Applications for hydrologic model parameterization and resilience assessment

Hydrologic recovery after wildfire is critical for restoring the ecosystem services of protecting of human lives and infrastructure from hazards and delivering water supply of sufficient quality and quantity. Recovery of soil‐hydraulic properties, such as field‐saturated hydraulic conductivity (K_fs), is a key factor for assessing the duration of watershed‐scale flash flood and debris flow risks after wildfire. Despite the crucial role of K_fs in parameterizing numerical hydrologic models to predict the magnitude of postwildfire run‐off and erosion, existing quantitative relations to predict K_fs recovery with time since wildfire are lacking. Here, we conduct meta‐analyses of 5 datasets from the literature that measure or estimate K_fs with time since wildfire for longer than 3‐year duration. The meta‐analyses focus on fitting 2 quantitative relations (linear and non‐linear logistic) to explain trends in K_fs temporal recovery. The 2 relations adequately described temporal recovery except for 1 site where macropore flow dominated infiltration and K_fs recovery. This work also suggests that K_fs can have low hydrologic resistance (large postfire changes), and moderate to high hydrologic stability (recovery time relative to disturbance recurrence interval) and resilience (recovery of hydrologic function and provision of ecosystem services). Future K_fs relations could more explicitly incorporate processes such as soil‐water repellency, ground cover and soil structure regeneration, macropore recovery, and vegetation regrowth.     

Watershed-scale vegetation,water quantity,and water quality responses to wildfire in the southern Appalachian mountain region,United States

Wildfires are landscape scale disturbances that can significantly affect hydrologic processes such as runoff generation and sediment and nutrient transport to streams. In Fall 2016, multiple large drought-related wildfires burned forests across the southern Appalachian Mountains. Immediately after the fires, we identified and instrumented eight 28.4–344 ha watersheds (four burned and four unburned) to measure vegetation, soil, water quantity, and water quality responses over the following two years. Within burned watersheds, plots varied in burn severity with up to 100% tree mortality and soil O-horizon loss. Watershed scale high burn severity extent ranged from 5% to 65% of total watershed area. Water quantity and quality responses among burned watersheds were closely related to the high burn severity extent. Total water yield (Q) was up to 39% greater in burned watersheds than unburned reference watersheds. Total suspended solids (TSS) concentration during storm events were up to 168 times greater in samples collected from the most severely burned watershed than from a corresponding unburned reference watershed, suggesting that there was elevated risk of localized erosion and sedimentation of streams. NO₃-N concentration, export, and concentration dependence on streamflow were greater in burned watersheds and increased with increasing high burn severity extent. Mean NO₃-N concentration in the most severely burned watershed increased from 0.087 mg L⁻¹ in the first year to 0.363 mg L⁻¹ (+317%) in the second year. These results suggest that the 2016 wildfires degraded forest condition, increased Q, and had negative effects on water quality particularly during storm events.     

Measurement Method Has a Larger Impact Than Spatial Scale For Plot-Scale Field-Saturated Hydraulic Conductivity (Kfs) After Wildfire and Prescribed Fire in Forests

Wildfires raise risks of floods, debris flows, major geomorphologic and sedimentologic change, and water quality and quantity shifts. A principal control on the magnitude of these changes is field-saturated hydraulic conductivity (K_fs), which dictates surface runoff generation and is a key input into numerical models. This work synthesizes 73 K_fs datasets from the literature in the first year following fire at the plot scale (≤ 10 m²). A meta-analysis using a random effects analysis showed significant differences between burned and unburned K_fs. The reductions in K_fs after fire, expressed by the ratio of K_{fs Burned}/K_{fs Unburned}, were 0.46 (95% confidence interval of 0.31-0.70) combining wildfire and prescribed fire and 0.3 (95% confidence interval of 0.13-0.71) for wildfire. No significant differences for K_fs were observed between wildfire and prescribed fire or moderate and high fire severity. Both K_fs magnitude and variability depended more on measurement method than measurement support area at the plot scale, with methods applying head ≥0.5 cm producing larger estimates of K_fs. It is recommended that post-fire efforts to characterize K_fs for modeling or process-based interpretations use methods that reflect the dominant infiltration processes: tension infiltrometers and simulated rainfall methods when soil matrix flow dominates and ponded head methods when macropore flow is critical. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.     

Effect of storms during drought on post‐wildfire recovery of channel sediment dynamics and habitat in the southern California chaparral,USA

Current global warming projections suggest a possible increase in wildfire and drought, augmenting the need to understand how drought following wildfire affects the recovery of stream channels in relation to sediment dynamics. We investigated post‐wildfire geomorphic responses caused by storms during a prolonged drought following the 2013 Springs Fire in southern California (USA), using multi‐temporal terrestrial laser scanning and detailed field measurements. After the fire, a dry‐season dry‐ravel sediment pulse contributed sand and small gravel to hillslope‐channel margins in Big Sycamore Creek and its tributaries. A small storm in WY 2014 generated sufficient flow to mobilize a portion of the sediment derived from the dry‐ravel pulse and deposited the fine sediment in the channel, totaling ~0.60 m³/m of volume per unit length of channel. The sediment deposit buried step‐pool habitat structure and reduced roughness by over 90%. These changes altered sediment transport characteristics of the bed material present before and after the storm; the ratio of available to critical shear stress (τ_o/τ_c) increased by five times. Storms during WY 2015 contributed additional fine sediment from tributaries and lower hillslopes and hyperconcentrated flow transported and deposited additional sediment in the channel. Together these sources delivered sediment on the order of six times that in 2014, further increasing τ_o/τ_c. These storms during multi‐year drought following wildfire transformed channel dynamics. The increased sediment transport capacity persisted during the drought period characterized by the longer residence time of relatively fine‐grained post‐fire channel sedimentation. This contrasts with wetter years, when post‐fire sediment is transported from the fluvial system during the same season as the post‐fire sediment pulse. Results of this short‐term study highlight the complex and substantial effects of multi‐year drought on geomorphic responses following wildfire. These responses influence pool habitat that is critical to longer‐term post‐wildfire riparian ecosystem recovery. Copyright © 2017 John Wiley & Sons, Ltd.     

Streamflow and suspended sediment yield following the 2000 Bobcat fire,Colorado

Postfire runoff and erosion are a concern, and more data are needed on the effects of wildfire at the watershed‐scale, especially in the Colorado Front Range. The goal of this study was to characterize and compare the streamflow and suspended sediment yield response of two watersheds (Bobcat Gulch and Jug Gulch) after the 2000 Bobcat fire. Bobcat Gulch had several erosion control treatments applied after the fire, including aerial seeding, contour log felling, mulching, and straw wattles. Jug Gulch was partially seeded. Study objectives were to: (1) measure precipitation, streamflow, and sediment yields; (2) assess the effect of rainfall intensity on peak discharges, storm runoff, and sediment yields; (3) evaluate short‐term hydrologic recovery. Two months after the fire, a storm with a maximum 30 min rainfall intensity I₃₀ of 42 mm h^?1 generated a peak discharge of 3900 l s^?1 km^?2 in Bobcat Gulch. The same storm produced less than 5 l s^?1 km^?2 in Jug Gulch, due to less rainfall and the low watershed response. In the second summer, storms with, I₃₀ of 23 mm h^?1 and 32 mm h^?1 generated peak discharges of 1100 l s^?1 km^?2 and 1700 l s^?1 km^?2 in the treated and untreated watersheds respectively. Maximum water yield efficiencies were 10% and 17% respectively, but 18 of the 23 storms returned ≤2% of the rainfall as runoff, effectively obscuring interpretation of the erosion control treatments. I₃₀ explained 86% of the variability in peak discharges, 74% of the variability in storm runoff, and >80% of the variability in sediment yields. Maximum single‐storm sediment yields in the second summer were 370 kg ha^?1 in the treated watershed and 950 kg ha^?1 in the untreated watershed. Copyright © 2005 John Wiley & Sons, Ltd.     

Synthesis of soil‐hydraulic properties and infiltration timescales in wildfire‐affected soils

We collected soil‐hydraulic property data from the literature for wildfire‐affected soils, ash, and unburned soils. These data were used to calculate metrics and timescales of hydrologic response related to infiltration and surface runoff generation. Sorptivity (S) and wetting front potential (Ψ_f) were significantly different (lower) in burned soils compared with unburned soils, whereas field‐saturated hydraulic conductivity (K_fs) was not significantly different. The magnitude and duration of the influence of capillarity during infiltration was greatly reduced in burned soils, causing faster ponding times in response to rainfall. Ash had large values of S and K_fs but moderate values of Ψ_f, compared with unburned and burned soils, indicating ash has long ponding times in response to rainfall. The ratio of S²/K_fs was nearly constant (~100 mm) for unburned soils but more variable in burned soils, suggesting that unburned soils have a balance between gravity and capillarity contributions to infiltration that may depend on soil organic matter, whereas in burned soils the gravity contribution to infiltration is greater. Changes in S and K_fs in burned soils act synergistically to reduce infiltration and accelerate and amplify surface runoff generation. Synthesis of these findings identifies three key areas for future research. First, short timescales of capillary influences on infiltration indicate the need for better measurements of infiltration at times less than 1 min to accurately characterize S in burned soils. Second, using parameter values, such as Ψ_f, from unburned areas could produce substantial errors in hydrologic modeling when used without adjustment for wildfire effects, causing parameter compensation and resulting underestimation of K_fs. Third, more thorough measurement campaigns that capture soil‐structural changes, organic matter impacts, quantitative water repellency trends, and soil‐water content along with soil‐hydraulic properties could drive the development of better techniques for numerically simulating infiltration in burned areas.     

Delineation of hydrologically similar units in a watershed based on fuzzy classification of soil hydraulic properties

Evaluation of flow and transport processes in a watershed‐scale requires that the watershed be divided into homogenous spatial units referred to as hydrologically similar units (HSUs). Although a few discretization schemes are already in use, a universally acceptable method of obtaining HSUs is yet to emerge. In this study, we developed a fuzzy inference system (FIS) to classify the saturated hydraulic conductivity (K_s) and two water‐retention parameters α and n into fuzzy logic‐based soil hydrologic classes (FSHCs). Analysis of these classes showed that soil properties within an FSHC have less variability and those between two FSHCs have large variability. This result suggested that soils belonging to a specific FSHC may be more similar than those across different FSHCs and may be grouped together to represent an HSU. Soils within a specific hydrologic class were aggregated to delineate HSUs within the watershed. For the Dengei Pahad micro‐watershed (DPW), this approach showed five distinct regions representing a discretized zone having similar soil hydraulic properties. Application of this approach on a larger international database of soil hydraulic properties revealed that the developed hydrologic classes are quite comparable across different databases. The delineated HSUs based on these FSHCs were also better than the soil series map of the watershed in maintaining the soil heterogeneity of the watershed. Moreover, this new discretization scheme using the SWAT modelling environment showed better performance than the soil series‐based discretization approach. Copyright © 2010 John Wiley & Sons, Ltd.     

Hydrological resilience to forest fire in the subarctic Canadian shield

Understanding the role of forest fires on water budgets of subarctic Precambrian Shield catchments is important because of growing evidence that fire activity is increasing. Most research has focused on assessing impacts on individual landscape units, so it is unclear how changes manifest at the catchment scale enough to alter water budgets. The objective of this study was to determine the water budget impact of a forest fire that partially burned a ~450 km² subarctic Precambrian Shield basin. Water budget components were measured in a pair of catchments: one burnt and another unburnt. Burnt and unburnt areas had comparable net radiation, but thaw was deeper in burned areas. There were deeper snow packs in burns. Differences in streamflow between the catchments were within measurement uncertainty. Enhanced winter streamflow from the burned watershed was evident by icing growth at the streamflow gauge location, which was not observed in the unburned catchment. Wintertime water chemistry was also clearly elevated in dissolved organics, and organic-associated nutrients. Application of a framework to assess hydrological resilience of watersheds to wildfire reveal that watersheds with both high bedrock and open water fractions are more resilient to hydrological change after fire in the subarctic shield, and resilience decreases with increasingly climatically wet conditions. This suggests significant changes in runoff magnitude, timing and water chemistry of many Shield catchments following wildfire depend on pre-fire land cover distribution, the extent of the wildfire and climatic conditions that follow the fire.     

Impacts of successive wildfire on soil hydraulic properties: Implications for debris flow hazards and system resilience

Climate and land use changes have led to recent increases in fire size, severity, and/or frequency in many different geographic regions and ecozones. Most post-wildfire geomorphology studies focus on the impact of a single wildfire but changing wildfire regimes underscore the need to quantify the effects of repeated disturbance by wildfire and the subsequent impacts on system resilience. Here, we examine the impact of two successive wildfires on soil hydraulic properties and debris flow hazards. The 2004 Nuttall-Gibson Complex and the 2017 Frye Fire affected large portions of the Pinaleño Mountains in southern Arizona, creating a mosaic of burn severity patterns that allowed us to quantify differences in wildfire-induced hydrologic changes as a function of burn severity and recent fire history (i.e. burned in only the Frye Fire or burned in both fires). Field observations after the 2017 Frye Fire indicated debris flow activity in areas burned predominantly at low severity. Many of these areas, however, were also affected by the 2004 Nuttall-Gibson Complex, suggesting that the relatively short recovery time between the two wildfires may have played a role in the geomorphic response to the most recent wildfire. Field measurements of soil hydraulic properties suggest that soils burned at moderate severity in 2004 and low severity in 2017 have a lower infiltration capacity relative to those that remained unburned in 2004 and burned at low severity in 2017. Simulations of runoff demonstrate that measured differences in infiltration capacity between once- and twice-burned soils are sufficient in some cases to influence the rainfall intensities needed to initiate runoff generated debris flows. Results quantify the impact of wildfire history and burn severity on runoff and debris flow activity in a landscape affected by successive wildfires and provide insight into how the resilience of geomorphic systems may be affected by successive wildfires. © 2019 John Wiley & Sons, Ltd.     

Parameter estimation for multiple post-wildfire hydrologic models

Predictions of post-wildfire flooding and debris flows are needed, typically with short lead times. Measurements of soil-hydraulic properties necessary for model parameterization are, however, seldom available. This study quantified soil-hydraulic properties, soil-water retention, and selected soil physical properties within the perimeter of the 2017 Thomas Fire in California. The Thomas Fire burn scar produced catastrophic debris flows in January 2018, highlighting the need for improved prediction capability. Soil-hydraulic properties were also indirectly estimated using relations tied to soil-water retention. These measurements and estimates are examined in the context of parameterizing post-wildfire hydrologic models. Tension infiltrometer measurements showed significant decreases (p < .05) in field-saturated hydraulic conductivity (K_fs) and sorptivity (S) in burned areas relative to unburned areas. Wildfire effects on soil water-retention were dominated by significant decreases in saturated soil-water content (θ_S). The van Genuchten parameters α, N, and residual water content did not show significant wildfire effects. The impacts of the wildfire on hydraulic and physical soil properties were greatest in the top 1 cm, emphasizing that measurements of post-fire soil properties should focus on the near-surface. Reductions in K_fs, θ_s, and soil-water retention in burned soils were attributed to fire-induced decreases in soil structure evidenced by increases in dry bulk density. Sorptivity reductions in burned soils were attributed to increases in soil-water repellency. Rapid post-fire assessments of flash flood and debris flow hazards using physically-based hydrologic models are facilitated by similarities between K_fs, S, and the Green–Ampt wetting front potential (ψ_f) with measurements at other southern CA burned sites. We suggest that ratios of burned to unburned K_fs (0.37), S (0.36), and ψ_f (0.66) could be used to scale unburned values for model parameterization. Alternatively, typical burned values (K_fs = 20 mm hr⁻¹; S = 6 mm hr^−0.5; ψ_f = 1.6 mm) could be used for model parameterization.     

Hydraulic redistribution and hydrological controls on aspen transpiration and establishment in peatlands following wildfire

In the sub‐humid Western Boreal Plains of Alberta, where evapotranspiration often exceeds precipitation, trembling aspen (Populus tremuloides Michx.) uplands often depend on adjacent peatlands for water supply through hydraulic redistribution. Wildfire is common in the Boreal Plains, so the resilience of the transfer of water from peatlands to uplands through roots immediately following wildfire may have implications for aspen succession. The objective of this research was to characterize post‐fire peatland‐upland hydraulic connectivity and assess controls on aspen transpiration (as a measure of stress and productivity) among landscape topographic positions. In May 2011, a wildfire affected 90,000 ha of north central Alberta, including the Utikuma Region Study Area (URSA). Portions of an URSA glacio‐fluval outwash lake catchment were burned, which included forests and a small peatland. Within 1 year after the fire, aspen were found to be growing in both the interior and margins of this peatland. Across recovering land units, transpiration varied along a topographic gradient of upland midslope (0.42 mm hr^?1) > upland hilltop (0.29 mm hr^?1) > margin (0.23 mm hr^?1) > peatland (0.10 mm hr^?1); similar trends were observed with leaf area and stem heights. Although volumetric water content was below field capacity, P. tremuloides were sustained through roots present, likely before fire, in peatland margins through hydraulic redistribution. Evidence for this was observed through the analysis of oxygen (δ¹⁸O) and hydrogen (δ²H) isotopes where upland xylem and peat core signatures were ?10.0‰ and ?117.8‰ and ?9.2‰ and ?114.0‰, respectively. This research highlights the potential importance of hydraulic redistribution to forest sustainability and recovery, in which the continued delivery of water may result in the encroachment of aspen into peatlands. As such, we suggest that through altering ecosystem services, peatland margins following fire may be at risk to aspen colonization during succession.     

Clearcutting and pine planting effects on nutrient concentrations and export in two mixed use headwater streams: Upper Coastal Plain,Southeastern USA

Timber harvest temporarily increases water yield; however, relationships between hydrologic and nutrient chemistry changes have not been consistent. This study quantified the effects of forest harvesting and site preparation without fertilization and with modern best management practices on nutrient concentrations and yields in small headwater streams of the Southeastern Coastal Plain. We monitored two watershed pairs for 2 years prior to and 1 year following timber harvest and for 2 more years following site preparation and planting. Treatment watersheds were clearcut, and downstream portions of streamside management zones were thinned in Fall 2003. Site preparation (herbicide application and burning) and planting followed a year later. All operations followed 1999 Georgia forestry best management practices. Previously published research revealed a large increase in water yield following harvest. Nutrient concentrations varied significantly within and between monitoring periods, even in reference watersheds. Silvicultural activities had no discernible effect on phosphorus and ammonium concentrations; however, statistically significant increases in nitrate/nitrite (67–340 µg L⁻¹) and total nitrogen concentrations (100–400 µg L⁻¹) in treatment watersheds followed stand re‐establishment. Nutrient yields increased after timber harvest largely as a result of increased water yields, although increased nutrient yields were small relative to inter‐annual and inter‐watershed variability and variability. Annual water yield largely explained the variability in annual nitrogen and phosphorus export from reference and treatment streams (r² values from 0.65 to 0.98). High NO_x concentrations coming from an upstream agricultural area decreased 1600–1800 µg L⁻¹ over several hundred metres in the treatment streams by dilution, uptake or denitrification. Copyright © 2013 John Wiley & Sons, Ltd.     

Initial impacts of a wildfire on hydrology and suspended sediment and nutrient export in California chaparral watersheds

Stream and rainfall gauging and runoff sampling were used to determine changes in hydrology and export of nutrients and suspended sediment from a June 2004 wildfire that burned 3010 ha in chaparral coastal watersheds of the Santa Ynez Mountains, California. Precipitation during water year 2005 exceeded average precipitation by 200–260%. Burned watersheds had order of magnitude higher peak discharge compared with unburned watersheds but similar annual runoff. Suspended sediment export of 181 mt ha^?1 from a burned watershed was approximately ten times greater than from unburned watersheds. Ammonium export from burned watersheds largely occurred during the first three storms and was 32 times greater than from unburned watersheds. Nitrate, dissolved organic nitrogen, and phosphate export from burned watersheds increased by 5.5, 2.8, and 2.2 times, respectively, compared with unburned chaparral watersheds. Storm runoff and peak discharge increase in burned compared with unburned sites were greatest during early season storms when enhanced runoff occurred. As the winter progressed, closely spaced storms and above average precipitation reduced the fire‐related impacts that resulted in significant increases in annual post‐fire runoff and export in other studies in southern California chaparral. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydrologic response of upland catchments to wildfires

To which extent do wildfires affect runoff production, soil erosion and sediment transport in upland catchments? This transient effect is investigated here by combining data of long term precipitation, sediment yield and wildfire records with a fine resolution spatially distributed modeling approach to flow generation and surface erosion. The model accounts for changes in the structure and properties of soil and vegetation cover by combining the tube-flux approach to topographic watershed partition with a parsimonious parametrization of hydrologic processes. This model is used to predict hydrologic and sediment fluxes for nine small catchments in Saint Gabriel mountains of southern California under control (pre-fire) and altered (post-fire) conditions. Simulation runs using a 45 years record of hourly precipitation show the passage of fire to significantly modify catchment response to storms with a major effect on erosion and flood flows. The probability of occurrence of major floods in the post-fire season is shown to increase up to an order of magnitude under same precipitation conditions. Also, the expected anomaly of sediment yield can increase dramatically the desertification hazard in upland wildfire prone areas. One should further consider the role of firefloods produced by the combined occurrence of wildfires and storms as a fundamental source of non-stationarity in the assessment of hydrologic hazard.     

Watersheds dynamics following wildfires: Nonlinear feedbacks and implications on hydrologic responses

In recent years, wildfires in the western United States have occurred with increasing frequency and scale. Climate change scenarios in California predict prolonged periods of droughts with even greater potential for conditions amenable to wildfires. The Sierra Nevada Mountains provide 70% of water resources in California, yet how wildfires will impact watershed-scale hydrology is highly uncertain. In this work, we assess the impacts of wildfires perturbations on watershed hydrodynamics using a physically based integrated hydrologic model in a high-performance-computing framework. A representative Californian watershed, the Cosumnes River, is used to demonstrate how postwildfire conditions impact the water and energy balance. Results from the high-resolution model show counterintuitive feedbacks that occur following a wildfire and allow us to identify the regions most sensitive to wildfires conditions, as well as the hydrologic processes that are most affected. For example, whereas evapotranspiration generally decreases in the postfire simulations, some regions experience an increase due to changes in surface water run-off patterns in and near burn scars. Postfire conditions also yield greater winter snowpack and subsequently greater summer run-off as well as groundwater storage in the postfire simulations. Comparisons between dry and wet water years show that climate is the main factor controlling the timing at which some hydrologic processes occur (such as snow accumulation) whereas postwildfire changes to other metrics (such as streamflow) show seasonally dependent impacts primarily due to the timing of snowmelt, illustrative of the integrative nature of hydrologic processes across the Sierra Nevada-Central Valley interface.     

Synergistic effects of water repellency and macropore flow on the hydraulic conductivity of a burned forest soil,south‐east Australia

Research shows that water repellency is a key hydraulic property that results in reduced infiltration rates in burned soils. However, more work is required in order to link the hydrological behaviour of water repellent soils to observed runoff responses at the plot and hillslope scale. This study used 5 M ethanol and water in disc infiltrometers to quantify the role of macropore flow and water repellency on spatial and temporal infiltration patterns in a burned soil at plot (<10 m²) scale in a wet eucalypt forest in south‐east Australia. In the first summer and winter after wildfire, an average of 70% and 60%, respectively, of the plot area was water repellent and did not contribute to infiltration. Macropores (r > 0·5 mm), comprising just 5·5% of the soil volume, contributed to 70% and 95%, respectively, of the field‐saturated and ponded hydraulic conductivity (K_p). Because flow occurred almost entirely via macropores in non‐repellent areas, this meant that less than 2·5% of the soil surface effectively contributed to infiltration. The hydraulic conductivity increased by a factor of up to 2·5 as the hydraulic head increased from 0 to 5 mm. Due to the synergistic effect of macropore flow and water repellency, the coefficient of variation (CV) in K_p was three times higher in the water‐repellent soil (CV = 175%) than under the simulated non‐repellent conditions (CV = 66%). The high spatial variability in K_p would act to reduce the effective infiltration rate during runoff generation at plot scale. Ponding, which tend to increase with increasing scale, activates flow through macropores and would raise the effective infiltration rates at larger scales. Field experiments designed to provide representative measurements of infiltration after fire in these systems must therefore consider both the inherent variability in hydraulic conductivity and the variability in infiltration caused by interactions between surface runoff and hydraulic conductivity. Copyright © 2010 John Wiley & Sons, Ltd.     

Infiltration and runoff generation processes in fire‐affected soils

Post‐wildfire runoff was investigated by combining field measurements and modelling of infiltration into fire‐affected soils to predict time‐to‐start of runoff and peak runoff rate at the plot scale (1 m²). Time series of soil‐water content, rainfall and runoff were measured on a hillslope burned by the 2010 Fourmile Canyon Fire west of Boulder, Colorado during cyclonic and convective rainstorms in the spring and summer of 2011. Some of the field measurements and measured soil physical properties were used to calibrate a one‐dimensional post‐wildfire numerical model, which was then used as a ‘virtual instrument’ to provide estimates of the saturated hydraulic conductivity and high‐resolution (1 mm) estimates of the soil‐water profile and water fluxes within the unsaturated zone. Field and model estimates of the wetting‐front depth indicated that post‐wildfire infiltration was on average confined to shallow depths less than 30 mm. Model estimates of the effective saturated hydraulic conductivity, K_s, near the soil surface ranged from 0.1 to 5.2 mm h^?1. Because of the relatively small values of K_s, the time‐to‐start of runoff (measured from the start of rainfall), t_p, was found to depend only on the initial soil‐water saturation deficit (predicted by the model) and a measured characteristic of the rainfall profile (referred to as the average rainfall acceleration, equal to the initial rate of change in rainfall intensity). An analytical model was developed from the combined results and explained 92–97% of the variance of t_p, and the numerical infiltration model explained 74–91% of the variance of the peak runoff rates. These results are from one burned site, but they strongly suggest that t_p in fire‐affected soils (which often have low values of K_s) is probably controlled more by the storm profile and the initial soil‐water saturation deficit than by soil hydraulic properties. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Performance of methods for estimating the time of concentration in a watershed of a tropical region

The time of concentration (T_c) is a fundamental parameter in the design of hydrological projects for watersheds. In this study a graphical methodology is described for estimating T_c in a watershed, and this is applied to 17 rainfall–runoff events from a rural watershed located near the capital city of Mato Grosso do Sul State, in the Brazilian Cerrado. The T_c values obtained through the graphical method were compared to T_c values estimated using 20 equations from various references. The equations were selected by considering those that were not developed using data for urban watersheds, and the results of the graphical method were compared to those derived by applying the equations to sub-basin data. The graphical method was reliable in determining T_c, and Ventura’s equation was found to present the best performance for a rural watershed in a tropical climate region.     

Effectiveness of introducing crop coefficient and leaf area index to enhance evapotranspiration simulations in hydrologic models

Evapotranspiration (ET) is an essential component of the hydrological cycle and plays a critical role in water resource management. However, ET is often overlooked in order to transform rainfall to runoff for better streamflow simulation. Hydrological models are commonly used to estimate areal actual evapotranspiration (AET) after calibration against observed discharge. However, classical approaches are often inadequate to appropriately simulate other hydrologic components. Hence, it is important to introduce natural heterogeneity to enhance hydrological processes and reduce water balance errors. In this study, the effectiveness of introducing a constant crop coefficient (K_c), flux tower‐based K_c, and leaf area index (LAI) to three hydrological models (Three‐Parametric Hydrologic Model [TPHM], Génie Rural à 4 paramètres Journalier [GR4J], and Catchment hydrologic cycle Assessment Tool [CAT]) is assessed for the simulation of daily streamflow and AET in a mountainous mixed forest watershed (8.54 km²) in South Korea. The results show that the streamflow simulations after introduction of K_c and LAI to the original models are quite similar. However, the effectiveness of the AET estimation was significantly enhanced after introduction of the flux tower‐based K_c and LAI. The information criterion computed to compare the models reveals that the flux tower‐based K_c‐simulated AET demonstrated better agreement with the observed AET. The Pearson's correlation coefficients (R²) of the TPHM (8%), GR4J (55%), and CAT (55%) models also showed improvements that were greater than the constant based K_c simulation. Similarly, the limitations of the three models with respect to capturing seasonal variation as well as high and low flows were enhanced after the introduction of the flux tower‐based K_c, which adequately reproduced hydrological processes with minimum water balance errors and bias. A regression analysis revealed the potential of estimating K_c as a linear function of LAI (R² = 0.84). The results of this study indicate that introduction of K_c and LAI to the conceptual rainfall–runoff models can be considered an effective approach to reduce water balance errors and uncertainties in hydrological models and improve the reliability of climate change studies and water resource management.     

Thresholds and relations for soil‐hydraulic and soil‐physical properties as a function of burn severity 4 years after the 2011 Las Conchas Fire,New Mexico,USA

Wildfire effects on soil‐physical and ‐hydraulic properties as a function of burn severity are poorly characterized, especially several years after wildfire. A stratified random sampling approach was used in 2015 to sample seven sites representing a spectrum of remotely sensed burn severity in the area impacted by the 2011 Las Conchas Fire in New Mexico, USA. Replicate samples from each site were analysed in the laboratory. Linear and linear indicator regression were used to assess thresholds in soil‐physical and ‐hydraulic properties and functional relations with remotely sensed burn severity. Significant thresholds were present for initial soil‐water content (θ_i) at 0–6 cm depth between the change in the Normalized Burn Ratio (dNBR) equal to 618–802, for bulk density (ρ_b) at 3–6 cm between dNBR equal to 416–533, for gravel fraction at 0–1 cm between dNBR equal to 416–533, for fines (the silt + clay fraction) at 0–1 cm for dNBR equal to 416–533, and for fines at 3–6 cm for dNBR equal to 293–416. Significant linear relations with dNBR were present between ρ_b at 0–1 cm, loss on ignition (LOI) at 0–1 cm, gravel fraction at 0–1 cm, and the large organic fraction at 1–3 cm. No thresholds or effects on soil‐hydraulic properties of field‐saturated hydraulic conductivity or sorptivity were observed. These results suggest that ρ_b and LOI at 0–1 cm have residual direct impacts from the wildfire heat impulse. The θ_i threshold is most likely from delayed groundcover/vegetation recovery that increases evaporation at the highest burn severity sites. Gravel and silt + clay thresholds at 0–1 cm at the transition to high burn severity suggest surface gravel lag development from hydraulic erosion. Thresholds in ρ_b from 3 to 6 cm and in silt + clay fraction from 3 to 6 cm appear to be the result of soil variability between sites rather than wildfire impacts. This work suggests that gravel‐rich soils may have increased resilience to sustained surface runoff generation and erosion following wildfire, with implications for assessments of postwildfire hydrologic and erosion recovery potential.