Detection of fluid migration pathways in seismic data: implications for fault seal analysis

A new and efficient method for fault seal analysis using seismic data is presented. It uses multiple seismic attributes and neural networks to enhance fluid migration pathways, including subtle features that are not detectable using single attributes only. The method may be used as a first estimate of fault seal or to calibrate results from other techniques. The results provide information about which faults and fault segments are sealing or leaking. Fluid flow along individual faults appears to be focused along zones of weakness, and fault seal research should thus be focused on finding such weak locations within fault zones, a task that is best done using three‐dimensional (3D) seismic data. Under certain conditions, it is suggested that fluids migrate along fault planes by a diapiric fluid flow mechanism. The results assist in calibrating the bulk hydraulic properties of faults and rock formations and can be used in basin modelling.     

Characterization of the depocenters and the basement structure,below the central Chile Andean Forearc: A 3D geophysical modelling in Santiago Basin area

Since the last century, several geological and geophysical studies have been developed in the Santiago Basin to understand its morphology and tectonic evolution. However, some uncertainties regarding sedimentary fill properties and possible density anomalies below the sediments/basement boundary remain. Considering that this is an area densely populated with more than 6 million inhabitants in a highly active seismotectonic environment, the physical properties of the Santiago Basin are important to study the geological and structural evolution of the Andean forearc and to characterize its seismic response and related seismic hazard. Two and three‐dimensional gravimetric models were developed, based on a database of 797 compiled and 883 newly acquired gravity stations. To produce a well‐constrained basement elevation model, a review of 499 wells and 30 transient electromagnetic soundings were used, which contribute with basement depth or minimum sedimentary thickness information. For the 2‐D modelling, a total of 49 gravimetric profiles were processed considering a homogeneous density contrast and independent regional trends. A strong positive gravity anomaly was observed in the centre of the basin, which complicated the modelling process but was carefully addressed with the available constrains. The resulting basement elevation models show complex basement geometry with, at least, eight recognizable depocenters with maximum sedimentary infill of ~ 500 m. The 3‐D density models show alignments in the basement that correlates well with important intrusive units of the Cenozoic and Mesozoic. Along with interpreted fault zones westwards and eastwards of the basin, the observations suggest a structural control of Santiago basin geometry, where recent deformation associated with the Andean contractional deformation front and old structures developed during the Cenozoic extension are superimposed to the variability of river erosion/deposition processes.     

Three-dimensional-seismic coherency signature of Niger Delta growth faults: integrating sedimentology and tectonics

This case study of growth faults and associated deltaic sedimentation in the shallow‐offshore Niger Delta uses an integrated analysis of three‐dimensional (3D)‐seismic coherence facies and wireline data that supports an evaluation of the sedimentary response to delta tectonics. The study area comprises four fault blocks bounded by a set of kilometre‐scale, basinward‐dipping, synsedimentary normal faults. Correlation of highly variable growth stratigraphy across faults was achieved by a systematic visualization and interpretation of series of coherence horizon‐slices: the detection and matching of erosive and depositional patterns (e.g. channels) across faults allowed the establishment of sedimentology‐controlled links between diverse footwall and hanging‐wall growth successions. At the same time, this interpretation approach helped to visualize seismic‐sedimentological and seismic‐geomorphological features survey‐wide at all depth levels. The integration of this extensive 3D database with lithology information from wireline logs provides a powerful tool for subsurface sedimentology interpretation. Synoptic analysis of the 3D‐seismic sedimentology interpretation with stratigraphy based fault‐kinematic analysis using throw vs. depth plots (Th–Z plots) enabled a discussion of the relation between delta tectonics and sedimentary‐system development, and the evaluation of the Th–Z method for subsurface‐lithology prediction. The interpretation results document that both motion analysis of synsedimentary deltaic faults and Th–Z‐based lithology prediction are only feasible when supported by detailed 3D information on palaeoenvironment and palaeotopography at and around studied fault systems. We therefore recommend the use of fast‐track fault‐kinematic and subsurface‐lithology predictions based on Th–Z plots only when supported by comprehensive 3D seismic‐sedimentological interpretations.     

Geological structure and kinematics of normal faults in the Otway Basin,Australia, based on quantitative analysis of 3‐D seismic reflection data

