The role of climate variation in delta architecture: lessons from analogue modelling

Sequence‐stratigraphic models for fourth to sixth order, glacio‐eustatic sequences based only on relative sea‐level variations result in simplified and potentially false interpretations. Glacio‐eustatic sea‐level variations form only one aspect of cyclic climate variation; other aspects, such as variations in fluvial water discharge, vegetation cover, weathering and sediment supply can lead to variable sediment yield, thus adding complexity to sequence‐stratigraphic patterns normally attributed to sea‐level variations. Analogue flume models show a significant impact of water discharge on the timing and character of sequence boundaries, and on changes in the relative importance of systems tracts, as expressed in sediment volumes. Four deltas, generated under the influence of an identical sea‐level curve, and affected by different water‐discharge cycles were generated in the Eurotank facility: (1) constant discharge; (2) high‐frequency discharge variations (HFD); (3) discharge leading sea level by a quarter phase; (4) discharge lagging sea level by a quarter phase. HFD shift the parasequence stacking pattern consistently but do not alter large‐scale delta architecture. Water‐discharge changes that lead sea‐level changes result in high sediment yield during sea‐level rise and in the poor development of maximum flooding surfaces. Delta‐front erosion during sea‐level fall is expressed by multiple, small channels related to upstream avulsions, and does not result in an incised valley that efficiently routs sediment to the shelf edge. When water‐discharge changes lag sea‐level changes, sediment yield is high during falling sea level and results in rapid progradation during forced regression. Erosion from incised valleys is strong on the proximal delta top and dissipates towards the delta front. The combination of high discharge and sea‐level fall provides the most efficient mode of valley incision and sediment transport to the shelf edge. During sea‐level rise, low water discharge results in sediment starvation and well‐developed maximum flooding surfaces. Water‐discharge variations thus alter sequence‐stratigraphic patterns and provide an alternative explanation to the amplitude of sea‐level fall for generating either type 1 or 2 erosional unconformities.     

Relationships between shelf‐edge trajectories and sediment dispersal along depositional dip and strike: a different approach to sequence stratigraphy

Analysis of shelf‐edge trajectories in prograding successions from offshore Norway, Brazil, Venezuela and West Africa reveals systematic changes in facies associations along the depositional dip. These changes occur in conjunction with the relative sea‐level change, sediment supply, inclination of the substratum and the relief of the margin. Flat and ascending trajectories generally result in an accumulation of fluvial and shallow marine sediments in the topset segment. Descending trajectories will generally result in erosion and bypass of the topset segment and deposition of basin floor fans. An investigation of incised valley fills reveals multiple stages of filling that can be linked to distinct phases of deepwater fan deposition and to the overall evolution of the margin. In the case of high sediment supply, like the Neogene Niger and Orinoco deltas, basin floor fans may develop systematically even under ascending trajectory styles. In traditional sequence stratigraphic thinking, this would imply the deposition of basin floor fans during a period of relative sea‐level highstand. Facies associations and sequence development also vary along the depositional strike. The width and gradient of the shelf and slope show considerable variations from south to north along the Brazilian continental margin during the Cenozoic. During the same time interval, the continental shelf may display high or low accommodation conditions, and the resulting stacking patterns and facies associations may be utilized to reconstruct palaeogeography and for prediction of lithology. Application of the trajectory concept thus reveals nuances in the rock record that would be lost by the application of traditional sequence stratigraphic work procedures. At the same time, the methodology simplifies the interpretation in that less importance is placed on interpretation and labelling of surface boundaries and systems tracts.     

Quantitative analogue flume-model study of river–shelf systems: principles and verification exemplified by the Late Quaternary Colorado river–delta evolution

ABSTRACT Physical modelling of clastic sedimentary systems over geological time spans has to resort to analogue modelling since full scaling cannot be achieved within the spatial and temporal restrictions that are imposed by a laboratory set‐up. Such analogue models are suitable for systematic investigation of a sedimentary system's sensitivity to allocyclic changes by isolating governing parameters. Until now, analogue models of landscape evolution were mainly qualitative in nature. In this paper, we present a quantitative approach. The quantitative experimental results are verified and discussed by comparison with high‐resolution data from the Colorado river–shelf system of the Texas shelf that we used as a prototype. The model's dimensions are proportionally scaled to the prototype, except for a vertical exaggeration. Time is scaled using a Basin Response factor to maintain a similar ratio between the period of change and the system's equilibrium time for model and prototype. A Basin Fill factor was used to compare the ratio between the time‐averaged sedimentation rate and the rate of change in accommodation space of model and prototype. The flume‐model results are in the form of sediment budgets that are related to shelf cannibalism and fluvial supply, which are compared with the ancestral Colorado river–delta evolution of the last 40 kyr. Model and prototype have similarities in delta evolution in response to one cycle of sea‐level change. With sea‐level change as the isolated variable, the flume model generates a significant supply pulse caused by headward erosion of the shelf in response to the sea‐level fall. This pulse adds to the yield of the hinterland. The supply induced by sea‐level change persists during the early rise, although its rate declines. A similar trend is observed on the east Texas shelf. We argue that shelfal and fluvial degradation cycles induced by sea‐level changes can significantly influence the timing and amount of sediment supply to basins and must therefore be taken into consideration.     

