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.     

Depression threshold control proxy to improve HEC-HMS modeling of depression-dominated watersheds

ABSTRACT

Many hydrologic models utilize delineation results from traditional methods which create a hydrologically connected drainage system. In depression-dominated areas, topographic characteristics of depressions are vital to modeling unique hydrologic processes associated with puddle-to-puddle (P2P) filling-spilling dynamics. The objective of this study is to evaluate the impacts of the P2P processes and dynamic changes in contributing area on outlet discharge. To do so, an improved HEC-HMS model is developed by incorporating a depression threshold control proxy (DTCP) and an improved conceptual framework. The DTCP uses a storage–discharge function to simulate the P2P dynamics. The improved conceptual framework counteracts the effect of full hydrologic connectivity created by traditional delineation methods by introducing depressional and non-depressional areas to each sub-basin. Application of the improved HEC-HMS model demonstrated that it was capable of accurately simulating outlet discharge and providing the details on surface connectivity and depression storage.     

Land cover controls on depression‐focused recharge on the Oak Ridges Moraine,southern Ontario,Canada

The Oak Ridges Moraine (ORM) is a key hydrogeologic feature in southern Ontario. Previous research has emphasized the importance of depression‐focused recharge (DFR) for the timing and location of water recharge to the ORM's aquifers. However, the significance of DFR has not been empirically demonstrated, and the ORM's permeable surficial deposits imply that rainfall and snowmelt will largely recharge vertically rather than move laterally to topographic depressions. The exception may be during winter and spring, when concrete soil frost limits infiltration and encourages overland flow. The potential for DFR was examined for closed depressions under forest and agricultural land covers with similar soils and surficial geology. Air temperatures, precipitation, snow depth and water equivalent, soil water contents, soil freezing, and depression surface‐water levels were monitored during the winter and spring of 2012–2013 and 2013–2014. Recharge (R) was estimated at the crest and base of each depression using a 1‐dimensional water balance approach and surface‐applied Br^? tracing. Both forest and agricultural land covers experienced soil freezing; however, forest soils did not develop concrete frost. Conversely, agricultural fields saw concrete frost, overland flow, episodic ponding, and subsequent drainage of rain‐on‐snow and snowmelt inputs in open depressions. Recharge at the base of open depressions exceeded that in surrounding areas by an order of magnitude, suggesting that DFR is a significant hydrologic process during winter and spring under agricultural land cover on the ORM. Closed topographic depressions under agricultural land cover on the ORM crest may serve as critical recharge “hot spots” during winter and spring, and the ability of the unsaturated zone beneath these depressions to modify the chemistry of recharging water deserves further attention.     

Modelling surface‐water depression storage in a Prairie Pothole Region

In this study, the Precipitation‐Runoff Modelling System (PRMS) was used to simulate changes in surface‐water depression storage in the 1,126‐km² Upper Pipestem Creek basin located within the Prairie Pothole Region of North Dakota, USA. The Prairie Pothole Region is characterized by millions of small water bodies (or surface‐water depressions) that provide numerous ecosystem services and are considered an important contribution to the hydrologic cycle. The Upper Pipestem PRMS model was extracted from the U.S. Geological Survey's (USGS) National Hydrologic Model (NHM), developed to support consistent hydrologic modelling across the conterminous United States. The Geospatial Fabric database, created for the USGS NHM, contains hydrologic model parameter values derived from datasets that characterize the physical features of the entire conterminous United States for 109,951 hydrologic response units. Each hydrologic response unit in the Geospatial Fabric was parameterized using aggregated surface‐water depression area derived from the National Hydrography Dataset Plus, an integrated suite of application‐ready geospatial datasets. This paper presents a calibration strategy for the Upper Pipestem PRMS model that uses normalized lake elevation measurements to calibrate the parameters influencing simulated fractional surface‐water depression storage. Results indicate that inclusion of measurements that give an indication of the change in surface‐water depression storage in the calibration procedure resulted in accurate changes in surface‐water depression storage in the water balance. Regionalized parameterization of the USGS NHM will require a proxy for change in surface‐storage to accurately parameterize surface‐water depression storage within the USGS NHM.     

