NUTRIENT LOADING ASSESSMENT IN THE ILLINOIS RIVER USING A SYNTHETIC APPROACH1

ABSTRACT: A synthetic relationship is developed between nutrient concentrations and discharge rates at two river gauging sites in the Illinois River Basin. Analysis is performed on data collected by the U.S. Geological Survey (USGS) on nutrients in 1990 through 1997 and 1999 and on discharge rates in 1988 through 1997 and 1999. The Illinois River Basin is in western Arkansas and northeastern Oklahoma and is designated as an Oklahoma Scenic River. Consistently high nutrient concentrations in the river and receiving water bodies conflict with recreational water use, leading to intense stakeholder debate on how best to manage water quality. Results show that the majority of annual phosphorus (P) loading is transported by direct runoff, with high concentrations transported by high discharge rates and low concentrations by low discharge rates. A synthetic relationship is derived and used to generate daily phosphorus concentrations, laying the foundation for analysis of annual loading and evaluation of alternative management practices. Total nitrogen (N) concentration does not have as clear a relationship with discharge. Using a simple regression relationship, annual P loadings are estimated as having a root mean squared error (RMSE) of 39.8 t/yr and 31.9 t/yr and mean absolute percentage errors of 19 percent and 28 percent at Watts and Tahlequah, respectively. P is the limiting nutrient over the full range of discharges. Given that the majority of P is derived from Arkansas, management practices that control P would have the most benefit if applied on the Arkansas side of the border.     

VALIDATION OF AGNPS FOR SMALL WATERSHEDS USING AN INTEGRATED AGNPS/GIS SYSTEM1

ABSTRACT: The AGNPS (AGricultural NonPoint Source) model was evaluated for predicting runoff and sediment delivery from small watersheds of mild topography. Fifty sediment yield events were monitored from two watersheds and five nested subwater-sheds in East Central Illinois throughout the growing season of four years. Half of these events were used to calibrate parameters in the AGNPS model. Average calibrated parameters were used as input for the remaining events to obtain runoff and sediment yield data. These data were used to evaluate the suitability of the AGNPS model for predicting runoff and sediment yield from small, mild-sloped watersheds. An integrated AGNPS/GIS system was used to efficiently create the large number of data input changes necessary to this study. This system is one where the AGNPS model was integrated with the GRASS (Geographic Resources Analysis Support System) GIS (Geographical Information System) to develop a decision support tool to assist with management of runoff and erosion from agricultural watersheds. The integrated system assists with the development of input GIS layers to AGNPS, running the model, and interpretation of the results.     

ADJUSTMENT OF STREAM CHANNEL CAPACITY FOLLOWING DAM CLOSURE,YEGUA CREEK,TEXAS1

ABSTRACT: In Yegua Creek, a principal tributary of the Brazos River in Texas, surveys of a 19 km channel reach downstream of Somerville Dam show that channel capacity decreased by an average of 65 percent in a 34 year period following dam closure. The decrease corresponds with an approximately 85 percent reduction in annual flood peaks. Channel depth has changed the most, decreasing by an average of 61 percent. Channel width remained stable with an average decrease of only 9 percent, reflecting cohesive bank materials along with the growth of riparian vegetation resulting from increased low flows during dry summer months. Although large changes in stream channel geometry are not uncommon downstream of dams, such pronounced reductions in channel capacity could have long‐term implications for sediment delivery through the system.     

SEDIMENT DISPERSION IN SALMON SPAWNING STREAMS: THE INFLUENCE OF FLOODS AND SALMON REDD CONSTRUCTION1

