Reverse migration of lithofacies boundaries and shoreline in response to sea‐level rise

The migration of the lithofacies boundaries preserved in the sedimentary record is key to interpreting changes in depositional environments. Grain size is one of the most recognizable physical characteristics of lithofacies. The advance and retreat of grain‐size breaks, as a proxy for lithofacies boundaries (e.g. gravel–sand transition), is commonly attributed to variations in external controls (e.g. climate, sea level and tectonic subsidence). While most models of fluviodeltaic systems focus on predicting the response of the shoreline to these forcings, none have thoroughly incorporated the migration of grain‐size transitions (GST) that coevolve with the shoreline. We present a numerical delta evolution model that treats both the shoreline and GST as moving boundaries to provide quantitative understanding of the dynamic interaction between the downstream boundary (shoreline) and the upstream lithofacies boundaries (GSTs) of the fluviodeltaic system under relative sea‐level rise. We tested a range of relative sea‐level rise rates in the model. The shoreline and GST gradually reduced their progradation rates and eventually retreated landward as the fluviodeltaic topset and foreset elongated. However, their timings of retreat were different, resulting in a counterintuitive case for a quicker retreat of GST while the shoreline still continued to advance. A series of scaled flume experiments with a sand and crushed walnut sediment mixture captured the same behaviours of these two moving boundaries. We found that GST experienced higher relative sea‐level rise (RSLR) rates than the shoreline. This additional RSLR rate scales with the downstream river slope and the shoreline progradation rate to cause earlier GST retreat in comparison to the shoreline. The fundamental understanding from this study of migration of both the GST and shoreline in fluviodeltaic systems will aid in accurately assessing the trajectories of GST in sedimentary strata as a proxy for environmental change.     

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.     

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.     

Evolution and Morphodynamics of a Prograded Beach‐Ridge Foreland,Northern Baffin Island,Canadian Arctic Archipelago

Landward retreat (marine transgression) is a common response of coastal systems to rising relative sea level. However, given sufficient sediment supply, the coast may advance seaward. The latter response of gravel barriers has been recorded in parts of southeastern and northwestern Canada, where seaward‐rising sets of beach ridges are observed in areas of Holocene RSL rise. Cape Charles Yorke, northern Baffin Island, is a 5 km long gravel foreland characterized by seaward‐rising beach‐ridge crest elevations. The prograded morphology of the Cape Charles Yorke foreland is a prime example of coastal response to a combination of rising RSL and abundant sediment supply, an unusual and little‐documented pattern in the Canadian Arctic. The main gravel supply to Cape Charles Yorke is likely from eroding bedrock and raised marine deposits southwest of the foreland. Although not the dominant sediment source, the Cape Charles Yorke delta contributed to the formation of the foreland by sheltering it from easterly storm waves and providing an anchor point for the prograding ridges. The truncation of relict ridges by the modern shoreline suggests a recent regime shift from continuous deposition to predominant erosion. The cause and timing of this shift are unknown but could result from a recent dwindling in sediment supply, increased accommodation space, increased wave energy, and/or an accelerated rise of relative sea level.     

Tectono‐sedimentary development of early syn‐rift deposits: the Abura Graben,Suez Rift,Egypt

The thickness and distribution of early syn‐rift deposits record the evolution of structures accommodating the earliest phases of continental extension. However, our understanding of the detailed tectono‐sedimentary evolution of these deposits is poor, because in the subsurface, they are often deeply buried and below seismic resolution and sparsely sampled by borehole data. Furthermore, early syn‐rift deposits are typically poorly exposed in the field, being buried beneath thick, late syn‐rift and post‐rift deposits. To improve our understanding of the tectono‐sedimentary development of early syn‐rift strata during the initial stages of rifting, we examined quasi‐3D exposures in the Abura Graben, Suez Rift, Egypt. During the earliest stage of extension, forced folding above blind normal fault segments, rather than half‐graben formation adjacent to surface‐breaking faults, controlled rift physiography, accommodation development and the stratigraphic architecture of non‐marine, early syn‐rift deposits. Fluvial systems incised into underlying pre‐rift deposits and were structurally focused in the axis of the embryonic depocentre, which, at this time, was characterized by a fold‐bound syncline rather than a fault‐bound half graben. During this earliest phase of extension, sediment was sourced from the rift shoulder some 3 km to the NE of the depocentre, rather than from the crests of the flanking, intra‐basin extensional forced folds. Fault‐driven subsidence, perhaps augmented by a eustatic sea‐level rise, resulted in basin deepening and the deposition of a series of fluvial‐dominated mouth bars, which, like the preceding fluvial systems, were structurally pinned within the axis of the growing depocentre, which was still bound by extensional forced folds rather than faults. The extensional forced folds were eventually locally breached by surface‐breaking faults, resulting in the establishment of a half graben, basin deepening and the deposition of shallow marine sandstone and fan‐delta conglomerates. Because growth folding and faulting were coeval along‐strike, syn‐rift stratal units deposited at this time show a highly variable along‐strike stratigraphic architecture, locally thinning towards the growth fold but, only a few kilometres along‐strike, thickening towards the surface‐breaking fault. Despite displaying the classic early syn‐rift stratigraphic motif recording net upward‐deepening, extensional forced folding rather than surface faulting played a key role in controlling basin physiography, accommodation development, and syn‐rift stratal architecture and facies development during the early stages of extension. This structural and stratigraphic observations required to make this interpretation are relatively subtle and may go unrecognized in low‐resolution subsurface data sets.     

