Bankfull Width Research Articles

An ephemeral meandering river system: Sediment dispersal processes in the Río Colorado, Southern Altiplano Plateau, Bolivia

The Rio Colorado meandering river system feeds the Salar de Uyuni, the World's largest salt pan in the southern Altiplano plateau (Bolivia). It is characterized by ephemerality due to the high aridity of the region, and a downstream decrease of river bankfull width and depth in the river terminus. Dry land meandering river systems in high-altitude regions have rarely been reported. We investigated the active processes in this system from the alluvial fan to the lower coastal plain. On-site surveys combined with high-precision GPS and high-resolution satellite imagery show that Rio Colorado area is typified by generally fine sediment and a low gradient in the coastal plain. In the alluvial fan, where channel meandering is weak but still recognizable, the deposits are characterized by fining-upward sequences with gravel-prone sediments at the bottom, paleosols in the middle and clay-dominated sediments at the top. In the coastal plain the channel deposits consist of coarse to fine sand with sparse basal gravel, while the river banks and thus much of the floodplain consist of clay and silt in the upper coastal plain and very fine sand, silt and clay with salt deposits in the lower coastal plain. Grain-size distribution analysis show a linear decrease in bed load sediments downstream and a concurrent increase in suspended load sediments. These longitudinal variations in the proportions of sediment load are found to be consistent with changes in the topographic slope. The river gradient shows six intervals with downstream decreasing but internally quite constant slopes and sediment composition. Characteristic geomorphologic features are found to be head-cutting channels, chute channels, avulsions and crevasse splays. They illustrate the processes of erosion and deposition in this low-gradient river system and help understand sediment dispersal and their preservation potentials.

Downstream Morphologic Characteristics of the Alluvial Section of Lower River Ogun, Nigeria

Abstract Rivers constitute an important focus of attention in surface water studies because of their dynamic nature. Therefore, natural rivers develop a wide range of channel forms whose characteristics vary as a function of the position within the fluvial systems. This study examined the river channel morphologic parameters along the alluvial section of River Ogun in South western Nigeria. Data on the channel morphologic variables were collected through field measurement of the bankfull cross sectional characteristics of the river from where the longitudinal characteristics were defined. 48 cross sections were randomly established at bankfull stage along the river channel stretch of 90 km. Bankfull depth and width at each of the cross sections were determined using sonar (electronic sounding machine) that was mounted to a boat. Velocity was measured with the aid of a current meter, while other morphological parameters were estimated from the field data. Analysis of variance revealed that downstream morphological characteristics of the river varies distinctively at each cross section with bedslope as the most significantly varied among all other morphologic parameters (F=91.18; P=0.00). Pearson product moment correlation technique revealed that bankfull width had a correlation of 0.8 and 0.9 with wetted perimeter and cross sectional area respectively while bankfull depth (maximum) had correlations of 0.9, 0.8 and 0.78 with hydraulic radius, wetted perimeter and cross sectional area respectively. The research also revealed that gradient affects the discharge with a positive correlation of 0.9. The study ascertains the extent of variability in the morphologic characteristic of River Ogun which provides scientific basis for river maintenance and management.

Bankfull widths in steep headwaters of the flysch Carpathians

Bankfull widths in steep headwaters of the flysch Carpathians

Assessment Approach for Identifying Compatibility of Restoration Projects with Geomorphic and Flooding Processes in Gravel Bed Rivers.

