Characterizing the Extent of Spatially Integrated Floodplain and Wetland Systems in the White River,Indiana, USA

Floodplain delineation may inform protection of wetland systems under local, state, or federal laws. Nationally available Federal Emergency Management Agency Flood Insurance Rate Maps (FIRMs, “100‐year floodplain” maps) focus on urban areas and higher‐order river systems, limiting utility at large scales. Few other national‐scale floodplain data are available. We acquired FIRMs for a large watershed and compared FIRMs to floodplain and integrated wetland area mapping methods based on (1) geospatial distance, (2) geomorphic setting, and (3) soil characteristics. We used observed flooding events (OFEs) with recurrence intervals of 25‐50 to >100 years to assess floodplain estimate accuracy. FIRMs accurately reflected floodplain areas based on OFEs and covered 32% of river length, whereas soil‐based mapping was not as accurate as FIRMs but characterized floodplain areas over approximately 65% of stream length. Geomorphic approaches included more areas than indicated by OFE, whereas geospatial approaches tended to cover less area. Overall, soil‐based methods have the highest utility in determining floodplains and their integrated wetland areas at large scales due to the use of nationally available data and flexibility for regional application. These findings will improve floodplain and integrated wetland system extent assessment for better management at local, state, and national scales.     

Demonstrating Floodplain Uncertainty Using Flood Probability Maps1

Abstract: The U.S. Federal Emergency Management Agency (FEMA) flood maps depict the 100‐year recurrence interval floodplain boundary as a single line. However, because of natural variability and model uncertainty, no floodplain extents can be accurately defined by a single line. This article presents a new approach to floodplain mapping that takes advantage of accepted methodologies in hydrologic and hydraulic analysis while including the effects of uncertainty. In this approach, the extents of computed floodplain boundaries are defined as a continuous map of flood probabilities, rather than as a single line. Engineers and planners can use these flood probability maps for viewing the uncertainty of a floodplain boundary at any recurrence interval. Such a flood probability map is a useful tool for visualizing the uncertainty of a floodplain boundary and represents greater honesty in engineering technologies that are used for flood mapping. While institutional barriers may prevent adoption of such definitions for use in graduated flood insurance rates (as most other insurance industries use to account for relative risks), the methods open the door technically to such a reality.     

Comparative Analysis of Inundation Mapping Approaches for the 2016 Flood in the Brazos River,Texas

Accurate and timely flood inundation maps serve as crucial information for hydrologists, first‐responders, and decision makers of natural disaster management agencies. In this study, two modeling approaches are applied to estimate the inundation area for a large flooding event that occurred in May 2016 in the Brazos River: (1) Height Above the Nearest Drainage combined with National Hydrograph Dataset Plus (NHDPlus‐HAND) and (2) International River Interface Cooperative — Flow and Sediment Transport with Morphological Evolution of Channels (iRIC‐FaSTMECH). The inundation extents simulated from these two modeling approaches are then compared against the observed inundation extents derived from a Landsat 8 satellite image. The simulated results from NHDPlus‐HAND and iRIC‐FaSTMECH show 56% and 70% of overlaps with the observed flood extents, respectively. A modified version of the NHDPlus‐HAND model, considering networked catchment behaviors, is also tested with an improved fitness of 67%. This study suggests that NHDPlus‐HAND has the potential for real‐time continental inundation forecast due to its low computational cost and ease to couple with the National Water Model. Better performance of NHDPlus‐HAND can be achieved by considering the inter‐catchment flows during extreme riverine flood events. Overall, this study presents a comprehensive examination made of remote sensing compared with HAND‐based inundation mapping in a region of complex topography.     

AutoRAPID: A Model for Prompt Streamflow Estimation and Flood Inundation Mapping over Regional to Continental Extents

This article couples two existing models to quickly generate flow and flood‐inundation estimates at high resolutions over large spatial extents for use in emergency response situations. Input data are gridded runoff values from a climate model, which are used by the Routing Application for Parallel computatIon of Discharge (RAPID) model to simulate flow rates within a vector river network. Peak flows in each river reach are then supplied to the AutoRoute model, which produces raster flood inundation maps. The coupled tool (AutoRAPID) is tested for the June 2008 floods in the Midwest and the April‐June 2011 floods in the Mississippi Delta. RAPID was implemented from 2005 to 2014 for the entire Mississippi River Basin (1.2 million river reaches) in approximately 45 min. Discretizing a 230,000‐km² area in the Midwest and a 109,500‐km² area in the Mississippi Delta into thirty‐nine 1° by 1° tiles, AutoRoute simulated a high‐resolution (~10 m) flood inundation map in 20 min for each tile. The hydrographs simulated by RAPID are found to perform better in reaches without influences from unrepresented dams and without backwater effects. Flood inundation maps using the RAPID peak flows vary in accuracy with F‐statistic values between 38.1 and 90.9%. Better performance is observed in regions with more accurate peak flows from RAPID and moderate to high topographic relief.     

