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.     

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.     

Identifying riparian sinks for watershed nitrate using soil surveys

The capacity of riparian zones to serve as critical control locations for watershed nitrogen flux varies with site characteristics. Without a means to stratify riparian zones into different levels of ground water nitrate removal capacity, this variability will confound spatially explicit source-sink models of watershed nitrate flux and limit efforts to target riparian restoration and management. We examined the capability of SSURGO (1:15 840 Soil Survey Geographic database) map classifications (slope class, geomorphology, and/or hydric soil designation) to identify riparian sites with high capacity for ground water nitrate removal. The study focused on 100 randomly selected riparian locations in a variety of forested and glaciated settings within Rhode Island. Geomorphic settings included till, outwash, and organic/alluvial deposits. We defined riparian zones with "high ground water nitrate removal capacity" as field sites possessing both >10 m of hydric soil width and an absence of ground water surface seeps. SSURGO classification based on a combination of geomorphology and hydric soil status created two functionally distinct sets of riparian sites. More than 75% of riparian sites classified by SSURGO as organic/alluviumhydric or as outwash-hydric had field attributes that suggest a high capacity for ground water nitrate removal. In contrast, >85% of all till sites and nonhydric outwash sites had field characteristics that minimize the capacity for ground water nitrate removal. Comparing the STATSGO and SSURGO databases for a 64000-ha watershed, STATSGO grossly under-represented critical riparian features. We conclude that the SSURGO database can provide modelers and managers with important insights into riparian zone nitrogen removal potential.     

WATERSHED SCALE MODELING OF CRITICAL SOURCE AREAS OF RUNOFF GENERATION AND PHOSPHORUS TRANSPORT1

ABSTRACT: A curve number based model, Soil and Water Assessment Tool (SWAT), and a physically based model, Soil Moisture Distribution and Routing (SMDR), were applied in a headwater watershed in Pennsylvania to identify runoff generation areas, as runoff areas have been shown to be critical for phosphorus management. SWAT performed better than SMDR in simulating daily streamflows over the four‐year simulation period (Nash‐Sutcliffe coefficient: SWAT, 0.62; SMDR, 0.33). Both models varied streamflow simulations seasonally as precipitation and watershed conditions varied. However, levels of agreement between simulated and observed flows were not consistent over seasons. SMDR, a variable source area based model, needs further improvement in model formulations to simulate large peak flows as observed. SWAT simulations matched the majority of observed peak flow events. SMDR overpredicted annual flow volumes, while SWAT underpredicted the same. Neither model routes runoff over the landscape to water bodies, which is critical to surface transport of phosphorus. SMDR representation of the watershed as grids may allow targeted management of phosphorus sources. SWAT representation of fields as hydrologic response units (HRUs) does not allow such targeted management.     

Hydrologic and atrazine simulation of the Cedar Creek Watershed using the SWAT model

One of the major factors contributing to surface water contamination in agricultural areas is the use of pesticides. The Soil and Water Assessment Tool (SWAT) is a hydrologic model capable of simulating the fate and transport of pesticides in an agricultural watershed. The SWAT model was used in this study to estimate stream flow and atrazine (2-chloro-4-(ethylamino)-6-(isopropylamino)-s-triazine) losses to surface water in the Cedar Creek Watershed (CCW) within the St. Joseph River Basin in northeastern Indiana. Model calibration and validation periods consisted of five and two year periods, respectively. The National Agricultural Statistics Survey (NASS) 2001 land cover classification and the Soil Survey Geographic (SSURGO) database were used as model input data layers. Data from the St. Joseph River Watershed Initiative and the Soil and Water Conservation Districts of Allen, Dekalb, and Noble counties were used to represent agricultural practices in the watershed which included the type of crops grown, tillage practices, fertilizer, and pesticide application rates. Model results were evaluated based on efficiency coefficient values, standard statistical measures, and visual inspection of the measured and simulated hydrographs. The Nash and Sutcliffe model efficiency coefficients (E(NS)) for monthly and daily stream flow calibration and validation ranged from 0.51 to 0.66. The E(NS) values for atrazine calibration and validation ranged from 0.43 to 0.59. All E(NS) values were within the range of acceptable model performance standards. The results of this study indicate that the model is an effective tool in capturing the dynamics of stream flow and atrazine concentrations on a large-scale agricultural watershed in the midwestern USA.     

