A GIS-BASED DISTRIBUTED SOIL EROSION MODEL： A CASE STUDY OF TYPICAL WATERSHED, SICHUAN BASIN

Based on the measuring data and Digital Elevation Data (DEM) in a typical watershed--Hemingguan Watershed, Nanbu County, Sichuan Province of China, a GIS-based distributed soil erosion model was developed particularly for the purple soil type. It takes 20 m × 20 m grid as calculating unit and operates at 10-minute time interval. The required input data to the model include DEM, soil, land use, and time-series of precipitation and evaporation loss. The model enables one to estimate runoff, erosion and sediment yield for each grid cell and route the flow along its flow path to the watershed outlet. Furthermore, the model is capable of calculating the total runoff; erosion and sediment yield for the entire watershed by recursion algorithm. The validation of the model demonstrated that it could quantitatively simulate the spatial distribution of hydrological variables in a watershed, such as runoff, vegetation entrapment, soil erosion, the degree of soil and water loss. Moreover, it can evaluate the effect of land use change on the runoff generation and soil erosion with an accuracy of 80% and 75% respectively. The application of this model to a neighboring watershed with similar conditions indicates that this distributed model could be extended to other similar regions in China.     

Discussion: Misrepresentation of the USLE in ‘Is sediment delivery a fallacy?’

The assertion that the application of the USLE to predicting soil losses within a catchment or watershed is not sound because the USLE provides an estimate of erosion that would be measured if the entire area were divided up into 22·1 m long plots, and the output from them all added together, is incorrect. The slope length factor was derived from data obtained using a wide range of plot lengths and included the 22·1 m length simply to force it to take on a value of 1·0 when the slope length is 22·1 m. The 22·1 m length has no physical significance but the USLE slope length factor has a physical basis when applied to planar and convex hillslopes. The use of sediment delivery ratios when the USLE is applied to concave areas attempts to correct for applying the USLE beyond its design criteria. It fails because, in using the sediment delivery ratios in the prediction sediment delivery, it is incorrectly assumed that sediment delivery ratios de not vary with the amount of sediment entering a zone of deposition.     

How does DEM resolution affect microtopographic characteristics,hydrologic connectivity,and modelling of hydrologic processes?

The resolution of a digital elevation model (DEM) is a crucial factor in watershed hydrologic and environmental modelling. DEM resolution can cause significant variability in the representation of surface topography, which further affects quantification of hydrologic connectivity and simulation of hydrologic processes. The objective of this study is to examine the effects of DEM resolution on (1) surface microtopographic characteristics, (2) hydrologic connectivity, and (3) the spatial and temporal variations of hydrologic processes. A puddle‐to‐puddle modelling system was utilized for surface delineation and modelling of the puddle‐to‐puddle overland flow dynamics, surface runoff, infiltration, and unsaturated flow for nine DEM resolution scenarios of a field plot surface. Comparisons of the nine modelling scenarios demonstrated that coarser DEM resolutions tended to eliminate topographic features, reduce surface depression storage, and strengthen hydrologic connectivity and surface runoff. We found that reduction in maximum depression storage and maximum ponding area was as high as 97.56% and 76.36%, respectively, as the DEM grid size increased from 2 to 80 cm. The paired t‐test and fractal analysis demonstrated the existence of a threshold DEM resolution (10 cm for the field plot), within which the DEM‐based hydrologic modelling was effective and acceptable. The effects of DEM resolution were further evaluated for a larger surface in the Prairie Pothole Region subjected to observed rainfall events. It was found that simulations based on coarser resolution DEMs (>10 m) tended to overestimate ponded areas and underestimate runoff discharge peaks. The simulated peak discharge from the Prairie Pothole Region surface reduced by approximately 50% as the DEM resolution changed from 2 to 90 m. Fractal analysis results elucidated scale dependency of hydrologic and topographic processes. In particular, scale analysis highlighted a unique constant–threshold–power relationship between DEM scale and topographic and hydrologic parameters/variables. Not only does this finding allow one to identify threshold DEM but also further develop functional relationships for scaling to achieve valid topographic characterization as well as effective and efficient hydrologic modelling. Copyright © 2016 John Wiley & Sons, Ltd.     

Spatial discretization of large watersheds and its influence on the estimation of hillslope sediment yield

