Numerical Modelling to Evaluate Sedimentation Effects on Heat Flow and Subsidence during Continental Rifting

Sedimentation impacts thermal and subsidence evolution in continental rifting. Estimating the blanketing effect of sediments is crucial to reconstructing the heat flow during rifting. The sedimentary load affects the basin subsidence rate. Numerical investigation of these effects requires active and complex simulations of the thermal structure, lithospheric stretching, and sedimentation. In this paper, we introduce a numerical model to quantify these effects, which was developed using the COMSOL Multiphysics® simulation software. Our numerical setting for the analytical and numerical solutions of thermal structure and subsidence is based on previous continental rifting studies. In our model, we accumulate a column of 5 m thick sediment layers with varied stretching factors and sedimentation rates, spanning the syn-rift to early post-rift phases over a period of 12 myr. Our results provide intuitive models to understand these sedimentation effects. The models show that an increase in sedimentation thickness significantly decreases surface heat flow, leading to lower geothermal temperature, and amplifies the subsidence magnitude. The findings also demonstrate that increases in the stretching factor and sedimentation rate enhance the blanketing effect and subsidence rate. Based on these results, we discuss key outcomes for geological applications and the possible limitations of our approach.

Victoria continental microplate dynamics controlled by the lithospheric strength distribution of the East African Rift

The Victoria microplate between the Eastern and Western Branches of the East African Rift System is one of the largest continental microplates on Earth. In striking contrast to its neighboring plates, Victoria rotates counterclockwise with respect to Nubia. The underlying cause of this distinctive rotation has remained elusive so far. Using 3D numerical models, we investigate the role of pre-existing lithospheric heterogeneities in continental microplate rotation. We find that Victoria’s rotation is primarily controlled by the distribution of rheologically stronger zones that transmit the drag of the major plates to the microplate and of the mechanically weaker mobile belts surrounding Victoria that facilitate rotation. Our models reproduce Victoria’s GPS-derived counterclockwise rotation as well as key complexities of the regional tectonic stress field. These results reconcile competing ideas on the opening of the rift system by highlighting differences in orientation of the far-field divergence, local extension, and the minimum horizontal stress.

One of the largest continental microplates on Earth is situated in the center of the East African Rift System, and oddly, the Victoria microplate rotates counterclockwise with respect to the neighboring African tectonic plate. Here, the authors' modelling results suggest that Victoria microplate rotation is caused by edge-driven lithospheric processes related to the specific geometry of rheologically weak and strong regions.

The importance of structural softening for the evolution and architecture of passive margins

Lithospheric extension can generate passive margins that bound oceans worldwide. Detailed geological and geophysical studies in present and fossil passive margins have highlighted the complexity of their architecture and their multi-stage deformation history. Previous modeling studies have shown the significant impact of coarse mechanical layering of the lithosphere (2 to 4 layer crust and mantle) on passive margin formation. We built upon these studies and design high-resolution (~100–300 m) thermo-mechanical numerical models that incorporate finer mechanical layering (kilometer scale) mimicking tectonically inherited heterogeneities. During lithospheric extension a variety of extensional structures arises naturally due to (1) structural softening caused by necking of mechanically strong layers and (2) the establishment of a network of weak layers across the deforming multi-layered lithosphere. We argue that structural softening in a multi-layered lithosphere is the main cause for the observed multi-stage evolution and architecture of magma-poor passive margins.

Rift migration explains continental margin asymmetry and crustal hyper-extension

When continents break apart, continental crust and lithosphere are thinned until break-up is achieved and an oceanic basin is formed. The most remarkable and least understood structures associated with this process are up to 200 km wide areas of hyper-extended continental crust, which are partitioned between conjugate margins with pronounced asymmetry. Here we show, using high-resolution thermo-mechanical modelling, that hyper-extended crust and margin asymmetry are produced by steady state rift migration. We demonstrate that rift migration is accomplished by sequential, oceanward-younging, upper crustal faults, and is balanced through lower crustal flow. Constraining our model with a new South Atlantic plate reconstruction, we demonstrate that larger extension velocities may account for southward increasing width and asymmetry of these conjugate magma-poor margins. Our model challenges conventional ideas of rifted margin evolution, as it implies that during rift migration large amounts of material are transferred from one side of the rift zone to the other.

Depth-dependent extension, two-stage breakup and cratonic underplating at rifted margins

Uniform lithospheric extension predicts basic properties of non-volcanic rifted margins but fails to explain other important characteristics. Significant discrepancies are observed at 'type I' margins (such as the Iberia-Newfoundland conjugates), where large tracts of continental mantle lithosphere are exposed at the sea floor, and 'type II' margins (such as some ultrawide central South Atlantic margins), where thin continental crust spans wide regions below which continental lower crust and mantle lithosphere have apparently been removed. Neither corresponds to uniform extension. Instead, either crust or mantle lithosphere has been preferentially removed. Using dynamical models, we demonstrate that these margins are opposite end members: in type I, depth-dependent extension results in crustal-necking breakup before mantle-lithosphere breakup and in type II, the converse is true. These two-layer, two-stage breakup behaviours explain the discrepancies and have implications for the styles of the associated sedimentary basins. Laterally flowing lower-mantle cratonic lithosphere may underplate some type II margins, thereby contributing to their anomalous characteristics.

Sequential faulting explains the asymmetry and extension discrepancy of conjugate margins

