Benthic macroinvertebrates assemblages and macrohabitat connectivity in Mediterranean-climate streams of northern California

Abstract. Drought leads to a loss of longitudinal and lateral hydrologic connectivity, which causes direct or indirect changes in stream ecosystem properties. Changes in macrohabitat availability from a rifflepool sequence to isolated pools are among the most conspicuous consequences of connectivity loss. Macroinvertebrate assemblages were compared among 3 distinct stream macrohabitats (riffles [R], pools connected to riffles [Pc], disconnected pools [Pd]) of 19 Mediterranean-climate sites in northern California to examine the influence of loss of habitat resulting from drought disturbance. At the time of sampling, 10 sites were perennial and included R and Pc macrohabitats, whereas 9 sites were intermittent and included only Pd macrohabitats. Taxa richness was more variable in Pd, and taxa richness was significantly lower in Pd than in Pc but not R. These results suggested a decline in richness between Pc and Pd that might be associated with loss of connectivity. Lower Ephemeroptera, Plecoptera, and Trichoptera (EPT) richness relative to Odonata, Coleoptera, and Heteroptera (OCH) richness was observed for Pd than R and Pc macrohabitats. Family composition was more similar between R and Pc than between R or Pc and Pd macrohabitats. This similarity may be associated with greater connectivity between R and Pc macrohabitats. Correspondence analysis indicated that macroinvertebrate composition changed along a gradient from R to Pc and Pd that was related to a perennialintermittent gradient across sites. High variability among macroinvertebrate assemblages in Pd could have been related to variability in the duration of intermittency. In cluster analysis, macroinvertebrate assemblages were grouped by macrohabitat first and then by site, suggesting that the macrohabitat filter had a greater influence on macroinvertebrate assemblages than did local site characteristics. Few taxa were found exclusively in Pc, and this macrohabitat shared numerous taxa with R and Pd, indicating that Pc may act as a bridge between R and Pd during drought. Drought is regarded as a ramp disturbance, but our results suggest that the response of macroinvertebrate assemblages to the loss of hydrological connectivity among macrohabitats is gradual, at least in Mediterranean-climate streams where drying is gradual. However, the changes may be more dramatic in arid and semiarid streams or in Mediterranean-climate streams if drying is rapid.

A new method to estimate the infilling of alluvial sediment of glacial valleys using a sloping local base level

A new method is used to estimate the volumes of sediments of glacial valleys. This method is based on the concept of sloping local base level and requires only a digital terrain model and the limits of the alluvial valleys as input data. The bedrock surface of the glacial valley is estimated by a progressive excavation of the digital elevation model (DEM) of the filled valley area. This is performed using an iterative routine that replaces the altitude of a point of the DEM by the mean value of its neighbors minus a fixed value. The result is a curved surface, quadratic in 2D. The bedrock surface of the Rhone Valley in Switzerland was estimated by this method using the free digital terrain model Shuttle Radar Topography Mission (SRTM) (~92 m resolution). The results obtained are in good agreement with the previous estimations based on seismic profiles and gravimetric modeling, with the exceptions of some particular locations. The results from the present method and those from the seismic interpretation are slightly different from the results of the gravimetric data. This discrepancy may result from the presence of large buried landslides in the bottom of the Rhone Valley.

Method for Estimating the Magnitude and Frequency of Floods at Ungaged Sites on Unregulated Rural Streams in Iowa, HR-268, 1987

This report provides techniques and procedures for estimating the probable magnitude and frequency of floods at ungaged sites on Iowa streams. Physiographic characteristics were used to define the boundaries of five hydrologic regions. Regional regression equations that relate the size of the drainage area to flood magnitude are defined for estimating peak discharges having specified recurrence intervals of 2, 5, 10, 25, 50, and 100 years. Regional regression equations are applicable to sites on streams that have drainage areas ranging from 0.04 to 5,150 square miles provided that the streams are not affected significantly by regulation upstream from the sites and that the drainage areas upstream from the sites are not mostly urban areas. Flood-frequency characteristics for the mainstems of selected rivers are presented in graphs as a function of drainage area.

Estimating Design-Flood Discharges for Streams in Iowa Using Drainage-Basin and Channel-Geometry Characteristics, HR-322, 1993

