Effects of calibration on L-THIA GIS runoff and pollutant estimation

Urbanization can result in alteration of a watershed's hydrologic response and water quality. To simulate hydrologic and water quality impacts of land use changes, the Long-Term Hydrologic Impact Assessment (L-THIA) system has been used. The L-THIA system estimates pollutant loading based on direct runoff quantity and land use based pollutant coefficients. The accurate estimation of direct runoff is important in assessing water quality impacts of land use changes. An automated program was developed to calibrate the L-THIA model using the millions of curve number (CN) combinations associated with land uses and hydrologic soil groups. L-THIA calibration for the Little Eagle Creek (LEC) watershed near Indianapolis, Indiana was performed using land use data for 1991 and daily rainfall data for six months of 1991 (January 1-June 30) to minimize errors associated with use of different temporal land use data and rainfall data. For the calibration period, the Nash-Sutcliffe coefficient was 0.60 for estimated and observed direct runoff. The calibrated CN values were used for validation of the model for the same year (July 1-December 31), and the Nash-Sutcliffe coefficient was 0.60 for estimated and observed direct runoff. The Nash-Sutcliffe coefficient was 0.52 for January 1, 1991 to December 31, 1991 using uncalibrated CN values. As shown in this study, the use of better input parameters for the L-THIA model can improve accuracy. The effects on direct runoff and pollutant estimation of the calibrated CN values in the L-THIA model were investigated for the LEC. Following calibration, the estimated average annual direct runoff for the LEC watershed increased by 34%, total nitrogen by 24%, total phosphorus by 22%, and total lead by 43%. This study demonstrates that the L-THIA model should be calibrated and validated prior to application in a particular watershed to more accurately assess the effects of land use changes on hydrology and water quality.     

EQUIVALENCE BETWEEN TOPMODEL AND THE NRCS CURVE NUMBER METHOD IN PREDICTING VARIABLE RUNOFF SOURCE AREAS1

An equivalence is proposed between two rainfall‐runoff methods with a long history of use in the United States and Europe. In watersheds where variable source areas dominate runoff, the two methods can have comparable probability distribution functions of moisture deficit, and therefore predict similar saturated runoff source areas. A novel approach is introduced to determine the S parameter in the Natural Resources Conservation Service (NRCS) method. This approach constrains S by the physical soil and topography characteristics of the catchment and depth to water table. The NRCS curve number method is at the core of many rainfall‐runoff models in hydrology. As a simple lumped parameter method, it is often scrutinized because it is not obvious how to derive S from catchment hydromorphological characteristics. The novel approach provides a clear physical meaning for S, allowing better estimation of this parameter in humid shallow water table environments where the variable source area can be the dominant runoff mechanism.     

Assessing Watershed-Scale, Long-Term Hydrologic Impacts of Land-Use Change Using a GIS-NPS Model

Land-use change, dominated by an increase in urban/impervious areas, has a significant impact on water resources. This includes impacts on nonpoint source (NPS) pollution, which is the leading cause of degraded water quality in the United States. Traditional hydrologic models focus on estimating peak discharges and NPS pollution from high-magnitude, episodic storms and successfully address short-term, local-scale surface water management issues. However, runoff from small, low-frequency storms dominates long-term hydrologic impacts, and existing hydrologic models are usually of limited use in assessing the long-term impacts of land-use change. A long-term hydrologic impact assessment (L-THIA) model has been developed using the curve number (CN) method. Long-term climatic records are used in combination with soils and land-use information to calculate average annual runoff and NPS pollution at a watershed scale. The model is linked to a geographic information system (GIS) for convenient generation and management of model input and output data, and advanced visualization of model results. The L-THIA/NPS GIS model was applied to the Little Eagle Creek (LEC) watershed near Indianapolis, Indiana, USA. Historical land-use scenarios for 1973, 1984, and 1991 were analyzed to track land-use change in the watershed and to assess impacts on annual average runoff and NPS pollution from the watershed and its five subbasins. For the entire watershed between 1973 and 1991, an 18% increase in urban or impervious areas resulted in an estimated 80% increase in annual average runoff volume and estimated increases of more than 50% in annual average loads for lead, copper, and zinc. Estimated nutrient (nitrogen and phosphorus) loads decreased by 15% mainly because of loss of agricultural areas. The L-THIA/NPS GIS model is a powerful tool for identifying environmentally sensitive areas in terms of NPS pollution potential and for evaluating alternative land use scenarios for NPS pollution management.     