Topographic Effects on Soil Organic Carbon in Louisiana Watersheds

Terrestrial carbon storage is influenced by a number of environmental factors, among which topographic and geomorphological features are of special significance. This study was designed to examine the relationships of soil organic carbon (SOC) density to various terrain parameters and watershed characteristics across Louisiana, USA. A polygon data set of 484 watersheds and 12 river drainage basins for Louisiana was used to form the landscape units. SOC densities were calculated for each soil map unit using the State Soil Geographic (STATSGO) database. Average drainage densities and average slopes at watershed and basin scales were quantified with the 1:24 K Digital Elevation Models (DEM) data, and the Louisiana hydrographic water features. Correlation and regression analyses were performed to determine relationships among drainage density, slope, elevation, and SOC. The study found an average watershed drainage density of 1.6 km/km² and an average watershed slope of 2.9 degrees in Louisiana. The results revealed that SOC density at both watershed and basin scales was closely related to drainage density, slope, and elevation. SOC density was positively correlated with watershed drainage density, but negatively correlated with watershed slope gradient and elevation. Regression models were developed for predicting SOC density at watershed and basin scales, obtaining regression coefficients (r ²) ranging from 0.43 to 0.83. The study showed that estimation of SOC at watershed and drainage basin scales combining DEM data can be a feasible approach to improve the understanding of the relationships among SOC, topographic, and geomorphological features.     

EPIC modeling of soil organic carbon sequestration in croplands of Iowa

Depending on management, soil organic carbon (SOC) is a potential source or sink for atmospheric CO(2). We used the EPIC model to study impacts of soil and crop management on SOC in corn (Zea mays L.) and soybean (Glycine max L. Merr.) croplands of Iowa. The National Agricultural Statistics Service crops classification maps were used to identify corn-soybean areas. Soil properties were obtained from a combination of SSURGO and STATSGO databases. Daily weather variables were obtained from first order meteorological stations in Iowa and neighboring states. Data on crop management, fertilizer application and tillage were obtained from publicly available databases maintained by the NRCS, USDA-Economic Research Service (ERS), and Conservation Technology Information Center. The EPIC model accurately simulated state averages of crop yields during 1970-2005 (R(2) = 0.87). Simulated SOC explained 75% of the variation in measured SOC. With current trends in conservation tillage adoption, total stock of SOC (0-20 cm) is predicted to reach 506 Tg by 2019, representing an increase of 28 Tg with respect to 1980. In contrast, when the whole soil profile was considered, EPIC estimated a decrease of SOC stocks with time, from 1835 Tg in 1980 to 1771 Tg in 2019. Hence, soil depth considered for calculations is an important factor that needs further investigation. Soil organic C sequestration rates (0-20 cm) were estimated at 0.50 to 0.63 Mg ha(-1) yr(-1) depending on climate and soil conditions. Overall, combining land use maps with EPIC proved valid for predicting impacts of management practices on SOC. However, more data on spatial and temporal variation in SOC are needed to improve model calibration and validation.     

ARCGIS‐SWAT: A GEODATA MODEL AND GIS INTERFACE FOR SWAT1

This paper presents ArcGIS‐SWAT, a geodata model and geographic information system (GIS) interface for the Soil and Water Assessment Tool (SWAT). The ArcGIS‐SWAT data model is a system of geodatabases that store SWAT geographic, numeric, and text input data and results in an organized fashion. Thus, it is proposed that a single and comprehensive geodatabase be used as the repository of a SWAT simulation. The ArcGIS‐SWAT interface uses programming objects that conform to the Component Object Model (COM) design standard, which facilitate the use of functionality of other Windows‐based applications within ArcGIS‐SWAT. In particular, the use of MS Excel and MATLAB functionality for data analysis and visualization of results is demonstrated. Likewise, it is proposed to conduct hydrologic model integration through the sharing of information with a not‐model‐specific hub data model where information common to different models can be stored and from which it can be retrieved. As an example, it is demonstrated how the Hydrologic Modeling System (HMS) ‐ a computer application for flood analysis ‐ can use information originally developed by ArcGIS‐SWAT for SWAT. The application of ArcGIS‐SWAT to the Seco Creek watershed in Texas is presented.     