The Otway Basin in the south of Victoria, Australia underwent three phases of deformation during breakup of the southern Australian margin. We assess the geometry and kinematics of faulting in the basin by analysing a 3‐D reflection seismic volume. Eight stratigraphic horizons and 24 SW‐dipping normal faults as well as subordinate antithetic faults were interpreted. This resulted in a high‐resolution geological 3‐D model (ca. 8 km × 7 km × 4 km depth) that we present as a supplementary 3‐D PDF (Data S1). We identified hard‐ and soft‐linking fault connections over the entire area, such as antithetic faults and relay ramps, respectively. Most major faults were continuously active from Early to Late Cretaceous, with two faults in the northern part of the study area active until at least the Oligocene. Allan maps of faults show tectonic activity continuously waned over this time period. Isopach maps of stratigraphic volumes quantify the amount of syn‐sedimentary movement that is characteristic of passive margins, such as the Otway Basin. We show that the faults possess strong corrugations (with amplitudes above the seismic resolution), which we illustrated by novel techniques, such as cylindricity and curvature. We argue that the corrugations are produced by sutures between sub‐vertical fault segments and this morphology was maintained during fault growth. Thus, they can be used to indicate the kinematics vector of the fault movement. This evidences, together with left‐stepping relay ramps, that 40% of the faults had a small component (up to 25°) of dextral oblique slip as well as normal (dip‐slip) movement.     

Fault‐controlled fluid flow inferred from hydrothermal vents imaged in 3D seismic reflection data,offshore NW Australia

Fluid migration pathways in the subsurface are heavily influenced by pre‐existing faults. Although studies of active fluid‐escape structures can provide insights into the relationships between faults and fluid flow, they cannot fully constrain the geometry of and controls on the contemporaneous subsurface fluid flow pathways. We use 3D seismic reflection data from offshore NW Australia to map 121 ancient hydrothermal vents, likely related to magmatic activity, and a normal fault array considered to form fluid pathways. The buried vents consist of craters up to 264 m deep, which host a mound of disaggregated sedimentary material up to 518 m thick. There is a correlation between vent alignment and underlying fault traces. Seismic‐stratigraphic observations and fault kinematic analyses reveal that the vents were emplaced on an intra‐Tithonian seabed in response to the explosive release of fluids hosted within the fault array. We speculate that during the Late Jurassic the convex‐upwards morphology of the upper tip‐lines of individual faults acted to channelize ascending fluids and control where fluid expulsion and vent formation occurred. This contribution highlights the usefulness of 3D seismic reflection data to constraining normal fault‐controlled subsurface fluid flow.     

Geometry of half-grabens containing a mid-level viscous décollement

In this work, we explore by means of analogue models how different basin-bounding fault geometries and thickness of a viscous layer within the otherwise brittle pre-rift sequence influence the deformation and sedimentary patterns of basins related to extension. The experimental device consists of a rigid wooden basement in the footwall to simulate a listric fault. The hangingwall consists of a sequence of pre-rift deposits, including the shallow interlayered viscous layer, and a syn-rift sequence deposited at constant intervals during extension. Two different geometries exist of listric normal faults, dip at 30 and 60° at surface. This imposes different geometries in the hangingwall anticlines and their associated sedimentary basins. A strong contrast exists between models with and without a viscous layer. With a viscous décollement, areas near the main basement fault show a wide normal drag and the hangingwall basin is gently synclinal, with dips in the fault side progressively shallowing upwards. A secondary roll-over structure appears in some of the models. Other structures are: (1) reverse faults dipping steeply towards the main fault, (2) antithetic faults in the footwall, appearing only in models with the 30° dipping fault and silicone-level thicknesses of 1 and 1.5 cm and (3) listric normal faults linked to the termination of the detachment level opposite to the main fault, with significant thickness changes in the syn-tectonic units. The experiments demonstrate the importance of detachment level in conditioning the geometry of extensional sedimentary basins and the possibility of syncline basin geometries associated with a main basement fault. Comparison with several basins with half-graben geometries containing a mid-level décollement supports the experimental results and constrains their interpretation.     

Listric extensional fault systems - results of analogue model experiments

