Integrated Modular Modeling of Water and Nutrients From Point and Nonpoint Sources in the Patuxent River Watershed1

Abstract: We present a simple modular landscape simulation model that is based on a watershed modeling framework in which different sets of processes occurring in a watershed can be simulated separately with different models. The model consists of three loosely coupled submodels: a rainfall‐runoff model (TOPMODEL) for runoff generation in a subwatershed, a nutrient model for estimation of nutrients from nonpoint sources in a subwatershed, and a stream network model for integration of point and nonpoint sources in the routing process. The model performance was evaluated using monitoring data in the watershed of the Patuxent River, a tributary to the Chesapeake Bay in Maryland, from July 1997 through August 1999. Despite its simplicity, the landscape model predictions of streamflow, and sediment and nutrient loads were as good as or better than those of the Hydrological Simulation Program‐Fortran model, one of the most widely used comprehensive watershed models. The landscape model was applied to predict discharges of water, sediment, silicate, organic carbon, nitrate, ammonium, organic nitrogen, total nitrogen, organic phosphorus, phosphate, and total phosphorus from the Patuxent watershed to its estuary. The predicted annual water discharge to the estuary was very close to the measured annual total in terms of percent errors for both years of the study period (≤2%). The model predictions for loads of nutrients were also good (20‐30%) or very good (<20%) with exceptions of sediment (40%), phosphate (36%), and organic carbon (53%) for Year 1.     

Watershed Regressions for Pesticides (warp) Models for Predicting Atrazine Concentrations in Corn Belt Streams1

Stone, Wesley W. and Robert J. Gilliom, 2012. Watershed Regressions for Pesticides (WARP) Models for Predicting Atrazine Concentrations in Corn Belt Streams. Journal of the American Water Resources Association (JAWRA) 48(5): 970‐986. DOI: 10.1111/j.1752‐1688.2012.00661.x Abstract: Watershed Regressions for Pesticides (WARP) models, previously developed for atrazine at the national scale, are improved for application to the United States (U.S.) Corn Belt region by developing region‐specific models that include watershed characteristics that are influential in predicting atrazine concentration statistics within the Corn Belt. WARP models for the Corn Belt (WARP‐CB) were developed for annual maximum moving‐average (14‐, 21‐, 30‐, 60‐, and 90‐day durations) and annual 95th‐percentile atrazine concentrations in streams of the Corn Belt region. The WARP‐CB models accounted for 53 to 62% of the variability in the various concentration statistics among the model‐development sites. Model predictions were within a factor of 5 of the observed concentration statistic for over 90% of the model‐development sites. The WARP‐CB residuals and uncertainty are lower than those of the National WARP model for the same sites. Although atrazine‐use intensity is the most important explanatory variable in the National WARP models, it is not a significant variable in the WARP‐CB models. The WARP‐CB models provide improved predictions for Corn Belt streams draining watersheds with atrazine‐use intensities of 17 kg/km² of watershed area or greater.     

A Simple Process‐Based Snowmelt Routine to Model Spatially Distributed Snow Depth and Snowmelt in the SWAT Model1

Abstract: We present a method to integrate a process‐based (PB) snowmelt model that requires only daily temperature and elevation information into the Soil and Water Assessment Tool (SWAT) model. The model predicts the spatiotemporal snowpack distribution without adding additional complexity, and in fact reduces the number of calibrated parameters. To demonstrate the utility of the PB model, we calibrate the PB and temperature‐index (TI) SWAT models to optimize agreement with stream discharge on a 46‐km² watershed in northwestern Idaho, United States, for 10 individual years and use the calibrated parameters for the year with the best agreement to run the model for 15 remaining years. Stream discharge predictions by the PB and TI model were similar, although the PB model simulated snowmelt more accurately than the TI model for the remaining 15‐year period. Spatial snow distributions predicted by the PB model better matched observations from LandSat imagery and a SNOTEL station. Results for this watershed show that including PB snowmelt in watershed models is feasible, and calibration of TI‐based watershed models against discharge can incorrectly predict snow cover.     

POTENTIAL PHOSPHORUS LOAD REDUCTIONS UNDER THE LAKE OKEECHOBEE REGULATORY PROGRAM1

