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 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 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 in a basin, both perennial stream and should be considered. If a drainage basin contained only ephemeral streams, then the drainage density by the equation above would be calculated to be zero if 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, and relates 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.
where is the length of overland flow and is the drainage density of the catchment.
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.
Bras et al. (1991) 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 is the sediment flux, is the slope of the source area, and 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, then 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, then Bras et al. 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. Therefore, the characteristics of the source area, or potential source area, influence the drainage density and evolution of a drainage basin.
where is the change in reservoir storage, is precipitation, is evapotranspiration, and are the respective groundwater flux into and out of the basin, is the groundwater discharge into streams, and is groundwater discharge from the basin through . 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 is slower due to being thinned out and obstructed by vegetation or pores in the ground.
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. More of the water entering the drainage basin during and 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. 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.
These two equations agree with each other and follow the water balance equation. According to the equations, in 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.
The timing of the hydrograph in relation to the peak of the hyetograph is influenced by the drainage density. 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.
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. 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. According to the proportionality put forth by Gregory and Walling, 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, the decrease in drainage density results in an increase in baseflow to the channels and a more gradual decrease in the hydrograph.
where is the mean source width, is the density of water, is the average precipitation rate, is the width of the channel head, is the saturated bulk density of the soil, is the vertical saturated hydraulic conductivity, is the slope at the channel head, and is the soil angle of internal friction.
, 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. 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, 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.
where is in units of cubic feet per second per square mile and 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.
Vegetation stabilizes the unstable source area in basin and prevents channel initiation. Plants stabilize the hillslope that they grow in, which results in physical erosion processes such as rain splash, dry ravel, or freezing and thawing processes. 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. Hillslope diffusion was found to decrease exponentially with vegetation cover. By stabilizing the hillslope in the source area of the basins, 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.
At the basin scale, there are fewer channels in the basin, and the drainage density is lower than an unvegetated system. However, the effect of the vegetation on decreasing the drainage density is not unbounded. 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.
Vegetation also narrows the range of drainage density values for basins of various soil composition. 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 sediments 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.
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 hillslopes in the basin are then susceptible to channel formation processes, and the 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 simulations have validated that drainage density will be higher in regions that have more frequent forest fires.
where is sediment yield, is the average effective rainfall, , , and and vary depending on units. The graph of this equation has a maximum between 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 greater effective rainfalls, the decrease in sediment yield is interpreted to be the result of increasing vegetation cover. 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.
where is a coefficient of diffusivity of the hillslope, is the elevation of the hillslope, and is horizontal distance.
The range of drainage densities in the Caineville Badlands illustrates the complicated nature of drainage densities in low-precipitation environments. 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. 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 an increase the slope ratio is accompanied by a decrease in drainage density. 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. The smaller number of channels results in a smaller drainage density for the basin.
This qualitative topographic map 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, 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 of precipitation per year. 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. This behavior is more consistent with Langbein and Schumm's expression of erosion rate as a function of rainfall. 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 at which vegetation can be supported.
|
|