Ebro Basin continental sedimentation associated with late compressional Pyrenean tectonics(north‐eastern Iberia): controls on basin margin fans and fluvial systems

ABSTRACT This contribution deals with the External Sierras and a part of the foreland Ebro Basin related to the southern Pyrenean thrust front. The structure of the External Sierras consists of a south‐verging thrust system developed from middle Eocene to early Miocene times. Since the end of the early Oligocene, a regional‐scale detachment anticline (the Santo Domingo anticline) developed, folding the original thrust system and creating new thrust units. The molassic fill in this part of the Ebro Basin (Uncastillo Formation) mainly corresponds to an extensive, composite distributary fluvial system, termed the Luna system, which drained the uplifted Gavarnie Unit to the north. Small, marginal alluvial fans originated along the External Sierras and coalesced in the proximal‐middle portions of the Luna system. Three tecto‐sedimentary units (TSU), late Oligocene to early Miocene in age, comprise the Uncastillo Formation. Lateral relationships and areal distribution of lithofacies through time have been used to establish sedimentary models for the marginal alluvial fans and the Luna fluvial system. Their sedimentary evolution was controlled by tectonics affecting the drainage basins, and based on mapping and stratigraphic relationships of the TSU, the temporal succession of the marginal alluvial fans and their relationships with each thrust system in the south Pyrenean front can be shown. Alluvial fan formation evolved through time from west to east, in accord with the progressive eastward growth of the Santo Domingo anticline as a conical fold. The fluvial network of the Luna system appears to have been mainly radial, but near the basin margin its architecture was influenced by the syndepositional Fuencalderas and Uncastillo anticlines developed within the Ebro Basin. These low‐amplitude folds originated by layer‐parallel shearing caused by rotation of the southern flank of the Santo Domingo anticline. Progressive uplift of these anticlines constrained part of the fluvial discharge to synclinal areas parallel to the basin margin; these areas where characterized by meandering sandy channels. At the peripheral tips of the anticlines the channel system flowed basinward.     

Low-amplitude, synsedimentary folding of a deltaic complex: Roda Sandstone (lower Eocene), South-Pyrenean Foreland Basin

The South‐Pyrenean Foreland Basin is a sedimentary trough developed during Palaeogene times in response Pyrenean orogenesis. The structural development of the chain progressively involved the foreland basin deposits resulting in synsedimentary thrusting and growth of folds through the basin. The lower Eocene Roda Sandstone was deposited in a shallow‐marine environment whose topography (bathymetry) was modified by a series of growing gentle folds. This synsedimentary folding is evidenced in the field (i) thickening of sedimentary units above synclinal structures and thinning over anticlines; (ii) carbonate platform deposits growing on top of anticlines; (iii) the areal distribution of benthic foraminifera in transgressive facies assemblages determined by an irregular, fold‐influenced palaeobathymetry; (iv) variation of sandstone palaeocurrents related to the presence of a sedimentary trough formed by the synsedimentary growth of a syncline. In addition, synsedimentary folding has been evidenced from seismic data. In the Roda Sandstone example, the growth of gentle folds occurred in an area with high sedimentation rates (～0.21 ± 0.06 mm y^?1). Due to the high sedimentation rates, exceeding the folds uplift rate (～0.10 ± 0.01 mm y^?1), there are no noticeable unconformities in the growth strata at outcrop scale. However, the effects on the sedimentation are very significant, because the sediments were deposited close to sea level, and thus were very sensitive to fauna and facies distribution.     

Syn‐orogenic fluid flow in the Jaca basin (south Pyrenean fold and thrust belt) from fracture and vein analyses