AN IRRIGATION SCHEDULING MODEL WHICH INCORPORATES RAINFALL PREDICTIONS1

ABSTRACT In humid areas appreciable amounts of rainfall complicate irrigation scheduling. This rainfall tends to give supplemental water application a low priority. As a result irrigation may be delayed until there is not enough time to cover the crop area before some drought damage occurs. To improve the management of irrigation systems, a scheduling model has been developed. The model's water application decisions incorporate climatological records, soil-plant data, current pan evaporation and rainfall, the number of fields to be irrigated, and 5-day weather forecasts. The model updates the soil moisture conditions, predicts impending water depletion, and if supplemental water is needed both the field priority and amount to be applied is indicated for each of the next 5 days. Errors introduced through the use of forecasts and long-term pan evaporation records have been slight because of the tri-weekly updating. Also natural rains which restore the root zone to maximum water holding capacity prevent long-term bias.     

The effects of large-scale afforestation and climate change on water allocation in the Macquarie River catchment,NSW, Australia

Widespread afforestation has been proposed as one means of addressing the increasing dryland and stream salinity problem in Australia. However, modelling results presented here suggest that large-scale tree planting will substantially reduce river flows and impose costs on downstream water users if planted in areas of high runoff yield. Streamflow reductions in the Macquarie River, NSW, Australia are estimated for a number of tree planting scenarios and global warming forecasts. The modelling framework includes the Sacramento rainfall-runoff model and IQQM, a streamflow routing tool, as well as various global climate model outputs from which daily rainfall and potential evaporation data files have been generated in OzClim, a climate scenario generator. For a 10% increase in tree cover in the headwaters of the Macquarie, we estimate a 17% reduction in inflows to Burrendong Dam. The drying trend for a mid-range scenario of regional rainfall and potential evaporation caused by a global warming of 0.5 degree C may cause an additional 5% reduction in 2030. These flow reductions will decrease the frequency of bird-breeding events in Macquarie Marshes (a RAMSAR protected wetland) and reduce the security of supply to irrigation areas downstream. Inter-decadal climate variability is predicted to have a very significant influence on catchment hydrologic behaviour. A further 20% reduction in flows from the long-term historical mean is possible, should we move into an extended period of below average rainfall years, such as occurred in eastern Australia between 1890 and 1948. Because current consumptive water use is largely adapted to the wetter conditions of post 1949, a return to prolonged dry periods would cause significant environmental stress given the agricultural and domestic water developments that have been instituted.     

Assessing Surface Water Quality and Its Relation with Urban Land Cover Changes in the Lake Calumet Area, Greater Chicago

Urban land use and land cover change significantly affect spatial and temporal patterns of runoff, which in turn impacts surface water quality. With the exponential growth in urban areas over the past three decades, changes in land use and land cover to cater for the growth of cities has been a conspicuous spectacle in urban spaces. The main goal of this study was to assess the impacts of land cover change on runoff and surface water quality using a partial area hydrology framework. The study employed ArcHydro GIS extension and a modified version of Long-Term Hydrologic and Nonpoint Source Pollution model (L-THIA-NPS) in estimating runoff and nonpoint source pollutant concentration around Lake Calumet between 1992 and 2001. Data employed include National Land Cover Data set, rainfall data, digital elevation model (DEM), Soil Survey Geographic (SSURGO) data, and The United States Environmental Protection Agency’s STORET (storage and retrieval) water quality data. The model was able to predict surface water quality reasonably well over the study period. Sensitivity analysis facilitated a manual calibration of the model. Model validation was executed by comparing simulated results following calibration and observed water quality data for the study area. The study demonstrates that the level of concentration of nonpoint source pollutants in surface water within an urban watershed heavily depends on the spatiotemporal variations in areas that contribute towards runoff compared to the spatial extent of change in major land use/land cover.     

Design and risk assessment tool for vegetative treatment areas receiving agricultural wastewater: Preliminary results

