Challenging established narratives on soil erosion and shifting cultivation in Laos

The official environmental discourse in Laos describes a “chain of degradation” stretching from upland shifting cultivation, increased runoff and soil erosion to the siltation of wetlands and reservoirs. This perspective has had wide‐ranging impacts on rural development policy which, in the uplands, has long favoured forest conservation over agriculture. Integrating soil erosion and water sediment data with local perceptions of land degradation in an upland village of northern Laos, this study tests the validity of the official environmental discourse. Biophysical measurements made in a small agricultural catchment indicate a significant correlation between the spatial extent of cultivation and soil erosion rates. However, sediment yields recorded at the outlet of the catchment highlight relatively low levels of off‐site sediment exportation. Furthermore, farmers' perceptions suggest that local land degradation issues and crop yield declines could be less related to soil erosion than to agricultural land shortage, increased weed competition, and fertility losses resulting from the intensification of shifting cultivation. The study concludes that a better understanding and management of land degradation issues can be achieved by developing more inclusive and scientifically‐informed approaches to environmental perceptions and narratives.     

Use of Streamflow Data to Estimate Base Flow/Ground‐Water Recharge For Wisconsin1

Abstract: The average annual base flow/recharge was determined for streamflow‐gaging stations throughout Wisconsin by base‐flow separation. A map of the State was prepared that shows the average annual base flow for the period 1970‐99 for watersheds at 118 gaging stations. Trend analysis was performed on 22 of the 118 streamflow‐gaging stations that had long‐term records, unregulated flow, and provided aerial coverage of the State. The analysis found that a statistically significant increasing trend was occurring for watersheds where the primary land use was agriculture. Most gaging stations where the land cover was forest had no significant trend. A method to estimate the average annual base flow at ungaged sites was developed by multiple‐regression analysis using basin characteristics. The equation with the lowest standard error of estimate, 9.5%, has drainage area, soil infiltration and base flow factor as independent variables. To determine the average annual base flow for smaller watersheds, estimates were made at low‐flow partial‐record stations in 3 of the 12 major river basins in Wisconsin. Regression equations were developed for each of the three major river basins using basin characteristics. Drainage area, soil infiltration, basin storage and base‐flow factor were the independent variables in the regression equations with the lowest standard error of estimate. The standard error of estimate ranged from 17% to 52% for the three river basins.     

URBANIZATION INFLUENCES ON AQUATIC COMMUNITIES IN NORTHEASTERN ILLINOIS STREAMS1

ABSTRACT: Biotic indices and sediment trace element concentrations for 43 streams in northeastern Illinois (Chicago area) from the 1980s and 1990s were examined along an agricultural to urban land cover gradient to explore the relations among biotic integrity, sediment chemistry, and urbanization. The Illinois fish Alternative Index of Biotic Integrity (AIBI) ranged from poor to excellent in agricultural/rural streams, but streams with more than 10 percent watershed urban land (about 500 people/mi²) had fair or poor index scores. A macroinvertebrate index (MBI) showed similar trends. A qualitative habitat index (PIBI) did not correlate to either urban indicator. The AIBI and MBI correlated with urban associated sediment trace element concentrations. Elevated copper concentrations in sediment occurred in streams with greater than 40 percent watershed urban land. The number of intolerant fish species and modified index of biotic integrity scores increased in some rural, urbanizing, and urban streams from the 1980s to 1990s, with the largest increases occurring in rural streams with loamy/sandy surficial deposits. However, smaller increases also occurred in urban streams with clayey surficial deposits and over 50 percent watershed urban land. These data illustrate the potentially complex spatial and temporal relations among biotic integrity, sediment chemistry, watershed urban land, population density, and regional and local geologic setting.     

RIPARIAN LAND USES AND PRECIPITATION INFLUENCES ON STREAM BANK EROSION IN CENTRAL IOWA1