This study aims at understanding the origin and nature of syn‐orogenic fluid flow in the Jaca basin from the South Pyrenean fold and thrust‐belt, as recorded in calcite and quartz veins of the Sierras Interiores (Spain) and the turbiditic basin, which cover upper Cretaceous to Late Eocene syntectonic deposits. The fracture network consists of a classical pattern of transverse and longitudinal fractures related to Layer Parallel Shortening (LPS) and folding respectively. Veins filled equally about the third of fractures in the carbonate shelf and turbidites. Carbon and oxygen isotopes of calcite veins mostly indicate precipitation from isotopically buffered water, consistent with high water‐rock interaction. In the Sierras Interiores, petrographical observations and fluid inclusion microthermometry are consistent with two distinct stages of precipitation. The first stage is characterized by relatively low Th and low salinities (155–205 °C and 0.5–3.2 wt% eq. NaCl). The second stage, which was characterized both by the formation of mode‐I joints and by mode‐I reactivation of pre‐existing veins, shows higher Th and salinities (215–270 °C and 2.2–5.7 wt% eq. NaCl). Waters recorded in the second stage are interpreted to have interacted with underlying Triassic evaporites and flowed along major thrusts before vein precipitation, which are locally in thermal disequilibrium with host‐rocks. We suggest the transition from a rather closed hydrological system during the first stage of vein formation, interpreted to have occurred during Eaux‐chaudes thrusting (upper Lutetian‐Bartonian), to a more open hydrological system during the second stage, which likely occurred during Gavarnie thrusting (Priabonian‐early Rupelian). Finally, we also document the migration in space and time of hydrothermal pulses along the South Pyrenean Foreland Basin, related to the westward propagation of major thrusts during the Pyrenean orogeny.     

Magnetostratigraphy and detrital apatite fission track thermochronology in syntectonic conglomerates: constraints on the exhumation of the South‐Central Pyrenees

The syntectonic continental conglomerates of the South‐Central Pyrenees record the late stages of thin‐skinned transport of the South‐Pyrenean Central Units and the onset of exhumation of the Pyrenean Axial Zone (AZ) in the core of the orogen. New magnetostratigraphic data of these syntectonic continental conglomerates have established their age as Late Lutetian to Late Oligocene. The data reveal that these materials were deposited during intense periods of tectonic activity of the Pyrenean chain and not during the cessation of the deformation as considered previously. The magnetostratigraphic ages have been combined with new detrital apatite fission track (AFT) thermochronology from AZ‐derived granite cobbles within the syntectonic conglomerates. Distribution of the granitic cobbles with different AFT ages and track lengths combined with their depositional ages reveal information on the timing and rate of episodes of exhumation in the orogen. Some AFT ages are considerably older than the AFT ages of the outcropping AZ granitic massifs, indicating erosion from higher crustal levels within the massifs than presently exposed or from completely eroded plutons. Inverse thermal modelling reveals two well‐defined periods of rapid cooling in the hinterland at ca. 50–40 and ca. 30–25 Ma, with another poorly defined cooling episode at ca. 70–60 Ma. The lowest stratigraphic samples experienced postburial annealing caused by the deposition of younger syntectonic sediments during progressive burial of the south Pyrenean thrust and fold belt. Moreover, samples from the deeper stratigraphic levels also reveal postorogenic cooling during the Late Miocene as a response to the excavation of the Ebro River towards the Mediterranean Sea. Our data strongly support previous ideas about the burial of the South Pyrenean fold and thrust belt by Late Palaeogene syntectonic conglomerates and their subsequent re‐excavation and are consistent with other thermochronological data and thermal modelling from the interior part of the orogen.     

Discrete element modelling of extensional,growth, fault‐propagation folds

Extensional fault‐propagation folds are now recognised as being an important part of basin structure and development. They have a very distinctive expression, often presenting an upward‐widening monocline, which is subsequently breached by an underlying, propagating fault. Growth strata, if present, are thought to provide a crucial insight into the manner in which such structures grow in space and time. However, interpreting their stratigraphic signal is neither straightforward nor unique. Both analogue and numerical models can provide some insight into fold growth. In particular, the trishear kinematic model has been widely adopted to explain many aspects of the evolution and geometry of such fault‐propagation folds. However, in some cases the materials/rheologies used to represent the cover do not reproduce the key geometric/stratigraphic features of such folds seen in nature. This appears to arise from such studies not addressing adequately the very heterogenous mechanical stratigraphy seen in many sedimentary covers. In particular, flexural slip between beds/layers is often not explicitly modelled but, paradoxically, it appears to be an important deformation mechanism operative in such settings. Here, I present a 2D discrete element model of extensional fault‐propagation folding which explicitly includes flexural slip between predefined sedimentary units or layers in the cover. The model also includes growth strata and shows how they may reflect the various evolutionary stages of fold and fault growth. When flexural slip is included in the modelling scheme, the resultant breached monoclines and their growth strata are strikingly similar to some of those seen in nature. Results are also compared with those obtained using simple, homogeneous, frictional‐cohesive and elastic cover materials. Both un‐lithified and lithified growth strata are considered and clearly show that, rather than just being passive recorders of structural evolution, growth strata can themselves have an important effect on fault‐related fold growth. Implications for the evolution of and strain within, the resultant growth structures are discussed. A final focus of this study is the relationship that trishear might have with the upward‐widening zone of flexural slip activation away from a fault tip singularity.     

