Plan
Comptes Rendus

Article original - Tectonique, Tectonophysique
Understanding the role of décollement thickness on the evolution of décollement folds: insights from discrete element models
Comptes Rendus. Géoscience, Volume 354 (2022), pp. 75-91.

Résumé

We explored the influence of a variety of geometric and mechanical factors in the evolution of décollement folds above a weak décollement level in a sedimentary sequence. Under an applied overburden pressure, we construct nine discrete element models (DE-models) to test the effects of mechanical stratigraphy, the thickness of the décollement layer, and the number of strong and weak layers within the sedimentary cover sequence. The effects of the aforementioned parameters on the structural style of fold belts and the mechanical evolution of décollement folds that develop during a progressive deformation are what we are searching for. The one-way lateral motion of a rigid, vertical boundary wall was used to shorten the particle assemblage. The models presented in this study were then compared to a natural structure in Iran and demonstrated that the combined effects of mechanical behaviour, the thickness of décollement layer, and the number of strong and weak layers within the sedimentary cover sequence are playing the decisive role in the structural style, kinematic and mechanical evolution of décollement folds.

Métadonnées
Reçu le :
Révisé le :
Accepté le :
Publié le :
DOI : 10.5802/crgeos.103
Mots clés : Structural style, Kinematics, Mechanical evolution, Fault-related folding, Décollement folding, Discrete element modelling

Nazanin Eslamirezaei 1 ; Seyed Ahmad Alavi 1 ; Seyed Tohid Nabavi 1 ; Mohammad R. Ghassemi 2

