Drainage density: Difference between revisions
the article is listed in the navbox per WP:NAV-WITHIN |
Open access status updates in citations with OAbot #oabot |
||
(14 intermediate revisions by 11 users not shown) | |||
Line 1: | Line 1: | ||
{{Short description|Physical parameter of a drainage basin}} |
|||
'''Drainage density''' is a quantity used to describe physical parameters of a [[drainage basin]]. First described by Robert E. Horton, drainage density is defined as the total length of channel in a drainage basin divided by the total area, represented by the following equation: |
'''Drainage density''' is a quantity used to describe physical parameters of a [[drainage basin]]. First described by [[Robert E. Horton]], drainage density is defined as the total length of channel in a drainage basin divided by the total area, represented by the following equation: |
||
<center><math>D_{d}=\frac{\sum{L}}{A_{basin}}</math></ |
<div class="center"><math>D_{d}=\frac{\sum{L}}{A_{basin}}</math></div><ref name="Horton1932">{{cite journal |last1=Horton |first1=Robert E. |title=Drainage-basin characteristics |journal=Transactions of the American Geophysical Union |date=June 1932 |volume=13 |issue=1 |pages=350–361 |doi=10.1029/TR013i001p00350|bibcode=1932TrAGU..13..350H }}</ref> |
||
The quantity represents the average length of channel per unit area of [[catchment]] and has units <math>\frac{\left[L\right]}{\left[L^2\right]}</math>, which is often reduced to <math>\left[L^{-1}\right]</math>. |
The quantity represents the average length of channel per unit area of [[catchment]] and has units <math>\frac{\left[L\right]}{\left[L^2\right]}</math>, which is often reduced to <math>\left[L^{-1}\right]</math>. |
||
Drainage density depends upon both [[climate]] and physical characteristics of the drainage basin. Soil [[Permeability (earth sciences)|permeability (infiltration difficulty)]] and underlying rock type affect the runoff in a watershed; [[Permeability (earth sciences)|impermeable]] ground or exposed [[bedrock]] will lead to an increase in [[Surface runoff|surface water runoff]] and therefore to more frequent streams. Rugged regions or those with high [[relief]] will also have a higher drainage density than other drainage basins if the other characteristics of the basin are the same. |
Drainage density depends upon both [[climate]] and physical characteristics of the drainage basin. Soil [[Permeability (earth sciences)|permeability (infiltration difficulty)]] and underlying rock type affect the runoff in a watershed; [[Permeability (earth sciences)|impermeable]] ground or exposed [[bedrock]] will lead to an increase in [[Surface runoff|surface water runoff]] and therefore to more frequent streams. Rugged regions or those with high [[Terrain#Relief|relief]] will also have a higher drainage density than other drainage basins if the other characteristics of the basin are the same. |
||
When determining the total length of [[stream]]s in a basin, both [[perennial stream|perennial]] and [[ephemeral stream]]s should be considered.<ref name="Horton1945">{{cite journal |last1=Horton |first1=Robert E. |title=Erosional Development of Streams and Their Drainage Basins; Hydrophysical Approach to Quantitative Morphology |journal=GSA Bulletin |date=1945-03-01 |volume=56 |issue=3 |pages=275–370 |doi=10.1130/0016-7606(1945)56[275:EDOSAT]2.0.CO;2|url=http://revele.uncoma.edu.ar/htdoc/revele/index.php/geografia/article/view/1761 }}</ref> If a drainage basin contained only ephemeral streams, the drainage density by the equation above would be calculated to be zero if only the total length of streams was calculated using only perennial streams. Ignoring ephemeral streams in the calculations does not consider the behavior of the basin during flood events and is therefore not completely representative of the drainage characteristics of the basin. |
When determining the total length of [[stream]]s in a basin, both [[perennial stream|perennial]] and [[ephemeral stream]]s should be considered.<ref name="Horton1945">{{cite journal |last1=Horton |first1=Robert E. |title=Erosional Development of Streams and Their Drainage Basins; Hydrophysical Approach to Quantitative Morphology |journal=GSA Bulletin |date=1945-03-01 |volume=56 |issue=3 |pages=275–370 |doi=10.1130/0016-7606(1945)56[275:EDOSAT]2.0.CO;2|bibcode=1945GSAB...56..275H |s2cid=129509551 |url=http://revele.uncoma.edu.ar/htdoc/revele/index.php/geografia/article/view/1761 }}</ref> If a drainage basin contained only ephemeral streams, the drainage density by the equation above would be calculated to be zero if only the total length of streams was calculated using only perennial streams. Ignoring ephemeral streams in the calculations does not consider the behavior of the basin during flood events and is therefore not completely representative of the drainage characteristics of the basin. |
||
Drainage density is indicative of [[Infiltration (hydrology)|infiltration]] and [[Hydraulic permeability|permeability]] of a drainage basin, as well as relating to the shape of the [[hydrograph]]. Drainage density depends upon both [[climate]] and physical characteristics of the drainage basin. |
Drainage density is indicative of [[Infiltration (hydrology)|infiltration]] and [[Hydraulic permeability|permeability]] of a drainage basin, as well as relating to the shape of the [[hydrograph]]. Drainage density depends upon both [[climate]] and physical characteristics of the drainage basin. |
||
Line 13: | Line 14: | ||
== Inverse of drainage density as a physical quantity == |
== Inverse of drainage density as a physical quantity == |
||
Drainage density can be used to approximate the average length of [[overland flow]] in a catchment. Horton (1945) used the following equation to describe the average length of overland flow as a function of drainage density:<ref name="Horton1945" /> |
Drainage density can be used to approximate the average length of [[overland flow]] in a catchment. Horton (1945) used the following equation to describe the average length of overland flow as a function of drainage density:<ref name="Horton1945" /> |
||
<center><math>l_O=\frac{1}{2D_d}</math></ |
<div class="center"><math>l_O=\frac{1}{2D_d}</math></div> |
||
Where <math>l_0</math> is the length of overland flow with units of length and <math>D_d</math> is the drainage density of the catchment, expressed in units of inverse length. |
Where <math>l_0</math> is the length of overland flow with units of length and <math>D_d</math> is the drainage density of the catchment, expressed in units of inverse length. |
||
Considering the geometry of channels on the hillslope, Horton also proposed the following equation: |
Considering the geometry of channels on the hillslope, Horton also proposed the following equation: |
||
<center><math>l_O=\frac{1}{2D_d\sqrt{1-\left(\frac{s_c}{s_g}\right)^2}}</math></ |
<div class="center"><math>l_O=\frac{1}{2D_d\sqrt{1-\left(\frac{s_c}{s_g}\right)^2}}</math></div><ref name="Horton1945" /> |
||
Where <math>s_c</math> is the channel slope and <math>s_g</math> is the average slope of the ground in the area. |
Where <math>s_c</math> is the channel slope and <math>s_g</math> is the average slope of the ground in the area. |
||
== Elementary components of drainage basins == |
== Elementary components of drainage basins == |
||
A drainage basin can be defined by three elementary quantities: channels, the hillslope area associated with those channels, and the source areas.<ref name="MontgomeryDietrich" /> The channels are the well-defined segments that efficiently carry water through the catchment. Labeling these features as “channels” rather than “streams” indicates that there need not be a continuous flow of water to capture the behavior of this region as a conduit of water. According to Arthur Strahler’s stream ordering system,<ref name="Strahler1957">{{cite journal |last1=Strahler |first1=Arthur N. |title=Quantitative Analysis of Watershed Geomorphology |journal=Transactions, American Geophysical Union |date=December 1957 |volume=38 |issue=6 |pages=913–920|doi=10.1029/TR038i006p00913 }}</ref> the channels are not defined to be any single order or range of orders. Channels of lower orders combine to form higher order channels. The associated hillslope areas are the hillslopes that slope directly into the channels.<ref name="MontgomeryDietrich" /> Precipitation that enters the system on the hillslopes areas and is not lost to infiltration or evapotranspiration enters the channels. The source areas are concave regions of hillslope that are associated with a single channel.<ref name="Strahler1957" /> Precipitation entering a source area that is not lost to infiltration or evapotranspiration flows through the source area and enters the channel at the channel’s head. Source areas and the hillslope areas associated with channels are differentiated by source areas draining through the channel head, while the associated hillslope areas drain into the rest of the stream.<ref name="MontgomeryDietrich" /> According to Strahler’s stream ordering system, all source areas drain into a primary channel, by the definition of a primary channel.<ref name="Strahler1957" /> |
A drainage basin can be defined by three elementary quantities: channels, the hillslope area associated with those channels, and the source areas.<ref name="MontgomeryDietrich" /> The channels are the well-defined segments that efficiently carry water through the catchment. Labeling these features as “channels” rather than “streams” indicates that there need not be a continuous flow of water to capture the behavior of this region as a conduit of water. According to Arthur Strahler’s stream ordering system,<ref name="Strahler1957">{{cite journal |last1=Strahler |first1=Arthur N. |title=Quantitative Analysis of Watershed Geomorphology |journal=Transactions, American Geophysical Union |date=December 1957 |volume=38 |issue=6 |pages=913–920|doi=10.1029/TR038i006p00913 |bibcode=1957TrAGU..38..913S }}</ref> the channels are not defined to be any single order or range of orders. Channels of lower orders combine to form higher order channels. The associated hillslope areas are the hillslopes that slope directly into the channels.<ref name="MontgomeryDietrich" /> Precipitation that enters the system on the hillslopes areas and is not lost to infiltration or evapotranspiration enters the channels. The source areas are concave regions of hillslope that are associated with a single channel.<ref name="Strahler1957" /> Precipitation entering a source area that is not lost to infiltration or evapotranspiration flows through the source area and enters the channel at the channel’s head. Source areas and the hillslope areas associated with channels are differentiated by source areas draining through the channel head, while the associated hillslope areas drain into the rest of the stream.<ref name="MontgomeryDietrich" /> According to Strahler’s stream ordering system, all source areas drain into a primary channel, by the definition of a primary channel.<ref name="Strahler1957" /> |
||
Bras et al. (1991)<ref name="PhysicalBasis">{{cite journal |last1=Tarboton |first1=David G. |last2=Bras |first2=Rafael L. |last3=Rodriguez-Iturbe |first3=Ignacio |title=A physical basis for drainage density |journal=Geomorphology |date=October 19, 1991 |volume=5 |issue=1–2 |pages=59–76 |doi=10.1016/0169-555X(92)90058-V}}</ref> describe the conditions that are necessary for channel formation. Channel formation is a concept intimately tied to the formation and evolution of a drainage system and influence the drainage density of catchment. The relation they propose determines the behavior of a given hillslope in response to a small perturbation. They propose the following equation as a relation between source area, source slope, and the sediment flux through this source area: |
Bras et al. (1991)<ref name="PhysicalBasis">{{cite journal |last1=Tarboton |first1=David G. |last2=Bras |first2=Rafael L. |last3=Rodriguez-Iturbe |first3=Ignacio |title=A physical basis for drainage density |journal=Geomorphology |date=October 19, 1991 |volume=5 |issue=1–2 |pages=59–76 |doi=10.1016/0169-555X(92)90058-V}}</ref> describe the conditions that are necessary for channel formation. Channel formation is a concept intimately tied to the formation and evolution of a drainage system and influence the drainage density of catchment. The relation they propose determines the behavior of a given hillslope in response to a small perturbation. They propose the following equation as a relation between source area, source slope, and the sediment flux through this source area: |
||
<center><math>a\frac{\partial F}{\partial S}\frac{dS}{da}=F-a\frac{\partial F}{\partial a}</math></ |
<div class="center"><math>a\frac{\partial F}{\partial S}\frac{dS}{da}=F-a\frac{\partial F}{\partial a}</math></div><ref name="PhysicalBasis" /> |
||
Where '''F''' is the sediment flux, '''S''' is the slope of the source area, and '''a''' is the source area. The right-hand side of this relation determines channel stability or instability. If the right-hand side of the equation is greater than zero, the hillslope is stable, and small perturbations such as small erosive events do no develop into channels. Conversely, if the right-hand side of the equation is less than zero, Bras et al.<ref name="PhysicalBasis" /> determine the hillslope to be unstable, and small erosive structures, such as rills, will tend to grow and form a channel and increase the drainage density of a basin. In this sense, "unstable" is not used in the sense of the gradient of the hillslope being greater than the angle of repose and therefore susceptible to mass wasting, but rather fluvial erosive processes such as sheet flow or channel flow tend to incise and erode to form a singular channel.<ref name="PhysicalBasis" /> Therefore, the characteristics of the source area, or potential source area, influence the drainage density and evolution of a drainage basin.<ref name="PhysicalBasis" /> |
Where '''F''' is the sediment flux, '''S''' is the slope of the source area, and '''a''' is the source area. The right-hand side of this relation determines channel stability or instability. If the right-hand side of the equation is greater than zero, the hillslope is stable, and small perturbations such as small erosive events do no develop into channels. Conversely, if the right-hand side of the equation is less than zero, Bras et al.<ref name="PhysicalBasis" /> determine the hillslope to be unstable, and small erosive structures, such as rills, will tend to grow and form a channel and increase the drainage density of a basin. In this sense, "unstable" is not used in the sense of the gradient of the hillslope being greater than the angle of repose and therefore susceptible to mass wasting, but rather fluvial erosive processes such as sheet flow or channel flow tend to incise and erode to form a singular channel.<ref name="PhysicalBasis" /> Therefore, the characteristics of the source area, or potential source area, influence the drainage density and evolution of a drainage basin.<ref name="PhysicalBasis" /> |
||
Line 31: | Line 32: | ||
== Relation to water balance == |
== Relation to water balance == |
||
Drainage density is tied to the water balance equation: |
Drainage density is tied to the water balance equation: |
||
<center><math>\frac{dV}{dt}=R+G_i-G_o-G_s-ET-Q_w</math></ |
<div class="center"><math>\frac{dV}{dt}=R+G_i-G_o-G_s-ET-Q_w</math></div><ref name="Fitts">{{cite book |last1=Fitts |first1=Charles R. |title=Groundwater Science |date=2013 |publisher=Elsevier Inc. |location=Waltham, MA |isbn=978-0-12-384705-8 |page=14 |edition=2}}</ref> |
||
Where <math>\frac{dV}{dt}</math> is the change in [[reservoir]] storage, '''R''' is [[precipitation]], '''ET''' is [[evapotranspiration]], '''G<sub>i</sub>''' and '''G<sub>o</sub>''' are the respective [[groundwater]] flux into and out of the basin, '''G<sub>s</sub>''' is the [[groundwater discharge]] into streams, and '''Q<sub>w</sub>''' is groundwater discharge from the basin through [[well]]s. Drainage density relates to the storage and runoff terms. Drainage density relates to the efficiency by which water is carried over the landscape. Water is carried through channels much faster than over hillslopes, as [[overland flow|saturated overland flow]] is slower due to being thinned out and obstructed by [[vegetation]] or pores in the ground.<ref name="GeomorphologyTextbook">{{cite book |last1=Anderson |first1=Robert S. |last2=Anderson |first2=Suzanne P. |title=Geomorphology: The Mechanics and Chemistry of Landscapes |date=2010 |publisher=Cambridge University Press |location=Cambridge, UK |isbn=978-0-521-51978-6}}</ref> Consequently, a drainage basin with a relatively higher drainage density will be more efficiently drained than a higher density one.<ref name="GeomorphologyTextbook" /> Because of the more extensive drainage system in a higher density basin, precipitation entering the basement will, on average, travel a shorter distance over the slower hillslopes before reaching the faster-flowing channels and exit the basin through the channels in less time. Conversely, precipitation entering a lower drainage density basin will take longer to exit the basin due to travelling over the slower hillslope longer.<ref name="GeomorphologyTextbook" /> |
Where <math>\frac{dV}{dt}</math> is the change in [[reservoir]] storage, '''R''' is [[precipitation]], '''ET''' is [[evapotranspiration]], '''G<sub>i</sub>''' and '''G<sub>o</sub>''' are the respective [[groundwater]] flux into and out of the basin, '''G<sub>s</sub>''' is the [[groundwater discharge]] into streams, and '''Q<sub>w</sub>''' is groundwater discharge from the basin through [[well]]s. Drainage density relates to the storage and runoff terms. Drainage density relates to the efficiency by which water is carried over the landscape. Water is carried through channels much faster than over hillslopes, as [[overland flow|saturated overland flow]] is slower due to being thinned out and obstructed by [[vegetation]] or pores in the ground.<ref name="GeomorphologyTextbook">{{cite book |last1=Anderson |first1=Robert S. |last2=Anderson |first2=Suzanne P. |title=Geomorphology: The Mechanics and Chemistry of Landscapes |date=2010 |publisher=Cambridge University Press |location=Cambridge, UK |isbn=978-0-521-51978-6}}</ref> Consequently, a drainage basin with a relatively higher drainage density will be more efficiently drained than a higher density one.<ref name="GeomorphologyTextbook" /> Because of the more extensive drainage system in a higher density basin, precipitation entering the basement will, on average, travel a shorter distance over the slower hillslopes before reaching the faster-flowing channels and exit the basin through the channels in less time. Conversely, precipitation entering a lower drainage density basin will take longer to exit the basin due to travelling over the slower hillslope longer.<ref name="GeomorphologyTextbook" /> |
||
In his 1963 paper on drainage density and streamflow, Charles Carlston found that [[baseflow]] into streams is inversely related to the drainage density of the drainage basin: |
In his 1963 paper on drainage density and streamflow, Charles Carlston found that [[baseflow]] into streams is inversely related to the drainage density of the drainage basin: |
||
<center><math>Q_{baseflow}\propto{D_d}^{-2}</math></ |
<div class="center"><math>Q_{baseflow}\propto{D_d}^{-2}</math></div><ref name="Carlston">{{cite book |last1=Carlston |first1=Charles |title=Drainage Density and Streamflow |date=1963 |publisher=US Government Printing Office}}</ref> |
||
This equation represents the effect of drainage density on infiltration. As drainage density increases, baseflow discharge into a stream decreases for a given basin because there is less infiltration to contribute to baseflow.<ref name = "Carlston"/> More of the water entering the drainage basin during a precipitation immediately following a rainfall event exits quickly through streams and does not become infiltration to contribute to baseflow discharge. |
This equation represents the effect of drainage density on infiltration. As drainage density increases, baseflow discharge into a stream decreases for a given basin because there is less infiltration to contribute to baseflow.<ref name = "Carlston"/> More of the water entering the drainage basin during a precipitation immediately following a rainfall event exits quickly through streams and does not become infiltration to contribute to baseflow discharge. |
||
Gregory and Walling (1968) found that the average discharge through a drainage basin is proportional to the square of drainage density: |
Gregory and Walling (1968) found that the average discharge through a drainage basin is proportional to the square of drainage density: |
||
<center><math>Q_{runoff}\propto D_d^2</math></ |
<div class="center"><math>Q_{runoff}\propto D_d^2</math></div><ref name="GregoryAndWalling">{{cite journal |last1=Gregory |first1=K.J. |last2=Walling |first2=D.E. |title=The Variation of Drainage Density Within a Catchment |journal=Hydrological Sciences Journal |date=1968 |volume=13 |issue=2 |pages=61–68 |doi=10.1080/02626666809493583|doi-access=free }}</ref> |
||
This relation illustrates that a higher drainage density environment transports water more efficiently through the basin.<ref name="GeomorphologyTextbook" /> In a relatively low drainage density environment, the lower average discharge results predicted by this relation would be the result of the [[surface runoff]] spending more time travelling over hillslope and having a larger time for infiltration to occur. The increased infiltration results in a decreased surface |
This relation illustrates that a higher drainage density environment transports water more efficiently through the basin.<ref name="GeomorphologyTextbook" /> In a relatively low drainage density environment, the lower average discharge results predicted by this relation would be the result of the [[surface runoff]] spending more time travelling over hillslope and having a larger time for infiltration to occur. The increased infiltration results in a decreased surface runoff according to the water balance equation.<ref name="Fitts" /> |
||
These two equations agree with each other and follow the water balance equation. According to the equations, a basin with high drainage density, the contribution of surface runoff to stream discharge will be high, while that from baseflow will be low. Conversely, a stream in a low drainage density system will have a larger contribution from baseflow and a smaller contribution from [[overland flow]].<ref name="Carlston" /><ref name="GregoryAndWalling" /> |
These two equations agree with each other and follow the water balance equation. According to the equations, a basin with high drainage density, the contribution of surface runoff to stream discharge will be high, while that from baseflow will be low. Conversely, a stream in a low drainage density system will have a larger contribution from baseflow and a smaller contribution from [[overland flow]].<ref name="Carlston" /><ref name="GregoryAndWalling" /> |
||
Line 53: | Line 54: | ||
=== Montgomery and Dietrich (1989) === |
=== Montgomery and Dietrich (1989) === |
||
Montgomery and Dietrich (1989)<ref name="MontgomeryDietrich" /> determined the following equation for drainage density by observing drainage basins in the Tennessee Valley, California: |
Montgomery and Dietrich (1989)<ref name="MontgomeryDietrich" /> determined the following equation for drainage density by observing drainage basins in the Tennessee Valley, California: |
||
<center><math>D=\frac{w_s\rho_wR_0}{W^\ast\rho_sK_z\ sin\ \theta cos\ \theta}\left[1-\left(\frac{tan\theta}{tan\psi}\right)\right]^{-1}</math></ |
<div class="center"><math>D=\frac{w_s\rho_wR_0}{W^\ast\rho_sK_z\ sin\ \theta cos\ \theta}\left[1-\left(\frac{tan\theta}{tan\psi}\right)\right]^{-1}</math></div><ref name="MontgomeryDietrich">{{cite journal |last1=Montgomery |first1=David R. |last2=Dietrich |first2=William E. |title=Source Areas, Drainage Density, and Channel Initiation |journal=Water Resources Research |date=August 1989 |volume=25 |issue=8 |pages=1907–1918|doi=10.1029/WR025i008p01907 |bibcode=1989WRR....25.1907M |citeseerx=10.1.1.658.9871 }}</ref> |
||
Where '''w<sub>s</sub>''' is the mean source width, '''ρ<sub>w</sub>''' is the density of water, '''R<sub>0</sub>''' is the average precipitation rate, '''W*''' is the width of the channel head, '''ρ<sub>s</sub>''' is the saturated bulk density of the soil, '''K<sub>z</sub>''' is the vertical saturated hydraulic conductivity, '''θ''' is the slope at the channel head, and '''φ''' is the soil angle of internal friction. |
Where '''w<sub>s</sub>''' is the mean source width, '''ρ<sub>w</sub>''' is the density of water, '''R<sub>0</sub>''' is the average precipitation rate, '''W*''' is the width of the channel head, '''ρ<sub>s</sub>''' is the saturated bulk density of the soil, '''K<sub>z</sub>''' is the vertical saturated hydraulic conductivity, '''θ''' is the slope at the channel head, and '''φ''' is the soil angle of internal friction. |
||
Line 60: | Line 61: | ||
== Relation to the mean annual flood == |
== Relation to the mean annual flood == |
||
Charles Carlston (1963)<ref name="Carlston" /> determined an equation to express the mean annual flood runoff, Q2.33, for a given drainage basin as a function of drainage density. Carlston found a correlation between the two quantities when plotting data from 15 drainage basins and determined the following equation: |
Charles Carlston (1963)<ref name="Carlston" /> determined an equation to express the mean annual flood runoff, Q2.33, for a given drainage basin as a function of drainage density. Carlston found a correlation between the two quantities when plotting data from 15 drainage basins and determined the following equation: |
||
<center><math>Q_{2.33}=1.3D^2</math></ |
<div class="center"><math>Q_{2.33}=1.3D^2</math></div><ref name="Carlston" /> |
||
Where '''Q''' is in units of cubic feet per second per square mile and '''D<sub>d</sub>''' is in units of inverse miles. From that equation, it is concluded that a drainage basin will adjust itself through erosion such that this equation is satisfied. |
Where '''Q''' is in units of cubic feet per second per square mile and '''D<sub>d</sub>''' is in units of inverse miles. From that equation, it is concluded that a drainage basin will adjust itself through erosion such that this equation is satisfied. |
||
== Effect of vegetation on drainage density == |
== Effect of vegetation on drainage density == |
||
The presence of [[vegetation]] in a drainage basin has multiple effects on the drainage density. Vegetation prevents landslides<ref name="EffectOfVegetation">{{cite journal |last1=Istanbulluoglu |first1=Erkan |last2=Bras |first2=Rafael L. |title=Vegetation-modulated landscape evolution: Effects of vegetation on landscape processes, drainage density, and topography |journal=Journal of Geophysical Research |date=2005 |volume=110 |doi=10.1029/2004JF000249|doi-access=free }}</ref> in the source area of a basin that would result in channel formation as well as decrease the range of drainage density values regardless of soil composition.<ref name="EffectOfVegetation" /> |
The presence of [[vegetation]] in a drainage basin has multiple effects on the drainage density. Vegetation prevents landslides<ref name="EffectOfVegetation">{{cite journal |last1=Istanbulluoglu |first1=Erkan |last2=Bras |first2=Rafael L. |title=Vegetation-modulated landscape evolution: Effects of vegetation on landscape processes, drainage density, and topography |journal=Journal of Geophysical Research |date=2005 |volume=110 |issue=F2 |doi=10.1029/2004JF000249|bibcode=2005JGRF..110.2012I |doi-access=free }}</ref> in the source area of a basin that would result in channel formation as well as decrease the range of drainage density values regardless of soil composition.<ref name="EffectOfVegetation" /> |
||
Vegetation stabilizes the unstable source area in basin and prevents [[channel initiation]].<ref name="EffectOfVegetation" /> Plants stabilize the hillslope that they grow in, which results in physical [[erosion]] processes such as rainsplash, dry |
Vegetation stabilizes the unstable source area in basin and prevents [[channel initiation]].<ref name="EffectOfVegetation" /> Plants stabilize the hillslope that they grow in, which results in physical [[erosion]] processes such as rainsplash, dry ravel, or freezing and thawing processes.<ref name="EffectOfVegetation" /> While there is significant variation between species, plant roots grow in underground networks that holds the soil in place. Because the soil is held in place, it is less prone to erosion from those physical methods.<ref name="EffectOfVegetation" /> Hillslope diffusion was found to decrease exponentially with vegetation cover.<ref name="EffectOfVegetation" /> By stabilizing the hillslope in the source area of the basins, channel initiation, channel initiation is less likely. The erosional processes that may lead to channel initiation are prevented. The increased soil strength also protects against surface runoff erosion, which hinders channel evolution once it has begun.<ref name="EffectOfVegetation" /> |
||
At the basin scale, there are fewer channels in the basin and the drainage density is lower than an unvegetated system. The effect of the vegetation on decreasing the drainage density is not unbounded though. At high vegetative coverage, the effect of increasing the coverage diminishes. This effect imposes an upper limit to the total reduction in drainage density that vegetation can result in.<ref name="EffectOfVegetation" /> |
At the basin scale, there are fewer channels in the basin and the drainage density is lower than an unvegetated system. The effect of the vegetation on decreasing the drainage density is not unbounded though. At high vegetative coverage, the effect of increasing the coverage diminishes. This effect imposes an upper limit to the total reduction in drainage density that vegetation can result in.<ref name="EffectOfVegetation" /> |
||
Line 83: | Line 84: | ||
== Effect of climate change on drainage density == |
== Effect of climate change on drainage density == |
||
Drainage density may also be influenced by climate change. Langbein and Schumm (1958)9 propose an equation for the rate of sediment discharge through catchment as a function of precipitation rate: |
Drainage density may also be influenced by climate change. Langbein and Schumm (1958)9 propose an equation for the rate of sediment discharge through catchment as a function of precipitation rate: |
||
<center><math>P=\frac{aR^\alpha}{1+bR^\gamma}</math></ |
<div class="center"><math>P=\frac{aR^\alpha}{1+bR^\gamma}</math></div><ref name="EffectOfClimateChange">{{cite journal|author1-link=Walter B. Langbein |last1=Langbein |first1=Walter B. |last2=Schumm |first2=Stanley A. |title=Yield of sediment in relation to mean annual precipitation |journal=Eos, Transactions American Geophysical Union |date=1958 |volume=39 |issue=6 |pages=1076–1084|doi=10.1029/TR039i006p01076 |bibcode=1958TrAGU..39.1076L }}</ref> |
||
Where '''P''' is sediment yield, '''R''' is the average effective rainfall, '''α''' ~ 2.3, '''γ''' ~ 3.33, and '''a''' and '''b''' vary depending on units. The graph of this equation has a maximum between 10 and 14 inches and sharp declines on either side of the peak. At lower effective rainfalls, sediment discharge is lower because there is less rainfall to erode the hillslope. At effective rainfalls of greater than 10-14 inches, the decrease in sediment yield is interpreted to be the result of increasing vegetation cover.<ref name="EffectOfClimateChange" /> Increasing precipitation supports denser vegetation coverage and prevents overland flow and other methods of physical erosion. This finding is consistent with the Istanbulluoglu and Bras’ findings on the effect of vegetation on erosion and channel formation.<ref name="EffectOfVegetation" /> |
Where '''P''' is sediment yield, '''R''' is the average effective rainfall, '''α''' ~ 2.3, '''γ''' ~ 3.33, and '''a''' and '''b''' vary depending on units. The graph of this equation has a maximum between 10 and 14 inches and sharp declines on either side of the peak. At lower effective rainfalls, sediment discharge is lower because there is less rainfall to erode the hillslope. At effective rainfalls of greater than 10-14 inches, the decrease in sediment yield is interpreted to be the result of increasing vegetation cover.<ref name="EffectOfClimateChange" /> Increasing precipitation supports denser vegetation coverage and prevents overland flow and other methods of physical erosion. This finding is consistent with the Istanbulluoglu and Bras’ findings on the effect of vegetation on erosion and channel formation.<ref name="EffectOfVegetation" /> |
Revision as of 06:12, 12 January 2024
Drainage density is a quantity used to describe physical parameters of a drainage basin. First described by Robert E. Horton, drainage density is defined as the total length of channel in a drainage basin divided by the total area, represented by the following equation:
The quantity represents the average length of channel per unit area of catchment and has units , which is often reduced to .
Drainage density depends upon both climate and physical characteristics of the drainage basin. Soil permeability (infiltration difficulty) and underlying rock type affect the runoff in a watershed; impermeable ground or exposed bedrock will lead to an increase in surface water runoff and therefore to more frequent streams. Rugged regions or those with high relief will also have a higher drainage density than other drainage basins if the other characteristics of the basin are the same.
When determining the total length of streams in a basin, both perennial and ephemeral streams should be considered.[2] If a drainage basin contained only ephemeral streams, the drainage density by the equation above would be calculated to be zero if only the total length of streams was calculated using only perennial streams. Ignoring ephemeral streams in the calculations does not consider the behavior of the basin during flood events and is therefore not completely representative of the drainage characteristics of the basin.
Drainage density is indicative of infiltration and permeability of a drainage basin, as well as relating to the shape of the hydrograph. Drainage density depends upon both climate and physical characteristics of the drainage basin.
High drainage densities also mean a high bifurcation ratio.
Inverse of drainage density as a physical quantity
Drainage density can be used to approximate the average length of overland flow in a catchment. Horton (1945) used the following equation to describe the average length of overland flow as a function of drainage density:[2]
Where is the length of overland flow with units of length and is the drainage density of the catchment, expressed in units of inverse length.
Considering the geometry of channels on the hillslope, Horton also proposed the following equation:
Where is the channel slope and is the average slope of the ground in the area.
Elementary components of drainage basins
A drainage basin can be defined by three elementary quantities: channels, the hillslope area associated with those channels, and the source areas.[3] The channels are the well-defined segments that efficiently carry water through the catchment. Labeling these features as “channels” rather than “streams” indicates that there need not be a continuous flow of water to capture the behavior of this region as a conduit of water. According to Arthur Strahler’s stream ordering system,[4] the channels are not defined to be any single order or range of orders. Channels of lower orders combine to form higher order channels. The associated hillslope areas are the hillslopes that slope directly into the channels.[3] Precipitation that enters the system on the hillslopes areas and is not lost to infiltration or evapotranspiration enters the channels. The source areas are concave regions of hillslope that are associated with a single channel.[4] Precipitation entering a source area that is not lost to infiltration or evapotranspiration flows through the source area and enters the channel at the channel’s head. Source areas and the hillslope areas associated with channels are differentiated by source areas draining through the channel head, while the associated hillslope areas drain into the rest of the stream.[3] According to Strahler’s stream ordering system, all source areas drain into a primary channel, by the definition of a primary channel.[4]
Bras et al. (1991)[5] describe the conditions that are necessary for channel formation. Channel formation is a concept intimately tied to the formation and evolution of a drainage system and influence the drainage density of catchment. The relation they propose determines the behavior of a given hillslope in response to a small perturbation. They propose the following equation as a relation between source area, source slope, and the sediment flux through this source area:
Where F is the sediment flux, S is the slope of the source area, and a is the source area. The right-hand side of this relation determines channel stability or instability. If the right-hand side of the equation is greater than zero, the hillslope is stable, and small perturbations such as small erosive events do no develop into channels. Conversely, if the right-hand side of the equation is less than zero, Bras et al.[5] determine the hillslope to be unstable, and small erosive structures, such as rills, will tend to grow and form a channel and increase the drainage density of a basin. In this sense, "unstable" is not used in the sense of the gradient of the hillslope being greater than the angle of repose and therefore susceptible to mass wasting, but rather fluvial erosive processes such as sheet flow or channel flow tend to incise and erode to form a singular channel.[5] Therefore, the characteristics of the source area, or potential source area, influence the drainage density and evolution of a drainage basin.[5]
Relation to water balance
Drainage density is tied to the water balance equation:
Where is the change in reservoir storage, R is precipitation, ET is evapotranspiration, Gi and Go are the respective groundwater flux into and out of the basin, Gs is the groundwater discharge into streams, and Qw is groundwater discharge from the basin through wells. Drainage density relates to the storage and runoff terms. Drainage density relates to the efficiency by which water is carried over the landscape. Water is carried through channels much faster than over hillslopes, as saturated overland flow is slower due to being thinned out and obstructed by vegetation or pores in the ground.[7] Consequently, a drainage basin with a relatively higher drainage density will be more efficiently drained than a higher density one.[7] Because of the more extensive drainage system in a higher density basin, precipitation entering the basement will, on average, travel a shorter distance over the slower hillslopes before reaching the faster-flowing channels and exit the basin through the channels in less time. Conversely, precipitation entering a lower drainage density basin will take longer to exit the basin due to travelling over the slower hillslope longer.[7]
In his 1963 paper on drainage density and streamflow, Charles Carlston found that baseflow into streams is inversely related to the drainage density of the drainage basin:
This equation represents the effect of drainage density on infiltration. As drainage density increases, baseflow discharge into a stream decreases for a given basin because there is less infiltration to contribute to baseflow.[8] More of the water entering the drainage basin during a precipitation immediately following a rainfall event exits quickly through streams and does not become infiltration to contribute to baseflow discharge. Gregory and Walling (1968) found that the average discharge through a drainage basin is proportional to the square of drainage density:
This relation illustrates that a higher drainage density environment transports water more efficiently through the basin.[7] In a relatively low drainage density environment, the lower average discharge results predicted by this relation would be the result of the surface runoff spending more time travelling over hillslope and having a larger time for infiltration to occur. The increased infiltration results in a decreased surface runoff according to the water balance equation.[6]
These two equations agree with each other and follow the water balance equation. According to the equations, a basin with high drainage density, the contribution of surface runoff to stream discharge will be high, while that from baseflow will be low. Conversely, a stream in a low drainage density system will have a larger contribution from baseflow and a smaller contribution from overland flow.[8][9]
Relation to hydrographs
The discharge through the central stream draining a catchment reflects the drainage density, which makes it a useful diagnostic for predicting the flooding behavior of a catchment following a storm event due to being intimately tied to the hydrograph.[7] The material that overland flow travels over is one factor that influences the speed that water can flow out of a catchment. Water flows significantly slower over hillslopes compared to channels that form to efficiently carry water and other flowing material.[7] According to Horton’s interpretation of half of the inverse of drainage density as the average length of overland flow[2] implies that overland flow in high-drainage environments will reach a fast-flowing channel faster over a shorter range. On the hydrograph, the peak is higher and occurs over a shorter range. This more compact and higher peak is often referred to as being “flashy”.[7]
The timing of the hydrograph in relation to the peak of the hyetograph is influenced by the drainage density.[7] The water that enters a high-drainage watershed during a storm will reach a channel relatively fast and travel in the high-velocity channels to the outlet of the watershed in a relatively short time. Conversely, the water entering a low drainage density basin will, on average, have to travel a longer distance over the low velocity hillslope to reach the channels. As a result, the water will require more time to reach the exit of the catchment. The lag time between the peak of the hyetograph and the hydrograph is then inversely related to drainage density; as drainage density increases, water is more efficiently drained from the basin and the lag time decreases.[7]
Another impact on the hydrograph that drainage density has is a steeper falling limb following the storm event due to its impact on both overland flow and baseflow.[7][10] The falling limb occurs after the peak of the hydrograph curve and is when overland flow is decreasing back to ambient levels. In higher drainage systems, the overland flow reaches the channels quicker resulting in a narrower spread in the falling limb. Baseflow is the other contributor to the hydrograph. The peak of baseflow to the channels will occur after the quick-flow peak because groundwater flow is much slower than quick-flow. Because the baseflow peak occurs after the quick-flow peak, the baseflow peak influences the shape of the falling limb.[10] According to the proportionality put forth by Gregory and Walling,[9] as drainage density increases, the contribution of baseflow to the falling limb of the hydrograph diminishes. During a storm event in a high drainage density basin, there is little water that infiltrates into the ground as infiltration because water spends less time flowing over the surface in the catchment before exiting through the central channel. Because there is little water that enters the water as infiltration, baseflow will contribute only a small part to falling limb. The falling limb is thus quite steep. Conversely, a low drainage system will have a shallower falling limb. According to Gregory and Walling’s relation,[9] the decrease in drainage density results in an increase in baseflow to the channels and a more gradual decrease in the hydrograph.
Formula for drainage density
Montgomery and Dietrich (1989)
Montgomery and Dietrich (1989)[3] determined the following equation for drainage density by observing drainage basins in the Tennessee Valley, California:
Where ws is the mean source width, ρw is the density of water, R0 is the average precipitation rate, W* is the width of the channel head, ρs is the saturated bulk density of the soil, Kz is the vertical saturated hydraulic conductivity, θ is the slope at the channel head, and φ is the soil angle of internal friction.
R0, the average precipitation term, shows the dependence of drainage density on climate. With all other factors being constant, an increase in precipitation in the drainage basin results in an increase in drainage density.[3] A decrease in precipitation, such as in an arid environment, results in a lower drainage density. The equation also shows the dependence on the physical characteristics and lithology of the drainage basin. Materials with a low hydraulic conductivities, such as clay or solid rock,[6] would result in a higher-drainage density system. Because of the low hydraulic conductivity, there is little water lost to infiltration and that water exits the system as runoff and can contribute to erosion. In a basin with a higher vertical hydraulic conductivity, water more effectively infiltrates into the ground and does not contribute to saturated overland flow erosion, resulting in a less developed channel system and therefore lower drainage density.[3]
Relation to the mean annual flood
Charles Carlston (1963)[8] determined an equation to express the mean annual flood runoff, Q2.33, for a given drainage basin as a function of drainage density. Carlston found a correlation between the two quantities when plotting data from 15 drainage basins and determined the following equation:
Where Q is in units of cubic feet per second per square mile and Dd is in units of inverse miles. From that equation, it is concluded that a drainage basin will adjust itself through erosion such that this equation is satisfied.
Effect of vegetation on drainage density
The presence of vegetation in a drainage basin has multiple effects on the drainage density. Vegetation prevents landslides[11] in the source area of a basin that would result in channel formation as well as decrease the range of drainage density values regardless of soil composition.[11]
Vegetation stabilizes the unstable source area in basin and prevents channel initiation.[11] Plants stabilize the hillslope that they grow in, which results in physical erosion processes such as rainsplash, dry ravel, or freezing and thawing processes.[11] While there is significant variation between species, plant roots grow in underground networks that holds the soil in place. Because the soil is held in place, it is less prone to erosion from those physical methods.[11] Hillslope diffusion was found to decrease exponentially with vegetation cover.[11] By stabilizing the hillslope in the source area of the basins, channel initiation, channel initiation is less likely. The erosional processes that may lead to channel initiation are prevented. The increased soil strength also protects against surface runoff erosion, which hinders channel evolution once it has begun.[11]
At the basin scale, there are fewer channels in the basin and the drainage density is lower than an unvegetated system. The effect of the vegetation on decreasing the drainage density is not unbounded though. At high vegetative coverage, the effect of increasing the coverage diminishes. This effect imposes an upper limit to the total reduction in drainage density that vegetation can result in.[11]
Vegetation also narrows the range of drainage density values for basins of various soil composition.[11] Unvegetated basins can have a large range in drainage densities, from low to high. Drainage density is related to the ease at which channels can form. According to Montgomery and Dietrich’s equation, drainage density is a function of vertical hydraulic conductivity. Coarse-grained sediment like sand would have a higher hydraulic conductivity and are predicted by the equation to form a relatively higher drainage density system than a system formed by finer silt with a lower hydraulic conductivity.[6]
Forest fires play an indirect role in a basin’s drainage density. Forest fires, both natural and unnatural, destroy some or all of the existing vegetation, which removes the stability that the plants and their roots provide. Newly destabilized hillslope in the basin is then susceptible to channel formation processes, and drainage density of the basin may increase until the vegetation grows back to the previous state. The type of plants and the associated depth and density of the plant roots determine how strongly the soil is held in place as well as the intensity of the forest fire in killing and removing the vegetation. Computer simulation experiments have validated that drainage density will be higher in regions that have more frequent forest fires.[11]
Relation to flood hydrograph
The discharge through the central stream draining a catchment reflects the drainage density, which makes it a useful diagnostic for predicting the flooding behavior of a catchment following a storm event due to being intimately tied to the hydrograph.[7] The material that overland flow travels over is one factor that influences the speed that water can flow out of a catchment. Water flows significantly slower over hillslopes compared to channels that form to efficiently carry water and other flowing material. According to Horton’s interpretation of half of the inverse of drainage density as the average length of overland flow[2] implies that overland flow in high-drainage environments will reach a fast-flowing channel faster over a shorter range. On the hydrograph, the peak is higher and occurs over a shorter range. This more compact and higher peak is often referred to as being “flashy”.[7]
The timing of the hydrograph in relation to the peak of the hyetograph is influenced by the drainage density.[7] The water that enters a high-drainage watershed during a storm will reach a channel relatively fast and travel in the high-velocity channels to the outlet of the watershed in a relatively short time. Conversely, the water entering a low drainage density basin will, on average, have to travel a longer distance over the low velocity hillslope to reach the channels. As a result, the water will require more time to reach the exit of the catchment. The lag time between the peak of the hyetograph and the hydrograph is then inversely related to drainage density; as drainage density increases, water is more efficiently drained from the basin and the lag time decreases.[7]
Another impact on the hydrograph that drainage density has is a steeper falling limb following the storm event due to its impact on both overland flow and baseflow.[7][10] The falling limb occurs after the peak of the hydrograph curve and is when overland flow is decreasing back to ambient levels. In higher drainage systems, the overland flow reaches the channels quicker resulting in a narrower spread in the falling limb. Baseflow is the other contributor to the hydrograph. The peak of baseflow to the channels will occur after the quick-flow peak because groundwater flow is much slower than quick-flow.[10] Because the baseflow peak occurs after the quick-flow peak, the baseflow peak influences the shape of the falling limb.4 According to the proportionality put forth by Gregory and Walling,[9] as drainage density increases, the contribution of baseflow to the falling limb of the hydrograph diminishes. During a storm event in a high drainage density basin, there is little water that infiltrates into the ground as infiltration because water spends less time flowing over the surface in the catchment before exiting through the central channel. Because there is little water that enters the water as infiltration, baseflow will contribute only a small part to falling limb. The falling limb is thus quite steep. Conversely, a low drainage system will have a shallower falling limb. According to Gregory and Walling’s relation,[9] the decrease in drainage density results in an increase in baseflow to the channels and a more gradual decrease in the hydrograph.
Effect of climate change on drainage density
Drainage density may also be influenced by climate change. Langbein and Schumm (1958)9 propose an equation for the rate of sediment discharge through catchment as a function of precipitation rate:
Where P is sediment yield, R is the average effective rainfall, α ~ 2.3, γ ~ 3.33, and a and b vary depending on units. The graph of this equation has a maximum between 10 and 14 inches and sharp declines on either side of the peak. At lower effective rainfalls, sediment discharge is lower because there is less rainfall to erode the hillslope. At effective rainfalls of greater than 10-14 inches, the decrease in sediment yield is interpreted to be the result of increasing vegetation cover.[12] Increasing precipitation supports denser vegetation coverage and prevents overland flow and other methods of physical erosion. This finding is consistent with the Istanbulluoglu and Bras’ findings on the effect of vegetation on erosion and channel formation.[11]
The Caineville Badlands
The badlands of Caineville, Utah are often cited as a region of extremely high drainage density. The region features steep slopes, high relief, an arid climate, and a complete absence of vegetation.[13][7] Because the slopes of hillslopes are often greater than the angle of repose, the dominant erosional process in the Caineville badlands is mass wasting.[13] There is no vegetation to provide stability to the slopes and increase the angle of repose and prevent mass-wasting. The regions below the angle of repose, however, are still generally at a significant angle, and hillslope diffusion, according to the following relation, is still a significant source of erosion:
Where Ks is a coefficient of diffusivity of the hillslope, z is the elevation of the hillslope, and x is horizontal distance.
The range of drainage densities in the Caineville Badlands illustrates the complicated nature of drainage densities in low-precipitation environments.[11] In a study on the region, Alan Howard (1996) found that the effect of increasing relief angles in different basins did not have a constant effect on the drainage density.[13] For regions of relatively low relief, drainage density and relief are positively correlated. This occurs until a threshold is reached at a higher relief ratio, when increase the slope ratio is accompanied by a decrease in drainage density.[13] This is interpreted by Howard to be a result of the critical source area needed to support a channel increasing. At a higher slope, the erosion is faster and more efficiently funneled through fewer channels.[13] The smaller number of channels results in a smaller drainage density for the basin.
This qualitative topographic map of a section of a section of the Caineville Badlands shows the extensive drainage network in the arid environment. Relating to Montgomery and Dietrich’s definition of the elementary parts of a drainage basin,[3] the source area for each of the channels is relatively very small, resulting in a large number of channels forming. The image of the Caineville Badlands displays the lack of vegetation and numerous channels. The Caineville Badlands are located in an arid environment, receiving an average of 125mm of precipitation per year.[13] This low precipitation contrasts with Montgomery and Dietrich’s equation of drainage density, which predicts that drainage density should be low where rainfall is low.[3] This behavior is more consistent with Langbein and Schumm’s expression of erosion rate as a function of rainfall.[12] According to the equation, erosion will increase with precipitation up to a point where the precipitation can support stabilizing vegetation. The lack of vegetation present in the image of the Caineville Badlands implies that the rainfall rate of this region is below the critical rainfall amount vegetation can be supported.
References
- ^ Horton, Robert E. (June 1932). "Drainage-basin characteristics". Transactions of the American Geophysical Union. 13 (1): 350–361. Bibcode:1932TrAGU..13..350H. doi:10.1029/TR013i001p00350.
- ^ a b c d e Horton, Robert E. (1945-03-01). "Erosional Development of Streams and Their Drainage Basins; Hydrophysical Approach to Quantitative Morphology". GSA Bulletin. 56 (3): 275–370. Bibcode:1945GSAB...56..275H. doi:10.1130/0016-7606(1945)56[275:EDOSAT]2.0.CO;2. S2CID 129509551.
- ^ a b c d e f g h i Montgomery, David R.; Dietrich, William E. (August 1989). "Source Areas, Drainage Density, and Channel Initiation". Water Resources Research. 25 (8): 1907–1918. Bibcode:1989WRR....25.1907M. CiteSeerX 10.1.1.658.9871. doi:10.1029/WR025i008p01907.
- ^ a b c Strahler, Arthur N. (December 1957). "Quantitative Analysis of Watershed Geomorphology". Transactions, American Geophysical Union. 38 (6): 913–920. Bibcode:1957TrAGU..38..913S. doi:10.1029/TR038i006p00913.
- ^ a b c d e Tarboton, David G.; Bras, Rafael L.; Rodriguez-Iturbe, Ignacio (October 19, 1991). "A physical basis for drainage density". Geomorphology. 5 (1–2): 59–76. doi:10.1016/0169-555X(92)90058-V.
- ^ a b c d Fitts, Charles R. (2013). Groundwater Science (2 ed.). Waltham, MA: Elsevier Inc. p. 14. ISBN 978-0-12-384705-8.
- ^ a b c d e f g h i j k l m n o p Anderson, Robert S.; Anderson, Suzanne P. (2010). Geomorphology: The Mechanics and Chemistry of Landscapes. Cambridge, UK: Cambridge University Press. ISBN 978-0-521-51978-6.
- ^ a b c d e Carlston, Charles (1963). Drainage Density and Streamflow. US Government Printing Office.
- ^ a b c d e f Gregory, K.J.; Walling, D.E. (1968). "The Variation of Drainage Density Within a Catchment". Hydrological Sciences Journal. 13 (2): 61–68. doi:10.1080/02626666809493583.
- ^ a b c d Erickson, Timothy O.; Stefan, Heinz G. (June 2008). "Baseflow Analysis of the Upper Vermillion River, Dakota County, Minnesota". Minnesota Pollution Control Agency. Project Report No. 507.
- ^ a b c d e f g h i j k l Istanbulluoglu, Erkan; Bras, Rafael L. (2005). "Vegetation-modulated landscape evolution: Effects of vegetation on landscape processes, drainage density, and topography". Journal of Geophysical Research. 110 (F2). Bibcode:2005JGRF..110.2012I. doi:10.1029/2004JF000249.
- ^ a b c Langbein, Walter B.; Schumm, Stanley A. (1958). "Yield of sediment in relation to mean annual precipitation". Eos, Transactions American Geophysical Union. 39 (6): 1076–1084. Bibcode:1958TrAGU..39.1076L. doi:10.1029/TR039i006p01076.
- ^ a b c d e f Howard, Alan D. (1996). "Badland Morphology and Evolution: Interpretation Using a Simulation Model". Earth Surface Processes and Landforms. 22 (3): 211–227. doi:10.1002/(SICI)1096-9837(199703)22:3<211::AID-ESP749>3.0.CO;2-E.