Ephemeral rollover points and clinothem evolution in the modern Po Delta based on repeated bathymetric surveys

Reconstructions of ancient delta systems rely typically on a two-dimensional (2D) view of prograding clinothems but may miss their three-dimensional (3D) stratigraphic complexity which can, instead, be best documented on modern delta systems by integrating high-resolution geophysical data, historical cartography, core data and geomorphological reconstructions offshore. We quantitatively compare three precisely positioned, high-resolution multi-beam bathymetry maps in the delta front and pro delta sectors (0.3 to 10 m water depth) of Po di Pila, the most active of the modern Po Delta five branches. By investigating the detailed morphology of the prograding modern Po Delta, we shed new light on the mechanisms that control the topset to foreset transition in clinothems and show the temporal and spatial complexity of a delta and its pro delta slope, under the impact of oceanographic processes. This study documents the ephemeral nature of the rollover point at the transition between sandy topset (fluvial, delta plain to mouth-bar) and muddy seaward-dipping foreset deposits advancing, in this case, in >20 m of water depth. Three multibeam surveys, acquired between 2013 and 2016, document the complexity in space and time of the topset and foreset regions and their related morphology, a diagnostic feature that could not be appreciated using solely 2D, even very high-resolution, seismic profiles. In addition, the comparison of bathymetric surveys gathered with one-year lapses shows the migration of subaqueous sand dunes on the clinothem topset, the formation of ephemeral cut-and-fill features at the rollover point (few m below mean sea level), the presence of collapse depressions derived by sagging of sediments and fluid expulsion (possibly induced by storm waves) on the foreset, and splays of sand likely reflecting gravity flows on the lower foreset. Though the modern Po Delta is anthropogenic in many respects, its subaqueous clinothem can be studied as a scale model for ancient clinothems that are less resolved geometrically and far less constrained chronologically.     

Shallow-water deltaic clinoforms and process regime

Utilizing two outcrop data sets with dip direction exposures of shallow-water (tens of meters) deltaic clinoforms, this paper quantifies sedimentary facies proportions and clinoform lengths and gradients, and links process regimes to delta clinoform dimensions. Both data sets are from foreland basins, the Cretaceous Chimney Rock Sandstone of the Rock Springs Formation from the US Western Interior, and the Eocene Brogniartfjellet Clinoform Complex 8 of the Battfjellet Formation from the Central Basin of Spitsbergen. Sedimentary facies indicate presence of both river- and wave-dominated clinothems in each data set. Facies characteristics and distribution implies that river-dominated clinothem progradation was primarily driven by deposition from weak hyperpycnal flow turbidity currents across the clinoforms, and minor slumps. Wave-dominated clinothems were constructed by wave processes rather than alongshore currents, and are also progradational subaerial clinoforms, with one exception, where the formation of a compound subaqueous clinoform set indicates erosion and sediment bypass above the wave base. Sediment distribution and lithological heterogeneity in the river-dominated clinothems is controlled by individual hyperpycnal flow events or mouth-bar collapse events, and thus by self-organization and minimal reworking that results in a heterogeneity that is difficult to predict (high entropy). The efficient reworking of river-derived sediments in wave-dominated clinothems results in predictable lithological sediment partitioning (low entropy). Clinoform dimension analyses show that although of similar sediment caliber, river-dominated clinoforms in both data sets are on average 3–4 times steeper and 3–4 times shorter than the wave-dominated clinoforms, with mean gradients of ca 4 degrees and ca 1 degree, respectively, and mean lengths of 150–230 m and 640–760 m. These results require corroboration from additional data sets, but do suggest that river- and wave-dominated delta clinoforms are likely to have distinct downdip extents (lengths) and gradients for given clinoform heights. Clinoform shape can thus be a method for differentiating ancient river- vs. wave-dominated deltaic clinoforms, in addition to their sedimentary facies, biogenic features and sandstone maturity, and helpful when incorporated into reservoir models.     

Upslope-climbing shelf-edge clinoforms and the stepwise evolution of the northern European glaciation (lower Pleistocene Eridanos Delta system,U.K. North Sea): When sediment supply overwhelms accommodation