ABSTRACT: Nutrient loading from beef pastures located within the northern Lake Okeechobee watershed in Florida, has been identified as a source of phosphorus contributing to the accelerated eutrophication of the lake. Since 1989 within the watershed, 557 agricultural drainage sites, mainly beef pasture, have been monitored for compliance under a regulatory program. Of those sites, 154 were actively monitored for phosphorus concentrations from October 1, 1998, to September 30, 1999. Of these 154 sites, 77 were considered to be out of compliance (OOC). An OOC site is defined as having runoff with a 12‐month average phosphorus concentration exceeding the permitted discharge limit. The average annual phosphorous load from the 77 OOC sites for an eight‐year study period from October 1, 1991, to September 30, 1999, was estimated using measured concentration values and simulated runoff obtained from an agricultural nonpoint source pollution model, CREAMS‐WT. The 77 OOC sites produced an estimated average annual 46 metric tonnes of phosphorus load, of which an estimated 22 tonnes of phosphorus reached Lake Okeechobee on an average annual basis. The remaining estimated average annual 24 tonnes of phosphorus load was retained by streams and wetlands in the discharge transport system between the sites and the lake. The estimated average annual load reaching Lake Okeechobee from the OOC sites represented 11 percent of the phosphorus load above a five‐year average annual target load for the lake. However, the OOC site drainage areas represented only 3 percent of the northern watershed that drains into the lake. Of the 77 OOC sites, 12 sites had an average annual phosphorus loading rate equal to or greater than 3.0 kg/ha and were placed on the priority list for the Critical Restoration Project in the Lake Okeechobee watershed. To estimate the possible phosphorus load reductions from the 77 sites, two scenarios were modeled. The first scenario reduced phosphorus concentrations in runoff to the permitted discharge limits under the Lake Okeechobee regulatory program. The second scenario changed current land uses to native rangeland with an estimated annual offsite total phosphorus areal loading rate of 0.114 kg/ha. These two scenarios are hypothetical with assumed concentration values and loading rate. Model results showed that the first management scenario reduced the average annual phosphorus load to the lake by an estimated 15 tonnes. The second scenario reduced the average annual phosphorus load to the lake by an estimated 21 tonnes.     

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.     

IMPACTS OF CALIFORNIA'S CLIMATIC REGIMES AND COASTAL LAND USE CHANGE ON STREAMFLOW CHARACTERISTICS1

ABSTRACT: To investigate the impacts of urbanization and climatic fluctuations on stream flow magnitude and variability in a Mediterranean climate, the HEC‐HMS rainfall/runoff model is used to simulate stream flow for a 14‐year period (October 1, 1988, to September 30, 2002) in the Atascadero Creek watershed located along the southern coast of California for 1929, 1998, and 2050 (estimated) land use conditions (8, 38 and 52 percent urban, respectively). The 14‐year period experienced a range of climatic conditions caused mainly by El Nino‐Southern Oscillation variations. A geographic information system is used to delineate the watershed and parameterize the model, which is calibrated using data from two stream flow and eight rainfall gauges. Urbanization is shown to increase peak discharges and runoff volume while decreasing stream flow variability. In all cases, the annual and 14‐year distributions of stream flow are shown to be highly skewed, with the annual maximum 24 hours of discharge accounting for 22 to 52 percent of the annual runoff and the maximum ten days of discharge from an average El Nino year producing 10 to 15 percent of the total 14‐year discharge. For the entire period of urbanization (1929 to 2050), the average increase in annual maximum discharges and runoff was 45 m³/s (300 percent) and 15 cm (350 percent), respectively. Additionally, the projected increase in urbanization from 1998 to 2050 is half the increase from 1929 to 1998; however, increases in runoff (22 m³/s and 7 cm) are similar for both scenarios because of the region's spatial development pattern.     

Addressing data challenges in riverine nutrient load modeling of an intensively managed agro-industrial watershed

Data limitations often challenge the reliability of water quality models, especially in intensively managed watersheds. While numerous studies report successful hydrological model setup and calibration, few have addressed in detail the data challenges for multisite and multivariable model calibration to an intensively managed watershed. In this study, we address some of these challenges based on our reflective experience calibrating the Soil and Water Assessment Tool (SWAT) to the Upper Sangamon River Watershed in central Illinois based on daily flow, annual crop yield, and monthly sediment, nitrate, and total phosphorus loads. We highlight some challenges in SWAT calibration processes due to data errors and inconsistencies, and insufficient precipitation and water quality observations. Following, we demonstrate the merits of additional weather and water quality observations that could help reduce input uncertainties, and we provide suggestions for selecting appropriate observations for the model calibration. After dealing with the data issues, we show that the SWAT model could be calibrated with acceptable results for the case study watershed.     

