NETWORK AND SUBWATERSHED PARAMETERS EXTRACTED FROM DIGITAL ELEVATION MODELS: THE BILLS CREEK EXPERIENCE1

ABSTRACT: An automated extraction of channel network and sub-watershed characteristics from digital elevation models (DEM) is performed by model DEDNM. This model can process DEM data of limited vertical resolution representing low relief terrain. Such representations often include ill-defined drainage boundaries and indeterminate flow paths. The application watershed is an 84 km² low relief watershed in southwestern Oklahoma. The standard for validation is the network and subwatershed parameters defined by the blue line method on USGS 7.5–minute topographic maps. Evaluation of the generated and validation networks by visual comparisons shows a high degree of correlation. Comparison of selected network parameters (channel length, slope, drainage density, etc.) and of drainage network composition (bifurcation, length, slope, and area ratios) shows that, on the average, the generated parameters are within 5 percent of those derived from the validation network. The largest discrepancies were found for the channel slope values. The results of this application demonstrate that DEDNM effectively addresses network definition problems often encountered in low relief terrain and that it can generate accurate network and subwatershed parameters under those conditions.     

A GEOGRAPHIC INFORMATION SYSTEM PROCEDURE TO QUANTIFY DRAINAGE-BASIN CHARACTERISTICS1

ABSTRACT: The Basin Characteristics System (BCS) has been developed to quantify characteristics of a drainage basin. The first of four main BCS processing steps creates four geographic information system (GIS) digital maps representing the drainage divide, the drainage network, elevation contours, and the basin length. The drainage divide and basin length are manually digitized from 1:250,000-scale topographic maps. The drainage network is extracted using GIS software from 1:100,000-scale digital line graph data. The elevation contours are generated using GIS software from 1:250,000-scale digital elevation model data. The second and third steps use software developed to assign attributes to specific features in three of the four digital maps and analyze the four maps to quantify 24 morphometric basin characteristics. The fourth step quantifies two climatic characteristics from digitized State maps of precipitation data. Compared to manual methods of measurement, the BCS provides a reduction in the time required to quantify the 26 basin characteristics. Comparison tests indicate the BCS measurements are not significantly different from manual topographic-map measurements for 11 of 12 primary drainage-basin characteristics. Tests indicate the BCS significantly underestimates basin slope. Comparison-measurement differences for basin slope, main channel slope, and basin relief appear to be due to limitations in the digital elevation model data.     

CHANNEL NETWORK EXTENSION BY LOGGING ROADS IN TWO BASINS,WESTERN CASCADES,OREGON1

ABSTRACT: Based on field surveys and analysis of road networks using a geographic information system (GIS), we assessed the hydrologic integration of an extensive logging-road network with the stream network in two adjacent 62 and 119 km² basins in the western Cascades of Oregon. Detailed surveys of road drainage for 20 percent of the 350 km road network revealed two hydrologic flow paths that link roads to stream channels: roadside ditches draining to streams (35 percent of the 436 culverts examined), and roadside ditches draining to culverts with gullies incised below their outlets (23 percent of culverts). Gully incision is significantly more likely below culverts on steep (< 40 percent) slopes with longer than average contributing ditch length. Fifty-seven percent of the surveyed road length is connected to the stream network by these surface flowpaths, increasing drainage density by 21 to 50 percent, depending on which road segments are assumed to be connected to streams. We propose a conceptual model to describe the hydrologic function of roads based on two effects: (1) a volumetric effect, increasing the volume of water available for quickflow and (2) a timing effect, altering flow-routing efficiency through extensions to the drainage network. This study examines the second of these two effects. Future work must quantify discharge along road segments connected to the stream network in order to more fully explain road impacts on basin hydrology.     

DEM AGGREGATION FOR WATERSHED MODELING1

