Next Article in Journal
Characteristics of Stormwater Quality in Singapore Catchments in 9 Different Types of Land Use
Next Article in Special Issue
Climatic Trends of West Virginia: A Representative Appalachian Microcosm
Previous Article in Journal
Methane Emissions Driven by Adding a Gradient of Ethanol as Carbon Source in Integrated Vertical-Flow Constructed Wetlands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterizing Land Use Impacts on Channel Geomorphology and Streambed Sedimentological Characteristics

by
Sean J. Zeiger
1,* and
Jason A. Hubbart
2,3
1
School of Natural Resources, University of Missouri, 203-T ABNR Building, Columbia, MO 65211, USA
2
Institute of Water Security and Science, West Virginia University, 3109 Agricultural Sciences Building, Morgantown, WV 26506, USA
3
Davis College of Agriculture, Natural Resources and Design, Schools of Agriculture and Food, and Natural Resources, West Virginia University, 3109 Agricultural Sciences Building, Morgantown, WV 26506, USA
*
Author to whom correspondence should be addressed.
Water 2019, 11(5), 1088; https://doi.org/10.3390/w11051088
Submission received: 3 April 2019 / Revised: 10 May 2019 / Accepted: 12 May 2019 / Published: 24 May 2019

Abstract

:
Land use can radically degrade stream physical habitat via alterations to channel geomorphology and sedimentological characteristics. However, independent and combined influences such as those of agricultural and urban land use practices on channel geomorphology and substrate composition remain poorly understood. To further understanding of mixed land use influence on stream physical habitat, an intensive, 56 km hydrogeomorphological assessment was undertaken in a representative mixed land use watershed located in Midwestern USA. Sub-objectives included quantitative characterization of (1) channel geomorphology, (2) substrate frequency and embeddedness, and (3) relationships between land use, channel geomorphology, and substrate frequency and embeddedness. Channel geomorphology, and stream substrate data were directly measured at survey transects (n = 561) every 100 m of the entire 56 km distance of the reference stream. Observed data were averaged within five sub-basins (Sites #1 to #5) nested across an agricultural-urban land use gradient. Multiple regression results showed agricultural and urban land use explained nearly all of the variance in average width to depth ratios (R2 = 0.960; p = 0.020; n = 5), and maximum bank angle (R2 = 0.896; p = 0.052; n = 5). Streambed substrate samples of pools indicated significantly (p < 0.001) increased substrate embeddedness at agricultural Site #1 (80%) located in the headwaters and urban Site #5 (79%) located in the lower reaches compared to rural-urban Sites #2 to #4 (39 to 57%) located in the mid-reaches of the study stream. Streambed substrate embeddedness samples of riffles that ranged from 51 to 72% at Sites #1 and #5, and 27 to 46% at Sites #2 to #4 were significantly different between sites (p = 0.013). Percent embeddedness increased with downstream distance by 5% km−1 with the lower urban reaches indicating symptoms of urban stream syndrome linked to degraded riffle habitat. Collectively, observed alterations to channel morphology and substrate composition point to land use alterations to channel geomorphology metrics correlated with increased substrate embeddedness outside of mid-reaches where bedrock channel constraints accounted for less than 3% of substrate frequency. Results from this study show how a hydrogeomorphological assessment can help elucidate casual factors, target critical source areas, and thus, guide regional stream restoration efforts of mixed-land-use watersheds.

Graphical Abstract

1. Introduction

The importance of channel geomorphology and stream substrate composition with regards to stream habitat has been well-documented [1,2,3,4,5]. Previous studies indicated that channel geomorphology and stream substrate are important aspects of stream physical habitat [6,7,8]. For example, substrate is a medium for settlement by propagules [9]. Successful propagules that mature into aquatic plants provide substrate stability, food, dissolved oxygen, physical habitat, and refugia beneficial for healthy stream ecosystems [7]. There is therefore also a direct link between suitable stream substrate composition and the abundance and distribution of macroinvertebrates [8,10], mussels [6], and fishes [11]. Given the influence on stream health, it is important to understand the dominant physical processes that control channel geomorphology and substrate composition.
Process-based understanding indicates channel geomorphology and stream substrate composition is dependent on sediment supply, and transport [12]. Surface runoff and streamflow are largely the erosive forces and transport mechanisms for overland soil erosion and stream bed and bank erosion, respectively [12]. However, a myriad of independent and interacting natural and anthropogenic factors influence surface runoff, streamflow, sediment supply and transport (e.g., meteorological conditions, hydrology, topography, soils and underlying geology, land cover, and human activities) [12]. Thus, regions associated with physiographic characteristics known to increase sediment supply and transport are areas of concern susceptible to alterations of channel geomorphology and stream substrate composition, and therefore physical habitat.
Streams of the Midwestern USA are commonly associated with increased sediment supply and degraded stream health [13]. For example, Gellis et al. [13] found channel sources of sediment accounted for the majority (>50%) of bed sediment in 79% of 99 Midwestern watersheds sampled. Claypan soils overlain by loess are especially subject to increased surface runoff and sediment supply to nearby streams [14,15]. A study by Willett et al. [15] indicated bank sediment accounted for 79% to 96% of the total in-stream sediment in two claypan watersheds of the central USA. A study by Huang [16] in mid-Missouri, showed bank erosion contributed approximately 67% of in-stream suspended sediment load in Hinkson Creek Watershed located in central claypan region of the USA. Given previous findings, it is not surprising that Midwestern streams are commonly characterized by increased frequency of fine bed sediment textures and increased substrate embeddedness related to degraded algal, macroinvertebrate, and fish communities [17,18,19]. Increased substrate embeddedness can influence stream community composition by reducing available riffle-pool habitat when increased deposition of fine sediments and sand fill interstitial spaces of course gravel substrate [20,21]. In Midwestern streams and other regions associated with increased sediment supply, human activity (i.e., land use) can exacerbate problems and further reduce suitable stream physical habitat and aquatic refugia.
Previous studies showed agricultural land use activities can increase sediment supply potentially leading to problems with sedimentation and substrate embeddedness [13,14,22,23]. For example, conventional tillage operations and livestock presence have been shown to increase erosion rates and sediment supply to nearby stream networks [12]. Given stream access to riparian areas, livestock can trample biomass and erode stream banks thereby increasing sediment supply [22]. Near channel removal of streambank stabilizing riparian vegetation has also been shown to increase bank erosion and alter channel geomorphology in agricultural watersheds [23].
Urban land use has also been shown to increase stream sediment supply. A literature review by Walsh et al. [24] presented consistent confirmation of urban land use influence on watershed hydrology, water chemistry, channel geomorphology, organic matter, algae, macroinvertebrates, fishes, and ecosystem processes. The term ‘urban stream syndrome’ is widely accepted to indicate consistent urban land use influence on stream health [24,25]. More specifically, previous studies indicated urban land use can increase sediment supply and pollutant transport when impervious surfaces and engineered water infrastructure increase the volume and velocity of surface runoff, stream power, and transport capacity [24]. Increased stream power associated with large storms can alter channel geomorphology via channel widening, channel incision, bank erosion and mass wasting [24]. For example, Gellis et al. [13] showed Tropical Storm Lee (a recent 100-year storm) caused stream bank erosion that accounted for 70% of suspended sediment concentration in the suburban and urban Upper Difficult Run, Virginia, USA. Results from Gellis et al. [13] emphasized the importance of streambank sources of sediment supply in suburban and urban watersheds. Thus, results from previous studies indicate streambank stabilization efforts are often necessary in urbanizing watersheds [13,14,24].
The combined influence of agricultural and urban land use practices on channel geomorphology and substrate composition is not well-understood. Investigations that quantitatively characterize land use influence on channel geomorphology and stream substrate composition in representative watersheds could help guide regional stream restoration efforts [4,26]. Therefore, the overall objective of this study was to quantify general trends and relationships between land use, stream hydromorphology, and substrate composition using observed data collected during a physical habitat assessment (PHA) in a representative experimental watershed located in the Midwestern USA. Sub-objectives were to (1) present observed spatial variability of channel geomorphology metrics, substrate frequency, and substrate embeddedness, and (2) quantify relationships between land use, channel geomorphology, and substrate frequency and embeddedness. Additionally, results are discussed with implications for watershed management in similar watersheds where agricultural and urban land uses have exacerbated problems with channel geomorphology alteration and streambed sedimentation.

2. Methods

2.1. Study Site

Hinkson Creek Watershed (HCW), is a rapidly urbanizing mixed-land use (forest, agriculture, urban) watershed located in the Lower Missouri Mississippi River Basin (LMMRB) (Table 1, Figure 1). Drainage area of HCW is approximately 230 km2. Elevation ranges from 274 m above mean sea level (AMSL) in the headwaters to 177 m AMSL near the watershed outlet. At the time of this study, agricultural land use accounted for 36.7% HCW (28.4% grazing pasture and 8.3% row crop). Urban and suburban areas accounted for 29.0% of the total land use. Approximately 28 km2 of impervious surfaces associated with urban land use were located within the municipal boundary of The City of Columbia (population 121,717; USCB [27]).
Soils in the upper elevations of HCW are generally loamy loess with a well-developed underling claypan in the argillic horizon of smectitic mineralization corresponding to the Mexico-Leonard association [28]. Claypan soils have increased surface runoff potential [29]. Soils that consist of a cherty clay solution residuum corresponding to the Weller-Bardley-Clinkenbeard (CBC) association are found in the lower reaches (upland of the alluvium). Soils in HCW are underlain by Burlington formation Mississippian series limestone in the lower reaches [28,29,30].
Climate in HCW is dominated by continental polar air masses in the winter and maritime and continental tropical air masses in the summer [31]. A 20-year climate record (2000–2019) obtained from the Sanborn Field climate station located on University of Missouri campus showed mean annual total precipitation was 962 mm and mean annual air temperature was 13.5 °C. A wet season occurs primarily during March through June [31,32]. The HCW main drainage, Hinkson Creek, is primarily stormflow dominated with a base flow index (ratio of base flow to total streamflow) of approximately 0.25 calculated using daily discharge data collected at a U.S. Geological Survey (USGS) gaging station (USGS 06910230, Site #4 in Figure 1) that has been intermediately monitoring stage since 1967.
During winter 2008, HCW was instrumented with a nested-scale experimental study design to investigate land use alterations to water quantity and quality [31]. Gauging sites (n = 5) partitioned the catchment into five sub-basins, each with different dominant land uses [31]. Site #5 was located near the watershed outlet, and Sites #1 to #4 were nested within (Figure 1). Site #1 (57.0% agricultural, 35.9% forested, 4.7% urban) was located in the agricultural land use dominated headwaters. Moving downstream from agricultural Site #1, agricultural land use decreased as urban land use increased. Site #3 (49.5% agricultural; 35.4% forested; 13.0% urban) was located at the rural-urban interface of HCW. Site #5 (39.4% agricultural; 33.1% forested; 26.5% urban) is located in the lower most urbanized reaches. Results from previous studies showed land use alterations to the flow regime [28,33,34,35], environmental flows [36], stream temperature [37,38], suspended sediment [14,39,40], nutrients (i.e., total inorganic nitrogen species, total phosphorus) [41], and chloride [42] in HCW.

2.2. PHA Data Collection and Analysis

A physical habitat assessment (PHA) was performed in HCW during the study period (2013–2014) [43,44,45]. During a single survey, channel geomorphology and stream substrate data were collected at survey points (n = 561) every 100 m along the entire 56 km length of Hinkson Creek. At each survey point, data were collected (using a clinometer, extension pole, laser level, a laser range finder and/or meter stick) at principal transects that spanned from stream bank to bank perpendicular to flow. Additionally, data were collected at transects located 5 m upstream and downstream of the principal transect. Upstream and downstream transects were parallel to the principal transect. At any confluence of the main channel, three transects (upstream, downstream, and upstream in the tributary) were located equidistant from the center of the confluence.

2.2.1. Observed Channel Geomorphology

Channel width, wetted width of the stream, bankfull width, bank angle, bank height, and channel depth were measured using a laser level and/or laser range finder at each principal transect (n = 561). Following methods suggested by Harrelson et al. [46], physical indicators of bankfull level included the top of pointbars, changes in vegetation from aquatic to terrestrial, changes in slope, changes in bank material (e.g., from coarse gravel to sand), bank undercuts, or stain lines on bedrock or boulders. Bankfull width was measured parallel of the stream surface and perpendicular to stream flow from the lowest bank (i.e., bankfull bank) to the opposite bank.
Observed channel geomorphology data (collected every 100 m) dependent on drainage area were normalized as per methods used by Yanites and Tucker [47] using the following equation:
C h * = C h m k D A b
where Ch* is area-normalized channel geomorphology, Chm is the measured channel morphology metric, DA is drainage area, and k and b are fitted power regression coefficients.
Channel geomorphology data were reduced by averaging within five nested HCW sub-basins (Figure 1) to quantify the change in channel geomorphology with downstream distance across a rural-urban land use gradient. Sub-basins were delineated between nested gauging sites, and thus, were not cumulative watersheds. For example, channel width at Site #2 reflected the average of all the channel width measurements collected every 100 m upstream of Site #2 and downstream of Site #1. One-way Analysis of Variance (ANOVA) and Tukey Kramer post-hoc multiple comparison tests were used to test for significant differences (CI = 95%, p < 0.05) in average channel geomorphology metrics (i.e., channel width, bank height, bank angle, etc.) between five HCW sub-basins [48]. Tukey Kramer post-hoc multiple comparison test was selected to elucidate site differences with a narrow confidence interval [49].

2.2.2. Observed Streambed Substrate Frequency and Embeddedness

Substrate particle size-class (Table 2) was estimated at each principal, upstream, downstream transect following methods suggested by Peck et al. [50], and Wolman [51]. At each transect, particles were sampled in five locations from bank to bank (left bank, left center, center, right center, and right bank) perpendicular to flow for a total of 8415 stream substrate samples. Substrate frequency was quantified as a percent of each individual substrate type relative to the total number of substrate samples collected within the area of interest. Additionally, substrate embeddedness of each particle was estimated as the percent vertical entrainment of a streambed substrate sample. Channel unit type (i.e., trench pool, plunge pool, impoundment pool, pool, split channel, riffle, glide, dry channel, etc.) was also recorded at each principal transect as per methods used by Peck et al. [50]. While the influence of water depth on channel unit identification was not quantified in the current study, adverse effects were assumed negligible considering the water depth recorded in the thalweg ranged from 0 to 280 cm with a median of 40 cm at the time of channel unit identification. Channel units were grouped into general categories of pool, riffle, and glide for further analysis. Observed substrate particle size frequency and percent embeddedness data were reduced by averaging in channel units (i.e., pools, riffles, glides, and dry channels), from bank to bank, and within each HCW sub-basin to quantify the change in stream substrate composition associated with riffle-pool habitat with downstream distance across a rural-urban land use gradient. Sub-basins were delineated between nested gauging sites and were thus not cumulative.

3. Results

3.1. Observed Geomorphology

Results from power regression analysis showed drainage area explained 41.2% of the variance in channel width, 58.6% of the variance in bankfull width, 29.9% of the variance in bank height, and 42.1% of the variance in bankfull depth (p < 0.001; n = 561) (Figure 2). After channel width and height and bankfull width and depth metrics were area-normalized, results showed significant differences in channel morphology metrics averaged within five sub-basins of HCW (Table 3). For example, area-normalized average channel width ranged from 0.12 to 0.19 m km−2. Area-normalized average bankfull width ranged from 0.21 to 0.27 m km−2. Area-normalized average bank height ranged from 0.38 to 0.56 m km−2. Area-normalized average bankfull depth ranged from 0.34 to 0.49 m km−2. Average minimum and maximum bank angle ranged from 28.0 to 47.8 degrees. Average width to depth ratios and bed slope ranged from 6.3 to 10.3, and 0.14 to 0.22%, respectively. These results quantitatively characterized channel geomorphology metrics important for stream restoration efforts in HCW (Table 3).
No significant relationships were observed between any one land use and land cover (LULC) index shown in Table 1 and average channel morphology metrics at five sub-basins in HCW (R2 ≤ 0.593; p ≥ 0.080; n = 5). However, multiple linear regression (MLR) results showed agricultural and urban land use combined explained nearly all of the variance in average width to depth ratios (R2 = 0.960; p = 0.020; n = 5), and maximum bank angle (R2 = 0.896; p = 0.052; n = 5). Thus, results implied the expected relationships between drainage area and channel morphology were altered, at least in part, by the combined influence of agricultural and urban land use in HCW.

3.2. Observed Streambed Substrate Frequency and Embeddedness

Results from stream substrate frequency analysis quantitatively characterized bank to bank variability in streambed substrate composition important for stream physical habitat and aquatic refugia in HCW. Results showed percent fine substrate (FN) summarized for each bank to bank position along the entire drainage ranged from 11.6% at the center channel to 62.0% on the right descending bank. Sand substrate (SA) ranged from 17.3% on the left bank to 33.9% center channel. Percent fine gravel (GF), course gravel (GC), cobble (CB), and small boulder (SB) followed a similar trend as sand substrate from bank to bank. For example, GC ranged from 3.3% at the right bank to 21.3% at the center channel position. There was no obvious bank to bank trend in substrate frequency large boulders (XB), rough bedrock (RR), smooth bedrock (RS), or the rest of the substrate types presented in Table 4.
Examination of substrate frequency at five sub-basins showed the frequency of substrate types smaller than GC decreased with downstream distance from agricultural Site #1 in the headwaters to urban Site #3 at the rural-urban interface of the watershed. Continuing downstream, substrate frequency of substrate types smaller than GC increased from Site #3 to urban Site #5 located near the watershed outlet. For example, FN, SA, and GF decreased from Site #1 to Site #3 by 17%, 19%, and 3%, respectively. Continuing downstream, FN, SA, and GF increased from Site #3 to Site #5 by 19%, 12%, and 3%, respectively. Conversely, some of the larger substrate types (i.e., GC, CB, SB, and XB) increased with downstream distance from Site #1 to Site #3, and then, decreased from Site #3 to Site #5. Observed substrate composition data showed increased bed rock channel constraints at Sites #2 to #4 (9 to 11% RR and RS) located in the mid-reaches relative to Site #1 (3% RR) located in the agricultural headwaters and Site #5 (1% RR and RS) in the lower urban reaches (Table 5). These results quantifiably characterized the change of stream substrate frequency across an agricultural-urban land use gradient, and in synthesis, indicated increased substrate embeddedness in the agricultural headwaters and lower urban reaches located outside of the observed bedrock channel constraints in the mid-reaches of the study catchment.
Substrate embeddedness also varied between different stream habitats (i.e., riffle, glide, and pool). For example, results showed that pools (74.9% embeddedness) were associated with significantly greater substrate embeddedness than riffles (58.5% embeddedness) and glides (60.1% embeddedness) (p = 0.0014; CI = 95%). However, no significant differences were observed between riffles and glides (p = 0.9367; CI = 95%). Significant differences in substrate embeddedness were especially apparent at Sites #1 and #5 compared to the mid-reaches (i.e., Sites #2 to #4) (p = 0.013; CI = 95%) (Table 6). Stream banks were embedded with FN and SA in more than 88% of pools habitat at Sites #1 and #5. Even greater were differences in substrate embeddedness at center of channel samples (i.e., left center, center, and right center) of pools at Sites #1 and #5 (70 to 79%) compared to Sites #2 to #4 (39 to 57%). Center of channel samples of riffles ranged from 51 to 72% at Sites #1 and #5, and 27 to 46% at Sites #2 to #4. Due to the observed widespread substrate embeddedness, results generally showed reduced spatial complexity of substrate composition in riffles, glides, and pools of HCW.
Substrate embeddedness was significantly correlated at the 0.05 level with 4 of 8 channel morphology metrics considered in this study. For example, minimum bank angle (R2 = 0.8511; p = 0.0164; n = 5), maximum bank angle (R2 = 0.8400; p = 0.0183; n = 5), width to depth ratio (R2 = 0.7225; p = 0.0424; n = 5), and bed slope (R2 = 0.7424; p = 0.0384; n = 5) explained a substantial amount of variance in substrate embeddedness between sites in HCW. Minimum and maximum bank angle were positive correlates. Width to depth ratios and bed slope were negative correlates. Relationships between drainage area and land use were not significant at a 0.05 level. However, there was a general trend for increased substrate embeddedness in the agricultural headwaters at Site #1 and in the lower urban reaches at Site #5 (80% embedded) relative to the mid-reaches at Sites #2 to #4 (56 to 60% embedded) of HCW.
There was also a general trend for increased frequency of pool habitat coupled to decreased frequency of riffle and glide habitats in the agricultural headwaters and in the lower urban reaches of HCW (Table 7). The frequency of pools ranged from 56.1% at the urban Site #4 of HCW to 85.7% at urban Site #5. Conversely, the frequency of riffle habitat increased from 12.2% at Site #5 in the lower urban reaches to 28.0% at urban Site #4. The frequencies of pool, riffle, and glide habitats were not significantly correlated with substrate embeddedness (R2 ≥ 0.420; p ≤ 0.140; n = 5). However, results showed strong correlations between channel unit frequency and frequency of fine substrate type (0.927 ≤ R2 ≥ 0.854; p ≤ 0.016; n = 5) highlighting observed streambed siltation in Hinkson Creek.

4. Discussion

4.1. Observed Geomorphology

Results from this study show how expected longitudinal change in channel geomorphology can be altered independently or by combined bedrock channel constraints and human activity. For example, the mid-reaches of HCW were associated with increased bedrock channel constraints (Figure 3), which according to Montgomery et al. [52] indicates transport capacity in excess of sediment supply. Outside of the mid-reaches (i.e., Sites #1 and #5) where less than 3% presence of channel bedrock constraints were observed, signs of channel incision were apparent. Site #1 in the agricultural headwaters and Site #5 located the lower urban reaches of HCW (Figure 3). Agricultural and urban land use explained nearly all the variance in average width to depth ratios, and maximum bank angle (0.896 ≤ R2 ≥ 0.960; 0.020 ≤ p ≥ 0.052; n = 5). At Site #1 located in the agricultural headwaters (drainage area = 79.0 km2), width to depth ratios increased at a rate of about 0.3 km−1, and maximum bank angle decreased by about 0.5 degrees km−1 as agricultural land use decreased from 100 to 56.9%, forested land use increased from 0 to 35.9%, and urban land use increased from 0 to 4.7% over 22.8 km stream distance (Figure 4). With stream distance from agricultural Site #1 to sub-urban Site #3 located at the rural-urban interface of HCW, (drainage area = 116.2.0 km2), width to depth ratios increased by a rate of about 0.01 m m−1 km−1, and maximum bank angle decreased by about 0.1 degrees km−1 as agricultural land use decreased to 49.5%, forested land use decreased by 0.5%, and urban land use increased to 13.0% over 12.6 km stream distance (Figure 4). Continuing downstream from sub-urban Site #3 to urban Site #5 located in the lower urban reaches (drainage area = 207.5 km2), width to depth ratios decreased by about 0.2 km−1, and maximum bank angle increased by about 0.2 degrees km−1 as agricultural land use decreased to 38.4%, forested land use decreased to 33.1%, and urban land use increased to 26.5% over 17.6 km stream distance (Figure 4). Bankfull depth increased by 0.13 m km−1 at Site #1, 0.02 m km−1 between Sites #1 and #3, and 0.15 m km−1 between Sites #3 and #5. In combination, these results indicated increased channel incision at Site #1 in the agricultural headwaters where agricultural land use accounted for greater than 50% of total catchment area, and Site #5 the lower urban reaches where urban land use accounted for greater than 20% of total catchment area in HCW. Thus, Sites #1 and #5 were considered hot spots of channel incision due to increased rates of change in channel geomorphology metrics.
Previous studies have shown agricultural land use can result in channel incision through various mechanisms including, but not limited to, channelization, deforestation, and alterations to soil and channel hydraulics [53,54,55,56]. Simon et al. [56] noted increased pore water pressure coupled to decreased sheer strength can lead to streambank erosion. Given the influence of soil hydraulic forces on streambank stability, agricultural areas are subject to increased bank failure where crop irrigation lowers water table levels in the vicinity of a stream. A study by Zaimes et al. [54] showed incised stream channels were associated with increased bank mass wasting, streambank erosion and sediment load in Beak Creek, an agricultural land use dominated Midwestern stream located in central Iowa, USA. Streambank erosion rates varied from meandering row crop fields (387 mm year−1), cattle and horse pastures (295 mm year−1), and meandering forest buffered reaches (142 mm year−1). Results indicated that forested riparian buffers would reduce streambank erosion by 72%. Streambank stabilizing root systems associated with riparian vegetation add roughness, reduce stream power, and enhance bank accretion [56]. Results from the aforementioned studies are in agreement with results from the current work showing agricultural land use alterations to channel morphology of Midwestern streams.
Previous studies have also shown urban land use can cause channel incision via alterations to streamflow regimes [24,25,57,58,59]. A literature review by Walsh et al. [24] showed increased impervious surfaces associated with urban land use can cause increased volume and velocity of surface runoff and a flashy hydrologic streamflow response linked to increased bank wasting, channel incision, and scouring. For example, Jordan et al. [60] showed urban land use alterations to flow caused a 9 to 61% increase of sediment yield due to channel incision and bank erosion in Berryessa Creek, California, USA. A study by Shields et al. [59] quantified differences in channel incision, streamflow, water quality, and stream physical habitat between rural and urban catchments located in the Yazoo River basin, Mississippi, USA. Results showed urban land use was associated with decreased physical aquatic habitat, 6.4 times median rate of rise, 1.8 times channel depth, 3.5 times channel width, 2 to 3 times turbidity and suspended solids, 2 times fish species, and 4 times the amount of fish biomass per unit of effort [59]. Results from the current work are a novel addition to previous studies considering the intensive sampling regimen (n = 561) that made possible the estimation of the rate of change in channel morphology and substrate composition across an agricultural-urban land use gradient.

4.2. Observed Streambed Substrate Frequency and Embeddedness

No channel morphology variable significantly explained the variance in substrate embeddedness of pools, riffles, and glides for each bank to bank sampling position (Figure 5) highlighting the spatial complexity of streambed substrate composition in this study. Spatial complexity of substrate composition was, in part, shaped by the thalweg which meandered from bank to bank. The thalweg was generally associated with decreased substrate embeddedness due to increased stream velocity, and thus, increased sediment transport capacity. While much of the observed bank to bank variability in substrate embeddedness was attributed to thalweg position, results from the current work in combination with previous research in the region indicated longitudinal variability in substrate composition was attributed to the presence claypan soils and agricultural land use in the headwaters, increased bed slope and bedrock channel constraints in the mid-reaches, and the influence of urban land use associated with increased impervious surfaces in the lower reaches (Figure 3).
Previous studies in HCW and elsewhere have shown claypan soils and agricultural land use are associated with increased surface runoff, soil erosion, and channel sediment supply [14,15,29]. Lerch et al [29] noted claypan soils corresponding to the Mexico-Leonard association consisting of an argillic soil horizon of smectitic mineralogy with clay content of 450–650 g kg−1 formed at 10 to 50 cm depth are characterized by increased surface runoff. Willett et al. [15] showed claypan soils were associated with increased bank sediment supply that accounted for 88% of total sediment supply. Streambank erosion was particularly high during winter months attributed to a combination of increased frequency and large magnitude flow events, freeze/thaw cycles, high antecedent moisture conditions, and lack of vegetation. In the current work, increased bank angle and substrate embeddedness were apparent during field sampling in the agricultural headwaters of HCW where claypan soils are present. Results from the current work are among the first to quantitatively characterize agricultural land use influence on substrate frequency and substrate embeddedness in the Central Claypan Region and point to a need to mitigate the influence claypan soils and agricultural land use on degradation of stream hydrogeomorphology in HCW. Thus, results from this study in combination with previous research have implications in other agricultural watersheds where near surface soil features (e.g., claypans, argillic horizons, or fragipans) have increased sediment supply.
While substrate embeddedness was greater in the agricultural headwaters, substrate embeddedness decreased as bed slope increased in the mid-reaches. Agricultural Site #1 was associated with 26 to 39% greater frequency of substrate smaller than GC, and 20 to 25% greater substrate embeddedness compared to Sites #2 to #4 in the mid-reaches of HCW (Figure 6). Bed slope, width to depth ratio, and the frequency of substrate greater than GC were inversely related to substrate embeddedness (Figure 4). Process-based understanding of the control of channel morphology on streamflow and sediment transport [52] indicated the aforementioned general trends in the observed data were physically meaningful. Given that bed slope controls velocity of streamflow, stream capacity, and stream competence, it makes sense that Sites #2 to #4 located in the mid-reaches were generally associated with less substrate embeddedness and frequency of substrate less than GC in diameter compared to Sites #1 and #5. The mid-reaches of HCW were also associated with increased bedrock substrate and channel constraints (Figure 3), which indicated transport capacity in excess of sediment supply [52]. However, outside of the mid-reaches, increased substrate embeddedness was observed at agricultural Site # 1 and urban Site #5 (Figure 3).
The greatest rate of change in substrate embeddedness with downstream distance was observed in riffles located in the lower urban reaches where bank height exceeded 7 m (Figure 7). For example, percent embeddedness of riffle habitat decreased by at a rate of about 2% km−1 as agricultural land use decreased from 100 to 56.9%, forested land use increased from 0 to 35.9%, and urban land use increased from 0 to 4.7% over 22.8 km stream distance at Site #1 (Figure 7). However, percent embeddedness of riffles remained relatively constant (i.e., negligible rate of change) from agricultural Site #1 to sub-urban Site #3 located at the rural-urban interface where increased bedrock channel constraints were observed in HCW (Figure 7). Continuing downstream from sub-urban Site #3 to urban Site #4, substrate embeddedness of riffles began to increase by about 1.4% km−1 over approximately 8 km of stream distance. Further downstream, percent embeddedness of riffles increased rapidly by 5.3% km−1 between urban Sites #4 and #5 (Figure 7). Similar trends were observed in glides and pools as well (Figure 7). Thus, these results showed increased rate of change in percent embeddedness linked to degraded physical habitat (riffles, glides and pools) at agricultural Site #1 and urban Site #5, with a disproportionate rate of increase of substrate embeddedness in riffle habitat of the lower urban reaches pointing to symptoms of urban stream syndrome in HCW.
The observed channel incision in the lower reaches of the current work was corroborated by other authors who also reported urban land use influence on channel morphology [24,60,61]. Blann et al. [53] discussed how increased channel incision disrupts hydrologically and ecologically important stream-floodplain-riparian connectivity. In urban catchments, channel incision is often attributed to increased impervious surfaces and engineered waterways that connect impervious surfaces directly to stream channels [24,25]. Increased impervious surfaces have been shown to increase the volume and velocity of surface runoff in HCW [36] and elsewhere [24,25,62]. Increased surface runoff commonly translates to alterations to flow regimes (i.e., flow frequency, magnitude, timing, duration and rate of change), sediment transport regimes, water quality, and aquatic ecosystem health [24,25]. Clearly, there is a need to mitigate land use alterations to channel morphology via channel stabilization efforts in HCW and similar catchments globally. However, Vietz et al. [58] noted channel stabilization methods may not be sufficient to combat channel degradation in urbanized catchments. There is also a critical need for flow regime management efforts to reduce flow magnitude, frequency, and rate of change associated with alterations to channel morphology.
It was evident in the current work that simultaneously occurring agricultural and urban land uses exacerbated problems with substrate embeddedness in HCW particularly in the headwaters and lower urban reaches were bedrock constraints were less abundant. Previous studies also showed a general trend for suspended sediment and nutrients to decrease from Site #1 in the agricultural headwaters to Site #3 at the rural-urban interface of the watershed, and then, increase from Site #3 to urban Site #5 located near the watershed outlet in HCW [14,40,41]. In fact, significant relationships were observed between substrate embeddedness values reported in the current work and concentrations of suspended sediment (R2 = 0.798; p = 0.026), nitrite-N (R2 = 0.975; p = 0.001), and ammonia-N (R2 = 0.956; p = 0.002) reported in previous studies in HCW [14,41]. Notably, suspended sediment and total phosphorous yields were particularly high compared to other regions within the Mississippi River Basin. For example, observed suspended sediment yields exceeding 300 Mg km−2 year−1 were 54%, 80%, and 87% greater than sediment yields from the Ohio River, the Upper Mississippi River, and the Missouri River basins, respectively [14]. Total phosphorous yields (0.979 kg ha−1 year−1) were also high for the region [41]. In combination with results of LULC alterations to substrate embeddedness from the current work, results indicate increased sediment supply has degraded water quality, physical habitat, and aquatic refugia in HCW, especially in the agricultural headwaters and the lower urban reaches.
Collectively, previous studies showed streams of the Midwestern USA are commonly associated with increased sediment supply and degraded stream health [13,15,57]. For example, Gellis et al. [57] noted channel sources of sediment accounted for the majority (>50%) of bed sediment in 79% of 99 Midwestern watersheds sampled. Increased channel sediment supply and subsequent bed sedimentation can bury riffle habitat, reduce egg and fry survivorship, and lower prey densities [63,64]. Results from previous studies often focused on agricultural influence on stream bank erosion and sedimentation which was observed in the current work. However, it should be noted that urban land use may cause a greater influence on channel morphology and streambed sedimentation relative to agricultural land use. For example, substrate embeddedness was observed to decrease between Site #1 and #2 where agricultural land use accounted for 55% of total catchment area. Continuing downstream, substrate embeddedness began to increase between Sites #2 and #3 at the rural—urban interface of HCW which was associated with about 7 to 13% urban land use in this study. Thus, results from this study were in agreement with Paul and Meyer [65] who noted urban land use can account for a small percentage of total catchment area while causing a disproportionate influence on water quality and stream health relative to other land uses.

4.3. Management Implications and Future Research

The PHA and nested-scale experimental watershed study design presented in the current work is a useful model for managers that need to elucidate casual factors, target critical source areas, and thus, guide regional stream restoration efforts of mixed land use watersheds globally. The intensive field data collection associated with this assessment made it possible to estimate downstream changes in metrics important for understanding available stream physical habitat and biological refugia across an agricultural-urban land use gradient. Assuming that reduction of sediment supply below sediment transport capacity will help to decrease substrate embeddedness, management efforts might focus on reducing land use alterations to (1) erosive forces of surface runoff and streamflow, and (2) channel incision and bank losses in watersheds prone to increased sediment supply.
A complete decoupling of the simultaneously occuring natural and anthropogenic factors (legacy and ongoing) that influence channel geomorphology and streambed sedimentological characteristics was beyond the scope of the current work. Additionally, results from this study do not show temporal variablity of channel geomorphology and substrate composition. Thus, results may not reflect the total influence of past and present land use alterations that are expected to continue to alter channel geomorphology and substrate composition over long time periods. Future work focused on re-measuring channel geomorphology, substrate frequency and embeddedness may elucidate relationship between land use change and subsequent alterations to stream physical habitat.

5. Conclusions

Key findings from the current work point to (1) agricultural and urban land use alterations on channel geomorphology and stream substrate embeddedness, and (2) channel geomorphology as an indicator of stream substrate embeddedness in a mixed land use Midwestern stream. Expected relationships between drainage area and channel geomorphology were altered. Agricultural and urban land use explained nearly all of the variance in average width to depth ratios (R2 = 0.960; p = 0.020; n = 5), and maximum bank angle (R2 = 0.896; p = 0.052; n = 5). Also observed were reduced frequencies of riffle habitat at Site #1 in the agricultural headwaters (13.1%) and Site #5 (12.2%) in the lower urban reaches. Increased rate of change in percent embeddedness of riffle habitat exceeding 5% km−1 in the lower urban reaches indicated a disproportionate influence of urban land use on hydromorphology. Results showed nearly all the variability in channel unit frequency was explained by increased frequency of fine substrate type (0.927 ≤ R2 ≥ 0.854; p ≤ 0.016; n = 5). These results highlighted observed streambed siltation in Hinkson Creek, especially in the agricultural dominated headwaters and the lower urban reaches. Thus, results from this study point to a critical need to mitigate observed agricultural and urban land use impacts on stream hydrogeomorphology.
Given the influence of stream physical habitat on stream ecosystem health, stream physical habitat assessments are an integral component of regional stream restoration efforts. The robust observed data set collected for this study can provide critical information needed to guide regional policy development and watershed management efforts. Results of this study are particularly important for regional management efforts considering (1) rigorous hydrogeomorphic data sets are rare, and (2) the magnitude of alterations to channel geomorphology, and streambed composition in Midwestern streams of the USA. Results from this study quantitatively characterize channel geomorphology, substrate frequency, and embeddedness across 56 km stream length in the Midwestern USA where streams are commonly burdened with increased channel sediment supply, increased suspended sediment and total phosphorous, and problems associated with substrate siltation. While results from the current work are regionally applicable, key findings may also be useful to guide policy development and management decisions in physiographically similar watersheds globally.

Author Contributions

Conceptualization, J.A.H. and S.J.Z.; Formal Analysis, S.J.Z.; Investigation, J.A.H.; Data Curation, S.J.Z.; Writing-Original Draft Preparation, S.J.Z.; Writing-Review & Editing, J.A.H.; Supervision, J.A.H.; Project Administration, J.A.H.; Funding Acquisition, J.A.H.

Funding

This research was funded by the Missouri Department of Conservation and the U.S. Environmental Protection Agency Region 7 through the Missouri Department of Natural Resources (P.N: G08-NPS-17) under Section 319 of the Clean Water Act and through joint agreement of the University of Missouri, the City of Columbia, and Boone County Public works and partners of the Hinkson Creek Collaborative Adaptive Management (CAM) program. Additional funding was provided by the National Science Foundation under Award Number OIA-1458952, the USDA National Institute of Food and Agriculture, Hatch project accession number 1011536, and the West Virginia Agricultural and Forestry Experiment Station. Results presented may not reflect the views of the sponsors and no official endorsement should be inferred.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Van Looy, K.; Tonkin, J.D.; Floury, M.; Leigh, C.; Soininen, J.; Larsen, S.; Datry, T. The three Rs of river ecosystem resilience: Resources, recruitment, and refugia. River Res. Appl. 2019, 35, 107–120. [Google Scholar] [CrossRef] [Green Version]
  2. Keppel, G.; Van Niel, K.P.; Wardell-Johnson, G.W.; Yates, C.J.; Byrne, M.; Mucina, L.; Franklin, S.E. Refugia: Identifying and understanding safe havens for biodiversity under climate change. Glob. Ecol. Biogeogr. 2012, 21, 393–404. [Google Scholar] [CrossRef]
  3. Scherrer, D.; Körner, C. Topographically controlled thermal-habitat differentiation buffers alpine plant diversity against climate warming. J. Biogeogr. 2011, 38, 406–416. [Google Scholar] [CrossRef]
  4. Maddock, I. The importance of physical habitat assessment for evaluating river health. Freshw. Biol. 1999, 41, 373–391. [Google Scholar] [CrossRef]
  5. Poff, N.L.; Ward, J.V. Physical habitat template of lotic systems: Recovery in the context of historical pattern of spatiotemporal heterogeneity. Environ. Manag. 1990, 14, 629. [Google Scholar] [CrossRef]
  6. Klos, P.Z.; Rosenberry, D.O.; Nelson, G.R. Influence of hyporheic exchange, substrate distribution, and other physically linked hydrogeomorphic characteristics on abundance of freshwater mussels. Ecohydrology 2015, 8, 1284–1291. [Google Scholar] [CrossRef]
  7. Duan, X.; Wang, Z.; Zang, K. Effect of streambed sediment on benthic ecology. Int. J. Sediment Res. 2009, 24, 325–338. [Google Scholar] [CrossRef]
  8. Duan, X.; Wang, Z.; Tian, S. Effect of streambed substrate on macroinvertebrate biodiversity. Front. Environ. Sci. Eng. China 2008, 2, 122–128. [Google Scholar] [CrossRef]
  9. Erman, D.C.; Erman, N.A. The response of stream macroinvertebrates to substrate size and heterogeneity. Hydrobiologia 1984, 108, 75–82. [Google Scholar] [CrossRef]
  10. Nichols, J.; Hubbart, J.A.; Poulton, B.C. Using macroinvertebrate assemblages and multiple stressors to infer urban stream system condition: A case study in the central US. Urban Ecosyst. 2016, 19, 679–704. [Google Scholar] [CrossRef]
  11. Blann, K.L.; Anderson, J.L.; Sands, G.R.; Vondracek, B. Effects of Agricultural Drainage on Aquatic Ecosystems: A Review. Crit. Rev. Environ. Sci. Technol. 2009, 39, 909–1001. [Google Scholar] [CrossRef]
  12. Neitsch, S.L.; Arnold, J.G.; Kiniry, J.R.; Williams, J.R. Soil and Water Assessment Tool Theoretical Documentation Version 2009; Texas Water Resources Institute: College Station, TX, USA, 2011; Available online: https://swat.tamu.edu/documentation/ (accessed on 26 August 2018).
  13. Gellis, A.C.; Myers, M.K.; Noe, G.B.; Hupp, C.R.; Schenk, E.R.; Myers, L. Storms, channel changes, and a sediment budget for an urban-suburban stream, Difficult Run, Virginia, USA. Geomorphology 2017, 278, 128–148. [Google Scholar] [CrossRef]
  14. Zeiger, S.J.; Hubbart, J.A. Quantifying suspended sediment flux in a mixed land use urbanizing watershed using a nested-scale study design. Sci. Total Environ. 2015, 242, 315–323. [Google Scholar] [CrossRef] [PubMed]
  15. Willett, C.D.; Lerch, R.N.; Schultz, R.C.; Berges, S.A.; Peacher, R.D.; Isenhart, T.M. Streambank erosion in two watersheds of the Central Claypan Region of Missouri, United States. J. Soil Water Conserv. 2012, 67, 249–263. [Google Scholar] [CrossRef] [Green Version]
  16. Huang, D. Quantifying Stream Bank Erosion Deposition Rates in a Central, U.S. Urban Watershed. Master’s Thesis, University of Missouri, Columbia, MO, USA, 2012. [Google Scholar]
  17. Waite, I.R.; Van Metre, P.C. Multistressor predictive models of invertebrate condition in the Corn Belt, USA. Freshw. Sci. 2017, 36, 901–914. [Google Scholar] [CrossRef]
  18. Meador, M.R.; Frey, J.W. Relative Importance of Water-Quality Stressors in Predicting Fish Community Responses in Midwestern Streams. JAWRA J. Am. Water Resour. Assoc. 2018, 54, 708–723. [Google Scholar] [CrossRef]
  19. Munn, M.D.; Waite, I.; Konrad, C.P. Assessing the influence of multiple stressors on stream diatom metrics in the upper Midwest, USA. Ecol. Indic. 2018, 85, 1239–1248. [Google Scholar] [CrossRef]
  20. Poff, N.L.; Allan, J.D.; Bain, M.B.; Karr, J.R.; Prestegaard, K.L.; Richter, B.D.; Stromberg, J.C. The natural flow regime. BioScience 1997, 47, 769–784. [Google Scholar] [CrossRef]
  21. Berkman, H.E.; Rabeni, C.F. Effect of siltation on stream fish communities. Environ. Biol. Fishes 1987, 18, 285–294. [Google Scholar] [CrossRef]
  22. Tufekcioglu, M.; Isenhart, T.M.; Schultz, R.C.; Bear, D.A.; Kovar, J.L.; Russell, J.R. Stream bank erosion as a source of sediment and phosphorus in grazed pastures of the Rathbun Lake Watershed in southern Iowa, United States. J. Soil Water Conserv. 2012, 67, 545–555. [Google Scholar] [CrossRef] [Green Version]
  23. Dosskey, M.G.; Vidon, P.; Gurwick, N.P.; Allan, C.J.; Duval, T.P.; Lowrance, R. The role of riparian vegetation in protecting and improving chemical water quality in streams1. JAWRA J. Am. Water Resour. Assoc. 2010, 46, 261–277. [Google Scholar] [CrossRef]
  24. Walsh, C.J.; Roy, A.H.; Feminella, J.W.; Cottingham, P.D.; Groffman, P.M.; Morgan, I.I.R.P. The urban stream syndrome: Current knowledge and the search for a cure. J. N. Am. Benthol. Soc. 2005, 24, 706–723. [Google Scholar] [CrossRef]
  25. Meyer, J.L.; Paul, M.J.; Taulbee, W.K. Stream ecosystem function in urbanizing landscapes. J. N. Am. Benthol. Soc. 2005, 24, 602–612. [Google Scholar] [CrossRef]
  26. Poff, N.L.; Richter, B.D.; Arthington, A.H.; Bunn, S.E.; Naiman, R.J.; Kendy, E.; Henriksen, J. The ecological limits of hydrologic alteration (ELOHA): A new framework for developing regional environmental flow standards. Freshw. Biol. 2010, 55, 147–170. [Google Scholar] [CrossRef]
  27. United States Census Bureau (USCB). U.S. Census Bureau: State and County Quick Facts. 2018. Available online: http://quickfacts.census.gov/qfd/states/29000.html (accessed on 13 May 2019).
  28. Hubbart, J.A.; Zell, C. Considering streamflow trend analyses uncertainty in urbanizing watersheds: A base flow case study in the central United States. Earth Interact. 2013, 17, 1–28. [Google Scholar] [CrossRef]
  29. Lerch, R.N.; Kitchen, N.R.; Baffaut, C.; Vories, E.D. Long-term agroecosystem research in the Central Mississippi River Basin: Goodwater Creek Experimental Watershed and regional nutrient water quality data. J. Environ. Qual. 2015, 44, 37–43. [Google Scholar] [CrossRef] [PubMed]
  30. Miller, D.; Vandike, J. Groundwater Resources of Missouri; Missouri Department of Natural Resources Division of Geology and Land Survey Water Resources: Jefferson City, MO, USA, 1997; Volume 46, p. 136. [Google Scholar]
  31. Hubbart, J.A.; Holmes, J.; Bowman, G. Integrating science based decision making and TMDL allocations in urbanizing watersheds. Watershed Sci. Bull. 2010, 1, 19–24. [Google Scholar]
  32. Hubbart, J.A.; Kellner, E.; Hooper, L.; Lupo, A.R.; Market, P.S.; Guinan, P.E.; Stephan, K.; Fox, N.I.; Svoma, B.M. Localized Climate and Surface Energy Flux Alterations across an Urban Gradient in the Central U.S. Energies 2014, 7, 1770–1791. [Google Scholar] [CrossRef]
  33. Zeiger, S.J.; Hubbart, J.A. Rainfall-Stream Flow Responses in a Mixed land use Municipal Watershed of the Central, U.S.A. Environ. Earth Sci. 2018, 77, 438. [Google Scholar] [CrossRef]
  34. Zeiger, S.J.; Hubbart, J.A. Quantifying land use influences on event-based flow frequency, timing, magnitude, and rate of change in an urbanizing watershed of the central USA. Environ. Earth Sci. 2018, 77, 107. [Google Scholar] [CrossRef]
  35. Wei, L.; Hubbart, J.A.; Zhou, H. Variable streamflow contributions in nested subwatersheds of a US midwestern urban watershed. Water Resour. Manag. 2018, 32, 213–228. [Google Scholar] [CrossRef]
  36. Zeiger, S.J.; Hubbart, J.A. Assessing environmental flow targets using pre-settlement land cover: A SWAT modeling application. Water 2018, 10, 791. [Google Scholar] [CrossRef]
  37. Zeiger, S.J.; Hubbart, J.A. Quantifying urban stormwater temperature surges: A central US watershed study. Hydrology 2015, 2, 193–209. [Google Scholar] [CrossRef]
  38. Zeiger, S.J.; Hubbart, J.A.; Anderson, S.H.; Stambaugh, M.C. Quantifying and modelling urban stream temperature: A central US watershed study. Hydrol. Process. 2015, 30, 503–514. [Google Scholar] [CrossRef]
  39. Kellner, E.; Hubbart, J.A. Improving understanding of mixed land use watershed suspended sediment regimes: Mechanistic progress through high-frequency sampling. Sci. Total Environ. 2017, 598, 228–238. [Google Scholar] [CrossRef] [PubMed]
  40. Kellner, E.; Hubbart, J.A. Quantifying Urban Land-Use Impacts on Suspended Sediment Particle Size Class Distribution: A Method and Case Study. Stormwater J. 2014, 15, 40–50. [Google Scholar]
  41. Zeiger, S.J.; Hubbart, J.A. Nested-Scale Nutrient Yields from a Mixed land use Urbanizing Watershed. Hydrol. Process. 2016, 30, 1475–1490. [Google Scholar] [CrossRef]
  42. Hubbart, J.A.; Kellner, E.; Hooper, L.W.; Zeiger, S.J. Quantifying loading, toxic concentrations, and systemic persistence of chloride in a contemporary mixed land use watershed using an experimental watershed approach. Sci. Total Environ. 2017, 581, 822–832. [Google Scholar] [CrossRef]
  43. Hubbart, J.A.; Kellner, E.; Kinder, P.; Stephan, K. Challenges in aquatic physical habitat assessment: Improving conservation and restoration decisions for contemporary watersheds. Challenges 2017, 8, 31. [Google Scholar] [CrossRef]
  44. Hooper, L.; Hubbart, J.A. A rapid physical habitat assessment of wadeable streams for mixed land use watersheds. Hydrology 2016, 3, 37. [Google Scholar] [CrossRef]
  45. Hooper, L. A Stream Physical Habitat Assessment of an Urbanizing Stream of the Central U.S.A. Master’s Thesis, University of Missouri, Columbia, SC, USA, 2015. [Google Scholar]
  46. Harrelson, C.C.; Rawlins, C.L.; Potyondy, J.P. Stream Channel Reference Sites: An Illustrated Guide to Field Technique; Gen. Tech. Rep. RM-245; US Department of Agriculture, Forest Service, Rocky Mountain Forest and Range Experiment Station: Fort Collins, CO, USA, 1994. [Google Scholar]
  47. Yanites, B.J.; Tucker, G.E. Controls and limits on bedrock channel geometry. J. Geophys. Res. 2010, 115, F04019. [Google Scholar] [CrossRef]
  48. Monk, W.A.; Wood, P.J.; Hannah, D.M.; Wilson, D.A. Macroinvertebrate community response to inter-annual and regional river flow regime dynamics. River Res. Appl. 2008, 24, 988–1001. [Google Scholar] [CrossRef] [Green Version]
  49. Stoline, M.R. The status of multiple comparisons: Simultaneous estimation of all pairwise comparisons in one-way ANOVA designs. Am. Stat. 1981, 35, 134–141. [Google Scholar]
  50. Peck, D.V.; Herlihy, A.T.; Hill, B.H.; Hughes, R.M.; Kaufmann, P.R.; Klemm, D.J.; Lazorchak, J.M.; Cappaert, M. Environmental Monitoring and Assessment Program-Surface Waters Western Pilot Study: Field Operations Manual for Wadeable Streams; EPA/620/R-06/003; U.S. Environmental Protection Agency, Office of Research and Development: Washington, DC, USA, 2006; pp. 128–130.
  51. Wolman, M.G. A method of sampling coarse river-bed material. Trans. Am. Geophys. Union 1954, 35, 951–956. [Google Scholar] [CrossRef]
  52. Montgomery, D.R.; Buffington, J.M.; Peterson, N.P.; Schuett-Hames, D.; Quinn, T.P. Stream-bed scour, egg burial depths, and the influence of salmonid spawning on bed surface mobility and embryo survival. Can. J. Fish. Aquat. Sci. 1996, 53, 1061–1070. [Google Scholar] [CrossRef]
  53. Diana, M.; Allan, J.D.; Infante, D. The influence of physical habitat and land use on stream fish assemblages in southeastern Michigan. Am. Fish. Soc. Symp. 2006, 48, 359–374. [Google Scholar]
  54. Zaimes, G.N.; Schultz, R.C.; Isenhart, T.M. Stream bank erosion adjacent to riparian forest buffers, row-crop fields, and continuously-grazed pastures along Bear Creek in central Iowa. J. Soil Water Conserv. 2004, 59, 19–27. [Google Scholar]
  55. Urban, M.L.; Rhoads, B.L. Catastrophic Human-Induced Change in Stream-Channel Planform and Geometry in an Agricultural Watershed, Illinois, USA. Ann. Assoc. Am. Geogr. 2003, 93, 783–796. [Google Scholar] [CrossRef]
  56. Simon, A.; Curini, A.; Darby, S.E.; Langendoen, E.J. Bank and near-bank processes in an incised channel. Geomorphology 2000, 35, 193–217. [Google Scholar] [CrossRef]
  57. Gellis, A.C.; Fuller, C.C.; Van Metre, P.C. Sources and ages of fine-grained sediment to streams using fallout radionuclides in the Midwestern United States. J. Environ. Manag. 2017, 194, 73–85. [Google Scholar] [CrossRef] [Green Version]
  58. Vietz, G.J.; Walsh, C.J.; Fletcher, T.D. Urban hydrogeomorphology and the urban stream syndrome: Treating the symptoms and causes of geomorphic change. Prog. Phys. Geogr. 2016, 40, 480–492. [Google Scholar] [CrossRef]
  59. Shields, F.D.; Lizotte, R.E.; Knight, S.S.; Cooper, C.M.; Wilcox, D. The stream channel incision syndrome and water quality. Ecol. Eng. 2010, 36, 78–90. [Google Scholar] [CrossRef]
  60. Jordan, B.A.; Annable, W.K.; Watson, C.C.; Sen, D. Contrasting stream stability characteristics in adjacent urban watersheds: Santa Clara Valley, California. River Res. Appl. 2010, 26, 1281–1297. [Google Scholar] [CrossRef]
  61. Fletcher, T.D.; Vietz, G.; Walsh, C.J. Protection of stream ecosystems from urban stormwater runoff: The multiple benefits of an ecohydrological approach. Prog. Phys. Geogr. 2014, 38, 543–555. [Google Scholar] [CrossRef]
  62. Trimble, S.W. Fluvial processes, morphology and sediment budgets in the Coon Creek Basin, WI, USA, 1975–1993. Geomorphology 2009, 108, 8–23. [Google Scholar] [CrossRef]
  63. Walters, D.M.; Leigh, D.S.; Bearden, A.B. Urbanization, sedimentation, and the homogenization of fish assemblages in the Etowah River Basin, USA. In The Interactions between Sediments and Water; Springer: Dordrecht, the Netherlands, 2003; pp. 5–10. [Google Scholar]
  64. Waters, T.F. Sediment in Streams: Sources, Biological Effects, and Control; American Fisheries Society: Bethesda, MA, USA, 1995. [Google Scholar]
  65. Paul, M.J.; Meyer, J.L. Streams in the urban landscape. Annu. Rev. Ecol. Syst. 2001, 32, 333–365. [Google Scholar] [CrossRef]
Figure 1. Land use and land cover of Hinkson Creek Watershed, Missouri, USA. Five nested gauging sites (numbered 1 to 5) and corresponding sub-basins are shown. Site #5 was located near the watershed outlet, Sites #1 to #4 nested within.
Figure 1. Land use and land cover of Hinkson Creek Watershed, Missouri, USA. Five nested gauging sites (numbered 1 to 5) and corresponding sub-basins are shown. Site #5 was located near the watershed outlet, Sites #1 to #4 nested within.
Water 11 01088 g001
Figure 2. Power relationships between drainage area and observed channel morphological metrics in Hinkson Creek Watershed, Missouri, USA.
Figure 2. Power relationships between drainage area and observed channel morphological metrics in Hinkson Creek Watershed, Missouri, USA.
Water 11 01088 g002
Figure 3. Examples of channel characteristics associated with the agricultural headwaters (a), rural-urban mid-reaches (b,c), and lower urban reaches (d) of Hinkson Creek Watershed, USA.
Figure 3. Examples of channel characteristics associated with the agricultural headwaters (a), rural-urban mid-reaches (b,c), and lower urban reaches (d) of Hinkson Creek Watershed, USA.
Water 11 01088 g003
Figure 4. Trend lines associated with channel morphology variables with downstream distance (km) in Hinkson Creek Watershed, Missouri, USA. Vertical lines show location of five nested gauging sites.
Figure 4. Trend lines associated with channel morphology variables with downstream distance (km) in Hinkson Creek Watershed, Missouri, USA. Vertical lines show location of five nested gauging sites.
Water 11 01088 g004
Figure 5. Explained variance between channel morphology metrics (x-axis) and substrate embeddedness from left bank (lft. bnk.) to right bank (rt. bnk.) associated with riffles, glides, and pools (y-axis) at five sub-basins of Hinkson Creek Watershed, Missouri, USA. No glide channel units were observed at Site #1, instead, percent embeddedness of dry channel units is shown.
Figure 5. Explained variance between channel morphology metrics (x-axis) and substrate embeddedness from left bank (lft. bnk.) to right bank (rt. bnk.) associated with riffles, glides, and pools (y-axis) at five sub-basins of Hinkson Creek Watershed, Missouri, USA. No glide channel units were observed at Site #1, instead, percent embeddedness of dry channel units is shown.
Water 11 01088 g005
Figure 6. Substrate frequency and embeddedness from left bank (lft. bnk.) to right bank (rt. bnk.) associated with riffles, glides, and pools located within five nested sub-basins (numbered 1 to 5) of Hinkson Creek Watershed, Missouri, USA. Sub-basin #1 was located in the headwaters. Sub-basin #5 was located near the watershed outlet. No glide channel units were observed at Site #1, instead, percent embeddedness of dry channel units is shown. Substrate definitions are presented in Table 2.
Figure 6. Substrate frequency and embeddedness from left bank (lft. bnk.) to right bank (rt. bnk.) associated with riffles, glides, and pools located within five nested sub-basins (numbered 1 to 5) of Hinkson Creek Watershed, Missouri, USA. Sub-basin #1 was located in the headwaters. Sub-basin #5 was located near the watershed outlet. No glide channel units were observed at Site #1, instead, percent embeddedness of dry channel units is shown. Substrate definitions are presented in Table 2.
Water 11 01088 g006
Figure 7. Watershed characteristics and trend lines associated with substrate embeddedness in pools, riffles, and glides with downstream distance (km) in Hinkson Creek Watershed, Missouri, USA. Vertical lines show location of five nested gauging sites.
Figure 7. Watershed characteristics and trend lines associated with substrate embeddedness in pools, riffles, and glides with downstream distance (km) in Hinkson Creek Watershed, Missouri, USA. Vertical lines show location of five nested gauging sites.
Water 11 01088 g007
Table 1. Cumulative land use and land cover (LULC), drainage area, and stream length corresponding to each gauging site located in Hinkson Creek Watershed (HCW), Missouri, USA. Percent cumulative LULC is shown parenthetically.
Table 1. Cumulative land use and land cover (LULC), drainage area, and stream length corresponding to each gauging site located in Hinkson Creek Watershed (HCW), Missouri, USA. Percent cumulative LULC is shown parenthetically.
VariableSite #1Site #2Site #3
[km2 (%)]
Site #4Site #5HCW
Agricultural45.0 (57.0)56.4 (54.9)57.6 (49.5)78.5 (44.1)79.7 (39.4)85.4 (36.7)
Forested28.4 (35.9)37.5 (36.4)41.1 (35.4)62.8 (34.5)68.6 (33.1)74.9 (32.2)
Urban3.7 (4.7)6.6 (6.4)15 (13.0)37.1 (20.4)54.9 (26.5)67.6 (29.0)
Wetland1.9 (2.4)2.4 (2.3)2.5 (2.1)3.6 (2.0)4.4 (2.0)4.9 (2.1)
Total area79.0102.9116.2182.0207.5232.8
Stream length §22.829.835.443.653.056.1
§ Stream length is shown in km.
Table 2. Definition of streambed substrate variables used in this study.
Table 2. Definition of streambed substrate variables used in this study.
CodeDescription
LWD PILELarge woody debris
VGVegetation
SHShell
CUCulvert
BRBridge
RPRip-rap
WDWood
RSSmooth Bedrock (>4000 mm)
RRRough Bedrock (>4000 mm)
RCConcrete/asphalt (>4000 mm)
XBLarge boulder (1000 to 4000 mm)
SBSmall boulder (256 to 1000 mm)
CBCobble (64 to 256 mm)
GCCoarse gravel (16 to 64 mm)
GFFine gravel (2 to 16 mm)
SASand (0.06 to 2 mm)
FNSilt/clay/muck (<0.06 mm)
Table 3. Channel morphology metrics measured at survey sites (n = 561) and averaged within five nested sub-basins located in Hinkson Creek Watershed, Missouri, USA. Sub-basins were delineated between nested gauging sites and were thus not cumulative.
Table 3. Channel morphology metrics measured at survey sites (n = 561) and averaged within five nested sub-basins located in Hinkson Creek Watershed, Missouri, USA. Sub-basins were delineated between nested gauging sites and were thus not cumulative.
Variable §Site #1Site #2Site #3Site #4Site #5
Channel width (m km−2) 0.14 abc0.19 ade0.15 df0.16 bg0.12 cefg
Bankfull width (m km−2) 0.23 a0.27 abcd0.23 b0.23 c0.21 d
Bank height (m km−2) 0.47 abc0.39 ad0.44 e0.38 bf0.56 cdef
Bankfull depth (m km−2) 0.42 abcd0.36 ae0.37 bf0.34 cg0.49 defg
Min bank angle (°)28.0 ab22.5 ac22.0 bd25.428.0 cd
Max bank angle (°)47.8 abc38.7 a40.3 b39.9 c44.4
Width: depth (--)6.3 abc10.2 ade8.5 bdfg10.3 cfh6.6 egh
Bed slope (%)0.14 abcd0.22 aef0.24 agh0.18 aeg0.16 dfh
§ Averages with corresponding letters (a–h) indicate significant differences (α = 0.05). Channel geomorphic variables normalized by drainage area.
Table 4. Bank to bank variability of streambed substrate frequency in Hinkson Creek Watershed, Missouri, USA. Each datum is a percent of a substrate type observed at a bank position along the entire 56km stream length of the reference stream. Substrate definitions are presented in Table 2.
Table 4. Bank to bank variability of streambed substrate frequency in Hinkson Creek Watershed, Missouri, USA. Each datum is a percent of a substrate type observed at a bank position along the entire 56km stream length of the reference stream. Substrate definitions are presented in Table 2.
Substrate TypeLeft BankLeft CenterCenter (%)Right CenterRight Bank
Silt/clay/muck (FN)57.917.511.617.662.0
Sand (SA)17.328.833.934.317.5
Fine gravel (GF)0.25.16.14.10.4
Coarse gravel (GC)3.520.121.316.83.3
Cobble (CB)5.613.013.815.94.1
Small boulder (SB)3.74.34.34.33.0
Large boulder (XB)1.40.61.20.91.2
Rough Bedrock (RR)4.25.93.72.51.2
Smooth Bedrock (RS)1.52.01.41.11.1
Wood (WD)2.90.40.60.42.5
Rip-rap (RP)0.80.30.10.11.7
Bridge (BR)----0.1----
Culvert (CU)----------
Concrete/asphalt (RC)0.1------0.3
Shell (SH)----0.1----
Vegetation (VG)0.92.01.62.11.8
Large woody debris (LWD PILE)0.10.10.10.1--
Table 5. Streambed substrate frequency observed within five nested sub-basins in Hinkson Creek Watershed, Missouri, USA. Each datum is a percent of a substrate type. Sub-basins were delineated between nested gauging sites, and thus, were not cumulative watersheds excepting Site #1 located in the headwaters. Substrate definitions are presented in Table 2.
Table 5. Streambed substrate frequency observed within five nested sub-basins in Hinkson Creek Watershed, Missouri, USA. Each datum is a percent of a substrate type. Sub-basins were delineated between nested gauging sites, and thus, were not cumulative watersheds excepting Site #1 located in the headwaters. Substrate definitions are presented in Table 2.
Substrate TypeSite #1Site #2Site #3 (%)Site #4Site #5
Silt/clay/muck (FN)402623942
Sand (SA)3218132925
Fine gravel (GF)41134
Coarse gravel (GC)919161814
Cobble (CB)71621175
Small boulder (SB)259102
Large boulder (XB)01331
Rough Bedrock (RR)310911
Smooth Bedrock (RS)00280
Wood (WD)22001
Rip-rap (RP)01101
Bridge (BR)----0----
Culvert (CU)----------
Concrete/asphalt (RC)0------0
Shell (SH)----0----
Vegetation (VG)22104
Large woody debris (LWD PILE)0000--
Table 6. Bank to bank variability of streambed substrate embeddedness in riffles, glides, and pools measured at survey sites (n = 561) and averaged within five nested sub-basins located in Hinkson Creek Watershed, Missouri, USA. Each datum shows substrate embeddedness as a percent. Sub-basins were delineated between nested gauging sites, and thus, were not cumulative watersheds excepting Site #1 located in the headwaters.
Table 6. Bank to bank variability of streambed substrate embeddedness in riffles, glides, and pools measured at survey sites (n = 561) and averaged within five nested sub-basins located in Hinkson Creek Watershed, Missouri, USA. Each datum shows substrate embeddedness as a percent. Sub-basins were delineated between nested gauging sites, and thus, were not cumulative watersheds excepting Site #1 located in the headwaters.
Bank PositionSite (#)Riffles Glides (%) *Pools
Left bankSite #18310088
Site #2766771
Site #3768167
Site #4727172
Site #5856588
Left centerSite #1569175
Site #2272342
Site #3395644
Site #4465954
Site #5593270
CenterSite #1518076
Site #2311846
Site #3403939
Site #4435152
Site #5617075
Right centerSite #1529379
Site #2302757
Site #3344347
Site #4415952
Site #5724276
Right bankSite #17010090
Site #2777880
Site #3706986
Site #4657377
Site #5887890
* Percent embeddedness of glides at Site #1 indicate dry channel units.
Table 7. Channel unit frequency measured at survey sites (n = 561) and averaged within five nested sub-basins located in Hinkson Creek Watershed, Missouri, USA.
Table 7. Channel unit frequency measured at survey sites (n = 561) and averaged within five nested sub-basins located in Hinkson Creek Watershed, Missouri, USA.
Site (#)PoolsRiffles (%)GlidesDry
Site #178.613.10.08.3
Site #272.123.64.30.0
Site #364.325.010.70.0
Site #456.128.015.90.0
Site #585.712.22.10.0

Share and Cite

MDPI and ACS Style

Zeiger, S.J.; Hubbart, J.A. Characterizing Land Use Impacts on Channel Geomorphology and Streambed Sedimentological Characteristics. Water 2019, 11, 1088. https://doi.org/10.3390/w11051088

AMA Style

Zeiger SJ, Hubbart JA. Characterizing Land Use Impacts on Channel Geomorphology and Streambed Sedimentological Characteristics. Water. 2019; 11(5):1088. https://doi.org/10.3390/w11051088

Chicago/Turabian Style

Zeiger, Sean J., and Jason A. Hubbart. 2019. "Characterizing Land Use Impacts on Channel Geomorphology and Streambed Sedimentological Characteristics" Water 11, no. 5: 1088. https://doi.org/10.3390/w11051088

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop