Applications of the random model of drainage basin composition

The random model of drainage basin composition is founded on the assumptions that (a) natural channels are topologically random in the absence of geological controls and (b) for channel networks developed in similar environments, the exterior and interior link lengths are independent random variables with a common distribution for each type. The effectiveness of this model in estimating the values of geomorphic variables and in explaining and predicting geomorphic relationships is illustrated by several examples. The data required for these examples were obtained from map studies of 30 channel networks, comprising a total of about 8700 links, in eastern Kentucky. A common factor in the success of all three applications of the model is the way in which the planimetric features of drainage basins are determined by their underlying topologic structure.     

The analysis of drainage network composition

The Horton method of analyzing drainage network composition is reviewed, with the conclusion that it has not been very effective either in improving understanding or in developing useful methods of characterizing drainage basins. New methods which are based on the link rather than the Horton (or Strahler) stream are described. A number of detailed examples of the application of these new methods to the topologic and geometric properties of networks are provided. The results are compared with the predictions of the random model. Data used in the analysis were obtained from 1:24,000 U.S.G.S. topographic maps of eastern Kentucky. Thirty drainage basins were selected and their channel networks were outlined first by the contour-crenulation (CC) method and then by another, more objective, method (SC) in which stream sources were identified by a quantitative slope criterion. The CC and SC samples comprise about 8,700 and 1,700 links, respectively. The three most important results of the analysis are: (1) the channel networks are slightly but significantly more elongated than predicted by the random model, (2) there are fewer second magnitude links than predicted, and (3) the length distribution for interior links depends upon the kind of link (interior or exterior) joined downstream. These features are found in both CC and SC networks.     

DRAINAGE NETWORKS AFTER WILDFIRE

Predicting runoff and erosion from watersheds burned by wildfires requires an understanding of the three-dimensional structure of both hillslope and channel drainage networks. We investigate the small- and large-scale structures of drainage networks using field studies and computer analysis of 30- m digital elevation model. Topologic variables were derived from a composite 30-m DEM, which included 14 order 6 watersheds within the Pikes Peak batholith. Both topologic and hydraulic variables were measured in the field in two smaller burned watersheds （3.7 and 7.0 hectares） located within one of the order 6 watersheds burned by the 1996 Buffalo Creek Fire in Central Colorado. Horton ratios of topologic variables （stream number, drainage area, stream length, and stream slope） for small-scale and large-scale watersheds are shown to scale geometrically with stream order （i.e., to be scale invariant）. However, the ratios derived for the large-scale drainage networks could not be used to predict the rill and gully drainage network structure. Hydraulic variables （width, depth, cross- sectional area, and bed roughness） for small-scale drainage networks were found to be scale invariant across 3 to 4 stream orders. The relation between hydraulic radius and cross-sectional area is similar for fills and gullies, suggesting that their geometry can be treated similarly in hydraulic modeling. Additionally, the rills and gullies have relatively small width-to-depth ratios, implying sidewall friction may be important to the erosion and evolutionary process relative to main stem channels.     

Gravel sediment routing from widespread,low- intensity landscape disturbance,Current River Basin,Missouri

During the last 160 years, land-use changes in the Ozarks have had the potential to cause widespread, low-intensity delivery of excess amounts of gravel-sized sediment to stream channels. Previous studies have indicated that this excess gravel bedload is moving in wave-like forms through Ozarks drainage basins. The longitudinal, areal distribution of gravel bars along 160 km of the Current River, Missouri, was evaluated to determine the relative effects of valley-scale controls, tributary basin characteristics, and lagged sediment transport in creating areas of gravel accumulations. The longitudinal distribution of gravel-bar area shows a broad scale wave-like form with increases in gravel-bar area weakly associated with tributary junctions. Secondary peaks of gravel area with 1·8–4·1 km spacing (disturbance reaches) are superimposed on the broad form. Variations in valley width explain some, but not all, of the short-spacing variation in gravel-bar area. Among variables describing tributary drainage basin morphometry, present-day land use and geologic characteristics, only drainage area and road density relate even weakly to gravel-bar areal inventories. A simple, channel network-based sediment routing model shows that many of the features of the observed longitudinal gravel distribution can be replicated by uniform transport of sediment from widespread disturbances through a channel network. These results indicate that lagged sediment transport may have a dominant effect on the synoptic spatial distribution of gravel in Ozarks streams; present-day land uses are only weakly associated with present-day gravel inventories; and valley-scale characteristics have secondary controls on gravel accumulations in disturbance reaches. Copyright © 1999 John Wiley & Sons, Ltd.     