Thresholds for run‐off generation in a drained closed depression

Hydrological threshold behaviour has been observed across hillslopes and catchments with varying characteristics. Few studies, however, have evaluated rainfall–run‐off response in areas dominated by agricultural land use and artificial subsurface drainage. Hydrograph analysis was used to identify distinct hydrological events over a 9‐year period and examine rainfall characteristics, dynamic water storage, and surface and subsurface run‐off generation in a drained and farmed closed depression in north‐eastern Indiana, USA. Results showed that both surface flow and subsurface tile flow displayed a threshold relationship with the sum of rainfall amount and soil moisture deficit (SMD). Neither surface flow nor subsurface tile flow was observed unless rainfall amount exceeded the SMD. Timing of subsurface tile flow relative to soil moisture response on the shoulder slope of the depression indicated that the formation and drainage of perched water tables on depression hillslopes were likely the main mechanism that produced subsurface connectivity. Surface flow generation was delayed compared with subsurface tile flow during rainfall events due to differences in soil water storage along depression hillslopes and run‐off generation mechanisms. These findings highlight the substantial impact of subsurface tile drainage on the hydrology of closed depressions; the bottom of the depression, the wettest area prior to drainage installation, becomes the driest part of the depression after installation of subsurface drainage. Rapid connectivity of localized subsurface saturation zones during rainfall events is also greatly enhanced because of subsurface drainage. Thus, less fill is required to generate substantial spill. Understanding hydrologic processes in drained and farmed closed depressions is a critical first step in developing improved water and nutrient management strategies in this landscape.     

Effect of permafrost thaw on the dynamics of lakes recharged by ice‐jam floods: case study of Yukon Flats,Alaska

Large river floods are a key water source for many lakes in fluvial periglacial settings. Where permeable sediments occur, the distribution of permafrost may play an important role in the routing of floodwaters across a floodplain. This relationship is explored for lakes in the discontinuous permafrost of Yukon Flats, interior Alaska, using an analysis that integrates satellite‐derived gradients in water surface elevation, knowledge of hydrogeology, and hydrologic modelling. We observed gradients in water surface elevation between neighbouring lakes ranging from 0.001 to 0.004. These high gradients, despite a ubiquitous layer of continuous shallow gravel across the flats, are consistent with limited groundwater flow across lake basins resulting from the presence of permafrost. Permafrost impedes the propagation of floodwaters in the shallow subsurface and constrains transmission to ‘fill‐and‐spill’ over topographic depressions (surface sills), as we observed for the Twelvemile‐Buddy Lake pair following a May 2013 ice‐jam flood on the Yukon River. Model results indicate that permafrost table deepening of 1–11 m in gravel, depending on watershed geometry and subsurface properties, could shift important routing of floodwater to lakes from overland flow (fill‐and‐spill) to shallow groundwater flow (‘fill‐and‐seep’). Such a shift is possible in the next several hundred years of ground surface warming and may bring about more synchronous water level changes between neighbouring lakes following large flood events. This relationship offers a potentially useful tool, well suited to remote sensing, for identifying long‐term changes in shallow groundwater flow resulting from thawing of permafrost. Copyright © 2015 John Wiley & Sons, Ltd.     

Analysis of geomorphologic and hydrological characteristics in watershed saturated areas using topographic‐index threshold and geomorphology‐based runoff model

