Unravelling the conundrum of river response to rising sea‐level from laboratory to field. Part I: Laboratory experiments

The most recent deglaciation resulted in a global sea‐level rise of some 120 m over approximately 12 000 years. In this Part I of two parts, a moving boundary numerical model is developed to predict the response of rivers to this rise. The model was motivated by experiments at small scale, which have identified two modes describing the transgression of a river mouth: autoretreat without abandonment of the river delta (no sediment starvation at the topset–foreset break) and sediment‐starved autoretreat with abandonment of the delta. In the latter case, transgression is far more rapid and its effects are felt much further upstream of the river mouth. The moving boundary numerical model is checked against experiments. The generally favourable results of the check motivate adaptation of the model to describe the response of the much larger Fly‐Strickland River system, Papua New Guinea to Holocene sea‐level rise; this is done in the companion paper, Part II.     

Sedimentary facies and geomorphic evolution of a blocked‐valley lake: Lake Futululu,northern Kwazulu‐Natal,South Africa

Blocked‐valley lakes are formed when tributaries are impounded by the relatively rapid aggradation of a large river and its floodplain. These features are common in the landscape, and have been identified in the floodplains of the Solimões‐Amazon (Brazil) and Fly‐Strickland Rivers (Papua New Guinea), for example, but their inaccessibility has resulted in studies being limited to remotely sensed image analysis. This paper documents the sedimentology and geomorphic evolution of a blocked‐valley lake, Lake Futululu on the Mfolozi River floodplain margin, in South Africa, while also offering a context for the formation of lakes and wetlands at tributary junctions. The study combines aerial photography, elevation data from orthophotographs and field survey, and longitudinal sedimentology determined from a series of cores, which were sub‐sampled for organic content and particle size analysis. Radiocarbon dating was used to gauge the rate and timing of peat accumulation. Results indicate that following the last glacial maximum, rising sea‐levels caused aggradation of the Mfolozi River floodplain. By 3980 years bp , aggradation on the floodplain had impounded the Futululu drainage line, creating conditions suitable for peat formation, which has since occurred at a constant average rate of 0·13 cm year^?1. Continued aggradation on the Mfolozi River floodplain has raised the base level of the Futululu drainage line, resulting in a series of back‐stepping sedimentary facies with fluvially derived sand and silt episodically prograding over lacustrine peat deposits. Blocked‐valley lakes form where the trunk river has a much larger sediment load and catchment than the tributary stream. Similarly, when the relative difference in sediment loads is less, palustrine wetlands, rather than lakes, may be the result. In contrast, where tributaries drain a steep, well‐connected catchment, they may impound much larger trunk rivers, creating lakes or wetlands upstream.     

Millennial‐scale stratigraphy of a tide‐dominated incised valley during the last 14 kyr: Spatial and quantitative reconstruction in the Tokyo Lowland,central Japan

Spatial and quantitative analysis of infilling processes of the tide‐dominated incised valleys beneath the Tokyo Lowland during the last 14 kyr was undertaken by using data from 18 sediment cores, 467 radiocarbon dates and 6100 borehole logs. The post‐Last Glacial Maximum valley fills consist of braided river, meandering river, estuary, spit and delta systems in ascending order. The boundary between the estuary and delta systems is regarded as the maximum flooding surface. The maximum flooding surface beneath the Tokyo Lowland is dated at 8 ka in the Arakawa Valley and 7 ka in the Nakagawa Valley. This age difference is due to the migration of the Tone River from the Arakawa Valley to the Nakagawa Valley at 5 ka, and suggests that the widely held view that the global initiation of deltas coincided with the abrupt rise of sea‐level at 9 to 8 ka is true only where there has been steady sediment supply from major rivers. The meandering river system is dominated by sheet‐like sands that were deposited during lateral migration of channels during the Younger Dryas and isolated vertical sands within muds that reflect vertical aggradation of channels before and after the Younger Dryas. The transition between these channel geometries is controlled by a threshold sea‐level rise of 4 to 7 mm yr^?1. Before migration of the Tone River at 5 ka, the tide‐dominated bay in the Nakagawa Valley was filled by upward‐fining laterally accreting muds. The muds accreted from the margin to the axis of the bay. Such lateral accretion of suspended particles derived from outside the bay has been documented in other tide‐dominated coastal environments and is probably common in other similar settings. After the migration of the Tone River, the bay was filled by upward‐coarsening deltaic sediments.     