Basal sediment evacuation by subglacial meltwater: suspended sediment transport from Haut Glacier d'Arolla,Switzerland

Proglacial suspended sediment transport was monitored at Haut Glacier d'Arolla, Switzerland, during the 1998 melt season to investigate the mechanisms of basal sediment evacuation by subglacial meltwater. Sub‐seasonal changes in relationships between suspended sediment transport and discharge demonstrate that the structure and hydraulics of the subglacial drainage system critically influenced how basal sediment was accessed and entrained. Under hydraulically inefficient subglacial drainage at the start of the melt season, sediment availability was generally high but sediment transport increased relatively slowly with discharge. Later in the melt season, sediment transport increased more rapidly with discharge as subglacial meltwater became confined to a spatially limited network of channels following removal of the seasonal snowpack from the ablation area. Flow capacity is inferred to have increased more rapidly with discharge within subglacial channels because rapid changes in discharge during highly peaked diurnal runoff cycles are likely to have been accommodated largely by changes in flow velocity. Basal sediment availability declined during channelization but increased throughout the remainder of the monitored period, resulting in very efficient basal sediment evacuation over the peak of the melt season. Increased basal sediment availability during the summer appears to have been linked to high diurnal water pressure variation within subglacial channels inferred from the strong increase in flow velocity with discharge. Basal sediment availability therefore appears likely to have been increased by (1) enhanced local ice‐bed separation leading to extra‐channel flow excursions and[sol ]or (2) the deformation of basal sediment towards low‐pressure channels due to a strong diurnally reversing hydraulic gradient between channels and areas of hydraulically less‐efficient drainage. Copyright © 2005 John Wiley & Sons, Ltd.     

Understanding Amazonian fluvial rias based on a Late Pleistocene–Holocene analog

Fluvial rias are elongated lakes at tributary mouths that can reach dozens of kilometers in length, constituting one of the most remarkable features in the Amazonian landscape. Thus far, definitive data which documents the genesis of fluvial rias have not been published. The main goal of this work was to integrate morphological, sedimentological and chronological information in order to characterize fluvial paleorias in the interfluve of the Purus and Madeira Rivers and discuss the most likely hypothesis for their genesis. These paleorias were first observed through remote sensing imagery as several elongated and interconnecting belts of open vegetation that are in sharp contact with the surrounding dense forest. The belts are branched and form a dendritic pattern similar to many modern drainage networks. The sedimentary record of these belts revealed the prevalence of sharp‐based sandstones and mudstones arranged into fining‐upward successions, which are compatible with deposition within channels. Active channel and abandoned channel deposits were recognized. These are topped by continuous mudstones related to rapid channel abandonment and formation of a low energy basin or ria environment. Radiocarbon dating of these deposits recorded only Late Pleistocene and Holocene ages ranging from 21 547–22 285 cal yr bp to 5928–6124 cal yr bp . This chronology for sediment deposition is not compatible with the hypothesis of Amazonian rias being formed by fluvial erosion during the Last Glaciation Maximum low sea level, with sediment accumulation during the subsequent Holocene transgression. Instead, the studied paleorias record previous tributaries of the Madeira River that became abandoned as the position of this river shifted southeastward and its interfluve tilted northward, inverting the drainage systems. Therefore, a neotectonic origin of some Amazonian paleorias seems most likely. This hypothesis should be considered in further investigations aiming at understanding the origin of numerous modern fluvial rias that typify the Amazonian landscape. Copyright © 2014 John Wiley & Sons, Ltd.     

The Quaternary development of tributary channels to the Nile River at Kom Ombo area,Eastern Desert of Egypt,and their implication for groundwater resources

