Next Article in Journal
Fate of Heavy Metals in the Surface Water-Dump Rock System of the Mine Lupikko I (Karelia): Field Observations and Geochemical Modeling
Previous Article in Journal
Optimization Method to Determine the Kinetic Rate Constants for the Removal of Benzo[a]pyrene and Anthracene in Water through the Fenton Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Streambed Median Grain Size (D50) across the Contiguous U.S.

1
Department of Civil, Architectural and Environmental Engineering, North Carolina A&T State University, Greensboro, NC 27411, USA
2
Department of Geoscience, Baylor University, Waco, TX 76798, USA
3
Grassland Soil and Water Research Laboratory USDA-ARS, Temple, TX 76502, USA
4
Department of Ecoscience, Aarhus University, Vejlsøvej 25, DK-8600 Silkeborg, Denmark
*
Author to whom correspondence should be addressed.
Water 2022, 14(21), 3378; https://doi.org/10.3390/w14213378
Submission received: 15 September 2022 / Revised: 16 October 2022 / Accepted: 20 October 2022 / Published: 25 October 2022
(This article belongs to the Section Water Erosion and Sediment Transport)

Abstract

:
The streambed median grain size (D50) has been an integral part of many sediment transport and stream power equations seeking to characterize stream channel stability conditions. However, its previous usage is constrained by regional applicability, localization of datasets, and a limited number of data points. This study uses a large and geographically diverse data set (n > 2400), from five published sources, to present quantitative information and assess the distribution of D50 data across the contiguous U.S. Spatial distribution was analyzed based on the three regional frameworks: Physiographic Provinces, Level III Ecoregions, and Hydrologic Landscape Regions (HLRs). Gravel was found to be the dominant streambed material in most Physiographic Provinces. Regions with a humid climate, permeable soil, and plateaus exhibit a higher average D50 than regions with other climate, geologic texture, and landscape forms. Further analysis of all data across the U.S. using smoothed spatial maps showed the dominance of sand and fine gravel in streams located in the central and southern U.S., and the dominance of coarse gravel and cobbles in the northeastern U.S.

1. Introduction

Streams are composed of active channels and floodplains, which play a vital role in the water cycle, ensuring the sustainability of the natural water balance. They are characterized by dynamic processes that maintain a state of equilibrium under variable sediment and water throughputs [1,2]. A stream can maintain its viability so long as it is under a stable equilibrium condition. [3] indicated a stable river can sustain sufficient transport capacity to convey the sediment and water flow coming from its watershed while preserving its channel pattern and dimension. However, natural and anthropogenic disturbances alter sediment and water flow [4], forcing streams and rivers to readjust their shape, size, and slope. Such perturbations result in reduced water storage capacity, diminished channel roughness, and shortened flow paths ensuring the prevalence of higher sediment transport and lower sediment storage capacity [5]. Consequently, the extra energy acquired through such modifications will be invested in scouring the streambed.
The resulting channel alterations can have a pronounced negative impact on nearby infrastructures [6,7], the benthic environment [8,9], and on the sediment transport [10]. The severity of these changes depends upon the landscape character and geological conditions [11], which mainly influence the magnitude of geomorphic driving and resisting forces [4]. This imbalance can be expressed by the stream power proportionality proposed by [1], which correlates to available stream power (the discharge-slope product) and the discharge of streambed material sediment. [12] employed parameters such as top width, average flow depth, flow velocity, longitudinal slope, and bankfull discharge to assess channel stability. The median grain size of mobile streambed sediments (D50) has been investigated as an indicator of channel stability [13,14,15].
D50 of streambed sediments is a valuable parameter in assessing watershed management strategies, analyzing and predicting stream behavior, and improving overall understanding of stream processes [16]. Previous studies recognized the influence of streambed material on channel stability and channel geometry [17,18]. The potential to assess channel stability and channel processes through median grain size data justifies its ubiquitous use in hydrological and geomorphological studies.
Golden and Springer [19] and Hack [20] used D50 in combination with drainage area to establish a regression relation for capturing the variation in the streambed slope. Grain-size data has also been used to predict the onset of scouring and streambed instability. Ref [16,21,22,23] used D50 as one of the inputs for estimating the critical specific stream power. Parker et al. [24] proposed a simplified expression of specific stream power using D50 as the sole dependent variable. Similarly, multiple authors [8,25,26] have used D50 as an integral parameter in developing formulas for critical shear stress and bedload. Kiel [27] developed a map of D50 using data collected from more than 400 sites to estimate the horizontal hydraulic conductivity of riverbanks. These studies demonstrated the extensive applicability of D50 in streambed stability analysis, but they also showed its limitations for its application in large regional scales.
This study developed a comprehensive database of D50 across the contiguous U.S. The data is presented in a stratified format using three regional frameworks: Physiographic Provinces, Level III Ecoregions, and Hydrologic Landscape Regions (HLRs). These frameworks use different sets of variables to define areas with similar hydrological processes. Many hydrological studies use these regional stratifications because they effectively characterize processes [28,29,30,31,32]. We prepared maps that depicted the distribution of streambed material over the contiguous U.S. Our dataset and analysis are applicable to diverse hydrological studies, but particularly those that address streambed stability.