Conceptual framework for assessing the response of delta channel networks to Holocene sea level rise

Recent research has identified two fundamental unit processes that build delta distributary channels. The first is mouth-bar deposition at the shoreline and subsequent channel bifurcation, which is driven by progradation of the shoreline; the second is avulsion to a new channel, a result of aggradation of the delta topset. The former creates relatively small, branching networks such as Wax Lake Delta; the latter generates relatively few, long distributaries such as the Mississippi and Atchafalaya channels on the Mississippi Delta. The relative rate of progradation to aggradation, and hence the creation of accommodation space, emerges as a controlling parameter on channel network form. Field and experimental research has identified sea level as the dominant control on Holocene delta growth worldwide, and has empirically linked channel network changes to changes in the rate of sea level rise. Here I outline a simple modeling framework for distributary network evolution, and use this to explore large-scale changes in Holocene channel pattern that have been observed in deltas such as the Rhine-Meuse and Mississippi. Rapid early- to mid-Holocene sea level rise forced many deltas into an aggradational mode, where I hypothesize that avulsion and the generation of large-scale branches should dominate. Slowing of sea level rise in the last ～6000 yr allowed partitioning of sediment into progradation, facilitating the growth of smaller-scale distributary trees at the shorelines of some deltas, and a reduction in the number of large-scale branches. Significant antecedent topography modulates delta response; the filling of large incised valleys, for example, caused many deltas to bypass the aggradational phase. Human effects on deltas can be cast in terms of geologic controls affecting accommodation: constriction of channels forces rapid local progradation and mouth-bar bifurcation, while accelerated sea level rise increases aggradation and induces more frequent channel avulsion.     

长江中游泥沙淤积与全新世长江三角洲的发育——兼论长江中游防洪策略

全新世海侵以来,巨量泥沙在长江河口湾一带堆积,形成了自西向东雁列的六个亚三角洲河口沙坝。长江口随之从镇江、扬州向东南方向推进到长兴岛一带,加之百余年来长江中游荆江河曲的演化和簰洲湾弯曲河道的发展,长江中下游河道累计延长620 km。长江中下游河床纵比降的自动调整,促使了中游一带河道泥沙淤积达12.9 m。长江河道的不断延伸加长和长江中下冲积河床对纵比降的自动调整,是长江中游河道淤涨和洪水位上升的内因;而荆江一带人类工程活动的参与,强化了长江泥沙向下游输移,加速了长江三角洲的发展和长江河道的延伸,构成了长江中游河道淤涨和洪水位上升的正反馈机制。长江中游防洪减灾策略、长江流域防洪规划和流域国土空间规划,应统筹考虑长江口浚海清淤和长江三角洲泥沙资源化利用,或开辟新的防洪入海通道,如开辟南通-如东入海新运河暨防洪新河。     

Fluvial response to rapid high‐amplitude lake‐level changes during the Late Weichselian and early Holocene,Ain River valley,Jura, France