The combined use of water erosion models and geographic information systems has facilitated soil loss estimation at the watershed scale. Tools such as the Geo‐spatial interface for the Water Erosion Prediction Project (GeoWEPP) model provide a convenient spatially distributed soil loss estimate but require discretization to identify hillslopes and channels. In GeoWEPP, the TOpographic PArameteriZation (TOPAZ) model is used as an automated procedure to extract a watershed boundary, hillslopes and channels from a digital elevation model (DEM). Previous studies in small watersheds have shown that the size of the hillslopes and the channel distribution affect the model estimates, but in large watersheds, the effects on the soil loss estimates have yet to be tested. Therefore, the objective of this study was to evaluate the effect of discretization on the hillslope sediment yield estimates using GeoWEPP in two large watersheds (>10 km²). The watersheds were selected and discretized varying the TOPAZ parameters [critical source area (CSA) and minimum source channel length (MSCL)] in a 30‐m resolution digital elevation model. The drainage networks built with TOPAZ were compared with each other using the drainage density index. The results showed that the discretization affected hillslope sediment yield estimates and their spatial distribution more than the total runoff. The drainage density index and the hillslope sediment yield were proportional but inversely related; thus, soil loss estimates were highly affected by the spatial discretization. As a result of this analysis, a method to choose the CSA and MSCL values that generates the greatest fraction of hillslopes having profile lengths less than 200 m was developed. This slope length condition is particularly crucial when using the WEPP and GeoWEPP models, in order for them to produce realistic estimates of sheet and rill erosion. Finally, and as a result of this analysis, a more reliable method was developed for selecting the TOPAZ channel network parameters (CSA and MSCL). Copyright © 2015 John Wiley & Sons, Ltd.     

How long is a hillslope?

Hillslope length is a fundamental attribute of landscapes, intrinsically linked to drainage density, landslide hazard, biogeochemical cycling and hillslope sediment transport. Existing methods to estimate catchment average hillslope lengths include inversion of drainage density or identification of a break in slope–area scaling, where the hillslope domain transitions into the fluvial domain. Here we implement a technique which models flow from point sources on hilltops across pixels in a digital elevation model (DEM), based on flow directions calculated using pixel aspect, until reaching the channel network, defined using recently developed channel extraction algorithms. Through comparisons between these measurement techniques, we show that estimating hillslope length from plots of topographic slope versus drainage area, or by inverting measures of drainage density, systematically underestimates hillslope length. In addition, hillslope lengths estimated by slope–area scaling breaks show large variations between catchments of similar morphology and area. We then use hillslope length–relief structure of landscapes to explore nature of sediment flux operating on a landscape. Distinct topographic forms are predicted for end‐member sediment flux laws which constrain sediment transport on hillslopes as being linearly or nonlinearly dependent on hillslope gradient. Because our method extracts hillslope profiles originating from every ridgetop pixel in a DEM, we show that the resulting population of hillslope length–relief measurements can be used to differentiate between linear and nonlinear sediment transport laws in soil mantled landscapes. We find that across a broad range of sites across the continental United States, topography is consistent with a sediment flux law in which transport is nonlinearly proportional to topographic gradient. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Lumped surface and sub‐surface runoff for erosion modeling within a small hilly watershed in northern Vietnam

Developing models to predict on‐site soil erosion and off‐site sediment transport at the agricultural watershed scale represent an on‐going challenge in research today. This study attempts to simulate the daily discharge and sediment loss using a distributed model that combines surface and sub‐surface runoffs in a small hilly watershed (< 1 km²). The semi‐quantitative model, Predict and Localize Erosion and Runoff (PLER), integrates the Manning–Strickler equation to simulate runoff and the Griffith University Erosion System Template equation to simulate soil detachment, sediment storage and soil loss based on a map resolution of 30 m × 30 m and over a daily time interval. By using a basic input data set and only two calibration coefficients based, respectively, on water velocity and soil detachment, the PLER model is easily applicable to different agricultural scenarios. The results indicate appropriate model performance and a high correlation between measured and predicted data with both Nash–Sutcliffe efficiency (Ef) and correlation coefficient (r²) having values > 0.9. With the simple input data needs, PLER model is a useful tool for daily runoff and soil erosion modeling in small hilly watersheds in humid tropical areas. Copyright © 2013 John Wiley & Sons, Ltd.     

A comparison of hillslope drainage area estimation methods using high-resolution DEMs with implications for topographic studies of gullies