2. Data

A large and geographically diverse data set (n > 2400) was collected for D50 from five different sources (Table 1). Some of the sources were comprehensive in terms of the hydrological parameters that were included within their dataset along with D50. The majority of the data points were obtained from the U.S. EPA [33]. Figure 1 shows the distribution of data per source over the contiguous U.S. An abundance of the data was found to be located in the Northeast, the Central U.S., and to some extent on the West Coast.
Preliminary data quality check revealed several duplicate and erroneous data points, which were removed from the analysis. Some data points had the same latitude and longitude but different D50 values. The USEPA dataset [33] contained 101 duplicate points, but the other four sources contained only 5 duplicates.
Box plot of D50 (Figure 2a) for each of the dataset revealed outliers that ranged from 6 to 84. The Interquartile range, defined as the difference between the 75th and 25th percentile of a dataset, was used as a measure of the statistical dispersion. Table 2 presents the descriptive statistics for each dataset after removing the outliers. According to [37], outliers can either be removed, used with weight adjustment, or be replaced with statistical estimators. We merely removed all outliers because there were relatively few of them. 139 data points were removed in total, which reduced the total number of points for further analyses to about 2300. Figure 2b shows the box plots after the removal of outlier. It shows the new outliers, but the interquartile range is much smaller and these remaining outliers will have minimal influence on the overall analysis. All D50 data and related information are provided in Tables S1 and S2 in the supplementary material of this paper.

3. Distribution of D50 by Regions

The spatial and temporal changes of climatic, geographic, and ecological variables influence hydrological processes. Despite their diverse and complex nature, these variables share common features [38]. Conceptual and analytical frameworks are therefore developed on a regional basis where similarities in these factors are used to define consistent hydrological characteristics of the region [32]. In this study, the D50 data were stratified using three major regionalization scheme in the U.S.: Level III Ecoregions [39], Physiographic Provinces [40], and Hydrologic Landscape Regions (HLR) [38]. The spatial distribution of the D50 data varies based on the geographic regions.
Level III Ecoregions. Ecoregions are geographical divisions characterized by similar ecological systems, organisms and environments [41]. The contiguous U.S. has been divided into 85 Level III Ecoregions [39], which have already been employed to regionalize hydraulic geometry parameters [42]. Figure 3 shows the distribution of D50 data points across the Level III Ecoregions. Of the 85 regions, 13 contain no data points. Moreover, the number of data points also varied considerably among the other 72 ecoregions. For example, the Ridge and Valley ecoregion contained 151 points, but the Central California Valley ecoregion contained only two data points. The overall distribution of the data points was representative except for the spatial gaps observed in some parts of the U.S., especially the central and western regions. The details of data distribution per the ecoregion is provided in Appendix A Table A1.
Figure 4 and Figure 5 present the spatial distribution of all D50 data following the classification schemes suggested by [43,44] respectively. The former categorizes unconsolidated sediments into six classes: silt (<0.0625 mm), sand (0.0625–2 mm), fine gravel (2–16 mm), coarse gravel (16–64 mm), cobble (64–512 mm), and boulder (>512 mm). The spatial map (Figure 4) depicts the distribution of all classes across the U.S. with no geographical bias. The map classifies the streambed material based on their susceptibility to degradation and armoring: labile bed (<16 mm), intermediate bed (16–128 mm), and coarse or armored bed (>128 mm). The labile bed class has the least resistance to erosion, whereas the coarse or armored bed has the highest. A box plot analysis indicated that the interquartile range for all 72 ecoregions lies either within the intermediate or labile bed class. This interpretation is supported by the spatial map (Figure 5) which shows the dominance of those two streambed classes. The two maps illustrate the streambed material composition and erosion susceptibility of streams based on location across the contiguous U.S.
Physiographic Provinces. Fenneman and Johnson [40] devised a method that divides the contiguous U.S. into twenty-five Physiographic Provinces. These provinces assume overall simplicity in geologic and geomorphic conditions, creating an ideal platform for the gross characterization of stream reaches. Figure 6 shows the distribution of D50 over the Physiographic Provinces. The Central Lowland province contained the highest number of data points (n = 312) whereas the St. Lawrence Valley contained the lowest number (n = 6). The D50 data were further analyzed for streambed material composition using the sediment classification scheme suggested by [43]. The percentage of silt, sand, gravel, cobbles, and boulders was calculated based on the D50 data in each of the 25 Physiographic Provinces. Analysis of data presented in Table 3 and Figure 7 allowed the identification of dominant bed materials within the divisions. The analysis indicates that gravel comprises the highest percentage of streambeds in the majority of the provinces. The Appalachian Highlands and Rocky Mountains are dominated by gravel and cobble-bedded channels. Sand and gravel were found to be abundant in the beds of stream channels located in the Atlantic Plains, Intermontane Plateaus and Pacific Mountains. Furthermore, the Interior Highlands regions saw an abundance of gravel while the Interior plains exhibited a mixture of sand, gravel and cobble in the beds of their stream channels.
Hydrologic Landscape Region (HLR). Using an integrated GIS technique and statistical tool, [35] subdivided the contiguous U.S. into twenty non-contiguous HLRs. This categorization was based on similarities in land-surface form, geologic texture, and climate characteristics. Unlike the two other regional frameworks, a single HLR can exist in different parts of the U.S. [45]. Figure 8 shows the distribution of D50 over the HLRs. All twenty HLRs are represented in the dataset. HLR 18, semiarid mountains with permeable soils and impermeable bedrock, contained the highest number of data points with a total of 352 points, whereas HLR 5 and HLR 10 contained the lowest, only 30 data points each. Box plot (Figure 9) analysis indicated HLRs with humid climate conditions (HLR 2, 4, 7, 9, 11 and 16) have higher average D50 values compared to regions with arid and sub-humid conditions. Additionally, HLRs with permeable soil conditions have higher average D50 values than those with impermeable soil.

