Quantitative stratigraphic analysis in a source-to-sink numerical framework

The sedimentary architecture at continental margins reflects the interplay between the rate of change of accommodation creation (δA) and the rate of change of sediment supply (δS). Stratigraphic interpretation increasingly focuses on understanding the link between deposition patterns and changes in δA=δS, with an attempt to reconstruct the contributing factors. Here, we use the landscape modelling code pyBadlands to (1) investigate the development of stratigraphic sequences in a source-to-sink context; (2) assess the respective performance of two well-established stratigraphic interpretation techniques: the trajectory analysis method and the accommodation succession method; and (3) propose quantitative stratigraphic interpretations based on those two techniques. In contrast to most stratigraphic forward models (SFMs), pyBadlands provides self-consistent sediment supply to basin margins as it simulates erosion, sediment transport and deposition in a source-to-sink context. We present a generic case of landscape evolution that takes into account periodic sea level variations and passive margin thermal subsidence over 30 million years, under uniform rainfall. A set of post-processing tools are provided to analyse the predicted stratigraphic architecture. We first reconstruct the temporal evolution of the depositional cycles and identify key stratigraphic surfaces based on observations of stratal geometries and facies relationships, which we use for comparison to stratigraphic interpretations. We then apply both the trajectory analysis and the accommodation succession methods to manually map key stratigraphic surfaces and define sequence units on the final model output. Finally, we calculate shoreline and shelf-edge trajectories, the temporal evolution of changes in relative sea level (proxy for δA) and sedimentation rate (proxy for δS) at the shoreline, and automatically produce stratigraphic interpretations. Our results suggest that the analysis of the presented model is more robust with the accommodation succession method than with the trajectory analysis method. Stratigraphic analysis based on manually extracted shoreline and shelf-edge trajectory requires calibrations of time-dependent processes such as thermal subsidence or additional constraints from stratal terminations to obtain reliable interpretations. The 3-D stratigraphic analysis of the presented model reveals small lateral variations of sequence formations. Our work provides an efficient and flexible quantitative sequence stratigraphic framework to evaluate the main drivers (climate, sea level and tectonics) controlling sedimentary architectures and investigate their respective roles in sedimentary basin development.