In the present paper the effects of rapid, high‐amplitude base‐level changes during the last glacial‐interglacial transition were studied for the Ain River in eastern France. During the Würm glacial maximum (MIS 2) rapid aggradation by deep‐water Gilbert‐type deltas and shallow‐water fan deltas occurred at the margins of a 20 to 50 m deep proglacial lake. A temporal high‐amplitude lake‐level fall of 60 m resulted in gravel deposition by forced‐regressive deltas, followed by rapid lake‐level rise and fine‐grained glaciolacustrine deposition. During the final deglaciation, a rapid base‐level fall of 40 m resulted in a complex fluvial response. Knickpoint formation and headward incision of the highstand deltas and concomitant deposition of gravel sheets by forced‐regressive deltas and braided systems occurred in several depocentres on the former glacial lake floor. Preservation of highstand and falling‐stage deposits and terrace formation in the incised valley depended on vertical incision and lateral channel migration. Terraces are well developed in the former lake‐floor depressions, whereas vertical incision was dominant in the higher lake‐floor areas. The Ain terrace staircase was likely formed by autogenic processes during a single allogenic base‐level fall. This case study possibly offers an analogue for the preservation of interglacial highstand coastal deltas during sea‐level fall at warm‐to‐cold climate transitions, although the rates of base‐level fall are different.     

Holocene transgression of the Rhine river mouth area,The Netherlands/Southern North Sea: palaeogeography and sequence stratigraphy

This study presents a detailed reconstruction of the palaeogeography of the Rhine valley (western Netherlands) during the Holocene transgression with systems tracts placed in a precise sea‐level context. This approach permits comparison of actual versus conceptual boundaries of the lowstand, transgressive and highstand systems tracts. The inland position of the highstand Rhine river mouth on a wide, low‐gradient continental shelf meant that base‐level changes were the dominant control on sedimentation for a relatively short period of the last glacial cycle. Systems in such inland positions predominantly record changes in the balance between river discharge and sediment load, and preserve excellent records of climatic changes or other catchment‐induced forcing. It is shown here that the transgressive systems tract‐part of the coastal prism formed in three stages: (i) the millennium before 8·45 ka bp , when the area was dominated by fluvial environments with extensive wetlands; (ii) the millennium after 8·45 ka, characterized by strong erosion, increasing tidal amplitudes and bay‐head delta development; and (iii) the period between 7·5 and 6·3 ka bp when the Rhine avulsed multiple times and the maximum flooding surface formed. The diachroneity of the transgressive surface is strongly suppressed because of a pulse of accelerated sea‐level rise at 8·45 ka bp . That event not only had a strong effect on preservation, but has circum‐oceanic stratigraphical relevance as it divides the early and middle Holocene parts of coastal successions worldwide. The palaeogeographical reconstruction offers a unique full spatial–temporal view on the coastal and fluvial dynamics of a major river mouth under brief rapid forced transgression. This reconstruction is of relevance for Holocene and ancient transgressive systems worldwide, and for next‐century natural coasts that are predicted to experience a 1 m sea‐level rise.     

Morphological evolution of the Guadiana estuary and intertidal zone in response to projected sea‐level rise and sediment supply scenarios

Analysis of Holocene sediment accumulation in the Guadiana estuary (southern Portugal) during sea‐level rise since ca. 13 cal. ka BP was used to simulate the long‐term morphological evolution of the lower Guadiana estuary and the associated intertidal zone for 21st‐century predicted sea‐level rises. Three sea‐level rise scenarios given by the IPCC ( 2007 ) were used in the simulations of morphology using a large‐scale behaviour‐oriented modelling approach. Sedimentation rate scenarios were derived both from the Holocene evolution of the estuary and from a semi‐empirical estimation of present‐day sediment aggradation. Our results show that the net lateral expansion of the intertidal zone area would be about 3–5% of the present intertidal zone area for each 10 cm rise in sea level. Under constraints imposed by the lack of fluvial sediment supply, the lateral expansion of the landward boundary of the intertidal zone will occur mainly in the Portuguese margin of the Guadiana estuary, while submergence of the salt marshes will occur in the Spanish margin. Therefore the Spanish margin is highly vulnerable to both sea‐level rise and lack of sediment supply. Copyright © 2011 John Wiley & Sons, Ltd.     

From outwash to coastal systems in the Portneuf–Forestville deltaic complex (Québec North Shore): Anatomy of a forced regressive deglacial sequence