Clinoforms are basinward-dipping and accreting palaeo-bathymetric profiles that record palaeo-environmental conditions and processes; thus, clinothems represent natural palaeo-archives. Here, we document shelf-edge scale clinoform sets which prograded through the entire width of an epicontinental marine basin (ca. 400 km), eventually encroaching onto the opposite basin flank, where they started to prograde upslope and landward, in defiance of gravity (“upslope-climbing clinoforms”). The giant westward-prograding Eridanos muddy shelf-edge clinothem originated from the Baltic hinterland in the Oligocene and achieved maximum regression in the Early Pleistocene, on the UK Central Graben (CG) and Mid North Sea High (MNSH), after crossing the whole North Sea mesopelagic depocentre and causing near complete basin infill. Here we integrate well and seismic data through the MNSH and CG and examine the Eridanos final heyday and demise, identifying five clinothem complexes (A1, A2, A3, B and C) and six depositional sequence boundaries (SB1 to SB6) in the Miocene-Recent section. Tectonic and climatic events drove the recent evolution of this system. Early Pleistocene climate cooling, in particular, resulted in a stepwise increase in sediment supply. This climaxed in the earliest Calabrian, following a likely Eburonian eustatic fall (=SB3) when the Eridanos clastic wedge was restructured from a 100–300 m thick compound shelf-edge and delta system to a “hybrid” shelf-edge delta at sequence boundary SB3 (ca. 1.75 Ma). In the ca. 40 kyr that followed SB3, a progradation rate peak (>1,000 m/kyr) is associated with clinoforms starting to accrete upslope, onto the east-dipping slope between CG and MNSH. This “upslope-climbing clinoform” phase was quickly followed by the maximum regression and final retreat of the Eridanos system in the Early Calabrian (=SB4), likely as the result of climate-driven changes in the Baltic hinterland and/or delta auto-retreat. To our knowledge, this contributions represents the first documentation of “upslope-climbing clinoforms” recorded in the stratigraphic record.     

Climatically controlled river terrace staircases: A worldwide Quaternary phenomenon

A comparison of fluvial terrace sequences from around the world, based on data collected as part of International Geoscience Programme (IGCP) Project No. 449, has revealed significant patterns. River terraces provide important records of uplift, which is essential for their formation, and of landscape evolution. Their cyclic formation, however, almost invariably seems to have been a response to climatic fluctuation. Sequences in the European core area of IGCP 449, which has the longest and most extensive research history, have been used as templates for worldwide comparison. There is evidence for a global acceleration of uplift at the time of, and perhaps in response to, the Mid-Pleistocene Revolution, when climatic fluctuation switched to 100 kyr Milankovitch cycles. Terraces formed prior to this generally consist of wide aggradational sheets that probably each represent formation over several 41 kyr cycles. Subsequently, river valleys became more steeply entrenched and terraces formed in response to the stronger 100 kyr climatic forcing, in many cases at approximately one per cycle. This paper uses the new data resource to explore differences between records in different climate zones, between sequences with variable numbers of Middle–Late Pleistocene terraces and between systems in which the all-important incision event has occurred in different parts of climatic cycles. Key records discussed include European examples from the Rhine, Thames, Somme, Dniester, Dnieper, Don, Volga and Aguas; from Asia the Gediz (Turkey) and Orontes (Syria); from North America, the South Platte and Colorado; from South Africa the Vaal and Sundays; from Australia the Shoalhaven; and from South America, the Amazon, Paraguay and tributaries of the Colorado and Negro.     

The environmental significance of magnetic measurements of Late Pleistocene and Holocene sediments from the Grand Lac d'Annecy, eastern France