Deciphering the origin of cyclical gravel front and shoreline progradation and retrogradation in the stratigraphic record

Nearly all successions of the near‐shore strata exhibit cyclical movements of the shoreline, which have commonly been attributed to cyclical oscillations in relative sea level (combining eustasy and subsidence) or, more rarely, to cyclical variations in sediment supply. It has become accepted that cyclical change in sediment delivery from source catchments may lead to cyclical movement of boundaries such as the gravel front, particularly in the proximal segments of sediment‐routing systems. In order to quantitatively assess how variations in sediment transport as a consequence of change in relative sea‐level and surface run‐off control stratigraphic architecture, we develop a simple numerical model of sediment transport and explore the sensitivity of moving boundaries within the sediment‐routing system to change in upstream (sediment flux, precipitation rate) and downstream (sea level) controls. We find that downstream controls impact the shoreline and sand front, while the upstream controls can impact the whole system depending on the amplitude of change in sediment flux and precipitation rate. The model implies that under certain conditions, the relative movement of the gravel front and shoreline is a diagnostic marker of whether the sediment‐routing system experienced oscillations in sea level or climatic conditions. The model is then used to assess the controls on stratigraphic architecture in a well‐documented palaeo‐sediment‐routing system in the Late Cretaceous Western Interior Seaway of North America. Model results suggest that significant movement of the gravel front is forced by pronounced (±50%) oscillations in precipitation rate. The absence of such movement in gravel front position in the studied strata implies that time‐equivalent movement of the shoreline was driven by relative sea‐level change. We suggest that tracking the relative trajectories of internal boundaries such as the gravel front and shoreline is a powerful tool in constraining the interpretation of stratigraphic sequences.     

Shelf-margin clinoforms and prediction of deepwater sands

Early Eocene successions from Spitsbergen and offshore Ireland, showing well‐developed shelf‐margin clinoforms and a variety of deepwater sands, are used to develop models to predict the presence or absence of turbidite sands in clinoform strata without significant slope disturbance/ponding by salt or mud diapers. The studied clinoforms formed in front of narrow to moderate width (10–60 km) shelves and have slopes, 2–4°, that are typical of accreting shelf margins. The clinoforms are evaluated in terms of both shelf‐transiting sediment‐delivery systems and the resultant partitioning of the sand and mud budget along their different segments. Although this sediment‐budget partitioning is controlled by sediment type and flux, shelf width and gradient, process regime on the shelf and relative sea‐level behaviour, the most tell‐tale or predictive signs in the stratigraphic record appear to be (1) sediment‐delivery system type, (2) degree of shelf‐edge channelling and (3) character of shelf‐edge trajectory through time. The clinoform data sets from the Porcupine Basin (wells and 3‐D seismic) and from the Central Basin on Spitsbergen (outcrops) suggest that river‐dominated deltas are the most efficient delivery systems for dispersing sand into deep water beyond the shelf‐slope break. In addition, low‐angle or flat, channelled shelf‐edge trajectories associate with co‐eval deepwater slope and basin‐floor sands, whereas rising trajectories tend to associate with muddy slopes and basin floors. Characteristic features of the shelf‐edge, slope and basin‐floor segments of clinoforms for these trajectory types are documented. Seismic lines along the slope to basin‐floor transects tend to show apparent up‐dip sandstone pinchouts, but most of these are likely to be simply sidelap features. Dip lines aligned along the axes of sandy fairways show that stratigraphic traps are unlikely, unless slope channels become mud‐filled or are structurally partitioned. Another feature that is prominent in the data sets examined is the lack of slope onlap. During the relative rise of sea level back up to the shelf, the clinoform slopes are generally mud‐prone and they are characteristically aggradational.     

Allogenic forcing of the late Quaternary Rhine–Meuse fluvial record: the interplay of sea-level change, climate change and crustal movements