Human alterations to the Iowa landscape, such as elimination of native vegetation for row crop agriculture and grazing, channelization of streams, and tile and ditch drainage, have led to deeply incised channels with accelerated streambank erosion. The magnitude of streambank erosion and soil loss were compared along Bear Creek in central Iowa. The subreaches are bordered by differing land uses, including reestablished riparian forest buffers, row crop fields, and continuously grazed riparian pastures. Erosion pins were measured from June 1998 to July 2002 to estimate the magnitude of streambank erosion. Total streambank soil loss was estimated by using magnitude of bank erosion, soil bulk density, and severely eroded bank area. Significant seasonal and yearly differences in magnitude of bank erosion and total soil loss were partially attributed to differences in precipitation and associated discharges. Riparian forest buffers had significantly lower magnitude of streambank erosion and total soil loss than the other two riparian land uses. Establishment of riparian forest buffers along all of the nonbuffered subreaches would have reduced stream‐bank soil loss by an estimated 77 to 97 percent, significantly decreasing sediment in the stream, a major water quality problem in Iowa.     

Research Gaps Related to Forest Management and Stream Sediment in the United States

Water quality from forested landscapes tends to be very high but can deteriorate during and after silvicultural activities. Practices such as forest harvesting, site preparation, road construction/use, and stream crossings have been shown to contribute sediment, nutrients, and other pollutants to adjacent streams. Although advances in forest management accompanied with Best Management Practices (BMPs) have been very effective at reducing water quality impacts from forest operations, projected increases in demand for forest products may result in unintended environmental degradation. Through a review of the pertinent literature, we identified several research gaps related to water yield, aquatic habitat, sediment source and delivery, and BMP effectiveness that should be addressed for streams in the United States to better understand and address the environmental ramifications of current and future levels of timber production. We explored the current understanding of these topics based on relevant literature and the possible implications of increased demand for forest products in the United States.     

Determination of Soil Erosion Risk in the Mustafakemalpasa River Basin, Turkey, Using the Revised Universal Soil Loss Equation, Geographic Information System, and Remote Sensing

Sediment transport from steep slopes and agricultural lands into the Uluabat Lake (a RAMSAR site) by the Mustafakemalpasa (MKP) River is a serious problem within the river basin. Predictive erosion models are useful tools for evaluating soil erosion and establishing soil erosion management plans. The Revised Universal Soil Loss Equation (RUSLE) function is a commonly used erosion model for this purpose in Turkey and the rest of the world. This research integrates the RUSLE within a geographic information system environment to investigate the spatial distribution of annual soil loss potential in the MKP River Basin. The rainfall erosivity factor was developed from local annual precipitation data using a modified Fournier index: The topographic factor was developed from a digital elevation model; the K factor was determined from a combination of the soil map and the geological map; and the land cover factor was generated from Landsat-7 Enhanced Thematic Mapper (ETM) images. According to the model, the total soil loss potential of the MKP River Basin from erosion by water was 11,296,063?Mg?year(-1) with an average soil loss of 11.2?Mg?year(-1). The RUSLE produces only local erosion values and cannot be used to estimate the sediment yield for a watershed. To estimate the sediment yield, sediment-delivery ratio equations were used and compared with the sediment-monitoring reports of the Dolluk stream gauging station on the MKP River, which collected data for >41?years (1964-2005). This station observes the overall efficiency of the sediment yield coming from the Orhaneli and Emet Rivers. The measured sediment in the Emet and Orhaneli sub-basins is 1,082,010?Mg?year(-1) and was estimated to be 1,640,947?Mg?year(-1) for the same two sub-basins. The measured sediment yield of the gauge station is 127.6?Mg?km(-2)?year(-1) but was estimated to be 170.2?Mg?km(-2) year(-1). The close match between the sediment amounts estimated using the RUSLE-geographic information system (GIS) combination and the measured values from the Dolluk sediment gauge station shows that the potential soil erosion risk of the MKP River Basin can be estimated correctly and reliably using the RUSLE function generated in a GIS environment.     

Landscape Planning for Agricultural Nonpoint Source Pollution Reduction III: Assessing Phosphorus and Sediment Reduction Potential