Abstract Analogue models are a powerful tool for investigating progressive deformation in extensional fault systems. This paper presents exciting new insights into the progressive evolution of hanging wall structures in listric extensional terranes. Analogue models, scaled to simulate deformation in a sedimentary sequence, were constructed for simple listric and ramp/flat listric extensional detachments. For each detachment geometry homogeneous sand, sand/mica and sand/clay models were used to simulate respectively, deformation of isotropic sediments, of anisotropic sediments and of sedimentary sequences with competency contrasts. Roll-over anticlines with geometrically necessary crestal collapse graben structures are characteristic of the steepening-upwards segments of listric extensional fault systems in all of our models. With progressive deformation, crestal collapse grabens show hanging wall nucleation of new faults. Variations in graben size, amount of fault rotation and throw, are dependent on detachment curvature and amount of extension. Individual faults and associated fault blocks may significantly change shape during extension. Complex and apparently conjugate fault arrays are the result of superposition of successive crestal collapse grabens. Ramp/flat listric extensional fault systems are characterized by a roll-over anticline and a crestal collapse graben system associated with each steepening-upwards segment of the detachment and a ramp zone consisting of a hanging wall syncline and a complex deformation zone with local reverse faults. The roll-over anticlines and crestal collapse graben are similar in geometry to those formed in simple listric extensional systems. The models demonstrate that the geometry of the detachments exerts a fundamental control on the evolution of hanging wall structures. Analysis of particle displacement paths for these experiments provides new insights into the mechanical development of roll-over anticlines. Two general models for deformation above simple listric and ramp/flat listric extensional detachments have been erected.     

Growth‐faults from delta collapse – structural and sedimentological investigation of the Last Chance delta,Ferron Sandstone,Utah

Investigations of syn‐sedimentary growth faults in the Last Chance delta (Ferron Sandstone, Utah, USA) show that fault‐bounded half‐grabens arrested high amounts of sand in the mouth bar and/or distributary channel areas. Fault‐controlled morphology causes changes in routing of the delta top to delta front drainage towards the long axis of half‐grabens. Faulting was spatially and temporally non‐systematic, and polyphase, with 3D cusp/listric fault geometries instigated by linkage of variously oriented segments. Hanging wall rollover folds consisting of wedge‐shaped syn‐kinematic sand attest to rapid <1‐m slip increments on faults followed by mild erosion along crests of fault blocks and sedimentary infill of adjacent accommodation. Triangle‐zones in prodelta to delta front muds are located underneath steeper faults and interconnected rotated fault‐flats. Their geometry is that of antiformal stack duplexes, in an arrangement of low‐angle‐to‐bedding normal faults at the base, replaced by folded thrusts upwards. These faults show a brittle, frictional flow deformation mechanism ascribed to early compaction of mud. For syn‐kinematic sand, there is a change from general granular/hydroplastic flow in shear zones to later brittle failure and cataclasis, a transition instigated by precipitation of calcite cement. Extensional faulting in the Last Chance delta was likely controlled by gravity driven collapse towards the delta slope and prodelta, as is commonly observed in collapsing deltas. The trigger and driving mechanism is envisioned as localized loads from sand deposited within distributary channels/mouth bars and fault‐controlled basins along the delta top. A regional tilt and especially displacement of compacted mud below sand bodies towards less compacted muds also contributed to the faulting.     

The effect of fault relay and clay smearing on groundwater flow patterns in the Lower Rhine Embayment

Faults strongly impact groundwater flow in the unconsolidated sediments of the Lower Rhine Embayment. Hydraulic head maps show that many individual faults form a barrier to fluid flow whereas relay structures in these faults are sites of hydraulic contact between otherwise separated aquifers. The fluid flow patterns around the Rurrand Fault close to the largest open‐pit mine in the Lower Rhine Embayment is one of the first well‐documented examples of fluid flow around a fault relay zone. The effect of clay smearing could be quantified using the Shale Gouge Ratio (SGR) method that is common in hydrocarbon‐related studies but has not been applied to groundwater flow data so far. The effect of fault relay zones on groundwater flow is analysed using numerical simulations. It is concluded that fault relay needs special consideration in the evaluation of the sealing capacities of faults in sedimentary basins. Moreover, it is demonstrated that the SGR methodology is a promising tool for the estimation of fault zone hydraulic properties in hydrogeological modelling.     

Structure of the footwall of a listric fault system revealed by 3D seismic data from the Niger Delta