Error Assessment for Height Above the Nearest Drainage Inundation Mapping

Real‐time flood inundation mapping is vital for emergency response to help protect life and property. Inundation mapping transforms rainfall forecasts into meaningful spatial information that can be utilized before, during, and after disasters. While inundation mapping has traditionally been conducted on a local scale, automated algorithms using topography data can be utilized to efficiently produce flood maps across the continental scale. The Height Above the Nearest Drainage method can be used in conjunction with synthetic rating curves (SRCs) to produce inundation maps, but the performance of these inundation maps needs to be assessed. Here we assess the accuracy of the SRCs and calculate statistics for comparing the SRCs to rating curves obtained from hydrodynamic models calibrated against observed stage heights. We find SRCs are accurate enough for large‐scale approximate inundation mapping while not as accurate when assessing individual reaches or cross sections. We investigate the effect of terrain and channel characteristics and observe reach length and slope predict divergence between the two types of rating curves, and SRCs perform poorly for short reaches with extreme slope values. We propose an approach to recalculate the slope in Manning’s equation as the weighted average over a minimum distance and assess accuracy for a range of moving window lengths.     

USING GEOLOGY TO IMPROVE FLOOD HAZARD MANAGEMENT ON ALLUVIAL FANS ‐ AN EXAMPLE FROM LAUGHLIN,NEVADA1

A study of the piedmont of the Newberry Mountains near Laughlin, Nevada, demonstrates that geologic information can improve the scientific basis of flood‐hazard management on alluvial fans in desert areas. Comparison of geologic information against flood insurance rate maps (FIRMs) reveals flaws in conventional methods for flood hazard delineation in this setting. Geologic evidence indicates that large parts of the Newberry piedmont have been isolated from significant flooding for at least the past 10,000 years. This contrasts with existing FIRMs that include large tracts of nonflood prone land in the 100‐year and 500‐year flood hazard zones and exclude areas of indisputably flood prone land from the regulatory flood plain. From the basis of the geology, flood hazards on at least one‐third of the piedmont are mischaracterized on the regulatory maps. The formal incorporation of geologic data into flood hazard studies on desert piedmonts could significantly reduce this type of discrepancy and substantially reduce the scope, hence cost, of more elaborate engineering studies and hazard mitigation strategies. The results of this study affirm the value of new Federal Emergency Management Agency (FEMA) recommendations for characterizing alluvial fan flood hazards and support an argument for mandating geological studies in the regulatory process.     

Probabilistic Flood Inundation Forecasting Using Rating Curve Libraries

One approach for performing uncertainty assessment in flood inundation modeling is to use an ensemble of models with different conceptualizations, parameters, and initial and boundary conditions that capture the factors contributing to uncertainty. However, the high computational expense of many hydraulic models renders their use impractical for ensemble forecasting. To address this challenge, we developed a rating curve library method for flood inundation forecasting. This method involves pre‐running a hydraulic model using multiple inflows and extracting rating curves, which prescribe a relation between streamflow and stage at various cross sections along a river reach. For a given streamflow, flood stage at each cross section is interpolated from the pre‐computed rating curve library to delineate flood inundation depths and extents at a lower computational cost. In this article, we describe the workflow for our rating curve library method and the Rating Curve based Automatic Flood Forecasting (RCAFF) software that automates this workflow. We also investigate the feasibility of using this method to transform ensemble streamflow forecasts into local, probabilistic flood inundation delineations for the Onion and Shoal Creeks in Austin, Texas. While our results show water surface elevations from RCAFF are comparable to those from the hydraulic models, the ensemble streamflow forecasts used as inputs to RCAFF are the largest source of uncertainty in predicting observed floods.     

A Hydraulic MultiModel Ensemble Framework for Visualizing Flood Inundation Uncertainty