The results of a series of magnetic measurements on sediment samples are illustrated by means of examples selected from a more specialized parallel publication (Hu & Oldfield, submitted). They are interpreted in terms of magnetic components, notably, 'magnetite' of varying grain size/domain state, haematite, goethite and greigite. Wherever possible, some indication is given of changes in the relative importance of the various magnetic grain size and mineral phases through time. The main components of the magnetic assemblage are described for each of the main lithological units identified in the long core which is considered elsewhere in this volume. These assemblages are related to source areas and to depositional processes, especially in the pre-Holocene parts of the core. In the Holocene sediments, the magnetic properties cannot be interpreted solely in terms of shifts in detrital mineral sources. Evidence is presented for the presence of both biogenic magnetite and authigenic greigite. The most recent sediments point to renewed catchment erosion, probably as a result of human activity during historical times.     

On the use of geochronology of detrital grains in determining the time of deposition of clastic sedimentary strata

Placing geologic events in a temporal framework is essential to telling the story of Earth history. However, clastic sedimentary rocks can be difficult to date in an absolute reference because they are made up of grains that are older than the rock in which they are now found, and some clastic rocks do not contain fossils that allow precise reference to the Geologic Timescale. For such rocks, the isotopic dating of detrital minerals can be used to estimate the time of deposition; the clastic rock must be younger than the youngest grain analysed. However, many researchers eschew this simple and straightforward approach in favour of schemes that estimate the maximum allowable depositional age as the weighted mean of the age of several grains, chosen by a variety of selection criteria. This is a mistake; in the absence of a geochemical resemblance apart from the similarity of their age, detrital grains should not be assumed to have originated in the same system and therefore any averaging or other manipulation of such data is statistically invalid and produces results without geologic significance. In the absence of interbedded volcanic rocks or index fossils, dating of detrital minerals can be an important aid in understanding the time of deposition of clastic rocks, but the best estimate will come from taking note of the youngest single grain and not by inappropriately averaging data.     

Metals and models: Diagenic modelling in freshwater lacustrine sediments

This review aims to distill and synthesize the existing information on the use of models to describe and predict the distribution and movement of metals in lacustrine sediments. As such it examines the causes of metal diagenesis, the origin and form of the equations that govern these phenomena, and the predictability or measurability of the parameters that appear in the models. The paper concludes by highlighting some seminal results from modelling studies, including the determination of the factors controlling the formation of surficial Mn-Fe-enriched zone or layers, the substantial contribution possible from metal reduction to organic matter regeneration, the calculation of mixing-corrected metal input histories to lakes, and the prediction of growth rates and morphologies for both deep-sea and lacustrine ferromanganese nodules.     

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.     

Stratigraphic complexity in fluvial fans: Lower Eocene Green River Formation,Uinta Basin,USA

In this study, measured outcrop sections and geolocated photomosaics are integrated with areal mapping of channel dimensions, degree of amalgamation, calculations of channel‐to‐floodplain ratios and sedimentary facies variability to study and quantify the channel and floodplain deposits in the Sunnyside Delta Interval of the Lower Eocene Green River Formation in the Uinta Basin, Utah. Vertically, sand content and bed thickness increases, due to an increase in the channel‐to‐floodplain ratio, channel size and the degree of channel amalgamation. Laterally, the channel‐to‐floodplain ratio, channel size, the degree of channel amalgamation and the sand content in channel facies decreases in the paleo‐downstream direction. Such vertical and lateral transitions identify the Sunnyside Delta Interval as a fluvial fan (or distributive fluvial system). However, the vertical and lateral transitions occur at multiple spatial scales, demonstrating considerable stratigraphic complexity as compared to the existing facies and architectural models suggested for fluvial megafans and distributive fluvial systems. The smallest‐scale transitions are identified as avulsion‐related packages that form the building blocks of the stratigraphy, whereas the intermediate‐ and largest‐scale transitions are suggested to be related to lobe and whole fan progradation respectively. This documented complexity indicates the significance of self‐organization in building fluvial fan stratigraphy, and demonstrates that changes in the degree of channel amalgamation or in channel‐to‐floodplain ratio are not linked to accommodation changes. On facies scale, an abundance of Froude supercritical‐flow and high‐deposition‐rate facies, in‐channel mud deposits, and in‐channel bioturbation and desiccation indicate deposition in rivers with highly variable discharge. Such discharge conditions suggest seasonally and inter‐annually variable precipitation conditions in the US Western Interior in the Early Eocene.     

