Relationships between landscape morphology,climate and surface erosion in northern Peru at 5°S latitude

The northern segment of the Peruvian Andes is affected by a twofold climate with measurable implications on landscapes and landscape dynamics. During ‘normal’ or ‘neutral’ years easterly winds bring rain from the Atlantic and the Amazon Basin to the Sierras, which results in a seasonal climate with rather low-intensity precipitations. In contrast, during the large-scale warm phase of the ENSO cycle, El Niños transfer moisture from the Pacific to the Peruvian coast by westerly winds and result in high-intensity precipitation. We investigate the effects of this twofold climate for the case of the Piura drainage basin at ca. 5°S latitude (northern Peru). In the headwaters that have been under the influence of the easterlies, the landscape is mantled by a thick regolith cover and dissected by a network of debris flow channels that are mostly covered by a thick layer of unconsolidated sediment. This implies that in the headwaters of the Piura River sediment discharge has been limited by the transport capacity of the sediment transfer system. In the lower segment that has been affected by high-intensity rainfall in relation to the westerlies (El Niños), the hillslopes are dissected by debris flow channels that expose the bedrock on the channel floor, implying a supply-limited sediment discharge. Interestingly, measurements at the Piura gauging station near the coast reveal that, during the last decades, sediment was transferred to the lower reaches only in response to the 1982–1983 and 1997–1998 El Niño periods. For the latter period, synthetic aperture radar (SAR) intensity images show that the locations of substantial erosion are mainly located in areas that were affected by higher-than-average precipitation rates. Most important, these locations are coupled with the network of debris flow channels. This implies that the seasonal easterlies are responsible for the production of sediment through weathering in the headwaters, and the highly episodic El Niños result in export of sediment through channelized sediment transport down to the coastal segment. Both systems overlap showing a partially coupled sediment production–delivery system.     

Limits of sediment transfer in an alpine debris-flow catchment,Illgraben, Switzerland

The 9.5 km² Illgraben catchment, located in the Rhône valley in the Central Alps of Switzerland, is one of the most active debris flow torrents in the Alps. In this paper we present sediment yield data collected in 2006 for segments where hillslopes and channels form a fully connected network and contrast these with sediment yields measured for disconnected hillslopes. The data reveal that sediment yields are 1–2 orders of magnitude larger in segments where hillslopes are connected with the channel network than on disconnected hillslopes. Support for this conclusion is provided by observations made on 1959, 1999 and 2004 aerial photographs that the vegetation cover in the disconnected segments is still intact, whereas denudation rates of several centimeters per year in the connected segments have inhibited the establishment of a stable vegetation cover. Furthermore, sediment supplied from hillslopes during the past 40 years has temporarily accumulated along the Illgraben channel, indicating that the channel aggraded over this period and has not yet recovered. An implication of this observation is that initiation of debris flows in the Illgraben catchment is limited more by the availability of intense rainfall than sediment. In contrast, on disconnected hillslopes, sediment flux does not appear to be driven by precipitation.The petrographic composition of the Illgraben fan deposits indicates two distinct sediment sources, one related to rockfall and the other to landslides and debris flows. The presence of clasts from both sources implies multiple processes of erosion, deposition, mixing and re-entrainment in the catchment before the material is exported to the Illgraben fan and to the Rhône River. In addition, delivery of large amounts of coarse-grained sediment to the Rhône causes a modification of the flow pattern from meandering or anastomosing upstream to braided downstream. Hence, the direct connectivity between hillslope and channelized processes in the Illgraben catchment causes not only rapid topographic modifications in the catchment, but also morphologic adjustment in the Rhône valley downstream.     

山坡表层关键带结构与水文连通性研究进展

山丘区是洪水的"策源地",山丘区坡地、沟谷及间歇性河道为洪水的形成提供了通道,同时也是水文连通时空变化最为强烈的地带。然而,对流域表层关键带结构特征及其水文连通机制等的认识尚存不足,限制了产汇流理论及模型方法的发展和应用。通过对比国内外山坡水文实验,发现山坡物理结构连通性控制并深刻影响着水流的连通过程,现有水文连通实验侧重孔隙等微观尺度的规律研究,与水文模型理论存在尺度上的巨大偏差。为此,提出水文连通性应侧重揭示水流在山坡地表、地下的宏观表象通道及分布特征,探索径流连通的动力学机制,即山坡水文连通性研究重在剖析其结构特征的水文累积效应,应保持关键带结构特征合理概化与产汇流理论适度复杂之间的平衡。     

