Impacts of Urbanization on Stream Habitat and Fish Across Multiple Spatial Scales

We analyzed the relation of the amount and spatial pattern of land cover with stream fish communities, in-stream habitat, and baseflow in 47 small southeastern Wisconsin, USA, watersheds encompassing a gradient of predominantly agricultural to predominantly urban land uses. The amount of connected impervious surface in the watershed was the best measure of urbanization for predicting fish density, species richness, diversity, and index of biotic integrity (IBI) score; bank erosion; and base flow. However, connected imperviousness was not significantly correlated with overall habitat quality for fish. Nonlinear models were developed using quantile regression to predict the maximum possible number of fish species, IBI score, and base flow for a given level of imperviousness. At watershed connected imperviousness levels less than about 8%, all three variables could have high values, whereas at connected imperviousness levels greater than 12% their values were inevitably low. Connected imperviousness levels between 8 and 12% represented a threshold region where minor changes in urbanization could result in major changes in stream condition. In a spatial analysis, connected imperviousness within a 50-m buffer along the stream or within a 1.6-km radius upstream of the sampling site had more influence on stream fish and base flow than did comparable amounts of imperviousness further away. Our results suggest that urban development that minimizes amount of connected impervious surface and establishes undeveloped buffer areas along streams should have less impact than conventional types of development.     

PREDICTION OF MEAN ANNUAL FLOWS IN THE FRASER RIVER CATCHMENT,BRITISH COLUMBIA1

ABSTRACT: A method of predicting annual flows is presented and is applied to the Fraser River catchment. Statistical tests show the annual flow records to be stationary and aerially independent and can be adequately approximated by Gaussian distributions. Estimates are made of the Gaussian parameters for each subbasin. The spatial variations of these parameters are described by third order trend surfaces. The fitted surfaces can then be used to predict parameters of ungaged basins using the latitude and longitude of the basin centroids. The predicted parametric values are substituted into the Gaussian distribution to generate flows of various return periods.     

BASE FLOW TRENDS IN URBANIZING WATERSHEDS OF THE DELAWARE RIVER BASIN1

Rapid land development is raising concern regarding the ability of urbanizing watersheds to sustain adequate base flow during periods of drought. Long term streamflow records from unregulated watersheds of the lower to middle Delaware River basin are examined to evaluate the impact of urbanization and imperviousness on base flow. Trends in annual base flow volumes, seven‐day low flows, and runoff ratios are determined for six urbanizing watersheds and four reference watersheds across three distinct physiographic regions. Hydrograph separation is used to determine annual base flow and stormflow volumes, and nonparametric trend tests are conducted on the resulting time series. Of the watersheds examined, the expected effects of declining base flow volumes and seven‐day low flows and increasing stormflows are seen in only one watershed that is approximately 20 percent impervious and has been subject to a net water export over the past 15 years. Both interbasin transfers and hydrologic mechanisms are invoked to explain these results. The results show that increases in impervious area may not result in measurable reductions in base flow at the watershed scale.     

A NEW FLASHINESS INDEX: CHARACTERISTICS AND APPLICATIONS TO MIDWESTERN RIVERS AND STREAMS1

ABSTRACT: The term flashiness reflects the frequency and rapidity of short term changes in streamflow, especially during runoff events. Flashiness is an important component of a stream's hydrologic regime. A variety of land use and land management changes may lead to increased or decreased flashiness, often to the detriment of aquatic life. This paper presents a newly developed flashiness index, which is based on mean daily flows. The index is calculated by dividing the pathlength of flow oscillations for a time interval (i.e., the sum of the absolute values of day‐to‐day changes in mean daily flow) by total discharge during that time interval. This index has low interannual variability, relative to most flow regime indicators, and thus greater power to detect trends. Index values were calculated for 515 Midwestern streams for the 27‐year period from 1975 through 2001. Statistically significant increases were present in 22 percent of the streams, primarily in the eastern portion of the study area, while decreases were present in 9 percent, primarily in the western portion. Index values tend to decrease with increasing watershed area and with increasing unit area ground water inputs. Area compensated index values often shift at ecoregion boundaries. Potential index applications include evaluation of programs to restore more natural flow regimes.     

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.     

ASSESSING INSTREAM FLOWS AND RESERVOIR OPERATIONS ON AN EASTERN IDAHO RWER1