Sea‐level and ocean‐current control on carbonate‐platform growth,Maldives, Indian Ocean

Multichannel high‐resolution seismic and multibeam data were acquired from the Maldives‐isolated carbonate platform in the Indian Ocean for a detailed characterization of the Neogene bank architecture of this edifice. The goal of the research is to decipher the controlling factors of platform evolution, with a special emphasis on sea‐level changes and changes of the oceanic currents. The stacking pattern of Lower to Middle Miocene depositional sequences, with an evolution of a ramp geometry to a flat‐topped platform, reflects variations of accommodation, which here are proposed to be primarily governed by fluctuations of relative sea level. Easterly currents during this stage of bank growth controlled an asymmetric east‐directed progradation of the bank edge. During the late middle Miocene, this system was replaced by a twofold configuration of bank development. Bank growth continued synchronously with partial bank demise and associated sediment‐drift deposition. This turnover is attributed to the onset and/or intensification of the Indian monsoon and related upwelling and occurrence of currents, locally changing environmental conditions and impinging upon the carbonate system. Mega spill over lobes, shaped by reversing currents, formed as large‐scale prograding complexes, which have previously been interpreted as deposits formed during a forced regression. On a regional scale, a complex carbonate‐platform growth can occur, with a coexistence of bank‐margin progradation and aggradation, as well as partial drowning. It is further shown that a downward shift of clinoforms and offlapping geometries in carbonate platforms are not necessarily indicative for a sea‐level driven forced regression. Findings are expected to be applicable to other examples of Cenozoic platforms in the Indo‐Pacific region.     

Response of a proglacial delta to rapid high‐amplitude lake‐level change: an integration of outcrop data and high‐resolution shear wave seismics

In this paper, we will present the stratigraphic evolution, internal facies architecture and geomorphology of the Middle Pleistocene Emme delta, controlled by rapid high‐amplitude lake‐level change. The Emme delta was deposited on the northern margin of glacial Lake Weser, located in north‐west Germany. Rates of lake‐level rise were probably >50 mm year^?1 and rates of lake‐level fall 30–50 m within a few days or weeks, due to the opening of lake outlets. We use digital elevation models, sedimentology and shear wave seismics to improve earlier reconstructions and investigate the influence of rapid base‐level change on delta development. Shear wave seismic data resolve architectural elements in the range of metres and bridge the common gap between outcrop and conventional compression wave seismic data. The radial delta complex is about 2 km long, 1.8 km wide and up to 70 m thick, overlying a concave, up to 13° steep dipping ramp surface. It consists of vertically and laterally stacked delta lobes, caused by lobe switching during base‐level change. During the lake‐level rise, vertically stacked (Gilbert‐type) delta systems formed. The decrease in thickness and lateral extent indicates a rapid upslope shift of depocentres. A high rate and magnitude of lake‐level fall (50 m) promoted the development of a single incised valley and the deposition of forced regressive coarse‐grained delta lobes in front of the valley. The incised valley was filled during decreasing rates of lake‐level fall and low base‐level, because the alluvial gradient was larger than the emergent lake profile. Attached sand‐rich forced regressive aprons formed during lower magnitudes of lake‐level falls in the range of 30–35 m. Valley incision occurred, but was limited to the uppermost portion of the delta, controlled by the steep slope. The incised valley related to the final lake drainage is associated with long‐wavelength (60–90 m) bedforms at the downslope end, attributed to the formation of standing waves as a result of a hydraulic jump. Estimated palaeoflow depth during standing wave formation was ～9–14 m and flow velocity was 10–12 m s^?1. Because subsidence, waves or tides did not play a major role, the Emme delta can be used as an analogue‐based predictive stratigraphical and sedimentological model for steep glacigenic deltas controlled by rapid base‐level change and can help to understand better the facies distribution and three‐dimensional geometry of these depositional systems.     

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.     

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.     