The dry wadis ‘ephemeral channel’ constituting the main tributaries to the Nile River in Kom Ombo are structurally and tectonically controlled and exhibit complex drainage pattern. This complicated drainage pattern is inherited from the morpho‐tectonic evolution of the ancestral Nile River (‘Protonile’), which drained the Eastern Desert during the Middle Pleistocene. A digital elevation model derived from Shuttle Radar Topography Mission data is used to delineate the contemporary drainage networks and their catchments. Satellite images acquired during a flash flood event were used to validate the delineated watershed divides and flow pathways, particularly where the courses of dry wadis are interlocked. Currently, the westward flow of Wadi Abu‐Suberah is derived from a small area in the Eastern Desert, as the palaeo‐upper reaches of this wadi were captured due to tectonic movements along NW/SE and N/S faults by wadis in the Kharit and Elewa areas. The influence of these tectonic movements on groundwater distribution is also shown where the deep Nubian aquifer discharges its water into the Quaternary aquifer through fault planes. The northward flowing ‘Protonile’ main course has tectonically shifted from the Gallaba plain in the Western Desert, eastward to the current Nile River course. This shift has produced several cut‐off segments of the palaeo‐tributary drainage that was originally flowing westward towards the old ‘Protonile’ main course in the Gallaba plain. However, these segments have reversed their surface run‐off flow directions eastward towards the current Nile course; they could include potential groundwater resources, as their alluvium may be still recharged by the Nile River. Copyright © 2010 John Wiley & Sons, Ltd.     

Controls on valley spacing in landscapes subject to rapid base‐level fall

What controls the architecture of drainage networks is a fundamental question in geomorphology. Recent work has elucidated the mechanisms of drainage network development in steadily uplifting landscapes, but the controls on drainage‐network morphology in transient landscapes are relatively unknown. In this paper we exploit natural experiments in drainage network development in incised Plio‐Quaternary alluvial fan surfaces in order to understand and quantify drainage network development in highly transient landscapes, i.e. initially unincised low‐relief surfaces that experience a pulse of rapid base‐level drop followed by relative base‐level stasis. Parallel drainage networks formed on incised alluvial‐fan surfaces tend to have a drainage spacing that is approximately proportional to the magnitude of the base‐level drop. Numerical experiments suggest that this observed relationship between the magnitude of base‐level drop and mean drainage spacing is the result of feedbacks among the depth of valley incision, mass wasting and nonlinear increases in the rate of colluvial sediment transport with slope gradient on steep valley side slopes that lead to increasingly wide valleys in cases of larger base‐level drop. We identify a threshold magnitude of base‐level drop above which side slopes lengthen sufficiently to promote increases in contributing area and fluvial incision rates that lead to branching and encourage drainage networks to transition from systems of first‐order valleys to systems of higher‐order, branching valleys. The headward growth of these branching tributaries prevents the development of adjacent, ephemeral drainages and promotes a higher mean valley spacing relative to cases in which tributaries do not form. Model results offer additional insights into the response of initially unincised landscapes to rapid base‐level drop and provide a preliminary basis for understanding how varying amounts of base‐level change influence valley network morphology. Copyright © 2015 John Wiley & Sons, Ltd.     

Watershed modeling of rainfall excess transformation into runoff

In this paper an attempt is made to present a distributed physiographic conceptual model that uses the principles of flow continuity and momentum. For this purpose, the watershed under study is divided into subwatersheds keeping in view the drainage patterns and characteristics. Then the main tributaries are identified and their drainage areas are delineated to form tributary subwatersheds. The main channel subwatersheds have taken care of the remaining area in the vicinity of the main channel. The kinematic wave theory is applied for the overland runoff computations from these subwatersheds. Further, the overland flows are superimposed onto the main channel. The dynamic wave theory is used to route the flows through the main channel to compute the watershed responses at the outlet. The proposed model is tested onto a natural watershed. The computations were performed for few storm events. Comparison of the significant parameters of the computed and the observed hydrographs shows that the maximum relative error in prediction is 5.8%. Thus, the results are satisfactory. Better results can be obtained when measured rainfall-excess data are available or a more realistic loss index is adopted for rainfall-excess separation.     

Hillslope–channel sediment transfer in a slope failure event: Wet Swine Gill,Lake District,northern England