Augmenting Watershed Model Calibration with Incorporation of Ancillary Data Sources and Qualitative Soft Data Sources

Watershed simulation models such as the Soil & Water Assessment Tool (SWAT) can be calibrated using “hard data” such as temporal streamflow observations; however, users may find upon examination of model outputs, that the calibrated models may not reflect actual watershed behavior. Thus, it is often advantageous to use “soft data” (i.e., qualitative knowledge such as expected denitrification rates that observed time series do not typically exist) to ensure that the calibrated model is representative of the real world. The primary objective of this study is to evaluate the efficacy of coupling SWAT‐Check (a post‐evaluation framework for SWAT outputs) and IPEAT‐SD (Integrated Parameter Estimation and Uncertainty Analysis Tool‐Soft & hard Data evaluation) to constrain the bounds of soft data during SWAT auto‐calibration. IPEAT‐SD integrates 59 soft data variables to ensure SWAT does not violate physical processes known to occur in watersheds. IPEAT‐SD was evaluated for two case studies where soft data such as denitrification rate, nitrate attributed from subsurface flow to total discharge ratio, and total sediment loading were used to conduct model calibration. Results indicated that SWAT model outputs may not satisfy reasonable soft data responses without providing pre‐defined bounds. IPEAT‐SD provides an efficient and rigorous framework for users to conduct future studies while considering both soft data and traditional hard information measures in watershed modeling.     

INTEGRATION OF A BASIN-SCALE WATER QUALITY MODEL WITH GIS1

ABSTRACT: Geographic Information Systems (GIS) have been successfully integrated with distributed parameter, single-event, water quality models such as AGNPS (AGricultural NonPoint Source) and ANSWERS (Areal Nonpoint Source Watershed Environmental Response Simulation). These linkages proved to be an effective way to collect, manipulate, visualize, and analyze the input and output date of water quality models. However, for continuous-time, basin large-scale water quality models, collecting and manipulating the input data are more time-consuming and cumbersome due to the method of disaggregation (subdivisions are based on topographic boundaries). SWAT (Soil and Water Assessment Tool), a basin-scale water quality model, was integrated with a GIS to extract input data for modeling a basin. This paper discusses the detailed development of the integration of the SWAT water quality model with GRASS (Geographic Resources Analysis Support System) GIS, along with an application and advantages. The integrated system was applied to simulated a 114 sq. km upper portion of the Seco Creek Basin by subdividing it into 37 subbasins. The average monthly predicted streamflw is in agreement with measured monthly streamflw values.     

Spatial and Temporal Evaluation of Hydrological Response to Climate and Land Use Change in Three South Dakota Watersheds

This study analyzed changes in hydrology between two recent decades (1980s and 2010s) with the Soil and Water Assessment Tool (SWAT) in three representative watersheds in South Dakota: Bad River, Skunk Creek, and Upper Big Sioux River watersheds. Two SWAT models were created over two discrete time periods (1981‐1990 and 2005‐2014) for each watershed. National Land Cover Datasets 1992 and 2011 were, respectively, ingested into 1981‐1990 and 2005‐2014 models, along with corresponding weather data, to enable comparison of annual and seasonal runoff, soil water content, evapotranspiration (ET), water yield, and percolation between these two decades. Simulation results based on the calibrated models showed that surface runoff, soil water content, water yield, and percolation increased in all three watersheds. Elevated ET was also apparent, except in Skunk Creek watershed. Differences in annual water balance components appeared to follow changes in land use more closely than variation in precipitation amounts, although seasonal variation in precipitation was reflected in seasonal surface runoff. Subbasin‐scale spatial analyses revealed noticeable increases in water balance components mostly in downstream parts of Bad River and Skunk Creek watersheds, and the western part of Upper Big Sioux River watershed. Results presented in this study provide some insight into recent changes in hydrological processes in South Dakota watersheds. Editor's note: This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

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.     

