Basin Research - Volume 29, Issue 6, 2017

Multichannel high‐resolution seismic data along the northwestern margin of the Great Bahama Bank (GBB), Bahamas, detail the internal geometry and depositional history of a Neogene‐Quaternary carbonate slope‐to‐basin area. The stratigraphic architecture through this period evolves from (i) a mud‐dominated slope apron during the Miocene, (ii) a debris‐dominated base‐of‐slope apron during the Late Pliocene and then (iii) return to a slope apron with very short prograding clinoformal aprons during the Pleistocene. This geometric evolution was broadly constrained by the development of the Santaren Drift by bottom current since the Langhian. The drift expands along the northwestern GBB slope, forming a continuous correlative massive feature that shows successive phases of growth and retreat and influenced the downslope sediments distribution. Indeed, Late Pliocene deposits are confined into the moat, forming a strike‐continuous coarse debrites belt along the mid‐slope, preventing their free expansion into the basin. The occurrence of basinal drift that operated since 15 Ma showed a significant upslope growth around 3.6 Ma and is interpreted as resulting from the closure of the Central American Seaway which also coincides with a global oceanographic re‐organization and climate changes in the Northern Hemisphere.

Determining both short‐ and long‐term sedimentation rates is becoming increasingly important in geomorphic (exhumation and sediment flux), structural (subsidence/flexure) and natural resource (predictive modelling) studies. Determining sedimentation rates for ancient sedimentary sequences is often hampered by poor understanding of stratigraphic architecture, long‐term variability in large‐scale sediment dispersal patterns and inconsistent availability of absolute age data. Uranium–Lead (U‐Pb) detrital zircon (DZ) geochronology is not only a popular method to determine the provenance of siliciclastic sedimentary rocks but also helps delimit the age of sedimentary sequences, especially in basins associated with protracted volcanism. This study assesses the reliability of U‐Pb DZ ages as proxies for depositional ages of Upper Cretaceous strata in the Magallanes‐Austral retroarc foreland basin of Patagonia. Progressive younging of maximum depositional ages (MDAs) calculated from young zircon populations in the Upper Cretaceous Dorotea Formation suggests that the MDAs are potential proxies for absolute age, and constrain the age of the Dorotea Formation to be ca. 82–69 Ma. Even if the MDAs do not truly represent ages of contemporaneous volcanic eruptions in the arc, they may still indicate progressive‐but‐lagged delivery of increasingly younger volcanogenic zircon to the basin. In this case, MDAs may still be a means to determine long‐term (≥1–2 Myr) average sedimentation rates. Burial history models built using the MDAs reveal high aggradation rates during an initial, deep‐marine phase of the basin. As the basin shoaled to shelfal depths, aggradation rates decreased significantly and were outpaced by progradation of the deposystem. This transition is likely linked to eastward propagation of the Magallanes fold‐thrust belt during Campanian‐Maastrichtian time, and demonstrates the influence of predecessor basin history on foreland basin dynamics.

Recent advances in our understanding of palaeovalleys are largely guided by examples from passive margins, in which accommodation increases down depositional dip. This study tests these models against a dataset from the Pennsylvanian Breathitt Group of the central Appalachian foreland basin, USA. This fluvio‐deltaic succession contains extensive erosionally based fluvio‐estuarine sand bodies that can be tracked over 80 km down depositional dip from a proximal zone of high accommodation close to the orogenic margin to a distal, lower accommodation zone close to the cratonic margin of the basin. The sand bodies are up to 25 m thick, multi‐storey and characterized in their lower parts by strongly amalgamated storeys containing sandy fluvial to estuarine bar accretion elements, and in their middle to upper parts by more fully preserved storeys up to 10 m thick and laterally extensive over 100s of metres. The upper storeys include abandonment channel‐fills of heterolithic marine or marginal marine deposits or muddy to sandy point‐bar elements. Three major regional‐scale architectures include: (i) Tabular sand bodies that everywhere incise open marine prodelta and mouth bar facies and are interpreted as palaeovalleys formed during falling stage and lowstand systems tracts, when eustatic sea‐level fall outpaced tectonic subsidence across the entire study area. (ii) Sand bodies that incise genetically related floodplain lake and/or bay‐fill minor mouth bar deposits up depositional dip and open marine prodelta and mouth bar facies down dip. These stacked distributary channel deposits map down dip into palaeovalleys and formed when up dip subsidence rate resulted in positive, but reduced rate of accommodation creation, while lower tectonic subsidence rate down‐dip resulted in incision. (iii) Sand bodies that incise genetically related floodplain, lake and/or bay‐fill minor mouth bars up dip and pass down‐dip into genetically related unconfined floodplain, prodelta and mouth bar deposits. These sand bodies represent stacked distributary channel fills and channel amalgamation was the product of high rates of lateral migration, typical of the behaviour of channels above their backwater reach. Case (2) sand bodies demonstrate that in rapidly subsiding foreland basins, cross‐shelf palaeovalleys may form down depositional dip from aggradational, distributive fluival strata. Additionally, the genetic relationship between stacked distributary channels and palaeovalleys supports recent models for palaeovalley formation that emphasize diachronous, cut‐and‐fill during falling stage and lowstands of relative sea level.