Vegetative treatment areas (VTAs) are commonly being used as an alternative method of agricultural process wastewater treatment. However, it is also apparent that to completely prevent discharge of pollutants to the surrounding environment, settling of particulates and bound constituents from overland flow through VTAs is not sufficient. For effective remediation of dissolved agricultural pollutants, VTAs must infiltrate incoming wastewater. A simple water balance model for predicting VTA soil saturation and surface discharge in landscapes characterized by sloping terrain and a shallow restrictive layer is presented and discussed. The model accounts for the cumulative effect of successive rainfall events and wastewater input on soil moisture status and depth to water table. Nash–Sutcliffe efficiencies ranged from 0.65 to 0.81 for modeled and observed water table elevations after calibration of saturated hydraulic conductivity. Precipitation data from relatively low, average, and high annual rainfall years were used with soil, site, and contributing area data from an example VTA for simulations and comparisons. Model sensitivity to VTA width and contributing area (i.e. barnyard, feedlot, silage bunker, etc.) curve number was also investigated. Results of this analysis indicate that VTAs should be located on steeper slopes with deeper, more-permeable soils, which effectively lowers the shallow water table. In sloping landscapes (>2%), this model provides practitioners an easy-to-use VTA design and/or risk assessment tool that is more hydrological process-based than current methods.     

MODELING THE EFFECT OF BRUSH CONTROL ON RANGELAND WATER YIELD1

ABSTRACT: With the increase in water demand in Texas, attention has turned to improving water yield by brush control on rangeland watersheds. Several hydrologic models have been developed for either farmland or rangeland. However, none of the models were specifically developed to assess the impact of brush control on rangeland water yield. Yet, modeling the impact of brush control on water yield needs to be considered if alternative techniques are to be compared. Two models, Ekalaka Rangeland Hydrology and Yield Model (ERHYM-II) and Simulator for Water Resources on Rural Basins (SWRRB) were selected. The Soil Conservation Service curve number (SCS-CN) method is used in both models to predict surface runoff from each rainfall event. The major differences between the ERHYM-II and SWRRB models are the evapotranspiration, soil water routing, and plant growth components. The models were evaluated on brush-dominated and chemically and mechanically brush-controlled range watersheds in Texas. Results indicated that both models were capable of simulating soil water and water yield from brush dominated and chemically brush-controlled range watersheds. The models were not able to predict water yield from the mechanically brush-controlled (root plowed) watershed with acceptable accuracy. The depressions that were caused by root plowing stored surface runoff and reduced water yield from the watershed. Information about the size of depressions was not available for further model evaluation.     

Modeling the Effect of Summertime Heating on Urban Runoff Temperature1

Abstract: The summertime heating of runoff in urban areas is recognized as a common and consistent urban climatological phenomenon. In this study, a simple thermal urban runoff model (TURM) is presented for the net energy flux at the impervious surfaces of urban areas to account for the heat transferred to runoff. The first step in developing TURM consists of calculating the various factors that control how urban impervious areas absorb heat and transfer it to moving water on the surface. The runoff temperature is determined based on the interactions of the physical characteristics of the impervious areas, the weather, and the heat transfer between the moving film of runoff and the impervious surface common in urban areas. Key surface and weather factors that affect runoff temperature predictions are type of impervious surface, air temperature, humidity, solar radiation before and during rain, rainfall intensity, and rainfall temperature. Runoff from pervious areas is considered separately and estimated using the Green‐Ampt Mein‐Larson rainfall excess method. Pervious runoff temperature is estimated as the rainfall temperature. Field measurements indicate that wet bulb temperature can be used as a surrogate for rainfall temperature and that runoff temperatures from sod average just 2°C higher than rainfall temperatures. Differences between measured and predicted impervious runoff temperature average approximately 2°C, indicating that TURM is a useful tool for determining runoff temperatures for typical urban areas.     

A MODELING APPROACH FOR OPTIMUM SEDIMENT DETENTION POND DESIGN1

ABSTRACT: A design procedure to determine optimum size for a sediment detention pond is presented. The procedure is based on simulating the sediment removal efficiency of the pond in conjunction with temporal variations in rainfall and potential land use and/or management options. The simulation procedure is based on a combined probabilistic-deterministic modeling approach. The probabilistic model generates daily rainfall with hourly increments for a selected site. The deterministic model simulates sediment yield and concentration for drainage area (pond inflow) and sediment trapping efficiency of the pond. The sediment yield and concentration in pond effluent is estimated from the difference between sediment inflow to the pond and sediment trapped by the pond. As an example, the procedure is applied to determine optimum design for a sediment detention pond in a surface mined area using several pond design options and alternative mining operation/land reclamation strategies.     