Topographic models provide a useful tool for understanding gully occurrence in the landscape but require reliable estimates of gully head drainage areas. Modern high-resolution topography data (collected using structure from motion photogrammetry or light detection and ranging) is increasingly used for topographic studies of gullies, but little work has been done to assess the variability of gully head drainage area estimates using different methods. This study evaluated alternative approaches to using high-resolution digital elevation models (DEMs) so that gully topographic models can be more readily applied to any area with suitably high-resolution data. Specifically, we investigated the impact of single- or multiple-direction flow routing algorithms, DEM hydrologic-enforcement procedures and spatial resolution on gully head drainage area estimation. We tested these methods on a 40 km² site centred on Weany Creek, a low-relief semi-arid landscape draining towards the Great Barrier Reef, Australia. Using a subroutine to separate gully heads into those with divergent or convergent flow patterns upslope, we found that divergent flow conditions occurred at half of 484 studied gullies. Drainage areas estimated by different flow routing algorithms were more variable in these divergent cases than for convergent cases. This variation caused a significant difference between topographic threshold parameters (slope b and intercept k) derived from single- or multiple-direction flow routing algorithms, respectively. Different methods of hydrologic enforcement (filling or breaching) also affected threshold analysis, resulting in estimates of the exponent b being ~188% higher if the DEM was filled than if breached. The testing of the methods to date indicates that a finer resolution (≤2 m) DEM and a multiple-direction flow routing algorithm achieve the most realistic drainage area estimates in low-relief landscapes. For Weany Creek we estimated threshold parameters k = 0.033 and b = 0.189, indicating that it is highly susceptible to gully erosion.     

Water flux and sediment transport within a forested landscape: the role of connectivity,subsurface flow,and slope length scale on transport mechanism

The connectivity and upscaling of overland runoff and sediment transport are important issues in hillslope hydrology to identify water flux and sediment transport within landscape. These processes are highly variable in time and space with regard to their interactions with vegetation and soil surface conditions. The generation of overland runoff and its spatial connectivity were examined along a slope to determine the variations in the transport mechanism of runoff and soil particles by rain splash and overland runoff. Field experiments were conducted by erosion plots on a steep hillslope at lengths of 5, 10, and 15 m. The overland runoff connectivity and flow transport distance decreased with the slope length, while spatial variability of infiltration increased significantly with the slope length. Observation of subsurface flow revealed that surface soil and litter layer could have important role in water transport. However, the surface soil water content and water flux transport along the slope was highly variable for different storm events; the variability was related to the complexity of the system, mainly by way of the initial wetness conditions and infiltration characteristics. Only net rain‐splashed soil was measurable, but examination of the water flux, overland runoff and sediment transport connectivity, characteristics of sheetwash, and the variability in spatial infiltration indicated an increase in the contribution of the rain splash transport mechanism along the slope. Copyright © 2013 John Wiley & Sons, Ltd.     

ESTIMATION OF SOIL EROSION IN A RESERVOIR WATERSHED USING ^137CS FALLOUT RADIONUCLIDE

Sedimentation from soil erosion is a critical reservoir watershed management issue. Due to the difficulty of field investigations, empirical formulas are commonly used to estimate the soil erosion rate. However, these estimations are often far from accurate. An effective alternative to estimating soil erosion is to analyze the spatial variation of 137Cs inventory in the soil. 137Cs can be adsorbed by the soil and is widely assumed to change its distribution only when disturbed by rainfall and human activities. Thus, 137Cs distributed in soils can be a useful environmental tracer to estimate soil erosion. In this study, the net soil loss estimate is 108,346 t/yr and the gross erosion and net erosion rates are 10.1 and 9 t/ha yr respectively. The sediment delivery ratio is therefore estimated to be 0.9 based on the two erosion rates. Because of the steep hillsides in the watershed, only 10% of the sediment yield stayed in the deposition sites and 90% was transported to the river as the sediment output. Soil erosion estimates from spatial variations of the 137Cs activity in the Baishi river watershed showed satisfactory accuracy when compared to sediment yield data. Using soil 137Cs concentrations is therefore a feasible method for estimating soil loss or deposition in Taiwan. Data sampling, analysis and result of this approach are given in this paper.     

‘Out with the Old?’ Why coarse spatial datasets are still useful for catchment‐scale investigations of sediment (dis)connectivity