A critical requirement for a successful river restoration project in a dynamic gravel bed river is that it be compatible with natural hydraulic and sediment transport processes operating at the reach scale. The potential for failure is greater at locations where the influence of natural processes is inconsistent with intended project function and performance. We present an approach using practical GIS, hydrologic, hydraulic, and sediment transport analyses to identify locations where specific restoration project types have the greatest likelihood of working as intended because their function and design are matched with flooding and morphologic processes. The key premise is to identify whether a specific river analysis segment (length ~1-10 bankfull widths) within a longer reach is geomorphically active or inactive in the context of vertical and lateral stabilities, and hydrologically active for floodplain connectivity. Analyses involve empirical channel geometry relations, aerial photographic time series, LiDAR data, HEC-RAS hydraulic modeling, and a time-integrated sediment transport budget to evaluate trapping efficiency within each segment. The analysis segments are defined by HEC-RAS model cross sections. The results have been used effectively to identify feasible projects in a variety of alluvial gravel bed river reaches with lengths between 11 and 80 km and 2-year flood magnitudes between ~350 and 1330 m(3)/s. Projects constructed based on the results have all performed as planned. In addition, the results provide key criteria for formulating erosion and flood management plans.

하천 복원을 위한 과거 및 현재 자료 기반의 하천지형학적 특성 분석: 미호천과 내성천을 중심으로

As a basic work for river restoration, analysis on fluvial geomorphological characteristics is made using past and present data to understand close-to-nature geomorphic status. The Miho and the Naesung Rivers are targets of this study. Fluvial geomorphic variables including valley-floor width, sinuosity, bankfull width, channel gradient, bed material size, bankfull discharge and unit stream power are evaluated with dominant processes. Though common sand-bed rivers with similar catchment area, the Miho and the Naesung Rivers are different in terms of valley-floor width, channel shape variables and dominant processes related with longitudinal location. In addition, analyses on interrelationship among the geomorphological variables are carried. Bankfull width is shown to be proportional to bankfull discharge, as is in a rough agreement with the previous studies. Relationship of bankfull discharge and channel gradient shows meandering and braiding are prevalent in the Miho River, whereas the most of the sub-reaches of the Naesung River fall to braiding. Relationship of channel gradient with width-depth ratio indicates dune-ripple processes are dominant in the Miho River, while the Naesung River shows longitudinal diversity from braiding in the downstream sub-reaches to riffle-pool and plane-bed along the upper ones. Analyses based on the past data on a river in a close-to-nature status are thought to be rather reasonable in comparison with those on the same river in a engineered condition.

Development and Evaluation of Bankfull Hydraulic Geometry Relationships for the Physiographic Regions of the United States

AbstractBankfull hydraulic geometry relationships are used to estimate channel dimensions for streamflow simulation models, which require channel geometry data as input parameters. Often, one nationwide curve is used across the entire United States (U.S.) (e.g., in Soil and Water Assessment Tool), even though studies have shown that the use of regional curves can improve the reliability of predictions considerably. In this study, regional regression equations predicting bankfull width, depth, and cross‐sectional area as a function of drainage area are developed for the Physiographic Divisions and Provinces of the U.S. and compared to a nationwide equation. Results show that the regional curves at division level are more reliable than the nationwide curve. Reliability of the curves depends largely on the number of observations per region and how well the sample represents the population. Regional regression equations at province level yield even better results than the division‐level models, but because of small sample sizes, the development of meaningful regression models is not possible in some provinces. Results also show that drainage area is a less reliable predictor of bankfull channel dimensions than bankfull discharge. It is likely that the regional curves can be improved using multiple regression models to incorporate additional explanatory variables.

Instream wood loads in montane forest streams of the Colorado Front Range, USA