4. Analysis of D50 Data over the Contiguous U.S.

Based on the D50 data distribution across the Contiguous U.S., a series of spatial maps of D50 were developed. Figure 10 shows an interpolated map of D50 values in mm, which was achieved through the implementation of a function in ArcMap using the technique of inverse distance weighted interpolation. This method calculates the raster value for the map by taking the average of nearby weighted observations where the weight is dependent on the distance between the observation and the raster cell. Since the D50 data over such a large geographic region exhibits a broad spatial variation, the D50 values were converted into logarithmic values (log10D50) for the spatial interpolation, which reduced the impact of large variations between the values. The logarithmic values were then converted back to mm for interpretation and visualization.
As can be inferred from the map, streams located in the central and southern parts of the U.S. are dominated by bed material composed of sand and fine gravel. On the other hand, the northeastern part of the U.S. is dominated by coarse gravel and cobble. The map shows some degree of resemblance with the smoothed D50 map developed by Kiel [27], which used 400 D50 data points in his analysis of river shaping factors and their effect on groundwater quality. Both maps depicted the coarse nature of streambeds located in the northeastern part, but there are also some differences between the two maps. A coarse material majorly represented the southeastern part of the U.S. on the smoothed map developed by [27]; contrastingly, the region was designated by fine materials on the map shown in Figure 10. Part of the Midwest U.S. that was depicted as coarse bedded in Figure 10 was largely portrayed as being fine bedded on the map developed by [27]. These discrepancies could be attributed to the number of data points used in the map development. Though there still are spatial gaps in our database of over 2200 data points, the map developed here presents a more comprehensive and reliable representation of the streambed material compositions of streams located in the U.S.
We employed a simplified method of [46] to develop a map (Figure 11) that delineates the sediment types. The interpolated log10D50 value was converted back to “mm” and classified into the three streambed groups. The majority of the U.S. falls into the first class, labile bed width sand-dominated gravel. In accordance with Figure 10, this map also captures the coarse nature of stream channels located in the northeastern part of the U.S. The map depicts a more simplified representation of streambed material distribution that serves to identify the erosion susceptibility of stream channels in different geographical regions. Similarly, Figure 12 presented the interpolated D50 map overlaid by Physiographic Provinces and Level III Ecoregions of the U.S.

5. Summary and Conclusions

This study presents a comprehensive dataset that contains information on the median grain size (D50) of the streambed and other hydrological parameters. The data were collected from five publications and spanned the contiguous U.S. Regional frameworks (Physiographic Provinces, Level III Ecoregions, and HLRs) that were used in this study to stratify the data and identify underlying trends exhibited by the streambed material. We found that Gravel is the dominant streambed material in most Physiographic Provinces. HLRs with a humid climate, permeable soil, and plateaus exhibited higher average D50 values than regions with other climates, geologic structures, and dominant landforms. Interpolated maps of D50 indicated that the central and southern parts of the U.S. are dominated by streambed material composed of sand and fine gravel, whereas the northeastern part of the U.S. is dominated by coarse gravel and cobbles. The data presented in this study can serve as a guide for future hydrological studies that rely on D50 and encompass large study areas. For application in localized scales, caution should be taken as local variations in channel slope, geology, dams, and tributary input could create site-specific conditions that deviate widely from the D50 spatial maps presented in this study. In parameterizing the streambed for a hydraulic or hydrologic model, these findings provide the user with an expectation of the regional conditions, but actual parametrization would likely be more accurately achieved with a method that accounts for the local channel slope [47].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w14213378/s1. Table S1. Sources of D50 data. Table S2. D50 data with information on latitude, longitude, and value for each.

Author Contributions

Conceptualization, J.G.A., P.M.A., M.K.J.; methodology, P.M.A.; software, D.M.A.; formal analysis, D.M.A., M.K.J.; investigation, M.K.J., P.M.A.; resources, M.J.W., J.G.A.; data curation, D.M.A.; writing—original draft preparation, D.M.A., M.K.J.; writing—review and editing, M.J.W., K.B., J.G.A., P.M.A., M.K.J.; supervision, M.K.J.; funding acquisition, J.G.A., M.J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research is based upon work that is supported by the Agricultural Research Service (ARS), U.S. Department of Agriculture (USDA), under the agreement number 58-3098-8-003. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s) and do not necessarily reflect the view of the USDA. USDA is an equal opportunity provider and employer.