The increasing popularity of remote sensing techniques has created numerous options for researchers seeking spatial datasets, especially digital elevation models (DEMs), for geomorphic investigations. This yields an important question regarding what DEM resolution is most appropriate when answering questions of geomorphic significance. The highest possible resolution is not always the best choice for a particular research aim, and DEM resolution should be tailored to fit both the scale of investigation and the simplicity/complexity of modelling processes applied to the dataset. We find that DEM resolution has a significant effect on a simple model of bed load sediment connectivity in the Lockyer Valley, Queensland. We apply a simple bed load transport threshold to catchment DEMs at three different resolutions – 1 m, 5 m, and 25 m. We find that using a 1 m resolution DEM generates numerous disconnections along tributary channel networks that underestimates the sediment contributing area, i.e. effective catchment area (ECA), of seven tributary basins of Lockyer Creek. Utilizing a coarser (lower‐resolution) DEM helps eliminate erroneous disconnections, but can reduce the detail of stream network definition. We find that the 25 m resolution DEM provides the best measure of ECA for comparing sediment connectivity between tributary catchments. The utility of simple models and coarse‐resolution datasets is important for undertaking large, catchment‐scale geomorphic investigations. As catchment‐scale investigations are becoming increasingly entwined with river management and rehabilitation efforts, scientists need not embrace an ‘out with the old’ philosophy. Simple models and coarse‐resolution datasets can help better integrate geomorphic research with management strategies and provide inexpensive and quick first‐order insights into catchment‐scale processes that can help focus future management efforts. Copyright © 2017 John Wiley & Sons, Ltd.     

Impact of grid resolution on the integrated and distributed response of a coupled surface–subsurface hydrological model for the des Anglais catchment,Quebec

Digital elevation models (DEMs) at different resolutions (180, 360, and 720 m) are used to examine the impact of different levels of landscape representation on the hydrological response of a 690‐km² catchment in southern Quebec. Frequency distributions of local slope, plan curvature, and drainage area are calculated for each grid size resolution. This landscape analysis reveals that DEM grid size significantly affects computed topographic attributes, which in turn explains some of the differences in the hydrological simulations. The simulations that are then carried out, using a coupled, process‐based model of surface and subsurface flow, examine the effects of grid size on both the integrated response of the catchment (discharge at the main outlet and at two internal points) and the distributed response (water table depth, surface saturation, and soil water storage). The results indicate that discharge volumes increase as the DEM is coarsened, and that coarser DEMs are also wetter overall in terms of water table depth and soil water storage. The reasons for these trends include an increase in the total drainage area of the catchment for larger DEM cell sizes, due to aggregation effects at the boundary cells of the catchment, and to a decrease in local slope and plan curvature variations, which in turn limits the capacity of the watershed to transmit water downslope and laterally. The results obtained also show that grid resolution effects are less pronounced during dry periods when soil moisture dynamics are mostly controlled by vertical fluxes of evaporation and percolation. Copyright © 2010 John Wiley & Sons, Ltd.     

Bedform effect on the reorganization of surface and subsurface grain size distribution in gravel bedded channels

Quantification of river bedform variability and complexity is important for sediment transport modeling as well as for characterization of river morphology. Alluvial bedforms are shown to exhibit highly nonlinear dynamics across a range of scales, affect local bed roughness, and vary with local hydraulic, hydrologic, and geomorphic properties. This paper examines sediment sorting on the crest and trough of gravel bedforms and relates it to bed elevation statistics. The data analysed here are the spatial and temporal series of bed elevation, grain size distribution of surface and subsurface bed materials, and sediment transport rates from flume experiments. We describe surface topography through bedform variability in height and wavelength and multiscale analysis of bed elevations as a function of discharge. We further relate bedform migration to preferential distribution of coarse and fine sediments on the troughs and crests, respectively, measuring directly surface and subsurface grain size distributions, and indirectly the small scale roughness variations as estimated from high resolution topographic scans.     

Distributed evaluation of the contribution of soil erosion to the sediment yield from a watershed

The correct determination of the sediment yield from a basin is of paramount importance in several hydraulic and environmental applications, such as the evaluation of the storage reduction of artificial reservoirs. However, due to the highly episodic nature of sediment supply and transport in many environments and to the extreme complexity of the processes involved, the evaluation of the sediment load in a river is still highly uncertain. When the time scale of interest is sufficiently long, and when the primary sediment source comes from distributed erosion in the watershed, the problem can be tackled in an indirect fashion, by computing the contribution to the annual suspended yield from soil erosion. In order to accomplish this task, we propose a distributed application of the widely used USLE formula. The formula is automatically applied along drainage networks derived from a digital elevation model and properly modified to take into account the presence of deposition zones in the watershed. © 1997 John Wiley & Sons, Ltd.     