PRISM vs. CFSR Precipitation Data Effects on Calibration and Validation of SWAT Models

Precipitation is one of the most important drivers in watershed models. Our objective was to compare two sources of interpolated precipitation data in terms of their effect on calibration and validation of two Soil and Water Assessment Tool (SWAT) models. One model was a suburban watershed in metropolitan Atlanta, Georgia. The precipitation sources were Parameter‐elevation Relationships on Independent Slopes Model (PRISM) data on a 4‐km grid and climate forecast system reanalysis (CFSR) data on a 38‐km grid. The PRISM data resulted in a better fit to the calibration data (Nash Sutcliffe efficiency [NSE] = 0.64, Kling‐Gupta efficiency [KGE] = 0.74, p‐factor = 0.84, and r‐factor = 0.43) than the CFSR data (NSE = 0.47, KGE = 0.53, p‐factor = 0.67, and r‐factor = 0.39). Validation results were similar. Sensitive parameters were similar in both the PRISM and CFSR models, but fitted values indicated more rapid groundwater flow to the streams with the PRISM data. The same comparison was made in the Big Creek watershed located approximately 1,000 km away, in central Louisiana. Results were similar with a more responsive groundwater system indicating PRISM data may produce better predictions of streamflow because of a more accurate estimate of rainfall within a watershed or because of a denser grid. Our study implies PRISM is providing a better estimate than CFSR of precipitation within a watershed when rain gauge data are not available, resulting in more accurate simulations of streamflows at the watershed outlet. Editor's note: This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Development of a SWAT Patch for Better Estimation of Sediment Yield in Steep Sloping Watersheds1

Abstract: The watershed scale Soil and Water Assessment Tool (SWAT) model divides watersheds into smaller subwatersheds for simulation of rainfall‐runoff and sediment loading at the field level and routing through stream networks. Typically, the SWAT model first needs to be calibrated and validated for accurate estimation through adjustment of sensitive input parameters (i.e., Curve Number values, USLE P, slope and slope‐length, and so on). However, in some instances, SWAT‐simulated results are greatly affected by the watershed delineation and Digital Elevation Models (DEM) cell size. In this study, the SWAT ArcView GIS Patch II was developed for steep sloping watersheds, and its performance was evaluated for various threshold values and DEM cell size scenarios when delineating subwatersheds using the SWAT model. The SWAT ArcView GIS Patch II was developed using the ArcView GIS Avenue program and Spatial Analyst libraries. The SWAT ArcView GIS Patch II improves upon the SWAT ArcView GIS Patch I because it reflects the topographic factor in calculating the field slope‐length of Hydrologic Response Units in the SWAT model. The simulated sediment value for 321 subwatersheds (watershed delineation threshold value of 25 ha) is greater than that for 43 subwatersheds (watershed delineation threshold value of 200 ha) by 201% without applying the SWAT ArcView GIS Patch II. However, when the SWAT ArcView GIS Patch II was applied, the difference in simulated sediment yield decreases for the same scenario (i.e., difference in simulated sediment with 321 subwatersheds and 43 subwatersheds) was 12%. The simulated sediment value for DEM cell size of 50 m is greater than that for DEM cell size of 10 m by 19.8% without the SWAT ArcView GIS Patch II. However, the difference becomes smaller (3.4% difference) between 50 and 10 m with the SWAT ArcView GIS Patch II for the DEM scenarios. As shown in this study, the SWAT ArcView GIS Patch II can reduce differences in simulated sediment values for various watershed delineation and DEM cell size scenarios. Without the SWAT ArcView GIS Patch II, variations in the SWAT‐simulated results using various watershed delineation and DEM cell size scenarios could be greater than those from input parameter calibration. Thus, the results obtained in this study show that the SWAT ArcView GIS Patch II should be used when simulating hydrology and sediment yield for steep sloping watersheds (especially if average slope of the subwatershed is >25%) for more accurate simulation of hydrology and sediment using the SWAT model. The SWAT ArcView GIS Patch II is available at http://www.EnvSys.co.kr/~swat for free download.     