ABSTRACT: The widely available USGS 7.5‐minute Digital Elevation Model (DEM) has a cell size of approximately 30 m × 30 m. This high resolution topographic information is impractical for many applications of distributed hydrologic and water quality models. In this study, cells were aggregated into coarse‐resolution areal units, termed grids, and a method to approximate flow direction for coarse‐resolution grids from 30 m DEM cells was developed. The method considers the flow path defined from the fine‐resolution DEM in determining a grid's flow direction and makes flow directions for grids closely follow the flow pattern suggested by the DEM. The aggregation method was applied to a DEM of Goodwater Creek, a nearly flat watershed that is located in central Missouri. The drainage networks derived for different levels of cell aggregations showed that grid aggregates of the Goodwater Creek watershed provided an adequate representation of the landscape topography.     

SENSITIVITY CONSIDERATIONS WHEN MODELING HYDROLOGIC PROCESSES WITH DIGITAL ELEVATION MODEL1

ABSTRACT: The purpose of this paper is to investigate the sensitivity of a hydrologic models to the type of DEM used. This was done while modeling basin water quality with 1:24,000 and 1:250,000 U.S. Geological Survey DEMs as input to model hydro‐logical processes. The manner in which the model results were sensitive to the choice of raster cell size (scale) is investigated in this study. The Broadhead watershed, located in New Jersey, USA, was chosen as a study area. Curve numbers were estimated by a trial and error to match simulated and observed total discharge. Monthly runoff for the watershed was used in the calibration process. Higher runoff volumes were simulated by the model when the 1:24,000 DEM were used as input data, probably due to the finer resolution which simulated increased average slope and hence higher estimated runoff from the watershed. As the simulated slope of the watershed is flatten with the 1:250,000 DEM, the response of stream flow was delayed and simulated less runoff volume.     

Evaluating the Slope‐Area Method to Accurately Identify Stream Channel Heads in Three Physiographic Regions

Estimation of stream channel heads is an important task since ephemeral channels play a significant role in the transport of sediment and materials to perennial streams. The slope‐area method utilizes digital elevation model (DEM) and related information to develop slope‐area threshold relationships used to estimate the position of channel heads in the watershed. A total of 162 stream channel heads were mapped across the three physiographic regions of Alabama, including the Southwestern Appalachians (51), Piedmont/Ridge and Valley (61), and Coastal Plains (51). Using Geographic Information System and DEM, the local slope and drainage area for each mapped channel head was calculated and region‐specific models were developed and evaluated. Results demonstrated the local slope and drainage area had an inverse and strong correlation in the Piedmont/Ridge and Valley region (r² = 0.71) and the Southwestern Appalachian region (r² = 0.61). Among three physiographic regions, the weakest correlation was observed in the Coastal Plain region (r² = 0.45). By comparing the locations of modeled channel heads to those located in the field, calculated reliability and sensitivity indices indicated model accuracy and reliance were weak to moderate. However, the slope‐area method helped define the upstream boundaries of a more detailed channel network than that derived from the 1:24,000‐scale National Hydrography Dataset, which is commonly used for planning and regulatory purposes.     

River pollution from non-point sources: a new simplified method of assessment

Assessment of the pollution of water bodies from non-point sources is a complex data- and time-consuming task. The potential non-point pollution index (PNPI), is a new tool designed to assess the global pressure exerted on rivers and other surface water bodies by different land uses. The main feature of PNPI is the wide availability of its input data. Very detailed input maps, often lacking over many areas, are not needed for PNPI calculation. As a consequence of the input data used, the modelling of physical reality and of processes is heavily simplified. The authors counterbalanced such a simplification using an 'expert system' approach. The system bypasses the accurate representation of the physical reality to assess globally the pollution potential of different land uses according to the judgement of scientists. The scientific community proposes many models for depicting the dynamics of pollutants coming from diffuse sources. Most of them can be grouped into two broad categories: statistical models and physically based models. PNPI belongs to neither of the above-mentioned groups. PNPI is a GIS-based, watershed-scale tool designed to inform decision makers and public opinion about the potential environmental impacts of different land management scenarios. PNPI applies the multicriteria technique to pollutant dynamics and water quality. The pressure exerted on water bodies by diffuse pollution coming from land units is expressed as a function of three indicators: land use, run-off and distance from the river network. They are calculated from land use data, geological maps and a digital elevation model (DEM). The weights given to different land uses and to the three indicators were set according to experts' evaluations and allow calculation of the value of the PNPI for each node of a grid representing the watershed; the higher the PNPI of the cell, the greater the potential impact on the river network. The output of the calculation is presented in the form of maps that highlight areas that are more likely to produce pollution. Last, possibilities, strategies and results of the validation of the PNPI are described. In the authors' view, the explicit link between land use and potential pollution on which PNPI is based, together with its high communication potential, make it particularly interesting for a participatory and integrated approach to land management and environmental protection.     