We present an interpolation model that describes Holocene groundwater level rise and the creation of accommodation space in 3D in the Rhine‐Meuse delta – the Netherlands. The model area (ca. 12 400 km²) covers two palaeovalleys of Late Pleistocene age (each 30 km wide) and the overlying Holocene deposits of the Rhine‐Meuse delta, the Holland coastal plain, and the Zuiderzee former lagoon. Water table rise is modelled from 10 800 to 1000 cal. BP, making use of age‐depth relations based on 384 basal peat index points, and producing output in the form of stacked palaeo groundwater surfaces, groundwater age‐depth curves, and voxel sets. These products allow to resolve (i) regional change and variations of inland water table slopes, (ii) spatial differences in the timing and pacing of transgression, and (iii) analysis of interplay of coastal, fluvial and subsidence controls on the provision of accommodation space. The interpolation model is a multi‐parameter trend function, to which a 3D‐kriging procedure of the residuals is added. This split design deploys a generic approach for modelling provision of accommodation space in deltas and coastal lowlands, aiming to work both in areas of intermediate data availability and in the most data‐rich environments. Major provision of accommodation space occurred from 8500 cal BP onwards, but a different evolution occurred in each of the two palaeovalleys. In the northern valley, creation of accommodation space began to stall at 7500 cal BP, while in the southern valley provision of new accommodation space in considerable quantities continued longer. The latter is due to the floodplain gradient that was maintained by the Rhine, which distinguishes the fluvial deltaic environment from the rest of the back‐barrier coastal plain. The interpolation results allow advanced mapping and investigation of apparent spatial differences in Holocene aggradation in larger coastal sedimentary systems. Furthermore, they provide a means to generate first‐order age information with centennial precision for 3D geological subsurface models of Holocene deltas and valley fills. As such, the interpolation is of use in studies into past and present land subsidence and into low land sedimentation.

Subduction zones provide direct insight into plate boundary deformation and by studying these areas we better understand tectonic processes and variability over time. We studied the structure of the offshore subduction zone system of the Pampean flat‐slab segment (ca. 29–33°S) of the Chilean margin using seismic and bathymetric constraints. Here, we related and analysed the structural styles of the offshore and onshore western fore‐arc. Overlying the acoustic top of the continental basement, two syn‐extensional seismic sequences were recognised and correlated with onshore geological units and the Valparaíso Forearc Basin seismic sequences: (SII) Pliocene‐Pleistocene and (SI) Miocene‐Pliocene (Late Cretaceous (?) to Miocene‐Pliocene) syn‐extensional sequences. These sequences are separated by an unconformity (i.e. Valparaíso Unconformity). Seismic reflection data reveal that the eastward dipping extensional system (EI) recognised at the upper slope can be extended to the middle slope and controlled the accumulation of the older seismic package (SI). The westward dipping extensional system (EII) is essentially restricted to the middle slope. Here, EII cuts the eastward dipping extensional system (EI), preferentially parallel to the inclination of the older sequences (SI), and controlled a series of middle slope basins which are filled by the Pliocene‐Pleistocene seismic sequence (SII). At the upper slope and in the western Coastal Cordillera, the SII sequence is controlled by eastward dipping faults (EII) which are the local reactivation of older extensional faults (EI). The tectonic boundary between the middle (eastern outermost forearc block) and upper continental slope (western coastal block) is a prominent system of trenchward dipping normal fault scarps (ca. 1 km offset) that resemble a major trenchward dipping extensional fault system. This prominent structural feature can be readily detected along the Chilean erosive margin as well as the two extensional sets (EI and EII). Evidence of slumping, thrusting, reactivated faults and mass transport deposits, were recognised in the slope domain and locally restricted to some eastern dipping faults. These features could be related to gravitational effects or slope deformation due to coseismic deformation. The regional inclination of the pre‐Pliocene sequences favoured the gravitational collapse of the outermost forearc block. We propose that the structural configuration of the study area is dominantly controlled by tectonic erosion as well as the uplift of the Coastal Cordillera, which is partially controlled by pre‐Pliocene architecture.