HYDROLOGICAL MODELING OF THE IROQUOIS RIVER WATERSHED USING HSPF AND SWAT1

ABSTRACT: The performance of two popular watershed scale simulation models — HSPF and SWAT — were evaluated for simulating the hydrology of the 5,568 km² Iroquois River watershed in Illinois and Indiana. This large, tile drained agricultural watershed provides distinctly different conditions for model comparison in contrast to previous studies. Both models were calibrated for a nine‐year period (1987 through 1995) and verified using an independent 15‐year period (1972 through 1986) by comparing simulated and observed daily, monthly, and annual streamflow. The characteristics of simulated flows from both models are mostly similar to each other and to observed flows, particularly for the calibration results. SWAT predicts flows slightly better than HSPF for the verification period, with the primary advantage being better simulation of low flows. A noticeable difference in the models' hydrologic simulation relates to the estimation of potential evapotranspiration (PET). Comparatively low PET values provided as input to HSPF from the BASINS 3.0 database may be a factor in HSPF's overestimation of low flows. Another factor affecting baseflow simulation is the presence of tile drains in the watershed. HSPF parameters can be adjusted to indirectly account for the faster subsurface flow associated with tile drains, but there is no specific tile drainage component in HSPF as there is in SWAT. Continued comparative studies such as this, under a variety of hydrologic conditions and watershed scales, provide needed guidance to potential users in model selection and application.     

Modeling the Spatially Varying Water Balance Processes in a Semiarid Mountainous Watershed of Idaho1

Stratton, Benjamin T., Venakataramana Sridhar, Molly M. Gribb, James P. McNamara, and Balaji Narasimhan, 2009. Modeling the Spatially Varying Water Balance Processes in a Semiarid Mountainous Watershed of Idaho. Journal of the American Water Resources Association (JAWRA) 45(6):1390‐1408. Abstract: The distributed Soil Water Assessment Tool (SWAT) hydrologic model was applied to a research watershed, the Dry Creek Experimental Watershed, near Boise Idaho to investigate its water balance components both temporally and spatially. Calibrating and validating SWAT is necessary to enable our understanding of the water balance components in this semiarid watershed. Daily streamflow data from four streamflow gages were used for calibration and validation of the model. Monthly estimates of streamflow during the calibration phase by SWAT produced satisfactory results with a Nash Sutcliffe coefficient of model efficiency 0.79. Since it is a continuous simulation model, as opposed to an event‐based model, it demonstrated the limited ability in capturing both streamflow and soil moisture for selected rain‐on‐snow (ROS) events during the validation period between 2005 and 2007. Especially, soil moisture was generally underestimated compared with observations from two monitoring pits. However, our implementation of SWAT showed that seasonal and annual water balance partitioning of precipitation into evapotranspiration, streamflow, soil moisture, and drainage was not only possible but closely followed the trends of a typical semiarid watershed in the intermountain west. This study highlights the necessity for better techniques to precisely identify and drive the model with commonly observed climatic inversion‐related snowmelt or ROS weather events. Estimation of key parameters pertaining to soil (e.g., available water content and saturated hydraulic conductivity), snow (e.g., lapse rates, melting), and vegetation (e.g., leaf area index and maximum canopy index) using additional field observations in the watershed is critical for better prediction.     

Assessing the Impact of Areal Precipitation Input on Streamflow Simulations Using the SWAT Model1