Along‐strike evolution of folding,stretching and breaching of supra‐salt strata in the Plataforma Burgalesa extensional forced fold system (northern Spain)

The Plataforma Burgalesa is a partly exposed extensional forced fold system with an intermediate salt layer, which has developed along the southern portion of the Basque‐Cantabrian Basin from Malm to Early Cretaceous as part of the Bay of Biscay‐Pyrenean rift system. Relationships between syn‐ and pre‐rift strata of the supra‐salt cover sequence and distribution of intra‐cover second‐order faults are observed both along seismic sections and at the surface. These relationships indicate an along‐strike variability of the extensional structural style. After a short period of salt mobilization and forced folding, high slip rates in the central portion of the major basement faults have rapidly promoted brittle behaviour of the salt layer, preventing further salt mobilization and facilitating the propagation of the fault across the salt layer. In contrast, at the tip regions of basement faults, slower slip rates have facilitated ductile salt behaviour, ensuring its further evaporite evacuation, preventing fault propagation across the salt layer and, in essence, allowing for a long‐living forced folding process. Our results indicate the important effect of along‐strike variation in displacement and displacement rates in controlling evaporite behaviour in extensional basins. Amount of displacement and displacement rates are key factors controlling the propagation of basement faults across evaporite layers. In addition, growth strata patterns are recognized as a powerful tool for constraining the up‐dip propagation history of basement faults in extensional fault‐related fold systems with intermediate décollement levels.     

Detrital record of Mesozoic shortening in the Yanshan belt, NE China: testing structural interpretations with basin analysis

The Yanshan fold‐thrust belt is an exposed portion of a major Mesozoic orogenic system that lies north of Beijing in northeast China. Structures and strata within the Yanshan record a complex history of thrust faulting characterized by multiple deformational events. Initially, Triassic thrusting led to the erosion of a thick sequence of Proterozoic and Palaeozoic sedimentary strata from northern reaches of the thrust belt; Triassic–Lower Jurassic strata that record this episode are deposited in a thin belt south of this zone of erosion. This was followed by postulated Late Jurassic emplacement of a major allochthon (the Chengde thrust plate), which is thought to have overridden structures and strata associated with the Triassic event and is cut by two younger thrusts (the Gubeikou and Chengde County thrusts). The Chengde allochthon is now expressed as a major east–west trending, thrust‐bounded synform (the Chengde synform), which has been interpreted as a folded klippe 20 km wide underlain by a single, north‐vergent thrust fault. Two sedimentary basins, defined on the basis of provenance, geochronology and palaeodispersal trends, developed within the Yanshan belt during Late Jurassic–Early Cretaceous time and are closely associated with the Chengde thrust and allied structures. Shouwangfen basin developed in the footwall of the Gubeikou thrust and records syntectonic unroofing of the hanging wall of that fault. Chengde basin developed in part atop Proterozoic strata interpreted as the upper plate of the Chengde allochthon and records unroofing of the adjacent Chengde County thrust. Both the Chengde County thrust and the Gubeikou thrust are younger than emplacement of the postulated Chengde allochthon, and structurally underlie it, yet neither Shouwangfen basin nor Chengde basin contain a detrital record of the erosion of this overlying structure. In addition, facies, palaeodispersal patterns and geochronology of Upper Jurassic strata that are cut by the Chengde thrust suggest only limited (ca. 5 km) displacement along this fault. We suggest that the units forming the Chengde synform are autochthonous, and that the synform is bounded by two limited‐displacement faults of opposing north and south vergence, rather than a single large north‐directed thrust. This conclusion implies that the Yanshan belt experienced far less Late Jurassic shortening than was previously thought, and has major implications for the Mesozoic evolution of the region. Specifically, we argue that the bulk of shortening and uplift in the Yanshan belt was accomplished during Triassic–Early Jurassic time, and that Late Jurassic structures modified and locally ponded sediments from a well‐developed southward drainage system developed atop this older orogen. Although Upper Jurassic strata are widespread throughout the Yanshan belt, it is clear that these strata developed within several discrete intermontane basins that are not correlable across the belt as a single entity. Thus, the Yanshan has no obvious associated foreland basin, and determining where the Mesozoic erosional products of this orogen ultimately lie is one of the more intriguing unresolved questions surrounding the palaeogeography of North China.     