The Rhine–Meuse system in the west‐central Netherlands is a continental‐scale fluvial system bordered by an extremely wide continental shelf. Consequently, late Quaternary eustatic sea‐level changes have resulted in dramatic shoreline displacements, by as much as 800 km. In addition, changes in climate have been severe, given the latitudinal and palaeogeographic setting of the Rhine–Meuse system. We investigated the relative importance of these allogenic controls on fluvial aggradation and incision during the last two glacial–interglacial cycles. We used optical dating of quartz from ～30 samples in a cross‐section perpendicular to the palaeoflow direction, allowing us to correlate periods of aggradation and incision with independent records of sea‐level change, climate change and glacio‐isostatic crustal movements. We found the long‐term aggradation rate to be ～8 cm kyr^?1, a value similar to previous estimates of tectonic subsidence rates in the study area. Several excursions from this long‐term aggradation trend could be identified for the last glacial–interglacial cycle. Dry climatic conditions with relatively high sediment supply induced aggradation during oxygen‐isotope stages (OIS) 4 and 3. Build‐up of a glacio‐isostatic forebulge during OIS 2 is a likely cause of incision around the Last Glacial Maximum, followed by an aggradation phase during forebulge collapse. Sea‐level highstands during OIS 5 have likely resulted in the aggradation of coastal prisms, but only minor, basal estuarine deposits have been preserved because these coastal prisms were prone to erosion during ensuing sea‐level falls. Overall, the sedimentary record is dominated by strata formed during time intervals when the study area was completely unaffected by sea‐level control, and our evidence shows that the falling‐stage systems tract has the highest preservation potential. Our study highlights the importance of considering the complex interplay of both upstream and downstream controls to obtain a comprehensive understanding of the evolution of basin‐margin successions.     

Low‐accommodation and backwater effects on sequence stratigraphic surfaces and depositional architecture of fluvio‐deltaic settings (Cretaceous Mesa Rica Sandstone,Dakota Group,USA)

The adequate documentation and interpretation of regional‐scale stratigraphic surfaces is paramount to establish correlations between continental and shallow marine strata. However, this is often challenged by the amalgamated nature of low‐accommodation settings and control of backwater hydraulics on fluvio‐deltaic stratigraphy. Exhumed examples of full‐transect depositional profiles across river‐to‐delta systems are key to improve our understanding about interacting controlling factors and resultant stratigraphy. This study utilizes the ~400 km transect of the Cenomanian Mesa Rica Sandstone (Dakota Group, USA), which allows mapping of down‐dip changes in facies, thickness distribution, fluvial architecture and spatial extent of stratigraphic surfaces. The two sandstone units of the Mesa Rica Sandstone represent contemporaneous fluvio‐deltaic deposition in the Tucumcari sub‐basin (Western Interior Basin) during two regressive phases. Multivalley deposits pass down‐dip into single‐story channel sandstones and eventually into contemporaneous distributary channels and delta‐front strata. Down‐dip changes reflect accommodation decrease towards the paleoshoreline at the Tucumcari basin rim, and subsequent expansion into the basin. Additionally, multi‐storey channel deposits bound by erosional composite scours incise into underlying deltaic deposits. These represent incised‐valley fill deposits, based on their regional occurrence, estimated channel tops below the surrounding topographic surface and coeval downstepping delta‐front geometries. This opposes criteria offered to differentiate incised valleys from flood‐induced backwater scours. As the incised valleys evidence relative sea‐level fall and flood‐induced backwater scours do not, the interpretation of incised valleys impacts sequence stratigraphic interpretations. The erosional composite surface below fluvial strata in the continental realm represents a sequence boundary/regional composite scour (RCS). The RCS’ diachronous nature demonstrates that its down‐dip equivalent disperses into several surfaces in the marine part of the depositional system, which challenges the idea of a single, correlatable surface. Formation of a regional composite scour in the fluvial realm throughout a relative sea‐level cycle highlights that erosion and deposition occur virtually contemporaneously at any point along the depositional profile. This contradicts stratigraphic models that interpret low‐accommodation settings to dominantly promote bypass, especially during forced regressions. Source‐to‐sink analyses should account for this in order to adequately resolve timing and volume of sediment storage in the system throughout a complete relative sea‐level cycle.     

Middle Miocene–Pliocene siliciclastic influx across a carbonate shelf and influence of deltaic sedimentation on shelf construction,Northern Carnarvon Basin,Northwest Shelf of Australia