Contrasting sediment flux in Val Lumnezia (Graubünden, Eastern Swiss Alps), and implications for landscape development

This paper presents qualitative estimates of sediment discharge from opposite valley flanks in the S–N-oriented Val Lumnezia, eastern Swiss Alps, and relates inferred differences in sediment flux to the litho-tectonic architecture of bedrock. The valley flank on the western side hosts the deep-seated Lumnezia landslide where an area of ca. 30 km² has experienced slip rates of several centimetres per year, potentially resulting in high sediment discharge to the trunk stream (i.e. the Glogn River). High slip rates have resulted in topographic changes that are detectable on aerial photographs and measurable with geodetic tools. In contrast, a network of tributary channels dissects the valley flank on the eastern side. There, an area of approximately 18 km² corresponding to < 30% of the surface has experienced a change in the landscape mainly by rock avalanche and rock fall, and the magnitudes of changes are below the calibration limit of digital photogrammetry. We thus infer lower magnitudes of sediment discharge on the eastern tributaries than on the western valley side, where landsliding has been the predominant erosional process. These differences are interpreted to be controlled by the dip-slope situation of bedrock on the western side that favours down-slope slip of material. Morphometric investigations reveal that the western valley side is characterized by a low topographic roughness because this valley flank has not been dissected by a channel network. It appears that high sediment discharge of the Lumnezia landslide has inhibited the establishment of a stable channel network and has largely controlled the overall evolution of the landscape. This contrasts to the general notion that channelized processes exert the first-order control on landscape evolution and formation of relief and needs to be considered in future studies about landscape architecture, drainage network and sediment discharge.     

Hillslope coupled stream morphology,flow conditions,and their effects on detrital sedimentology in Garnet Canyon,Teton Range,Wyoming

Characterizing stream erosion in any steep mountain landscape is arduous, but the challenge level increases when the stream flows through a glaciated catchment frequently modified by hillslope debris.Glacial landforms and stochastic mass wasting in alpine systems may interfere with sediment delivery to downstream sites where detrital sediments are often collected to represent upstream bedrock sources.To use detrital sediments as indicators of erosion, we need to understand potential sediment accumulation in flat glaciated reaches or behind rockfall barriers. This study investigates the stream channel in Garnet Canyon, a glaciated catchment located in the central Teton Range, to describe hillslope coupled channel morphology and the subsequent effects on sediment transport throughout the catchment.Stream cross-section surveys and sediment size measurements of the surface bedload were collected in the field within a glacially flattened segment of Garnet Canyon. Calculations of shear stress conditions allowed evaluation of the importance of mineral densities on potential grain entrainment. The length of the Garnet Canyon stream observed in this study was coupled with hillslope deposits. Critical shear stresses were sufficient to move gravel-sized sediments through all sections when calculated with quartz mineral density and through most sections when applying apatite mineral density. These results verify the application of detrital sediments to evaluate erosion rates or spatial bedrock sources because snowmelt stream flow efficiently moves entrained sediment past glacially reduced slopes and potential talus barriers.     

Pleistocene geomorphology and geochronology of eastern Grand Canyon: linkages of landscape components during climate changes

We report new mapping, soils, survey, and geochronologic (luminescence, U-series, and cosmogenic-nuclide) data from Pleistocene deposits in the arid setting of eastern Grand Canyon. The result is a stratigraphic framework of inset fill gravels and associated terraces that provide a record of the responses of hillslopes, tributary streams, and the Colorado River to the last 400 kyr of glacial–interglacial climate change. The best-preserved last 80 kyr of this record indicates a stratigraphic–chronologic disconnect between both deposition and incision along the Colorado River versus along the trunks of local tributaries. For example, the Colorado River finished aggrading and had already begun incising before the main pulse of aggradation in the trunks of local catchments during Marine Isotope Stage 3, and then tributary incision followed during the millennial-scale fluctuations of the last glacial epoch, potentially concurrent with mainstem aggradation. The mainstem record appears to broadly correlate with regional paleoclimate and upstream geomorphic records and thus may be responding to climatic–hydrologic changes in its mountain headwaters, with aggradation beginning during full-glacial times and continuing into subsequent interglacials. The contrasting lag time in responses of the dryland catchments within Grand Canyon may be largely a function of the weathering-limited nature of hillslope sediment supply.     