Use of Predictive Weather Uncertainties in an Irrigation Scheduling Tool Part II: An Application of Metrics and Adjoints

We apply predictive weather metrics and land model sensitivities to improve the Colorado State University Water Irrigation Scheduler for Efficient Application (WISE). WISE is an irrigation decision aid that integrates environmental and user information for optimizing water use. Rainfall forecasts and verification performance metrics are used to estimate predictive rainfall probabilities that are used as input data within the irrigation decision aid. These input data errors are also used within a land model sensitivity study to diagnose important prognostic water movement behaviors for irrigation tool development purposes simultaneously performing the analysis in space and time. Thus, important questions such as “how long can a crop water application be delayed while maintaining crop yield production?” are addressed by evaluating crop growth stage interactions as a function of soil depth (i.e., space), rainfall events (i.e., time), and their probabilistic uncertainties. Editor’s note : This paper is part of the featured series on Optimizing Ogallala Aquifer Water Use to Sustain Food Systems. See the February 2019 issue for the introduction and background to the series.     

PARAMETER UNCERTAINTY IN RAINFALL-RUNOFF MODELING FOR POLDER SYSTEM DESIGN1

ABSTRACT: An approach is developed for incorporating the uncertainty of parameters for estimating runoff in the design of polder systems in ungaged watersheds. Monte Carlo Simulation is used to derive a set of realizations of streamflow hydrographs for a given design rainstorm using the U. S. Soil Conservation Service (SCS) unit hydrograph model. The inverse of the SCS curve number, which is a function of the antecedent runoff condition in the SCS model, is the random input in the Monte Carlo Simulation. Monte Carlo realizations of streamfiow hydrographs are used to simulate the performance of a polder flood protection system. From this simulation the probability of occurrence of flood levels for a particular hydraulic design may be used to evaluate its effectiveness. This approach is demonstrated for the Pluit Polder flood protection system for the City of Jakarta, Indonesia. While the results of the application indicate that uncertainty in the antecedent runoff condition is important, the effects of uncertainty in rainfall data, in additional runoff parameters, such as time to peak, in the hydraulic design, and in the rainfall-runoff model selected should also be considered. Although, the SCS model is limited to agricultural conditions, the approach presented herein may be applied to other flood control systems if appropriate storm runoff models are selected.     

Rainfall conditions and rainwater harvesting potential in the urban area of Khartoum

Runoff water management is among the inherent challenges which face the sustainability of the development of arid urban centers. These areas are particularly at risk from flooding due to rainfall concentration in few heavy showers. On the other hand, they are susceptible to drought. The capital of Sudan (Khartoum) stands as exemplary for these issues. Hence, this research study aims at investigating the potential of applying rainwater harvesting (RWH) in Khartoum City Center as a potential urban runoff management tool. Rapid urbanization coupled with the extension of impervious surfaces has intensified the heat island in Khartoum. Consequently, increased frequency of heat waves and dust storms during the dry summer and streets flooding during the rainy season have led to environmental, economical, and health problems. The study starts with exposing the rainfall behavior in Khartoum by investigating rainfall variability, number of raindays, distribution of rain over the season, probability of daily rainfall, maximum daily rainfall and deficit/surplus of rain through time. The daily rainfall data show that very strong falls of >30 mm occur almost once every wet season. Decreased intra- and inter-annual rainfall surpluses as well as increased rainfall concentration in the month of August have been taking place. The 30-year rainfall variability is calculated at decade interval since 1941. Increasing variability is revealed with 1981–2010 having coefficients of variation of 66.6% for the annual values and 108.8–118.0% for the wettest months (July–September). Under the aforementioned rainfall conditions, this paper then explores the potential of RWH in Khartoum City Center as an option for storm water management since the drainage system covers only 40% of the study area. The potential runoff from the 6.5 km² center area is computed using the United States Natural Resources Conservation Services method (US-NRCS), where a weighted Curve Number (CN) of 94% is found, confirming dominant imperviousness. Rainfall threshold for runoff generation is found to be 3.3 mm. A 24,000 m³ runoff generated from a 13.1 mm rainfall (with 80% probability and one year return period) equals the drainage system capacity. An extreme rainfall of 30 mm produces a runoff equivalent to fourfold the drainage capacity. It is suggested that the former and latter volumes mentioned above could be harvested by applying the rational method from 18% and 80% rooftops of the commercial and business district area, respectively. Based on the above results, six potential sites can be chosen for RWH with a total roof catchment area of 39,558 m² and potential rooftop RWH per unit area of 0.033 m³. These results reflect the RWH potential for effective urban runoff management and better water resources utilization. RWH would provide an alternative source of water to tackle the drought phenomenon.     

