Modeling Phosphorus Transport in the Blue River Watershed, Summit County, Colorado
Total Page:16
File Type:pdf, Size:1020Kb
Modeling Phosphorus Transport in the Blue River Watershed, Summit County, Colorado Paula Jo Lemonds and John E. McCray Abstract Introduction Lake Dillon in Summit County, Colorado, is a primary Numerical models are useful tools because they allow a drinking-water reservoir for Denver. Eutrophication of quantitative assessment of the environmental impacts Lake Dillon is a concern, primarily due to phosphorus of wastewater pollutants and improve understanding of (P) loading. There is little agriculture in the watershed. watershed-scale pollutant transport. Projecting future Thus, many local officials attribute the P loading to water quantity and quality is especially important in onsite wastewater systems (OWS). A watershed developing communities that rely on shallow modeling effort using the SWAT model is underway to groundwater as a source of drinking water while understand the potential influence of various point and disposing of wastewater in the shallow subsurface. nonpoint sources of P in the Blue River watershed (the Some models capable of simulating watershed-scale most developed of three watersheds that supply Lake pollutant transport include Soil and Water Assessment Dillon). The watershed model was calibrated to Tool (SWAT) (Arnold 1998), MIKESHE (Danish measured flow rates and P concentrations. The Hydraulic Institute 1999), Watershed Analysis Risk hydrologic model results are most sensitive to the Management Framework (WARMF) (Chen et al. 1999) physical parameters of snowmelt, and orographic and Hydrologic Simulation Program – Fortran (HSPF) effects on precipitation and evapotranspiration. (Bicknell et al. 1996). SWAT was used for this effort However, uncertainties in chemical-hydrologic because it is a public-domain model that can parameters preclude a rigorous assignment of relative incorporate large amounts of data and simulate many contributions of various P sources. Rather, the effort hydrologic processes. has resulted in a better understanding of P chemical parameters required to simulate watershed-scale Several watershed-scale models have been developed transport. The model was most sensitive to the P using SWAT (Arnold et al. 1999, Arnold and Allen sorption coefficient, the P availability index, and the P 1996, Fontaine et al. 2002, Manguerra and Engel 1998, enrichment ratio (a measure of P in runoff sediments Santhi et al. 2001, Srinivasan et al. 1998). However, compared to immobile sediments). Modeling results these projects did not specifically address the indicate that OWS are not significant sources of P to watershed-scale impacts of wastewater pollutants from Lake Dillon. onsite wastewater systems (OWS). The goals of this study are to accurately simulate mountain watershed Keywords: watershed, modeling, phosphorus, hydrology and to quantify the impacts of OWS-derived wastewater, hydrology, SWAT phosphorus (P) in the Lake Dillon Watershed. The study area is the Lake Dillon watershed located in Summit County, Colorado (Figure 1). Lake Dillon is the main drinking-water storage reservoir for Denver. Towns in the watershed include Frisco, Dillon, Lemonds is a research associate, Colorado School of Silverthorne, Breckenridge, and Blue River. Mines, Golden, CO 80401 (currently a water resources engineer at HDR, Denver, CO 80203). E-mail: Model Setup [email protected]. McCray is an associate professor at Colorado School of Mines, Department of The ArcView Interface for SWAT was used in model Geology and Geological Engineering, Golden, CO development. Subwatersheds were delineated using 80401. SWAT and a USGS 300-m resolution, 1-degree Digital 272 Elevation Model (DEM). Information extracted and watershed. Daily precipitation and minimum/maximum calculated from the DEM includes overland slope, temperature values were incorporated into the model. slope length, and elevation corrections for precipitation and evapotranspiration. The subwatershed delineation Other information defined in initial model setup is illustrated in Figure 1. included wastewater treatment plant point-source discharges into the Blue River and its tributaries, Lake Lake Dillon watershed Dillon water levels, reservoir outflow, and surface area of the reservoir. Information specific to each r Colorado 6 8 subwatershed, including stream-water chemical 1 3 2 13 10 9 34 4 r 7 properties, groundwater-flow properties, stream-routing 5 16 11 36 35 1433 12 parameters, consumptive water use, and agricultural 30 32 17 15 diversions were included. For details on the model 19 22 18 N 23 20 25 setup and simulation parameters, the reader is referred r 26 r 31 24 to Lemonds (2003). 27 21 29 28 Incorporation of OWS r r Weather Stations Currently, no algorithms exist in SWAT to specifically r Streams 5 0 5 10 Kilometers Subbasins simulate OWS. Therefore, a fertilizer management practice was used to simulate OWS input. The mass Figure 1. The study area is the Lake Dillon watershed input rate of OWS pollutants was set equal to the mass in Summit County, Colorado. of nutrient input by the fertilizer. OWS inputs were established from reviews of OWS effluent flow rates The land-use information was derived from 1:250,000- and water-quality parameters completed by Kirkland scale Landuse/Landcover Geographic Information (2001). Retrieval Analysis System (GIRAS) spatial data. The land use/land cover digital data were collected by the Fertilizer application input parameters were adjusted in USGS and converted to ARC/INFO by the USEPA. SWAT to achieve the appropriate inorganic P mass This information was used when simulating infiltration, input rate to the subsurface based on the number of runoff, ET, and natural sources of nutrients. OWS in each subwatershed. It was assumed that implementing the fertilizer management practice in The soil attributes were taken from the State Soil seven-day intervals would adequately represent OWS Geographic (STATSGO) Database, which was effluent processes. SWAT allows the user to apply the developed by the National Cooperative Soil Survey fertilizer into the first soil layer. As a result, the (USDA-NRCS Soil Survey 2002). In the STATSGO simulated OWS nutrients are not affected by runoff and database, soil attributes are stored in polygon format. are allowed to percolate through the vadose zone. The Each polygon includes multiple soil series with source of percolation is the natural precipitation, which information on its areal percentage of the polygon. In is orders of magnitude greater than OWS effluent input. the SWAT ArcView Interface, the dominant soil series is selected, and the interface extracts properties for the Results model from a relational database. Examples of the properties extracted include soil texture, bulk density, Prior to simulating nutrient transport, physical hydraulic conductivity, available water capacity, hydrologic input parameters were adjusted to calibrate organic carbon, and total depth of soil. These the model to stream flow rates. Adequate calibration of parameters are used in computations for infiltration, the physical hydrologic system was critical to runoff, groundwater flow, and P transport. simulating nutrient transport. Precipitation and temperature data were available from Hydrology Simulation and Calibration the National Climatic Data Center (NCDC) and the Natural Resources Conservation Service’s (NRCS) Measured streamflow was obtained from two USGS Snowpack Telemetry (SNOTEL) Data Network. Six gaging stations on the Blue River: one near the stations were available within the Lake Dillon headwaters and the other located approximately one 273 half mile upstream of Lake Dillon. The model precipitation are adjusted for each elevation band in a simulation was executed for 11 years (1990-2000). The subbasin as a function of the lapse rate and the first two years were not used for model evaluation difference between elevation of the meteorological because parameters such as soil water content and gaging station and the average elevation specified for residue cover are initially not in equilibrium with actual the band (Neitsch et al. 2000). physical conditions (Fontaine et al. 2002; Santhi et al. 2001). Prior to calibration, comparison of annual- Adjustment of Snowmelt/Snow Accumulation average streamflow data to simulated values show an Parameters under-prediction of flow (Figure 2). Parameters in SWAT that simulate snowmelt processes and control the formation of snow were also adjusted to Because SWAT was developed for watersheds in non- create a better match to observed streamflow data. The mountainous terrains, special adjustments were parameters that were modified include a factor that necessary to accurately simulate hydrologic processes accounts for snow pack characteristics and two that are strongly affected by elevation changes empirical factors that account for the melting rate of characteristic of this watershed. Fontaine et al. (2002), snow. who applied SWAT to the mountainous Wind River Basin in Wyoming, discovered that orographic A “lagging factor” accounts for temperature processes were very important. Processes that are characteristics of the snow pack that influence the affected by elevation include evapotranspiration, snow-pack density, snow-pack depth, exposure, and precipitation and snowmelt/snow-formation processes. other factors (Neitsch et al. 2001). As the lagging factor approaches 1.0, the mean air temperature on the current 20 Observed day exerts an increasingly greater influence on the 15 Simulated snow pack temperature, and the snow pack temperature /s) 3 from