Clinoforms as paleogeographic tools: Development of the Danube catchment above the deep Paratethyan basins in Central and Southeast Europe

The Miocene marine basins of Central and Southeast Europe, once comprising the Paratethys Sea, were gradually filled with sediments during the Neogene and turned to be the catchment area of the proto-Danube and finally that of the modern Danube. Seismic data from various parts of the large Danube catchment area show that these several hundred meter deep basins were filled by lateral accretion of river-transported sediments, appearing as shelf edge scale clinoform sets in seismic profiles. The direction of shelf edge progradation is NW to SE (N to S, W to E) in each basin, except for the Dacian basin where NE to SW direction prevails. The age of the clinoform sets is generally younging downstream: 19–18 Ma in the North Alpine Foreland basin, 14–13 Ma in the Vienna basin, 10–9 Ma in the Danube (Kisalföld) basin, 8.6–4 Ma in the Central Pannonian basin (Alföld), ?9–5 Ma in the Dacian basin, and 6–0 Ma in the Euxinian (Black Sea) basin. In spite of this geographical and temporal pattern, only the Danube (Kisalföld) and the western and central part of the Central Pannonian basin were filled by the proto-Danube shelf accretion. Formation of the Danube, as a longitudinal river of the Alpine foreland that gradually elongated to the east and followed the retreating shoreline of the Paratethys, most probably took place at the beginning of the Late Miocene, ca. 11 Ma ago, thus the Early and Middle Miocene shelf advance in the North Alpine Foreland and Vienna basins, respectively, cannot be attributed to a „paleo-Danube”. The clinoform systems of the Dacian basin are coeval with those of the upstream Central Pannonian basin, indicating that by the time the Danube sedimentary system reached the Dacian basin, it was already a shallow basin. The vast clinoforms of the northwestern Euxinian shelf also significantly overlap in age with the Pannonian basin ones; only the <4 Ma part of the shelf accretion can be attributed to the Danube sensu stricto.     

Evaluating alluvial stratigraphic response to cyclic and non‐cyclic upstream forcing through process‐based alluvial architecture modelling

Formation of alluvial stratigraphy is controlled by autogenic processes that mix their imprints with allogenic forcing. In some alluvial successions, sedimentary cycles have been linked to astronomically‐driven, cyclic climate changes. However, it remains challenging to define how such cyclic allogenic forcing leads to sedimentary cycles when it continuously occurs in concert with autogenic forcing. Accordingly, we evaluate the impact of cyclic and non‐cyclic upstream forcing on alluvial stratigraphy through a process‐based alluvial architecture model, the Karssenberg and Bridge (2008) model (KB08). The KB08 model depicts diffusion‐based sediment transport, erosion and deposition within a network of channel belts and associated floodplains, with river avulsion dependent on lateral floodplain gradient, flood magnitude and frequency, and stochastic components. We find cyclic alluvial stratigraphic patterns to occur when there is cyclicity in the ratio of sediment supply over water discharge (Q_s/Q_w ratio), in the precondition that the allogenic forcing has sufficiently large amplitudes and long, but not very long, wavelengths, depending on inherent properties of the modelled basin (e.g. basin subsidence, size, and slope). Each alluvial stratigraphic cycle consists of two phases: an aggradation phase characterized by rapid sedimentation due to frequent channel shifting and a non‐deposition phase characterized by channel belt stability and, depending on Q_s/Q_w amplitudes, incision. Larger Q_s/Q_w ratio amplitudes contribute to weaker downstream signal shredding by stochastic components in the model. Floodplain topographic differences are found to be compensated by autogenic dynamics at certain compensational timescales in fully autogenic runs, while the presence of allogenic forcing clearly impacts the compensational stacking patterns.     

Drainage integration and sediment dispersal in active continental rifts: A numerical modelling study of the central Italian Apennines