Middle Miocene to Pliocene siliciclastics of the Bare Formation represent a long‐lived (ca. 11 Myr) break in the otherwise carbonate‐dominated shelf of the Northern Carnarvon Basin, Northwest Shelf of Australia. The quartz‐sandstone interval is correlated with the appearance of spectacular clinoform sets mapped on 3D and dense 2D seismic data. Twenty‐seven clinoform sets are interpreted as delta lobes primarily based on their plan‐view morphology (strike‐elongate to lobate features) and their 40–100‐m‐high clinoform amplitudes. The delta lobes were deposited on outer‐shelf to shelf‐edge positions, and the older deltas show evidence of a higher degree wave reworking than the younger deltas. Measurements of the along‐strike (migration) and down‐dip (progradation) movement of these deltas are compared with relative sea‐level behaviour inferred from shelf‐edge trajectory analysis. Delta lobes exhibit greater lateral shifting during relative sea‐level rise, whereas delta lobes are more restricted to dip‐oriented fairways during sea‐level fall, although no major incised valleys have been identified. Long‐term (cumulative) progradation of this delta system and subsequent backstepping correlates with long‐term sea‐level fall and rise during the late middle and late Miocene. In addition, a long‐term northeastward migration trend for these delta lobes was likely a result of localized uplift of an inversion anticline in the Rosemary–Legendre Trend; the growth of this anticline probably steered the fluvial source for the delta system towards the northeast. The Bare Formation siliciclastic influx correlates with other middle Miocene increases in siliciclastic sediment supply worldwide. Global cooling and a shift to more arid conditions, negatively influencing vegetation cover, may have combined with more seasonally variable rainfall to generate the high sediment supply that built the deltas. Retreat of the siliciclastics could correlate with ice‐sheet growth in the Northern Hemisphere and/or increase in the Indonesian Throughflow and Leeuwin Current (ca. 1.6 Ma), which might have modified climate regionally.     

Role of autoretreat and A/S changes in the understanding of deltaic shoreline trajectory: a semi-quantitative approach

ABSTRACT There is continued interest in how the rate of relative sea‐level rise [A ( > 0)] and the rate of sediment supply [S] function during the growth and evolution of deltaic shorelines. The theory of shoreline autoretreat, recently corroborated in flume experiments, claims that (1) A( > 0) and S can never be in equilibrium, and (2) shoreline or shelf‐edge progradation inevitably turns to retrogradation, when relative sea level is rising even modestly and even if A/S = const (> 0). Autoretreat arises because the area of the clinoform surface of the delta (or shelf edge) per kilometer of shoreline must increase as the relative sea level rises, and the delta (or shelf edge) progrades into deeper water. A finite sediment supply rate is thus liable to become inadequate to sustain progradation. The problem increases further as a rising sea level also greatly increases the delta‐plain volume that needs to be filled, further limiting the progradation of the system. The fundamental trajectory of shoreline migration is thus one characterized by a concave‐landward shape, even under the steady forcing of the basin. The magnitudes of A (> 0) and S, or A/S do not determine whether the landward turnaround of the shoreline is realized or not, but affect merely the length and height of the fundamental trajectory curve. Thus, any attempt to detect and interpret temporal changes in A and S from the observed stratigraphic record of shoreline trajectory needs first to take full account of the inbuilt autoretreat mechanism. We develop here a simple, semi‐quantitative method of reconstructing the basin conditions (A and S) from the stratigraphic record of prograding deltaic shorelines (or prograding shelf‐margin clinoforms) on the basis of the theory of shoreline autoretreat. The deterministic nature of the autoretreat theory is advantageous in managing this latter issue, because any expected or unexpected change emerges as some discrepancy from a trajectory that was predicted for the initial conditions. The autoretreat theory also provides a convenient graphical method of dealing with the uncertainty of the field data, and with evaluating the accuracy of any reconstruction. Our methodology has been developed to deal with the behaviour of deltaic shorelines, but is basically applicable to any clinoform system, the development of which is affected by relative sea level. The suggested method is applied to an Early Eocene (Ypresian) regressive shoreline succession in the Central Tertiary Basin on Spitsbergen. The studied regressive wedge developed as a delta‐driven, progradational shelf‐margin system under a regime of overall (i.e. long‐term) rise of relative sea level, but also suffered short‐term sea‐level falls associated with valley incisions on the coastal plain and shelf. On the assumption that S was constant or was steadily decreasing, the analysis of field data obtained from three sites within the basin suggests that the initial water depth in the basin was around 0.45 km, and that the overall relative sea‐level rise (c. 0.80 km) happened largely during an early time period and was followed by a longer period of much lower rate of rise. This pattern of relative sea‐level rise is consistent with the Palaeogene tectonic subsidence trend of the basin which was determined independently through a geohistory analysis. The uncertainty of the field data does not negate our reconstruction. The combined effects of autoretreat and A/S changes on a deltaic shoreline trajectory are confirmed through the development of an autoretreat‐based methodology. Conventional sequence stratigraphic models that assume a possible equilibrium condition between A and S are both conceptually misleading and insufficient to analyse basin conditions quantitatively. Sequence stratigraphic analyses of shorelines need to incorporate the autoretreat concept.     