Deglacial sequences typically include backstepping grounding zone wedges and prevailing glaciomarine depositional facies. However, in coastal domains, deglacial sequences are dominated by depositional systems ranging from turbiditic to fluvial facies. Such deglacial sequences are strongly impacted by glacio‐isostatic rebound, the rate and amplitude of which commonly outpaces those of post‐glacial eustatic sea‐level rise. This results in a sustained relative sea‐level fall covering the entire depositional time interval. This paper examines a Late Quaternary, forced regressive, deglacial sequence located on the North Shore of the St. Lawrence Estuary (Portneuf Peninsula, Québec, Canada) and aims to decipher the main controls that governed its stratigraphic architecture. The forced regressive deglacial sequence forms a thick (>100 m) and extensive (>100 km²) multiphased deltaic complex emplaced after the retreat of the Laurentide Ice Sheet margin from the study area ca 12 500 years ago. The sedimentary succession is composed of ice‐contact, glaciomarine, turbiditic, deltaic, fluvial and coastal depositional units. A four‐stage development is recognized: (i) an early ice‐contact stage (esker, glaciomarine mud and outwash fan); (ii) an in‐valley progradational stage (fjord head or moraine‐dammed lacustrine deltas) fed by glacigenics; (iii) an open‐coast deltaic progradation, when proglacial depositional systems expanded beyond the valley outlets and merged together; and (iv) a final stage of river entrenchment and shallow marine reworking that affected the previously emplaced deltaic complex. Most of the sedimentary volume (10 to 15 km³) was emplaced during the three‐first stages over a ca 2 kyr interval. In spite of sustained high rates of relative sea‐level fall (50 to 30 mm·year^?1), delta plain accretion occurred up to the end of the proglacial open‐coast progradational stage. River entrenchment only occurred later, after a significant decrease in the relative sea‐level fall rates (<30 mm·year^?1), and was concurrent with the formation and preservation of extensive coastal deposits (raised beaches, spit platform and barrier sands). The turnaround from delta plain accretion to river entrenchment and coastal erosion is interpreted to be a consequence of the retreat of the ice margin from the river drainage basins that led to the drastic drop of sediment supply and the abrupt decrease in progradation rates. The main internal stratigraphic discontinuity within the forced regressive deglacial sequence does not reflect changes in relative sea‐level variations.     

晚更新世江苏海岸带沉积分布模拟研究

中国边缘海大陆架在晚更新世时期曾是海岸平原,在古长江、古黄河泥沙填充下形成了陆架堆积体,并在全新世发育了南黄海辐射沙脊群、废黄河三角洲和长江三角洲。根据点状的地质钻孔分析和重建,对南黄海-江苏海岸的沉积体系的分布和变化机制尚不明了。作为动力机制探讨,基于气候-海面-沉积系统,根据气候水文学、沉积学原理以及泥沙沉积面的动态高程计算,构建了气候冰川驱动-东黄海地海系统响应-河流沉积建造的数值模式,模拟了14万年、4万年和1万年不同时间尺度江苏海岸线和长江三角洲沉积的变化过程和分布,进而对冰川气候、构造沉降、沉积压实等复杂效应下的海面特征、陆源泥沙沉积和海岸线进行分析。模拟结果与地质钻孔资料揭示的层序和埋深能够进行对比。     

Autostratigraphic modelling of the growth of alluvial‐shelf systems during steady base‐level cycles: Two‐dimensional tank experiments