ABSTRACT: Magnetically tagged particles were used to investigate the effects of sockeye salmon (Oncorhynchus nerka) and floods on the dispersion of coarse bed material in the Stuart‐Takla region, British Columbia, Canada. The dominant annual sediment transporting event in the channels is the snowmelt flood events, with lesser activity usually accomplished during summer floods. Annually in August, the channel bed material is reworked by the Early Stuart sockeye salmon spawning, as the fish excavate the streambed to deposit and bury their eggs. These nesting cavities are called redds. Results from 67 tracer recovery experiments over five years were highly variable, subject to the magnitude of floods and the returning population of salmon. Overall, the depositional pattern from nival flood events usually demonstrated a high degree of clast mobilization, long travel distances (up to 150 m), and mean depths of burial up to 18 cm. Summer flood events showed somewhat lower rates of mobilization, distances of travel, and depths of burial. Although the fish did not move the tracers very far, their effect on the bed was generally quite pervasive ‐ up to 100 percent of the clasts were mobilized, and the depth of burial was considerable (mean burial depths up to 14 cm). The amount of vertical mixing of sediment by salmon was often on the same order of magnitude as flood events. The significant vertical mixing of sediments by the fish has important implications for the mobility of sediment in the stream. Since any armoring layer formed during high flows throughout the year is subject to the bioturbation of salmon, this suggests that the transport threshold in the creeks remains relatively low. Salmon likely play an integral role in the sediment transport dynamics and annual sediment budget of the lower reaches of these creeks.     

EFFECTS OF DAM IMPOUNDMENT ON THE FLOOD REGIME OF NATURAL FLOODPLAIN COMMUNITIES IN THE UPPER CONNECTICUT RIVER1

ABSTRACT: Understanding the effects of dams on the inundation regime of natural floodplain communities is critical for effective decision making on dam management or dam removal. To test the implications of hydrologic alteration by dams for floodplain natural communities, we conducted a combined field and modeling study along two reaches in the Connecticut River Rapids Macrosite (CRRM), one of the last remaining flowing water sections of the Upper Connecticut River. We surveyed multiple channel cross sections at both locations and concurrently identified and surveyed the elevations of important natural communities, native species of concern, and nonnative invasive species. Using a hydrologic model, HEC‐RAS, we routed estimated pre‐and post‐impoundment discharges of different design recurrence intervals (two year through 100 year floods) through each reach to establish corresponding reductions in elevation and effective wetted perimeter following post‐dam discharge reductions. By comparing (1) the frequency and duration of flooding of these surfaces before and after impoundment and (2) the total area flooded at different recurrence intervals, our goal was to derive a spatially explicit assessment of hydrologic alteration, directly relevant to natural floodplain communities. Post‐impoundment hydrologic alteration profoundly affected the subsequent inundation regime, and this impact was particularly true of higher floodplain terraces. These riparian communities, which were flooded, on average, every 20 to 100 years pre‐impoundment, were predicted to flood at 100 ? 100 year intervals, essentially isolating them completely from riverine influence. At the pre‐dam five to ten year floodplain elevations, we observed smaller differences in predicted flood frequency but substantial differences in the total area flooded and in the average flood duration. For floodplain forests in the Upper Connecticut River, this alteration by impoundment suggests that even if other stresses facing these communities (human development, invasive exotics) were alleviated, this may not be sufficient to restore intact natural communities. More generally, our approach provides a way to combine site specific variables with long term gage records in assessing the restorative potential of dam removal.     

HYDROLOGICAL MODELING OF THE IROQUOIS RIVER WATERSHED USING HSPF AND SWAT1

ABSTRACT: The performance of two popular watershed scale simulation models — HSPF and SWAT — were evaluated for simulating the hydrology of the 5,568 km² Iroquois River watershed in Illinois and Indiana. This large, tile drained agricultural watershed provides distinctly different conditions for model comparison in contrast to previous studies. Both models were calibrated for a nine‐year period (1987 through 1995) and verified using an independent 15‐year period (1972 through 1986) by comparing simulated and observed daily, monthly, and annual streamflow. The characteristics of simulated flows from both models are mostly similar to each other and to observed flows, particularly for the calibration results. SWAT predicts flows slightly better than HSPF for the verification period, with the primary advantage being better simulation of low flows. A noticeable difference in the models' hydrologic simulation relates to the estimation of potential evapotranspiration (PET). Comparatively low PET values provided as input to HSPF from the BASINS 3.0 database may be a factor in HSPF's overestimation of low flows. Another factor affecting baseflow simulation is the presence of tile drains in the watershed. HSPF parameters can be adjusted to indirectly account for the faster subsurface flow associated with tile drains, but there is no specific tile drainage component in HSPF as there is in SWAT. Continued comparative studies such as this, under a variety of hydrologic conditions and watershed scales, provide needed guidance to potential users in model selection and application.     