While deterministic forecasts provide a single realization of potential inundation, the inherent uncertainty associated with forecasts also needs to be conveyed for improved decision support. The objective of this study was to develop an ensemble framework for the quantification and visualization of uncertainty associated with flood inundation forecast maps. An 11‐member ensemble streamflow forecast at lead times from 0 to 48 hr was used to force two hydraulic models to produce a multimodel ensemble. The hydraulic models used are (1) the International River Interface Cooperative along with Flow and Sediment Transport with Morphological Evolution of Channels solver and (2) the two‐dimensional Hydrologic Engineering Center‐River Analysis System. Uncertainty was quantified and augmented onto flood inundation maps by calculating statistical spread among the ensemble members. For visualization, a series of probability flood maps conveying the uncertainty in forecasted water extent, water depth, and flow velocity was disseminated through a web‐based decision support tool. The results from this study offer a framework for quantifying and visualizing model uncertainty in forecasted flood inundation maps.     

Coastal Flood Inundation Monitoring with Satellite C‐band and L‐band Synthetic Aperture Radar Data

Satellite Synthetic Aperture Radar (SAR) was evaluated as a method to operationally monitor the occurrence and distribution of storm‐ and tidal‐related flooding of spatially extensive coastal marshes within the north‐central Gulf of Mexico. Maps representing the occurrence of marsh surface inundation were created from available Advanced Land Observation Satellite (ALOS) Phased Array type L‐Band SAR (PALSAR) (L‐band) (21 scenes with HH polarizations in Wide Beam [100 m]) data and Environmental Satellite (ENVISAT) Advanced SAR (ASAR) (C‐band) data (24 scenes with VV and HH polarizations in Wide Swath [150 m]) during 2006‐2009 covering 500 km of the Louisiana coastal zone. Mapping was primarily based on a decrease in backscatter between reference and target scenes, and as an extension of previous studies, the flood inundation mapping performance was assessed by the degree of correspondence between inundation mapping and inland water levels. Both PALSAR‐ and ASAR‐based mapping at times were based on suboptimal reference scenes; however, ASAR performance seemed more sensitive to reference‐scene quality and other types of scene variability. Related to water depth, PALSAR and ASAR mapping accuracies tended to be lower when water depths were shallow and increased as water levels decreased below or increased above the ground surface, but this pattern was more pronounced with ASAR. Overall, PALSAR‐based inundation accuracies averaged 84% (n = 160), while ASAR‐based mapping accuracies averaged 62% (n = 245).     

EFFECTS OF DAM IMPOUNDMENT ON THE FLOOD REGIME OF NATURAL FLOODPLAIN COMMUNITIES IN THE UPPER CONNECTICUT RIVER1

ABSTRACT: Understanding the effects of dams on the inundation regime of natural floodplain communities is critical for effective decision making on dam management or dam removal. To test the implications of hydrologic alteration by dams for floodplain natural communities, we conducted a combined field and modeling study along two reaches in the Connecticut River Rapids Macrosite (CRRM), one of the last remaining flowing water sections of the Upper Connecticut River. We surveyed multiple channel cross sections at both locations and concurrently identified and surveyed the elevations of important natural communities, native species of concern, and nonnative invasive species. Using a hydrologic model, HEC‐RAS, we routed estimated pre‐and post‐impoundment discharges of different design recurrence intervals (two year through 100 year floods) through each reach to establish corresponding reductions in elevation and effective wetted perimeter following post‐dam discharge reductions. By comparing (1) the frequency and duration of flooding of these surfaces before and after impoundment and (2) the total area flooded at different recurrence intervals, our goal was to derive a spatially explicit assessment of hydrologic alteration, directly relevant to natural floodplain communities. Post‐impoundment hydrologic alteration profoundly affected the subsequent inundation regime, and this impact was particularly true of higher floodplain terraces. These riparian communities, which were flooded, on average, every 20 to 100 years pre‐impoundment, were predicted to flood at 100 ? 100 year intervals, essentially isolating them completely from riverine influence. At the pre‐dam five to ten year floodplain elevations, we observed smaller differences in predicted flood frequency but substantial differences in the total area flooded and in the average flood duration. For floodplain forests in the Upper Connecticut River, this alteration by impoundment suggests that even if other stresses facing these communities (human development, invasive exotics) were alleviated, this may not be sufficient to restore intact natural communities. More generally, our approach provides a way to combine site specific variables with long term gage records in assessing the restorative potential of dam removal.     

Estimating Floodwater Depths from Flood Inundation Maps and Topography