The objective of this paper is to investigate the variation of geomorphology and runoff characteristics in saturated areas under different partial contributing area (PCA) conditions. Geomorphologic information and hydrologic records from two mid‐size watersheds in northern Taiwan were selected for analysis. The PCA ratio in the watershed during a storm was assumed equal to the ratio of the surface‐flow volume to the direct runoff volume from measured hydrologic data. The extents of PCA regions were then determined by using a topographic‐index threshold. Consequently, the geomorphologic factors in saturated and unsaturated areas could be calculated using a digital elevation model, and these factors could then be linked to a geomorphology‐based IUH model for runoff simulation, which can consider both the surface‐ and subsurface‐flow processes in saturated and unsaturated areas, respectively. The results show that geomorphologic characteristics in the saturated areas vary significantly with different PCA ratios especially for higher order streams. A large PCA ratio results in a sharp hydrograph because the quick surface flow dominates the runoff process, whereas the hydrologic response in a low PCA case is dominated by the delayed subsurface flow. Copyright © 2007 John Wiley & Sons, Ltd.     

The bed morphology of upland single‐thread channels in semi‐arid environments: evidence of repeating bedforms and their wider implications for gravel‐bed rivers

Single‐thread, gravel‐bed streams of moderate slope in the northern Negev are characterized by three channel units: bars exhibit steeper than average slopes and poorly sorted mixtures of small–medium cobbles and coarse–very coarse pebbles; flats are associated with more gentle slopes and well‐sorted medium–fine pebbles and granules; and transitional units have intermediate slopes and grain size. In general, all three units are planar, span the full channel width and have well‐defined boundaries. Bars and flats are more common than the transitional units and alternate downstream for distances of several hundred metres, forming sequences that are reminiscent of the riffle–pool structure commonly observed in humid‐temperate gravel‐bed rivers. A notable contrast is the absence of significant bed relief: bars lack crests and flats lack depressions. The relative lack of bed relief in bar–flat sequences is attributed to the high rate of sediment supply from the sparsely vegetated hillslopes which promotes the infilling of depressions and to the erosion of crests under conditions of intense transport. This reduction of bed relief lowers channel roughness, which in turn increases flow velocity and, therefore, the ability of the channel to transmit the large sediment loads it receives. Although our analyses pertain to a semi‐arid river system, the results have wider implications for understanding the adjustment of channel bedform to high sediment loads in other fluvial environments. Copyright © 2012 John Wiley & Sons, Ltd.     

Delineating floodplain and upland areas for hydrologic models: a comparison of methods

A spatially distributed representation of basin hydrology and transport processes in hydrologic models facilitates the identification of critical source areas and the placement of management and conservation measures. Floodplains are critical landscape features that differ from neighbouring uplands in terms of their hydrological processes and functions. Accordingly, an important step in watershed modelling is the representation of floodplain and upland areas within a watershed. The aim of this study is (1) to evaluate four floodplain–upland delineation methods that use readily available topographic data (topographic wetness index, slope position, uniform flood stage, and variable flood stage) with regard to their suitability for hydrological models and (2) to introduce an evaluation scheme for the delineated landscape units. The methods are tested in three U.S. watersheds ranging in size from 334 to 629 km² with different climatic, hydrological, and geomorphological characteristics. Evaluation of the landscape delineation methods includes visual comparisons, error matrices (i.e. cross‐tabulations of delineated vs reference data), and geometric accuracy metrics. Reference data were obtained from the Soil Survey Geographic (SSURGO) database and Federal Emergency Management Agency (FEMA) flood maps. Results suggest that the slope position and the variable flood stage method work very well in all three watersheds. Overall percentages of floodplain and upland areas allocated correctly were obtained by comparing delineated and reference data. Values range from 83 to 93% for the slope position and from 80 to 95% for the variable flood stage method. Future studies will incorporate these two floodplain–upland delineation methods into the subwatershed‐based hydrologic model Soil and Water Assessment Tool (SWAT) to improve the representation of hydrological processes within floodplain and upland areas. Copyright © 2016 John Wiley & Sons, Ltd.     

Supraglacial melt channel networks in the Jakobshavn Isbræ region during the 2007 melt season