NUTRIENT CONTRIBUTION OF NONPOINT SOURCE RUNOFF IN THE LAS VEGAS VALLEY1

ABSTRACT: A Geographic Information System (GIS) based non‐point source runoff model is developed for the Las Vegas Valley, Nevada, to estimate the nutrient loads during the years 2000 and 2001. The estimated nonpoint source loads are compared with current wastewater treatment facilities loads to determine the non‐point source contribution of total phosphorus (TP), total nitrogen (TN), and total suspended solids (TSS) on a monthly and annual time scale. An innovative calibration procedure is used to estimate the pollutant concentrations for different land uses based on available water quality data at the outlet. Results indicate that the pollutant concentrations are higher for the Las Vegas Valley than previous published values for semi‐arid and arid regions. The total TP and TN loads from nonpoint sources are approximately 15 percent and 4 percent, respectively, of the total load to the receiving water body, Lake Mead. The TP loads during wet periods approach the permitted loads from the wastewater treatment plants that discharge into Las Vegas Wash. In addition, the GIS model is used to track pollutant loads in the stream channels for one of the subwatersheds. This is useful for planning the location of Best Management Practices to control nonpoint pollutant loads.     

PARAMETER IDENTIFICATION IN A TWO‐MULTIPLIER SEDIMENT YIELD MODEL1

ABSTRACT: A process based, distributed runoff erosion model (KINEROS2) was used to examine problems of parameter identification of sediment entrainment equations for small watersheds. Two multipliers were used to reflect the distributed nature of the sediment entrainment parameters: one multiplier for a raindrop induced entrainment parameter, and one multiplier for a flow induced entrainment parameter. The study was conducted in three parts. First, parameter identification was studied for simulated error free data sets where the parameter values were known. Second, the number of data points in the simulated sedigraphs was reduced to reflect the effect of temporal sampling frequency on parameter identification. Finally, event data from a small range‐land watershed were used to examine parameter identifiability when the parameter values are unknown. Results demonstrated that whereas unique multiplier values can be obtained for simulated error free data, unique parameter values could not be obtained for some event data. Unique multiplier values for raindrop induced entrainment and flow induced entrainment were found for events with greater than a two‐year return period (～25 mm) that also had at least 10 mm of rain in ten minutes. It was also found that the three‐minute sampling frequency used for the sediment sampler might be inadequate to identify parameters in some cases.     

APPLICATION OF THE BASINS DATABASE AND NPSM MODEL ON A SMALL OHIO WATERSHED1

ABSTRACT: The purpose of this study was to evaluate the Better Assessment Science Integrating Point and Nonpoint Sources (BASINS) watershed management system. BASINS data were used with the NPSM model to predict discharge and sediment concentrations at the outlet of a 103 km² Ohio watershed. It was concluded that the NPSM model should always be calibrated but only a few of the parameters provided with BASINS needed to be calibrated. For a three‐year study period, there was a 2 percent underestimation of discharge using area weighted precipitation values and a 25 percent overestimation using the single station data in BASINS. A comparison of observed and predicted monthly discharge resulted in an r² of 0.86 with area‐weighted precipitation and an r² of 0.74 with the single station data. Calibrating the model to substantially improve sediment predictions was unsuccessful and we concluded that a calibration period of one year was too short. For the three‐year study period, the r² for sediment was 0.36 with a slope of 0.37 and an intercept of 18.8 mg/l. The mean observed and predicted sediment concentrations were 27.1 mg/l and 22.6 mg/l, respectively.     

PROCESS DOMAINS AND THE RIVER CONTINUUM1

