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.     

GIS BASED LONG TERM HYDROLOGIC IMPACT EVALUATION FOR WATERSHED URBANIZATION1

ABSTRACT: To adequately manage impacts of ongoing or future land use changes in a watershed, the magnitude of their hydrologic impacts needs to be assessed. A grid based daily streamflow model was calibrated with two years of observed streamflow data, using time periods when land use data are available and verified by comparison of model predictions with observed streamflow data. Streamflow data were separated into direct runoff and baseflow to estimate the impacts of urbanization on each hydrologic component. Analysis of the ratio between direct runoff and total runoff from 30 years of simulation results and the change in these ratios with urbanization shows that estimated annual direct runoff increased from 49.2 percent (1973) to 63.1 percent (1984) and 65.0 percent (1991), indicating the effects of urbanization are greater on direct runoff than on total runoff. The direct runoff ratio also varies with annual rainfall, with dry year ratios larger than those for wet years. This suggests that the impact of urbanization on areas that are sensitive to runoff ratios, such as stream ecosystems, might be more serious during drier years than in wetter years in terms of water quality and water yield. This indicates that sustainable base‐flow is important to maintaining sound stream ecosystems.     

Evaluation of the MIKE SHE Model for Application in the Loess Plateau,China1

Abstract: Quantifying the hydrologic responses to land use/land cover change and climate variability is essential for integrated sustainable watershed management in water limited regions such as the Loess Plateau in Northwestern China where an adaptive watershed management approach is being implemented. Traditional empirical modeling approach to quantifying the accumulated hydrologic effects of watershed management is limited due to its complex nature of soil and water conservation practices (e.g., biological, structural, and agricultural measures) in the region. Therefore, the objective of this study was to evaluate the ability of the distributed hydrologic model, MIKE SHE to simulate basin runoff. Streamflow data measured from an overland flow‐dominant watershed (12 km²) in northwestern China were used for model evaluation. Model calibration and validation suggested that the model could capture the dominant runoff process of the small watershed. We found that the physically based model required calibration at appropriate scales and estimated model parameters were influenced by both temporal and spatial scales of input data. We concluded that the model was useful for understanding the rainfall‐runoff mechanisms. However, more measured data with higher temporal resolution are needed to further test the model for regional applications.     

HYDROLOGIC CALIBRATION OF THE SWAT MODEL IN A WATERSHED CONTAINING FRAGIPAN SOILS1

ABSTRACT: The Soil and Water Assessment Tool (SWAT) model, designed for use on rural ungaged basins and incorporating a GRASS GIS interface, was used to model the hydrologic response of the Ariel Creek watershed of northeastern Pennsylvania. Model evaluation of daily flow prior to calibration revealed a deviation of runoff volumes (D_v) of 68.3 percent and a Nash-Sutcliffe coefficient of-0.03. Model performance was affected by unusually large observed snowmelt events and the inability of the model to accurately simulate baseflow, which was influenced by the presence of fragipans. Seventy-five percent of the soils in the watershed contain fragipans. Model calibration yielded a D_v of 39.9 percent and a Nash-Sutcliffe coefficient of 0.04, when compared on a daily basis. Monthly comparisons yielded a Nash-Sutcliffe coefficient of 0.14. Snowmelt events in the springs of 1993 and 1994, which were unusually severe, were not adequately simulated. Neglecting these severe events, which produced the largest and third largest measured flows for the period of record, a D_v of 4.1 percent and Nash-Sutcliffe coefficient of 0.20 were calculated on a daily comparison, while on a monthly basis the Nash-Sutciffe coefficient was 0.55. These results suggest that the SWAT model is better suited to longer period simulations of hydrologic yields. Baseflow volumes were accurately simulated after calibration (D_v= -0.2 percent). Refinements made to the algorithms controlling subsurface hydrology and snowmelt, to better represent the presence of fragipans and snowmelt events, would likely improve model performance.     

EFFECTS OF INITIAL ABSTRACTION AND URBANIZATION ON ESTIMATED RUNOFF USING CN TECHNOLOGY1