ABSTRACT: A methodology for assessing reservoir management was applied to the historical conflict between winter fish and wildilife flows below Island Park Reservoir on Henrys Fork of the Snake River and the fulfillment of storage water rights. The methodology consists of (1) identifying impacts of flow regulation, (2) quantifying relationships among variables affecting physical reservoir fill, and (3) assessing effects of these discharges on the fulfillment of water rights in the context of a larger system of interrelated reservoirs. Winter (storage season) flows are critical to management of fish and wildlife populations below Island Park Dam, but flow regulation has resulted in decreased winter discharge. Allowable winter flows are a function of inflow, length of storage season, reservoir content at the start of storage season, and potential for downstream capture of excess storage season water discharged at Island Park. Modeling results indicate that winter flows in the range of those recommended for fish and wildlife management are attainable during average years but not during years when initial reservoir content is low. The methodology was successful in quantifying information useful to decision makers in a variety of agencies and disciplines and could be applied to solve water management problems on other regulated river systems.     

QUANTIFICATION OF INSTREAM FLOW NEEDS OF A WILD AND SCENIC RIVER FOR WATER RIGHTS LITIGATION1

ABSTRACT: The lower 4 miles of the Red River, a tributary of the Rio Grande in northern New Mexico, was designated as one of the “instant” components of the National Wild and Scenic River System in 1968. The Bureau of Land Management (BLM), as the managing agency of the wild and scenic river, was a participant in a general water rights adjudication of the Red River stream system. The BLM sought a federal reserved water right and asserted a claim to the instream flows necessary to protect and maintain the values of the river. Instream flows are not recognized under New Mexico water law. Instream flow requirements were determined by several methods to quantify the claims made by the United States for a federal reserved water right under the Wild and Scenic Rivers Act. The scenic (aesthetic), recreational, and fish and wildlife values are the purposes for which instream flow requirements were claimed. Since water quality is related to these values, instream flows for waste transport and protection of water quality were also included in the claim. The U.S. Fish and Wildlife Service's Instream Flow Incremental Methodology was used to quantify the relationship between various flow regimes and fish habitat. Experience in this litigation indicates the importance of using state-of-the-art methodologies in quantifying instream flow claims. The incremental methodology held up well under technical and legal scrutiny and is an example of the latest methodology that was applied successfully in an adjudication. On February 23, 1984, the parties involved in the adjudication entered a precedential stipulation recognizing a federal reserved right to instream flows for the Red River component of the National Wild and Scenic River System.     

WATER QUALITY SAMPLING SCHEMES FOR VARIABLE FLOW CANALS AT REMOTE SITES1

ABSTRACT: Growing interest in water quality has resulted in the development of monitoring networks and intensive sampling for various constituents. Common purposes are regulatory, source and sink understanding, and trend observations. Water quality monitoring involves monitoring system design; sampling site instrumentation; and sampling, analysis, quality control, and assurance. Sampling is a process to gather information with the least cost and least error. Various water quality sampling schemes have been applied for different sampling objectives and time frames. In this study, a flow proportional composite sampling scheme is applied to variable flow remote canals where the flow rate is not known a priori. In this scheme, historical weekly flow data are analyzed to develop high flow and low flow sampling trigger volumes for auto‐samplers. The median flow is used to estimate low flow sampling trigger volume and the five percent exceedence probability flow is used for high flow sampling trigger volume. A computer simulation of high resolution sampling is used to demonstrate the comparative bias in load estimation and operational cost among four sampling schemes. Weekly flow proportional composite auto‐sampling resulted in the least bias in load estimation with competitive operational cost compared to daily grab, weekly grab sampling and time proportional auto‐sampling.     

Modeling potential herbicide loss to surface waters on the Swiss plateau