Autostratigraphy is the stratigraphy generated by large‐scale autogenesis, developed based on the full recognition of the non‐equilibrium behaviour of depositional systems in response to steady external forcing. The existing autostratigraphic concepts were derived mostly from studies of river deltas growing during a single rise or fall of base level (or relative sea level). The present study challenges to extend the autostratigraphic framework to the alluvial‐shelf system growing through steady base‐level cycles by two‐dimensional tank experiments. During each experimental run, the base level was changed symmetrically, wherein through cycles, the rise and fall had the same constant rate (|R_bl|) and period (T_bl), and thus the same constant amplitude (A_bl), but with no basin tectonism. In total, nine runs with different combinations of |R_bl| and A_bl were performed. The experimental results brought the following implications. (i) The shelf‐transiting active depositional system takes non‐equilibrium responses in earlier base‐level cycles, during each of which the system experiences episodes of degradation with base‐level fall. (ii) After the system has sufficiently grown through cycles, non‐equilibrium responses change into equilibrium responses, whereby the shelf‐transiting alluvial system, whether retrogradational or progradational, is free from degradation and continues to aggrade but with a gradually decreasing rate of aggradation. (iii) The alluvial topset river tentatively but autogenically attains a graded state during the falling limb of an intermediate cycle, which separates the earlier degradation‐inclusive and later aggradation‐sustainable cycles. (iv) The number (or duration) of cycles elapsed prior to this phase‐transition is linearly proportional to the amplitude (or the square of the period) of base‐level cycles, with a coefficient defined by the rates of base‐level change and sediment supply. Such a growth pattern does not necessarily hold when considering long‐term tectonic subsidence or uplift. These notions help to understand the stratigraphic architectures of natural alluvial‐shelf systems evolved through base level cycles.     

Sea-level and tectonic control of middle to late Pleistocene turbidite systems in Santa Monica Basin, offshore California

Small turbidite systems offshore from southern California provide an opportunity to track sediment from river source through the turbidity‐current initiation process to ultimate deposition, and to evaluate the impact of changing sea level and tectonics. The Santa Monica Basin is almost a closed system for terrigenous sediment input, and is supplied principally from the Santa Clara River. The Hueneme fan is supplied directly by the river, whereas the smaller Mugu and Dume fans are nourished by southward longshore drift. This study of the Late Quaternary turbidite fill of the Santa Monica Basin uses a dense grid of high‐resolution seismic‐reflection profiles tied to new radiocarbon ages for Ocean Drilling Program (ODP) Site 1015 back to 32 ka. Over the last glacial cycle, sedimentation rates in the distal part of Santa Monica Basin averaged 2–3 mm yr^?1, with increases at times of extreme relative sea‐level lowstand. Coarser‐grained mid‐fan lobes prograded into the basin from the Hueneme, Mugu and Dume fans at times of rapid sea‐level fall. These pulses of coarse‐grained sediment resulted from river channel incision and delta cannibalization. During the extreme lowstand of the last glacial maximum, sediment delivery was concentrated on the Hueneme Fan, with mean depositional rates of up to 13 mm yr^?1 on the mid‐ and upper fan. During the marine isotope stage (MIS) 2 transgression, enhanced rates of sedimentation of > 4 mm yr^?1 occurred on the Mugu and Dume fans, as a result of distributary switching and southward littoral drift providing nourishment to these fan systems. Longer‐term sediment delivery to Santa Monica Basin was controlled by tectonics. Prior to MIS 10, the Anacapa ridge blocked the southward discharge of the Santa Clara River into the Santa Monica Basin. The pattern and distribution of turbidite sedimentation was strongly controlled by sea level through the rate of supply of coarse sediment and the style of initiation of turbidity currents. These two factors appear to have been more important than the absolute position of sea level.     

Holocene depocenter shift in the middle-lower Changjiang River basins and coastal area in response to sea level change