Information on flood inundation extent is important for understanding societal exposure, water storage volumes, flood wave attenuation, future flood hazard, and other variables. A number of organizations now provide flood inundation maps based on satellite remote sensing. These data products can efficiently and accurately provide the areal extent of a flood event, but do not provide floodwater depth, an important attribute for first responders and damage assessment. Here we present a new methodology and a GIS‐based tool, the Floodwater Depth Estimation Tool (FwDET), for estimating floodwater depth based solely on an inundation map and a digital elevation model (DEM). We compare the FwDET results against water depth maps derived from hydraulic simulation of two flood events, a large‐scale event for which we use medium resolution input layer (10 m) and a small‐scale event for which we use a high‐resolution (LiDAR; 1 m) input. Further testing is performed for two inundation maps with a number of challenging features that include a narrow valley, a large reservoir, and an urban setting. The results show FwDET can accurately calculate floodwater depth for diverse flooding scenarios but also leads to considerable bias in locations where the inundation extent does not align well with the DEM. In these locations, manual adjustment or higher spatial resolution input is required.     

A Framework to Develop Nationwide Flooding Extents Using Climate Models and Assess Forecast Potential for Flood Resilience

The methods used to simulate flood inundation extents can be significantly improved by high‐resolution spatial data captured over a large area. This paper presents a hydraulic analysis methodology and framework to estimate national‐level floodplain changes likely to be generated by climate change. The hydraulic analysis was performed using existing published Federal Emergency Management Agency 100‐year floodplains and estimated 100‐ and 10‐year return period peak flow discharges. The discharges were estimated using climate variables from global climate models for two future growth scenarios: Representative Concentration Pathways 2.6 and 8.5. River channel dimensions were developed based on existing regional United States Geological Survey publications relating bankfull discharges with channel characteristics. Mathematic relationships for channel bankfull topwidth, depth, and side slope to contributing drainage area measured at model cross sections were developed. The proposed framework can be utilized at a national level to identify critical areas for flood risk assessment. Existing hydraulic models at these “hot spots” could be repurposed for near–real‐time flood forecasting operations. Revitalizing these models for use in simulating flood scenarios in near–real time through the use of meteorological forecasts could provide useful information for first responders of flood emergencies.     

INFLUENCES OF SOIL DATASET RESOLUTION ON HYDROLOGIC MODELING1

ABSTRACT: Soils represent a fundamental abiotic parameter in defining the characteristics of an ecosystem. The U.S. Department of Agriculture (USDA) Natural Resources Conservation Service (NRCS) produces the most detailed digital spatial soil datasets that are publicly available. The Soil Survey Geographic (SSURGO) database contains basic attributes for the continuous coverage of soils across the United States. In its standard format, the SSURGO database is incompatible for use within the ArcView Soil and Water Assessment Tool (SWAT). A modified version of the State Soil and Geographic (STATSGO) database is the template soils dataset used by ArcView SWAT. This paper presents the methodology and development of a SSURGO database preprocessor extension for the ArcView SWAT model. A case study for the Upper Sabinal River Watershed near Uvalde, Texas, is given. Results indicate that hydro‐logic output parameter differences occur when comparing the STATSGO and SSURGO database information in the ArcView SWAT model under identical modeling conditions. Specifically, the SSURGO model produced a greater daily mean water yield with evapotranspiration and surface runoff being found consistently lower across the watershed. The most likely causes assigned to this phenomenon were higher percolation and resulting ground water return flow values due to significantly larger saturated hydraulic conductivity values associated with the SSURGO 2.x database.     

A CyberGIS Integration and Computation Framework for High‐Resolution Continental‐Scale Flood Inundation Mapping

We present a Digital Elevation Model‐based hydrologic analysis methodology for continental flood inundation mapping (CFIM), implemented as a cyberGIS scientific workflow in which a 1/3rd arc‐second (10 m) height above nearest drainage (HAND) raster data for the conterminous United States (CONUS) was computed and employed for subsequent inundation mapping. A cyberGIS framework was developed to enable spatiotemporal integration and scalable computing of the entire inundation mapping process on a hybrid supercomputing architecture. The first 1/3rd arc‐second CONUS HAND raster dataset was computed in 1.5 days on the cyberGIS Resourcing Open Geospatial Education and Research supercomputer. The inundation mapping process developed in our exploratory study couples HAND with National Water Model forecast data to enable near real‐time inundation forecasts for CONUS. The computational performance of HAND and the inundation mapping process were profiled to gain insights into the computational characteristics in high‐performance parallel computing scenarios. The establishment of the CFIM computational framework has broad and significant research implications that may lead to further development and improvement of flood inundation mapping methodologies.     