During early extension, cold continental lithosphere thins and subsides, creating rift basins. If extension continues to final break-up, the split and greatly thinned plates subside deep below sea level to form a conjugate pair of rifted margins. Although basins and margins are ubiquitous structures, the deformation processes leading from moderately extended basins to highly stretched margins are unclear, as studies consistently report that crustal thinning is greater than extension caused by brittle faulting. This extension discrepancy might arise from differential stretching of brittle and ductile crustal layers, but that does not readily explain the typical asymmetric structure of conjugate margins-in cross-section, one margin displays gradual thinning accompanied by large faults, and the conjugate margin displays abrupt thinning but smaller-scale faulting. Whole-crust detachments, active from early in the rifting, could in theory create both thinning and asymmetry, but are mechanically problematical. Furthermore, the extension discrepancy occurs at both conjugate margins, leading to the apparent contradiction that both seem to be upper plates to a detachment fault. Alternative models propose that much brittle extension is undetected because of seismic imaging limitations caused either by subseismic-resolution faulting, invisible deformation along top-basement 100-km-scale detachments or the structural complexity of cross-cutting arrays of faults. Here we use depth-migrated seismic images to accurately measure fault extension and compare it with crustal thinning. The observations are used to create a balanced kinematic model of rifting that resolves the extension discrepancy by producing both fault-controlled crustal thinning which progresses from a rift basin to the asymmetric structure, and extreme thinning of conjugate rifted margins. Contrary to current wisdom, the observations support the idea that thinning is to a first degree explained by simple Andersonian faulting that is unambiguously visible in seismic data.

The formation of non-volcanic rifted margins by the progressive extension of the lithosphere: the example of the West Iberian margin

Non-volcanic margins such as the West Iberian margin exhibit certain characteristics, such as a deficit of synrift igneous rock, a zone of exhumed subcontinental mantle in the continent–ocean transition and an apparent extension discrepancy. These observations can be explained as a consequence of the progressive extension of the lithosphere above relatively cool mantle. The evolving rheological stratification of the lithosphere controls the style of extension at different lithospheric levels at different times; extension is probably heterogeneous at all stages, with lower crustal and upper mantle boudinage controlling the patterns of thinning and mantle upwelling early in the rift history, and complete crustal embrittlement and mantle serpentinization controlling the formation of late-stage detachment faults. Extension in the brittle crust is via multiple phases of faulting, with a general focusing of extension towards the incipient ocean.

The lack of melt is explained by a combination of heterogeneous extension of the lower lithosphere and a cool subcontinental geotherm. The extension discrepancy may in places be controlled by depth-dependent stretching of the crust through lower crustal boudinage, but may also simply be the result of incomplete recognition of the entire polyphase faulting history. The latter seems to be the case for West Iberia.

Evidence for all these processes can be found at the West Iberian rifted margins as well as those preserved and partially exposed in the Alps.

A mechanism to thin the continental lithosphere at magma-poor margins

Where continental plates break apart, slip along multiple normal faults provides the required space for the Earth's crust to thin and subside. After initial rifting, however, the displacement on normal faults observed at the sea floor seems not to match the inferred extension. Here we show that crustal thinning can be accomplished in such extensional environments by a system of conjugate concave downward faults instead of multiple normal faults. Our model predicts that these concave faults accumulate large amounts of extension and form a very thin crust (< 10 km) by exhumation of mid-crustal and mantle material. This transitional crust is capped by sub-horizontal detachment surfaces over distances exceeding 100 km with little visible deformation. Our rift model is based on numerical experiments constrained by geological and geophysical observations from the Alpine Tethys and Iberia/Newfoundland margins. Furthermore, we suggest that the observed transition from broadly distributed and symmetric extension to localized and asymmetric rifting is directly controlled by the existence of a strong gabbroic lower crust. The presence of such lower crustal gabbros is well constrained for the Alpine Tethys system. Initial decoupling of upper crustal deformation from lower crustal and mantle deformation by progressive weakening of the middle crust is an essential requirement to reproduce the observed rift evolution. This is achieved in our models by the formation of weak ductile shear zones.

Crustal flow modes in large hot orogens

Crustal-scale channel flow numerical models support recent interpretations of Himalayan—Tibetan tectonics proposing that gravitationally driven channel flows of low-viscosity, melt-weakened, middle crust can explain both outward growth of the Tibetan Plateau and ductile extrusion of the Greater Himalayan Sequence. We broaden the numerical model investigation to explore three flow modes: homogeneous channel flow (involving laterally homogeneous crust); heterogeneous channel flow (involving laterally heterogeneous lower crust that is expelled and incorporated into the mid-crustal channel flow); and the hot fold nappes style of flow (in which mid-/lower crust is forcibly expelled outward over a lower crustal indentor to create fold nappes that are inserted into the mid-crust). The three flow modes are members of a continuum in which the homogeneous mode is driven by gravitational forces but requires very weak channel material. The hot fold nappe mode is driven tectonically by, for example, collision with a strong crustal indentor and can occur in crust that is subcritical for homogeneous flows. The heterogeneous mode combines tectonic and gravitationally driven flows. Preliminary results also demonstrate the existence and behaviour of mid-crustal channels during advancing and retreating dynamical mantle lithosphere subduction. An orogen temperature—magnitude (T-M) diagram is proposed and the positions of orogens in T-M space that may exhibit the flow modes are described, together with the characteristic positions of a range of other orogen types.