Authors collected 38 sedimentary boreholes and numerous seismic profiles from previous publications to delineate the Holocene sedimentation rate of six major depositional sinks in the middle-lower Changjiang River basins and its river coast. The results demonstrate that the highest sedimentation rate of ca. 15 m/ka occurred in the mono-depositional sink of the former Changjiang River mouth during 10 000–8000 aBP, when post-glacial transgression happened and the Changjiang water level remained at lower stand. With the rising of the Changjiang water level in response to sea level rise, Jianghan Basin of the middle Changjiang River becomes the other important depositional sink with highest sedimentation rate of 10 m/ka since 7000 aBP. As Jianghan Basin was mostly filled up at ca. 4000 aBP, Dongting Basin and the lower Changjiang valley trapped sediments in great amounts like in the river mouth. A considerable amount of Changjiang sediments has been delivered, both eastward and southward, to the inner continental shelf of the East China Sea, especially after 2000 aBP. This indicates reduced sediment storage capacity of the middle-lower Changjiang valley and the river mouth. In total, ca. 1307.4 billion tons of sediment have deposited in the middle-lower Changjiang floodplain since 7000 aBP. In the meantime, ca. 947 billion tons of sediment have been deposited in the river coast to form the Changjiang subaqueous delta and the Zhejiang-Fujian along-shelf mud wedge. Our result also reveals two time stages with lower sedimentation rates(< 4 m/ka) in all basins during 8000–7000 aBP and in the estuarine area during 4000–2000 aBP, probably owing to stengthened chemical weathering of decline of monsoon precipitation. __________ Translated from Journal of Palaeogeography, 2007, 9(4):419–429 [译自: 古地理学报]     

Spatial variability in fluvial style and likely responses to sea‐level change,Herbert River,Queensland

In broad terms, fluvial systems can be considered as comprising two basic geomorphologic features, a channel and its floodplain (overbank), each of which may accumulate sediment or undergo erosion. The sedimentary relationships between channels and floodplains, the resultant sedimentary architecture and the form of the dependent landscape may all be considered in terms of the relative rates of channel and floodplain aggradation and/or erosion. Using this approach, the Herbert River in north Queensland can be divided into seven ‘fluvial fields’. By considering the likely migration directions of field boundaries in the lower floodplain we conclude that, contrary to many sequence‐stratigraphic models, lowering sea‐level would drive a general aggradation of the system on the Great Barrier Reef shelf, whereas a sea‐level rise would cause further incision of the modern coastal plain.     

Sedimentological and ichnological implications of rapid Holocene flooding of a gently sloping mud‐dominated incised valley – an example from the Red River (Gulf of Tonkin)

The Gulf of Tonkin coastline migrated at an average rate of ca 60 m year^?1 landward during Holocene sea‐level rise (20 to 8 ka). Due to a combination of rapid coastline migration and undersupply of sand, neither coastal barriers nor tidal sand bars developed at the mouth of the Red River incised valley. Only a 30 to 80 cm thick sandy interval formed at the base of full‐marine deposits. Thus, the river mouth represented a mud‐dominated open funnel‐shaped estuary during transgression. At the base of the valley fill, a thin fluvial lag deposit marks a period of lowered sea‐level when the river did not reach geomorphic equilibrium and was thus prone to erosion. The onset of base‐level rise is documented by non‐bioturbated to sparsely bioturbated mud that occasionally contains pyrite indicating short‐term seawater incursions. Siderite in overlying deposits points to low‐salinity estuarine conditions. The open funnel‐shaped river mouth favoured upstream incursion of seawater that varied inversely to the seasonal strongly fluctuating discharge: several centimetres to a few tens of centimetres thick intervals showing marine or freshwater dominance alternate, as indicated by bioturbational and physical sedimentary structures, and by the presence of Fe sulphides or siderite, respectively. Recurrent short‐term seawater incursions stressed the burrowing fauna. The degree of bioturbation increases upward corresponding to increasing marine influence. The uppermost estuarine sediments are completely bioturbated. The estuarine deposits aggraded on average rapidly, up to several metres kyr^?1. Siphonichnidal burrows produced by bivalves, however, document recurrent episodes of enhanced deposition (>0·5 m) and pronounced erosion (<1 m) that are otherwise not recorded. The slope of the incised valley affected the sedimentary facies. In steep valley segments, the marine transgressive surface (equivalent to the onset of full‐marine conditions) is accentuated by the Glossifungites ichnofacies, whereas in gently sloped valley segments the marine transgressive surface is gradational and bioturbated. Marine deposits are completely bioturbated.     