Amalgamation surfaces, bed thicknesses, and dish structures in sand-rich submarine fans: numeric differences in channelized and unchannelized deposits and their diagnostic value

Hydraulic differences between channelized and unchannelized flows in sand-rich submarine fans result in different distributions of amalgamation surfaces, bed thicknesses, and dish structures in successions of these two different environments. Distribution trends of these fabrics were quantified for the sand-rich fans of the Reiselsberger Sandstein (Cenomanian–Turonian). These trends can be used as criteria to distinguish channelized from unchannelized paleoenvironments of sand-rich submarine fans.

Amalgamation surfaces in the studied fans' channelized regions are considerably more abundant than in the unchannelized fan areas. In unchannelized deposits, tabular amalgamation surfaces outnumber nontabular ones, whereas the opposite occurs in channelized successions. These results indicate a higher degree of erosive power of gravity-driven sediment flows in channels as a result of a greater flow thickness, higher flow velocity, and turbulence.

The average turbidite layer thickness in channelized successions is markedly greater than in unchannelized deposits (“layer” as defined herein). This is mainly attributed to the combined effects of differences in sediment fall-out rate and the inefficiency of sand-rich suspensions to transport sand. In the proximal and channelized fan areas, more sediment is deposited from a flow in the form of a layer than in distal unchannelized fan regions despite a higher degree of erosion in channels. The greater average bed thickness in channel fills is a function of layer thickness and more frequent amalgamations (“bed” as defined herein).

Dish structures seem to be considerably more common in midfan than in outer-fan successions. This may indicate a higher sedimentation rate from individual suspension currents in midfan areas.     

黄土丘陵区土壤侵蚀链垂直带水沙流空间分布

采用多坡段组合模型、人工模拟降雨实验研究了土壤侵蚀链垂直带水沙流空间分布变化特征。结果表明:愈向下坡方向的坡段产流量愈大,单位面积、单位时间的产流量按坡段排列为:谷坡>梁峁坡下部>梁峁坡中部>梁峁坡上部。雨强为29.7mm/h时,坡面上没有沟蚀发育;雨强为60.5mm/h时,细沟主要分布在梁峁坡下部和谷坡处;雨强为90 2mm/h时,各种侵蚀形态在坡面均有较好发育,细沟的出现部位一直伸展到梁峁坡中部和上部之间。由于上坡来水来沙的作用,梁峁坡的产沙量增大了20.2%～63.5%,谷坡的产沙量增大了42.9%～74.5%。     

Triggering upstream channel incision in the fluvial systems of the trunk valley of Wadi Watir drainage basin, South Sinai Peninsula, Egypt

Dredging the alluvial fans for repaving the international road located in the bottom of the Wadi Watir valley produced vertical cliff faces of different heights, and at different locations of the fans. The heights of the cliff faces resulted in considerable elevation differences between the surface of the dredged alluvial fans and the local base level provided by the Watir trunk valley. The principal geomorphic response to this anthropogenic intervention is triggering upstream channel incision waves at different intensities in the fluvial systems of the downstream reaches of the Watir drainage basin. The channel incision processes resulted in subsequent geomorphic adjustment scenarios that vary from widening the active channels on the surface of the dredged fans, triggering rockfalls from the adjacent hillslopes, and transporting coarse alluvial deposits from the main sediment sources of the fluvial systems, and eventually re-depositing them as sheetform gravel, channelform gravel, and new fan lobes. The major outcome of the various geomorphic adjustment processes was changing the role of the alluvial fans within the fluvial systems from buffer zones where fan aggradation was dominant into dynamic coupled zones. Being coupled zones, the dredged alluvial fans allowed high potential of mass transmission from the feeder catchment areas into the Watir trunk valley. Under such conditions, it could be stipulated that considerable changes in the morphology of landscapes are highly anticipated in response to flash flood events that intermittently occur in the Watir drainage basin.     