Lack of sufficiently detailed data often limits the applicability of complex transport-reaction models for estimating potential herbicide loss to surface waters. Therefore, there is also a need for simple models that are easy to apply but still capture the main features of the underlying processes.In this study, a simple regression model was developed to assess the vulnerability of catchments in the Swiss Plateau to diffuse herbicide loss to surface waters. The model is designed as a screening tool to rank the catchments in a relative sense and not to calculate Predicted Environmental Concentrations (PEC) of pesticides. The main goal is to capture two dominating factors controlling diffuse herbicide transport into streams and rivers. These factors are herbicide application and fast flow processes that are mainly responsible for herbicide transport. In a first step vulnerability of sites to herbicide loss is estimated based on site-specific conditions irrespective of actual herbicide application. In the second step, this vulnerability assessment is combined with actual herbicide application data to estimate the potential herbicide loss.The fast flow index (FFI), derived from discharge data using a base flow separation method, was applied as a proxy for the amount of fast flow occurring. The influence of catchment attributes (including topographic, climatic and soil data) on the FFI was analyzed using a multiple regression approach based on data from 57 catchments of the Swiss Plateau. By combining regression analysis with mechanistic knowledge, a two factor non-linear model based on river density and soil permeability as dominant input factors was selected as the best model for FFI prediction given the available data. Higher dimensional models had to be excluded because the strong correlation between the potential input factors led to unrealistic dependences while only minimally improving the quality of the fit.The spatial pattern of the predicted FFI as a measure for the vulnerability to diffuse herbicide losses shows a clearly increasing trend from the western to the eastern part of the Swiss Plateau and towards the pre-alpine/alpine regions in the south.In general the pattern of herbicide use corresponds to site conditions typical of a low FFI. However, the spatial analysis revealed exceptions, namely areas in which high actual herbicide use coincides with a high FFI.Despite the uncertainties in the model, this simple approach seems to be useful for supporting site-adapted agricultural practice whenever the higher accuracy of more detailed models is not required or too expensive to achieve. In addition, in combination with data on actual herbicide application, it can support the design of monitoring strategies by identifying critical areas of actual herbicide loss.     

CORRELATION OF ANNUAL PEAK FLOWS FOR PENNSYLVANIA STREAMS1

A common assumption in flood frequency analysis is that annual peak flows are independent events. This study was undertaken to investigate the validity of this assumption with regard to Pennsylvania streams by statistically analyzing the dependence between annual peak flows and to determine if basin carryover effects relate to the degree of dependence. Five tests of dependence, the autocorrelation test, the median crossing test, the turning points test, the rank difference test, and the Spearman rank order serial correlation coefficient test were applied to the series of annual peak flows for 57 streams. Of the 57 streams analyzed, only two exhibited signs of dependence by at least two of the tests performed, and the baseflow component of annual peak flows was found to be unrelated to the degree of dependence exhibited between annual peak flows. It was concluded that the assumption of independence of annual peak flows is valid in flood frequency analysis for Pennsylvania streams.     

EFFECT OF SURFACE MINING ON STORM FLOW AND PEAK FLOW FROM SIX SMALL BASINS IN EASTERN KENTUCKY1

ABSTRACT: Hydrologic records from six small Eastern Kentucky watersheds were analyzed to determine the effect of surface mining on storm flows and peak flows. Average storm flow volumes were not changed by surface mining, whereas average peak flows were increased 36 percent. Peak flow increases were only in the summer. Smaller peak flows are doubled; moderate ones are increased by about a third; peak flows around 100 csm seem to be largely unaffected; and the larger peak flows may have been reduced by surface mining. The maximum annual storm flows, usually in winter or spring, appeared slightly reduced. No time trend in either storm flows or peak flows could be detected in five years of postmining record. Surface mining is not a serious floodwater discharge problem.     

EPIC tile flow and nitrate loss predictions for three Minnesota cropping systems

Subsurface tile drains are a key source of nitrate N (NO3-N) losses to streams in parts of the north central USA. In this study, the Erosion Productivity Impact Calculator (EPIC) model was evaluated by comparing measured vs. predicted tile flow, tile NO3-N loss, soil profile residual NO3-N, crop N uptake, and yield, using 4 yr of data collected at a site near Lamberton, MN, for three crop rotations: continuous corn (Zea mays L.) or CC, corn-soybean [Glycine max (L.) Merr.] or CS, and continuous alfalfa (Medicago sativa L.) or CA. Initially, EPIC was run using standard Soil Conservation Service (SCS) runoff curve numbers (CN2) for CC and CS; monthly variations were accurately tracked for tile flow (r2 = 0.86 and 0.90) and NO3-N loss (r2 = 0.69 and 0.52). However, average annual CC and CS tile flows were underpredicted by -32 and -34%, and corresponding annual NO3-N losses were underpredicted by -11 and -52%. Predicted average annual tile flows and NO3-N losses generally improved following calibration of the CN2; tile flow underpredictions were -9 and - 12%, whereas NO3-N losses were 0.6 and -54%. Adjusting a N parameter further improved predicted CS NO3-N losses. Predicted monthly tile flows and NO3-N losses for the CA simulation compared poorly with observed values (r2 values of 0.27 and 0.19); the annual drainage volumes and N losses were of similar magnitude to those measured. Overall, EPIC replicated the relative impacts of the three cropping systems on N fate.     