Extensional fault and fold growth: Impact on accommodation evolution and sedimentary infill

Extensional faults and folds exert a fundamental control on the location, thickness and partitioning of sedimentary deposits on rift basins. The connection between the mode of extensional fault reactivation, resulting fault shape and extensional fold growth is well‐established. The impact of folding on accommodation evolution and growth package architecture, however, has received little attention; particularly the role‐played by fault‐perpendicular (transverse) folding. We study a multiphase rift basin with km‐scale fault displacements using a large high‐quality 3D seismic data set from the Fingerdjupet Subbasin in the southwestern Barents Sea. We link growth package architecture to timing and mode of fault reactivation. Dip linkage of deep and shallow fault segments resulted in ramp‐flat‐ramp fault geometry, above which fault‐parallel fault‐bend folds developed. The folds limited the accommodation near their causal faults, leading to deposition within a fault‐bend synclinal growth basin further into the hangingwall. Continued fold growth led to truncation of strata near the crest of the fault‐bend anticline before shortcut faulting bypassed the ramp‐flat‐ramp structure and ended folding. Accommodation along the fault‐parallel axis is controlled by the transverse folds, the location and size of which depends on the degree of linkage in the fault network and the accumulated displacement on causal faults. We construct transverse fold trajectories by tracing transverse fold hinges through space and time to highlight the positions of maximum and minimum accommodation and potential sediment entry points to hangingwall growth basins. The length and shape of the constructed trajectories relate to the displacement on their parent faults, duration of fault activity, timing of transverse basin infill, fault linkage and strain localization. We emphasize that the considerable wavelength, amplitudes and potential periclinal geometry of extensional folds make them viable targets for CO₂ storage or hydrocarbon exploration in rift basins.     

Geometry of growth strata in a transpressive fold belt in field and analogue model: Gosau Group at Muttekopf,Northern Calcareous Alps,Austria

The thrust sheets of the Northern Calcareous Alps were emplaced during Late Cretaceous thrust‐dominated transpression expressed by thrust sheets segmented by closely spaced tear faults. Thrust sheet‐top sediments were deposited during thrusting and associated fold growth and were controlled by active folding and tearing. We observe two types of angular unconformities: (1) Angular unconformities above folds between tear faults conform with the model of progressive unconformities. Across these unconformities dip decreases upsection. (2) Here, we define progressive unconformities that are related to tear faults and are controlled by both folding and tearing. Across these unconformities both strike and dip change. In growth strata overlying folds dissected by high‐angle faults, such unconformities are expected to be common. We used analogue modelling to define the geometry of the tear faults and related unconformities. Within the syn‐tectonic sediments, a steep, upward flattening thrust within a broader, roughly tulip‐shaped drag zone develops. The thrust roots in the tear fault in pre‐tectonic deposits and is curved upward toward the downthrown block. Vertical offset on the thrust is related to differential vertical uplift caused by, for example, growth of folds with different wavelength and amplitude on either side of the tear fault. Formation of progressive unconformities is governed by the relationship between the rates of deposition and vertical growth of a structure. Fault‐related progressive unconformities are additionally controlled by the growth of the vertical step across the tear fault. When the rates of vertical growth of two neighbouring folds separated by a tear fault are similar, the rate of growth across the tear fault is small; if the first differ, the latter is high. Episodic tear fault activity may create several angular unconformities attached to a tear fault or allow the generation of angular unconformities near tear faults in sedimentary systems that have a rate of deposition too high to generate classical progressive unconformities between the tear faults.     

Influence of the mechanical behaviour of brittle–ductile fold–thrust belts on the development of foreland basins