Supraglacial channels are an important mechanism for surface water transport over the ablation zone of western Greenland. The first assessment of the spatio‐temporal distribution of surface melt channels and their relationship to supraglacial lakes over the Jakobshavn Isbræ region of Western Greenland was analysed using Landsat Enhanced Thematic Mapper Plus panchromatic images during the 2007 melt season. A total of 1188 melt channels were delineated and show an increase in the number of melt channels throughout the season, reaching a peak on 9 August. Water‐filled melt channels advanced to a maximum elevation of 1647 m on 9 August and attained a minimum average slope of 0.009 on 8 July. The ablation zone demonstrates two hydrologic modes, where crevasse and moulin terminating channels dominate at elevations <800 m and higher‐order channel networks >800 m. Development of higher‐order networks is interrupted by flow divergence due to partitioning of melt water into vertical infiltration through moulins and crevasse fields prevalent at lower elevations. Tributary and connector networks between 800 and 1200 m in elevation are correlated with fewer lake occurrences, lower surface velocities (~50 m a^?1), and ice flow dominated by internal deformation over basal sliding. High‐order channels are associated with lake basins that exceed melt water storage capacity. Evolution of channel networks is coupled to changes in melt water production, runoff, and ice dynamics with implication for the englacial and subglacial environments. © 2013 The Authors. Hydrological Processes Published by John Wiley & Sons, Ltd.     

Spatial variability of soils developing on basalt flows in the Potrillo volcanic field,southern New Mexico: prelude to a chronosequence study

If the systematic spatial variability of soils in a chronosequence is identified and accounted for, the accuracy of quantitative data derived from soil chronosequence studies will be increased. A sample design using landscape positions with minimal variability could result in more accurate chronofunctions from these studies. Four basalt flows in the Potrillo volcanic field, southern New Mexico, with ages ranging between 20 ka and 260 ka (⁴⁰Ar/³⁹Ar and/or cosmogenic ³He methods) provide a sound basis for a soil chronosequence study. Basalt flow surface relief in the Potrillos reduces with time as depressions fill with basalt rubble and aeolian dust. Soil variability is primarily a function of landscape position with respect to ridges and swales in the original basalt flow topography. Soils developing over original topographic lows (swale soils) form primarily in aeolian dust, have larger amounts of total carbonate and soluble salts, and display greater variability than soils developing over original topographic highs (ridge soils). It is thus concluded that ridge soils, which have minimal variability, should be employed for a soil chronosequence study of basalt surfaces in the Potrillo volcanic field. The spatial variability of swale soils results in part from the significant hydrologic variability of low-lying landscape positions. Depth profiles of chloride concentrations suggest that hydrologic variability systematically correlates with the size and shape of depressions in which soils are forming. The infilling of depressions with aeolian material results in increasingly arid hydrologic conditions both by increasing the volume of aeolian material that is being drained and by reducing the catchment area for runoff into the depression. Depressions fill at different rates, however, depending on their size, shape and catchment area. Small, narrow depressions fill quickly, and their associated soils form under more arid conditions and have stronger development than soils in large depressions. Therefore, a number of geomorphic surfaces of varying age may develop on a single isochronous basalt flow. Each of these surfaces will have unique hydrologic characteristics and consequently different degrees of soil development. The pre-burial high water flux evident in depressions suggests that basalt flows may play an important role in aquifer recharge in this area of New Mexico. Copyright © 1999 John Wiley & Sons, Ltd.     

Influence of river channel morphology and bank characteristics on water surface boundary delineation using high‐resolution passive remote sensing and template matching