The use of landscape evolution models in mining rehabilitation design

Landscape evolution models can be useful tools for the evaluation of rehabilitation designs for post-mining landscapes. When calibrated for the erodible material, landscape evolution models can predict sediment loss over entire landscapes (i.e. tonnes/hectare/year), method of erosion (i.e. slope wash, gullying) and also where on a hillslope erosion is likely to occur. The models provide the ability to examine simple hillslopes through to complex whole landscapes. These models can also be used for a probabilistic risk assessment of rehabilitation design for high-risk situations such as tailings dams. Importantly, unlike other erosion models they allow the eroded landscape to be visualised. This paper outlines the capabilities of the SIBERIA landscape evolution model for the rehabilitation of mining landscapes and proposes a probabilistic approach for risk assessment and site stability.     

The influence of shallow landslides on sediment supply: A flume-based investigation using sandy soil

The impact of rainfall-induced shallow landslides on hillslope sediment discharge is not well understood. This paper reports experimental measurements of sediment discharge after water-induced shallow landslides are triggered on sandy soil in a flume under simulated rainfall. The principal aim of the research was to investigate how varying soil depth affects the location and occurrence of shallow slope failures, as well as how it affects sediment yields downslope. Four experiments were conducted using the same sandy soil and a 30° and 10° compound slope configuration under average rainfall intensity of 50 mm h^− 1 for up to 390 min. Soil depths were set to 200, 300, 400 and 500 mm. Engineering and geotechnical properties of the soil were examined. Sediment discharge and runoff were collected from the flume outlet at 15 minute intervals. Changes in the soil slope profiles after landslides and soil physical properties resulted from soil armouring, under continuous rainfall were also recorded. Results showed that sediment yields at the flume outlet, before landslides occurred, were very low and limited to the finer soil particles as would be expected for a sandy soil. However subsequent variations in sediment discharge were strongly related to failure events and their proximity to the outlet. Sediment yield was also affected by the original soil depth; the greater the depth, the higher the sediment yields. Post-failure reductions in sediment discharge were observed and attributed to post-failure slope stabilization under continuing rainfall and extensive soil armouring near the flume outlet. The results provide a clear linkage between landslides and sediment discharge due to hydrological processes occurring in the hillslope. This knowledge is being used to develop a model to predict sediment discharges from hillslopes following shallow landslide events.     

Fluvial to tidal transition in proximal,mixed tide‐influenced and wave‐influenced deltaic deposits: Cretaceous lower Sego Sandstone,Utah, USA

Facies models for regressive, tide‐influenced deltaic systems are under‐represented in the literature compared with their fluvial‐dominated and wave‐dominated counterparts. Here, a facies model is presented of the mixed, tide‐influenced and wave‐influenced deltaic strata of the Sego Sandstone, which was deposited in the Western Interior Seaway of North America during the Late Cretaceous. Previous work on the Sego Sandstone has focused on the medial to distal parts of the outcrop belt where tides and waves interact. This study focuses on the proximal outcrop belt, in which fluvial and tidal processes interact. Five facies associations are recognized. Bioturbated mudstones (Facies Association 1) were deposited in an offshore environment and are gradationally overlain by hummocky cross‐stratified sandstones (Facies Association 2) deposited in a wave‐dominated lower shoreface environment. These facies associations are erosionally overlain by tide‐dominated cross‐bedded sandstones (Facies Association 4) interbedded with ripple cross‐laminated heterolithic sandstones (Facies Association 3) and channelized mudstones (Facies Association 5). Palaeocurrent directions derived from cross‐bedding indicate bidirectional currents which are flood‐dominated in the lower part of the studied interval and become increasingly ebb‐directed/fluvial‐directed upward. At the top of the succession, ebb‐dominated/fluvial‐dominated, high relief, narrow channel forms are present, which are interpreted as distributary channels. When distributary channels are abandoned they effectively become estuaries with landward sediment transport and fining trends. These estuaries have sandstones of Facies Association 4 at their mouth and fine landward through heterolithic sandstones of Facies Association 3 to channelized mudstones of Facies Association 5. Therefore, the complex distribution of relatively mud‐rich and sand‐rich deposits in the tide‐dominated part of the lower Sego Sandstone is attributed to the avulsion history of active fluvial distributaries, in response to a subtly expressed allogenic change in sediment supply and relative sea‐level controls and autocyclic delta lobe abandonment.     