Allogenic controls on autogenic variability in fluvio‐deltaic systems: inferences from analysis of synthetic stratigraphy

Fluvio‐deltaic stratigraphy develops under continuous morphodynamic interactions of allogenic and autogenic processes, but the role and relative contribution of these processes to the stratigraphic record are poorly understood. We analysed synthetic fluvio‐deltaic deposits of several accommodation‐to‐supply cycles (sequences) with the aim to relate stratigraphic variability to autogenic and allogenic controls. The synthetic stratigraphy was produced in a series of long time‐scale (10⁵ years) numerical experiments with an aggregated process‐based model using a typical passive‐margin topography with constant rates of liquid and solid river discharge subjected to sinusoidal sea‐level fluctuation. Post‐processing of synthetic stratigraphy allowed us to quantify stratigraphic variability by means of local and regional net sediment accumulation over equally spaced time intervals (1–10 kyr). The regional signal was subjected to different methods of time‐series analysis. In addition, major avulsion sites (>5 km from the coastline) were extracted from the synthetic stratigraphy to confirm the interpretations of our analyses. Regional stratigraphic variability as defined in this study is modulated by a long‐term allogenic signal, which reflects the rate of sea‐level fluctuation, and it preserves two autogenic frequency bands: the intermediate and high‐frequency components. The intermediate autogenic component corresponds to major avulsions with a median inter‐avulsion period of ca. 3 kyr. This component peaks during time intervals in which aggradation occurs on the delta plain, because super‐elevation of channel belts is a prerequisite for large‐scale avulsions. Major avulsions occur occasionally during early stages of relative sea‐level fall, but they are fully absent once the coast line reaches the shelf edge and incision takes place. These results are consistent with a number of field studies of falling‐stage deposition in fluvial systems. The high‐frequency autogenic component (decadal to centennial time scales) represents mouthbar‐induced bifurcations occurring at the terminal parts of the system, and to a lesser extent, partial or small‐scale avulsions (<5 km from the coastline). Bifurcation intensity correlates strongly with the rate of progradation, and thus reaches its maximum during forced regression. However, its contribution to overall stratigraphic variability is much less than that of the large‐scale avulsions, which affect the entire area downstream of avulsion nodes. The results of this study provide guidelines for predicting fluvio‐deltaic stratigraphy in the context of co‐existing autogenic and allogenic processes and underscore the fact that the relative importance and the type of autogenic processes occurring in fluvio‐deltaic systems are governed by allogenic forcing.     

100 Myr record of sequences, sedimentary facies and sea level change from Ocean Drilling Program onshore coreholes, US Mid-Atlantic coastal plain

We analyzed the latest Early Cretaceous to Miocene sections (～110–7 Ma) in 11 New Jersey and Delaware onshore coreholes (Ocean Drilling Program Legs 150X and 174AX). Fifteen to seventeen Late Cretaceous and 39–40 Cenozoic sequence boundaries were identified on the basis of physical and temporal breaks. Within‐sequence changes follow predictable patterns with thin transgressive and thick regressive highstand systems tracts. The few lowstands encountered provide critical constraints on the range of sea‐level fall. We estimated paleowater depths by integrating lithofacies and biofacies analyses and determined ages using integrated biostratigraphy and strontium isotopic stratigraphy. These datasets were backstripped to provide a sea‐level estimate for the past ～100 Myr. Large river systems affected New Jersey during the Cretaceous and latest Oligocene–Miocene. Facies evolved through eight depositional phases controlled by changes in accommodation, long‐term sea level, and sediment supply: (1) the Barremian–earliest Cenomanian consisted of anastomosing riverine environments associated with warm climates, high sediment supply, and high accommodation; (2) the Cenomanian–early Turonian was dominated by marine sediments with minor deltaic influence associated with long‐term (10⁷ year) sea‐level rise; (3) the late Turonian through Coniacian was dominated by alluvial and delta plain systems associated with long‐term sea‐level fall; (4) the Santonian–Campanian consisted of marine deposition under the influence of a wave‐dominated delta associated with a long‐term sea‐level rise and increased sediment supply; (5) Maastrichtian–Eocene deposition consisted primarily of starved siliciclastic, carbonate ramp shelf environments associated with very high long‐term sea level and low sediment supply; (6) the late Eocene–Oligocene was a starved siliciclastic shelf associated with moderately high sea‐level and low sediment supply; (7) late early–middle Miocene consisted of a prograding shelf under a strong wave‐dominated deltaic influence associated with major increase in sediment supply and accommodation due to local sediment loading; and (8) over the past 10 Myr, low accommodation and eroded coastal systems were associated with low long‐term sea level and low rates of sediment supply due to bypassing.     