Although several studies examine instream wood loads and associated geomorphic effects in streams of subalpine forests in the U.S. Southern Rocky Mountains, little is known of instream wood loads in lower elevation, montane forests of the region. We compare instream wood loads and geomorphic effects between streams draining montane forest stands of differing age (old growth versus younger) and disturbance history (healthy versus infested by mountain pine beetles). We examined forest stand characteristics, instream wood load, channel geometry, pool volume, and sediment storage in 33 pool–riffle or plane-bed stream reaches with objectives of determining whether (i) instream wood and geomorphic effects differed significantly among old-growth, younger, healthy, and beetle-infested forest stands and (ii) wood loads correlated with valley and channel characteristics. Wood loads were standardized to drainage area, stream gradient, reach length, bankfull width, and floodplain area. Streams flowing through old-growth forests had significantly larger wood loads and logjam volumes (pairwise t-tests), as well as logjam frequencies (Kruskal–Wallis test), residual pool volume, and fine sediment storage around wood than streams flowing through younger forests. Wood loads in streams draining beetle-infested forest did not differ significantly from those in healthy forest stands, but best subset regression models indicated that elevation, stand age, and beetle infestation were the best predictors of wood loads in channels and on floodplains, suggesting that beetle infestation is affecting instream wood characteristics. Wood loads are larger than values from subalpine streams in the same region and jams are larger and more closely spaced. We interpret these differences to reflect greater wood piece mobility in subalpine zone streams. Stand age appears to exert the dominant influence on instream wood characteristics within pool–riffle streams in the study area rather than beetle infestation, although this may reflect the relatively recent nature (<10years) of the infestation.

Downstream effects of stream flow diversion on channel characteristics and riparian vegetation in the Colorado Rocky Mountains, USA

AbstractFlow diversions are widespread and numerous throughout the semi‐arid mountains of the western United States. Diversions vary greatly in their structure and ability to divert water, but can alter the magnitude and duration of base and peak flows, depending upon their size and management. Channel geometry and riparian plant communities have adapted to unique hydrologic and geomorphic conditions existing in the areas subject to fluvial processes. We use geomorphic and vegetation data from low‐gradient (≤3%) streams in the Rocky Mountains of north‐central Colorado to assess potential effects of diversion. Data were collected at 37 reaches, including 16 paired upstream and downstream reaches and five unpaired reaches. Channel geometry data were derived from surveys of bankfull channel dimensions and substrate. Vegetation was sampled using a line‐point intercept method along transects oriented perpendicular to the channel, with a total of 100 sampling points per reach. Elevation above and distance from the channel were measured at each vegetation sampling point to analyze differences in lateral and vertical zonation of plant communities between upstream and downstream reaches.Geomorphic data were analyzed using mixed effects models. Bankfull width, depth, and cross‐sectional area decreased downstream from diversions. Vegetation data were analyzed using biological diversity metrics, richness, evenness and diversity, as well as multivariate community analysis. Evenness increased downstream from diversions, through reduced frequency of wetland indicator species and increased frequency of upland indicator species. Probability of occurrence for upland species downstream of a diversion increases at a greater rate beginning around 0·5 m above active channel. The results suggest that channel morphology and riparian plant communities along low‐gradient reaches in montane environments in the Colorado Rocky Mountains are impacted by diversion‐induced flow alteration, with the net effect of simplifying and narrowing the channel and homogenizing and terrestrializing riparian plant communities. Copyright © 2014 John Wiley &amp; Sons, Ltd.

An Evaluation of the Stream Simulation Culvert Design Method in Washington State

AbstractStream simulation has become an increasingly common culvert design method around the world. It is based on the assumption that geologic and hydraulic conditions in natural channels define passage characteristics for migrating fish and that water crossing structures that imitate these conditions can then achieve those same passage characteristics. This study expands on an initial evaluation of 19 culverts in 2003 to 50 culverts and includes methods and analyses comparing hydraulic characteristics based on cross sections, profile variation, and bed texture between each culvert and its paired reference reach situated in an adjacent section of the natural channel of each stream. Taken as a group, these culverts simulate bed texture, 100‐year recurrence interval flood velocity and 2‐year flood width but did not simulate thalweg complexity or other hydraulic metrics. Culvert span, relative to the bankfull width of the stream, does not by itself determine whether the culvert simulates the reference reach. Of the 50 culverts, many of which experienced record floods, only one showed significant bed degradation. Copyright © 2014 John Wiley &amp; Sons, Ltd.

Prediction of recent bank retreat processes at typical sections in the Jingjiang Reach