WATERSHED CONFIGURATION AND GEOGRAPHIC INFORMATION SYSTEM PARAMETERIZATION FOR SPUR MODEL HYDROLOGIC SIMULATIONS1

ABSTRACT: A grid cell geographic information system (GIS) is used to parameterize SPUR, a quasi-physically based surface runoff model in which a watershed is configured as a set of stream segments and contributing areas. GIS analysis techniques produce various watershed configurations by progressive simplification of a stream network delineated from digital elevation models (DEM). We used three watershed configurations: ≥ 2nd, ≥ 4th, and ≥ 13th Shreve order networks, where the watershed contains 28, 15, and 1 channel segments with 66, 37, and 3 contributing areas, respectively. Watershed configuration controls simulated daily and monthly sums of runoff volumes. For the climatic and topographic setting in southeastern Arizona the ≥ 4th order configuration of the stream network and contributing areas produces results that are typically as good as the ≥ 2nd order network. However both are consistently better than the ≥ 13th order configuration. Due to the degree of parameterization in SPUR, model simulations cannot be significantly improved by increasing watershed configuration beyond the ≥ 4th order network. However, a range of Soil Conservation Service curve numbers derived from rainfall/runoff data can affect model simulations. Higher curve numbers yield better results for the ≥ 2nd order network while lower curve numbers yield better results for the ≥ 4th order network.     

CHANGES IN HYDROLOGIC RESPONSE DUE TO STREAM NETWORK EXTENSION VIA LAND DRAINAGE ACTIVITIES1

ABSTRACT: The objective of this work is to determine the effects of extension of a stream network through land drainage activities during the late 1800s on the hydrologic response of a watershed. The Mackinaw River Basin in Central Illinois was chosen as the focus and the pre‐land and post‐land drainage activity hydrologic responses were obtained through convolution of the hill slope and channel responses and compared. The hill slope response was computed using the kinematic wave model and the channel response was determined using the geomorphologic instantaneous unit hydrograph method. Our hypothesis was that the hydrologic response of the basin would exhibit the characteristic effects of settlement (i.e., increases in peak discharges and decreases in times to peak). This, indeed, is what occurred; however, the increase in peak discharges diminishes as scale increases, leaving only the decrease in times to peak. At larger scales, the dispersive effects of the longer hill slope lengths in the pre‐settlement scenario seem to balance the depressive effects of the longer path lengths in the post‐settlement scenario, thus the pre‐settlement and post‐settlement peak discharges are approximately equivalent. At small scales, the dispersion caused by the hill slope is larger in the pre‐settlement case; thus, the post‐settlement peak discharges are greater than the pre‐settlement.     

Development and Comparison of Multiple Regression Models to Predict Bankfull Channel Dimensions for Use in Hydrologic Models

Channel dimensions are important input variables for many hydrologic models. As measurements of channel geometry are not available in most watersheds, they are often predicted using bankfull hydraulic geometry relationships. This study aims at improving existing equations that relate bankfull width, depth, and cross‐sectional area to drainage area (DA) without limiting their use to well‐gauged watersheds. We included seven additional variables in the equations that can be derived from data that are generally required by hydrologic models anyway and conducted several multiple regression analyses to identify the ideal combination of additional variables for nationwide and regional models for each Physiographic Division of the United States (U.S.). Results indicate that including the additional variables in the regression equations generally improves predictions considerably. The selection of relevant variables varies by Physiographic Division, but average annual precipitation (PCP) and temperature (TMP) were generally found to improve the models the most. Therefore, we recommend using regression equations with three independent variables (DA, PCP, and TMP) to predict bankfull channel dimensions for hydrologic models. Furthermore, we recommend using the regional equations for watersheds within regions from which data were used for model development, whereas in all other parts of the U.S. and the rest of the world, the nationwide equations should be given preference.     