A two‐dimensional mathematical model considering coupling between a deforming elasto‐visco‐plastic fold–thrust belt, flexural subsidence and diffusional surface processes is solved using the Finite Element Method to investigate how the mechanical behaviour of brittle–ductile wedges influences the development of foreland basins. Results show that, depending mainly on the strength of the basal décollement, two end‐member types of foreland basin are possible. When the basal detachment is relatively strong, the foreland basin system is characterised by: (1) Highly asymmetrical orogen formed by thrusts concentrated in the incoming pro‐wedge. (2) Sedimentation on retro‐side takes place in one major foredeep basin which grows throughout orogen evolution. (3) Deposition on the pro‐side occurs initially in the foredeep, and continues in the wedge‐top before isolated basins are advected towards the orogen core where they become uplifted and exhumed. (4) Most pro‐wedge basins show an upward progression from low altitude, foredeep deposits at the base to high altitude, wedge‐top deposits near the surface. In contrast, when the basal detachment behaves weakly due to the presence of low viscosity material such as salt, the foreland basin system is characterised by (1) Broad, low relief orogen showing little preferential vergence and predominance of folding relative to faulting. (2) Deposition mainly in wedge‐top basins showing growth strata. (3) Many basins are initiated contemporaneously but form discontinuously due to the locus of active deformation jumping back and forth between different structures. Model results successfully reproduce first order observations of deforming brittle–ductile wedges and foreland basins. Moreover, the results support and provide a framework for understanding the existence of two main end‐member foreland basin types, simple and complex, associated with fold–thrust belts whose detachments are relatively strong and weak, respectively.     

Basin modelling of the Limón Back-arc Basin (Costa Rica): burial history and temperature evolution of an island arc-related basin-system

The Limón back‐arc basin belongs to the southern Central American arc‐trench system and is situated at the east coast of Costa Rica. The basin‐fill consists of Late Cretaceous to Pleistocene sedimentary rocks. A northern and a southern sub‐basin can be defined, separated by the E–W‐trending Trans Isthmic Fault System. The North Limón Basin is nearly undeformed, whereas the South Limón Basin is characterized by a fold‐and‐thrust belt. Both sub‐basins have a very similar sedimentary fill and can act as a natural laboratory for reconstructing controlling factors of arc‐related sedimentary basins as well as the influence of deformation on a basin system. Modelling focused on burial history and temperature evolution. Two‐dimensional simulations were carried out with the software PetroMod^®. The geohistory curve of the North Limón Basin is overall linear, indicating continuous subsidence. The South Limón Basin is also characterized by continuous subsidence, but rates strongly increased at the beginning of the Neogene. Despite a rapid Plio‐Pleistocene deformation of the fold‐and‐thrust belt, the present‐day temperature field is not disturbed in that area. The modelling results indicate a mean heat flow of 60 mW m^?2 for the North Limón Basin and 41 mW m^?2 for the South Limón Basin. These values are low compared with other back‐arc basins. The lower values are attributed to the following effects: (1) underlying basaltic crust, (2) the lack of an initial rift phase, (3) the low extension rates, (4) absence of volcanic activity and (5) insulation effects of a thick sediment pile. The reasons for the locally lower heat flow in the southern sub‐basin can be found in the low‐angle subduction of the Cocos Ridge. Owing to the low subduction angle, the cool fore‐arc mantle‐wedge below the island‐arc is pushed backwards increasing the cooled area.     

Stratigraphic expression of the lateral propagation and growth of isolated fault-related uplifts

The lateral propagation of faults and folds is known to be an important process during the development of mountain belts, but little is known about the manner in which along‐strike fault–fold growth is expressed in pre‐ and syntectonic (growth) strata. We use a coupled tectonic and stratigraphic model to investigate the along‐strike stratigraphic expression of fault‐related folds/uplifts that grow in both the transport and strike directions. We consider faults that propagate following a quadratic (nonself‐similar evolution) or linear (self‐similar evolution) scaling law, using different slip distributions per episode of fault propagation, under general background sedimentation. We find that the long‐strike geometry of pre‐ and syntectonic strata and the geometry of growth axial surfaces reflect the mode of fault propagation. The geometry of strata observed in the model is similar to that observed in natural contractional structures when: (1) the evolution of the fault is nonself‐similar, or (2) the fault grows as a result of thrust faulting events with similar displacements along strike that are terminated abruptly at the fault tips.     

Closing and continentalization of the South Pyrenean foreland basin (NE Spain): magnetochronological constraints