ABSTRACT: Few studies have been conducted to explore the effects of initial abstraction on estimated direct runoff despite the widespread use of the curve number (CN) method in many hydrologic models to estimate direct runoff. In this study, use of a 5 percent ratio of initial abstraction (I_a) to storage (S) to estimate daily direct runoff with modified CN values for a 5 percent I_a/S value was investigated using the Long‐Term Hydrologic Impact Assessment (L‐THIA) geographic information system (GIS). In addition, the effects on estimated runoff of altering the hydrologic soil group due to urbanization were investigated. The L‐THIA model was applied to the Indiana Little Eagle Creek watershed with 5 percent and 20 percent I_a/S values, considering hydrologic soil group alteration due to urbanization. The results indicate that uses of a 5 percent l_a/S and modified CN values and Hydrologic Soil Group D for urbanized areas in model runs can improve long term direct runoff prediction.     

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.     

Modeling Hydrology in a Small Rocky Mountain Watershed Serving Large Urban Populations1

Abstract: This article describes the development of a calibrated hydrologic model for the Blue River watershed (867 km²) in Summit County, Colorado. This watershed provides drinking water to over a third of Colorado’s population. However, more research on model calibration and development for small mountain watersheds is needed. This work required integration of subsurface and surface hydrology using GIS data, and included aspects unique to mountain watersheds such as snow hydrology, high ground‐water gradients, and large differences in climate between the headwaters and outlet. Given the importance of this particular watershed as a major urban drinking‐water source, the rapid development occurring in small mountain watersheds, and the importance of Rocky Mountain water in the arid and semiarid West, it is useful to describe calibrated watershed modeling efforts in this watershed. The model used was Soil and Water Assessment Tool (SWAT). An accurate model of the hydrologic cycle required incorporation of mountain hydrology‐specific processes. Snowmelt and snow formation parameters, as well as several ground‐water parameters, were the most important calibration factors. Comparison of simulated and observed streamflow hydrographs at two U.S. Geological Survey gaging stations resulted in good fits to average monthly values (0.71 Nash‐Sutcliffe coefficient). With this capability, future assessments of point‐source and nonpoint‐source pollutant transport are possible.     

CALIBRATION OF STORM LOADS IN THE SOUTH PRONG WATERSHED,FLORIDA, USING BASINS/HSPF1

ABSTRACT: The South Prong watershed is a major tributary system of the Sebastian River and adjacent Indian River Lagoon. Continued urbanization of the Sebastian River drainage basin and other watersheds of the Indian River Lagoon is expected to increase runoff and nonpoint source pollutant loads. The St. Johns River Water Management District developed watershed simulation models to estimate potential impacts on the ecological systems of receiving waters and to assist planners in devising strategies to prevent further degradation of water resources. In the South Prong system, a storm water sampling program was carried out to calibrate the water quality components of the watershed model for total suspended solids (TSS), total phosphorous (TP), and total nitrogen (TN). During the period of May to November 1999, water quality and flow data were collected at three locations within the watershed. Two of the sampling stations were located at the downstream end of major watercourses. The third station was located at the watershed outlet. Five storm events were sampled and measured at each station. Sampling was conducted at appropriate intervals to represent the rising limb, peak, and recession limb of each storm event. The simulations were handled by HSPF (Hydrologic Simulation Program‐Fortran). Results include calibration of the hydrology and calibration of the individual storm loads. The hydrologic calibration was continuous over the period 1994 through 1999. Simulated storm runoff, storm loads, and event mean concentrations were compared with their corresponding observed values. The hydrologic calibration showed good results. The outcome of the individual storm calibrations was mixed. Overall, however, the simulated storm loads agreed reasonably well with measured loads for a majority of the storms.     

Incorporating Variable Source Area Hydrology into a Spatially Distributed Direct Runoff Model1