Riparian buffers have the potential to improve stream water quality in agricultural landscapes. This potential may vary in response to landscape characteristics such as soils, topography, land use, and human activities, including legacies of historical land management. We built a predictive model to estimate the sediment and phosphorus load reduction that should be achievable following the implementation of riparian buffers; then we estimated load reduction potential for a set of 1598 watersheds (average 54 km²) in Wisconsin. Our results indicate that land cover is generally the most important driver of constituent loads in Wisconsin streams, but its influence varies among pollutants and according to the scale at which it is measured. Physiographic (drainage density) variation also influenced sediment and phosphorus loads. The effect of historical land use on present-day channel erosion and variation in soil texture are the most important sources of phosphorus and sediment that riparian buffers cannot attenuate. However, in most watersheds, a large proportion (approximately 70%) of these pollutants can be eliminated from streams with buffers. Cumulative frequency distributions of load reduction potential indicate that targeting pollution reduction in the highest 10% of Wisconsin watersheds would reduce total phosphorus and sediment loads in the entire state by approximately 20%. These results support our approach of geographically targeting nonpoint source pollution reduction at multiple scales, including the watershed scale.     

Evaluation of Deposited Sediment and Macroinvertebrate Metrics Used to Quantify Biological Response to Excessive Sedimentation in Agricultural Streams

The objective of this study was to evaluate which macroinvertebrate and deposited sediment metrics are best for determining effects of excessive sedimentation on stream integrity. Fifteen instream sediment metrics, with the strongest relationship to land cover, were compared to riffle macroinvertebrate metrics in streams ranging across a gradient of land disturbance. Six deposited sediment metrics were strongly related to the relative abundance of Ephemeroptera, Plecoptera and Trichoptera and six were strongly related to the modified family biotic index (MFBI). Few functional feeding groups and habit groups were significantly related to deposited sediment, and this may be related to the focus on riffle, rather than reach-wide macroinvertebrates, as reach-wide sediment metrics were more closely related to human land use. Our results suggest that the coarse-level deposited sediment metric, visual estimate of fines, and the coarse-level biological index, MFBI, may be useful in biomonitoring efforts aimed at determining the impact of anthropogenic sedimentation on stream biotic integrity.     

Modeling of Sediment Yield From Anjeni-Gauged Watershed,Ethiopia Using SWAT Model1

Setegn, Shimelis G., Bijan Dargahi, Ragahavan Srinivasan, and Assefa M. Melesse, 2010. Modeling of Sediment Yield From Anjeni-Gauged Watershed, Ethiopia Using SWAT Model. Journal of the American Water Resources Association (JAWRA) 46(3):514-526. DOI: 10.1111/j.1752-1688.2010.00431.x Abstract: The Soil and Water Assessment Tool (SWAT) was tested for prediction of sediment yield in Anjeni-gauged watershed, Ethiopia. Soil erosion and land degradation is a major problem on the Ethiopian highlands. The objectives of this study were to evaluate the performance and applicability of SWAT model in predicting monthly sediment yield and assess the impacts of subbasin delineation and slope discretization on the prediction of sediment yield. Ten years monthly meteorological, flow and sediment data were used for model calibration and validation. The annual average measured sediment yield was 24.6 tonnes/ha. The annual average simulated sediment yield was 27.8 and 29.5 tones/ha for calibration and validation periods, respectively. The study found that the observed values showed good agreement with the simulated sediment yield with Nash-Sutcliffe efficiency (NSE) = 0.81, percent bias (PBIAS) = 28%, RMSE-observations standard deviation ratio (RSR) = 0.23, and coefficient of determination (R²) = 0.86 for calibration and NSE = 0.79, PBIAS = 30%, RSR = 0.29, and R² = 0.84 for validation periods. The model can be used for further analysis of different management scenarios that could help different stakeholders to plan and implement appropriate soil and water conservation strategies.     

ASSESSING STREAM CHANNEL STABILITY THRESHOLDS USING FLOW COMPETENCE ESTIMATES AT BANKFULL STAGE1

ABSTRACT: Stream channel stability is affected by peak flows rather than average annual water yield. Timber harvesting and other land management activities that contribute to soil compaction, vegetation removal, or increased drainage density can increase peak discharges and decrease the recurrence interval of bankfull discharges. Increased peak discharges can cause more frequent movement of large streambed materials, leading to more rapid stream channel change or instability. This study proposes a relationship between increased discharge and channel stability, and presents a methodology that can be used to evaluate stream channel stability thresholds on a stream reach basis. Detailed surveys of the channel cross section, water surface slope, streambed particle size distribution, and field identification of bankfull stage are used to estimate existing bankfull flow conditions. These site specific stream channel characteristics are used in bed load movement formulae to predict critical flow conditions for entrainment of coarse bed material (D₈₄ size fraction). The “relative bed stability” index, defined as the ratio of critical flow condition to the existing condition at bankfull discharge, can predict whether increased peak discharges will exceed stream channel thresholds.     