We use three‐dimensional (3D) seismic reflection data to analyse the architecture of the footwall of a listric fault, in a gravitationally driven extensional system, in the north‐western Niger Delta. In contrast to conventional listric normal fault models with a single master listric fault plane the level of detachment switches from a deeper to shallower level. The footwall evolves through the generation of new master detachment faults and detachments, which transfers hanging wall rocks into the footwall. New detachments form by branching off pre‐existing detachment levels, cutting‐up through stratigraphy to the next mechanical weakness, separating discrete sections of extended strata. As a consequence a deeper, older array of seaward‐dipping, tilted extensional fault blocks is now located in the footwall beneath the master listric detachment fault. The structural complexity located below the master detachment fault highlights extensional episodes on separate detachment faults that are not captured in conventional listric models. We speculate that changes in the level of the detachment are caused by mechanical weaknesses controlled by lithology, pore pressure and episodes of sediment loading related to deltaic progradation.     

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.     

The physical scale modelling of braided alluvial architecture and estimation of subsurface permeability

The quantitative modelling of fluvial reservoirs, especially in the stages of enhanced oil recovery, requires detailed three‐dimensional data at both the scale of the channel belt and within‐channel. Although studies from core, analogue outcrop and modern environments may partially meet these needs, they often cannot provide detail on the smaller‐scale (i.e. channel‐scale) heterogeneity, frequently suffer from limited three‐dimensional exposure and cannot be used to examine the influence of different variables on the process–deposit relationship. Physical modelling offers a complementary technique that can address many of these quantitative requirements and holds great future potential for integration with reservoir modelling. Physical modelling provides the potential to upscale results and derive reservoir information on three‐dimensional facies geometry, connectivity and permeability. This paper describes the development and use of physical modelling, which employs generic Froude‐scaling principles, in an experimental basin that permits aggradation in order to model the morphology and subsurface depositional stratigraphy of coarse‐grained braided rivers. An example is presented of a 1:50 scale model based on the braided Ashburton River, Canterbury Plains, New Zealand and the adjacent late Quaternary braided alluvium exposed in the coastal cliffs. Critically, a full, bimodal grain size distribution (20% sand and 80% gravel) was used to replicate the prototype, which allows the realistic reproduction of the surface morphology and importantly permits grain size sorting during deposition. Uncertainties associated with the compression of time, sediment mass balance and the hydrodynamics of the finest particle sizes do not appear to affect the reproducibility of stratigraphy between experimental and natural environments. Sectioning of the preserved sedimentary sequence in the physical model allows quantification of the geometry, shape, spatial distribution and internal sedimentary structure of the coarse‐ and fine‐grained facies. A six‐fold facies scheme is proposed for the model braided alluvium and a direct link is established between the grain size distribution and facies type: this allows permeability to be estimated for each facies, which can be mapped onto two‐dimensional vertical cross‐sections of the preserved stratigraphy. Results demonstrate the dominance of four facies based on permeability that range over three orders of magnitude in hydraulic conductivity. Quantification of such variability, and linkage to both vertical proportion curves for facies distribution and connectivity presents significant advantages over other methodologies and offers great potential for the modelling of heterogeneous braided river sediments at the within channel‐belt scale. This paper outlines how physical models may be used to develop high‐resolution, geologically‐accurate, object‐based reservoir simulation models.     

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.     

Dynamic modelling of passive margin salt tectonics: effects of water loading, sediment properties and sedimentation patterns

We investigate the evolution of passive continental margin sedimentary basins that contain salt through two‐dimensional (2D) analytical failure analysis and plane‐strain finite‐element modelling. We expand an earlier analytical failure analysis of a sedimentary basin/salt system at a passive continental margin to include the effects of submarine water loading and pore fluid pressure. Seaward thinning sediments above a weak salt layer produce a pressure gradient that induces Poiseuille flow in the viscous salt. We determine the circumstances under which failure at the head and toe of the frictional–plastic sediment wedge occurs, resulting in translation of the wedge, landward extension and seaward contraction, accompanied by Couette flow in the underlying salt. The effects of water: (i) increase solid and fluid pressures in the sediments; (ii) reduce the head to toe differential pressure in the salt and (iii) act as a buttress to oppose failure and translation of the sediment wedge. The magnitude of the translation velocity upon failure is reduced by the effects of water. The subsequent deformation is investigated using a 2D finite‐element model that includes the effects of the submarine setting and hydrostatic pore pressures. The model quantitatively simulates a 2D approximation of the evolution of natural sedimentary basins on continental margins that are formed above salt. Sediment progradation above a viscous salt layer results in formation of landward extensional basins and listric normal growth faults as well as seaward contraction. At a later stage, an allochthonous salt nappe overthrusts the autochthonous limit of the salt. The nature and distribution of major structures depends on the sediment properties and the sedimentation pattern. Strain weakening of sediment favours landward listric growth faults with formation of asymmetric extensional depocentres. Episodes of low sediment influx, with partial infill of depocentres, produce local pressure gradients in the salt that result in diapirism. Diapirs grow passively during sediment aggradation.     

