16The Paleozoic intracratonic North African Platform is characterized by an association of 17 arches (ridges, domes, swells or paleo-highs) and low subsidence rate syncline basins of

Abstract. The Paleozoic intracratonic North African Platform is characterized by an association of arches (ridges, domes, swells, or paleo-highs) and low subsidence rate syncline basins of different wavelengths (75–620 km). The Reggane, Ahnet, Mouydir and Illizi basins are successively delimited from east to west by the Amguid El Biod, Arak-Foum Belrem, and Azzel Matti arches. Through the analysis of new unpublished geological data (i.e., satellite images, well logs, seismic lines), the deposits associated with these arches and syncline basins exhibit thickness variations and facies changes ranging from continental to marine environments. The arches are characterized by thin amalgamated deposits with condensed and erosional surfaces, whereas the syncline basins exhibit thicker and well-preserved successions. In addition, the vertical facies succession evolves from thin Silurian to Givetian deposits into thick Upper Devonian sediments. Synsedimentary structures and major unconformities are related to several tectonic events such as the Cambrian–Ordovician extension, the Ordovician–Silurian glacial rebound, the Silurian–Devonian Caledonian extension/compression, the late Devonian extension/compression, and the Hercynian compression. Locally, deformation is characterized by near-vertical planar normal faults responsible for horst and graben structuring associated with folding during the Cambrian–Ordovician–Silurian period. These structures may have been inverted or reactivated during the Devonian (i.e., Caledonian, Mid–Late Devonian) compression and the Carboniferous (i.e., pre-Hercynian to Hercynian). Additionally, basement characterization from geological and geophysics data (aeromagnetic and gravity maps), shows an interesting age-dependent zonation of the terranes which are bounded by mega-shear zones within the arches–basins framework. The old terranes are situated under arches while the young terranes are located under the basins depocenter. This structural framework results from the accretion of Archean and Proterozoic terranes inherited from former orogeny (e.g., Pan-African orogeny 900–520 Ma). Therefore, the sedimentary infilling pattern and the nature of deformation result from the repeated slow Paleozoic reactivation of Precambrian terranes bounded by subvertical lithospheric fault systems. Alternating periods of tectonic quiescence and low-rate subsidence acceleration associated with extension and local inversion tectonics correspond to a succession of Paleozoic geodynamic events (i.e., far-field orogenic belt, glaciation).

Intracratonic basins tend to subside much longer than the timescale predicted by thermal relaxation of the lithosphere. Many hypotheses have been suggested to explain their longevity, yet few have been tested using quantitative thermo-mechanical numerical models, which capture the dynamic of the lithosphere. Lithospheric scale geodynamic modelling preserving the tectono-stratigraphic architecture of these basins is challenging because they display only few kilometres of subsidence over 1000 of km during time periods exceeding 250 Myr. Here we present simulations that are designed to examine the relative role of thermal anomaly, tectonics and heterogeneity of the lithosphere on the dynamics of intracratonic basins. Our results demonstrate that initial heterogeneity of accretionary continental lithosphere explains long-term subsidence and the arches-basins architecture of Saharan type intracratonic basins at first order. The simulations show that initially warm and heterogeneous lithospheres inherited from accretion are strong enough to resist local isostatic re-equilibration for very long period of time. Indeed, the lateral density variations store potential gravitational energy that is then slowly dissipated by differential erosion and slow vertical movements. For relatively well-accepted coefficient of erosion of 10-6 m2/s, the subsidence last longer than 250 Myr. Extensional tectonic forcing and thermal anomalies both results in an effective strength drop of the lithosphere, which allows a temporarily acceleration of local isostatic re-equilibration. Periodic changes in far field tectonic forcing from extension to compression complicate the tectono-stratigraphic architecture (intra-basin arches, sub-basins) introducing stratigraphic unconformities between different neighbouring basins such as the ones observed in North Africa.

<p>For a long time, the complexity of the lithosphere was ignored by numerical modelling because the inherited structural and compositional complexity of the “real” lithosphere is indeed mainly unknown to geologists so modeler preferred to understand first order parameters such as rate of extension, lithospheric thickness, mechanical coupling or decoupling at the Moho. These models were not representative of any particular region but they were helpful. As a wider community of geologist became interested in numerical modelling, a growing number of numerical models have attempted to account for a major player in structural geology: inheritance. However, the complexity of “real” Earth has been simplified and “idealized” where inherited “anomalies” (e.g. fault, pluton, craton) or a combination of them has been added without really knowing the exact initial conditions which are the unknown of the problem. Yet another approach has been to add a lot of them in a more or less random mater or to replace them by initial noise in the parameters. None of these approaches actually fulfil the need for end-users community to have predictive models.</p><p>Realizing that structural inheritance is some kind of kinematic forcing in the solution of the models but also that it is not possible to anticipate and identify all the geological structures that can be inherited in rifted margin lithospheres, we have developed a new approach, through the integration of a new kinematic module to pTatin2D thermomechanical code, permitting to understand the kinematics of deformation of the continental lithosphere and asthenosphere through time leading to the establishment of rifted margins. The method is settled and validated by fitting the architecture (i.e. basement, Moho, LAB, Tmax) and by solving the kinematics of a random unknow 2D cross section extracted from 3D thermomechanical rifted margin model.</p><p>This new tool aims to help geologists to better constrain and draw on their 2D geological cross sections the position of the Moho, the Lithosphere-Asthenosphere boundary (LAB), the temperature isotherms and the heat flux.</p><p>Key words: Kinematic thermomechanical modelling, asymmetric rifted margin architecture, modelling method.</p>