Data Availability Statement

The data used in this study are available in the supplementary materials.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. D50 data distribution per Level III Ecoregions.
Table A1. D50 data distribution per Level III Ecoregions.
Level III EcoregionNo. of PointsMinQ1MedianQ3Max
Acadian Plains and Hills375.2406218.7739.3950.8287.90
Arizona/New Mexico Mountains460.007750.350.3532.00126.49
Arizona/New Mexico Plateau160.007750.350.3512.6632.00
Arkansas Valley40.3464124.0932.0032.0032.00
Atlantic Coastal Pine Barrens65.6569032.0032.00126.49126.49
Blue Mountains50.3464124.0943.0072.12126.49
Blue Ridge220.007755.6632.00126.49126.49
Boston Mountains60.007751.425.665.6632.00
Canadian Rockies110.007753.005.6662.00126.49
Cascades70.007750.3516.1770.5090.00
Central Appalachians500.0077529.0067.20115.82216.41
Central Basin and Range150.007750.015.6632.00126.49
Central California Valley232.0055.6279.25102.87126.49
Central Corn Belt Plains190.007753.0032.0079.25126.49
Central Great Plains300.130000.9632.00126.49126.49
Central Irregular Plains180.007750.505.6679.25126.49
Coast Range710.007751.6732.0053.30200.00
Colorado Plateaus350.007750.482.135.85165.00
Driftless Area130.007755.6632.0032.00126.49
Eastern Cascades Slopes and Foothills432.0032.0079.25126.49126.49
Eastern Corn Belt Plains330.007757.3722.8040.66126.49
Erie Drift Plains290.0077510.1228.0452.30126.49
Huron/Erie Lake Plains60.007750.284.8032.00126.49
Idaho Batholith950.007758.3337.5091.00190.00
Interior Plateau600.0077528.7537.9778.53194.00
Lake Agassiz Plain250.007750.350.875.66126.49
Madrean Archipelago70.007750.353.0025.4132.00
Middle Atlantic Coastal Plain190.007750.331.0021.50126.49
Middle Rockies850.007750.3519.0040.00268.00
Mississippi Alluvial Plain160.007751.6718.8332.00126.49
North Central Appalachians240.0077543.4081.90119.00133.00
North Central Hardwood Forests110.007755.6632.00126.49126.49
Northeastern Coastal Zone250.0077510.7332.0086.26171.00
Northeastern Highlands910.007755.6632.0084.96179.45
Northern Allegheny Plateau420.6520.7051.3068.20116.00
Northern Basin and Range210.007750.3532.0034.00126.49
Northern Glaciated Plains150.007750.350.3532.00126.49
Northern Lakes and Forests380.007750.3532.0032.95126.49
Northern Piedmont670.007751.1314.0032.00132.81
Northern Rockies410.007755.3332.0087.50315.00
Northwestern Glaciated Plains110.007755.6632.00126.49126.49
Northwestern Great Plains620.007750.355.6632.00178.89
Ouachita Mountains100.0077532.0032.0036.20126.49
Ozark Highlands150.007750.355.6632.00126.49
Piedmont570.007755.6632.00126.49126.49
Ridge and Valley1510.0077532.0050.6091.44222.00
S. Michigan/N. Indiana Drift Plains110.007750.3532.00126.49126.49
Sierra Nevada220.007750.3532.00126.49126.49
Sonoran Basin and Range60.346410.3516.1755.62126.49
South Central Plains340.007755.6632.00126.49126.49
Southeastern Plains580.007750.351.119.98126.49
Southeastern Wisconsin Till Plains90.007750.1832.00126.49126.49
Southern and Central California Chaparral and Oak320.007750.355.66126.49126.49
Southern California Mountains470.007750.350.3532.00126.49
Southern California/Northern Baja Coast570.125000.903.3032.93151.80
Southern Coastal Plain150.007750.590.715.66126.49
Southern Rockies990.007755.6632.0076.00190.00
Southwestern Appalachians230.007755.6628.0046.30135.00
Southwestern Tablelands260.007750.010.351.67126.49
Wasatch And Uinta Mountains60.007750.938.0827.36126.49
Western Allegheny Plateau620.007755.6631.2354.12131.06
Western Corn Belt Plains440.007755.6632.00126.49126.49
Western Gulf Coastal Plain120.007751.1418.8355.62126.49
Western High Plains120.007750.353.0032.0032.00
Willamette Valley645.0048.0054.0060.0066.00
Wyoming Basin240.0077534.0068.55102.88210.00