ABSTRACT: The concept of process domains is proposed as an alternative to the River Continuum Concept for the influence of geomorphic processes on aquatic ecosystems. Broadly defined, the Process Domain Concept is a multi-scale hypothesis that spatial variability in geomorphic processes governs temporal patterns of disturbances that influence ecosystem structure and dynamics. At a coarse scale, regional climate, geolog vegetation, and topography control the suite of geomorphic processes that are distributed over a landscape. Within the broad context so defined, stream channel classification can guide identification of functionally similar portions of a channel network, but the response of otherwise similar reaches can depend upon their geologic and geomorphic context. Within geomorphic provinces defined by differences in topography, climate history, and tectonic setting, areas with generally similar geology and topography define lithotopo units, which are useful for stratifying different suites of dominant geomorphic processes. Process domains are spatially identifiable areas characterized by distinct suites of geomorphic processes, and the Process Domain Concept implies that channel networks can be divided into discrete regions in which community structure and dynamics respond to distinctly different disturbance regimes. The concepts of process domains and lithotopo units provide both a framework for the application of patch dynamics concepts to complex landscapes and a context for addressing the effects of watershed processes on the ecology of mountain drainage basins.     

RELATION OF NITRATE CONCENTRATIONS TO BASEFLOW IN THE RACCOON RIVER,IOWA1

ABSTRACT: Excessive nitrate‐nitrogen (nitrate) export from the Raccoon River in west central Iowa is an environmental concern to downstream receptors. The 1972 to 2000 record of daily streamflow and the results from 981 nitrate measurements were examined to describe the relation of nitrate to streamflow in the Raccoon River. No long term trends in streamflow and nitrate concentrations were noted in the 28‐year record. Strong seasonal patterns were evident in nitrate concentrations, with higher concentrations occurring in spring and fall. Nitrate concentrations were linearly related to streamflow at daily, monthly, seasonal, and annual time scales. At all time scales evaluated, the relation was improved when baseflow was used as the discharge variable instead of total streamflow. Nitrate concentrations were found to be highly stratified according to flow, but there was little relation of nitrate to streamflow within each flow range. Simple linear regression models developed to predict monthly mean nitrate concentrations explained as much as 76 percent of the variability in the monthly nitrate concentration data for 2001. Extrapolation of current nitrate baseflow relations to historical conditions in the Raccoon River revealed that increasing baseflow over the 20th century could account for a measurable increase in nitrate concentrations.     

WAVE INDUCED FLOW AND TRANSPORT IN SEDIMENT BEDS1

ABSTRACT: The fate of contaminants in large water bodies is highly influenced by the transfer of flow and solutes across the water sediment interface. In this paper, an analytical model is presented where flow in both sediment bed and open channel is coupled at the interface through a boundary layer occupying the upper part of the sediment bed. The presence of this layer allows not only the capture of the inertia effects through a drag term in the generalized Darcy's equation, but also the specification of different soil parameters for the two porous zones. The flow is advective and driven by wave action along the water surface. The resulting system is solved for the pressure and flux in each sediment layer. The generated transport velocity fields are linked to a random walk simulation that is used to examine the trajectories of solute particles. Comparison of these trajectories against experimental tracer tests suggests a pattern very similar to the one attributed to the presence of surface mounds. The results clearly show the significance of the boundary layer and the drag term for soil with high permeability and the impact of both the thickness of the boundary layer and the length of the gravity wave relative to the depth of the water channel on the transport and exchange across the interface. The paper also examines the sensitivity of the mass exchange to the permeability of the two porous zones.     

THE CONCEPT OF HYDROLOGIC LANDSCAPES1

ABSTRACT: Hydrologic landscapes are multiples or variations of fundamental hydrologic landscape units. A fundamental hydrologic landscape unit is defined on the basis of land‐surface form, geology, and climate. The basic land‐surface form of a fundamental hydrologic landscape unit is an upland separated from a lowland by an intervening steeper slope. Fundamental hydrologic landscape units have a complete hydrologic system consisting of surface runoff, ground‐water flow, and interaction with atmospheric water. By describing actual landscapes in terms of land‐surface slope, hydraulic properties of soils and geologic framework, and the difference between precipitation and evapotranspiration, the hydrologic system of actual landscapes can be conceptualized in a uniform way. This conceptual framework can then be the foundation for design of studies and data networks, syntheses of information on local to national scales, and comparison of process research across small study units in a variety of settings. The Crow Wing River watershed in central Minnesota is used as an example of evaluating stream discharge in the context of hydrologic landscapes. Lake‐research watersheds in Wisconsin, Minnesota, North Dakota, and Nebraska are used as an example of using the hydrologic‐land‐scapes concept to evaluate the effect of ground water on the degree of mineralization and major‐ion chemistry of lakes that lie within ground‐water flow systems.     