UPLAND SOIL EROSION SIMULATION FOR AGRICULTURAL WATERSHEDS1

ABSTRACT: A soil erosion simulation model that considered the physical conditions of agricultural watersheds and that interfaced with the modified USDAHL-74 watershed hydrology model was developed. The erosion model simulates the detachment and transport of soil particles caused by raindrop impact and overland flow from rill and interrill areas. The model considers temporal and spatial variation of plant residue, crop canopy cover, snow cover, and the moisture content of surface soil as modifying factors of the erosive forces of raindrop impact and overland flow. The hydrology model simulates overland flow and some of the physical parameters that are used in the erosion model. The simulation is executed in the time interval determined by the rainfall rate or snowmelt rate. The erosion model compares the transport capacity of the overland flow and the sediment loaded in the overland flow to determine the fate account for the free soil particles that have already been detached and are readily available to be transported by the overland flow. The model was tested with data from two small agricultural watersheds in the Palouse region of the Pacific Northwest dryland. The model was calibrated by trial-and-error to determine the coefficients of the model.     

USE OF A HYDROLOGIC MODEL IN A BASIN-WIDE WATER ALLOCATION PROCEEDING1

ABSTRACT: The Montana Department of Natural Resources and Conservation developed a hydrologic model to help analyze the effects of allocating water for consumptive and instream uses in the upper Missouri River basin of Montana. The model, a PC-based FORTRAN program, uses a mass-balance approach to compute monthly streamflows, reservoir operations, hydropower production, and irrigation and municipal water uses throughout the 54,000 square mile basin for a 59-year base period. Simulation results are presented as monthly mean and percentile-exceedence values. The model was run for baseline conditions and six hypothetical water-allocation alternatives. Results were used by staff resource area specialists to assess potential impacts to water quantity and distribution, water rights, water quality, stream channel form, fisheries, wildlife, recreation, hydropower production, and economics. These analyses were presented to the public and the decision-making board in an environmental impact statement (EIS). Though, in many instances, the model did not allow for detailed, site-specific analyses, the model was an important tool and its simulation results formed the hydrologic basis for the EIS.     

QUANTIFYING MODEL OUTPUT UNCERTAINTY DUE TO SPATIAL VARIABILITY OF RAINFALL1

Traditionally in the application of hydrologic/water quality (H/WQ) models, rainfall is assumed to be spatially homogeneous and is considered not to contribute to output uncertainty. The objective of this study was to assess the uncertainty induced in model outputs solely due to rainfall spatial variability. The study was conducted using the AGNPS model and the rainfall pattern captured by a network of 17 rain gauges. For each rainfall event, the model was run using the rainfall captured by each rain gauge, one at a time, under the assumption of rainfall spatial homogeneity. A large uncertainty in the modeled outputs resulted from the rainfall spatial variability. The uncertainty in the modeled outputs exceeded the input rainfall uncertainty. Results of this study indicate that spatial variability of rainfall should be captured and used in H/WQ models in order to accurately assess the release and transport of pollutants. A large uncertainty in the model outputs can be expected if this rainfall property is not taken into account.     

GROUND WATER QUALITY IMPLICATIONS OF SOIL CONSERVATION MEASURES: AN ECONOMIC PERSPECTWE1