References

  1. Lane, E.W. The Importance of Fluvial Morphology in Hydraulic Engineering. In Proceedings of the American Society of Civil Engineers, San Diego, CA, USA, 22–26 October 1995. [Google Scholar]
  2. Ritter, J.B.; Rumschlag, J.H.; Zaleha, M.J. Evaluating recent stream channel and pattern changes for stream resource protection and restoration: An example from west-central Ohio. J. Great Lakes Res. 2007, 33, 154–166. [Google Scholar] [CrossRef]
  3. Rosgen, D.L. A Stream Channel Stability Assessment Methodology. In Proceedings of the Seventh Federal Interagency Sedimentation Conference, Reno, Nevada, 25–29 March 2001. [Google Scholar]
  4. Simon, A.; Rinaldi, M. Channel instability in the loess area of the midwestern United States. J. Am. Water Res. Assoc. 2000, 36, 133–150. [Google Scholar] [CrossRef]
  5. Morét, S. Predicting Channel Stability in Colorado Mountain Streams using Hydrobiogeomorphic and Land Use Data: A Cost-sensitive Machine Learning Approach to Modeling Rapid Assessment Protocols. Ph.D. Thesis, Oregon State University, Corvallis, OR, USA, 16 March 2001. [Google Scholar]
  6. Johnson, P.A. Assessing Stream Channel Stability at Bridges in Physiographic Regions. A Report of the U.S. Department of Transportation. 2006. Available online: https://www.fhwa.dot.gov/publications/research/infrastructure/hydraulics/05072 (accessed on 1 August 2022).
  7. Lagasse, P.F.; Schall, J.D.; Johnson, F.; Richardson, E.V.; Chang, F. Stream Stability at Highway Structures. A Report of the U.S Department of Transportation. 1995. Available online: https://rosap.ntl.bts.gov/view/dot/25666 (accessed on 1 August 2022).
  8. Lammers, R.W.; Bledsoe, B.P. Parsimonious sediment transport equations based on Bagnold’s stream power approach. Earth Surf. Process. Landf. 2018, 43, 242–258. [Google Scholar] [CrossRef]
  9. Schwendel, A.C.; Death, R.G.; Fuller, I.C. The assessment of shear stress and bed stability in stream ecology. Freshw. Biol. 2010, 55, 261–281. [Google Scholar] [CrossRef] [Green Version]
  10. Lorang, M.S.; Hauer, F.R. Flow competence and streambed stability: An evaluation of technique and application. J. N. Am. Benthol. Soc. 2003, 22, 475–491. [Google Scholar] [CrossRef]
  11. Kaufmann, P.R.; Larsen, D.P.; Faustini, J.M. Bed stability and sedimentation associated with human disturbances in Pacific Northwest streams. J. Am. Water Res. Assoc. 2009, 45, 434–459. [Google Scholar] [CrossRef]
  12. Afzalimehr, H.; Abdolhosseini, M.; Singh, V.P. Hydraulic geometry relations for stable channel design. J. Hydrol. Eng. 2010, 15, 859–864. [Google Scholar] [CrossRef]
  13. Andrews, E. Bed-material entrainment and hydraulic geometry of gravel-bed rivers in Colorado. Geol. Soc. Am. Bull. 1984, 95, 371. [Google Scholar] [CrossRef]
  14. MacKenzie, L.G.; Eaton, B.C.; Church, M. Breaking from the average: Why large grains matter in gravel-bed streams. Earth Surf. Process. Landf. 2018, 43, 3190–3196. [Google Scholar] [CrossRef]
  15. Pfeiffer, A.; Finnegan, N.; Willenbring, J. Sediment supply controls equilibrium channel geometry in gravel rivers. Proc. Natl. Acad. Sci. USA 2017, 114, 3346–3351. [Google Scholar] [CrossRef]
  16. Bunte, K.; Abt, S.R. Sampling Surface and Subsurface Particle-Size Distributions in Wadable Gravel- and Cobble-Bed Streams for Analyses in Sediment Transport, Hydraulics, and Streambed Monitoring; US Department of Agriculture, Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2001. [Google Scholar] [CrossRef]
  17. Parker, G.; Wilcock, P.R.; Paola, C.; Dietrich, W.E.; Pitlick, J. Physical basis for quasi-universal relations describing bankfull hydraulic geometry of single-thread gravel bed rivers. J. Geophys. Res. 2007, 112. [Google Scholar] [CrossRef]
  18. Wilkerson, G.V.; Parker, G. Physical basis for quasi-universal relationships describing bankfull hydraulic geometry of sand-bed rivers. J. Hydrol. Eng. 2011, 137, 739–753. [Google Scholar] [CrossRef]
  19. Golden, L.A.; Springer, G.S. Channel geometry, median grain size, and stream power in small mountain streams. Geomorphology 2006, 78, 64–76. [Google Scholar] [CrossRef]
  20. Hack, J.T. Studies of Longitudinal Stream Profiles in Virginia and Maryland; Geological Survey Professional Paper, 245-B; U.S. Government Publishing Office: Washington, DC, USA, 1957. [Google Scholar]
  21. Bagnold, R.A. An Empirical Correlation of Bedload Transport Rates in Flumes and Natural Rivers. Proc. R. Soc. Lond. Math. Phys. Sci. 1980, 372, 453–473. [Google Scholar]
  22. Camenen, B. Discussion of “Understanding the influence of slope on the threshold of coarse grain motion: Revisiting critical stream power by C. Parker, N.J. Clifford, and C.R. Thorne. Geomorphology 2012, 126, 51–65. [Google Scholar] [CrossRef] [Green Version]
  23. Ferguson, R.I. Hydraulics and hydraulic geometry. Prog. Phys. Geogr. Earth Environ. 1986, 10, 1–31. [Google Scholar] [CrossRef]
  24. Parker, C.; Clifford, N.J.; Thorne, C.R. Understanding the influence of slope on the threshold of coarse grain motion: Revisiting critical stream power. Geomorphology 2011, 126, 51–65. [Google Scholar] [CrossRef]
  25. Parker, G.; Klingeman, P.C.; McLean, D.G. Bedload and size distribution in paved gravel-bed streams. J. Hydraul. Div.—ASCE 1982, 108, 544–571. [Google Scholar] [CrossRef]
  26. Wilcock, P.R.; Crowe, J.C. Surface-based transport model for mixed-size sediment. J. Hydraul. Eng. 2003, 129, 120–128. [Google Scholar] [CrossRef]
  27. Kiel, B.A. Measurements of U.S. Rivers Clarify River-Shaping Factors and Interaction with Groundwater. Ph.D. Thesis, University of Texas at Austin, Austin, TX, USA, May 2015. [Google Scholar]
  28. Bieger, K.; Rathjens, H.; Allen, P.M.; Arnold, J.G. Development and evaluation of bankfull hydraulic geometry relationships for the physiographic regions of the United States. J. Am. Water Res. Assoc. 2015, 51, 842–858. [Google Scholar] [CrossRef]
  29. Blackburn-Lynch, W.; Agouridis, C.T.; Barton, C.D. Development of regional curves for hydrologic landscape regions (HLR) in the contiguous United States. J. Am. Water Res. Assoc. 2017, 53, 903–928. [Google Scholar] [CrossRef]
  30. Faustini, J.M.; Kaufmann, P.R.; Herlihy, A.T. Downstream variation in bankfull width of wadeable streams across the conterminous United States. Geomorphology 2009, 108, 292–311. [Google Scholar] [CrossRef] [Green Version]
  31. Johnson, P.A.; Fecko, B.J. Regional channel geometry equations: A statistical comparison for physiographic provinces in the eastern US. River Res. Appl. 2008, 24, 823–834. [Google Scholar] [CrossRef]
  32. Simon, A.; Dickerson, W.; Heins, A. Suspended-sediment transport rates at the 1.5-year recurrence interval for ecoregions of the United States: Transport conditions at the bankfull and effective discharge? Geomorphology 2004, 58, 243–262. [Google Scholar] [CrossRef] [Green Version]
  33. USEPA. Wadeable Stream Assessment: A Collaborative Summary of the Nation’s Streams. A Report of the United States Environmental Protection Agency. 2006. Available online: https://www.epa.gov/sites/production/files/2014-10/documents/2007_5_16_streamsurvey_wsa_assessment_may2007.pdf (accessed on 1 August 2022).
  34. Bledsoe, B.; Baker, D.; Nelson, P.; Rosburg, T.; Sholtes, J.; Stroth, T. Design Hydrology for Stream Restoration and Channel Stability at Stream Crossings; Project Report for the Transportation Research Board of The National Academics of Sciences, Engineering, and Medicine; Colorado State University: Fort Collins, CO, USA, 2016. [Google Scholar]
  35. Hawley, R.J.; Bledsoe, B.P. Channel enlargement in semiarid suburbanizing watersheds: A southern California case study. J. Hydrol. 2013, 496, 17–30. [Google Scholar] [CrossRef]
  36. Slater, L.J.; Singer, M.B. Imprint of climate and climate change in alluvial riverbeds: Continental United States, 1950–2011. Geology 2013, 41, 595–598. [Google Scholar] [CrossRef] [Green Version]
  37. Vic, B.; Toby, L. Outliers in Statistical Data, 3rd ed.; Wiley: Hoboken, NJ, USA, 1994; ISBN 978-0-471-93094-5. [Google Scholar]
  38. Wolock, D.M.; Winter, T.C.; McMahon, G. Delineation and evaluation of hydrologic-landscape regions in the United States using geographic information system tools and multivariate statistical analyses. Environ. Manag. 2004, 34, S71–S88. [Google Scholar] [CrossRef] [PubMed]
  39. Omernik, J.M.; Griffith, G.E. Ecoregions of the conterminous United States: Evolution of a hierarchical spatial framework. Environ. Manag. 2014, 54, 1249–1266. [Google Scholar] [CrossRef] [PubMed]
  40. Fenneman, N.M.; Johnson, D.W. Physiographic Divisions of the Conterminous, U.S.; U.S. Geological Survey Special Map Series, Scale 1:7,000,000; U.S. Geological Survey: Reston, VA, USA, 1946. Available online: http://water.usgs.gov/GIS/metadata/usgswrd/XML/physio.xml (accessed on 1 August 2022).
  41. Wright, R.G.; Murray, M.P.; Merrill, T. Ecoregions as a level of ecological analysis. Biol. Conserv. 1998, 86, 207–213. [Google Scholar] [CrossRef]
  42. Castro, J.M.; Jackson, P.L. Bankfull discharge recurrence intervals and regional hydraulic geometry relationships: Patterns in the Pacific Northwest, USA. J. Am. Water Res. Assoc. 2001, 37, 1249–1262. [Google Scholar] [CrossRef]
  43. Konrad, C.; Gellis, A. Factors influencing fine sediment on stream beds in the midwestern United States. J. Environ. Qual. 2018, 47, 1214–1222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Bledsoe, B.P.; Hawley, R.J.; Booth, D.B. Hydromodification Screening Tools: Field Manual for Assessing Channel Susceptibility; Colorado State University: Fort Collins, CO, USA, 2010; Available online: http://ftp.sccwrp.org/pub/download/DOCUMENTS/TechnicalReports/606_HydromodScreeningTools_FieldManual.pdf (accessed on 1 August 2022).
  45. Santhi, C.; Allen, P.M.; Muttiah, R.S.; Arnold, J.G.; Tuppad, P. Regional estimation of base flow for the conterminous United States by hydrologic landscape regions. J. Hydrol. 2008, 351, 139–153. [Google Scholar] [CrossRef]
  46. Ferguson, R.I. River channel slope, flow resistance, and gravel entrainment thresholds. Water Resour. Res. 2012, 48, W05517. [Google Scholar] [CrossRef]
  47. Snyder, N.P.; Nesheim, A.O.; Wilkins, B.C.; Edmonds, D.A. Predicting grain size in gravel-bedded rivers using digital elevation models: Application to three Maine watersheds. Geol. Soc. Amer. Bull. 2013, 125, 148–163. [Google Scholar] [CrossRef]