Accurate mapping of water surface boundaries in rivers is an important step for monitoring water stages, estimating discharge, flood extent, and geomorphic response to changing hydrologic conditions, and assessing riverine habitat. Nonetheless, it is a challenging task in spatially and spectrally heterogeneous river environments, commonly characterized by high spatiotemporal variations in morphology, bed material, and bank cover. In this study, we investigate the influence of channel morphology and bank characteristics on the delineation of water surface boundaries in rivers using high spatial resolution passive remote sensing and a template‐matching (object‐based) algorithm, and compare its efficacy with that of Support Vector Machine (SVM) (pixel‐based) algorithm. We perform a detailed quantitative evaluation of boundary‐delineation accuracy using spatially explicit error maps in tandem with the spatial maps of geomorphic and bank classes. Results show that template matching is more successful than SVM in delineating water surface boundaries in river sections with spatially challenging geomorphic landforms (e.g. sediment bar structures, partially submerged sediment deposits) and shallow water conditions. However, overall delineation accuracy by SVM is higher than that of template matching (without iterative hierarchical learning). Vegetation and water indices, especially when combined with texture information, improve the accuracy of template matching, for example, in river sections with overhanging trees and shadows – the two most problematic conditions in water surface boundary delineation. By identifying the influence of channel morphology and bank characteristics on water surface boundary mapping, this study helps determine river sections with higher uncertainty in delineation. In turn, the most suitable methods and data sets can be selectively utilized to improve geomorphic/hydraulic characterization. The methodology developed here can also be applied to similar studies on other geomorphic landforms including floodplains, wetlands, lakes, and coastlines. Copyright © 2014 John Wiley & Sons, Ltd.     

The effects of topographic depressions on multiscale overland flow connectivity: A high‐resolution spatiotemporal pattern analysis approach based on connectivity statistics

In watershed modelling, the traditional practice of arbitrarily filling topographic depressions in digital elevation models has raised concerns. Advanced high‐resolution remote sensing techniques, including airborne scanning laser altimetry, can identify naturally occurring depressions that impact overland flow. In this study, we used an ensemble physical and statistical modelling approach, including a 2D hydraulic model and two‐point connectivity statistics, to quantify the effects of depressions on high‐resolution overland flow patterns across spatial scales and their temporal variations in single storm events. Computations for both models were implemented using graphic processing unit‐accelerated computing. The changes in connectivity statistics for overland flow patterns between airborne scanning laser altimetry‐derived digital elevation models with (original) and without (filled) depressions were used to represent the shifts of overland flow response to depressions. The results show that depressions can either decrease or increase (to a lesser degree and shorter duration) the probability that any two points (grid locations) are hydraulically connected by overland flow pathways. We used macro‐connectivity states (Φ) as a watershed‐specific indicator to describe the spatiotemporal thresholds of connectivity variability caused by depressions. Four states of Φ are identified in a studied watershed, and each state represents different magnitudes of connectivity and connectivity changes (caused by depressions). The magnitude of connectivity variability corresponds to the states of Φ, which depend on the topological relationship between depressions, the rising/recession limb, and the total rainfall amount in a storm event. In addition, spatial distributions of connectivity variability correlate with the density of depression locations and their physical structures, which cause changes in streamflow discharge magnitude. Therefore, this study suggests that depressions are “nontrivial” in watershed modelling, and their impacts on overland flow should not be neglected. Connectivity statistics at different spatial scales and time points within a watershed provide new insights for characterizing the distributed and accumulated effects of depressions on overland flow.     

Field‐Scale Relationships Among Soil Properties and Shallow Groundwater Quality

It is important to understand the link between land surface/soil properties and shallow groundwater quality. To that end, soil properties and near‐water‐table groundwater chemistry of a shallow, unconfined aquifer were measured on a 100‐m grid on a 64‐ha irrigated field in southeastern North Dakota. Soil properties and hydrochemistry were compared via multivariate analysis that included product‐moment correlations and factor analysis/principal component analysis. Topographic low areas where the water table was in close proximity to the soil surface generally had higher apparent electrical conductivity (EC_a) and higher percent silt and clay than higher positions on the landscape. The majority of the groundwater was characterized by Ca‐ and Mg‐HCO₃ type water and was associated with topographic high areas with lower EC_a and net groundwater recharge. Small topographic depressions were areas of higher EC_a (net groundwater discharge) where salts that precipitated via evapotranspiration and evaporative discharge dissolved and leached to the groundwater during short‐term depression‐focused recharge events. At this site, groundwater quality and soil EC_a were related to surface topography. High‐resolution topography and EC_a measurements are necessary to characterize the land surface/soil properties and surficial groundwater quality at the field‐scale and to delineate areas where the shallow groundwater is most susceptible to contamination.     