EFFECTS OF STATSGO AND SSURGO AS INPUTS ON SWAT MODEL'S SNOWMELT SIMULATION1

ABSTRACT: Soil data comprise a basic input of SWAT (Soil and Water Assessment Tool) for a watershed application. For watersheds where site specific soil data are unavailable, the two U.S. Department of Agriculture (USDA) soil databases, the State Soil Geographic (STATSGO) and Soil Survey Geographic (SSURGO) databases, may be the best alternatives. Although it has been noted that SWAT models using the STATSGO and SSURGO data may give different simulation results for water, sediment, and agricultural chemical yields, information is scarce on the effects of using these two databases in predicting streamflows that are predominantly generated from melting snow in spring. The objective of this study was to assess the effects of using STATSGO versus SSURGO as an input for the SWAT model's simulation of the streamflows in the upper 45 percent of the Elm River watershed in eastern North Dakota. Designating the model as SWAT‐STATSGO when the STATSGO data were used and SWAT‐SSURGO when the SSURGO data were used, SWAT‐STATSGO and SWAT‐SSURGO were separately calibrated and validated using the observed daily streamflows. The results indicated that SWAT‐SSURGO provided an overall better prediction of the discharges than SWAT‐STATSGO, although both did a good and comparable job of predicting the high streamflows. However, SWAT‐STATSGO predicted the low streamflows more accurately and had a slightly better performance during the validation period. In addition, the discrepancies between the discharges predicted by these two SWAT models tended to be larger at upstream locations than at those farther downstream within the study area.     

Understanding the Large‐Scale Influence of Levees on Floodplain Connectivity Using a Hydrogeomorphic Approach

The widespread construction of levees has reduced river–floodplain connectivity and altered associated fluvial processes in many river systems. Despite the recognition that levees can alter floodplain connectivity, few studies have examined the role of levees in reducing floodplain areas at large watershed scales. This paper explores the application of a hydrogeomorphic floodplain inundation model in the Wabash Basin, located in the Midwestern United States, to assess changes in floodplain area in levee‐protected areas. We evaluate 10‐ and 30‐m topographic resolutions and spatially examine the influence of levees on floodplain area in relation to river network attributes. Generally, floodplains in levee‐protected areas were influenced by topographic resolution, stream order, and elevation details of levees found in topography datasets. We show, when compared to Federal Emergency Management Agency maps, our approach underpredicts floodplain area when using 10‐m resolution topography data but only slightly overpredicts when using 30‐m resolution data. After removing details of levees from topography data, we found changes in floodplain area varied spatially, but basin‐aggregate results changed little compared to topography datasets that contain levees, though larger floodplain areas were produced in some regions where levees were removed. This work contributes to a growing research emphasis on using hydrogeomorphic floodplain models to understand floodplain disconnectivity.     

Using Steady‐State Backwater Analysis to Predict Inundated Area from National Water Model Streamflow Simulations

National Water Model (NWM) simulates the hydrologic cycle and produces streamflow forecasts for 2.7 million reaches in the National Hydrography Dataset for continental United States (U.S.). NWM uses Muskingum–Cunge channel routing, which is based on the continuity equation. However, the momentum equation also needs to be considered to obtain more accurate estimates of streamflow and stage in rivers, especially for applications such as flood‐inundation mapping. Here, we used a steady‐state backwater version of Simulation Program for River NeTworks (SPRNT) model. We evaluated SPRNT’s and NWM’s abilities to predict inundated area for the record flood of Hurricane Matthew in October 2016. The Neuse River experienced record‐breaking floods and was well‐documented by U.S. Geological Survey. Streamflow simulations from NWM retrospective analysis were used as input for the SPRNT simulation. Retrospective NWM discharge predictions were converted to stage. The stages (from both SPRNT and NWM) were utilized to produce flood‐inundation maps using the Height Above Nearest Drainage method which uses the local relative heights to find out the local draining potentials and provide spatial representation of inundated area. The inundated‐area accuracies for NWM and SPRNT (based on comparison to a remotely sensed dataset) were 65.1% and 67.6%, respectively. These results show using steady‐state SPRNT results in a modest improvement of inundation‐forecast accuracy compared to NWM.     