This paper presents new magnetostratigraphic results from a 1100‐m‐thick composite section across the marine to continental sediments of the central part of the SE margin of the Ebro basin (NE Spain). Integration with existing marine and continental biochronological data allows a robust correlation with the geomagnetic polarity time scale. The resulting absolute chronology ranges from 36.3 to 31.1 Ma (Priabonian to Rupelian), and yields an interpolated age of ～36.0 Ma (within chron C16n.2n) for the youngest marine sediments of the eastern Ebro basin. This age is in concordance with a reinterpretation of earlier magnetostratigraphic data from the western South Pyrenean foreland basin, and indicates that continentalization of the basin occurred as a rapid and isochronous event. The basin continentalization, determined by the seaway closure that resulted from the uplift of the western Pyrenees, was probably coincident with a mid‐amplitude eustatic sea level low with a maximum at 36.2 Ma. The base level drop that followed the basin closure and desiccation does not appear associated to a significant sedimentary hiatus along the margins, suggesting a late Eocene shallow marine basin that rapidly refilled and raised its base level after the seaway closing. Rapid basin filling following continentalization predates the phase of rapid exhumation of the Central Pyrenean Axial Zone from 35.0 to 32.0 Ma, determined from the thermochronology data. It is possible then that sediment aggradation at the front of the fold‐and‐thrust belt could have contributed to a decrease in the taper angle, triggering growth of the inner orogenic wedge through break‐back thrusting and underplating. Contrasting sedimentation trends between the western and eastern sectors of the South Pyrenean foreland indicate that basin closing preferentially affected those areas subjected to sediment bypass towards the ocean domain. As a result, sediment ponding after basin closure is responsible for a two‐fold increase of sedimentation rates in the western sector, while changes of sedimentation rates are undetected in the more restricted scenario of the eastern Ebro basin.     

Drainage evolution in response to fold growth in the hanging-wall of the Khazar fault, north-eastern Alborz, Iran

The Khazar fault is an active thrust fault in the northern part of the Alborz Range, which is associated with folding (the Khazar anticline) in its hanging‐wall. Regional geological studies indicate activity of the fault in Cenozoic time, and active propagation of the fault‐related folding towards the west. The Neka river drainage basin, which is a longitudinal river flowing mostly in the backlimb of the Khazar anticline, shows evidence for active folding and faulting influencing drainage evolution. Observations are made in different parts of the Neka river course, according to which a new morphotectonic feature is introduced within the river basin, termed as ‘tilted reach’. This feature is considered as a result of river course tilting in the backlimb of the growing fold, diversion, and capture of the river by other rivers. Consecutive episodes of similar events would have resulted in the development of a long drainage basin parallel to the growing fold structure.     

Structure and evolution of mass transport deposits in the South Caspian Basin,Azerbaijan

The Quaternary to late Pliocene sedimentary succession along the margin of the South Caspian Basin contains numerous kilometre‐scale submarine slope failures, which were sourced along the basin slope and from the inclined flanks of contemporaneous anticlines. This study uses three‐dimensional (3D) seismic reflection data to visualise the internal structure of 27 mass transport deposits and catalogues the syndepositional structures contained within them. These are used to interpret emplacement processes occurring during submarine slope failure. The deposits consist of three linked structural domains: extensional, translational and compressive, each containing characteristic structures. Novel features are present within the mass transport deposits: (1) a diverging retrogression of the headwall scarp; (2) the absence of a conventional headwall scarp around growth stratal pinch outs; (3) restraining bends in the lateral margin; (4) a downslope increase in the throw of thrust faults. The results of this study shed light on the deformation that occurred during submarine slope failure, and highlight an important geological process in the evolution of the South Caspian Basin margin.     

Three-dimensional restoration of original sedimentary geometries in deformed basin fill, onshore Brunei Darussalam, NW Borneo

The western flank of the Jerudong Anticline, onshore Brunei Darussalam, provides a rare opportunity to analyse the base of a major Miocene mud‐rich delta in outcrop, including kilometre‐scale prograding clinoforms, delta‐front turbidites and large‐scale syndepositional faults. The lateral continuation of this system in the subsurface of the Belait Syncline is documented on two‐dimensional (2D) reflection seismic data and wireline logs. In order to link geological observations at surface with corresponding geophysical subsurface signatures, we constructed a combined, quantitative 3D surface–subsurface model of onshore Brunei Darussalam. This 3D model is used to analyse and discuss the relation between field geology and geophysical subsurface interpretation, and provides the base for a quantitative kinematic restoration of the Miocene Belait delta to its original shape before folding. Final decompaction of the balanced rock volume allows the reconstruction of the palaeo‐relief of a series of Miocene clinoforms, indicating a close relation between delta‐lobe activity, clinoform morphology and the generation of delta‐front turbidites.     