This paper describes and analyses a hillslope–channel slope failure event that occurred at Wet Swine Gill, Lake District, northern England. This comprised a hillslope slide (180 m³, c. 203 ± 36 t), which coupled with the adjacent stream, resulting in a channelized debris flow and fluvial flood. The timing of the event is constrained between January and March 2002. The hillslope failure occurred in response to a rainfall/snowmelt trigger, on ground recently disturbed by a heather moorland fire and modified by artificial drainage. Slide and flow dynamics are estimated using reconstructed velocity and discharge values along the sediment transfer path. There is a rapid downstream reduction in both maximum velocity, from 9·8 to 1·3 m s^?1; and maximum discharge, ranging from 33·5 to 2·4 m³ s^?1. A volumetric sediment budget quantified a high degree of coupling between the hillslope and immediate channel (～92%: 167 m³), but virtually all of the sediment was retained in the first‐order tributary channel. Approximately 44% (81 m³) of the slide volume was retained in the run‐up deposit, and termination of the debris flow prior to the main river meant that the remainder did not discharge into the fluvial system downstream. These results suggest poor transmission of sediment to the main river at the time of the event, but importantly an increase in available material for post‐event sediment transfer processes within the small upland tributary. Copyright © 2007 John Wiley & Sons, Ltd.     

Martian oceans, valleys and climate

The new Mars Global Surveyor altimetry shows that the heavily cratered southern hemisphere of Mars is 5 km higher than the sparely cratered plains of the northern hemisphere. Previous suggestions that oceans formerly occupied the northern plains as evidenced by shorelines are partly supported by the new data. A previously identified outer boundary has a wide range of elevations and is unlikely to be a shoreline but an inner contact with a narrow range of elevations is a more likely candidate. No shorelines are visible in the newly acquired, 1.5 metre/pixel imaging. Newly imaged valleys provide strong support for sustained or episodic flow of water across the Martian surface. A major surprise, however, is the near absence of valleys less than 100 m across. Martian valleys seemingly do not divide into ever smaller valleys as terrestrial valleys commonly do. This could be due to lack of precipitation or lack of surface runoff because of high infiltration rates. High erosion rates and formation of valley networks supports warm climates and presence of large bodies of water during heavy bombardment. The climate history and fate of the water after heavy bombardment remain controversial.     

Rates and patterns of thermal mixing at a small stream confluence under variable incoming flow conditions

Confluences are important locations for river mixing within drainage networks, yet few studies have examined in detail the dynamics of mixing within confluences. This study examines the influence of momentum flux ratio, the scale of the flow (cross‐sectional area) and the density differences between incoming flows on thermal mixing at a small stream confluence. Results reveal that rates and patterns of thermal mixing depend on event‐specific combinations of the three factors. The mixing interface at this confluence is generally distorted towards the mouth of the lateral tributary by strong helical motion associated with curvature of flow from the lateral tributary as it aligns with the downstream channel. As the momentum flux from the lateral tributary increases, mixing is enhanced because helical motion from the curving tributary flow expands over the width of the downstream channel. The cross‐sectional area of the flow is negatively correlated with mixing rates, suggesting that the amount of mixing over a fixed distance downstream of the confluence is inversely related to the scale of the flow. Density differences are not strongly related to rates of mixing. Results confirm that mixing rates within the region of confluent flow interaction can be highly variable among flow events with different incoming flow conditions, but that, in general, length scales of mixing are short, and rates of mixing are high at this small confluence compared with those typically documented at large‐river confluences. Copyright © 2015 John Wiley & Sons, Ltd.     

How does network structure and complexity in river systems affect population abundance and persistence?

Models of habitat connectivity, and how network structure and connectivity affects resident populations, are increasingly being developed for terrestrial habitat networks. Rivers, unlike many terrestrial habitat networks, follow a rigid hierarchical branching structure to form dendritic networks. It has been hypothesised that this unique structure must have implications for population processes. We developed a theoretical model to relate local-scale habitat quality and connectivity to landscape-scale population dynamics of mobile organisms (e.g. fish, aquatic invertebrates). River networks were modelled as directed graphs, with nodes being habitat patches, and edges the connections between them. Using population simulation analyses, we investigated the effects of network structure on population abundance and persistence. Network structural complexity affected landscape-scale population abundance, but the apparent effect depended upon how structure was quantified. There were no noticeable effects of dendritic network structure on population persistence. Previous research on the effects of habitat network structure on population persistence has used metapopulation patch occupancy models, which do not directly consider population dynamics. Our results show that spatially-explicit population modelling is possible, and that it provides information beyond that possible with patch occupancy models (e.g. population abundance). More importantly, it calls into question whether metapopulation models provide an adequate representation of population dynamics in dendritic habitat networks.     