Effects of large sea‐level variations in connected basins: the Dacian–Black Sea system of the Eastern Paratethys

Sea‐level changes provide an important control on the interplay between accommodation space and sediment supply, in particular, for shallow‐water basins where the available space is limited. Sediment exchange between connected basins separated by a subaqueous sill (bathymetric threshold) is still not well understood. When sea‐level falls below the bathymetric level of this separating sill, the shallow‐water basin evolution is controlled by its erosion and rapid fill. Once this marginal basin is filled, the sedimentary depocenter shifts to the open marine basin (outward shift). With new accommodation space created during the subsequent sea‐level rise, sediment depocenter shifts backwards to the marginal basin (inward shift). This new conceptual model is tested here in the context of Late Miocene to Quaternary evolution of the open connection between Dacian and Black Sea basins. By the means of seismic sequence stratigraphic analysis of the Miocene‐Pliocene evolution of this Eastern Paratethys domain, this case study demonstrates these shifts in sedimentary depocenter between basins. An outward shift occurs with a delay that corresponds to the time required to fill the remaining accommodation space in the Dacian Basin below the sill that separates it from the Black Sea. This study provides novel insight on the amplitude and sedimentary geometry of the Messinian Salinity Crisis (MSC) event in the Black Sea. A large (1.3–1.7 km) sea‐level drop is demonstrated by quantifying coeval sedimentation patterns that change to mass‐flows and turbiditic deposits in the deep‐sea part of this main sink. The post‐MSC sediment routing continued into the present‐day pattern of Black Sea rivers discharge.     

Interactions between onshore bedrock-channel incision and nearshore wave-base erosion forced by eustasy and tectonics

We explore the response of bedrock streams to eustatic and tectonically induced fluctuations in base level. A numerical model coupling onshore fluvial erosion with offshore wave‐base erosion is developed. The results of a series of simulations for simple transgressions with constant rate of sea‐level change (SLR) show that response depends on the relative rates of rock uplift (U) and wave‐base erosion (?_w). Simple regression runs highlight the importance of nearshore bathymetry. Shoreline position during sea‐level fall is set by the relative rate of base‐level fall (U‐SLR) and ?_w, and is constant horizontally when these two quantities are equal. The results of models forced by a realistic Late Quaternary sea‐level curve are presented. These runs show that a stable shoreline position cannot be obtained if offshore uplift rates exceed ?_w. Only in the presence of a relatively stable shoreline position, fluvial profiles can begin to approximate a steady‐state condition, with U balanced by fluvial erosion rate (?_f). In the presence of a rapid offshore decrease in rock‐uplift rate (U), short (～5 km) fluvial channels respond to significant changes in rock‐uplift rate in just a few eustatic cycles. The results of the model are compared to real stream‐profile data from the Mendocino triple junction region of northern California. The late Holocene sea‐level stillstand response exhibited by the simulated channels is similar to the low‐gradient mouths seen in the California streams.     

Evolution of the late Cenozoic Chaco foreland basin, Southern Bolivia

Eastward Andean orogenic growth since the late Oligocene led to variable crustal loading, flexural subsidence and foreland basin sedimentation in the Chaco basin. To understand the interaction between Andean tectonics and contemporaneous foreland development, we analyse stratigraphic, sedimentologic and seismic data from the Subandean Belt and the Chaco Basin. The structural features provide a mechanism for transferring zones of deposition, subsidence and uplift. These can be reconstructed based on regional distribution of clastic sequences. Isopach maps, combined with sedimentary architecture analysis, establish systematic thickness variations, facies changes and depositional styles. The foreland basin consists of five stratigraphic successions controlled by Andean orogenic episodes and climate: (1) the foreland basin sequence commences between ～27 and 14 Ma with the regionally unconformable, thin, easterly sourced fluvial Petaca strata. It represents a significant time interval of low sediment accumulation in a forebulge‐backbulge depocentre. (2) The overlying ～14–7 Ma‐old Yecua Formation, deposited in marine, fluvial and lacustrine settings, represents increased subsidence rates from thrust‐belt loading outpacing sedimentation rates. It marks the onset of active deformation and the underfilled stage of the foreland basin in a distal foredeep. (3) The overlying ～7–6 Ma‐old, westerly sourced Tariquia Formation indicates a relatively high accommodation and sediment supply concomitant with the onset of deposition of Andean‐derived sediment in the medial‐foredeep depocentre on a distal fluvial megafan. Progradation of syntectonic, wedge‐shaped, westerly sourced, thickening‐ and coarsening‐upward clastics of the (4) ～6–2.1 Ma‐old Guandacay and (5) ～2.1 Ma‐to‐Recent Emborozú Formations represent the propagation of the deformation front in the present Subandean Zone, thereby indicating selective trapping of coarse sediments in the proximal foredeep and wedge‐top depocentres, respectively. Overall, the late Cenozoic stratigraphic intervals record the easterly propagation of the deformation front and foreland depocentre in response to loading and flexure by the growing Intra‐ and Subandean fold‐and‐thrust belt.     