Buchanan, Brian, Zachary M. Easton, Rebecca Schneider, and M. Todd Walter, 2011. Incorporating Variable Source Area Hydrology Into a Spatially Distributed Direct Runoff Model. Journal of the American Water Resources Association (JAWRA) 48(1): 43‐60. DOI: 10.1111/j.1752‐1688.2011.00594.x Abstract: Few hydrologic models simulate both variable source area (VSA) hydrology, and runoff‐routing at high enough spatial resolutions to capture fine‐scale hydrologic pathways connecting VSA to the stream network. This paper describes a geographic information system‐based operational model that simulates the spatio‐temporal dynamics of VSA runoff generation and distributed runoff‐routing, including through complex artificial drainage networks. The model combines the Natural Resource Conservation Service’s Curve Number (CN) equation for estimating storm runoff with the topographic index concept for predicting the locations of VSA and a runoff‐routing algorithm into a new spatially distributed direct hydrograph (SDDH) model (SDDH‐VSA). Using a small agricultural watershed in central New York, SDDH‐VSA results were compared to those from a SDDH model using the traditional land use assumptions for the CN (SDDH‐CN). The SDDH‐VSA model generally agreed better with observed discharge than the SDDH‐CN model (average, Nash‐Sutcliffe efficiency of 0.69 vs. 0.58, respectively) and resulted in more realistic spatial patterns of runoff‐generating areas. The SDDH approach did not correctly capture the timing of runoff from small storms in dry periods. Despite this type of limitation, SDDH‐VSA extends the applicability of the SDDH technique to VSA conditions, providing a basis for new tools to help identify critical management areas and assess water quality risks due to landscape alterations.     

Forecasting land use change and its environmental impact at a watershed scale

Urban expansion is a major driving force altering local and regional hydrology and increasing non-point source (NPS) pollution. To explore these environmental consequences of urbanization, land use change was forecast, and long-term runoff and NPS pollution were assessed in the Muskegon River watershed, located on the eastern coast of Lake Michigan. A land use change model, LTM, and a web-based environmental impact model, L-THIA, were used in this study. The outcomes indicated the watershed would likely be subjected to impacts from urbanization on runoff and some types of NPS pollution. Urbanization will slightly or considerably increase runoff volume, depending on the development rate, slightly increase nutrient losses in runoff, but significantly increase losses of oil and grease and certain heavy metals in runoff. The spatial variation of urbanization and its impact were also evaluated at the subwatershed scale and showed subwatersheds along the coast of the lake and close to cities would have runoff and nitrogen impact. The results of this study have significant implications for urban planning and decision making in an effort to protect and remediate water and habitat quality of Muskegon Lake, which is one of Lake Michigan's Areas of Concern (AOC), and the techniques described here can be used in other areas.     

SENSITIVITY CONSIDERATIONS WHEN MODELING HYDROLOGIC PROCESSES WITH DIGITAL ELEVATION MODEL1

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

Predicting soil erosion for alternative land uses

The APEX (Agricultural Policy-Environmental eXtender) model developed in the United States was calibrated for northwestern China's conditions. The model was then used to investigate soil erosion effects associated with alternative land uses at the ZFG (Zi-Fang-Gully) watershed in northwestern China. The results indicated that the APEX model could be calibrated reasonably well (+/-15% errors) to fit those areas with >50% slope within the watershed. Factors being considered during calibration include runoff, RUSLE (Revised Universal Soil Loss Equation) slope length and steepness factor, channel capacity flow rate, floodplain saturated hydraulic conductivity, and RUSLE C factor coefficient. No changes were made in the APEX computer code. Predictions suggest that reforestation is the best practice among the eight alternative land uses (the status quo, all grass, all grain, all grazing, all forest, half tree and half grass, 70% tree and 30% grain, and construction of a reservoir) for control of water runoff and soil erosion. Construction of a reservoir is the most effective strategy for controlling sediment yield although it does nothing to control upland erosion. For every 1 Mg of crop yield, 11 Mg of soil were lost during the 30-yr simulation period, suggesting that expanding land use for food production should not be encouraged on the ZFG watershed. Grass species are less effective than trees in controlling runoff and erosion on steep slopes because trees generally have deeper and more stable root systems.     

增江流域非点源污染模拟研究