EVALUATION OF RUNOFF,EROSION, AND PHOSPHORUS MODELING SYSTEM—SIMPLE1

ABSTRACT: The purpose of this study was to evaluate the performance of Spatially Integrated Models for Phosphorus Loading and Erosion (SIMPLE) in predicting runoff volume, sediment loss, and phosphorus loading from two watersheds. The modeling system was applied to the 334 ha QOD subwatershed, part of the Owl Run watershed, located in Fauquier County, Virginia, and to the 2240 ha watershed, Battle Branch, located in Delaware County, Oklahoma. Simulation runs were conducted at cell and field scales, and simulation results were compared with observed data. Runoff volume and dissolved phosphorus loading were measured at the Battle Branch watershed. Runoff volume, sediment yield, and total phosphorus loading were measured at the QOD site. SIMPLE tended to underestimate runoff volumes during the dormant period, from November to March. The comparison between observed and predicted dissolved phosphorus showed better correlation than for observed and predicted total phosphorus loading. Cell level simulations provided similar estimates of runoff volume and phosphorus loading when compared to field level simulations for both watersheds. However, observed sediment yields better compared with the values predicted from the cell level simulation when compared to field level simulation. Finally, results of model evaluation indicated that SIMPLE's predictive ability is acceptable for screening applications but not for site-specific quantitative predictions.     

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.     

NONPOINT SOURCE DISCHARGES OF NUTRIENTS FROM PIEDMONT WATERSHEDS OF CHESAPEAKE BAY1

ABSTRACT: We measured annual discharges of water, sediments, and nutrients from 10 watersheds with differing proportions of agricultural lands in the Piedmont physiographic province of the Chesapeake Bay drainage. Flow-weighted mean concentrations of total N, nitrate, and dissolved silicate in watershed discharges were correlated with the proportion of cropland in the watershed. In contrast, concentrations of P species did not correlate with cropland. Organic P and C correlated with the concentration of suspended particles, which differed among watersheds. Thus, the ratio of N:P:Si in discharges differed greatly among watersheds, potentially affecting N, P or Si limitation of phytoplankton growth in the receiving waters. Simple regression models of N discharge versus the percentage of cropland suggest that croplands discharge 29–42 kg N ha^-1 yr^-1 and other lands discharge 1.2–5.8 kg N ha^-1 yr^-1. We estimated net anthropogenic input of N to croplands and other lands using county level data on agriculture and N deposition from the atmosphere. For most of the study watersheds, N discharge amounted to less than half of the net anthropogenic N.     

SENSITIVITY ANALYSIS,CALIBRATION, AND VALIDATIONS FOR A MULTISITE AND MULTIVARIABLE SWAT MODEL1

The ability of a watershed model to mimic specified watershed processes is assessed through the calibration and validation process. The Soil and Water Assessment Tool (SWAT) watershed model was implemented in the Beaver Reservoir Watershed of Northwest Arkansas. The objectives were to: (1) provide detailed information on calibrating and applying a multisite and multivariable SWAT model; (2) conduct sensitivity analysis; and (3) perform calibration and validation at three different sites for flow, sediment, total phosphorus (TP), and nitrate‐nitrogen (NO₃‐N) plus nitrite‐nitrogen (NO₂‐N). Relative sensitivity analysis was conducted to identify parameters that most influenced predicted flow, sediment, and nutrient model outputs. A multi objective function was defined that consisted of optimizing three statistics: percent relative error (RE), Nash‐Sutcliffe Coefficient (R_NS²), and coefficient of determination (R²). This function was used to successfully calibrate and validate a SWAT model of Beaver Reservoir Watershed at multi‐sites while considering multivariables. Calibration and validation of the model is a key factor in reducing uncertainty and increasing user confidence in its predictive abilities, which makes the application of the model effective. Information on calibration and validation of multisite, multivariable SWAT models has been provided to assist watershed modelers in developing their models to achieve watershed management goals.     

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.     

NONPARAMETRIC STATISTICAL ANALYSIS OF ANNUAL PEAK FLOW DATA FROM A RECENTLY URBANIZED WATERSHED1