INTEGRATING LANDSCAPE ASSESSMENT AND HYDROLOGIC MODELING FOR LAND COVER CHANGE ANALYSIS1

ABSTRACT: Significant land cover changes have occurred in the watersheds that contribute runoff to the upper San Pedro River in Sonora, Mexico, and southeast Arizona. These changes, observed using a series of remotely sensed images taken in the 1970s, 1980s, and 1990s, have been implicated in the alteration of the basin hydrologic response. The Cannonsville subwatershed, located in the Catskill/Delaware watershed complex that delivers water to New York City, provides a contrast in land cover change. In this region, the Cannonsville watershed condition has improved over a comparable time period. A landscape assessment tool using a geographic information system (GIS) has been developed that automates the parameterization of the Soil and Water Assessment Tool (SWAT) and KINEmatic Runoff and EROSion (KINEROS) hydrologic models. The Automated Geospatial Watershed Assessment (AGWA) tool was used to prepare parameter input files for the Upper San Pedro Basin, a subwatershed within the San Pedro undergoing significant changes, and the Cannonsville watershed using historical land cover data. Runoff and sediment yield were simulated using these models. In the Cannonsville watershed, land cover change had a beneficial impact on modeled watershed response due to the transition from agriculture to forest land cover. Simulation results for the San Pedro indicate that increasing urban and agricultural areas and the simultaneous invasion of woody plants and decline of grasslands resulted in increased annual and event runoff volumes, flashier flood response, and decreased water quality due to sediment loading. These results demonstrate the usefulness of integrating remote sensing and distributed hydrologic models through the use of GIS for assessing watershed condition and the relative impacts of land cover transitions on hydrologic response.     

COMPARISON OF MODELS FOR ESTIMATING FLOOD QUANTILES IN NEW HAMPSHIRE AND VERMONT1

We develop and compare three regression models for estimating flood quantiles at ungaged stream reaches in New Hampshire and Vermont. These models emerge from systematic analysis and validation of relations between flood magnitude and six candidate predictors reflecting basin size, topography, and climate and channel size at 36 gaging stations with record lengths exceeding 20 years. Of the candidate predictors, bank full width is most highly correlated with flood magnitude and the best prediction equation is based on width. Thus channel geometry is closely related to the current hydrologic regime in spite of geologically recent glaciation and apparently non-alluvial bank materials. We also develop models that use information obtainable from maps or GIS. The best of these uses drainage area and drainage-basin elevation as predictors, but it is substantially less precise than the width-based relation. A third relation using only drainage area as a predictor is even less precise but may be useful for some purposes. No other single predictors or combinations yielded useful predictions, although some had been included in previously-established models for the region. Model comparison included examination of residuals generated by regression using one-at-a-time suppression of data points and comparison with precision obtainable with gaging records of varying lengths.     

INTEGRATION OF GEOGRAPHIC INFORMATION SYSTEMS AND A COMPUTER MODEL TO EVALUATE IMPACTS OF AGRICULTURAL RUNOFF ON WATER QUALITY1

ABSTRACT: This study integrates an Agricultural Non-Point Source Pollution Model (AGNPS), the Geographic Resource Analysis Support System (GRASS) (U.S. Army Corps of Engineers, 1987), and GRASS WATERWORKS (a hydrologic modeling tool box being developed at the Michigan State University Center for Remote Sensing) to evaluate the impact of agricultural runoff on water quality in the Cass River, a subwatershed of Saginaw Bay. AGNPS is used to estimate the amounts, origin, and distribution of sediment, nitrogen (N), and phosphorus (P) in the watershed. GRASS and GRASS WATERWORKS are used to generate parameters needed for AGNPS from digital maps, which include soil association, land use, watershed boundaries, water features, and digital elevation. Outputs of the model include spatially distributed estimates of volume and peak runoff, overland and channel erosion, sediment yields, and concentrations of nitrogen and phosphorus. Management scenarios are explored in the AGNPS model to minimize sedimentation and nutrient loading. Scenarios evaluated include variations in crop cover, tillage methods, and other agricultural management practices. In addition, areas vulnerable to erosion are identified for best management practices.     