山坡地形曲率分布特征及其水文效应分析——真实流域的野外实验及相关分析研究

研究了山坡地形曲率的空间分布特征,发现流域略呈现出微凹、发散的形状特征,平均的剖面和水平曲率分别为-4.62×10-4和3.49×10-4。通过划分源头型和边坡型山坡,发现源头山坡多是收敛的,边坡型山坡则多为发散的;在山坡内部收敛、发散、凹和凸等地形地貌类型是同时存在的,即山坡一般由相互组合的各种坡型所组成。基于野外采样观测的成果,分析了和睦桥子流域90个采样点的土壤含水量及其对应曲率的相互关系。结果显示:一般来说,凹形山坡、收敛形山坡对应的土壤含水量更高,而凸形和发散形山坡对应的土壤含水量较低;不论是收敛的凹坡,还是发散的凸坡,其采样点土壤含水量都有随高程下降而升高的趋势,即愈靠近沟谷土壤含水量愈大。     

Trench-slope channels from the New Zealand Jurassic: the Otekura Formation, Sandy Bay, South Otago

The Otekura Formation (Early Jurassic, Pseudaucella zone) at Sandy Bay comprises part of a 10+ km thick, regressive, forearc shelf and slope sequence, the Hokonui facies belt of the Rangitata Geosyncline. The Otekura Formation is dominantly fine grained, being mostly mudstone, silty mudstone and siltstone. The sediments are volcanogenic throughout. The upper 150 m of the formation contains two 20 m thick, channelized bodies of medium-thick bedded sandy flysch, each associated with thin bedded muddy flysch interpreted as overbank turbidites. Directional indicators within the channel sequence indicate emplacement from the south-southwest. In contrast, rare turbidites that occur below the channel sequence, within the background mudstone sediment, were emplaced from the east, i.e. at right angles to the channelized flows. The immediately overlying Omaru Formation contains more abundant macrofossils, intraclastic conglomerates, and appreciable amounts of traction-emplaced cross-bedded sand. Bioturbated calcareous siltstones with an in situ molluscan fauna follow (Boatlanding Formation), and are of shelf origin. The Omaru Formation is therefore interpreted as a shelf-slope break deposit, and the Otekura Formation as an upper slope facies. Reconnaissance studies indicate that the Otekura Formation is underlain by several kilometres of dominantly fine grained, deep water slope sediments, containing occasional sand and conglomerate filled channels similar to those here described in detail from the Otekura Formation. Such channels are inferred to form when a mass-transported sand, derived from failure higher on the slope, ploughs erosively into the sea floor. After their incision, the channels served for a short time as conduits for downslope transport of sediment, the redeposited deposits of which are found filling each channel. Both channel fills at Sandy Bay are capped by thin-bedded turbidites inferred to have overspilled from similar channels nearby on the slope.     

A Markov model for the description and classification of hillslopes

If considered as a sequence of morphological units, hillslopes can be treated as realizations of a stochastic process in which morphological units correspond to events and their attributes (e.g., slope)to values. A first-order, discrete-state, discrete-time Markov model, applied to a group of hillslopes measured in Calabria, explained a significant part of the hillslope morphology. A Markov property was determined to be present in the sequence of slopes and lengths of the morphological units. In the subsequent analysis of the Markov property in the sequence of slopes, an optimum transition scheme was developed. This indicated, inter alia,that concavity in the hillslope was most probable where the slope was in the range from 7° to <9°. The hillslopes were divided into classes based upon their geomorphological position. The resulting transition probability matrices and their derivatives were examined for significant differences.     

A structured approach to enhance flood hazard assessment in mountain streams