1 Faculty of Earth Sciences, Department of Geology, Shahid Beheshti University, Tehran, Iran
2 Research Institute for Earth Sciences, Geological Survey of Iran, Tehran, Iran
Licence : CC-BY 4.0
Droits d'auteur : Les auteurs conservent leurs droits
@article{CRGEOS_2022__354_G1_75_0,
     author = {Nazanin Eslamirezaei and Seyed Ahmad Alavi and Seyed Tohid Nabavi and Mohammad R. Ghassemi},
     title = {Understanding the role of d\'ecollement thickness on the evolution of d\'ecollement folds: insights from discrete element models},
     journal = {Comptes Rendus. G\'eoscience},
     pages = {75--91},
     publisher = {Acad\'emie des sciences, Paris},
     volume = {354},
     year = {2022},
     doi = {10.5802/crgeos.103},
     language = {en},
}
TY  - JOUR
AU  - Nazanin Eslamirezaei
AU  - Seyed Ahmad Alavi
AU  - Seyed Tohid Nabavi
AU  - Mohammad R. Ghassemi
TI  - Understanding the role of décollement thickness on the evolution of décollement folds: insights from discrete element models
JO  - Comptes Rendus. Géoscience
PY  - 2022
SP  - 75
EP  - 91
VL  - 354
PB  - Académie des sciences, Paris
DO  - 10.5802/crgeos.103
LA  - en
ID  - CRGEOS_2022__354_G1_75_0
ER  - 
%0 Journal Article
%A Nazanin Eslamirezaei
%A Seyed Ahmad Alavi
%A Seyed Tohid Nabavi
%A Mohammad R. Ghassemi
%T Understanding the role of décollement thickness on the evolution of décollement folds: insights from discrete element models
%J Comptes Rendus. Géoscience
%D 2022
%P 75-91
%V 354
%I Académie des sciences, Paris
%R 10.5802/crgeos.103
%G en
%F CRGEOS_2022__354_G1_75_0
Nazanin Eslamirezaei; Seyed Ahmad Alavi; Seyed Tohid Nabavi; Mohammad R. Ghassemi. Understanding the role of décollement thickness on the evolution of décollement folds: insights from discrete element models. Comptes Rendus. Géoscience, Volume 354 (2022), pp. 75-91. doi : 10.5802/crgeos.103. https://comptes-rendus.academie-sciences.fr/geoscience/articles/10.5802/crgeos.103/

Version originale du texte intégral (Proposez une traduction )

1. Introduction

The interaction of folding, thrust faulting, and mechanical behaviour in contractional settings, particularly fold-and-thrust belts, has been recognized in deformed sequence so that assessment of these relationships is fundamental for understanding deformation processes, kinematic and mechanical completion of the structures, and structural style [Brandes and Tanner 2014, Derikvand et al. (2018, Motamedi and Gharabeigli 201]. There are several assumptions to consider when analysing thrust belts, where one of the most important assumptions is that internal layer deformation occurs via simple shear parallel to bedding [e.g., Rowan et al. 2019]. As a result, a wide range of studies have been conducted to investigate the determinants based on these assumptions [Ridley and Casey 1989, Marshak et al. 1992, Corredor et al. 2005]. Many fold-and-thrust belts, for example, are the result of orthogonal or oblique convergence–subduction and/or collision. There are studies being conducted to investigate the influencing factors in the evaluation of thrust belts, such as those formed by convergence and the influence of convergence rate on structural style, which is a critical factor in the formation of fold-and-thrust belts [e.g., Bonini 2003, McClay et al. 2004, Farzipour-Saein and Koyi 2014, Schori et al. 2015, Ghanadian et al. 2017, Chen et al. 2018, Cui et al. 2020, Pla et al. 2019, Gu et al. 2021].

Structural terminologies of fold-and-thrust structures in contractional settings have been investigated using kinematic models of fault-related folds including detachment (or décollement), fault-bend, fault-propagation, and trishear folds [e.g., Brandes and Tanner 2014, Cosgrove 2015, Butler et al. 2018]. Moreover, fold-and-thrust belts often possess rich hydrocarbon reserves, leading to do wide attention to their structures and evolution [e.g., Kendall et al. 2019, Hammerstein et al. 2020]. Therefore, the evaluation of fault-related folds influences the development and exploration of a hydrocarbon-bearing basin by affecting the trap, seal capacity, reservoir quality, migration and also drilling costs and hazards. According to many studies on fault-related folds, folds geometry generally is controlled by fault geometry and fault slip/fault-propagation ratio in fault-propagation folds, and fault displacement and décollement layer thickness in detachment folds [e.g., Jamison 1987, Mitra 1990, Allmendinger 1998, Suppe et al. 2004, Brandes and Tanner 2014].

Among these models, décollement folds are considered as a prevalent structural style in fold-and-thrust belts and in most cases thin-skinned (décollement-dominated) tectonics. Décollement fold structures have formed from the shortening of a competent rock layers upper a décollement layer played as a lower competent layer such as over-pressured shales, salts and evaporates. These folds are not related to ramps and the propagating thrust segments is sub-horizontal. Their geometries range from concentric, chevron, to box-shape fold in geometry [Dean et al. 2013, Meng and Hodgetts 2019b]. Detachment folds can have a variety of structural styles, like: symmetric, asymmetric, disharmonic, lift-off, multi-detachment to faulted [e.g., Mitra and Namson 1989, Poblet and McClay 1996, Mitra 2003, Wallace and Homza 2004, Suppe 2011, Morley et al. 2017, Li and Mitra 2020]. Mitra [2003] demonstrates buckling of the above décollement levels where mobile material is placed beneath synclines and flows into anticline hinges. As contraction and distributed simple shear increases, décollement folds, in addition to amplification, can be asymmetric or overturned due to limb rotation and hinge migration. Faulting in their final phase of deformation is also common to form fault-bend and or fault-propagation folds [e.g., Mitra 2003, Morley et al. 2017, Nabavi et al. 2019].

Décollement folds are abundant in fold-and-thrust belts, for example, the Appalachian fold-and-thrust belt [e.g., Lammie et al. 2020], the Jura Mountains [e.g., Schori et al. 2015], the Zagros fold and thrust [e.g., Sherkati et al. 2005, 2006, Vergés et al. 2011], the Mexican Ridges fold belt [e.g., Yarbuh et al. 2018], and the deep-water Niger Delta [e.g., Briggs et al. 2006, Morley et al. 2011]. Studies on décollement folds have been successfully conducted using different kinds of approaches such as field data [e.g., Suppe 2011, Morley et al. 2017], subsurface seismic reflection data [e.g., Derikvand et al. (2018, Morley et al. 2017, Li and Mitra 2020], analogue [e.g., Farzipour-Saein 2017, Ghanadian et al. 2017, Li and Mitra 2017] and various numerical techniques such as FEM or DEM. For example, studies conducted by the discrete element technique [e.g., Schöpfer et al. 2006, 2017, Hardy and Finch 2005, Meng and Hodgetts 2019b].

Proper interpretation of the structural evolution and structural style are essential to build the accurate understanding of the surface and subsurface geometry of structures. In this way, numerical modelling using the discrete element method (DEM), provides means for simulating and tracking the spatial and temporal evolution of stress distribution, characteristic and sequential evolution of structures, strain localization and the prediction of onset of associated geological structures at different scales, for example, the studies of the development of normal and thrust faults [Saltzer and Pollard 1992, Donzé et al. (1994, Abe et al. 2011, Schöpfer et al. 2006, 2016, 2017, Hardy 2011, 2013, Smart et al. 2011, Deng et al. 2017, Finch and Gawthorpe 2017], relay structures [Imber et al. 2004], décollement and fault-related folds [Finch et al. 2003, Cardozo et al. 2005, Hardy and Finch 2005, 2006, 2007, Benesh et al. 2007, Hughes et al. 2014], accretionary wedges and FAT belts [Burbidge and Braun 2002, Naylor et al. 2005, Yamada et al. 2006, Hardy et al. 2009, Dean et al. 2013, Morgan 2015, Morgan and Bangs 2017, Meng and Hodgetts 2019a,b]. In this regard, we utilize results from a series of numerical DE-models using PFC3D® software package that can be followed during deformation process, for example, décollement slip, onset and amplification of décollement folds, and the initiation and propagation of discrete thrust faults.

This study seeks to address the importance and influence of (1) the thickness of décollement layer, and (2) the number of strong and weak layers within the sedimentary cover sequence under an applied horizontal pressure, on the structural and mechanical formation of décollement folds in a contractional setting. It should be noted that the presented models are simplistic. Geological processes such as diagenesis, erosion, fluid flow, temperature, etc., are not considered. The primary reason for this approach is that we seek to isolate the effects of a few selected parameters in order to keep the analysis traceable [Nabavi et al. 2020]. The outcome of this study is compared to a natural example from the central Iran structural zone as well as the Zagros fold-and-thrust belt and large-scale fold-and-thrust belt in South China.

2. The discrete element modelling method (DEM)

2.1. Basic principles

In this study, a series of DE-mechanical models and a particle-based numerical method is applied for simulation of deformation of rock layers during shortening situation. This approach uses a finite difference approach to solve Newton’s equations of motion. The rock mass is considered as an assemblage of soft circular, spheres, or disc, which interacts with each other in pairs as if connected by breakable elastic bonds (springs) and undergo motion relative to one other. These particles interact through elastic, gravitational, frictional, and viscous forces. At each discrete time step, the particles are displaced to their new positions within the model by integrating their equations of motion using Newtonian laws and a velocity-Verlet-based numerical scheme.1 The particle bonding is controlled by normal and shear bond stiffness, bond strength, and bond radius. In this regard, the normal and shear contact forces (Mohr–Coulomb criterion) are calculated between two overlapping or contacting particles using a rheological contact law. Bonds are broken when the inter-particle forces exceed the bond strength. In this regard, slip that is resisted by a frictional strength can occur between particles that are not bonded. Thus, discontinuities and heterogeneities form and evolve in response to changing stress condition and material properties [e.g., Cundall and Strack 1979, Hardy et al. 2009, Hughes et al. 2014, Morgan 2015, Smart et al. 2011, Hardy 2019, Meng and Hodgetts 2019b]. Many previous investigations have used the DE-models and have successfully produced realistic styles of geological and tectonic features such as normal faulting in layered sequences [e.g., Botter et al. 2014, Hardy 2019], fault-bend folding [e.g., Benesh et al. 2007, Hughes et al. 2014], fault-propagation folding [e.g., Hardy and Finch 2006, Hughes and Shaw 2015], and fold-and-thrust belts [e.g., Hardy et al. 2009, Morgan 2015, Meng and Hodgetts 2019b]. Hence, DEM, is an appropriate approach to study the initiation, propagation, interaction and linkage, and the evolution of folding, fracturing, and faulting in a realistic manner involving large deformation and large relative displacements.

2.2. Model setup

The mechanical modelling method we utilize is DEM that was developed using the Itasca Particle Flow Code (PFC3D®, version 6,13, Itasca Consulting Group) 3D software package. PFC® was initially developed for describing the mechanical behaviour of material (Soil or rocks), allows flexibility in simulations and is widely used among geologists as it permits modelling of large deformations.

Figure 1.

The initial design and boundary conditions for discrete element models. The models were made up of 21,403 elements with radii ranging from 25 to 41 m for the décollement and cover layers, the red lower décollement, and a series of weak and strong layers in a hexagonal box (10 km length, 1 km width, 0.5 km thickness). The boundary condition is considered free in surface. The left sidewall is moving towards right side. The fs  = shear stress, fn  = normal stress, and b is particle bond.

Table 1.

Summary of material properties and parameters used in DE-models

ParametersDécollement layerOverburden
Weak layersStrong layers
Particle radius25–41 m
Contact Young’s modulus0.1 × 1090.5 × 1091 × 109
Contact friction 00.5
Tensile bond strength 02.5 × 1065 × 106
Average density3000 kg/m3

We created a regular hexagonal model (Figure 1) consisting of 21,403 spherical particles in the part of upper crust model with dimensions of 10 km length × 1 km width × 0.5 km thickness (x,y,z). The particle radii range from 25 to 41 m both for the décollement and cover layers that follow a Gaussian distribution for particle size. The mechanically layered cover sequence consists of alternating relatively “strong” layer and relatively “weak” layers. The average rock density (𝜌) for all particles is about 3000 kg/m3. Particle stiffness values of 1 × 109, 0.5 × 109, and 0.1 × 109 N/m, for strong, weak, and décollement layer, respectively, were assigned for normal (kn) and shear (ks) stiffness (Table 1). The coefficient of sliding friction between particles (𝜇) was assigned as 0.5, whereas the décollement layer is frictionless to ensure its low strength [e.g., Morgan 2015, Meng and Hodgetts 2019b].

Table 2.

Presented various initial configurations of the models

Model No. 1 2 3 4 5 6 7 8 9
Td0.10.10.10.20.20.20.30.30.3
N543543543

Each model has unique initial setting as (Td, N) that Td is the thickness of basal décollement layer as a fraction of total thickness (0.1, 0.2, 0.3 km of total thickness or 100 m, 200 m, 300 m) and N is the number of weak layers.

We examined nine models (Table 2) with varying décollement layer thickness as a fraction of total initial thickness (0.1, 0.2, 0.3 km of total thickness or 100 m, 200 m, 300 m) and the number of strong and weak layers within the sedimentary cover (5 weak layers in models 1, 4, 7; 4 weak layers in models 2, 5, 7; and 3 weak layers in models 3, 6, 9). Coloured layers within the cover sequence above the décollement layer indicate contrast in material properties.

Our experimental setup consists of a horizontal surface and two vertical boundaries. All surfaces are defined in the model by rigid, frictional boundary walls (i.e., free-slip boundary conditions). Numerically, free surface boundary conditions are applied on the top model, thus permitting self-consistent evolution of topography [Nabavi et al. 2020]. Acceleration due to gravity (9.81 m/s2) is imposed throughout the deformation as a body force. In all models, the horizontal, vertical, compressional load was imposed horizontally [i.e., displacement-based boundary conditions; Nabavi et al. 2018, 2020]. All results are shown here after 10%, 20%, and 30% shortening, which is equal to 1000 m, 2000 m, after 10%, 20%, and 30% shortening, which is equal to 1000 m, 2000 m, and 3000 m displacement on the décollement and overburden layers. The shear strength 5 × 106 N/m and 2.5 × 106 N/m for strong and weak layers, respectively.

3. Results

We present the main features observed in a series of DE-models after 30% (3000 m displacement) horizontal shortening to investigate décollement folding above a flat décollement layer. There are large number of possible model setups that we could present, however, we discussed the effects of décollement layer thickness and the number of strong and weak layers within the sedimentary cover on the mechanical evolution and structural style of box-shaped décollement folds in a contractional environment through nine different models which is based on progressive deformation events. In all models, the sedimentary cover sequence overly a décollement horizon at the base of models.

3.3. Models with a thin décollement

Figure 2.

Evolution of DEM to show the folding in the models 1, 2, and 3 with thin décollement after 30% (3000 m displacement) orthogonal shortening in three steps of shortening 10%, 20%, and 30%. The models are the same in terms of the thickness of the décollement (Td = 0.1) and different in the sedimentary sequences (N). The value of N is equal to 5, 4, and 3 weak layers, respectively, for models 1, 2, and 3. The final geometries are denoted by a 30% shortening. The uplift is about 1.5 km for model 1, 1.65 km for model 2, and 1.5 km for model 3.

The structures observed in a set of numerical models exhibit a wide range of structural style that varies as a result and the function of model design, boundary conditions, and the material being deformed. Model 1 (Figure 2a) presents how a thin décollement layer affects deformation, the structural style and evolution of décollement folding, and its spatial and temporal geometric evolution is shown in three models and displacement steps at intervals of 100 m of shortening. In this model, the sedimentary cover layers consisted of a series of strong bonded layers (six layers thick) intercalated with weak bonded layers (five layers thick). The décollement layer in this model is 0.1 km thick. After the primary phase of homogeneous shortening 10%, the deformation shows pure-shear thickening (contraction or buckling) with the low-amplitude asymmetric décollement fold at the onset of deformation. Increasing deformation (20%) causes three major asymmetric anticlines at different locations of the model above the décollement layer including F1 in the hinterland (left sidewall), F2 at the model centre, and F3 in the foreland (right sidewall) (Model 1 in Figure 2). Anticlines have higher amplitude for those bounded by the lateral boundaries of the model. With increasing horizontal shortening and displacement, both fold limbs rotate about 10° and lengthen through time. After 30% shortening, the middle anticline gradually becomes symmetric and amplifying due to strain accumulation in the décollement layer (i.e., flexural slip) and localized strain zones. Finally, developing to the nearly symmetric sinusoidal, bulb-like to box-shaped décollement fold with a curved hinge. The curved hinge is bounded by a set of conjugate high-strain shear zones, that are closed to the strong layers, whereas its wavelength did not show significant changes through time. Figure 2(a) shows the final geometry that demonstrates most of the strain after 30% shortening and that propagates the fold structure upwards. At the end of shortening, the basal layers were folded isoclinally in the core of the structure, with limbs dipping sub-vertically in the lower layers. Upper layers did not dip as steeply and were extended and thinned in crestal areas as a result of buckling of the uppermost hanging-wall strata, whereas they were folded and thickened in adjacent synclinal areas. The crest of folds widened outwards and upwards from the core of them, which generally the resulting structure was upright and symmetrical box-shaped décollement fold. Furthermore, Figure 2(a) indicates flexural flow developed in the core of the box-shaped décollement folds, so that it is filled with weak material that evacuated beneath synclines and flowed into the anticline core. Some minor faults and shear bands develop in limbs of the box-shaped décollement fold to accommodate contraction as penetrative strain [e.g., Burberry 2015, Nabavi et al. 2019] that does not propagate across the entire multilayer but rather links the weak layers with the layered rock sequence.

Model 2 (Figure 2b) has the same décollement thickness of 0.1 km, but the number of weak layers inside the sedimentary cover is reduced to four, and the model is dominated by strong layers. In 10% of shortening, there are deformation in the surface that shaped the topography. Increasing deformation (20%) causes two major asymmetric anticlines at two parts of the model above the décollement layer including F1 in the hinterland (left sidewall) and F2 in the foreland (right sidewall). Anticlines have constant wavelength and amplitude for those bounded by the lateral boundaries of the model. As the deformation progresses (30% horizontal shortening), F3, an upright and asymmetric box-shaped décollement fold and low-amplitude anticline developing and amplifying at the middle part of the model due to continuous contraction, strain accumulation in the décollement layer (i.e., flexural slip) and localized strain zones. This low-amplitude anticline finally may bound by a set of conjugate faults or high-strain shear zones, which are more restricted to the strong layers to accommodate contraction as penetrative strain, and developed an asymmetric faulted décollement fold. Like Model 1, fold limbs both rotated and lengthened though time due to distributed strain throughout the fold limbs. As the limbs steepen, slip becomes more localized on flexural slip surfaces, then finally strain localizes at the synclinal axial zones where more material is progressively incorporated into the structure. At the end of the shortening, the basal layers were folded isoclinally in the core of the structure and limb dips were 80–90° in the lower layers. Upper layers did not dip as steeply and were extended and thinned in crestal areas as a result of buckling of the uppermost hanging-wall strata, whereas they were folded and thickened in adjacent synclinal areas. In contrast to model 1, the folds formed in model 2 have much higher amplitudes. Similar to Model 1, Model 2 (Figure 2b) shows the core of décollement folds filled with weak material that evacuated beneath synclines and flowed into the anticline core.

Model 3 (Figure 2c) has the same décollement thickness of 0.1 as Models 1 and 2, but there are three weak layers and the model is dominated by strong layers. The whole sedimentary cover sequence, as previously indicated, is deposited over a basal weak décollement level. Figure 2c depicts the evolution of this model. Model 3 developed two main folds (F1 and F2 in the middle part of the model) and two minor folds at two parts of the model above the décollement layer including F3 in the hinterland (left sidewall) and F4 in the foreland (right sidewall) to accommodate shortening. However, in addition to the symmetric box-shape décollement fold and differently from Models 1 and 2, fault-propagation fold takes over décollement fold as one of the main fold type in this model after 20–30% horizontal shortening (Figure 2c). The faults accommodated shortening by thrusting (both fore thrust and back thrust). The thrusts and accumulated fault displacement gradually and dip angles approximately keeping unchanged. However, after this step, the development of the structure changed, fault-propagation folding ceased, and instead the structure tightened and developed a classically upright and symmetric box-shaped décollement fold with continued displacement.

3.2. Models with an intermediate-thick décollement

Models 4, 5, and 6 show how an intermediate-thick décollement with a 0.2 km thickness affects deformation and structural styles. In these models, the sedimentary cover sequence consisted of a series of strong bonded layers (six, five, and four strong layers thick, respectively, for Model 4, Model 5, and Model 6) intercalated with weak bonded layers (six, five and three weak layers thick, respectively, for Model 4, Model 5 and Model 6). In Model 4 (Figure 3a), first, deformation is essentially pure-shear thickening (contraction or buckling) with the low-amplitude asymmetric décollement fold at the onset of deformation. Increasing deformation through time causes two major nearly symmetric anticlines F1 and F2 as sinusoidal, bulb-like to box-shaped geometries at different locations of the model above the décollement layer including.

Figure 3.

Evolution of DEM to show the folding in the models 4, 5, and 6 with an intermediate-thick décollement after 30% (3000 m displacement) orthogonal shortening in three steps of shortening 10%, 20%, and 30%. The models are the same in terms of the thickness of the décollement (Td = 0.2) and different in the sedimentary sequences (N). The value of N is equal to 5, 4, and 3 weak layers, respectively, for models 1, 2, and 3. The final geometries are denoted by a 30% shortening. The uplift is about 1.59 km for model 1, 1.75 km for model 2, and 1.64 km for model 3.

In Model 5 (Figure 3b), the deformation that is pure-shear thickening (contraction or buckling) started from 10% of shortening by changes in the surface with the low-amplitude asymmetric fold. Increasing deformation through time (deformation 20 to 30%) accommodated by a combination of the asymmetric box-shaped décollement fold F1 and the foreland fault-propagation fold F2, so that flexural flow with a distinct fault segment developed in the lower part and core of the fault-propagation fold and it is filled with weak material that evacuated beneath synclines and flowed into the anticline core (Figure 3b in 30%). High-strain shear zones or bands are more restricted to the strong layers of the left-hand limb of the structure. Upper layers dip as moderately to steeply (25° to 45°) as faulted décollement folds and fault-related folds develop and were extended and thinned in crestal areas with possible extensional faulting and collapse structures as a result of buckling of the uppermost hanging-wall strata, whereas they were folded and thickened in adjacent synclinal areas.

The final model in this series (Model 6 or c in Figure 3) responses to 10% shortening by deformation from left side, also at the model centre (20% in Model 6 in Figure 3c), and in the foreland (right sidewall). Anticlines are more cylindrical and have higher amplitude for those bounded by the lateral boundaries of the model. With increasing horizontal shortening and displacement, both fold limbs rotate and lengthen through time. As the deformation progresses and after 30% shortening, the middle anticline F1, gradually became asymmetric and amplifying due to strain accumulation in the décollement layer as both flexural slip and flexural flow, respectively in strong and weak layer. In addition, localized strain zones are developed in different domains of Model 6, where its wavelength increases through time but in a rather irregularly spaced with a steep to overturned left-hand limb and a shallower and longer right-hand limb (i.e., highly asymmetric décollement folds), and the contractional fault-propagation fold so that flexural flow with a distinct fault segment developed in the lower part and core of this fold and it is filled with weak material that evacuated beneath synclines and flowed into the anticline core. Some minor faults and shear bands develop in limbs of décollement folds to accommodate contraction as penetrative strain that does not propagate across the entire multilayer but rather links the weak layers with the layered rock sequence.

3.3. Models with a thick décollement

Figure 4.

Evolution of DEM to show the folding in the models 7, 8, and 9 with a thick décollement after 30% (3000 m displacement) orthogonal shortening in three steps of shortening 10%, 20%, and 30%. The models are the same in terms of the thickness of the décollement (Td = 0.3) and different in the sedimentary sequences (N). The value of N is equal to 5, 4, and 3 weak layers, respectively, for models 1, 2, and 3. The final geometries are denoted by a 30% shortening. The uplift is about 1.54 km for model 1, 1.62 km for model 2, and 1.62 km for model 3.

Model 7 (Figure 4a) depicts how a thick décollement layer influences deformation, the structural style, and evolution of décollement folding. The sedimentary cover sequence in these models was made up of a succession of strongly bonded layers (six strong layers deep) separated by weak layers (five weak layers thick). In this model, the décollement layer is 0.3 km thick. Overall, following the first stage of homogeneous shortening, the deformation is essentially pure-shear thickening (contraction or buckling) with the foreland-verging low-amplitude asymmetric décollement fold at the onset of 10% deformation. Increasing deformation through time (deformation 20%) causes three major nearly symmetric anticlines as sinusoidal, bulb-like to box-shaped geometries at the different location of the model above the décollement layer including F1 in the hinterland (left sidewall), F2 at the model centre, and F3 in the foreland (right sidewall). With increasing horizontal shortening and displacement, both fold limbs rotate and lengthen through time. As the deformation progresses and after 30% shortening, the middle anticline gradually became highly asymmetric and amplifying due to strain accumulation in the décollement layer as both flexural slip and flexural flow, respectively in strong and weak layers, with localized strain zones and a distinct fault segment developed in the lower part and core of this fold, which finally develops to the fault-propagation fold, whereas two box-like décollement folds in hinterland and foreland domains do not show any significant changes in their geometries through time. In the latter case, thrusts would preferentially develop to accommodate applied horizontal shortening rather than buckling. High-strain shear zones or bands are more restricted to the strong layers of the left-hand limb of the structure. Upper layers dip as moderately to steeply as fault-related folds develop and were extended and thinned in crestal areas with possible extensional faulting and collapse structures as a result of buckling of the uppermost hanging-wall strata, whereas they were folded and thickened in adjacent synclinal areas.

Model 8 (Figure 4b) illustrates how a thick décollement layer determines the structural style and evolution of a sedimentary overlay succession composed of a sequence of strong bonded layers (five strong layers thick) intercalated with weak bonded layers (four weak layers thick). In this model, the décollement layer is 0.3 km thick. The deformation is mainly pure-shear thickening (contraction or buckling) with the foreland-verging low-amplitude asymmetric décollement fold at the initiation of deformation, after an early phase of homogeneous shortening of the model. With the passing of time, the degree of deformation increases (deformation 20%) causes the foreland-verging fault-propagation fold with localized strain zones and a distinct fault segment developed in the lower part and core of this fold in the hinterland (left sidewall) and the low-amplitude asymmetric décollement fold at the middle domain of the model. In this case, thrusts would preferentially develop to accommodate applied horizontal shortening rather than buckling. As the deformation progresses and after 30% shortening, the F1 fault-propagation fold further develops and its asymmetry increases, and also F2 the middle anticline gradually became more symmetric and amplifying due to strain accumulation in the décollement layer as a box-like décollement fold with the low-amplitude asymmetric décollement fold in the foreland (right sidewall). Upper layers dip as moderately to steeply as fault-related folds develop and were extended and thinned in crestal areas with possible extensional faulting and collapse structures as a result of buckling of the uppermost hanging-wall strata, whereas they were folded and thickened in adjacent synclinal areas (Figure 4b).

Model 9 (Figure 4c) exhibits that a thick décollement layer impacts the structural style and evolution of a sedimentary cover sequence composed of four strong bonded layers intercalated with weak bonded layers (three weak layers thick). In this model, the décollement layer is 0.3 km thick. Overall, the deformation is mainly pure-shear thickening (contraction or buckling) with the hinterland-verging low-amplitude asymmetric décollement fold (F1) at the initiation of deformation following an initial step of homogeneous shortening of the model. With the passage of time, the deformity becomes more evident (deformation 20%) and develops two hinterland-verging fault-propagation folds with localized strain zones and distinct fault segments developed in lower parts and core of folds in the hinterland (left sidewall) and foreland (right sidewall). In this case, thrusts would preferentially develop to accommodate applied horizontal shortening rather than buckling. As the deformation progresses and after 30% shortening, the fault-propagation folds further develop and their asymmetries increase, however, a new well-developed fault-propagation fold F2, is also formed at the middle part of the model. Increasing shortening causes both flexural slip and flexural flow, respectively, in high-strain zones and in the lower parts of fault-related folds so that they are filled with weak material that evacuated beneath synclines and flowed into the anticline core. Upper layers dip as moderately to steeply (15° to 25°) as fault-related folds develop and were extended and thinned in crestal areas with possible extensional faulting and collapse structures as a result of buckling of the uppermost hanging-wall strata, whereas they were folded and thickened in adjacent synclinal areas (Figure 4c).

4. Discussion

4.1. Comparison to natural examples

Using a series of DE-models, we have model circumstances and structures in the contractional environment during progressive deformation. We make no attempt to simulate the exact structure of any natural example; however, in order to validate the final results, we compared one of the modelling results with a well-exposed décollement fold example developed in the Central Iran structural zone, near Isfahan city (Lat.: 33° 09′ 30′′ N, Long.: 51° 46′ 30′′ E) (Figure 5). In this area, the upper layers consisted of limestone and lower part included shale and silty shale so that there are interesting structures in this area such as asymmetric box-shape décollement fold (Figure 5d). The shape and characteristic of box-shape model that is depicted in the modelling results (e.g. Model 4 in Figure 3) are compatible with the mentioned box fold in the nature that consisted of weak and strong layers as overburden and the silty shale as the décollement layer (Figure 6). The nature of weak layer in lower part and the décollement thickness is similar to the structure in the field.

Figure 5.

(a) Simplified topographic map of central Iran. The rectangle indicates the natural example located at central Iran structural zone, near Isfahan city (Lat.: 33° 09′ 30′′ N, Long.: 51° 46′ 30′′ E). (b) Geological map that is part of 1:100,000 TARQ quadrangle map from geological survey of Iran showing the location of the natural example. The lithology of the structure consists of limestone for upper layers and shale for lower part. The rectangle shows the exact position of the studied structure. (c) Perspective of the structure in nature near the road between Natanz city and Esfahan city in Iran.

Figure 6.

(a) The geometry and characteristic of box-shape structure are compatible with (b) the box fold in model 4 (depicted in Figure 3) in final stage of shortening that consisted of weak and strong layers as overburden and the silty shale as the décollement layer.

In addition to the natural example mentioned above, we can also compare them to other regional and tectonic-scale examples in order to verify the applicability of numerical models, for example, the Zagros fold-and-thrust belt (Iran). The Zagros fold-and-thrust belt is a structural zone which is associated with the significant influence of mechanical stratigraphy through sedimentary sequences and multi-décollement layers [e.g., Sherkati et al. 2005, 2006, Alavi 2007, Vergés et al. 2011]. Furthermore, the Zagros fold-and-thrust belt is often recognized as one of the regions in the world showing the best examples of large-scale décollement folds as shown by their regular and well-rounded geometries [Sherkati et al. 2005, Meng and Hodgetts 2019b]. The result of this study showed that the thickness changes of décollement horizon is one of the controlling factors of variations in the geometry, also as seen in the other example named the Zeloi anticline in the Zagros belt (SW Iran) that is a symmetric and asymmetric faulted décollement fold.

Another good example is the analysis conducted by Gu et al. [2021] about large-scale fold-and-thrust belt in South China. The stratigraphic succession is composed of competent layers separated by three main incompetent layers which act as décollements. The basal décollement, the Cambrian evaporites, played a dominant role in the deformation; the intermediate and upper décollement accommodated the displacement during deformation. Structural styles are dominated by faulted décollement folds with the breakthroughs of the fore- and/or backlimbs.

Comparison of the present study with previous researches show that the décollement thickness is one of the main controlling factors of the structural evolution and style of the tectonic setting.

4.2. Variation in geometry and kinematics of models as a function of décollement thickness and mechanical stratigraphy

Our modelling results indicated that the geometry and kinematics of structures developed above a décollement level significantly changes in dependence on the décollement thickness and number of weak layers, which are considered as our main variables in the present study. Décollement folds can also develop ahead of the upper tip of a bedding-parallel thrust [Poblet and McClay 1996] in more competent rocks and are often filled by less competent rocks [Homza and Wallace 1995]. Distinct décollement folds can only form if there is enough ductile material present in the décollement layer to penetrate the core of the growing anticline [Stewart 1996]. If there is not enough material to fill the core of the fold, fold growth will finish and shortening may be completed by faulting [Stewart 1996] as we can see in the models with less décollement thickness such as Model 1 that faulting developed by increasing shortening. Mitra [2002, 2003] presented a kinematic model for décollement folds that explains why a large number of geometrically different décollement folds can occur as seen in Figures 23, and 4. Based on the modelling results it concluded that the fold develops in a competent rock layer that is underlain by a less competent unit. The fold mainly deforms by flexural slip with some additional faulting and fracturing. The wavelength of the fold is controlled by the thickness of the involved rock layers. Successive deformation leads to increasing amplitude and wavelength of the fold. This is motivated by the rotation of limb segments and hinge migration [Mitra 2003]. Looking into models such as Model 4 in Figure 3 revealed that in the models defined with more number of layers, the deformation formed by more complicated structures that is comparable with the findings presented by Mitra [2003] and it is possible to compare the décollement folding styles showed in (Figure 3) with the structural variations that is presented by Mitra [2003]. The similarity of the numerical models and kinematic models are significantly presented that the prepared numerical models are similar to the algorithms in nature (Figure 6).

As described in model results, contractional fault-related folds that formed as a result of horizontal shortening exhibit varied geometries especially wavelength and amplitude. Generally, it is demonstrated that the thicker basal décollement layer result in the higher the maximum flexural flow height. Similarly, surface uplift and fold amplification are mainly progressively affected by the basal décollement thickness so that uplift and fold amplitude increase from model with a thin décollement to models with intermediate to thick décollement. In this regard, the fold amplification (i.e., vertical growth of folds) with fold cores filled with basal materials through flexural flow led to a more significant uplift. Moreover, the number of ideal décollement folds and thrust-related folds and also their geometry developed in models influenced by the material strength and the décollement layer thickness, so that models with thick décollement layer result in more thrust-related folds and less ideal décollement folds at a specific shortening value. Box-like folds tend to show tight cores and broad, flat crests. Therefore, model results show that décollement thickness is able of effecting the fold geometries. With the growth of the fold, the fold limbs lengthen and become steeper, because of hinge migration and limb rotation [Mitra 2003]. As the fold tightens and the anticlinal area increases, the synclinal area decreases. Therefore, in the late stage of the deformation, the units almost return to their regional positions in the synclines [Li and Mitra 2020]. This mechanism is clearly visible in the models such as Models 6 and 8. It should be noted that this study considered just two main variables that are important in the formation of décollement folds, and the presented models are by no means perfect. Although considering more parameters led to hardness in distinguishing the determinants, for more understanding of the evolution of these kinds of folds, it is suggested that future models should incorporate other sets of controlling factors that are not considered in this paper.

5. Conclusions

In this research, we created a series of discrete element models (DEM) to assess the importance of geometric and mechanical factors in the evolution of décollement folds in sedimentary layers. The results reached from a series of modelling process showed us that shortening was mainly accommodated by décollement folds (sinusoidal, bulb- and box-like folds), fault-propagation folds, with thinning and thickening in the models. By changing the considered variables such as décollement thickness and number of weak and strong sedimentary upper layers, a range of models produced. The first groups or models with a thin décollement (Figure 2) accommodated by more distributed strain and décollement folds as box, lift-off and disharmonic folds, whereas more localized strain and contractional fault-related folds, such as faulted décollement and fault-propagation folds, are developed in intermediate to thick décollement models. Surface uplift and fold amplification are mainly progressively affected by the basal décollement thickness so that uplift and fold amplitude increase from model with a thin décollement to models with intermediate to thick décollement. Moreover, model results show that as the number of weak layers’ decreases, shortening was mainly accommodated by brittle deformation and faulting. Therefore, it could be concluded that with intermediate and thick décollement, the structural styles that developed are more complex. Finally, comprising the modelling results with the structures in the field such as central part of Iran or Zagros fold-and-thrust belt demonstrated that thickness changes of the décollement level and the overburden layers are both contributing factors in the structural styles of décollement folding during horizontal shortening.

Conflicts of interest

Authors have no conflict of interest to declare.

Acknowledgement

We greatly thank Professor Martin Schöpfer (University of Wien) for his careful help and comments regarding numerical modelling process.

1The velocity-Verlet integration is a numerical method used to integrate Newton’s equations of motion as this integrator strike a good balance between computational cost, stability, and accuracy.


Bibliographie

[Abe et al., 2011] S. Abe; H. Van Gent; J. L. Urai DEM simulation of normal faults in cohesive materials, Tectonophysics, Volume 512 (2011), pp. 12-21 | DOI

[Alavi, 2007] M. Alavi Structures of the Zagros fold-thrust belt in Iran, Am. J. Sci., Volume 307 (2007), pp. 1064-1095 | DOI

[Allmendinger, 1998] R. W. Allmendinger Inverse and forward numerical modeling of trishear fault propagation folds, Tectonics, Volume 17 (1998), pp. 640-656 | DOI

[Benesh et al., 2007] N. P. Benesh; A. Plesch; J. H. Shaw; E. K. Frost Investigation of growth fault bend folding using discrete element modeling: implications for signatures of active folding above blind thrust faults, J. Geophys. Res. Solid Earth, Volume 112 (2007), B03S04 | DOI

[Bonini, 2003] M. Bonini Detachment folding, fold amplification, and diapirism in thrust wedge experiments, Tectonics, Volume 22 (2003) no. 6, 1065 | DOI

[Botter et al., 2014] C. Botter; N. Cardozo; S. Hardy; I. Lecomte; A. Escalona From mechanical modeling to seismic imaging of faults: A synthetic workflow to study the impact of faults on seismic, Mar. Pet. Geol., Volume 57 (2014), pp. 187-207 | DOI

[Brandes and Tanner, 2014] C. Brandes; D. C. Tanner Fault-related folding: a review of kinematic models and their application, Earth Sci. Rev., Volume 138 (2014), pp. 352-370 | DOI

[Briggs et al., 2006] S. E. Briggs; R. J. Davies; J. A. Cartwright; R. Morgan Multiple detachment levels and their control on fold styles in the compressional domain of the deep-water west Niger Delta, Basin Res., Volume 18 (2006), pp. 435-450 | DOI

[Burberry, 2015] C. M. Burberry Spatial and temporal variation in penetrative strain during compression: Insights from analog models, Lithosphere, Volume 6 (2015), pp. 611-624 | DOI

[Burbidge and Braun, 2002] D. R. Burbidge; J. Braun Numerical models of the evolution of accretionary wedges and fold-and-thrust belts using the distinct-element method, Geophys. J. Int., Volume 148 (2002), pp. 542-561 | DOI

[Butler et al., 2018] R. Butler; C. E. Bond; M. A. Cooper; H. Watkins Interpreting structural geometry in fold-thrust belts: Why style matters, J. Struct. Geol., Volume 114 (2018), pp. 251-273 | DOI

[Cardozo et al., 2005] N. Cardozo; R. W. Allmendinger; J. K. Morgan Influence of mechanical stratigraphy and initial stress state on the formation of two fault propagation folds, J. Struct. Geol., Volume 27 (2005), pp. 1954-1972 | DOI

[Chen et al., 2018] Z. Chen; Y. Lei; Y. Hu; L. Wang; G. Yang Structural analysis of multi-level detachments and identification of deep-seated anticline, Pet. Explor. Dev., Volume 45 (2018) no. 2, pp. 281-289 | DOI

[Corredor et al., 2005] F. Corredor; J. H. Shaw; F Bilotti Structural styles in the deep-water fold and thrust belts of the Niger Delta, Am. Assoc. Pet. Geol. Bull., Volume 89 (2005) no. 6, pp. 753-780 | DOI

[Cosgrove, 2015] J. Cosgrove The association of folds and fractures and the link between folding, fracturing and fluid flow during the evolution of a fold–thrust belt: a brief review, Geol. Soc. Lond. Spec. Publ., Volume 421 (2015), pp. 41-68 | DOI

[Cui et al., 2020] J. Cui; D. Jia; H. Yin; Zh. Chen; Y. Li; M. Wang; F. Xiaogen; L. Shen; Ch. Sun; D. Ma; Y. Zhang The influence of a weak upper ductile detachment on the Longmen Shan fold-and-thrust belt (eastern margin of the Tibetan Plateau): Insights from sandbox experiments, J. Asian Earth Sci., Volume 200 (2020), 104507 | DOI

[Cundall and Strack, 1979] P. A. Cundall; O. D. L. Strack A discrete numerical model for granular assemblies, Geotechnique, Volume 29 (1979) no. 1, pp. 47-65 | DOI

[Dean et al., 2013] S. L. Dean; J. K. Morgan; T. Fournier Geometries of frontal fold and thrust belts: insights from discrete element simulations, J. Struct. Geol., Volume 53 (2013), pp. 43-53 | DOI

[Deng et al., 2017] C. Deng; R. L. Gawthorpe; E. Finch; H. Fossen Influence of a pre-existing basement weakness on normal fault growth during oblique extension: insights from discrete element modeling, J. Struct. Geol., Volume 105 (2017), pp. 44-61 | DOI

[Derikvand et al., 2018] B. Derikvand; S. A. Alavi; I. Abdollahie Fard; H. Hajialibeigi Folding style of the Dezful Embayment of Zagros Belt: Signature of detachment horizons, deep-rooted faulting and syn-deformation deposition, Mar. Pet. Geol., Volume 91 (2018), pp. 501-518 | DOI

[Donzé et al., 1994] F. Donzé; P. Mora; S.-A. Magnier Numerical simulation of faults and shear zones, Geophys. J. Int., Volume 116 (1994), pp. 46-52 | DOI

[Farzipour-Saein and Koyi, 2014] A. Farzipour-Saein; H. Koyi Effect of lateral thickness variation of an intermediate decollement on the propagation of deformation front in the Lurestan and Izeh zones of the Zagros fold-thrust belt, insights from analogue modeling, J. Struct. Geol., Volume 65 (2014), pp. 17-32 | DOI

[Farzipour-Saein, 2017] A. Farzipour-Saein Folding style controlled by intermediate decollement thickness change in the Lurestan region (NW of the Zagros fold-and-thrust belt), using analogue models, Int. J. Earth Sci., Volume 106 (2017), pp. 1525-1537 | DOI

[Finch and Gawthorpe, 2017] E. Finch; R. Gawthorpe Growth and interaction of normal faults and fault network evolution in rifts: insights from three-dimensional discrete element modelling, Geol. Soc. Lond. Spec. Publ., Volume 439 (2017), pp. 219-248 | DOI

[Finch et al., 2003] E. Finch; S. Hardy; R. Gawthorpe Discrete element modelling of contractional fault-propagation folding above rigid basement fault blocks, J. Struct. Geol., Volume 25 (2003), pp. 515-528 | DOI

[Ghanadian et al., 2017] M. Ghanadian; A. Faghih; I. Abdollahie Fard; T. Kusky; M. Maleki On the role of incompetent strata in the structural evolution of the Zagros Fold-Thrust Belt, Dezful embayment, Iran, Mar. Pet. Geol., Volume 81 (2017), pp. 320-333 | DOI

[Gu et al., 2021] Z. Gu; X. Wang; A. Nunns; B. Zhang; H. Jiang; L. Fu; X. Zhai Structural styles and evolution of a thin-skinned fold-and-thrust belt with multiple detachments in the eastern Sichuan Basin, South China, J. Struct. Geol., Volume 142 (2021), 104191 | DOI

[Hammerstein et al., 2020] A. Hammerstein; R. Di Cuia; M. A. Cottam; G. Zamora; R. W. H. Butler Fold and thrust belts: Structural style, evolution and exploration—An introduction, Geol. Soc. Lond. Spec. Publ., Volume 490 (2020), pp. 1-8 | DOI

[Hardy and Finch, 2005] S. Hardy; E. Finch Discrete-element modelling of detachment folding, Basin Res., Volume 17 (2005), pp. 507-520 | DOI

[Hardy and Finch, 2006] S. Hardy; E. Finch Discrete element modelling of the influence of cover strength on basement-involved fault-propagation folding, Tectonophysics, Volume 415 (2006), pp. 225-238 | DOI

[Hardy and Finch, 2007] S. Hardy; E. Finch Mechanical stratigraphy and the transition from trishear to kink-band fault-propagation fold forms above blind basement thrust faults: a discrete element study, Mar. Pet. Geol., Volume 24 (2007), pp. 75-90 | DOI

[Hardy et al., 2009] S. Hardy; K. McClay; J. A. Munoz Deformation and fault activity in space and time in high-resolution numerical models of doubly vergent thrust wedges, Mar. Pet. Geol., Volume 26 (2009), pp. 232-248 | DOI

[Hardy, 2011] S. Hardy Cover deformation above steep, basement normal faults: insights from 2D discrete element modeling, Mar. Pet. Geol., Volume 28 (2011), pp. 966-972 | DOI

[Hardy, 2013] S. Hardy Propagation of blind normal faults to the surface in basaltic sequences: insights from 2D discrete element modelling, Mar. Pet. Geol., Volume 48 (2013), pp. 149-159 | DOI

[Hardy, 2019] S. Hardy Discrete element modelling of extensional, growth, fault-propagation folds, Basin Res., Volume 31 (2019), pp. 584-599 | DOI

[Homza and Wallace, 1995] T. X. Homza; W. K. Wallace Geometric and kinematic models for detachment folds with fixed and variable depths, J. Struct. Geol., Volume 17 (1995), pp. 575-588 | DOI

[Hughes and Shaw, 2015] A. N. Hughes; J. H. Shaw Insights into the mechanics of fault-propagation folding styles, Geol. Soc. Am. Bull., Volume 127 (2015) no. 11-12, pp. 1752-1765 | DOI

[Hughes et al., 2014] A. N. Hughes; N. P. Benesh; J. H. Shaw Factors that control the development of fault-bend versus fault-propagation folds: insights from mechanical models based on the discrete element method (DEM), J. Struct. Geol., Volume 68 (2014), pp. 121-141 | DOI

[Imber et al., 2004] J. Imber; G. W. Tuckwell; C. Childs; J. J. Walsh; T. Manzocchi; A. E. Heath; C. G. Bonson; J. Strand Three-dimensional distinct element modelling of relay growth and breaching along normal faults, J. Struct. Geol., Volume 26 (2004), pp. 1897-1911 | DOI

[Jamison, 1987] W. Jamison Geometric analysis of fold development in overthrust terranes, J. Struct. Geol., Volume 9 (1987), pp. 207-219 | DOI

[Kendall et al., 2019] J. Kendall et al. Petroleum tectonic comparison of fold and thrust belts: The Zagros of Iraq and Iran, the Pyrenees of Spain, the Sevier of Western USA and the Beni Sub-Andean of Bolivia, Geol. Soc. Lond. Spec. Publ., Volume 490 (2019) no. 1, pp. 79-103 (2020) | DOI

[Lammie et al., 2020] D. Lammie et al. Quantifying shortening across the central Appalachian fold-thrust belt, Virginia and West Virginia, USA: Reconciling grain-, outcrop-, and map-scale shortening, Geosphere, Volume 16 (2020) no. 5, pp. 1276-1292 | DOI

[Li and Mitra, 2017] J. Li; S. Mitra Geometry and evolution of fold-thrust structures at the boundaries between frictional and ductile detachments, Mar. Pet. Geol., Volume 85 (2017), pp. 16-34 | DOI

[Li and Mitra, 2020] J. Li; S. Mitra Seismic models of detachment and faulted detachment folds, Mar. Pet. Geol., Volume 117 (2020), 104385 | DOI

[Marshak et al., 1992] S. Marshak; M. S. Wilkerson; A. T. Hsui Generation of curved fold-thrust belts: Insight from simple physical and analytical models, Thrust Tectonics (K. R. McClay, ed.), Springer, Dordrecht, 1992 | DOI

[McClay et al., 2004] K. R. McClay; P. S. Whitehouse; T. Dooley; M. Richards 3D evolution of fold and thrust belts formed by oblique convergence, Mar. Pet. Geol., Volume 21 (2004) no. 7, pp. 857-877 | DOI

[Meng and Hodgetts, 2019a] Q. Meng; D. Hodgetts Structural styles and decoupling in stratigraphic sequences with double décollements during thin-skinned contractional tectonics: Insights from numerical modelling, J. Struct. Geol., Volume 127 (2019a), 103862 | DOI

[Meng and Hodgetts, 2019b] Q. Meng; D. Hodgetts Combined control of décollement layer thickness and cover rock cohesion on structural styles and evolution of fold belts: A discrete element modelling study, Tectonophysics, Volume 757 (2019b), pp. 58-67 | DOI

[Mitra and Namson, 1989] S. Mitra; J. Namson Equal-area balancing, Am. J. Sci., Volume 289 (1989), pp. 563-599 | DOI

[Mitra, 1990] S. Mitra Fault-propagation folds: geometry, kinematic evolution, and hydrocarbon traps, Am. Assoc. Pet. Geol. Bull., Volume 74 (1990), pp. 921-945 | DOI

[Mitra, 2002] S. Mitra Structural models of faulted detachment folds, Am. Assoc. Petrol. Geol. Bull., Volume 86 (2002) no. 9, pp. 1673-1694 | DOI

[Mitra, 2003] S. Mitra A unified kinematic model for the evolution of detachment folds, J. Struct. Geol., Volume 25 (2003), pp. 1659-1673 | DOI

[Morgan and Bangs, 2017] J. K. Morgan; N. L. Bangs Recognizing seamount-forearc collisions at accretionary margins: insights from discrete numerical simulations, Geology, Volume 45 (2017), pp. 635-638 | DOI

[Morgan, 2015] J. K. Morgan Effects of cohesion on the structural and mechanical evolution of fold and thrust belts and contractional wedges: discrete element simulations, J. Geophys. Res. Solid Earth, Volume 120 (2015), pp. 3870-3896 | DOI

[Morley et al., 2011] C. K. Morley; M. Tingay; R. Hillis; R. King Deepwater fold and thrust belt classification, tectonics, structure and hydrocarbon prospectivity: a review, Earth Sci. Rev., Volume 104 (2011), pp. 41-91 | DOI

[Morley et al., 2017] C. K. Morley; C. von Hagke; R. L. Hansberry; A. S. Collins; W. Kanitpanyacharoen; R. King Review of major shale-dominated detachment and thrust characteristics in the diagenetic zone: Part I, meso and macro-scopic scale, Earth Sci. Rev., Volume 173 (2017), pp. 168-228 | DOI

[Motamedi and Gharabeigli, 2019] Hossein Motamedi; Gholamreza Gharabeigli Structural Style in the Fars Geological Province: Interaction of Diapirism and Multidetachment Folding, Dev. Struct. Geol. Tecton., Volume 3 (2019), pp. 145-160 | DOI

[Nabavi et al., 2018] S. T. Nabavi; S. A. Alavi; F. Maerten 2D finite-element elastic models of transtensional pull-apart basins, C. R. Geosci., Volume 350 (2018), pp. 222-230 | DOI

[Nabavi et al., 2019] S. T. Nabavi; S. A. Alavi; H. J. Javanbakht The Dinevar transtensional pull-apart basin, NW Zagros Mountains, Iran: a geological study and comparison to 2D finite element elastic models, Intl J. Earth Sci., Volume 108 (2019), pp. 329-346 | DOI

[Nabavi et al., 2020] S. T. Nabavi; S. A. Alavi; M. Díaz-Azpiroz; S. Mohammadi; M. R. Ghassemi; C. Fernández; L. Barcos; M. Frehner Deformation mechanics in inclined, brittle-ductile transpression zones: Insights from 3D finite element modelling, J. Struct. Geol., Volume 137 (2020), 104082 | DOI

[Naylor et al., 2005] M. Naylor; H. D. Sinclair; S. Willett; P. A. Cowie A discrete element model for orogenesis and accretionary wedge growth, J. Geophys. Res.: Solid Earth, Volume 110 (2005), B12403 | DOI

[Pla et al., 2019] O. Pla; E. Roca; H. W. Xie; E. Izquierdo-Llavall; J. A. Muñoz; M. G. Rowan; O. Ferrer; O. Gratacós; Y. Neng; S. Y. Huang Influence of syntectonic sedimentation and décollement rheology on the geometry and evolution of orogenic wedges: analog modeling of the Kuqa fold-and-thrust belt (NW China), Tectonics, Volume 38 (2019), pp. 2727-2755 | DOI

[Poblet and McClay, 1996] J. Poblet; K. McClay Geometry and kinematics of single-layer detachment folds, Am. Assoc. Pet. Geol. Bull., Volume 80 (1996), pp. 1085-1109 | DOI

[Ridley and Casey, 1989] J. Ridley; M. Casey Numerical modeling of folding in rotational strain histories: Strain regimes expected in thrust belts and shear zones, Geology, Volume 17 (1989) no. 10, pp. 875-878 | DOI

[Rowan et al., 2019] M. G. Rowan et al. Deformation of intrasalt competent layers in different modes of salt tectonics, Solid Earth, Volume 10 (2019), pp. 987-1013 | DOI

[Saltzer and Pollard, 1992] S. D. Saltzer; D. D. Pollard Distinct element modeling of structures formed in sedimentary overburden by extensional reactivation of basement normal faults, Tectonics, Volume 11 (1992), pp. 165-174 | DOI

[Schori et al., 2015] M. Schori; J. Mosar; G. Schreurs Multiple detachments during thin-skinned deformation of the Swiss Central Jura: a kinematic model across the Chasseral, Swiss J. Geosci., Volume 108 (2015), pp. 327-343 | DOI

[Schöpfer et al., 2006] M. P. J. Schöpfer; C. Childs; J. J. Walsh Localisation of normal faults in multilayer sequences, J. Struct. Geol., Volume 28 (2006), pp. 816-833 | DOI

[Schöpfer et al., 2016] M. P. J. Schöpfer; C. Childs; T. Manzocchi; J. J. Walsh Three-dimensional Distinct Element Method modelling of the growth of normal faults in layered sequences, Geol. Soc. Lond. Spec. Publ., Volume 439 (2016) no. 1, pp. 307-332 | DOI

[Schöpfer et al., 2017] M. P. J. Schöpfer; C. Childs; T. Manzocchi; J. J. Walsh; A. Nicol; B. Grasemann The emergence of asymmetric normal fault systems under symmetric boundary conditions, J. Struct. Geol., Volume 104 (2017), pp. 159-171 | DOI

[Sherkati et al., 2005] S. Sherkati; M. Molinaro; D. F. de Lamotte; J. Letouzey Detachment folding in the central and eastern Zagros fold-belt (Iran): salt mobility, multiple detachments and final basement control, J. Struct. Geol., Volume 27 (2005), pp. 1680-1696 | DOI

[Sherkati et al., 2006] S. Sherkati; J. Letouzey; D. Frizon de Lamotte Central Zagros fold-thrust belt (Iran): new insights from seismic data, field observation, and sandbox modeling, Tectonics, Volume 25 (2006), TC4007 | DOI

[Smart et al., 2011] K. J. Smart; D. Y. Wyrick; D. A. Ferrill Discrete element modeling of Martian pit crater formation in response to extensional fracturing and dilational normal faulting, J. Geophys. Res. Planets, Volume 116 (2011), E04005 | DOI

[Stewart, 1996] S. A. Stewart Influence of detachment layer thickness on style of thin-skinned shortening, J. Struct. Geol., Volume 18 (1996), pp. 1271-1274 | DOI

[Suppe et al., 2004] J. Suppe; C. D. Connors; Y. Zhang Shear fault-bend folding, Thrust Tectonics and Hydrocarbon Systems (K. R. McClay, ed.), Volume 82, American Association of Petroleum Geologists, Tulsa, OK, 2004 | DOI

[Suppe, 2011] J. Suppe Mass balance and thrusting in detachment folds, Thrust Fault-Related Folding (K. McClay; J. H. Shaw; J. Suppe, eds.) (AAPG Memoir 94), American Association of Petroleum Geologists, Tulsa, OK, 2011, pp. 21-37 | DOI

[Vergés et al., 2011] J. Vergés; E. Saura; E. Casciello; M. Fernandez; A. Villaseñor; I. Jimenez-Munt; D. García-Castellanos Crustal-scale cross-sections across the NW Zagros belt: implications for the Arabian margin reconstruction, Geol. Mag., Volume 148 (2011), pp. 739-761 | DOI

[Wallace and Homza, 2004] W. K. Wallace; T. Homza Detachment folds versus fault-propagation folds, and their truncation by thrust faults, Thrust Tectonics and Hydrocarbon Systems (K. R. McClay, ed.) (AAPG Memoir 82), American Association of Petroleum Geologists, Tulsa, OK, 2004, pp. 303-323

[Yamada et al., 2006] Y. Yamada; K. Baba; T. Matsuoka Analogue and numerical modelling of accretionary prisms with a décollement in sediments, Geol. Soc. Lond. Spec. Publ., Volume 253 (2006) no. 1, pp. 169-183 | DOI

[Yarbuh et al., 2018] I. Yarbuh; A. González; R. M. Spelz-Madero; R. Negrete-Aranda; J. Contreras Development of detachment folds in the Mexican Ridges foldbelt, western Gulf of Mexico basin, Tectonics, Volume 37 (2018), pp. 2013-2028 | DOI


Commentaires - Politique