DEVELOPMENT OF LARGE SCALE GRIDDED RIVER NETWORKS FROM VECTOR STREAM DATA1

ABSTRACT: The Network Tracing Method (NTM) has been developed to determine gridded coarse river networks for modeling large hydrologic systems. For a coarse resolution grid, the NTM determines the downstream cell of each cell and the distance along the actual meandering flow paths between them. Unlike previously developed methods, the NTM uses fine resolution vector river networks as the source of information of the flow patterns rather than digital elevation models. The main advantage of using vector river networks as input is that they capture the hydrologic terrain features better than topographic data do, particularly in areas of low topographic relief. The NTM was applied to South America with a grid resolution of 1 degree by 1 degree and to the globe with a resolution of 2.815 degrees by 2.8125 degrees. Overall, the method captured the flow patterns well. Generated digital river networks and drainage divides showed minor disagreement with those obtained from existing maps, and most of them were consistent with the resolution of the coarse river network. The majority of estimated basin areas were also close to documented values. River lengths calculated with the NTM, however, were consistently underpredicted.     

EXTRACTION AND UTILIZATION OF SPACE ACQUIRED PHYSIOGRAPHIC DATA FOR WATER RESOURCES DEVELOPMENT1

ABSTRACT: ERTS-1 satellite imagery has been evaluated as a means of providing useful watershed physiography information. From these data physiographic parameters such as drainage basin area and shape, drainage density, stream length and sinuosity, and the percentage of a watershed occupied by major land use types were obtained in three study areas. The study areas were: (1) Southwestern Wisconsin; (2) Eastern Colorado; and (3) portions of the Middle Atlantic States Using ERTS-1 imagery at 1:250,000 and 1:100,000 scales it was found that drainage basin area and shape and stream sinuosity were comparable (within 10%) in all study areas to physiographic measurements derived from conventional topographic maps at the same scales Land use information can be usefully extracted for watersheds as small as 30 mi²(78 km²) in area. Improved drainage network and density information is obtained from ERTS-1 imagery in dissected areas such as Southwestern Wisconsin, but in heavily vegetated areas (Middle Atlantic States) or areas with little physical relief (Eastern Colorado) low order streams are difficult to detect and the derived drainage densities are significantly smaller than those obtained from standard maps. It is concluded that ERTS-1 imagery can be employed to advantage in mean annual runoff prediction techniques and in providing or maintaining land use information used in the calibration and operation of watershed models.     

Mapping Outlets of Iowa Flood Center and National Water Center River Networks for Hydrologic Model Comparison

River networks based on Digital Elevation Model (DEM) data differ depending on the DEM resolution, accuracy, and algorithms used for network extraction. As spatial scale increases, the differences diminish. This study explores methods that identify the scale where networks obtained by different methods agree within some margin of error. The problem is relevant for comparing hydrologic models built around the two networks. An example is the need to compare streamflow prediction from the Hillslope Link Model (HLM) operated by the Iowa Flood Center (IFC) and the National Water Model (NWM) operated by the National Water Center of the National Oceanic and Atmospheric Administration. The HLM uses landscape decomposition into hillslopes and channel links while the NWM uses the NHDPlus dataset as its basic spatial support. While the HLM resolves the scale of the NHDPlus, the outlets of the latter do not necessarily correspond to the nodes of the HLM model. The authors evaluated two methods to map the outlets of NHDPlus to outlets on the IFC network. The methods compare the upstream areas of the channels and their spatial location. Both methods displayed similar performance and identified matches for about 80% of the outlets with a tolerance of 10% in errors in the upstream area. As the aggregation scale increases, the number of matches also increases. At the scale of 100 km², 90% of the outlets have matches with tolerance of 5%. The authors recommend this scale for comparing the HLM and NWM streamflow predictions.     

APPLICATIONS OF A GIS FOR MODELING THE SENSITIVITY OF WATER RESOURCES TO ALTERATIONS IN CLIMATE IN THE GUNNISON RIVER BASIN,COLORADO1