Timing of depth‐dependent lithosphere stretching on the S. Lofoten rifted margin offshore mid‐Norway: pre‐breakup or post‐breakup?

Depth‐dependent stretching, in which whole‐crustal and whole‐lithosphere extension is significantly greater than upper‐crustal extension, has been observed at both non‐volcanic and volcanic rifted continental margins. A key question is whether depth‐dependent stretching occurs during pre‐breakup rifting or during sea‐floor spreading initiation and early sea‐floor spreading. Analysis of post‐breakup thermal subsidence and upper‐crustal faulting show that depth‐dependent lithosphere stretching occurs on the outer part of the Norwegian volcanic rifted margin. For the southern Lofoten margin, large breakup lithosphere β stretching factors approaching infinity are required within 100 km of the continent–ocean boundary to restore Lower Eocene sediments and flood basalt surfaces (～54 Ma) to interpreted sub‐aerial depositional environments at sea level as indicated by well data. For the same region, the upper crust shows no significant Palaeocene and Late Cretaceous faulting preceding breakup with upper‐crustal β stretching factors <1.05. Further north on the Lofoten margin, reverse modelling of post‐breakup subsidence with a β stretching factor of infinity predicts palaeo‐bathymetries of ～1500 m to the west of the Utrøst Ridge and fails to restore Lower Eocene sediments and flood basalt tops to sea level at ～54 Ma. If these horizons were deposited in a sub‐aerial depositional environment, as indicated by well data to the south, an additional subsidence event younger than 54 Ma is required compatible with lower‐crustal thinning during sea‐floor spreading initiation. For the northern Vøring margin, breakup lithosphere β stretching factors of ～2.5 are required to restore Lower Eocene sediments and basalts to sea level at deposition, while Palaeocene and Late Cretaceous upper‐crustal β stretching factors for the same region are < 1.1. The absence of significant Palaeocene and late Cretaceous extension on the southern Lofoten and northern Vøring margins prior to continental breakup supports the hypothesis that depth‐dependent stretching of rifted margin lithosphere occurs during sea‐floor spreading initiation or early sea‐floor spreading rather than during pre‐breakup rifting.     

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.     

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.     

Non‐unique stratal geometries: implications for sequence stratigraphic interpretations

There is now strong evidence that stratal geometries on basin margins are most likely a consequence of multiple controls, not just variations in accommodation. Consequently, correct sequence stratigraphic interpretation of stratal geometries requires an understanding of how multiple different controls may generate similar geometries. Using a simple numerical stratigraphic forward model, we explore the impact of time variable sediment supply and different sediment transport rates on stratal geometries. We demonstrate how four common types of stratal geometry can form by more than one set of controlling parameter values and are thus likely to be non‐unique, meaning that there may be several sets of controlling factors that can plausibly explain their formation. For example, a maximum transgressive surface can occur in the model due to an increase in rate of relative sea‐level rise during constant sediment supply, and due to a reduction in rate of sediment supply during a constant rate of relative sea‐level rise. Sequence boundaries, topset aggradation and shoreline trajectories are also examples of non‐unique stratal geometries. If the model simulations in this work are sufficiently realistic, then the modelled stratal geometries are important examples of non‐uniqueness, suggesting the need for a shift towards sequence stratigraphic methods based on constructing and evaluating multiple hypotheses and scenarios.     

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.     

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.     

Tectonic controls on Cenozoic foreland basin development in the north‐eastern Andes,Colombia