INFLUENCE OF VEGETATION,ROCK COVER,AND TRAMPLING ON INFILTRATION RATES AND SEDIMENT PRODUCTION1

ABSTRACT: This study evaluated the impact of selected soil surface characteristics on infiltration rates and sediment production from interrill erosion from loam soil. Treatments were two different grass species (crested wheatgrass and intermediate wheatgrass), three levels of grass cover (30, 50, and 80 percent), four levels of rock cover (5, 10, 15, and 20 percent), and six levels of simulated trampling (10 to 60 percent of the respective plot area by 10 percent increments). Results indicated that plots with sod forming grass infiltrated only slightly more water than plots with bunchgrass, though the differences were significant. Trampling reduced infiltration rates significantly. On uncompacted soil, infiltration rates increased as percentage of rock cover increased. Trampling gradually destroyed this relationship however. Rock cover did not significantly affect sediment production. The tradeoff between vegetal cover and rock cover was affected by simulated trampling. Once trampling disturbance reached 20 percent, no relationship between vegetal cover and rock cover existed. Trampling was the most important factor influencing infiltration rates, explaining 35 to 48 percent of the variation in infiltration rates. The most important factor influencing sediment production was grass cover, which explained 40 to 62 percent of the variations associated with sediment yield at various trampling percentages. Results strongly suggest that, for slopes and soils as used here, adequate watershed protection may be obtained by maintaining 50 percent protective ground cover. Additional validation studies are recommended.     

INTERACTIONS BETWEEN GROUND WATER AND SURFACE WATER IN THE SUWANNEE RIVER BASIN,FLORIDA1

ABSTRACT: Ground water and surface water constitute a single dynamic system in most parts of the Suwannee River basin due to the presence of karst features that facilitate the interaction between the surface and subsurface. Low radon-222 concentrations (below background levels) and enriched amounts of oxygen-18 and deuterium in ground water indicate mixing with surface water in parts of the basin. Comparison of surface water and regional ground water flow patterns indicate that boundaries for ground water basins typically do not coincide with surface water drainage subbasins. There are several areas in the basin where ground water flow that originates outside of the Suwannee River basin crosses surface water basin boundaries during both low-flow and high-flow conditions. In a study area adjacent to the Suwannee River that consists predominantly of agricultural land use, 18 wells tapping the Upper Floridan aquifer and 7 springs were sampled three times during 1990 through 1994 for major dissolved inorganic constituents, trace elements, and nutrients. During a period of above normal rainfall that resulted in high river stage and high ground water levels in 1991, the combination of increased amounts of dissolved organic carbon and decreased levels of dissolved oxygen in ground water created conditions favorable for the natural reduction of nitrate by denitrification reactions in the aquifer. As a result, less nitrate was discharged by ground water to the Suwannee River.     

ANALYSIS AND MODELING OF LONG-TERM STREAM TEMPERATURES ON THE STEAMBOAT CREEK BASIN,OREGON: IMPLICATIONS FOR LAND USE AND FISH HABITAT1

ABSTRACT: Steamboat Creek basin is an important source of timber and provides crucial spawning and rearing habitat for anadromous steelhead trout (Oncorhynchus mykiss). Because stream temperatures are near the upper limit of tolerance for the survival of juvenile steelhead, the possible long-term effect of clear-cut logging on stream temperatures was assessed. Twenty-year (1969–1989) records of summer stream temperature and flow from four tributaries and two reaches of Steamboat Creek and Boulder Creek (a nearby unlogged watershed) were analyzed. Logging records for the Steamboat Creek basin and air temperature records also were used in the analysis. A time-series model of the components of stream temperature (seasonal cycle of solar radiation, air temperature, streamflow, an autoregressive term of order 1, and a linear trend variable) was fitted to the water-temperature data. The linear trend variable was significant in all the fitted models except Bend Creek (a tributary fed by cool ground-water discharge) and Boulder Creek. Because no trends in either climate (i.e., air temperature) or streamflow were found in the data, the trend variable was associated with the pre-1969 loss and subsequent regrowth of riparian vegetation and shading canopies.     