ABSTRACT. In order to demonstrate the feasibility of a nonparametric statistical application in investigating the hydrologic impact of the rapid land use change accompanying intense urbanization, annual maximum peak flow data from an actual example (the Northeast Branch Basin, a recently urbanized Washington, D.C., suburban watershed) were analyzed. Annual peak flow data from the Patuxent Basin above Unity, Maryland, a rural watershed in close proximity to the study area, were compared to data from the Northeast Branch for the Same period utilizing the Wilcoxon matched-pairs signed-ranks test. A change in central tendency of each series was noted at the 0.01 significance level; however, the change was negative in the rural basin and positive in the urbanized Northeast Branch Basin. This central tendency change was considered indicative of an average decrease in the size of rainstorms producing annual maximum peak discharges. Rainstorm data from the Northeast Branch Basin were divided into two equal periods (before urbanization and after) and the Wilcoxon test was applied. It was found that rainstorms producing maximum annual peak discharges in the urbanized period were indeed smaller than those in the prior period (0.01 level of significance); however, larger annual peaks were produced. It was concluded that nonparametric statistical methods can be used readily with conventional methods to isolate and clearly analyze the various problems in an actual urban hydrologic study.     

WATER BUDGET ANALYSIS FOR THE EVERGLADES AGRICULTURAL AREA DRAINAGE BASIN1

ABSTRACT: Water budget studies are essential for water resources and environmental management. In this study, a water budget analysis is presented for the Everglades Agricultural Area (EAA) in South Florida for the period from 1973 to 1991. The EAA is a highly productive irrigation/drainage basin that has a high water table and organic soils. Water quality problems are associated with the drainage discharge from the basin. During dry periods, supplemental water is used for irrigation and in rainy periods excess water with relatively higher phosphorus content is pumped out of the basin to Lake Okeechobee and the Everglades ecosystem. Elevated concentrations of phosphorus in the runoff/drainage that is discharged from the EAA basin have created water quality problems. The mean surface water inflow to the basin was 63,990 ha-m, and the outflow was 131,447 ha-m per year. On the average, supplemental surface water use was 47,411 ha-m, and runoff/drainage was 114,816 ha-m per year. The mean annual basin rainfall was 120.9 cm. A general trend in the decline of the wet season rainfall is observed.     

EVALUATING HSPF IN AN ARID,URBANIZED WATERSHED1

ABSTRACT: The Hydrologic Simulation Program‐FORTRAN (HSPF) is a powerful time variable hydrologic model that has rarely been applied in arid environments. Here, the performance of HSPF in southern California was assessed, testing its ability to predict annual volume, daily average flow, and hourly flow. The model was parameterized with eight land use categories and physical watershed characteristics. It was calibrated using rainfall and measured flow over a five‐year period in a predominantly undeveloped watershed and it was validated using a subsequent 4‐year period. The process was repeated in a separate, predominantly urbanized watershed over the same time span. Annual volume predictions correlated well with measured flow in both the undeveloped and developed watersheds. Daily flow predictions correlated well with measured flow following rain events, but predictions were poor during extended dry weather periods in the developed watershed. This modeling difficulty during dry‐weather periods reflects the large influence of, and the poor accounting in the model for, artificially introduced water from human activities, such as landscape overwatering, that can be important sources of water in urbanized arid environments. Hourly flow predictions mistimed peak flows, reflecting spatial and temporal heterogeneity of rainfall within the watershed. Model correlation increased considerably when predictions were averaged over longer time periods, reaching an asymptote after an 11‐hour averaging window.     

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.     

Use of SWAT to Estimate Spatial Scaling of Phosphorus Export Coefficients and Load Reductions Due to Agricultural BMPS

Phosphorus export coefficients (kg/ha/yr) from selected land covers, also called phosphorus yields, tend to get smaller as contributing areas get larger because some of the phosphorus mobilized on local fields gets trapped during transport to regional watershed outlets. Phosphorus traps include floodplains, wetlands, and lakes, which can then become impaired by eutrophication. The Sunrise River watershed in east central Minnesota, United States, has numerous lakes impaired by excess phosphorus. The Sunrise is tributary to the St. Croix River, whose much larger watershed is terminated by Lake St. Croix, also impaired by excess phosphorus. To support management of these impairments at both local and regional scales, a Soil and Water Assessment Tool (SWAT) model of the Sunrise watershed was constructed to estimate load reductions due to selected best management practices (BMPs) and to determine how phosphorus export coefficients scaled with contributing area. In this study, agricultural BMPs, including vegetated filter strips, grassed waterways, and reduction of soil‐phosphorus concentrations reduced phosphorus loads by 4‐20%, with similar percentage reductions at field and watershed spatial scales. Phosphorus export coefficients from cropland in rotation with corn, soybeans, and alfalfa decreased as a negative power function of contributing area, from an average of 2.12 kg/ha/yr at the upland field scale (~0.6 km²) to 0.63 kg/ha/yr at the major river basin scale (20,000 km²). 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.     