In order to evaluate the relationship between thrust loading and sedimentary facies evolution, we analyse the progradation of fluvial coarse‐grained deposits in the retroarc foreland basin system of the northern Andes of Colombia. We compare the observed sedimentary facies distribution with the calculated one‐dimensional (1D) Eocene to Quaternary sediment‐accumulation rates in the Medina wedge‐top basin and with a three‐dimensional (3D) sedimentary budget based on the interpretation of ～1800 km of industry‐style seismic reflection profiles and borehole data. Age constraints are derived from a new chronostratigraphic framework based on extensive fossil palynological assemblages. The sedimentological data from the Medina Basin reveal rapid accumulation of fluvial and lacustrine sediments at rates of up to ～500 m my^?1 during the Miocene. Provenance data based on gravel petrography and paleocurrents reveal that these Miocene fluvial systems were sourced from Upper Cretaceous and Paleocene sedimentary units exposed to the west in the Eastern Cordillera. Peak sediment‐accumulation rates in the upper Carbonera Formation and the Guayabo Group occur during episodes of coarse‐grained facies progradation in the early and late Miocene proximal foredeep. We interpret this positive correlation between sediment accumulation and gravel deposition as the direct consequence of thrust activity along the Servitá–Lengupá faults. This contrasts with one class of models relating gravel progradation in more distal portions of foreland basin systems to episodes of tectonic quiescence.     

Compound and hybrid clinothems of the last lowstand Mid-Adriatic Deep: Processes,depositional environments,controls and implications for stratigraphic analysis of prograding systems

Clinoforms with a range of scales are essential elements of prograding continental margins. Different types of clinoforms develop during margin growth, depending on combined changes in relative sea level, sediment supply and oceanographic processes. In studies of continental margin stratigraphy, trajectories of clinoform ‘rollover’ points are often used as proxies for relative sea-level variation and as predictors of the character of deposits beyond the shelf-break. The analysis of clinoform dynamics and rollover trajectory often suffers from the low resolution of geophysical data, the small scale of outcrops with respect to the dimensions of clinoform packages and low chronostratigraphic resolution. Here, through high-resolution seismic reflection data and sediment cores, we show how compound clinoforms were the most common architectural style of margin progradation of the late Pleistocene lowstand in the Adriatic Sea. During compound clinoform development, the shoreline was located landward of the shelf-break. It comprised a wave-dominated delta to the west and a barrier and back-barrier depositional system in the central and eastern area. Storm-enhanced hyperpycnal flows were responsible for the deposition of a sandy lobe in the river mouth, whereas a heterolithic succession formed elsewhere on the shelf. The storm-enhanced hyperpycnal flows built an apron of sand on the slope that interrupted an otherwise homogeneous progradational mudbelt. Locally, the late lowstand compound clinoforms have a flat to falling shelf-break trajectory. However, the main phase of shelf-break bypass and basin deposition coincides with a younger steeply rising shelf-break trajectory. We interpret divergence from standard models, linking shelf-break trajectory to deep-sea sand deposition, as resulting from a great efficiency of oceanographic processes in reworking sediment in the shelf, and from a high sediment supply. The slope foresets had a large progradational attitude during the late lowstand sea-level rise, showing that oceanographic processes can inhibit coastal systems to reach the shelf-edge. In general, our study suggests that where the shoreline does not coincide with the shelf-break, trajectory analysis can lead to inaccurate reconstruction of the depositional history of a margin.     

Continental margin syn‐rift salt tectonics at intermediate width margins

This paper examines interactions among syn‐rift continental margin extension, evaporites, particularly rocksalt (halite), deposited in the overlying sedimentary basins, and clastic sediment loading. We present dynamically evolving 2D numerical models that combine syn‐rift lithospheric extension, with salt (viscous halite, 10¹⁸–10¹⁹ Pa s) and clastic (frictional‐plastic) sediment deposition to investigate how salt is distributed and subsequently mobilized during syn‐rift extension. Example results are shown, contrasting salt deposition in the early, mid and late syn‐rift phases of a single lithospheric extension model. The lithospheric model is chosen to give depth‐dependent extension and intermediate width margins with proximal grabens and a hyperextended distal region. The models exhibit diachronous migration of extension towards the rift axis and this is reflected in the faulting of overlying sediments. The models illustrate the roles of timing of salt deposition, relative to rifting and subsequent sedimentation, in defining the location and deformation of syn‐rift salt, with post‐salt sediment progradation in some models. Late deposition of salt leads to increased lateral extent of the original salt body and decreased variation in salt thickness. Seaward flow of salt increases with later deposition; early syn‐rift salt is deposited and trapped in the grabens, whereas mid and late syn‐rift salt tends to flow towards the distal margin or even over the oceanic crust. Prograding clastic post‐salt sediments drive more substantial seaward movement of mid and late syn‐rift salt. A numerical model of the Red Sea with evaporite deposition during the mid to late syn‐rift period, preceded and followed by aggrading and prograding clastic sediment, shows reasonable agreement with observations from the central Red Sea.     