Farmyard point discharges and their influence on nutrient and labile carbon dynamics in a second order stream draining through a dairy unit

Two small piped sources deriving from a single farmyard together with the receiving second order stream above and below the farmyard region were sampled over a two-year period. Although not measured directly, observations at the time of sampling suggested that maximum drain flow was about 2% of downstream base flow. Both point sources were flowing on each sampling occasion (~62) and usually had concentrations of phosphorus (P), nitrate (NO(3)-N) and biological oxygen demand (BOD) well above those from the upstream site. Individual sample concentrations ranged over more than two orders of magnitude for most determinants and a large proportion of the total P was present as soluble (inorganic and organic) and therefore labile forms. More than 70% of samples collected at the downstream site had concentrations that were >1.2 times those of the corresponding upstream site. On certain sampling occasions >80% of total dissolved phosphorus (TDP) and >90% of the BOD and NO(3) instantaneous load appeared to originate from the farmyard region with the composition of downstream samples being completely overwhelmed after the passage through the farmyard. Extrapolations using instantaneous loads suggest that the farmyard and adjacent areas contributed on average 25-30% of the total and dissolved annual downstream P load of 3 kg P ha(-1) and 1.7 kg P ha(-1), respectively. There was no clear relationship between the relative proportion of the contaminant loading originating from the farmyard region and hydrological events. This emphasises the potential localised significance that small, highly concentrated, continuous or semi-continuous farmyard sources can impact headwater streams during periods of low stream flow.     

Evaluation of regression methodology with low-frequency water quality sampling to estimate constituent loads for ephemeral watersheds in Texas

Water quality regulation and litigation have elevated the awareness and need for quantifying water quality and source contributions in watersheds across the USA. In the present study, the regression method, which is typically applied to large (perennial) rivers, was evaluated in its ability to estimate constituent loads (NO(3)-N, total N, PO(4)-P, total P, sediment) on three small (ephemeral) watersheds with different land uses in Texas. Specifically, regression methodology was applied with daily flow data collected with bubbler stage recorders in hydraulic structures and with water quality data collected with four low-frequency sampling strategies: random, rise and fall, peak, and single stage. Estimated loads were compared with measured loads determined in 2001-2004 with an autosampler and high-frequency sampling strategies. Although annual rainfall and runoff volumes were relatively consistent within watersheds during the study period, measured annual nutrient and sediment concentrations and loads varied considerably for the cultivated and mixed watersheds but not for the pasture watershed. Likewise, estimated loads were much better for the pasture watershed than the cultivated and mixed landuse watersheds because of more consistent land management and vegetation type in the pasture watershed, which produced stronger correlations between constituent loads and mean daily flow rates. Load estimates for PO(4)-P were better than for other constituents possibly because PO(4)-P concentrations were less variable within storm events. Correlations between constituent concentrations and mean daily flow rate were poor and not significant for all watersheds, which is different than typically observed in large rivers. The regression method was quite variable in its ability to accurately estimate annual nutrient loads from the study watersheds; however, constituent load estimates were much more accurate for the combined 3-yr period. Thus, it is suggested that for small watersheds, regression-based annual load estimates should be used with caution, whereas long-term estimates can be much more accurate when multiple years of concentration data are available. The predictive ability of the regression method was similar for all of the low-frequency sampling strategies studied; therefore, single-stage or random strategies are recommended for low-frequency storm sampling on small watersheds because of their simplicity.     

Multiyear nutrient removal performance of three constructed wetlands intercepting tile drain flows from grazed pastures

Subsurface tile drain flows can be a major s ource of nurient loss from agricultural landscapes. This study quantifies flows and nitrogen and phosphorus yields from tile drains at three intensively grazed dairy pasture sites over 3- to 5-yr periods and evaluates the capacity of constructed wetlands occupying 0.66 to 1.6% of the drained catchments too reduce nutrient loads. Continuous flow records are combined with automated flow-proportional sampling of nutrient concentrations to calculate tile drain nutrient yields and wetland mass removal rates. Annual drainage water yields rangedfrom 193 to 564 mm (16-51% of rainfall) at two rain-fed sites and from 827 to 853 mm (43-51% of rainfall + irrigation) at an irrigated site. Annually, the tile drains exported 14 to 109 kg ha(-1) of total N (TN), of which 58 to 90% was nitrate-N. Constructed wetlands intercepting these flows removed 30 to 369 gTN m(-2) (7-63%) of influent loadings annually. Seasonal percentage nitrate-N and TN removal were negatively associated with wetland N mass loadings. Wetland P removal was poor in all wetlands, with 12 to 115% more total P exported annually overall than received. Annually, the tile drains exported 0.12 to 1.38 kg ha of total P, of which 15 to 93% was dissolved reactive P. Additional measures are required to reduce these losses or provide supplementary P removal. Wetland N removal performance could be improved by modifying drainage systems to release flows more gradually and improving irrigation practices to reduce drainage losses.     