An evidence-based flood hazard analysis in mountain streams requires the identification and the quantitative characterisation of multiple possible processes. These processes result from specific triggering mechanisms on the hillslopes (i.e. landslides, debris flows), in-channel morphodynamic processes associated with sudden bed changes and stochastic processes taking place at critical stream configurations (e.g. occlusion of bridges, failure of levees). From a hazard assessment perspective, such possible processes are related to considerable uncertainties underlying the hydrological cause-effect chains. Overcoming these uncertainties still remains a major challenge in hazard and risk assessment and represents a necessary condition for a reliable spatial representation of process intensities and the associated probabilities. As a result of an accurate analysis of the conceptual flaws present in the procedures currently employed for hazard mapping in South Tyrol (Italy) and Carinthia (Austria), we propose a structured approach as a means to enhance the integration of hillslope, morphodynamic and stochastic processes into conventional flood hazard prediction for mountain basins. To this aim, a functional distinction is introduced between prevailing one-dimensional and two-dimensional process propagation domains, i.e., between confined and semi- to unconfined stream segments. The former domains are mostly responsible for the generation of water, sediment and wood fluxes, and the latter are where flooding of inactive channel areas (i.e. alluvial fans, floodplains) can occur. For the 1D process propagation domain, we discuss how to carry out a process routing along the stream system and how to integrate numerical models output with expert judgement in order to derive consistent event scenarios, thus providing a consistent quantification of the input variables needed for the associated 2D domains. Within these latter domains, two main types of spatial sub-domains can be identified based on the predictability of their dynamics, i.e., stochastic and quasi-deterministic. Advantages and limitations offered by this methodology are finally discussed with respect to hazard and risk assessment in mountain basins.     

Multi‐scale analysis of a migrating submarine channel system in a tectonically‐confined basin: The Miocene Gorgoglione Flysch Formation,southern Italy

The Miocene Gorgoglione Flysch Formation records the stratigraphic product of protracted sediment transfer and deposition through a long‐lived submarine channel system developed in a narrow and elongate thrust‐top basin of the Southern Apennines (Italy). Channel‐fill deposits are exposed in an outcrop belt approximately 500 m thick and 15 km long, oriented oblique to the palaeoflow, which was roughly south‐eastward. These exceptional exposures of channel‐fill strata allow the stacking architectures and the evolution of the channel system to be analyzed at multiple scales, enabling the effects of syn‐sedimentary thrust tectonics and basin confinement on the depositional system development to be deciphered. Two end‐member types of elementary channel architecture have been identified: high‐aspect‐ratio, weakly‐confined channels, and low‐aspect‐ratio, incisional channels. Their systematic stacking results in a complex pattern of seismic‐scale depositional architectures that determines the stratigraphic framework of the deep‐water system. From the base of the succession, two prominent channel complex sets have been recognized, namely CS1 and CS2, consisting of amalgamated incisional channel elements and weakly‐confined channel elements. These channelized units are overlain by isolated incisional channels, erosional into mud‐prone slope deposits. The juxtaposition of different channel architectures is interpreted to have been governed by regional thrust‐tectonics, in combination with a high subsidence rate that promoted significant aggradation. In this scenario, the alternating ‘in sequence’ and ‘out of sequence’ tectonic pulses of the basin‐bounding thrusts controlled the activation of coarse‐clastic inputs in the basin and the resulting stacking architectures of channelized units. The tectonically‐driven confinement of the depositional system limited the lateral offset in channel stacking, preventing large‐scale avulsions. This study represents an excellent opportunity to analyze the stratigraphic evolution of a submarine channel system in tectonically‐active settings from an outcrop perspective. It should find wide applicability in analogous depositional systems, whose stratigraphic architecture has been influenced by tectonically‐controlled lateral confinement and associated lateral tilting.     

Pathways for sediment transport and retention in a vegetated,mid-channel island: Connecting sediment dynamics to morphology in Meinmahla Island,Ayeyarwady Delta,Myanmar