构建增江流域非点源污染数据库,包括DEM、土地利用,土壤类型,气象数据等,应用分布式流域水文模型SWAT(Soil andWater Assessment Tool,swat 2009版)对增江流域的非点源污染进行模拟。模型运行阶段为2000-2003年,分别应用2000-2001年和2002-2003年的实测月均流量及硝酸盐氮监测数据对模型的参数率定和验证,采用决定系数R2和Nash-Suttcliffe系数对模拟结果进行评定。其中水文模拟的R2均>0.9,Nash-Suttcliffe模型效率系数均>0.8;硝酸盐氮模拟的R2均>0.7,Nash-Suttcliffe模型效率系数均>0.6,表明SWAT模型在增江流域具有较好的适用性。     

ALAWAT: A spatially allocated watershed model for approximating stream,sediment, and pollutant flows in Hawaii,USA

The Ala Wai Canal Watershed Model (ALAWAT) is a planning-level watershed model for approximating direct runoff, streamflow, sediment loads, and loads for up to five pollutants. ALAWAT uses raster GIS data layers including land use, SCS soil hydrologic groups, annual rainfall, and subwatershed delineations as direct model parameter inputs and can use daily total rainfall from up to ten rain gauges and streamflow from up to ten stream gauges. ALAWAT uses a daily time step and can simulate flows for up to ten-year periods and for up to 50 subwatersheds. Pollutant loads are approximated using a user-defined combination of rating curve relationships, mean event concentrations, and loading/washoff parameters for specific subwatersheds, land uses, and times of year. Using ALAWAT, annual average streamflow and baseflow relationships and urban suspended sediment loads were approximated for the Ala Wai Canal watershed (about 10,400 acres) on the island of Oahu, Hawaii. Annual average urban suspended sediments were approximated using two methods: mean event concentrations and pollutant loading and washoff. Parameters for the pollutant loading and washoff method were then modified to simulate the effect of various street sweeping intervals on sediment loads.     

EXPORT COEFFICIENT MODELING TO ASSESS PHOSPHORUS LOADING IN AN URBAN WATERSHED1

ABSTRACT: An export coefficient modeling approach was used to assess the influence of land use on phosphorus loading to a Southern Ontario stream. A model was constructed for the 1995–1996 water year and calibrated within ± 3 percent of the observed mean concentration of total phosphorus. It was found that runoff from urban areas contributed most to the loading of phosphorus to the stream. When the model was assessed by running it for the 1977–1978 water year, using water quality and land use data collected independently, agreement within ± 7 percent was obtained. The model was then used to forecast the impact of future urban development proposed for the watershed, in terms of phosphorus loading, and to evaluate the reduction in loading resulting from several urban best management practices (BMP). It was determined that phosphorus removal will have to be applied to all the urban runoff from the watershed to appreciably reduce stream phosphorus concentration. Of the BMP designs assessed, an infiltration pond system resulted in the greatest phosphorus load reduction, 50 percent from the 1995–1996 baseline.     

An Optimization Method for Estimating Constituent Mean Concentrations in Base Flow‐Dominated Flow

Pollutant coefficients have been widely used to assess runoff nonpoint source pollution from individual land uses (e.g., agricultural, residential) of a watershed. Pollutant coefficients, known as event mean concentrations (EMCs), were developed by the U.S. Environmental Protection Agency's Nationwide Urban Runoff Program (NURP) to serve as a national measure for characterizing pollutant loading in a receiving water body. The term “baseflow pollutant coefficient (BPC)” is used in this study as a surrogate for EMC to describe mean concentration of pollutants in base flow‐dominated flow. A method for characterizing base flow quantity and quality for different land uses was explored using inverse modeling with two optimization techniques (a least square method and a genetic algorithm [GA] optimization), land use information, and streamflow quantity and quality data. The inverse model was formulated as a constrained minimization problem and demonstrated with data for 15 watersheds in Indiana. Results showed that estimated pollutant coefficients are comparable to the published literature. This indicates that the proposed method has the potential to effectively estimate constituent mean concentrations for pollutant load determination in gauged and ungauged watersheds, albeit more analysis with larger and more robust datasets is desirable to further refine and validate the accuracy of the approach.     