Estimating Forest Ecosystem Evapotranspiration at Multiple Temporal Scales With a Dimension Analysis Approach1

Abstract:  It is critical that evapotranspiration (ET) be quantified accurately so that scientists can evaluate the effects of land management and global change on water availability, streamflow, nutrient and sediment loading, and ecosystem productivity in watersheds. The objective of this study was to derive a new semi‐empirical ET modeled using a dimension analysis method that could be used to estimate forest ET effectively at multiple temporal scales. The model developed describes ET as a function of water availability for evaporation and transpiration, potential ET demand, air humidity, and land surface characteristics. The model was tested with long‐term hydrometeorological data from five research sites with distinct forest hydrology in the United States and China. Averaged simulation error for daily ET was within 0.5 mm/day. The annual ET at each of the five study sites were within 7% of measured values. Results suggest that the model can accurately capture the temporal dynamics of ET in forest ecosystems at daily, monthly, and annual scales. The model is climate‐driven and is sensitive to topography and vegetation characteristics and thus has potential to be used to examine the compounding hydrologic responses to land cover and climate changes at multiple temporal scales.     

Flow,Organic, and Inorganic Sediment Yields from a Channelized Watershed in the South Carolina Lower Coastal Plain

Many small streams in coastal watersheds in the southeastern United States are modified for agricultural, residential, and commercial development. In the South Carolina Lower Coastal Plain, low‐relief topography and a shallow water table make stream channelization ubiquitous. To quantify the impacts of urbanization and stream channelization, we measured flow and sediment from an urbanizing watershed and a small forested watershed. Flow and sediment export rates were used to infer specific yields from forested and nonforested regions of the urbanizing watershed. Study objectives were to: (1) quantify the range of runoff‐to‐rainfall ratios; (2) quantify the range of specific sediment yields; (3) characterize the quantity and quality of particulate matter exported; and (4) estimate sediment yield attributable to agriculture, development, and channelization activities in the urbanizing watershed. Our results showed that the urban watershed exported over five times more sediment per unit area compared with the forested watershed. Sediment concentration was related to flow flashiness in the urban watershed and to flow magnitude in the forested watershed. Sediments from the forested watershed were dominated by organic matter, whereas mineral matter dominated sediment from the urban stream. Our results indicated that a significant shift in sediment quality and quantity are likely to occur as forested watersheds are transformed by urbanization in coastal South Carolina.     

Magnitude,Frequency, and Duration Relations for Suspended Sediment in Stable (“Reference”) Southeastern Streams1