Shallow-marine sequences as the building blocks of stratigraphy: insights from numerical modelling

Two types of depositional sequences can be defined within the sequence stratigraphic framework: the parasequence and the high‐frequency sequence. Both sequences consist of stacked regressive and transgressive deposits. However, a parasequence forms under conditions of overall sea‐level rise, whereas a high‐frequency sequence forms as the sea level oscillates which results in typical forced regressive deposits during sea‐level fall. Both depositional sequences may develop over comparable temporal (10–100 kyr) and spatial (1–20 km wide and 1–40 m thick) scales. Numerical modelling is used to compare the architecture, preservation potential, internal volumes, bounding surfaces, condensed and expanded sections and facies assemblages of parasequences and high‐frequency sequences. Deposits originating from transgression are less pronounced than their regressive counterparts and consist of either preserved backbarrier deposits or shelf deposits. Shoreface deposits are not preserved during transgression. The second half of the paper evaluates in detail the preservation potential of backbarrier deposits and proposes a mechanism that explains the occurrence of both continuous and discontinuous barrier retreat in terms of varying rates of sea‐level rise and sediment supply. The key to this mechanism is the maximum washover capacity, which plays a part in both barrier shoreline retreat and backbarrier‐lagoonal shoreline retreat. If these two shorelines are not balanced, then the retreat of the coastal system as a whole is discontinuous and in time barrier overstep may take place.     

Large‐scale connectivity of fluvio‐deltaic stratigraphy: Inferences from simulated accommodation‐to‐supply cycles and automated extraction of chronosomes

Multiscale simulation of fluvio‐deltaic stratigraphy was used to quantify the elements of the geometry and architectural arrangement of sub‐seismic‐scale fluvial‐to‐shelf sedimentary segments. We conducted numerical experiments of fluvio‐deltaic system evolution by simulating the accommodation‐to‐sediment‐supply (A/S) cycles of varying wavelength and amplitude with the objective to produce synthetic 3‐D stratigraphic records. Post‐processing routines were developed in order to investigate delta lobe architecture in relation to channel‐network evolution throughout A/S cycles, estimate net sediment accumulation rates in 3‐D space, and extract chronostratigraphically constrained lithosomes (or chronosomes) to quantify large‐scale connectivity, that is, the spatial distribution of high net‐to‐gross lithologies. Chronosomes formed under the conditions of channel‐belt aggradation are separated by laterally continuous abandonment surfaces associated with major avulsions and delta‐lobe switches. Chronosomes corresponding to periods in which sea level drops below the inherited shelf break, that is, the youngest portions of the late falling stage systems tract (FSST), form in the virtual absence of major avulsions, owing to the incision in their upstream parts, and thus display purely degradational architecture. Detailed investigation of chronosomes within the late FSST showed that their spatial continuity may be disrupted by higher‐frequency A/S cycles to produce “stranded” sand‐rich bodies encased in shales. Chronosomes formed during early and late falling stage (FSST) demonstrate the highest large‐scale connectivity in their proximal and distal areas, respectively. Lower‐amplitude base level changes, representative of greenhouse periods during which the shelf break is not exposed, increase the magnitude of delta‐lobe switching and favour the development of system‐wide abandonment surfaces, whose expression in real‐world stratigraphy is likely to reflect the intertwined effects of high‐frequency allogenic forcing and differential subsidence.     

Sink‐to‐source: Influence of offshore dynamics on upstream processes

When we model fluvial sedimentation and the resultant alluvial stratigraphy, we typically focus on the effects of local parameters (e.g., sediment flux, water discharge, grain size) and the effects of regional changes in boundary conditions applied in the source region (i.e., climate, tectonics) and at the shoreline (i.e., sea level). In recent years this viewpoint has been codified into the “source‐to‐sink” paradigm, wherein major shifts in sediment flux, grain‐size fining trends, channel‐stacking patterns, floodplain deposition and larger stratigraphic systems tracts are interpreted in terms of (1) tectonic and climatic signals originating in the hinterland that propagate downstream; and (2) eustatic fluctuation, which affects the position of the shoreline and dictates the generation of accommodation. Within this paradigm, eustasy represents the sole means by which downstream processes may affect terrestrial depositional systems. Here, we detail three experimental cases in which coastal rivers are strongly influenced by offshore and slope transport systems via the clinoform geometries typical of prograding sedimentary bodies. These examples illustrate an underdeveloped, but potentially important “sink‐to‐source” influence on the evolution of fluvial‐deltaic systems. The experiments illustrate the effects of (1) submarine hyperpycnal flows, (2) submarine delta front failure events, and (3) deformable substrates within prodelta and offshore settings. These submarine processes generate (1) erosional knickpoints in coastal rivers, (2) increased river channel occupancy times, (3) rapid rates of shoreline movement, and (4) localized zones of significant offshore sediment accumulation. Ramifications for coastal plain and deltaic stratigraphic patterns include changes in the hierarchy of scour surfaces, fluvial sand‐body geometries, reconstruction of sea‐level variability and large‐scale stratal geometries, all of which are linked to the identification and interpretation of sequences and systems tracts.     