Drainage-basin and channel-geometry multiple-regression equations are presented for estimating design-flood discharges having recurrence intervals of 2, 5, 10, 25, 50, and 100 years at stream sites on rural, unregulated streams in Iowa. Design-flood discharge estimates determined by Pearson Type-III analyses using data collected through the 1990 water year are reported for the 188 streamflow-gaging stations used in either the drainage-basin or channel-geometry regression analyses. Ordinary least-squares multiple-regression techniques were used to identify selected drainage-basin and channel-geometry regions. Weighted least-squares multiple-regression techniques, which account for differences in the variance of flows at different gaging stations and for variable lengths in station records, were used to estimate the regression parameters. Statewide drainage-basin equations were developed from analyses of 164 streamflow-gaging stations. Drainage-basin characteristics were quantified using a geographic-information-system (GIS) procedure to process topographic maps and digital cartographic data. The significant characteristics identified for the drainage-basin equations included contributing drainage area, relative relief, drainage frequency, and 2-year, 24-hour precipitation intensity. The average standard errors of prediction for the drainage-basin equations ranged from 38.6% to 50.2%. The GIS procedure expanded the capability to quantitatively relate drainage-basin characteristics to the magnitude and frequency of floods for stream sites in Iowa and provides a flood-estimation method that is independent of hydrologic regionalization. Statewide and regional channel-geometry regression equations were developed from analyses of 157 streamflow-gaging stations. Channel-geometry characteristics were measured on site and on topographic maps. Statewide and regional channel-geometry regression equations that are dependent on whether a stream has been channelized were developed on the basis of bankfull and active-channel characteristics. The significant channel-geometry characteristics identified for the statewide and regional regression equations included bankfull width and bankfull depth for natural channels unaffected by channelization, and active-channel width for stabilized channels affected by channelization. The average standard errors of prediction ranged from 41.0% to 68.4% for the statewide channel-geometry equations and from 30.3% to 70.0% for the regional channel-geometry equations. Procedures provided for applying the drainage-basin and channel-geometry regression equations depend on whether the design-flood discharge estimate is for a site on an ungaged stream, an ungaged site on a gaged stream, or a gaged site. When both a drainage-basin and a channel-geometry regression-equation estimate are available for a stream site, a procedure is presented for determining a weighted average of the two flood estimates.

Automatic Tracer-Dilution Method Used for Stage-Discharge Ratings and Streamflow Hydrographs on Small Iowa Streams, HR-285, 1990

An automatic system was designed to concurrently measure stage and discharge for the purpose of developing stage-discharge ratings and high flow hydrographs on small streams. Stage, or gage height, is recorded by an analog-to-digital recorder and discharge is determined by the constant-rate tracer-dilution method. The system measures flow above a base stage set by the user. To test the effectiveness of the system and its components, eight systems, with a variety of equipment, were installed at crest-stage gaging stations across Iowa. A fluorescent dye, rhodamine-WT, was used as the tracer. Tracer-dilution discharge measurements were made during 14 flow periods at six stations from 1986 through 1988 water years. Ratings were developed at three stations with the aid of these measurements. A loop rating was identified at one station during rapidly-changing flow conditions. Incomplete mixing and dye loss to sediment apparently were problems at some stations. Stage hydrographs were recorded for 38 flows at seven stations. Limited data on background fluorescence during high flows were also obtained.

Submerged Vanes for Flow Control and Bank Protection in Streams, HR-255, 1984

This study was conducted for the purpose of evaluating a new concept for a bank-protection structure: The Iowa Vane . The underlying idea involves countering the torque exerted on the primary flow by its curvature and vertical velocity gradient, thereby eliminating or significantly reducing the secondary flow and thus reducing the undermining of the outer banks and the high-velocity attack on it. The new structure consists of an array of short, vertical, submerged vanes installed with a certain orientation on the channel bed. A relatively small number of vanes can produce bend flows which are practically uniform across the channel. The height of the vanes is less than half the water depth, and their angle with the flow direction is of the order of l0 degrees. In this study, design relations have been established. The relations, and the vanes' overall performance, have been tested in a laboratory model under different flow and sediment conditions. The results are used for the design of an Iowa-Vane bank protection structure for a section of East Nishnabotna River along U.S. Highway 34 at Red Oak, Iowa.

Impact of controlled changes in grain size and pore space characteristics on the hydraulic conductivity and spectral induced polarization response of "proxies" of saturated alluvial sediments