Abstract: Sediment is listed as one of the leading causes of water‐quality impairments in surface waters of the United States (U.S.). A water body becomes listed by a State, Territory or Tribe if its designated use is not being attained (i.e., impaired). In many cases, the prescribed designated use is aquatic health or habitat, indicating that total maximum daily loads (TMDL) targets for sediment should be functionally related to this use. TMDL targets for sediment transport have been developed for many ecoregions over the past several years using suspended‐sediment yield as a metric. Target values were based on data from “reference” streams or reaches, defined as those exhibiting geomorphic characteristics of equilibrium. This approach has proved useful to some states developing TMDLs for suspended sediment, although one cannot conclude that if a stream exceeds the target range, the aquatic ecosystem will be adversely impacted. To address this problem, historical flow‐transport and sediment‐transport data from hundreds of sites in the Southeastern U.S. were re‐examined to develop parameters (metrics) such as frequency and duration of sediment concentrations. Sites determined as geomorphically stable from field evaluations and from analysis of gauging‐station records were sorted by ecoregion. Mean‐daily flow data obtained from the U.S. Geological Survey were applied to sediment‐transport rating relations to determine suspended‐sediment load for each day of record. The frequency and duration that a given concentration was equaled or exceeded were then calculated to produce a frequency distribution for each site. “Reference” distributions were created using the stable sites in each ecoregion by averaging all of the distributions at specified exceedance intervals. As with the “reference” suspended‐sediment yields, there is a broad range of frequency and duration distributions that reflects the hydrologic and sediment‐transport regimes of the ecoregions. Ecoregions such as the Mississippi Valley Loess Plains (#74) maintain high suspended sediment concentrations for extended periods, whereas coastal plain ecoregions (#63 and 75) show much lower concentrations.     

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.     

A Sediment Budget for an Urbanizing Watershed, 1951‐1996, Montgomery County, Maryland,U.S.A.1

Abstract: Despite widespread interest, few sediment budgets are available to document patterns of erosion and sedimentation in developing watersheds. We assess the sediment budget for the Good Hope Tributary, a small watershed (4.05 km²) in Montgomery County, Maryland, from 1951‐1996. Lacking monitoring data spanning the period of interest, we rely on a variety of indirect and stratigraphic methods. Using regression equations relating sediment yield to construction, we estimated an upland sediment production of 5,700 m³ between 1951 and 1996. Regression equations indicate that channel cross‐sectional area is correlated with the extent of development; these relationships, when combined with historical land use data, suggest that upland sediment yield was augmented by 6,400 m³ produced by enlargement of first‐order and second‐order stream channels. We used dendrochronology to estimate that 4,000 m³ of sediment was stored on the floodplain from 1951‐1996. The sediment yield from the watershed, obtained by summing upstream contributions, totals 8,100 m³ of sediment, or 135 tons/km²/year. These results indicate that upland erosion, channel enlargement, and floodplain storage are all significant components of the sediment budget of our study area, and all three are approximately equal in magnitude. Erosion of “legacy” floodplain sediments originally deposited during poor agricultural practices of the 19th and early 20th Centuries has likely contributed between 0 and 20% of the total sediment yield, indicating that these remobilized deposits are not a dominant component of the sediment yield of our study area.     

“R1-R4” AND “BOISED” SEDIMENT PREDICTION MODEL TESTS USING FOREST ROADS IN GRANITICS1

ABSTRACT: Erosion and sedimentation data from research watersheds in the Silver Creek Study Area in central Idaho were used to test the prediction of logging road erosion using the R1-R4 sediment yield model, and sediment delivery using the “BOISED” sediment yield prediction model. Three small watersheds were instrumented and monitored such that erosion from newly constructed roads and sediment delivery to the mouths of the watersheds could be measured for four years following road construction. The errors for annual surface erosion predictions for the two standard road tests ranged from +31.2 t/ha/yr (+15 percent) to -30.3 t/ha/yr (-63 percent) with an average of zero t/ha/yr and a standard deviation of the differences of 18.7 t/ha/yr. The annual prediction errors for the three watershed scale tests had a greater range from -40.8 t/ha/yr (-70 percent) to +65.3 t/ha/yr (+38 percent) with a mean of -1.9 t/ha/yr and a standard deviation of the differences of 25.2 t/ha/yr. Sediment yields predicted by BOISED (watershed scale tests) were consistently greater (average of 2.5 times) than measured sediment yields. Hillslope sediment delivery coefficients in BOISED appear to be overly conservative to account for average site conditions and road locations, and thus over-predict sediment delivery. Mass erosion predictions from BOISED appear to predict volume well (465 tonnes actual versus 710 tonnes predicted, or a 35 percent difference) over 15 to 20 years, however mass wasting is more episodic than the model predicts.     

Applications of Turbidity Monitoring to Forest Management in California

Many California streams have been adversely affected by sedimentation caused by historic and current land uses, including timber harvesting. The impacts of timber harvesting and logging transportation systems on erosion and sediment delivery can be directly measured, modeled, or inferred from water quality measurements. California regulatory agencies, researchers, and land owners have adopted turbidity monitoring to determine effects of forest management practices on suspended sediment loads and water quality at watershed, project, and site scales. Watershed-scale trends in sediment discharge and responses to current forest practices may be estimated from data collected at automated sampling stations that measure turbidity, stream flow, suspended sediment concentrations, and other water quality parameters. Future results from these studies will provide a basis for assessing the effectiveness of modern forest practice regulations in protecting water quality. At the project scale, manual sampling of water column turbidity during high stream flow events within and downstream from active timber harvest plans can identify emerging sediment sources. Remedial actions can then be taken by managers to prevent or mitigate water quality impacts. At the site scale, manual turbidity sampling during storms or high stream flow events at sites located upstream and downstream from new, upgraded, or decommissioned stream crossings has proven to be a valuable way to determine whether measures taken to prevent post-construction erosion and sediment production are effective. Turbidity monitoring at the project and site scales is therefore an important tool for adaptive management. Uncertainty regarding the effects of current forest practices must be resolved through watershed-scale experiments. In the short term, this uncertainty will stimulate increased use of project and site-scale monitoring.     

Road Erosion and Delivery Index (READI): A Model for Evaluating Unpaved Road Erosion and Stream Sediment Delivery

The Road Erosion and Delivery Index (READI) is a new geographic information system–based model to assess erosion and delivery of water and sediment from unpaved road networks to streams. READI quantifies the effectiveness of existing road surfacing and drain placements in reducing road sediment delivery and guides upgrades to optimize future reductions. Roads are draped on a digital elevation model and parsed into hydrologically distinct segments. Segments are further divided by engineered drainage structures. For each segment, a kinematic wave approximation generates runoff hydrographs for specified storms, with discharge directly to streams at road–stream crossings and onto overland‐flow plumes at other discharge points. Plumes are attenuated by soil infiltration, which limits their length, with delivery occurring if plumes intersect streams. Sediment production and sediment delivery can be calculated as a relative dimensionless index. READI predicts only a small proportion of new drains and new surfacing results in the majority of sediment delivery reductions. The model illustrates how the spatial relationships between road and stream networks, controlled by topography and network geometries, influence patterns of road–stream connectivity. READI was applied in seven northern California basins. The model was also applied in a recent burn area to examine how reduced hillslope infiltration can result in increased hydrologic connectivity and sediment delivery.     

HYSTAR Sediment Model: Distributed Two‐Dimensional Simulation of Watershed Erosion and Sediment Transport Using Time‐Area Routing

An erosion and sediment transport component incorporated in the HYdrology Simulation using Time‐ARea method (HYSTAR) upland watershed model provides grid‐based prediction of erosion, transport and deposition of sediment in a dynamic, continuous, and fully distributed framework. The model represents the spatiotemporally varied flow in sediment transport simulation by coupling the time‐area routing method and sediment transport capacity approach within a grid‐based spatial data model. This avoids the common, and simplistic, approach of using the Universal Soil Loss Equation (USLE) to estimate erosion rates with a delivery ratio to relate gross soil erosion to sediment yield of a watershed, while enabling us to simulate two‐dimensional sediment transport processes without the complexity of numerical solution of the partial differential governing equations. In using the time‐area method for routing sediment, the model offers a novel alternative to watershed‐scale sediment transport simulation that provides detailed spatial representation. In predicting four‐year sediment hydrographs of a watershed in Virginia, the model provided good performance with R² of 0.82 and 0.78 and relative error of ?35% and 11% using the Yalin and Yang's sediment transport capacity equations, respectively. Prediction of spatiotemporal variation in sediment transport processes was evaluated using maps of sediment transport rates, concentrations, and erosion and deposition mass, which compare well with expected behavior of flow hydraulics and sediment transport processes.     

LAND USE AND AQUATIC BIOINTEGRITY IN THE BLACKFOOT RWER WATERSHED,MONTANA1

ABSTRACT: Benthic macroinvertebrate samples representing 151 taxa were collected in August 1995 to examine the linkage between land use, water quality, and aquatic biointegrity in seven tributaries of the Blackfoot River watershed, Montana. The tributaries represent silvicultural (timber harvesting), agricultural (irrigated alfalfa and hay and livestock grazing), and wilderness land uses. A 2.4 km (1.5 mile) reach of a recently restored tributary also was sampled for comparison with the other six sites. A geographic information system (GIS) was used to characterize the seven subwatersheds and estimate soil erosion, using the Modified Universal Soil Loss Equation, and sediment delivery. The wilderness stream had the highest aquatic biointegrity. Two agricultural streams had the largest estimated soil erosion and sediment delivery rates, the greatest habitat impairment from nonpoint source pollution, and the most impoverished macroinvertebrate communities. The silvicultural subwatersheds had greater rates of estimated soil erosion and sediment delivery and lower aquatic biointegrity than the wilderness reference site but evinced better conditions than the agricultural sites. A multiple-use (forestry, grazing, and wildlife management) watershed and the restored site ranked between the silvicultural and agricultural sites. This spectrum of land use and aquatic biointegrity illustrates both the challenges and opportunities that define watershed management.