Modelling soil erosion with a downscaled landscape evolution model

The measurement and prediction of soil erosion is important for understanding both natural and disturbed landscape systems. In particular numerical models of soil erosion are important tools for managing landscapes as well as understanding how they have evolved over time. Over the last 40 years a variety of methods have been used to determine rates of soil loss from a landscape and these can be loosely categorized into empirical and physically based models. Alternatively, physically based landscape evolution models (LEMs) have been developed that provide information on soil erosion rates at much longer decadal or centennial scales, over large spatial scales and examine how they may respond to environmental and climatic changes. Both soil erosion LEMs are interested in similar outcomes (landscape development and sediment delivery) yet have quite different methodologies and parameterizations. This paper applies a LEM (the CAESAR model) for the first time at time and space scales where soil erosion models have largely been used. It tests the ability of the LEM to predict soil erosion on a 30 m experimental plot on a trial rehabilitated landform in the Northern Territory, Australia. It then continues to discuss the synergies and differences between soil erosion and LEMs. The results demonstrate that once calibrated for the site hydrology, predicted suspended sediment and bedload yields from CAESAR show a close correspondence in both volume and timing of field measured data. The model also predicts, at decadal scales, sediment loads close to that of field measured data. Findings indicate that the small‐scale drainage network that forms within these erosion plots is an important control on the timing and magnitude of sediment delivery. Therefore, it is important to use models that can alter the DEM to reflect changing topography and drainage network as well as having a greater emphasis on channel processes. Copyright © 2012 John Wiley & Sons, Ltd. and Commonwealth of Australia     

Sediment delivery from hillslopes and the Universal Soil Loss Equation: some perceptions and misconceptions

The Universal Soil Loss Equation (USLE) or the revised USLE (RUSLE) are often used together with sediment delivery ratios in order to predict sediment delivery from hillslopes. In using sediment delivery ratios for this purpose, it is assumed that the sediment delivery ratio for a given hillslope does not vary with the amount of erosion occurring in the upslope area. This assumption is false. There is a perception that hillslope erosion is calculated on the basis that hillslopes are, in effect, simply divided into 22·1 m long segments. This perception fails to recognize the fact the inclusion of the 22·1 m length in the calculation has no physical significance but simply produces a value of 1·0 for the slope length factor when slopes have a length equal to that of the unit plot. There is a perception that the slope length factor is inappropriate because not all the dislodged sediment is discharged. This perception fails to recognize that the USLE and the RUSLE actually predict sediment yield from planar surfaces, not the total amount of soil material dislocated and removed some distance by erosion within an area. The application of the USLE/RUSLE to hillslopes also needs to take into account the fact that runoff may not be generated uniformly over that hillslope. This can be achieved by an equation for the slope length factor that takes account of spatial variations in upslope runoff on soil loss from a segment or grid cell. Several alternatives to the USLE event erosivity index have been proposed in order to predict event erosion better than can be achieved using the EI₃₀ index. Most ignore the consequences of changing the event erosivity index on the values for the soil, crop and soil conservation protection factors because there is a misconception that these factors are independent of one another. Copyright © 2007 John Wiley & Sons, Ltd.     

Development of a grid‐based version of the SWAT landscape model

Integrated river basin models should provide a spatially distributed representation of basin hydrology and transport processes to allow for spatially implementing specific management and conservation measures. To accomplish this, the Soil and Water Assessment Tool (SWAT) was modified by integrating a landscape routing model to simulate water flow across discretized routing units. This paper presents a grid‐based version of the SWAT landscape model that has been developed to enhance the spatial representation of hydrology and transport processes. The modified model uses a new flow separation index that considers topographic features and soil properties to capture channel and landscape flow processes related to specific landscape positions. The resulting model is spatially fully distributed and includes surface, lateral and groundwater fluxes in each grid cell of the watershed. Furthermore, it more closely represents the spatially heterogeneous distributed flow and transport processes in a watershed. The model was calibrated and validated for the Little River Watershed (LRW) near Tifton, Georgia (USA). Water balance simulations as well as the spatial distribution of surface runoff, subsurface flow and evapotranspiration are examined. Model results indicate that groundwater flow is the dominant landscape process in the LRW. Results are promising, and satisfactory output was obtained with the presented grid‐based SWAT landscape model. Nash–Sutcliffe model efficiencies for daily stream flow were 0.59 and 0.63 for calibration and validation periods, and the model reasonably simulates the impact of the landscape position on surface runoff, subsurface flow and evapotranspiration. Additional revision of the model will likely be necessary to adequately represent temporal variations of transport and flow processes in a watershed. Copyright © 2014 John Wiley & Sons, Ltd.     