Multiscale structure in sedimentary basins

Hierarchies of superimposed structures are found in maps of geological horizons in sedimentary basins. Mapping based on three‐dimensional (3D) seismic data includes structures that range in scale from tens of metres to hundreds of kilometres. Extraction of structures from these maps without a priori knowledge of scale and shape is analogous to pattern recognition problems that have been widely researched in disciplines outside of Geoscience. A number of these lessons are integrated and applied within a geological context here. We describe a method for generating multiscale representations from two‐dimensional sections and 3D surfaces, and illustrate how superimposed geological structures can be topologically analysed. Multiscale analysis is done in two stages – generation of scale‐space as a geometrical attribute, followed by identification of significant scale‐space objects. Results indicate that Gaussian filtering is a more robust method than conventional moving average filtering for deriving multiscale geological structure. We introduce the concept of natural scales for identifying the most significant scales in a geological cross section. In three dimensions, scale‐dependent structures are identified via an analogous process as discrete topological entities within a four‐dimensional scale‐space cube. Motivation for this work is to take advantage of the completeness of seismic data coverage to see ‘beyond the outcrop’ and yield multiscale geological structure. Applications include identifying artefacts, scale‐specific features and large‐scale structural domains, facilitating multiscale structural attribute mapping for reservoir characterisation, and a novel approach to fold structure classification.     

Normal fault growth influenced by basement fabrics: The importance of preferential nucleation from pre‐existing structures

Reactivation of pre‐existing intra‐basement structures can influence the evolution of rift basins, yet the detailed kinematic relationship between these structures and overlying rift‐related faults remains poorly understood. Understanding the kinematic as well as geometric relationship between intra‐basement structures and rift‐related fault networks is important, with the extension direction in many rifted provinces typically thought to lie normal to fault strike. We here investigate this problem using a borehole‐constrained, 3D seismic reflection dataset from the Taranaki Basin, offshore New Zealand. Excellent imaging of intra‐basement structures and a relatively weakly deformed, stratigraphically simple sedimentary cover allow us to: (a) identify a range of interaction styles between intra‐basement structures and overlying, Plio‐Pleistocene rift‐related normal faults; and (b) examine the cover fault kinematics associated with each interaction style. Some of the normal faults parallel and are physically connected to intra‐basement reflections, which are interpreted as mylonitic reverse faults formed during Mesozoic subduction and basement terrane accretion. These geometric relationships indicate pre‐existing intra‐basement structures locally controlled the position and attitude of Plio‐Pleistocene rift‐related normal faults. However, through detailed 3D kinematic analysis of selected normal faults, we show that: (a) normal faults only nucleated above intra‐basement structures that experienced late Miocene compressional reactivation, (b) despite playing an important role during subsequent rifting, intra‐basement structures have not been significantly extensionally reactivated, and (c) preferential nucleation and propagation of normal faults within late Miocene reverse faults and folds appears to be the key genetic relationship between contractionally reactivated intra‐basement structures and rift‐related normal faults. Our analysis shows that km‐scale, intra‐basement structures can control the nucleation and development of newly formed, rift‐related normal faults, most likely due to a local perturbation of the regional stress field. Because of this, simply inverting fault strike for causal extension direction may be incorrect, especially in provinces where pre‐existing, intra‐basement structures occur. We also show that a detailed kinematic analysis is key to deciphering the temporal as well as simply the spatial or geometric relationship between structures developed at multiple structural levels.     

Tectonic controls on Cenozoic foreland basin development in the north‐eastern Andes,Colombia