The Jingjiang Reach is experiencing continuous channel degradation owing to operation of the Three Gorges Project (TGP) in 2003. Significant retreat processes at composite riverbanks have occurred at local sites of this reach, which may influence the stability of the river regime and the effect of existing river training works. Therefore, bank retreat plays a key role in fluvial processes of the Jingjiang Reach, and it is necessary to predict the long-term processes of bank retreat at typical sections in the reach. In this study, various factors influencing bank retreat in the Jingjiang Reach are investigated, based on bank erosion processes at four sections and the corresponding flow and sediment conditions. It is discovered that fluvial erosion intensity is a dominant factor in controlling the processes of bank retreat in the reach, although other factors, such as bank soil properties, can also influence bank retreat. The bankfull width at a section with severe bank erosion since 2002 is defined as being equal to the sum of the bankfull width in 2002 and the accumulated bank retreat distance after 2002. The magnitude of the bankfull width can be expressed as an exponential function of the previous 5-year average fluvial erosion intensity during flood seasons. The accumulated distances of bank retreat at four sections over the period 2002–2012 are predicted using the proposed empirical relationships, with the calculated bank retreat processes agreeing well with observed data.

Channel-reach morphology controls of headwater streams based in flysch geologic structures: An example from the Outer Western Carpathians, Czech Republic

A detailed measurement of 93 channel reaches that were classified with an adjusted Montgomery–Buffington (1997) reach-scale system provided comprehensive information of approximately 9 at-a-reach parameters: the channel gradient, the bankfull width, the bankfull depth, d90, the percentage of resistant rocks in the bed sediment, the number of pieces of large woody debris, valley confinement, direct sediment inputs and the presence of fluvial accumulations in the stream channel. In addition, the quantified intensity of sediment transport (i.e. ratio between sediment supply and transport capacity in longitudinal stream profiles) during flood events has been estimated by the one-dimensional bedload transport model (TOMSED), which was validated in two local streams. The principal component analysis of the at-a-reach parameters did not reveal significant groups of channel-reach morphologies; thus, the selected parameters that exclude sediment transport dynamics within stream longitudinal profile cannot reliably distinguish or predict individual channel reach morphologies. Nevertheless, the channel gradient represented the most significant single explanatory variable for stepped-bed morphologies. The addition of bedload transport parameters demonstrated that limited sediment supply streams and streams with limited transport capacities featured different successions of the channel reach morphologies in terms of the channel gradient and, subsequently, the fluvial continuity. The bedrock-cascades and step-pools were significant for the first case, whereas cascade and step-rapid morphology often occurred in higher sediment supply conditions.

Morphological patterns of headwater streams based in flysch bedrock: Examples from the Outer Western Carpathians

A channel-reach morphology of mountain headwater streams results from the balance between fluvial and hillslope processes. A large dataset (n=102) of channel-reaches coming out from flysch mid-mountains was evaluated and compared with streams occurring in other environments. Two new channel-reach morphologies were distinguished with respect to their differences in the ratio of sediment supply to transport capacity: bedrock-cascades and step-rapids. The chaining of channel-reach morphologies by gradient criteria and the parameters of channel gradient and bankfull width related to the basin area show similar trends, when compared to other recent studies. By contrast, bankfull depth indicates its independence on increasing basin area. Significant downstream coarsening of bed material occurs only in alluvial stepped-bed morphologies (cascades, step-pools and step-rapids) mainly due to strong slope–channel coupling processes. Moreover, the specific ratio of bankfull width and d90 predicts the interlocking of boulder steps as well as well-developed step-pool morphology. An amount of large woody debris is closely related to the activity of strict local forest management in evaluated channels, as a significant decrease in large woody debris pieces is observed with increasing basin area.

Response of reach-scale bankfull channel geometry to the altered flow and sediment regime in the lower Yellow River