COMPARISON OF REGRESSION AND TIME-SERIES METHODS FOR SYNTHESIZING MISSING STREAMFLOW RECORDS1

ABSTRACT: Regression and time-series techniques have been used to synthesize and predict the stream flow at the Foresta Bridge gage from information at the upstream Pohono Bridge gage on the Merced River near Yosemite National Park. Using the available data from two time periods (calendar year 1979 and water year 1986), we evaluated the two techniques in their ability to model the variation in the observed flows and in their ability to predict stream flow at the Foresta Bridge gage for the 1979 time period with data from the 1986 time period. Both techniques produced reasonably good estimates and forecasts of the flow at the downstream gage. However, the regression model was found to have a significant amount of autocorrelation in the residuals, which the time-series model was able to eliminate. The time-series technique presented can be of great assistance in arriving at reasonable estimates of flow in data sets that have large missing portions of data.     

WATER QUALITY IN AGRICULTURAL WATERSHEDS: IMPACT OF RIPARIAN VEGETATION DURING BASE FLOW1

ABSTRACT: Water quality was monitored for 17 months during base flow periods in six agricultural watersheds to evaluate the impact of riparian vegetation on suspended solids and nutrient concentrations. In areas without riparian vegetation, both instream algal production and seasonal low flows appeared to be major determinants of suspended solids, turbidity, and phosphorus concentrations. Peak levels of all parameters were reached during the summer when flows were reduced and benthic algal production was high. Similar summer peaks were reached in streams receiving major point inputs but peaks occurred downstream from the input. Instream organic production was less important in regulating water quality in areas with riparian vegetation and permanent flows. Concentrations of suspended solids remained relatively constant, while phosphorus and turbidity increased in association with leaf fall in autumn. Intermittent flow conditions in summer increased the importance of instream organic production in controlling water quality, even when riparian vegetation was present. Efforts to improve water quality in agricultural watersheds during base flow should emphasize maintenance of riparian vegetation and stable flow conditions.     

A DECISION MODEL TO PREDICT SEDIMENT YIELD FROM FOREST PRACTICES1

ABSTRACT: In order to choose among “best management practices,” forest managers need to predict sediment yield to perennial streams following various forest land operations. The “universal soil loss equation” (USLE) is not directly applicable to forest operations because of the heterogenous soil surface conditions left by harvesting, site preparation, and planting. A sediment hazard index is proposed, to be based on the amount of exposed mineral soil and its proximity to streams. The model offered includes rainfall erosivity, soil erodibility and average land slope, together with the index W. A paired watershed experiment in the central Georgia Piedmont was used to estimate parameters in the model. The experimental basin (80 acres) was clearcut, drum roller chopped twice, and planted by machine. The standard error of estimate of sediment yield was computed to be about 50/lbs/ac per sampling period (four months). Use of William's erogivity index (storm flow times peak flow) reduced the standard error to 33/lbs/ac.     

Prediction Intervals for Estimates of Site Index Based on Ecosystem Type

/ British Columbia has an ecosystem classification system that classifies sites into site series. Foresters commonly measure the productivity of these sites by their site index. In British Columbia, site index is defined as the height of a stand at breast height age 50 and is usually estimated from height-age models. Biogeoclimatic site series/site index relationships are an increasingly popular method of estimating site index in British Columbia for stands where site index cannot be reliably estimated with height-age models. The precision of the predicted site index from these relationships can be evaluated with prediction intervals. This is done for the predicted site index of a single site, a group of sites, or the areally weighted site index of a group of sites. The methodology is also useful in determining the number of sites required to meet a specified precision. These prediction intervals will assist foresters in making sound forest management decisions.KEY WORDS: Biogeoclimatic Ecosystem Classification; Precision; Prediction interval; Site index; Site series