ABSTRACT: An evaluation of the intermedia movement of pesticides applied under various land management systems already in place, or to be implemented, under the Conservation Reserve and Conservation Compliance programs is presented. The simulation modeling approach followed in this analysis consists of a mathematical programming model and leaching/surface runoff, Pesticide Root Zone Model (PRZM) models. Special care was taken to ensure that the physical model was sensitive to the chemical characteristics of individual pesticides and the important physical changes brought about by different agricultural practices. Results show that, although these programs as now planned, increase farm income and achieve soil conservation goals, they may adversely affect ground water quality. Also, depending on soil and location characteristics, there are tradeoffs between surface and ground water quality implications. Hence, if these programs are to address water quality problems, the recommended practices must be evaluated for their impact on water quality, particularly in potentially vulnerable areas.     

An evaluation of the impact of persistent water level changes on the areal extent of Georgian Bay/North Channel marshlands

A conceptual mathematical model has recently been devised to assist environmental managers in predicting the impact on coastal marsh areas of long-term changes in water levels. The model considers such impact solely in terms of the geometry of the confining basin, the change in ambient water level, and the maximum depth for which bottom-rooted emergent vegetation is present. This model is applied to 17 shoreline marshes of various shapes in the Georgian Bay/North Channel region of the Great Lakes.Model outputs of predicted maximum and minimum marsh area subsequent to changes in long-term levels are compared to marsh areas measured from available historical air photos dating from 1935 to 1985. The results of such comparisons indicate that such a geometric model, despite its neglect of the biological complexities of marsh ecology, can serve as a valuable tool for assessing the range of impacts of both natural and man-made changes in long-term ambient water levels on shoreline marshes.     

AN INTEGRATED ONE‐DIMENSIONAL AND TWO‐DIMENSIONAL URBAN STORMWATER FLOOD SIMULATION MODEL1

ABSTRACT: Flash flooding is the rapid flooding of low lying areas caused by the stormwater of intense rainfall associated with thunderstorms. Flash flooding occurs in many urban areas with relatively flat terrain and can result in severe property damage as well as the loss of lives. In this paper, an integrated one‐dimensional (1‐D) and two‐dimensional (2‐D) hydraulic simulation model has been established to simulate stormwater flooding processes in urban areas. With rainfall input, the model simulates 2‐D overland flow and 1‐D flow in underground stormwater pipes and drainage channels. Drainage channels are treated as special flow paths and arranged along one or more sides of a 2‐D computational grid. By using irregular computation grids, the model simulates unsteady flooding and drying processes over urban areas with complex drainage systems. The model results can provide spatial flood risk information (e.g., water depth, inundation time and flow velocity during flooding). The model was applied to the City of Beaumont, Texas, and validated with the recorded rainfall and runoff data from Tropical Storm Allison with good agreement.     

EVALUATION OF HYDROLOGIC BENEFITS OF INFILTRATION BASED URBAN STORM WATER MANAGEMENT1

ABSTRACT: As watersheds are urbanized, their surfaces are made less pervious and more channelized, which reduces infiltration and speeds up the removal of excess runoff. Traditional storm water management seeks to remove runoff as quickly as possible, gathering excess runoff in detention basins for peak reduction where necessary. In contrast, more recently developed “low impact” alternatives manage rainfall where it falls, through a combination of enhancing infiltration properties of pervious areas and rerouting impervious runoff across pervious areas to allow an opportunity for infiltration. In this paper, we investigate the potential for reducing the hydrologic impacts of urbanization by using infiltration based, low impact storm water management. We describe a group of preliminary experiments using relatively simple engineering tools to compare three basic scenarios of development: an undeveloped landscape; a fully developed landscape using traditional, high impact storm water management; and a fully developed landscape using infiltration based, low impact design. Based on these experiments, it appears that by manipulating the layout of urbanized landscapes, it is possible to reduce impacts on hydrology relative to traditional, fully connected storm water systems. However, the amount of reduction in impact is sensitive to both rainfall event size and soil texture, with greatest reductions being possible for small, relatively frequent rainfall events and more pervious soil textures. Thus, low impact techniques appear to provide a valuable tool for reducing runoff for the events that see the greatest relative increases from urbanization: those generated by the small, relatively frequent rainfall events that are small enough to produce little or no runoff from pervious surfaces, but produce runoff from impervious areas. However, it is clear that there still needs to be measures in place for flood management for larger, more intense, and relatively rarer storm events, which are capable of producing significant runoff even for undeveloped basins.