MODELING FRAMEWORK FOR MANAGING COPPER RUNOFF IN URBAN WATERSHEDS1

ABSTRACT: A modeling framework was developed for managing copper runoff in urban watersheds that incorporates water quality characterization, watershed land use areas, hydrologic data, a statistical simulator, a biotic ligand binding model to characterize acute toxicity, and a statistical method for setting a watershed specific copper loading. The modeling framework is driven by export coefficients derived from water quality parameters and hydrologic inputs measured in an urban watershed's storm water system. This framework was applied to a watershed containing a copper roof built in 1992. A series of simulations was run to predict the change in receiving stream water chemistry caused by roof aging and to determine the maximum copper loading (at the 99 percent confidence level) a watershed could accept without causing acute toxicity in the receiving stream. Forecasting the amount of copper flux responsible for exceeding the assimilation capacity of a watershed can be directly related to maximum copper loadings responsible for causing toxicity in the receiving streams. The framework developed in this study can be used to evaluate copper utilization in urban watersheds.     

WATERSHED WEIGHTING OF EXPORT COEFFICIENTS TO MAP CRITICAL PHOSPHOROUS LOADING AREAS1

ABSTRACT: The Export Coefficient model (ECM) is capable of generating reasonable estimates of annual phosphorous loading simply from a watershed's land cover data and export coefficient values (ECVs). In its current form, the ECM assumes that ECVs are homogeneous within each land cover type, yet basic nutrient runoff and hydrological theory suggests that runoff rates have spatial patterns controlled by loading and filtering along the flow paths from the upslope contributing area and downslope dispersal area. Using a geographic information system (GIS) raster, or pixel, modeling format, these contributing area and dispersal area (CADA) controls were derived from the perspective of each individual watershed pixel to weight the otherwise homogeneous ECVs for phosphorous. Although the CADA‐ECM predicts export coefficient spatial variation for a single land use type, the lumped basin load is unaffected by weighting. After CADA weighting, a map of the new ECVs addressed the three fundamental criteria for targeting critical pollutant loading areas: (1) the presence of the pollutant, (2) the likelihood for runoff to carry the pollutant offsite, and (3) the likelihood that buffers will trap nutrients prior to their runoff into the receiving water body. These spatially distributed maps of the most important pollutant management areas were used within New York's West Branch Delaware River watershed to demonstrate how the CADA‐ECM could be applied in targeting phosphorous critical loading areas.     

INTEGRATING LANDSCAPE ASSESSMENT AND HYDROLOGIC MODELING FOR LAND COVER CHANGE ANALYSIS1

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

Long-Term Impacts of Land-Use Change on Non-Point Source Pollutant Loads for the St. Louis Metropolitan Area, USA

A land-use-change simulation model (LEAM) and a non-point-source (NPS) water quality model (L-THIA) were closely coupled as LEAMwq in order to determine the long-term implications of various degree of urbanization on NPS total nitrogen (TN), total suspended particles (TSP), and total phosphorus (TP) loads. A future land-use projection in the St. Louis metropolitan area from 2005 to 2030 using three economic growth scenarios (base, low, and high) and a long-term precipitation dataset were used to predict the mean annual surface runoff and mean annual NPS pollutant loads in the region. Results show mean annual TN increases of 0.21%, 0.13%, and 0.14% by 2030 compared to 2000 under the base, high, and low scenarios, respectively. TSP and TP showed similar trends with different magnitudes. Corresponding changes in annual mean surface runoff were shown to be lower than expected, which might be attributed to the small-scale conversion pattern of land uses. In the most dramatic change (high growth) scenario, the runoff would increase across time but at varying rates, and temporal pollutant loads would result in a more complicated pattern than in the other scenarios. This is attributed to the complex interactions between event mean concentrations of pollutants and the magnitude of changes in land-use acreages. By integrating L-THIA with LEAM, LEAMwq was found to be a useful planning tool to illustrate in a quick and simple manner how future water quality is connected to decision-making on future land-use change.