Progressive integration of drainage networks during active crustal extension is observed in continental areas around the globe. This phenomenon is often explained in terms of headward erosion, controlled by the distance to an external base‐level (e.g. the coast). However, conclusive field evidence for the mechanism(s) driving integration is commonly absent as drainage integration events are generally followed by strong erosion. Based on a numerical modelling study of the actively extending central Italian Apennines, we show that overspill mechanisms (basin overfilling and lake overspill) are more likely mechanisms for driving drainage integration in extensional settings and that the balance between sediment supply vs. accommodation creation in fault‐bounded basins is of key importance. In this area drainage integration is evidenced by lake disappearance since the early Pleistocene and the transition from internal (endorheic) to external drainage, i.e. connected to the coast. Using field observations from the central Apennines, we constrain normal faulting and regional surface uplift within the surface process model CASCADE (Braun & Sambridge, 1997, Basin Research, 9, 27) and demonstrate the phenomenon of drainage integration, showing how it leads to the gradual disappearance of lakes and the transition to an interconnected fluvial transport system over time. Our model results show that, in the central Apennines, the relief generated through both regional uplift and fault‐block uplift produces sufficient sediment to fill the extensional basins, enabling overspill and individual basins to eventually become fluvially connected. We discuss field observations that support our findings and throw new light upon previously published interpretations of landscape evolution in this area. We also evaluate the implications of drainage integration for topographic development, regional sediment dispersal and offshore sediment supply. Finally, we discuss the applicability of our results to other continental rifts (including those where regional uplift is absent) and the importance of drainage integration for transient landscape evolution.     

The lacustrine depression at Annecy (France), geological setting and Quaternary evolution

The Annecy lacustrine depression lies in a tectonic transverse valley of structural origin, linking the present day towns of Annecy and Ugine and bounded by the Bauges and Bornes massifs. Only the last two major glacial periods are represented in the stratigraphic and geomorphic evidence from the area, but its successive glaciations that have been largely responsible for the form of the valley within which the present day lake lies. Of particular importance were the lobes of the Arve glacier, one shaping the depression from the south-east; the other from the north. The retreat of the glaciers left behind a more extensive precursor of the present day Lac d'Annecy, which latter constitutes the residual body of open water after lake level lowering, partial infilling and drainage diversions within the catchment. The drainage basin of the extended lake formed shortly after the last glacial maximum (LGM) in the region. Its extent can be determined from field sampling and coring and its elevation can be established at 460 m above present day sea level. The present day lake is more isolated, since the Fier river, to the north, no longer drains into it. To the southeast, the course of the Chaise river has also been diverted, leaving the lake's drainage basin reduced to half its post-LGM size. The major paleogeographical events in the Late-Quaternary evolution of the basin are described in chronological sequence.     

Nested intrashelf platform clinoforms—Evidence of shelf platform growth exemplified by Lower Cretaceous strata in the Barents Sea

Two nested clinoform set types of different scales and steepness are mapped and analysed from high-resolution seismic data. Restoration of post-depositional faulting reveals a persistent pattern of small-scale, high-angle clinoforms contained within platform-scale, low-angle clinothems, showing a combined overall progradational depositional system. The large clinoforms lack a well-defined platform edge, and show a gradual increase in dip from topset to foreset. A consistent recurring stratal pattern is evident from the architecture, and is considered a result of interplay between relative sea-level change and autocyclic switching of sediment delivery focal points that brought sediment to the platform edge. This un-interrupted succession records how intra-shelf platforms prograde. Quantitative clinoform analysis may assist in determining the most influential depositional factors. Post-depositional uplift and erosion requires restoration with re-burial to maximum burial depth. Backstripping, decompaction and isostatic correction was performed assuming a range of lithologic compositions, as no wells test the lithology. Nearby wells penetrate strata basinward of the clinoforms, proving mudstone content above 50%, which in turn guide restoration values. Typical restored platform heights are 250–300 m, with correspondingly sized platform-scale clinoform heights. Typical large-scale clinoform foreset dip values are 1.3°–2.4°. Small-scale clinothems are typically 100 m thick, with restored foreset dip angles at 4.4° - > 10°. The results suggest that intrashelf platform growth occurs in pulses interrupted by draping of strata over its clinoform profile. The resultant architecture comprises small-scale clinoforms nested within platform-scale clinothems.