In order to evaluate the relationship between thrust loading and sedimentary facies evolution, we analyse the progradation of fluvial coarse‐grained deposits in the retroarc foreland basin system of the northern Andes of Colombia. We compare the observed sedimentary facies distribution with the calculated one‐dimensional (1D) Eocene to Quaternary sediment‐accumulation rates in the Medina wedge‐top basin and with a three‐dimensional (3D) sedimentary budget based on the interpretation of ～1800 km of industry‐style seismic reflection profiles and borehole data. Age constraints are derived from a new chronostratigraphic framework based on extensive fossil palynological assemblages. The sedimentological data from the Medina Basin reveal rapid accumulation of fluvial and lacustrine sediments at rates of up to ～500 m my^?1 during the Miocene. Provenance data based on gravel petrography and paleocurrents reveal that these Miocene fluvial systems were sourced from Upper Cretaceous and Paleocene sedimentary units exposed to the west in the Eastern Cordillera. Peak sediment‐accumulation rates in the upper Carbonera Formation and the Guayabo Group occur during episodes of coarse‐grained facies progradation in the early and late Miocene proximal foredeep. We interpret this positive correlation between sediment accumulation and gravel deposition as the direct consequence of thrust activity along the Servitá–Lengupá faults. This contrasts with one class of models relating gravel progradation in more distal portions of foreland basin systems to episodes of tectonic quiescence.     

Novel discrete element modelling of Gilbert‐type delta formation in an active tectonic setting—first results

Gilbert deltas are now recognised as an important stratigraphic component of many extensional basins. They are remarkable due to their coarse‐grained nature, large size and steep foresets (up to 30–35°) and may exhibit a variety of slope instability features (faulting, slump scars, avalanching, etc.). They are also often closely related to major, basin‐margin normal faults. There has been considerable research interest in Gilbert deltas, partly due to their economic significance as stratigraphic traps for hydrocarbons but also due to their sensitivity to relative base level changes, giving them an important role in basin analysis. In addition to field studies, numerical modelling has also been used to simulate such deltas, with some success. However, until now, such studies have typically employed continuum numerical techniques where the basic data elements created by simulations are stratigraphic volumes or timelines and the sediments themselves have no internal properties per se and merely represent areas/volumes of introduced coarse‐grained, clastic and sedimentary material. Faulting or folding (if present) are imposed externally and do not develop (naturally) within the modelled delta body itself. Here, I present first results from a novel 2D numerical model which simulates coarse‐grained (Gilbert‐type) deltaic sedimentation in an active extensional tectonic setting undergoing a relative base level rise. Sediment is introduced as packages of discrete elements which are deposited beneath sea level, from the shoreline, upon a pre‐existing basin or delta. These elements are placed carefully and then allowed to settle onto the system. The elements representing the coarse‐grained, deltaic sediments can have an intrinsic coefficient of friction, cohesion or other material properties appropriate to the system being considered. The spatial resolution of the modelling is of the order of 15 m and topsets, foresets, bottomsets, faults, slumps and collapse structures all form naturally in the modelled system. Examples of deltas developing as a result of sediment supply from both the footwall and hanging‐wall of a normal fault, and subject to changes in fault slip rate are presented. Implications of the modelling approach, and its application and utility in basin research, are discussed.     

Depositional systems in multiphase rifts: seismic case study from the Lofoten margin,Norway

The evolution of depositional systems in multiphase rifts is influenced by the selective reactivation of faults between subsequent rift phases. The Middle Jurassic to Palaeocene tectonic history of the Lofoten margin, a segment of the North Atlantic rift system, is characterised by three distinct rift phases separated by long (>20 Myr) inter‐rift periods. The initial rift phase comprised a distinct fault initiation and linkage stage, whereas the later rift phases were characterised by selective reactivation of previously linked through‐going faults which resulted in immediate rift climax. Using 2‐D and 3‐D seismic reflection data in conjunction with shallow core data we present a 100 Myr record of shallow to deep marine depositional environments that includes deltaic clinoform packages, slope aprons and turbidite fans. The rapid re‐establishment of major faults during the later rift phases impacts on drainage systems and sediment supply. Firstly, the immediate localisation of strain and accumulation of displacement on few faults results in pronounced footwall uplift and possible fault block rotation along those faults, which makes it more likely for any antecedent fault‐transverse depositional systems to become reversed. Secondly, any antecedent axially‐sourced depositional systems that are inherited from the foregoing rift phase(s) are likely to be sustained after reactivation because such axial systems have already been directed around fault tips. Hence, the immediate localisation of strain through selective reactivation in the later rift phases restricts fault‐transverse sediment supply more than axial sediment supply, which is likely to be a key aspect of the tectono‐sedimentary evolution of multiphase rifts.