Masih Ilyas, Shreedhar Maskey, Stefan Uhlenbrook, and Vladimir Smakhtin, 2011. Assessing the Impact of Areal Precipitation Input on Streamflow Simulations Using the SWAT Model. Journal of the American Water Resources Association (JAWRA) 47(1):179‐195. DOI: 10.1111/j.1752‐1688.2010.00502.x Abstract: Reduction of input uncertainty is a challenge in hydrological modeling. The widely used model Soil Water Assessment Tool (SWAT) uses the data of a precipitation gauge nearest to the centroid of each subcatchment as an input for that subcatchment. This may not represent overall catchment precipitation conditions well. This paper suggests an alternative – using areal precipitation obtained through interpolation. The effectiveness of this alternative is evaluated by comparing its simulations with those based on the standard SWAT precipitation input procedure. The model is applied to mountainous semiarid catchments in the Karkheh River basin, Iran. The model performance is evaluated at daily, monthly, and annual scales by using a number of performance indicators at 15 streamflow gauging stations each draining an area in the range of 590‐42,620 km². The comparison suggests that the use of areal precipitation improves model performance particularly in small subcatchments in the range of 600‐1,600 km². The modified areal precipitation input results in increased reliability of simulated streamflows in the areas of low rain gauge density. Both precipitation input methods result in reasonably good simulations for larger catchments (over 5,000 km²). The use of areal precipitation input improves the accuracy of simulated streamflows with spatial resolution and density of rain gauges having significant impact on results.     

Spatial extrapolation of soil characteristics using whole-soil particle size distributions

Soils support ecosystem functions such as plant growth and water quality because of certain physical, chemical, and biological properties. These properties have been studied at different spatial scales, including point scales to satisfy basic research needs, and regional scales to satisfy monitoring needs. Recently, soil property data for the entire USA have become available in the State Soil Geographic Data Base (STATSGO), which is appropriate for regional-scale research. We analyzed and created models of STATSGO data in this study to serve as a research tool, for example, for linking the soil to regional water quality monitoring data in our companion paper. Map units in STATSGO define geographic land areas by soil characteristics (SCs) of similar soil series. We selected 27 SCs that influenced water properties (in varying degrees), aggregated the layer and component SCs to map unit SCs, and used SCs to calculate relationships among map units. The relationships were defined by equations of conditional mean for the qth SC (SCq), while using the remaining 26 SCs as predictors. The relative standard errors for 22 of the 27 SCs were less than 10%, and less than 22% for the remaining five. We conclude that spatial extrapolation of SCs is feasible and the procedures are a first step toward extrapolating information across a region using SC-water property relationships. Although our procedure is for regional scale monitoring, it is also applicable to finer spatial scales commensurate with available soil data.     

WATER QUALITY MODEL OUTPUT UNCERTAINTY AS AFFECTED BY SPATIAL RESOLUTION OF INPUT DATA1

ABSTRACT: Resolution of the input GIS data used to parameterize distributed‐parameter hydrologic/water quality models may affect uncertainty in model outputs and impact the subsequent application of model results in watershed management. In this study we evaluated the impact of varying spatial resolutions of DEM, land use, and soil data (30 × 30 m, 100 × 100 m, 150 × 150 m, 200 × 200 m, 300 × 300 m, 500 × 500 m, and 1,000 × 1,000 m) on the uncertainty of SWAT predicted flow, sediment, NO₃‐N, and TP transport. Inputs included measured hydrologic, meteorological, and watershed characteristics as well as water quality data from the Moores Creek watershed in Washington County, Arkansas. The SWAT model output was most affected by input DEM data resolution. A coarser DEM data resolution resulted in decreased representation of watershed area and slope and increased slope length. Distribution of pasture, forest, and urban areas within the watershed was significantly affected at coarser resolution of land use and resulted in significant uncertainty in predicted sediment, NO₃‐N, and TP output. Soils data resolution had no significant effect on flow and NO₃‐N predictions; however, sediment was overpredicted by 26 percent, and TP was underpredicted by 26 percent at 1,000 m resolution. This may be due to change in relative distribution of various hydrologic soils groups (HSGs) in the watershed. Minimum resolution for input GIS data to achieve less than 10 percent model output error depended upon the output variable of interest. For flow, sediment, NO₃‐N, and TP predictions, minimum DEM data resolution should range from 30 to 300 m, whereas minimum land use and soils data resolution should range from 300 to 500 m.