Stereophotogrammetry of archive data and topographic approaches to debris‐flow torrent measurements: calculation of channel‐sediment states and a partial sediment budget for Manival torrent (Isère,France)

This paper focuses on a topographic methodology to characterize the amount of sediment stored in channels and the use of historical photographs for aerial survey by stereophotogrammetry, as part of wider research on debris‐flow magnitude prediction. The topographic methodology uses equidistant four‐point cross‐sections along the long profile of the channel. Each cross‐section is representative of a 50‐m reach of the channel. To calculate the volume of each reach, the difference is calculated between a reference level and the topographic surface. The reference level is the lowest level where the debris flow can erode, and in the current method this level is estimated from fixed points along the long profile of the channel. The accuracy of the method has been estimated by comparing results of a detailed topographic survey, with a standard deviation corresponding to about 6 per cent of the total calculated sediment volume. This topographic methodology has been used on aerial photographs by photogrammetry. This tool was applied to photographs taken on 12 past dates. The scales of the archive photographs used range from 1:3000 to 1:30 000, but results are consistent and permit us to calculate sediment states of the channel for different past dates with an uncertainty of about 6 per cent of the total volume. The application of the technique to the Manival debris‐flow torrent has permitted us to propose some partial sediment budgets and erosion‐rate estimates. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulated effect of a forest road on near‐surface hydrologic response: redux

In the work reported here the comprehensive physics‐based Integrated Hydrology Model (InHM) was employed to conduct both three‐ and two‐dimensional (3D and 2D) hydrologic‐response simulations for the small upland catchment known as C3 (located within the H. J. Andrews Experimental Forest in Oregon). Results from the 3D simulations for the steep unchannelled C3 (i) identify subsurface stormflow as the dominant hydrologic‐response mechanism and (ii) show the effect of the down‐gradient forest road on both the surface and subsurface flow systems. Comparison of the 3D results with the 2D results clearly illustrates the importance of convergent subsurface flow (e.g. greater pore‐water pressures in the hollow of the catchment for the 3D scenario). A simple infinite‐slope model, driven by subsurface pore‐water pressures generated from the 3D and 2D hydrologic‐response simulations, was employed to estimate slope stability along the long‐profile of the C3 hollow axis. As expected, the likelihood of slope failure is underestimated for the lower pore pressures from the 2D hydrologic‐response simulation compared, in a relative sense, to the higher pore pressures from the 3D hydrologic response simulation. The effort reported herein provides a firm quantitative foundation for generalizing the effects that forest roads can have on near‐surface hydrologic response and slope stability at the catchment scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Two‐dimensional continuous simulation of spatiotemporally varied hydrological processes using the time‐area method