LARGE AREA HYDROLOGIC MODELING AND ASSESSMENT PART I: MODEL DEVELOPMENT1

ABSTRACT: A conceptual, continuous time model called SWAT (Soil and Water Assessment Tool) was developed to assist water resource managers in assessing the impact of management on water supplies and nonpoint source pollution in watersheds and large river basins. The model is currently being utilized in several large area projects by EPA, NOAA, NRCS and others to estimate the off-site impacts of climate and management on water use, non-point source loadings, and pesticide contamination. Model development, operation, limitations, and assumptions are discussed and components of the model are described. In Part II, a GIS input/output interface is presented along with model validation on three basins within the Upper Trinity basin in Texas.     

DAM‐INDUCED MODIFICATIONS TO UPPER ALLEGHENY RIVER STREAMFLOW PATTERNS AND THEIR BIODIVERSITY IMPLICATIONS1

ABSTRACT: This study evaluates the streamflow characteristics of the upper Allegheny River during the periods preceding (1936 to 1965) and following (1966 to 1997) completion of the Kinzua Dam in northwestern Pennsylvania. Inter‐period trends in seasonal patterns of discharge and peak flow at three downstream sites are compared to those at two upstream sites to determine the influence of this large dam on surface water hydrology. Climatic records indicate that significant changes in annual total and seasonal precipitation occurred over the twentieth century. Increased runoff during the late summer through early winter led to increased discharge both upstream and downstream during these months, while slightly less early‐year rainfall produced minor reductions in spring flood peaks since 1966. The Kinzua Dam significantly enhanced these trends downstream, creating large reductions in peak flow, while greatly augmenting low flow during the growing season. This reduction in streamflow variability, coupled with other dam‐induced changes, has important biodiversity implications. The downstream riparian zone contains numerous threatened/endangered species, many of which are sensitive to the type of habitat modifications produced by the dam. Flood dynamics under the current post‐dam conditions are likely to compound the difficulties of maintaining their long‐term viability.     

APPLICATION OF THE SNOWMELT RUNOFF MODEL USING MULTIPLE-PARAMETER LANDSCAPE ZONES ON THE TOWANDA CREEK BASIN,PENNSYLVANIA1

ABSTRACT: The Snowmelt Runoff Model (SRM) is designed to compute daily stream discharge using satellite snow cover data for a basin divided into elevation zones. For the Towanda Creek basin, a Pennsylvania watershed with relatively little relief, analysis of snow cover images revealed that both elevation and land use affected snow accumulation and melt on the landscape. The distribution of slope and aspect on the watershed was also considered; however, these landscape features were not well correlated with the available snow cover data. SRM streamflow predictions for 1990, 1993 and 1994 snowmelt seasons for the Towanda Creek basin using a combination of elevation and land use zones yielded more precise streamflow estimates than the use of standard elevation zones alone. The use of multiple-parameter zones worked best in non-rain-on-snow conditions such as in 1990 and 1994 seasons where melt was primarily driven by differences in solar radiation. For seasons with major rain-on-snow events such as 1993, only modest improvements were shown since melt was dominated by rainfall energy inputs, condensation and sensible heat convection. Availability of GIS coverages containing satellite snow cover data and other landscape attributes should permit similar reformulation of multiple-parameter watershed zones and improved SRM streamflow predictions on other basins.     

THE VULNERABILITY OF WETLANDS TO CLIMATE CHANGE: A HYDROLOGIC LANDSCAPE PERSPECTIVE1

ABSTRACT: The vulnerability of wetlands to changes in climate depends on their position within hydrologic landscapes. Hydrologic landscapes are defined by the flow characteristics of ground water and surface water and by the interaction of atmospheric water, surface water, and ground water for any given locality or region. Six general hydrologic landscapes are defined; mountainous, plateau and high plain, broad basins of interior drainage, riverine, flat coastal, and hummocky glacial and dune. Assessment of these landscapes indicate that the vulnerability of all wetlands to climate change fall between two extremes: those dependent primarily on precipitation for their water supply are highly vulnerable, and those dependent primarily on discharge from regional ground water flow systems are the least vulnerable, because of the great buffering capacity of large ground water flow systems to climate change.