The channel geometry in an alluvial river usually varies in response to the altered flow and sediment regime from upstream damming. Bankfull channel dimensions have significantly adjusted along the lower Yellow River (LYR) because of recent channel degradation caused by the operation of the Xiaolangdi Reservoir, which has led to longitudinal variability in the cross-sectional geometry. Therefore, describing the bankfull channel geometry of the LYR is more representative using reach-averaged variables. In this study, an integrated approach is proposed for calculating reach-scale bankfull channel dimensions based on the measured section-scale bankfull geometry in the LYR, and the proposed approach integrates a geometric mean based on the log-transformation with a weighted average based on the spacing between two consecutive sections. The post-flood reach-scale bankfull channel dimensions in three different channel-pattern reaches of the LYR were then calculated annually during the period from 1999 to 2012 using the proposed approach and the surveyed post-flood profiles at 91 sedimentation sections. Calculated results indicate that the channel geometry changed in the components of bankfull width and depth in the braided reach, with an increase of about 400m in the reach-scale bankfull width over the past 13years; and the channel evolution occurred mainly in the component of bankfull depth in the transitional and meandering reaches, with the increased reach-scale bankfull depths of 2.4 and 1.8m, respectively. Finally, empirical relationships for the three reaches of the LYR were developed between the reach-scale bankfull channel dimensions and the previous four-year average discharge and incoming sediment coefficient during flood seasons, with higher correlations being obtained when calibrated by the measurements from 1999 to 2010. Furthermore, the developed relations were also verified by the observed data in 2011 and 2012, with encouraging results.

Continental-scale relationship between bankfull width and drainage area for single-thread alluvial channels

We explore the bankfull width (Wbf) versus drainage area (Ada) relationship across a range of climatic and geologic environments and ask (1) is the relationship between ln(Wbf) and ln(Ada) best described by a linear function and (2) can a reliable relationship be developed for predicting Wbf with Ada as the only independent variable. The principal data set for this study was compiled from regional curve studies and other reports that represent 1018 sites (1 m ≤ Wbf ≤ 110 m and 0.50 km2 ≤ Ada ≤ 22,000 km2) in the continental United States. Two additional data sets were used for validation. After dividing the data into small, medium, and large-size basins which, respectfully, correspond to Ada < 4.95 km2, 4.95 km2 ≤ Ada < 337 km2, and Ada ≥ 337 km2, regression lines from each data set were compared using one-way analysis of covariance (ANCOVA). A second ANCOVA was performed to determine if mean annual precipitation (P) is an extraneous factor in the Wbf versus Ada relationship. The ANCOVA results reveal that using Ada alone does not yield a reliable Wbf versus Ada relationship that is applicable across a wide range of environments and that P is a significant extraneous factor in the relationship. Considering data for very small basins (Ada ≤ 0.49 km2) and very large basins (Ada ≥ 1.0 × 105 km2) we conclude that a two-segment linear model is the most probable form of the ln(Wbf) versus ln(Ada) relationship. This study provides useful information for building complex multivariate models for predicting Wbf.

Direct and indirect drivers of instream wood in the interior Pacific Northwest, USA: decoupling climate, vegetation, disturbance, and geomorphic setting