Vegetated mid-channel islands play an important though poorly understood role in the sediment dynamics and morphology of tide-dominated deltas. Meinmahla Island is a mangrove-forest preserve at the mouth of the Bogale distributary channel, in the Ayeyarwady Delta, Myanmar. In this relatively unaltered mid-channel island, sediment dynamics can be directly connected to morphology. Field measurements from 2017 to 2019 provide insight into the pathways for sediment transport and resulting morphological evolution. Water depth, salinity and turbidity were monitored semi-continuously, and velocity profilers with turbidity and salinity sensors were deployed seasonally in single-entrance (dead-end/blind) and multi-entrance tidal channels of the island. The morphological evolution was evaluated using grain size, ²¹⁰Pb geochronology, remote sensing and channel surveys. The data show that ebb-dominant, single-entrance channels along the island exterior import sediment year-round to the land surface. However, these exterior channels do not deliver enough sediment to maintain the observed ca 0.8 cm/yr accretion rate, and most of the sediment import occurs via interior, multi-entrance channels. Interior channels retain water masses that are physically distinct from the water in the Bogale distributary, and estuarine processes at the tidal-channel mouths import sediment into the island. Sediment is sourced to the island from upriver in the wet season and from the Gulf of Mottoma in the dry season, as the location of the estuary shifts seasonally within the Bogale distributary. The salinity and biogeochemistry of the distributary water are affected by interactions with sediment and groundwater in the island interior. The largest interior channels have remained remarkably stable while the island has aggraded and prograded over decadal timescales. However, the studied multi-entrance channel is responding to a drainage-network change by narrowing and shoaling. Overall, mid-channel islands trap sediment and associated nutrients at the river–ocean interface, and these resilient landscape features evolve in response to changes in drainage-network connectivity.     

Storage of sediment-associated nutrients and contaminants in river channel and floodplain systems

Samples of fine-grained channel bed sediment and overbank floodplain deposits were collected along the main channels of the Rivers Aire (and its main tributary, the River Calder) and Swale, in Yorkshire, UK, in order to investigate downstream changes in the storage and deposition of heavy metals (Cr, Cu, Pb, Zn), total P and the sum of selected PCB congeners, and to estimate the total storage of these contaminants within the main channels and floodplains of these river systems. Downstream trends in the contaminant content of the <63 μm fraction of channel bed and floodplain sediment in the study rivers are controlled mainly by the location of the main sources of the contaminants, which varies between rivers. In the Rivers Aire and Calder, the contaminant content of the <63 μm fraction of channel bed and floodplain sediment generally increases in a downstream direction, reflecting the location of the main urban and industrialized areas in the middle and lower parts of the basin. In the River Swale, the concentrations of most of the contaminants examined are approximately constant along the length of the river, due to the relatively unpolluted nature of this river. However, the Pb and Zn content of fine channel bed sediment decreases downstream, due to the location of historic metal mines in the headwaters of this river, and the effect of downstream dilution with uncontaminated sediment. The magnitude and spatial variation of contaminant storage and deposition on channel beds and floodplains are also controlled by the amount of <63 μm sediment stored on the channel bed and deposited on the floodplain during overbank events. Consequently, contaminant deposition and storage are strongly influenced by the surface area of the floodplain and channel bed. Contaminant storage on the channel beds of the study rivers is, therefore, generally greatest in the middle and lower reaches of the rivers, since channel width increases downstream. Comparisons of the estimates of total storage of specific contaminants on the channel beds of the main channel systems of the study rivers with the annual contaminant flux at the catchment outlets indicate that channel storage represents <3% of the outlet flux and is, therefore, of limited importance in regulating that flux. Similar comparisons between the annual deposition flux of specific contaminants to the floodplains of the study rivers and the annual contaminant flux at the catchment outlet, emphasise the potential importance of floodplain deposition as a conveyance loss. In the case of the River Aire the floodplain deposition flux is equivalent to between ca. 2% (PCBs) and 36% (Pb) of the outlet flux. With the exception of PCBs, for which the value is ≌0, the equivalent values for the River Swale range between 18% (P) and 95% (Pb). The study emphasises that knowledge of the fine-grained sediment delivery system operating in a river basin is an essential prerequisite for understanding the transport and storage of sediment-associated contaminants in river systems and that conveyance losses associated with floodplain deposition exert an important control on downstream contaminant fluxes and the fate of such contaminants.