Evaluating foreland basin partitioning in the northern Andes using Cenozoic fill of the Floresta basin,Eastern Cordillera,Colombia

This paper addresses foreland basin fragmentation through integrated detrital zircon U–Pb geochronology, sandstone petrography, facies analysis and palaeocurrent measurements from a Mesozoic–Cenozoic clastic succession preserved in the northern Andean retroarc fold‐thrust belt. Situated along the axis of the Eastern Cordillera of Colombia, the Floresta basin first received sediment from the eastern craton (Guyana shield) in the Cretaceous–early Palaeocene and then from the western magmatic arc (Central Cordillera) starting in the mid‐Palaeocene. The upper‐crustal magmatic arc was replaced by a metamorphic basement source in the middle Eocene. This, in turn, was replaced by an upper‐crustal fold‐thrust belt source in the late Eocene which persisted until Oligocene truncation of the Cenozoic section by the eastward advancing thrust front. Sedimentary facies analysis indicates minimal changes in depositional environments from shallow marine to low‐gradient fluvial and estuarine deposits. These same environments are recorded in coeval strata across the Eastern Cordillera. Throughout the Palaeogene, palaeocurrent and sediment provenance data point to a uniform western or southwestern sediment source. These data show that the Floresta basin existed as part of a laterally extensive, unbroken foreland basin connected with the proximal western (Magdalena Valley) basin from mid‐Paleocene to late Eocene time when it was isolated by uplift of the western flank of the Eastern Cordillera. The Floresta basin was also connected with the distal eastern (Llanos) basin from the Cretaceous until its late Oligocene truncation by the advancing thrust front.     

Alluvial fan development and morpho-tectonic evolution in response to contractional fault reactivation (Late Cretaceous–Palaeocene), Provence, France

Along-strike variability within a Late Cretaceous to early Palaeocene contractional growth structure and associated alluvial fan deposits is documented at the northern margin of the Arc Basin (Provence, SE France). This contribution shows that alluvial fans can be used as high-resolution proxies to reconstruct structural segmentation and palaeo-geomorphological evolution of a source/basin margin system. Facies-based reconstruction allows the spatial and temporal distribution of alluvial fan bodies to be mapped. Relationships between fan area and catchment size from modern alluvial fan systems were used to estimate palaeo-catchment size. Combining alluvial fan morphologies with catchment area, pebble provenance analysis and growth structure reconstruction, we show that: (1) fan distribution and related depositional processes were strongly influenced by intrinsic parameters such as drainage basin evolution, local structural inheritance and lateral facies changes in source area lithologies; (2) Inherited structures trending N100 effectively controlled the first-order location of the fold and thrust structures (Montagne Sainte-Victoire Range) and adjacent depositional areas (Arc Basin); (3) Syn-sedimentary faults trending N010-030 influenced the source/basin margin development and interacted with developing growth structures; (4) Facies changes in Jurassic carbonates controlled fold development and consequently the structural evolution of the source area; and (5) the N010-030 faults and along-strike variability of the source/basin margin system were ultimately controlled by basement structures that controlled where Late Cretaceous deformation nucleated. The overall architecture of the source/basin margin system reflects segmentation and strain partitioning along strike, as demonstrated by diachronous alluvial fan distribution.     

Stratigraphic and structural expression of the lateral growth of thrust fault-propagation folds: results and implications from kinematic modelling