AbstractInstream wood is a driver of geomorphic change in low-order streams, frequently altering morphodynamic processes. Instream wood is a frequently measured component of streams, yet it is a complex metric, responding to ecological and geomorphic forcings at a variety of scales. Here we seek to disentangle the relative importance of physical and biological processes that drive wood growth and delivery to streams across broad spatial extents. In so doing, we ask two primary questions: (1) is riparian vegetation a composite variable that captures the indirect effects of climate and disturbance on instream wood dynamics? (2) What are the direct and indirect relationships between geomorphic setting, vegetation, climate, disturbance, and instream wood dynamics? We measured riparian vegetation composition and wood frequency and volume at 720 headwater reaches within the American interior Pacific Northwest. We used ordination to identify relationships between vegetation and environmental attributes, and subsequently built a structural equation model to identify how climate and disturbance directly affect vegetation composition and how vegetation and geomorphic setting directly affect instream wood volume and frequency. We found that large wood volume and frequency are directly driven by vegetation composition and positively correlated to wildfire, elevation, stream gradient, and channel bankfull width. Indicator species at reaches with high volumes of wood were generally long-lived, conifer trees that persist for extended durations once delivered to stream habitats. Wood dynamics were also indirectly mediated by factors that shape vegetation: wildfire, precipitation, elevation, and temperature. We conclude that wood volume and frequency are driven by multiple interrelated climatic, geomorphic, and ecological variables. Vegetation composition and geomorphic setting directly mediate indirect relationships between landscape environmental processes and instream large wood. Where climate or geomorphic setting preclude tree establishment, reaches may remain naturally depauperate of instream wood unless wood is transported from elsewhere in the stream network.

Conceptual Model for Simulating the Adjustments of Bankfull Characteristics in the Lower Yellow River, China

We present a conceptual model for simulating the temporal adjustments in the banks of the Lower Yellow River (LYR). Basic conservation equations for mass, friction, and sediment transport capacity and the Exner equation were adopted to simulate the hydrodynamics underlying fluvial processes. The relationship between changing rates in bankfull width and depth, derived from quasiuniversal hydraulic geometries, was used as a closure for the hydrodynamic equations. On inputting the daily flow discharge and sediment load, the conceptual model successfully simulated the 30-year adjustments in the bankfull geometries of typical reaches of the LYR. The square of the correlating coefficient reached 0.74 for Huayuankou Station in the multiple-thread reach and exceeded 0.90 for Lijin Station in the meandering reach. This proposed model allows multiple dependent variables and the input of daily hydrological data for long-term simulations. This links the hydrodynamic and geomorphic processes in a fluvial river and has potential applicability to fluvial rivers undergoing significant adjustments.

Predicting presence and absence of trout (Salmo trutta) in Iran

Species distribution modelling, as a central issue in freshwater ecology, is an important tool for conservation and management of aquatic ecosystems. The brown trout (Salmo trutta) is a sensitive species which reacts to habitat changes induced by human impacts. Therefore, the identification of suitable habitats is essential. This study explores the potential distribution of brown trout by a species distribution modelling approach for Iran. Furthermore, modelling results are compared to the distribution described in the literature. Areas outside the currently known distribution which may offer potential habitats for brown trout are identified. The species distribution modelling was based on five different modelling techniques: Generalised Linear Model, Generalised Additive Model, Generalised Boosting Model, Classification Tree Analysis and Random Forests, which are finally summarised in an ensemble forecasting approach. We considered four environmental descriptors at the local scale (slope, bankfull width, wetted width, and elevation) and three climatic parameters (mean air temperature, range of air temperature and annual precipitation) which were extracted on three different spatial extents (1/5/10 km). The performance of all models was excellent (≥0.8) according to the TSS (True Skill Statistic) criterion. Slope, mean and range of air temperature were the most important variables in predicting brown trout occurrence. Presented results deepen the knowledge about distribution patterns of brown trout in Iran. Moreover, this study gives a basic background for the future development of assessment methods for riverine ecosystems in Iran.

Scale-dependent interactions between wood and channel dynamics: Modeling jam formation and sediment storage in gravel-bed streams