ABSTRACT: The Gunnison River drains a mountainous basin in western Colorado, and is a large contributor of water to the Colorado River. As part of a study to assess water resource sensitivity to alterations in climate in the Gunnison River basin, climatic and hydrologic processes are being modeled. A geographic information system (GIS) is being used in this study as a link between data and modelers - serving as a common data base for project personnel with differing specialties, providing a means to investigate the effects of scale on model results, and providing a framework for the transfer of parameter values among models. Specific applications presented include: (1) developing elevation grids for a precipitation model from digital elevation model (DEM) point-elevation values, and visualizing the effects of grid resolution on model results; (2) using a GIS to facilitate the definition and parameterization of a distributed-parameters, watershed model in multiple basins; and (3) nesting atmospheric and hydrologic models to produce possible scenarios of climate change.     

Modeling Stream Channel Characteristics From Drainage‐Enforced DEMs in Puget Sound,Washington, USA1

Abstract: Mapping stream channels and their geomorphic attributes is an important step in many watershed research and management projects. Often insufficient field data exist to map hydromorphologic attributes across entire drainage basins, necessitating the application of hydrologic modeling tools to digital elevation models (DEMs) via a geographic information system (GIS). In this article, we demonstrate methods for deriving synthetic stream networks via GIS across large and diverse basins using drainage‐enforced DEMs, along with techniques for estimating channel widths and gradient on the reach scale. The two‐step drainage enforcement method we used produced synthetic stream networks that displayed a high degree of positional accuracy relative to the input streams. The accuracies of our estimated channel parameters were assessed with field data, and predictions of bankfull width, wetted width and gradient were strongly correlated with measured values (r² = 0.92, r² = 0.95, r² = 0.88, respectively). Classification accuracies of binned channel attributes were also high. Our methodology allows for the relatively rapid mapping of stream channels and associated morphological attributes across large geographic areas. Although initially developed to provide salmon recovery planners with important salmon habitat information, we suggest these methodologies are relevant to a variety of research and management questions.     

EFFECT OF WATERSHED SUBDIVISION ON SWAT FLOW,SEDIMENT, AND NUTRIENT PREDICTIONS1

ABSTRACT: The size, scale, and number of subwatersheds can affect a watershed modeling process and subsequent results. The objective of this study was to determine the appropriate level of subwatershed division for simulating flow, sediment, and nutrients over 30 years for four Iowa watersheds ranging in size from 2,000 to 18,000 km² with the Soil and Water Assessment Tool (SWAT) model. The results of the analysis indicated that variation in the total number of subwatersheds had very little effect on streamflow. However, the opposite result was found for sediment, nitrate, and inorganic P; the optimal threshold subwatershed sizes, relative to the total drainage area for each watershed, required to adequately predict these three indicators were found to be around 3, 2, and 5 percent, respectively. Decreasing the size of the subwatersheds below these threshold levels does not significantly affect the predicted levels of these environmental indicators. These threshold subwatershed sizes can be used to optimize input data preparation requirements for SWAT analyses of other watersheds, especially those within a similar size range. The fact that different thresholds emerged for the different indicators also indicates the need for SWAT users to assess which indicators should have the highest priority in their analyses.     

FACTORS AFFECTING FLOODS FROM WATERSHEDS IN HUMID REGIONS OF NORTHEASTERN OHIO1

This paper describes the results of a study of hydrologic factors affecting floods from humid region in northeastern Ohio. Statistical multiple correlation analysis was used to relate floods to hydrologic and basin characteristics. Results of the study emphasize that the characteristics of floods from small and large watersheds are so significantly different that the two problems cannot be combined into one solution. The studies show that the most important hydrologic characteristics in large watersheds were: drainage area size and main channel slope. For small watersheds the most important hydrologic characteristics were: drainage area size, rainfall intensity and soil index. For watershed effect by reservoir storage it was found that: (1) small drainage areas are relatively more affected by storage than large drainage areas; (2) storage of less than 25 acre feet per square mile will not have significant effect on the mean annual flood (for drainage area above 70 square miles).