*4.3. Past Geometries*

The restoration of past geometries for all models was performed with a backstripping method [83,84], implemented in the TemisFlowTM software (2016 version, Beicip-Franlab, Rueil-Malmaison, France) [85]. The uncompacting used a different Porosity vs. Depth table for each lithological type. Shalier facies have larger deposit porosities but compact more quickly with burial than the sandier ones (plots available in Figure A1 in the Appendix A).

The model bottom was used as a reference horizon. Its shape was arbitrarily input for each deposition event in order to model a progressive subsidence and tilting for the margin basement, following Table 2. The paleo-bathymetry profiles were then computed by the algorithm from the compaction tables and the basement shape. Lastly, absolute sea level variations were input by translating the full model upward or downward from the 0 mark on the vertical axis. These variations were set in accordance with the sedimentary track being deposited, and keeping the full model immersed during the whole simulation.