Figure 1. Distribution of D50 data points from five sources across the contiguous U.S. Dots with different colors in the figure represent different D50 data sources.
Figure 1. Distribution of D50 data points from five sources across the contiguous U.S. Dots with different colors in the figure represent different D50 data sources.
Water 14 03378 g001
Figure 2. Box plots using D50 data values from five datasets: (a) original data, (b) after removing outliers. The outliers are designated by shaded circles, which are the values beyond the whiskers (lines extended up to 1.5 times the interquartile range) in both directions from the box. The box represents the range of 25th and 75th percentile. The horizontal line within the box and the ‘x’ sign represent median and average values respectively. Outliers are designated by shaded circles.
Figure 2. Box plots using D50 data values from five datasets: (a) original data, (b) after removing outliers. The outliers are designated by shaded circles, which are the values beyond the whiskers (lines extended up to 1.5 times the interquartile range) in both directions from the box. The box represents the range of 25th and 75th percentile. The horizontal line within the box and the ‘x’ sign represent median and average values respectively. Outliers are designated by shaded circles.
Water 14 03378 g002
Figure 3. Distribution of D50 data points over Level III Ecoregions. Dots with different colors in the figure represent different D50 data sources.
Figure 3. Distribution of D50 data points over Level III Ecoregions. Dots with different colors in the figure represent different D50 data sources.
Water 14 03378 g003
Figure 4. Distribution of D50 data points over Level III Ecoregions.
Figure 4. Distribution of D50 data points over Level III Ecoregions.
Water 14 03378 g004
Figure 5. Streambed D50 values over Level III Ecoregions based on the channel susceptibility scheme.
Figure 5. Streambed D50 values over Level III Ecoregions based on the channel susceptibility scheme.
Water 14 03378 g005
Figure 6. Distribution of D50 data points over Physiographic Provinces.
Figure 6. Distribution of D50 data points over Physiographic Provinces.
Water 14 03378 g006
Figure 7. Box Plot showing the range of D50 data over Physiographic Provinces of the U.S.
Figure 7. Box Plot showing the range of D50 data over Physiographic Provinces of the U.S.
Water 14 03378 g007
Figure 8. Box Plot showing the range of D50 data over Hydrologic Landscape Regions of the U.S. Boxes indicate data range falling between the first and third quartile. Note: Horizontal lines within the boxes indicate median values whereas the ‘X’ sign designates averages.
Figure 8. Box Plot showing the range of D50 data over Hydrologic Landscape Regions of the U.S. Boxes indicate data range falling between the first and third quartile. Note: Horizontal lines within the boxes indicate median values whereas the ‘X’ sign designates averages.
Water 14 03378 g008
Figure 9. Distribution of D50 data points over HLR. Colors on the background map represent HLRs.
Figure 9. Distribution of D50 data points over HLR. Colors on the background map represent HLRs.
Water 14 03378 g009
Figure 10. Map of interpolated D50 over the Contiguous U.S. Note: The map was prepared using the log value of more than 2200 D50 data points collected from five sources; black dots on the map show the locations of D50 data points.
Figure 10. Map of interpolated D50 over the Contiguous U.S. Note: The map was prepared using the log value of more than 2200 D50 data points collected from five sources; black dots on the map show the locations of D50 data points.
Water 14 03378 g010
Figure 11. Map of interpolated D50 over the Contiguous U.S. using broad categories of streambed material. Note: The map was prepared using the log value of more than 2200 D50 data points collected from five sources; black dots on the map show the locations of D50 data points.
Figure 11. Map of interpolated D50 over the Contiguous U.S. using broad categories of streambed material. Note: The map was prepared using the log value of more than 2200 D50 data points collected from five sources; black dots on the map show the locations of D50 data points.
Water 14 03378 g011
Figure 12. Map of interpolated D50 over the Contiguous U.S. with three streambed classes (a) overlaid with a map of physiographic provinces (b) overlaid with a map of level III Ecoregions. The black dots on the map show the locations of D50 data points.
Figure 12. Map of interpolated D50 over the Contiguous U.S. with three streambed classes (a) overlaid with a map of physiographic provinces (b) overlaid with a map of level III Ecoregions. The black dots on the map show the locations of D50 data points.
Water 14 03378 g012aWater 14 03378 g012b
Table 1. D50 data sources with other hydrological parameters.
Table 1. D50 data sources with other hydrological parameters.
Data SourceHydrological ParameterNo. of Data Points
[28] Bieger et al. (2015)Q, D50, W, D, S, and DA642
[33] USEPA (2016)D501387
[34] Bledsoe et al. (2016)Q, D50, and DA103
[35] Hawley & Bledsoe (2013)D5066
[36] Slater & Singer (2013)D50, W, D, S, and DA255
Note: Q: bankfull discharge; D50: median grain size; W: bankfull width; D: bankfull depth; S: Slope; DA: Drainage area.
Table 2. Descriptive Statistics for the D50 data obtained from the five sources.
Table 2. Descriptive Statistics for the D50 data obtained from the five sources.
Data Characteristics[28] Bieger et al., 2015[33] USEPA 2006[34] Bledsoe et al., 2016[35] Hawley & Bledsoe 2013[36] Slater & Singer 2013
Total Number of Data60612029659245
Maximum (mm)152178.913469.7163.1
Minimum (mm))0.060.0080.060.1255.24
Q1 (mm)10.380.350.560.9031.4
Median (mm)355.665.2152.6051.6
Q3 (mm)60.93326314.5080.00
Mean (mm)56.4335.9231.2912.0260.29
Range151.94178.88133.9469.58157.83
Number of outliers68461033
Table 3. Bed material percentage composition of physiographic provinces of the U.S.
Table 3. Bed material percentage composition of physiographic provinces of the U.S.
No.Physiographic ProvinceBed Composition Based on D50 Data (%)No. of D50 Points
SiltSandGravelCobbleBoulder
1Superior Upland8.622.940.011.417.135
2Coastal Plain10.734.833.118.03.4171
3Piedmont5.619.054.817.53.2121
4Blue Ridge province7.715.438.534.63.826
5Valley and Ridge province3.23.857.630.45.1159
6St. Lawrence Valley16.70.050.033.30.06
7Appalachian Plateaus province5.15.850.034.34.7275
8New England Province8.93.063.025.20.0134
9Adirondack province7.715.430.846.20.014
10Interior Low Plateaus4.25.661.127.81.472
11Central Lowland9.722.837.823.46.3312
12Great Plains18.631.031.714.54.1138
13Ozark Plateaus20.015.045.015.05.020
14Ouachita province5.65.672.211.15.618
15Southern Rocky Mountains7.815.646.727.82.290
16Wyoming Basin6.93.448.341.40.031
17Middle Rocky Mountains20.76.924.144.83.420
18Northern Rocky Mountains9.311.749.227.82.0248
19Columbia Plateau12.50.037.531.318.816
20Colorado Plateaus7.540.343.37.51.570
21Basin and Range Province15.235.431.311.17.197
22Cascade-Sierra Mountains15.613.342.226.72.245
23Pacific Border province9.828.744.414.22.9272
24Lower California province15.035.050.00.00.020
25Continental shelf------
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jha, M.K.; Asamen, D.M.; Allen, P.M.; Arnold, J.G.; White, M.J.; Bieger, K. Streambed Median Grain Size (D50) across the Contiguous U.S. Water 2022, 14, 3378. https://doi.org/10.3390/w14213378

AMA Style

Jha MK, Asamen DM, Allen PM, Arnold JG, White MJ, Bieger K. Streambed Median Grain Size (D50) across the Contiguous U.S. Water. 2022; 14(21):3378. https://doi.org/10.3390/w14213378

Chicago/Turabian Style

Jha, Manoj K., Dawit M. Asamen, Peter M. Allen, Jeffrey G. Arnold, Michael J. White, and Katrin Bieger. 2022. "Streambed Median Grain Size (D50) across the Contiguous U.S." Water 14, no. 21: 3378. https://doi.org/10.3390/w14213378

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