Assessing the Sensitivity of a Model of Runoff Formation in the Ussuri River Basin

A physically based model of runoff formation with daily resolution has been developed for the upper part of the Ussuri basin with an area of 24400 km² based on ECOMAG hydrological modeling platform. Two versions of the hydrological model have been studied: (1) a crude version with the spatial schematization of the drainage area and river network based on DEM 1 × 1 km with the use of soil and landscape maps at a scale of 1: 2500000 and (2) a detailed version with DEM 80 × 80 m and soil and landscape maps of the scale of 1: 100000. Each version of the model has been tested for two variants of meteorological inputs: (1) meteorological forcing data (temperature, air humidity, precipitation) at eight weather stations and (2) with the involvement of additional data on precipitation collected at 15 gages in the basin. The model has been calibrated and validated over a 34-year period (1979–2012) with the use of runoff data for the Ussuri R. and its tributaries. The results of numerical experiments for assessing the sensitivity of model hydrological response to the spatial resolution of land surface characteristics and the density of precipitation gaging stations are discussed.     

How do sediment yields from post‐wildfire debris‐laden flows depend on terrain slope,soil burn severity class,and drainage basin area? Insights from airborne‐LiDAR change detection

We derived a high‐resolution, spatially continuous map of erosion and deposition associated with the debris‐laden flows triggered by the 2011 Las Conchas wildfire and subsequent rainstorms over a 197 km² area in New Mexico, USA. This map was produced using airborne‐LiDAR‐derived bare‐earth digital elevation models (DEMs) acquired approximately one year before and one year after the wildfire. Differencing of the pre‐wildfire and post‐wildfire‐and‐rainstorm bare‐earth DEMs yielded a DEM‐of‐difference (DoD) map that quantifies the magnitude of ground‐surface elevation changes due to erosion/deposition within each 1 m² pixel. We applied a 0.3 m threshold filter to our DoD to remove changes that could have been due to artifacts and/or imperfect georeferencing. The 0.3 m value for the threshold filter was chosen based on the stated accuracy of the LiDAR as well as a comparison of areas of significant topographic change mapped in aerial photographs with those predicted using a range of candidate threshold values for the DoD filter. We developed an automated procedure that accepts the DoD map as input and computes, for every pixel in the DEM, the net sediment volume exported through each pixel by colluvial and/or fluvial processes using a digital routing algorithm. An analysis of the resulting sediment volume map for the Las Conchas fire demonstrates that sediment volume is proportional to upstream contributing area. After normalized by contributing area, the average sediment yield (defined as the sediment volume divided by the contributing area) increases as a power‐law functions of the average terrain slope and soil burn severity class (SBSC) with exponents equal to approximately 1.5. Our analysis quantifies the relationships among sediment yield, average terrain slope, and average soil burn severity class at the watershed scale and should prove useful for predicting the geomorphic response of wildfire‐affected drainage basins. Copyright © 2014 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.     

Automated watershed subdivision for simulations using multi-objective optimization

A qualitative trial-and-error approach is commonly used to define watershed subdivisions through varying a single topographic threshold value. A methodology has been developed to quantitatively determine spatially variable threshold values using topography and a user-defined landscape reference layer. Optimization and topographic parameterization algorithms were integrated to create solutions that minimize the number of sub-watersheds and maximize the agreement between the discretized watershed and the reference layer. The system was evaluated using different reference datasets such as soil type, land management, and landscape form. Comparison of simulated results indicated that the scenario using land management as the reference layer yielded results closer to the scenario subdivided using a constant topographic threshold but with approximately 10 times more sub-catchments and therefore indicating customization of the watershed subdivision to the user-defined reference layer. The proposed optimization technology could be used in adequately applying watershed modeling technology in developing conservation practice implementation plans.