Understanding the influence of pore space characteristics on the hydraulic conductivity and spectral induced polarization (SIP) response is critical for establishing relationships between the electrical and hydrological properties of surficial unconsolidated sedimentary deposits, which host the bulk of the world's readily accessible groundwater resources. Here, we present the results of laboratory SIP measurements on industrial-grade, saturated quartz samples with granulometric characteristics ranging from fine sand to fine gravel, which can be regarded as proxies for widespread alluvial deposits. We altered the pore space characteristics by changing (i) the grain size spectra, (ii) the degree of compaction, and (iii) the level of sorting. We then examined how these changes affect the SIP response, the hydraulic conductivity, and the specific surface area of the considered samples. In general, the results indicate a clear connection between the SIP response and the granulometric as well as pore space characteristics. In particular, we observe a systematic correlation between the hydraulic conductivity and the relaxation time of the Cole-Cole model describing the observed SIP effect for the entire range of considered grain sizes. The results do, however, also indicate that the detailed nature of these relations depends strongly on variations in the pore space characteristics, such as, for example, the degree of compaction. The results of this study underline the complexity of the origin of the SIP signal as well as the difficulty to relate it to a single structural factor of a studied sample, and hence raise some fundamental questions with regard to the practical use of SIP measurements as site- and/or sample-independent predictors of the hydraulic conductivity.

Missouri River Studies: Alluvial Morphology and Engineering Soil Classification, HR-1, Progress Report, 1960

The primary purposes of this investigation are: 1) To delineate flood plain deposits with different geologic and engineering properties. 2) To provide basic data necessary for any attempt at stabilizing flood plain deposits. The alluvial valley of the Missouri River adjacent to Iowa was chosen as the logical place to begin this study. The river forms the western boundary of the state for an airline distance of approximately 139 miles; and the flood plain varies from a maximum width of approximately 18 miles (Plates 2 and 3, Sheets 75 and 75L) to approximately 4 miles near Crescent, Iowa (Plate 8, Sheet 66). The area studied includes parts of Woodbury, Monona, Harrison, Pottawattamie, Mills, and Fremont counties in Iowa and parts of Dakota, Thurston, Burt, Washington, Douglas, Sarpy, Cass and Otoe counties in Nebraska. Plate l is an index map of the area under consideration.

Techniques for Estimating Flood-Frequency Discharges for Streams in Iowa, HR-395A, 2001

A statewide study was conducted to develop regression equations for estimating flood-frequency discharges for ungaged stream sites in Iowa. Thirty-eight selected basin characteristics were quantified and flood-frequency analyses were computed for 291 streamflow-gaging stations in Iowa and adjacent States. A generalized-skew-coefficient analysis was conducted to determine whether generalized skew coefficients could be improved for Iowa. Station skew coefficients were computed for 239 gaging stations in Iowa and adjacent States, and an isoline map of generalized-skew-coefficient values was developed for Iowa using variogram modeling and kriging methods. The skew map provided the lowest mean square error for the generalized-skew- coefficient analysis and was used to revise generalized skew coefficients for flood-frequency analyses for gaging stations in Iowa. Regional regression analysis, using generalized least-squares regression and data from 241 gaging stations, was used to develop equations for three hydrologic regions defined for the State. The regression equations can be used to estimate flood discharges that have recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years for ungaged stream sites in Iowa. One-variable equations were developed for each of the three regions and multi-variable equations were developed for two of the regions. Two sets of equations are presented for two of the regions because one-variable equations are considered easy for users to apply and the predictive accuracies of multi-variable equations are greater. Standard error of prediction for the one-variable equations ranges from about 34 to 45 percent and for the multi-variable equations range from about 31 to 42 percent. A region-of-influence regression method was also investigated for estimating flood-frequency discharges for ungaged stream sites in Iowa. A comparison of regional and region-of-influence regression methods, based on ease of application and root mean square errors, determined the regional regression method to be the better estimation method for Iowa. Techniques for estimating flood-frequency discharges for streams in Iowa are presented for determining ( 1) regional regression estimates for ungaged sites on ungaged streams; (2) weighted estimates for gaged sites; and (3) weighted estimates for ungaged sites on gaged streams. The technique for determining regional regression estimates for ungaged sites on ungaged streams requires determining which of four possible examples applies to the location of the stream site and its basin. Illustrations for determining which example applies to an ungaged stream site and for applying both the one-variable and multi-variable regression equations are provided for the estimation techniques.

Methods for Estimating Annual Exceedance-Probability Discharges for Streams in Iowa, Based on Data through Water Year 2010, TR-519, 2013