Distributed, continuous hydrologic models promote better understanding of hydrology and enable integrated hydrologic analyses by providing a more detailed picture of water transport processes across the varying landscape. However, such models are not widely used in routine modelling practices, due in part to the extensive data input requirements, computational demands, and complexity of routing algorithms. We developed a two‐dimensional continuous hydrologic model, HYSTAR, using a time‐area method within a grid‐based spatial data model with the goal of providing an alternative way to simulate spatiotemporally varied watershed‐scale hydrologic processes. The model calculates the direct runoff hydrograph by coupling a time‐area routing scheme with a dynamic rainfall excess sub‐model implemented here using a modified curve number method with an hourly time step, explicitly considering downstream ‘reinfiltration’ of routed surface runoff. Soil moisture content is determined at each time interval based on a water balance equation, and overland and channel runoff is routed on time‐area maps, representing spatial variation in hydraulic characteristics for each time interval in a storm event. Simulating runoff hydrographs does not depend on unit hydrograph theory or on solution of the Saint Venant equation, yet retains the simplicity of a unit hydrograph approach and the capability of explicitly simulating two‐dimensional flow routing. The model provided acceptable performance in predicting daily and monthly runoff for a 6‐year period for a watershed in Virginia (USA) using readily available geographic information about the watershed landscape. Spatial and temporal variability in simulated effective runoff depth and time area maps dynamically show the areas of the watershed contributing to the direct runoff hydrograph at the outlet over time, consistent with the variable source area overland flow generation mechanism. The model offers a way to simulate watershed processes and runoff hydrographs using the time‐area method, providing a simple, efficient, and sound framework that explicitly represents mechanisms of spatially and temporally varied hydrologic processes. Copyright © 2015 John Wiley & Sons, Ltd.     

Efficient hybrid breaching‐filling sink removal methods for flow path enforcement in digital elevation models

Digital elevation models (DEMs) that are used in hydrological applications must be processed to remove sinks, mainly topographic depressions. Flow enforcement techniques include filling methods, which raise elevations within depressions, breaching, which carves channels through blockages, and hybrid methods. Despite previous research demonstrating the large impact to DEMs and subsequent analyses of depression filling, it is common practice apply this technique to flow enforcement. This is partly because of the greater efficiency of depression filling tools compared to breaching counterparts, which often limits breaching to applications of small‐ to moderate‐sized DEMs. A new hybrid flow enforcement algorithm is presented in this study. The method can be run in complete breaching, selective breaching (either breached or filled), or constrained breaching (partial breaching) modes, allowing for greater flexibility in how practitioners enforce continuous flow paths. Algorithm performance was tested with DEMs of varying topography, spatial extents, and resolution. The sites included three moderate sized DEMs (52 000 000 to 190 000 000 cells) and three massive DEMs of the Iberian Peninsula, and the Amazon and Nile River basins, the largest containing nearly one billion cells. In complete breaching mode, the new algorithm required 87% of the time needed by a filling method to process the test DEMs, while the selective breaching and constrained breaching modes, operating with maximum breach depth constraints, increased run times by 8% and 27% respectively. Therefore, the new algorithm offers comparable performance to filling and the ability to process massive topographic data sets, while giving practitioners greater flexibility and lowering DEM impact. Copyright © 2015 John Wiley & Sons, Ltd.     

Incorporating landscape depression heterogeneity into the Soil and Water Assessment Tool (SWAT) using a probability distribution

Modelling the hydrology of North American Prairie watersheds is complicated because of the existence of numerous landscape depressions that vary in storage capacity. The Soil and Water Assessment Tool (SWAT) is a widely applied model for long‐term hydrological simulations in watersheds dominated by agricultural land uses. However, several studies show that the SWAT model has had limited success in handling prairie watersheds. In past works using SWAT, landscape depression storage heterogeneity has largely been neglected or lumped. In this study, a probability distributed model of depression storage is introduced into the SWAT model to better handle landscape storage heterogeneity. The work utilizes a probability density function to describe the spatial heterogeneity of the landscape depression storages that was developed from topographic characteristics. The integrated SWAT–PDLD model is tested using datasets for two prairie depression dominated watersheds in Canada: the Moose Jaw River watershed, Saskatchewan; and the Assiniboine River watershed, Saskatchewan. Simulation results were compared to observed streamflow using graphical and multiple statistical criterions. Representation of landscape depressions within SWAT using a probability distribution (SWAT–PDLD) provides improved estimations of streamflow for large prairie watersheds in comparison to results using a lumped, single storage approach. Copyright © 2016 John Wiley & Sons, Ltd.