[1] A stochastic model is used to investigate how the geomorphic function of wood changes with watershed scale, assuming wood recruitment occurs due to the mortality of individual trees, not to mass recruitment events such as landslides or episodic bank erosion. The model replicates the downstream decline in total wood load observed in the field, but predicts that the functional wood load peaks in channels having bankfull widths about 33% of the characteristic riparian tree height. The model also predicts that the greatest potential impact of jams on channel pattern—both in terms of sediment stored behind individual jams and the potential for jams to trigger avulsions—will typically be associated with channel widths between 25% and 67% of the riparian tree height. The simulation results are used to refine the categories that describe wood in alluvial channels, and the equivalent terms that describe the size of streams with forested riparian areas: small channels (or channels with large wood) are associated with widths less than 25% of the tree height; large channels (or channels with small wood) are associated with widths greater than 67% of tree height; and medium channels (or channels with intermediate wood) have widths between 25% and 67% of the tree height. We surmise that large wood acts primarily to store bed material (in small channels); intermediate wood tends to form channel-spanning jams, which can induce channel avulsions and create anabranched channel patterns (in medium channels); and small wood may increase the morphologic diversity, but does not store significant quantities of bed material or form channel-spanning jams capable of inducing stream avulsions (in large channels).

Assessing stream bank condition using airborne LiDAR and high spatial resolution image data in temperate semirural areas in Victoria, Australia

Stream bank condition is an important physical form indicator for streams related to the environmental condition of riparian corridors. This research developed and applied an approach for mapping bank condition from airborne light detection and ranging (LiDAR) and high-spatial resolution optical image data in a temperate forest/woodland/urban environment. Field observations of bank condition were related to LiDAR and optical image-derived variables, including bank slope, plant projective cover, bank-full width, valley confinement, bank height, bank top crenulation, and ground vegetation cover. Image-based variables, showing correlation with the field measurements of stream bank condition, were used as input to a cumulative logistic regression model to estimate and map bank condition. The highest correlation was achieved between field-assessed bank condition and image-derived average bank slope (R2 1/4 0.60, n 1/4 41), ground vegetation cover (R2 1/4 0.43, n 1/4 41), bank width/height ratio (R2 1/4 0.41, n 1/4 41), and valley confinement (producer's accuracy 1/4 100%, n 1/4 9). Crossvalidation showed an average misclassification error of 0.95 from an ordinal scale from 0 to 4 using the developed model. This approach was developed to support the remotely sensed mapping of stream bank condition for 26,000 km of streams in Victoria, Australia, from 2010 to 2012.

Benthic macroinvertebrate communities in five rivers of the Coastal Hudson Bay Lowland

As a precursor to developing a biomonitoring program for rivers of the Coastal Hudson Bay Lowland, this study characterized and compared the benthic macroinvertebrate communities and water chemistry in 5 remote, previously undescribed, rivers near Fort Severn, Ontario, Canada. The pH of river water ranged from 8.1 to 8.7, total phosphorus from 11 to 26 μg L−1, dissolved organic carbon from 8 to 12 mg L−1, and chloride from 56 to 153 mg L−1. A total of 57 benthic macroinvertebrate taxa were represented, and the 10 most numerically dominant were the Chironominae (26 % of collected individuals), Orthocladiinae (16 %), oligochaetous clitellata (9 %), Hyalellidae (7 %), Hydropsychidae (6 %), Gammaridae (5 %), Elmidae (5 %), Sphaeriidae/Pisidiidae (4 %), Nemata (3 %), and Tanypodinae (3 %). Rivers’ positions in ordinations of chemical and biological datasets were similar, suggesting that water chemistry has a role in structuring riverine benthic communities in the study region. Correlations between water-chemistry or habitat predictors and site-scores in the ordination of benthic macroinvertebrate taxa counts suggested that biological community structure was most associated with river-water pH, nutrient concentrations (e.g., total phosphorus, nitrogenous compounds, dissolved organic carbon, calcium, and silicate), the relative abundance of submerged macrophytes, conductivity (i.e., the concentrations of chloride and various other dissolved ions), and several geomorphological variables (e.g., bank-full river width, current speed, and the size of the dominant inorganic particles in the pavement layer of the streambed). Interest in mineral extraction and other resource-based exploration in Ontario’s Far North is increasing. This study represents a start on baseline characterization for ecological monitoring and cumulative effects assessment that should proceed along with northern development.