A statewide study was performed to develop regional regression equations for estimating selected annual exceedance- probability statistics for ungaged stream sites in Iowa. The study area comprises streamgages located within Iowa and 50 miles beyond the State’s borders. Annual exceedanceprobability estimates were computed for 518 streamgages by using the expected moments algorithm to fit a Pearson Type III distribution to the logarithms of annual peak discharges for each streamgage using annual peak-discharge data through 2010. The estimation of the selected statistics included a Bayesian weighted least-squares/generalized least-squares regression analysis to update regional skew coefficients for the 518 streamgages. Low-outlier and historic information were incorporated into the annual exceedance-probability analyses, and a generalized Grubbs-Beck test was used to detect multiple potentially influential low flows. Also, geographic information system software was used to measure 59 selected basin characteristics for each streamgage. Regional regression analysis, using generalized leastsquares regression, was used to develop a set of equations for each flood region in Iowa for estimating discharges for ungaged stream sites with 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities, which are equivalent to annual flood-frequency recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years, respectively. A total of 394 streamgages were included in the development of regional regression equations for three flood regions (regions 1, 2, and 3) that were defined for Iowa based on landform regions and soil regions. Average standard errors of prediction range from 31.8 to 45.2 percent for flood region 1, 19.4 to 46.8 percent for flood region 2, and 26.5 to 43.1 percent for flood region 3. The pseudo coefficients of determination for the generalized leastsquares equations range from 90.8 to 96.2 percent for flood region 1, 91.5 to 97.9 percent for flood region 2, and 92.4 to 96.0 percent for flood region 3. The regression equations are applicable only to stream sites in Iowa with flows not significantly affected by regulation, diversion, channelization, backwater, or urbanization and with basin characteristics within the range of those used to develop the equations. These regression equations will be implemented within the U.S. Geological Survey StreamStats Web-based geographic information system tool. StreamStats allows users to click on any ungaged site on a river and compute estimates of the eight selected statistics; in addition, 90-percent prediction intervals and the measured basin characteristics for the ungaged sites also are provided by the Web-based tool. StreamStats also allows users to click on any streamgage in Iowa and estimates computed for these eight selected statistics are provided for the streamgage.

Evaluation of Design Flood Frequency Methods for Iowa Streams, TR-533, 2009

The objective of this project was to assess the predictive accuracy of flood frequency estimation for small Iowa streams based on the Rational Method, the NRCS curve number approach, and the Iowa Runoff Chart. The evaluation was based on comparisons of flood frequency estimates at sites with sufficiently long streamgage records in the Midwest, and selected urban sites throughout the United States. The predictive accuracy and systematic biases (under- or over-estimation) of the approaches was evaluated based on forty-six Midwest sites and twenty-one urban sites. The sensitivity of several watershed characteristics such as soil properties, slope, and land use classification was also explored. Recommendations on needed changes or refinements for applications to Iowa streams are made.

Biological Sampling and Physical Habitat Assessment Standard Operating Procedure for Iowa Wadeable Streams and Rivers, July 24, 2015

The Iowa Department of Natural Resources uses benthic macroinvertebrate and fish sampling data to assess stream biological condition and the support status of designated aquatic life uses (Wilton 2004; IDNR 2013). Stream physical habitat data assist with the interpretation of biological sampling results by quantifying important physical characteristics that influence a stream’s ability to support a healthy aquatic community (Heitke et al., 2006; Rowe et al. 2009; Sindt et al., 2012). This document describes aquatic community sampling and physical habitat assessment procedures currently followed in the Iowa stream biological assessment program. Standardized biological sampling and physical habitat assessment procedures were first established following a pilot sampling study in 1994 (IDNR 1994a, 1994b). The procedure documents were last updated in 2001 (IDNR 2001a; 2001b). The biological sampling and physical habitat assessment procedures described below are evaluated on a continual basis. Revision of this working document will occur periodically to reflect additional changes.

Biological Assessment of Iowa's Streams and Rivers, 2001

Summary of biological monitoring of Iowa's streams and rivers.

Floods of 2008: How Did Aquatic Life in Streams Fare?

Summary of biological water quality data collected during the Floods of 2008.

Fish Habitat Indicators for the Assessment of Wadeable, Warmwater Streams

Physical habitat characteristics such as stream width, depth, instream cover, and substrate composition are important environmental factors that shape Iowa’s stream fish species assemblages. The Iowa Department of Natural Resources (IDNR) stream biological assessment program collects physical habitat data to help interpret fish assemblage sampling results in order to assess stream health condition and the attainment status of designated aquatic life uses. The quantitative habitat indicators and interpretative guidelines developed in this study are designed for specific applications within the stream bioassessment program. These tools might also be useful to natural resource managers for purposes such as stream habitat improvement prioritization, goal-setting, and performance assessment.