Injectites sourced from base‐of‐slope and basin‐floor parent sandbodies are rarely reported in comparison to submarine slope channel systems. This study utilizes the well‐constrained palaeogeographic and stratigraphic context of three outcrop examples exposed in the Karoo Basin, South Africa, to examine the relationship between abrupt stratigraphic pinchouts in basin‐floor lobe complexes, and the presence, controls, and character of injectite architecture. Injectites in this palaeogeographic setting occur where there is: (i) sealing mudstone both above and below the parent sand to create initial overpressure; (ii) an abrupt pinchout of a basin‐floor lobe complex through steep confinement to promote compaction drive; (iii) clean, proximal sand beds aiding fluidization; and (iv) a sharp contact between parent sand and host lithology generating a source point for hydraulic fracture and resultant injection of sand. In all outcrop cases, dykes are orientated perpendicular to palaeoslope, and the injected sand propagated laterally beneath the parent sand, paralleling the base to extend beyond its pinchout. Understanding the mechanisms that determine and drive injection is important in improving the prediction of the location and character of clastic injectites in the subsurface. Here, we highlight the close association of basin‐floor stratigraphic traps and sub‐seismic clastic injectites, and present a model to explain the presence and morphology of injectites in these locations.

Pervasive fracture networks are common in many reservoir‐scale carbonate bodies even in the absence of large deformation and exert a major impact on their mechanical and flow behaviour. The Upper Cretaceous Jandaíra Formation is a few hundred meters thick succession of shallow water carbonates deposited during the early post‐rift stage of the Potiguar rift (NE Brazil). The Jandaíra Formation in the present onshore domain experienced <1.5 km thermal subsidence and, following Tertiary exhumation, forms outcrops over an area of >1000 km². The carbonates have a gentle, <5⁰, dip to the NE and are affected by few regional, low displacement faults or folds. Despite their simple tectonic history, carbonates display ubiquitous open fractures, sub‐vertical veins, and sub‐vertical as well as sub‐horizontal stylolites. Combining structural analysis, drone imaging, isotope studies and mathematical modelling, we reconstruct the fracturing history of the Jandaíra Formation during and following subsidence and analyse the impact fractures had on coeval fluid flow. We find that Jandaíra carbonates, fully cemented after early diagenesis, experienced negligible deformation during the first few hundreds of meters of subsidence but were pervasively fractured when they reached depths >400–500 m. Deformation was accommodated by a dense network of sub‐vertical mode I and hybrid fractures associated with sub‐vertical stylolites developed in a stress field characterised by a sub‐horizontal σ1 and sub‐vertical σ2. The development of a network of hybrid fractures, rarely reported in the literature, activated the circulation of waters charged in the mountainous region, flowing along the porous Açu sandstone underlying the Jandaíra carbonates and rising to the surface through the fractured carbonates. With persisting subsidence, carbonates reached depths of 800–900 m entering a depth interval characterised by a sub‐vertical σ1. At this stage, sub‐horizontal stylolites developed liberating calcite which sealed the sub‐vertical open fractures transforming them in veins and preventing further flow. During Tertiary exhumation, several of the pre‐existing veins and stylolites opened and became longer, and new fractures were created typically with the same directions of the older features. The simplicity of our model suggests that most rocks in passive margin settings might have followed a similar evolution and thus display similar structures.