Controls on fluvial sedimentary architecture and sediment‐fill state in salt‐walled mini‐basins: Triassic Moenkopi Formation,Salt Anticline Region,SE Utah,USA

The Triassic Moenkopi Formation in the Salt Anticline Region, SE Utah, represents the preserved record of a low‐relief ephemeral fluvial system that accumulated in a series of actively subsiding salt‐walled mini‐basins. Development and evolution of the fluvial system and its resultant preserved architecture was controlled by the following: (1) the inherited state of the basin geometry at the time of commencement of sedimentation; (2) the rate of sediment delivery to the developing basins; (3) the orientation of fluvial pathways relative to the salt walls that bounded the basins; (4) spatially and temporally variable rates and styles of mini‐basin subsidence and associated salt‐wall uplift; and (5) temporal changes in regional climate. Detailed outcrop‐based tectono‐stratigraphic analyses demonstrate how three coevally developing mini‐basins and their intervening salt walls evolved in response to progressive sediment loading of a succession of Pennsylvanian salt (the Paradox Formation) by the younger Moenkopi Formation, deposits of which record a dryland fluvial system in which flow was primarily directed parallel to a series of elongate salt walls. In some mini‐basins, fluvial channel elements are stacked vertically within and along the central basin axes, in response to preferential salt withdrawal and resulting subsidence. In other basins, rim synclines have developed adjacent to bounding salt walls and these served as loci for accumulation of stacked fluvial channel complexes. Neighbouring mini‐basins exhibit different styles of infill at equivalent stratigraphic levels: sand‐poor basins dominated by fine‐grained, sheet‐like sandstone fluvial elements, which are representative of nonchannelised flow processes, apparently developed synchronously with neighbouring sand‐prone basins dominated by major fluvial channel‐belts, demonstrating effective partitioning of sediment route‐ways by surface topography generated by uplifting salt walls. Reworked gypsum clasts present in parts of the stratigraphy demonstrate the subaerial exposure of some salt walls, and their partial erosion and reworking into the fill of adjoining mini‐basins during accumulation of the Moenkopi Formation. Complex spatial changes in preserved stratigraphic thickness of four members in the Moenkopi Formation, both within and between mini‐basins, demonstrates a complex relationship between the location and timing of subsidence and the infill of the generated accommodation by fluvial processes.     

Quantifying sand delivery to deep water during changing sea-level: Numerical models from the Quaternary Brazos Icehouse continental margin

Sequence stratigraphy for clastic continental margins predicts the development of sand-rich turbidite deposits during specific times in relation to base-level cycles. It is now widely understood that deltas can extend to the shelf-edge forced by high sediment flux and/or base level, providing a direct connection to transfer sediment and sand to the slope and basin floor even during high base level periods. Herein, we build a stratigraphic forward model for the last 120 kyr of the fluvio-deltaic to deep-water Brazos system (USA) where sediment partitioning along an Icehouse continental margin can be evaluated. The reduced-complexity stratigraphic forward model employs geologically constrained input parameters and mass balance. The modelled architecture is consistent with the location of depositional units previously mapped in the shelf. Sand bypasses the shelf and upper slope between 35 to 15 kyr before present and only about 20%–30% of all the sediment and sand supplied to the system is transferred to deep water. Several scenarios based on the initial Brazos model investigate the relationships between base level and deep-water sand ratio (DWSR). DWSR is defined as the relative amount of sand transferred to the deep-water portions of the system subdivided by the total sand input to the model. Linear correlations between DWSR and base level change rates or base level are very poor. Short-term variability due to local processes (for example avulsions) is superimposed to the long-term trends and mask the base level signal. DWSR for an entire base-level cycle is mainly controlled by the proportion of time the delta stays docked at the shelf-edge. Stratigraphic forward models are useful to complement field observations and quantify how different processes control stratigraphy, which is important for making predictions in areas with limited information.