Sediment storage and transport in Pancho Rico Valley during and after the Pleistocene‐Holocene transition,Coast Ranges of central California (Monterey County)

Factors influencing sediment transport and storage within the 156·6 km² drainage basin of Pancho Rico Creek (PRC), and sediment transport from the PRC drainage basin to its c. 11 000 km² mainstem drainage (Salinas River) are investigated. Numeric age estimates are determined by optically stimulated luminescence (OSL) dating on quartz grains from three sediment samples collected from a ‘quaternary terrace a (Qta)’ PRC terrace/PRC‐tributary fan sequence, which consists dominantly of debris flow deposits overlying fluvial sediments. OSL dating results, morphometric analyses of topography, and field results indicate that the stormy climate of the Pleistocene‐Holocene transition caused intense debris‐flow erosion of PRC‐tributary valleys. However, during that time, the PRC channel was backfilled by Qta sediment, which indicates that there was insufficient discharge in PRC to transport the sediment load produced by tributary‐valley denudation. Locally, Salinas Valley alluvial stratigraphy lacks any record of hillslope erosion occurring during the Pleistocene‐Holocene transition, in that the alluvial fan formed where PRC enters the Salinas Valley lacks lobes correlative to Qta. This indicates that sediment stripped from PRC tributaries was mostly trapped in Pancho Rico Valley despite the relatively moist climate of the Pleistocene‐Holocene transition. Incision into Qta did not occur until PRC enlarged its drainage basin by c. 50% through capture of the upper part of San Lorenzo Creek, which occurred some time after the Pleistocene‐Holocene transition. During the relatively dry Holocene, PRC incision through Qta and into bedrock, as well as delivery of sediment to the San Ardo Fan, were facilitated by the discharge increase associated with stream‐capture. The influence of multiple mechanisms on sediment storage and transport in the Pancho Rico Valley‐Salinas Valley system exemplifies the complexity that (in some instances) must be recognized in order to correctly interpret terrestrial sedimentary sequences in tectonically active areas. Copyright © 2009 John Wiley & Sons, Ltd.     

The effects of deforestation on slope and channel evolution in the tectonically active Darjeeling Himalaya

The main indicators of Quaternary tectonic uplift are the young mountain slopes of the Darjeeling Himalaya, rising straight above the Ganga–Brahmaputra foredeep, fragments of uplifted river terraces and fresh fault scarps. Evidence for the continuation of the uplift includes downcutting of the Tista and other straight rivers in the bedrock, continuing aggradation in the plains and overriding of the metamorphic rocks on the alluvia. Owing to deforestation and extensive land use, the earlier natural tendency of a dominance of channel incision over slope degradation has changed to prevailing aggradation, even in steep valley reaches, caused by intensive slope mass movements and the overloading of the mountain creeks. Aggradation progresses upstream along the rivers dissecting the mountain front.     

After horton

The divergent and yet related problems of post-Hortonian studies of drainage density and channel network geometry are viewed against the difficulties of defining first-order channels and basins. It is proposed that the junction of an unbranched perennial (or blue-line) channel with another perennial channel be taken as the starting point for definitions and that the entire contour-crenulation network tributary to that point be considered the first-order stream. It is shown that the concept of network diameter may be used to describe the networks so delimited and that it appears to provide a useful starting point for interregional comparisons. Finally, an analysis of Blyth and Rodda's (1973) data on channel lengths and discharge indicates that network diameter may be as closely related to discharge as is channel length itself.     

Drainage density on progressively tilted surfaces with different gradients,Wheeler Ridge,California