Evaluation of Three Watershed‐Scale Pesticide Environmental Transport and Fate Models1

Abstract: The U.S. Environmental Protection Agency (USEPA) Office of Pesticide Programs (OPP) has completed an evaluation of three watershed‐scale simulation models for potential use in Food Quality Protection Act pesticide drinking water exposure assessments. The evaluation may also guide OPP in identifying computer simulation tools that can be used in performing aquatic ecological exposure assessments. Models selected for evaluation were the Soil Water Assessment Tool (SWAT), the Nonpoint Source Model (NPSM), a modified version of the Hydrologic Simulation Program‐Fortran (HSPF), and the Pesticide Root Zone Model‐Riverine Water Quality (PRZM‐RIVWQ) model. Simulated concentrations of the pesticides atrazine, metolachlor, and trifluralin in surface water were compared with field data monitored in the Sugar Creek watershed of Indiana’s White River basin by the National Water Quality Assessment (NAWQA) program. The evaluation not only provided USEPA with experience in using watershed models for estimating pesticide concentration in flowing water but also led to the development of improved statistical techniques for assessing model accuracy. Further, it demonstrated the difficulty of representing spatially and temporally variable soil, weather, and pesticide applications with relatively infrequent, spatially fixed, point estimates. It also demonstrated the value of using monitoring and modeling as mutually supporting tools and pointed to the need to design monitoring programs that support modeling.     

Calibration and Verification of SWMM for Low Impact Development

The Storm Water Management Model was used to simulate runoff and nutrient export from a low impact development (LID) watershed and a watershed using traditional runoff controls. Predictions were compared to observed values. Uncalibrated simulations underpredicted weekly runoff volume and average peak flow rates from the multiple subcatchment LID watershed by over 80%; the single subcatchment traditional watershed had better predictions. Saturated hydraulic conductivity, Manning's n for swales, and initial soil moisture deficit were sensitive parameters. After calibration, prediction of total weekly runoff volume for the LID and traditional watersheds improved to within 12 and 5% of observed values, respectively. For the validation period, predicted total weekly runoff volumes for the LID and traditional watersheds were within 6 and 2% of observed values, respectively. Water quality simulation was less successful, Nash–Sutcliffe coefficients >0.5 for both calibration and validation periods were only achieved for prediction of total nitrogen export from the LID watershed. Simulation of a 100‐year, 24‐h storm resulted in a runoff coefficient of 0.46 for the LID watershed and 0.59 for the traditional watershed. Results suggest either calibration is needed to improve predictions for LID watersheds or expanded look‐up tables for Green–Ampt infiltration parameter values that account for compaction of urban soil and antecedent conditions are needed.     

Nutrient and Sediment Production,Watershed Characteristics,and Land Use in the Tahoe Basin,California‐Nevada1

Abstract: In efforts to control the degradation of water quality in Lake Tahoe, public agencies have monitored surface water discharge and concentrations of nitrogen, phosphorus, and suspended sediment in two separate sampling programs. The first program focuses on 20 watersheds varying in size from 162 to 14,000 ha, with continuous stream gaging and periodic sampling; the second focuses on small urbanized catchments, with automated sampling during runoff events. Using data from both programs, we addressed the questions (1) what are the fluxes and concentrations of nitrogen and phosphorus entering the lake from surface runoff; (2) how do the fluxes and concentrations vary in space and time; and (3) how are they related to land use and watershed characteristics? To answer these questions, we calculated discharge‐weighted average concentrations and annual fluxes and used multiple regression to relate those variable to a suite of GIS‐derived explanatory variables. The final selected regression models explain 47‐62% of the variance in constituent concentrations in the stormwater monitoring catchments, and 45‐72% of the variance in mean annual yields in the larger watersheds. The results emphasize the importance of impervious surface and residential density as factors in water quality degradation, and well‐developed soil as a factor in water quality maintenance.