Ian A. Lunt

Can we distinguish flood frequency and magnitude in the sedimentological record of rivers

Consideration of the origin of alluvial deposits and their paleoenvironmental interpretation has traditionally involved two schools of thought: that they are either the result of processes that, on average, have acted uniformly through time, or that they are related to exceptional events that occur infrequently. Despite the long-running debate of gradualism versus catastrophism within the Earth Sciences, there are surprisingly few quantitative data to assess the magnitude of events that produce alluvial sedimentary successions. This paper reports on a unique ‘natural experiment’ where surface (digital elevation model, DEM) and subsurface (ground penetrating radar, GPR) data were taken immediately prior to, and after, a large (1-in-40 yr) flood event on the sandy, braided, South Saskatchewan River, Canada. Results show that although this high-magnitude flood reworked the entire braidplain, the scale of scour and style of deposition was similar to that associated with lower-magnitude, annual, floods. The absence of a distinct imprint of this large flood within the deposits is related to the fact that as river discharge rises, and begins to flow overbank, channel width increases at a much faster rate than flow depth, and thus the rate of increase in channel bed shear stress declines. Hence, rather than being a product of either frequent or rare events, alluvial deposits will likely be created by a range of different magnitude floods, but discriminating between these different scale events in the rock record may be extremely difficult.

Source-to-sink systems on passive margins: theory and practice with an example from the Norwegian continental margin

Abstract Source-to-sink system analysis involves a complete, earth systems model approach from the ultimate onshore drainage point to the toe of related active deepwater sedimentary systems. Several methods and techniques have evolved in recent years, from experimental and numerical modelling through analysis of modern and recent systems, to analysis of ancient systems. A novel method has been developed, bridging between the previous approaches and dividing and analysing source-to-sink systems based on linked geomorphic segments along the source-to-sink profile. This approach builds on uniformitarian principles. The method is driven by the need to understand ancient, subsurface systems and still has high uncertainty but is an original, first-order approach to source-to-sink system analysis. In modern systems, entire onshore-to-offshore systems can be analysed with a higher degree of confidence than in ancient systems and semi-quantitative relationships can be established. Application in ancient systems is much more challenging but, in some cases, antecedent morphologies have been preserved onshore that can be matched with offshore known occurrences of, for instance, sandy submarine fan systems. Along the Norwegian North Sea and Norwegian Sea margins the Paleocene deep-marine reservoir of the giant Ormen Lange gas field is such an example. There, antecedent onshore drainage patterns which formed the feeder system to the offshore, deepwater fan system can be interpreted and aligned with onshore palaeogeomorphological evidence. Understanding the palaeogeomorphic development of basement regions such as the Fennoscandian shield is of high importance for understanding the offshore presence of deepwater sandstones.

Linking offshore stratigraphy to onshore paleotopography: The Late Jurassic–Paleocene evolution of the south Norwegian margin

The temporal link between onshore topography and offshore stratigraphy is of key importance for understanding the long-term development of continental margins. The Mesozoic and Cenozoic landscape topography of southern Norway has long been a matter of debate due to the absence of concrete data reflecting ancient onshore relief. We here present quantitative estimates of paleotopography in southern Norway during three key time intervals (Late Jurassic, Late Cretaceous, and Paleocene) based on the observed volume of sediment stored in point-sourced depocenters along the margin. The aim is to determine the syndepositional paleolandscapes that are in best agreement with the observed distribution, geometry, lithology, and volume of different offshore depositional units. Probability ranges of estimated paleorelief are based on Monte Carlo simulations of a relief prediction model. The results suggest that the late Phanerozoic sedimentary successions along the margin best reflect: (1) a Late Jurassic landscape with a maximum relief of ∼1.6 km; (2) a more subdued Late Cretaceous topography of <0.5 km in the south and up to ∼0.5 km in the north; and (3) a Paleocene topography of up to ∼1.1 km in the north and ∼0.5 km in the south. The temporal and spatial variability in offshore deposition and the linked onshore topography strongly suggest that relief must have been rejuvenated several times during the late Phanerozoic in response to tectonic activity along major fault systems. This quantitative assessment not only sheds new light on the topographic evolution of the south Norwegian margin, it also provides a potential tool with which to investigate the topographic evolution of continental margins and the close coupling between onshore relief and offshore stratigraphy in nearby sedimentary basins.

Deposits of the sandy braided South Saskatchewan River: Implications for the use of modern analogs in reconstructing channel dimensions in reservoir characterization

Estimation of the dimensions of fluvial geobodies from core data is a notoriously difficult problem in reservoir modeling. To try and improve such estimates and, hence, reduce uncertainty in geomodels, data on dunes, unit bars, cross-bar channels, and compound bars and their associated deposits are presented herein from the sand-bed braided South Saskatchewan River, Canada. These data are used to test models that relate the scale of the formative bed forms to the dimensions of the preserved deposits and, therefore, provide an insight as to how such deposits may be preserved over geologic time. The preservation of bed-form geometry is quantified by comparing the alluvial architecture above and below the maximum erosion depth of the modern channel deposits. This comparison shows that there is no significant difference in the mean set thickness of dune cross-strata above and below the basal erosion surface of the contemporary channel, thus suggesting that dimensional relationships between dune deposits and the formative bed-form dimensions are likely to be valid from both recent and older deposits. The data show that estimates of mean bankfull flow depth derived from dune, unit bar, and cross-bar channel deposits are all very similar. Thus, the use of all these metrics together can provide a useful check that all components and scales of the alluvial architecture have been identified correctly when building reservoir models. The data also highlight several practical issues with identifying and applying data relating to cross-strata. For example, the deposits of unit bars were found to be severely truncated in length and width, with only approximately 10% of the mean bar-form length remaining, and thus making identification in section difficult. For similar reasons, the deposits of compound bars were found to be especially difficult to recognize, and hence, estimates of channel depth based on this method may be problematic. Where only core data are available (i.e., no outcrop data exist), formative flow depths are suggested to be best reconstructed using cross-strata formed by dunes. However, theoretical relationships between the distribution of set thicknesses and formative dune height are found to result in slight overestimates of the latter and, hence, mean bankfull flow depths derived from these measurements. This article illustrates that the preservation of fluvial cross-strata and, thus, the paleohydraulic inferences that can be drawn from them, are a function of the ratio of the size and migration rate of bed forms and the time scale of aggradation and channel migration. These factors must thus be considered when deciding on appropriate length:thickness ratios for the purposes of object-based modeling in reservoir characterization.