HYSTAR: Hydrology and Sediment Transport Simulation using Time-Area Method

• Main Author: Her, Young Gu
• Other Authors: Biological Systems Engineering
• Format: Others
• Published: Virginia Tech 2014
• Subjects: margin of safety; critical source area; calibration; uncertainty analysis; sediment transport; flow routing; hydrology; distributed watershed modeling; time-area method
• Online Access: http://hdl.handle.net/10919/27503; http://scholar.lib.vt.edu/theses/available/etd-05022011-005324/

A distributed approach can improve functionality of H/WQ (Hydrology and Water Quality) modeling by facilitating a way to explicitly incorporate spatial characteristics of a watershed into the model. The time-area approach, with its intuitive and inherently distributed concept, provides a simple method to simulate runoff mechanisms. This study developed a distributed model based on the time-area approach with the goal of improved utility and efficiency in H/WQ modeling. Uncertainty is always introduced into watershed modeling because of imperfect knowledge and scale dependant spatial heterogeneity and temporal variability. Uncertainty analysis can provide a modeler, policy maker, and stakeholder with reliability information, better understanding, and better communication about the modeling results. This study quantified uncertainty of the model parameter and output through uncertainty analysis in order to assess risk in watershed management. The main goal of this study was to develop a hydrology and sediment transport model capable of routing overland flow using a time-area concept and providing reliability of the modeling results in a probabilistic manner through uncertainty analysis. The HYSTAR (HYdrology and Sediment transport simulation using Time-ARea method) model incorporates a modified Curve Number (CN) method and the newly devised time-area routing method to estimate runoff. HYSTAR is capable of simulating direct runoff, base flow, soil moisture, and sediment load in a distributed manner and in an hourly time step. In the model, the modified CN and a continuity equation are used to calculate infiltration of the routed runoff as well as rainfall on every overland cell. The effective direct runoff volume is distributed over downstream areas using the newly developed routing method. A direct runoff hydrograph is constructed directly through the discrete convolution of the time-area histogram and the effective direct runoff volume map without employing a unit hydrograph. In addition, sediment transport is simulated using the routing method and the sediment transport capacity approach without using a delivery ratio. The sensitivity analysis found that the CN and root zone depth were the most critical parameters for runoff simulation with HYSTAR. The model provided acceptable performance in predicting runoff and sediment load of a subwatershed of the Owl Run Watershed (ORD) with the Nash-Sutcliffe efficiency coefficient and coefficient of determination greater than 0.5. However, it failed to reproduce runoff for a subwatershed of Polecat Creek Watershed (PCA), where data show that runoff is not immediately responsive to rainfall. Uncertainty analysis revealed that the confidence intervals of the simulated monthly runoff and sediment load corresponded to 9.7 % and 10.2 % of their averages, respectively, at a significance level of 0.05. In addition, the average ranges of variation created by the Digital Elevation Model (DEM) and National Land Cover Data (NLCD) errors in the simulated monthly runoff and sediment load were equivalent to 7.5 % and 15.9 % of the average of their calibrated values, respectively. Based on the uncertainty analysis results, the Margin of Safety (MOS) of Total Maximum Daily Load (TMDL) were explicitly quantified as corresponding to 7.0 % and 21.3 % of the average of the simulated runoff and sediment load for ORD at significance level of 0.05. In conclusion, the HYSTAR model provided a new way to explicitly simulate runoff and sediment load of a watershed in a distributed manner. The approach developed here retains the simplicity of a unit hydrograph approach without employing numerical methods. Uncertainty analysis found that parameter uncertainty had greater impact on the model output than did expected Geographic Information System (GIS) data errors. In addition, the impact of the topographic data error on the model output was greater than was that of the land cover data error. Finally, this study provided a proof that a 5 to 10 % MOS that many TMDL studies consider underestimates modeling uncertainty. === Ph. D.