A simulation of continental basin margin sedimentation in response to crustal movements, eustatic sea level change, and sediment accumulation rates

As eustasy, subsidence, and sediment accumulation vary, a 2D computer-based graphical simulation generates on-lapping and off-lapping geometries of both marine and near coastal alluvial deposits, reproducing timelines within sediment-bodies at basin margins. In the simulation, deposition is expressed by creation of new surfaces above previous ones. Thicknesses of layers are reduced by both erosion and compaction while their surfaces move vertically in response to tectonic change and loading. Simulation is divided into a series of equal time steps in which sediment is deposited as an array of en-echelon columns that mark the top of the previous depositional surface. The volume of sediment deposited in each time step is expressed as a 2D cross section and is derived from two right-angle triangles (sand and shale), whose areas are a 2D expression of the quantity of sediment deposited at that time step and whose length matches the width of the offshore sediment wedge seaward of the shoreline. Each column in the array is filled by both marine sediments up to sea level, and alluvial sediments to a surface determined by an alluvial angle that is projected landward from the shore to its intersection with the previous surface. Each time the area representing the sediment column is subtracted from the triangles, the triangle heights are reduced correspondingly. This process is repeated until the triangle heights match the position of sea level above the sediment surface, at which time the remaining area of the sediment triangle is deposited seaward as a single wedge of offshore sediments. This simulation is designed to aid interpretation of stratigraphic sequences. It can be used as a complement to seismic stratigraphy or can be used alone as an inexpensive test of stratigraphic models.