Tidal amplification and along‐strike process variability in a mixed‐energy paralic system prograding onto a low accommodation shelf,Edgeøya,Svalbard

The study describes the depositional development and sediment partitioning in a prograding paralic Triassic succession. The deposits are associated with the advance of large prism‐scale clinoforms across a shallower platform area. Approaching the platform, the limited accommodation and associated relative higher rates of deposition generated straighter clinoforms with lower foreset angles. The vertical restriction across the platform is interpreted to have amplified the tidal signature. Sediment was redistributed from the coast into increasingly sandy delta‐front deposits, compared to offshore equivalents. The deposits comprise extensive compound dune fields of amalgamated and increasingly clean sandbodies up‐section. Rapid deposition of significant amounts of sand led to differential subsidence and growth‐faulting in the delta front, with downthrown fault blocks further amplifying the tidal energy through funnelling. A mixed‐energy environment created along‐strike variability along the delta front with sedimentation governing process‐regime. Areas of lower sedimentation were reworked by wave and storm‐action, whereas high sedimentation rates preserved fluvially dominated mouth bars. A major transgression, however, favoured tidally dominated deposits also in these areas, attributed to increasing rugosity of the coastline. Formation of an extensive subaqueous platform between the coast and delta front dampened incoming wave energy, and tidally dominated deposits dominate the near‐shore successions. Meanwhile formation of wave‐built sand‐bars atop the platform attest to continued wave influence. The strong tidal regime led to the development of a heterolithic near‐shore tidally dominated channel system, and sandier fluvial channels up‐river. The highly meandering tidal channels incising the subaqueous platform form kilometre wide successions of inclined heterolithic stratification. The fluvially dominated channels which govern deposition on the delta plain are narrower and slightly less deep, straighter, generally symmetric and filled with cleaner sands. This study provides important insight into tidal amplification and sand redistribution during shallowing on a wide shelf, along with along‐strike process‐regime variability resulting from variations in sediment influx.     

Microbial‐dominated carbonate platforms during the Ladinian rifting: sequence stratigraphy and evolution of accommodation in a fault‐controlled setting (Catalan Coastal Ranges,NE Spain)

The Upper Muschelkalk sedimentary record constitutes a major transgressive pulse of north‐eastern Iberia during the Ladinian. This record is arranged in two transgressive–regressive (T–R) sequences formed by two stepped microbial‐dominated carbonate ramp systems where accommodation was mainly controlled by extensional faults. This study seeks to gain new insights into how the evolution of syn‐rift subsidence controls the creation of accommodation space, the depositional styles and, especially, the palaeogeographical domains where specific microbialites developed (thrombolites and stromatolites). Thrombolite bodies (at least 40 m thick) display two types of architecture, biostromal and mud‐mounded and stromatolite bodies (at least 7 m thick) consist of tabular and domed, head‐shaped morphologies. Domed and mounded forms are usually developed during stages of increasing accommodation rates, low‐to flat‐nelief forms tend to grow in association with periods of low accommodation rates. A sea‐level fall of at least 50 m occurred at the end of the Early Ladinian leaving the platform subaerially exposed. As a result, a prominent karst with significant erosional incisions and profuse collapse breccia fillings was formed in the inner and middle ramp settings. The resultant subaerial unconformity bounds T–R sequences 1 and 2. Subsidence curves display two stages of rapid/decelerated total subsidence, constituting two discrete rift/post‐rift pulses in the large Triassic rifting period: (i) Buntsandstein – Middle Muschelkalk, and (ii) Late Muschelkalk – Imon Formation (Rhaetian). The second pulse is characterized by a rapid syn‐rift subsidence during the Late Muschelkalk, and a decelerated post‐rift subsidence throughout the deposition of Keuper facies and Imon Formation. The Late Muschelkalk rapid syn‐rift pulse of total subsidence produces gains in accommodation, which controls the development of the stromatolites and thrombolites (biostromes and mud‐mounds).