Wheeler Ridge in the Southern San Joaquin Valley, California, is an anticlinal fold which has been progressively uplifted during the last 250 ka. Drainage networks on the ridge become younger as the anticline's eastern tip is approached. Because of the fold's asymmetric shape, surfaces on opposite flanks of the ridge have similar ages but very different gradients. The ridge provides important insights into drainage development on progressively tilted surfaces, as existing studies are restricted to static topography. A surface gradient of between 4·8° and c.10° is needed to initiate channel networks. This gradient threshold is consistent with previous studies of the gradient and upslope area needed to incise a channel through overland flow. Comparison of coeval drainage networks on opposite flanks of the ridge allows the controversial relationship between drainage density and gradient to be investigated. A lower valley density is observed on the higher gradient flank of Wheeler Ridge. Field observations from the ridge indicate that this inverse relationship is associated with hillslope erosion by shallow mass-wasting, the rate of which increases rapidly as a threshold gradient is approached. Comparison of data from Wheeler Ridge with other field studies and numerical models, shows that the form of the relationship between gradient and drainage density is process-dependent. A positive correlation occurs when erosion is a result of overland flow, whilst a negative correlation occurs where erosion is dominated by shallow mass-wasting. Copyright © 1999 John Wiley & Sons, Ltd.     

Morphology of the Itapeva to Tramandai transgressive dunefield barrier system and mid‐ to late Holocene sea level change

The surface morphology of the transgressive dunefield barrier extending from Itapeva to Tramandai along the northern littoral of Rio Grande do Sul, is examined and an attempt is presented to link morphological changes across the barrier to Holocene sea level changes. The 4·5 km wide Holocene barrier displays two typical morphologies: an inner part dominated by large‐scale, continuous alongshore, overlapping dunefield phases comprising sand sheets, dunes, deflation plains and precipitation ridges; and an outer part dominated by discontinuous, medium‐ to small‐scale, triangular to lobate transgressive dunefield phases, cut by both active and relict (palaeo‐) creeks or washouts. Holocene sea level in the region rose to c. +1 to +3 m above present reaching a maximum around 5100 years bp and then progressively fell to the present level. We argue that the effect on barrier development was to suppress the development of a drainage system during the rising and maximum stages, and encourage the development of an organized drainage system in the form of regularly spaced washouts during the falling period, and that this change in sea level from rising to falling therefore produced the large‐scale differences in barrier morphology. Copyright © 2006 John Wiley & Sons, Ltd.     

Macroinvertebrate communities along the main stem and tributaries of a pre-Alpine river: composition responds to altitude,richness does not

We collected quantitative macroinvertebrate samples and measured environmental and geographical parameters at 13 sites: six along the main stem and seven in tributaries close to the main channel over a 700 m gradient in altitude and 22 km longitudinal distance along the River Kokra in the Slovenian Alps. Our objectives were 1) to compare longitudinal patterns in richness and community composition between main stem and tributary sites, and 2) to determine the relative importance of the replacement and richness difference component for overall beta diversity and of environmental versus spatial distance on beta diversity among main stem and tributary sites.In total 138 taxa were identified. There were no differences between main stem and tributary sites in mean abundance or taxon richness (67 and 58, respectively). A nMDS and ANOSIM based on Bray-Curtis similarities found no separation of main stem and tributary sites, but that upper (≥880 m a.s.l) and lower sites (≤680 m a.s.l.) formed two different groups. In both main stem and tributaries taxon richness increased only slightly going downstream while the community composition (DCA1) was much better explained by altitude and distance from source.Overall, beta diversity (Sørensen and Bray-Curtis dissimilarity) was similar for the two groups, and total Sørensen dissimilarity was driven mainly by replacement in main stem (78 %) and tributary sites (77 %). Mantel tests showed that main stem dissimilarities were significantly correlated to environmental PCA distance, watercourse distance, overland distance and altitudinal differences. Tributary dissimilarities were not correlated to any of these four factors. GLMs showed that dissimilarity among main stem sites was explained only by altitude difference, while no factors were significant among tributary sites, even though nearly so for environmental PCA distance.The study illustrates the importance of measuring beta diversity along ecological gradients, such as river continua and/or altitudinal gradients, where alpha diversity may fail to detect relatively minor changes in assemblage composition. Such changes are likely to occur due to present and future climate warming.