Sharp-based, tide-dominated deltas of the Sego Sandstone, Book Cliffs, Utah, USA

The lower part of the Cretaceous Sego Sandstone Member of the Mancos Shale in east‐central Utah contains three 10‐ to 20‐m thick layers of tide‐deposited sandstone arranged in a forward‐ and then backward‐stepping stacking pattern. Each layer of tidal sandstone formed during an episode of shoreline regression and transgression, and offshore wave‐influenced marine deposits separating these layers formed after subsequent shoreline transgression and marine ravinement. Detailed facies architecture studies of these deposits suggest sandstone layers formed on broad tide‐influenced river deltas during a time of fluctuating relative sea‐level. Shale‐dominated offshore marine deposits gradually shoal and become more sandstone‐rich upward to the base of a tidal sandstone layer. The tidal sandstones have sharp erosional bases that formed as falling relative sea‐level allowed tides to scour offshore marine deposits. The tidal sandstones were deposited as ebb migrating tidal bars aggraded on delta fronts. Most delta top deposits were stripped during transgression. Where the distal edge of a deltaic sandstone is exposed, a sharp‐based stack of tidal bar deposits successively fines upward recording a landward shift in deposition after maximum lowstand. Where more proximal parts of a deltaic‐sandstone are exposed, a sharp‐based upward‐coarsening succession of late highstand tidal bar deposits is locally cut by fluvial valleys, or tide‐eroded estuaries, formed during relative sea‐level lowstand or early stages of a subsequent transgression. Estuary fills are highly variable, reflecting local depositional processes and variable rates of sediment supply along the coastline. Lateral juxtaposition of regressive deltaic deposits and incised transgressive estuarine fills produced marked facies changes in sandstone layers along strike. Estuarine fills cut into the forward‐stepped deltaic sandstone tend to be more deeply incised and richer in sandstone than those cut into the backward‐stepped deltaic sandstone. Tidal currents strongly influenced deposition during both forced regression and subsequent transgression of shorelines. This contrasts with sandstones in similar basinal settings elsewhere, which have been interpreted as tidally influenced only in transgressive parts of depositional successions.     

Autocyclic behaviour of fan deltas: an analogue experimental study

Fan deltas are excellent recorders of fan‐building processes because of their high sedimentation rate, particularly in tectonically active settings. Although previous research focuses mainly on allogenic controls, there is clear evidence for autogenically produced storage and release of sediment by flume and numerical modelling that demands further definition of characteristics and significance of autogenically forced facies and stratigraphy. Analogue experiments were performed on fan deltas with constant extrinsic variables (discharge, sediment supply, sea‐level and basin relief) to demonstrate that fan‐delta evolution consists of prominent cyclic alternations of channellized flow and sheet flow. The channellized flow is initiated by slope‐induced scouring and subsequent headward erosion to form a channel that connected with the valley, while the removed sediment is deposited in a rapidly prograding delta lobe. The resulting decrease in channel gradient causes a reduction in flow strength, mouth‐bar formation, flow bifurcation and progressive backfilling of the channel. In the final stage of channel filling, sheet flow coexists for a while with channellized flow (semi‐confined flow), although in cycle 1 this phase of semi‐confined flow was absent. Subsequent autocyclic incisions are very similar in morphology and gradient. However, they erode deeper into the delta plain and, as a result, take more time to backfill. The duration of the semi‐confined flow increases with each subsequent cycle. During the period of sheet flow, the delta plain aggrades up to the ‘critical’ gradient required for the initiation of autocyclic incision. This critical gradient is dependent on the sediment transport capacity, defined by the input conditions. These autogenic cycles of erosion and aggradation confirm earlier findings that storage and release of sediment and associated slope variation play an important role in fan‐delta evolution. The erosional surfaces produced by the autocyclic incisions are well‐preserved by the backfilling process in the deposits of the fan deltas. These erosional surfaces can easily be misinterpreted as climate, sea‐level or tectonically produced bounding surfaces.