In order to better understand the development of thrust fault‐related folds, a 3D forward numerical model has been developed to investigate the effects that lateral slip distribution and propagation rate have on the fold geometry of pre‐ and syn‐tectonic strata. We consider a fault‐propagation fold in which the fault propagates upwards from a basal decollement and along‐strike normal to transport direction. Over a 1 Ma runtime, the fault reaches a maximum length of 10 km and accumulates a maximum displacement of 1 km. Deformation ahead of the propagating fault tip is modelled using trishear kinematics while backlimb deformation is modelled using kink‐band migration. The applicability of two different lateral slip distributions, namely linear‐taper and block‐taper, are firstly tested using a constant lateral propagation rate. A block‐taper slip distribution replicates the geometry of natural fold‐thrusts better and is then used to test the sensitivity of thrust‐fold morphology to varied propagation rates in a set of fault‐propagation folds that have identical final displacement to length (D_max/L_max) ratios. Two stratigraphic settings are considered: a model in which background sedimentation rates are high and no topography develops, and a model in which a topographic high develops above the growing fold and local erosion, transport and deposition occur. If the lateral propagation rate is rapid (or geologically instantaneous), the fault tips quickly become pinned as the fault reaches its maximum lateral extent (10 km), after which displacement accumulates. In both stratigraphic settings, this leads to strike‐parallel rotation of the syn‐tectonic strata near the fault tips; high sedimentation rates relative to rates of uplift result in along‐strike thinning over the structural high, while low sedimentation rates result in pinchout against it. In contrast, slower lateral propagation rates (i.e. up to one order of magnitude greater than slip rate) lead to the development of along‐strike growth triangles when sedimentation rates are high, whereas when sedimentation rates are low, offflap geometries result. Overall we find that the most rapid lateral propagation rates produce the most realistic geometries. In both settings, time‐equivalent units display both nongrowth and growth stratal geometries along‐strike and the transition from growth to nongrowth has the potential to delineate the time of fault/fold growth at a given location. This work highlights the importance of lateral fault‐propagation and fault tip pinning on fault and fold growth in three dimensions and the complex syn‐tectonic geometries that can result.     

Modern and ancient fluvial megafans in the foreland basin system of the central Andes, southern Bolivia: implications for drainage network evolution in fold-thrust belts

ABSTRACT Fluvial megafans chronicle the evolution of large mountainous drainage networks, providing a record of erosional denudation in adjacent mountain belts. An actualistic investigation of the development of fluvial megafans is presented here by comparing active fluvial megafans in the proximal foreland basin of the central Andes to Tertiary foreland‐basin deposits exposed in the interior of the mountain belt. Modern fluvial megafans of the Chaco Plain of southern Bolivia are large (5800–22 600 km²), fan‐shaped masses of dominantly sand and mud deposited by major transverse rivers (Rio Grande, Rio Parapeti, and Rio Pilcomayo) emanating from the central Andes. The rivers exit the mountain belt and debouch onto the low‐relief Chaco Plain at fixed points along the mountain front. On each fluvial megafan, the presently active channel is straight in plan view and dominated by deposition of mid‐channel and bank‐attached sand bars. Overbank areas are characterized by crevasse‐splay and paludal deposition with minor soil development. However, overbank areas also contain numerous relicts of recently abandoned divergent channels, suggesting a long‐term distributary drainage pattern and frequent channel avulsions. The position of the primary channel on each megafan is highly unstable over short time scales. Fluvial megafans of the Chaco Plain provide a modern analogue for a coarsening‐upward, > 2‐km‐thick succession of Tertiary strata exposed along the Camargo syncline in the Eastern Cordillera of the central Andean fold‐thrust belt, about 200 km west of the modern megafans. Lithofacies of the mid‐Tertiary Camargo Formation include: (1) large channel and small channel deposits interpreted, respectively, as the main river stem on the proximal megafan and distributary channels on the distal megafan; and (2) crevasse‐splay, paludal and palaeosol deposits attributed to sedimentation in overbank areas. A reversal in palaeocurrents in the lowermost Camargo succession and an overall upward coarsening and thickening trend are best explained by progradation of a fluvial megafan during eastward advance of the fold‐thrust belt. In addition, the present‐day drainage network in this area of the Eastern Cordillera is focused into a single outlet point that coincides with the location of the coarsest and thickest strata of the Camargo succession. Thus, the modern drainage network may be inherited from an ancestral mid‐Tertiary drainage network. Persistence and expansion of Andean drainage networks provides the basis for a geometric model of the evolution of drainage networks in advancing fold‐thrust belts and the origin and development of fluvial megafans. The model suggests that fluvial megafans may only develop once a drainage network has reached a particular size, roughly 10⁴ km²– a value based on a review of active fluvial megafans that would be affected by the tectonic, climatic and geomorphologic processes operating in a given mountain belt. Furthermore, once a drainage network has achieved this critical size, the river may have sufficient stream power to prove relatively insensitive to possible geometric changes imparted by growing frontal structures in the fold‐thrust belt.