Capturing LiDAR‐Derived Hydrologic Spatial Parameters to Evaluate Playa Wetlands

The digital elevation model data from traditional stereo photogrammetric methods are inadequate in providing accurate vertical parameters to feed hydrologic models for low‐lying, extremely flat areas. High‐resolution light detection and ranging (LiDAR) data provide the robust capability of capturing small variations in low‐relief playa wetlands. The Rainwater Basin in south‐central Nebraska includes a complex of seasonally shallow playa wetlands that attract millions of migratory waterfowl every spring and fall. This research focuses on the development of a procedure with applicable protocols to produce LiDAR‐derived three‐dimensional wetland maps and to extract the critical surface parameters (i.e., watershed boundaries, flow direction, flow accumulation, and drainage lines) for playa wetlands. The topo‐hydrologic conditions of playa wetlands were evaluated at the watershed level. The results show that in the Rainwater Basin, 70.7% of the historic hydric soil footprints identified in the Soil Survey Geographic (SSURGO) database were not functioning as topographically depressional wetlands. This finding was confirmed by a recent five‐year Annual Habit Survey showing that 69.8% of the historic hydric soil footprints did not function during the spring migratory bird seasons between 2004 and 2009. The majority of playa wetlands' topographic conditions have been substantially changed and the SSURGO data cannot fully reflect current topographic reality in the Rainwater Basin.     

The use of Landsat multispectral scanner data for the analysis and management of flooding on the river Severn,England

Remote sensing has emerged as one of the major techniques for the analysis and delineation of large floods. This analysis can provide data invaluable for the hydrological management of large river systems. A need for information on the extent of floodplain inundation for the lower reaches of the largest river in the UK was met by a search through Landsat images of floods and the analysis of the best example recorded. Automated classification of the Landsat imagery of this flood on the river Severn in 1977 was used to provide estimates of the extent and spatial distribution of inundation. Flood images were generated using the Plessey IDP 3000 image processor, and the maps derived accorded well with aerial photography and qualitative flood information. Three distinct floodplain environments were delineated and flood images produced by different spectral bands compared. Specific questions prompted by flood hazard management and concerning the processes and extent of flooding were answered by the Landsat data analysis. Management of the flood risk of large rivers is expensive and remote sensing data is a relatively cheap and effective way of monitoring control works and providing data for the prediction of the effects of future hydrological works. Remote sensing is a practical way in which spatial information concerning the behavior of large dynamic systems can be obtained both quickly and relatively cheaply.     

Effects of soil data resolution on SWAT model stream flow and water quality predictions

The prediction accuracy of agricultural nonpoint source pollution models such as Soil and Water Assessment Tool (SWAT) depends on how well model input spatial parameters describe the characteristics of the watershed. The objective of this study was to assess the effects of different soil data resolutions on stream flow, sediment and nutrient predictions when used as input for SWAT. SWAT model predictions were compared for the two US Department of Agriculture soil databases with different resolution, namely the State Soil Geographic database (STATSGO) and the Soil Survey Geographic database (SSURGO). Same number of sub-basins was used in the watershed delineation. However, the number of HRUs generated when STATSGO and SSURGO soil data were used is 261 and 1301, respectively. SSURGO, with the highest spatial resolution, has 51 unique soil types in the watershed distributed in 1301 HRUs, while STATSGO has only three distributed in 261 HRUS. As a result of low resolution STATSGO assigns a single classification to areas that may have different soil types if SSURGO were used. SSURGO included Hydrologic Response Units (HRUs) with soil types that were generalized to one soil group in STATSGO. The difference in the number and size of HRUs also has an effect on sediment yield parameters (slope and slope length). Thus, as a result of the discrepancies in soil type and size of HRUs stream flow predicted was higher when SSURGO was used compared to STATSGO. SSURGO predicted less stream loading than STATSGO in terms of sediment and sediment-attached nutrients components, and vice versa for dissolved nutrients. When compared to mean daily measured flow, STATSGO performed better relative to SSURGO before calibration. SSURGO provided better results after calibration as evaluated by R(2) value (0.74 compared to 0.61 for STATSGO) and the Nash-Sutcliffe coefficient of Efficiency (NSE) values (0.70 and 0.61 for SSURGO and STATSGO, respectively) although both are in the same satisfactory range. Modelers need to weigh the benefits before selecting the type of data resolution they are going to use depending on the watershed size and level of accuracy required because more effort is required to prepare and calibrate the model when a fine resolution soil data is used.