Crust‐scale 3D model of the Western Bredasdorp Basin (Southern South Africa): data‐based insights from combined isostatic and 3D gravity modelling

The southern South African continental margin documents a complex margin system that has undergone both continental rifting and transform processes in a manner that its present‐day architecture and geodynamic evolution can only be better understood through the application of a multidisciplinary and multi‐scale geo‐modelling procedure. In this study, we focus on the proximal section of the larger Bredasdorp sub‐basin (the westernmost of the five southern South African offshore Mesozoic sub‐basins), which is hereto referred as the Western Bredasdorp Basin. Integration of 1200 km of 2D seismic‐reflection profiles, well‐logs and cores yields a consistent 3D structural model of the Upper Jurassic‐Cenozoic sedimentary megasequence comprising six stratigraphic layers that represent the syn‐rift to post‐rift successions with geometric information and lithology‐depth‐dependent properties (porosities and densities). We subsequently applied a combined approach based on Airy's isostatic concept and 3D gravity modelling to predict the depth to the crust‐mantle boundary (Moho) as well as the density structure of the deep crust. The best‐fit 3D model with the measured gravity field is only achievable by considering a heterogeneous deep crustal domain, consisting of an uppermost less dense prerift meta‐sedimentary layer [ρ = 2600 kg m^?3] with a series of structural domains. To reproduce the observed density variations for the Upper Cenomanian–Cenozoic sequence, our model predicts a cumulative eroded thickness of ca. 800–1200 m of Tertiary sediments, which may be related to the Late Miocene margin uplift. Analyses of the key features of the first crust‐scale 3D model of the basin, ranging from thickness distribution pattern, Moho shallowing trend, sub‐crustal thinning to shallow and deep crustal extensional regimes, suggest that basin initiation is typical of a mantle involvement deep‐seated pull‐apart setting that is associated with the development of the Agulhas‐Falkland dextral shear zone, and that the system is not in isostatic equilibrium at present day due to a mass excess in the eastern domain of the basin that may be linked to a compensating rise of the asthenospheric mantle during crustal extension. Further corroborating the strike‐slip setting is the variations of sedimentation rates through time. The estimated syn‐rift sedimentation rates are three to four times higher than the post‐rift sedimentation, thereby indicating that a rather fast and short‐lived subsidence during the syn‐rift phase is succeeded by a significantly poor passive margin development in the post‐rift phase. Moreover, the derived lithospheric stretching factors [β = 1.5–1.75] for the main basin axis do not conform to the weak post‐rift subsidence. This therefore suggests that a differential thinning of the crust and the mantle‐lithosphere typical for strike‐slip basins, rather than the classical uniform stretching model, may be applicable to the Western Bredasdorp Basin.     

The architecture of submarine monogenetic volcanoes – insights from 3D seismic data

Many prospective sedimentary basins contain a variety of extrusive volcanic products that are ultimately sourced from volcanoes. However, seismic reflection‐based studies of magmatic rift basins have tended to focus on the underlying magma plumbing system, meaning that the seismic characteristics of volcanoes are not well understood. Additionally, volcanoes have similar morphologies to hydrothermal vents, which are also linked to underlying magmatic intrusions. In this study, we use high resolution 3D seismic and well data from the Bass Basin, offshore southern Australia, to document 34 cone‐ and crater‐type vents of Miocene age. The vents overlie magmatic intrusions and have seismic properties indicative of a volcanic origin: their moderate–high amplitude upper reflections and zones of “wash‐out” and velocity pull‐up beneath. The internal reflections of the vents are similar to those found in lava deltas, suggesting they are composed of volcaniclastic material. This interpretation is corroborated by data from exploration wells which penetrated the flanks of several vents. We infer that the vents we describe are composed of hyaloclastite and pyroclasts produced during submarine volcanic eruptions. The morphology of the vents is typical of monogenetic volcanoes, consistent with the onshore record of volcanism on the southern Australian margin. Based on temporal, spatial and volumetric relationships, we propose that submarine volcanoes can evolve from maars to tuff cones as a result of varying magma‐water interaction efficiency. The morphologies of the volcanoes and their links to the underlying feeder systems are superficially similar to hydrothermal vents. This highlights the need for careful seismic interpretation and characterization of vent structures linked to magmatic intrusions within sedimentary basins.     

Investigating the response of Cladocera to a major saltwater intrusion event in an Arctic lake from the outer Mackenzie Delta (NT, Canada)

An increase in the frequency and intensity of marine storm surges is a predicted consequence of climate warming, and therefore it is important to better understand the biological responses to such events in coastal regions. In late September 1999, a major storm surge resulted in a saltwater intrusion event over a large area of the Mackenzie Delta (NT, Canada) front, causing rapid salinization of lakes on the alluvial plain. Due to a lack of long-term ecological monitoring data in the region, the impacts that the saltwater intrusion event had on the biota of affected lakes were unknown. We used high-resolution paleolimnological approaches to reconstruct past assemblage changes in Cladocera from impacted Lake DZO-29 (unofficial name) in order to determine how different cladoceran species responded to a major increase in lake salinity following the 1999 storm surge. Camptocercus were extirpated from the lake following the saltwater intrusion and have not recovered. We also observed an initial decrease in Alona relative abundance following the marine flooding, likely reflecting a loss of A. quadrangularis, A. barbulata, and A. costata from the lake. A. circumfimbriata, Chydorus biovatus, C. brevilabris, and Bosmina spp. were abundant both before and after the saltwater intrusion, and Paralona pigra was present following the storm surge, but not prior to it. The most notable shift in Cladocera in the recent sedimentary record, however, occurred much earlier, with an increase in pelagic Bosmina taxa and a subsequent decrease in the benthic/littoral taxa Chydorus and Camptocercus, an assemblage shift that is consistent with a response to climate warming in this region, and strongly correlated to other changes in the lake inferred to be as a result of regional warming.     

Application of multi‐kinetic apatite fission track and (U‐Th)/He thermochronology to source rock thermal history: a case study from the Mackenzie Plain,NWT, Canada

Shale of the Upper Cretaceous Slater River Formation extends across the Mackenzie Plain of the Canadian Northwest Territories and has potential as a regional source rock because of the high organic content and presence of both oil‐ and gas‐prone kerogen. An understanding of the thermal history experienced by the shale is required to predict any potential petroleum systems. Our study integrates multi‐kinetic apatite fission track (AFT) and apatite (U‐Th)/He (AHe) thermochronometers from a basal bentonite unit to understand the timing and magnitude of Late Cretaceous burial experienced by the Slater River Formation along the Imperial River. We use LA‐ICP‐MS and EPMA methods to assess the chemistry of apatite, and use these values to derive the AFT kinetic parameter r_mr0. Our AFT dates and track lengths, respectively, range from 201.5 ± 36.9 Ma to 47.1 ± 12.3 Ma, and 16.8 to 10.2 μm, and single crystal AHe dates are between 57.9 ± 3.5 and 42.0 ± 2.5 Ma with effective uranium concentrations from 17 ppm to 36 ppm. The fission track data show no relationship with the kinetic parameter D_par and fail the χ²‐test indicating that the data do not comprise a single statistically significant population. However, when plotted against their r_mr0 value, the data are separated into two statistically significant kinetic populations with distinct track length distributions. Inverse thermal history modelling of both the multi‐kinetic AFT and AHe datasets, reveal that the Slater River Formation reached maximum burial temperatures of ~65–90 °C between the Turonian and Paleocene, indicating that the source rock matured to the early stages of hydrocarbon generation, at best. Ultimately, our data highlight the importance of kinetic parameter choice for AFT and AHe thermochronology, as slight variations in apatite chemistry may have significant implications on fission track and radiation damage annealing in apatite with protracted thermal histories through the uppermost crust.     

Sedimentary and thermal evolution of the Eocene‐Oligocene mudrocks from the southwestern Thrace Basin (NE Greece)

Paleothermal indicators based on clay mineral and organic matter analyses, were integrated with mudrock geochemistry and stratigraphic data to define the sedimentary evolution of the southwestern Thrace Basin during the Eocene to Oligocene. This multi‐method approach allowed us to reconstruct the burial evolution of the basin in Eocene and Oligocene times and to study the mudrock composition and relate this to their provenance and source area weathering. The studied mudrocks show similar chemical variations. The distribution of some major and trace elements for the studied samples reflect heterogeneous source areas containing both felsic to mafic rocks. In particular, the Light Rare Earth Elements/Transition elements (LREEs/TEs) ratios are very high for the Avdira and Organi samples (on the average between 1.5 and 2.2 for (La + Ce)/Cr and 3.5–8 for (La + Ce)/Ni), suggesting a felsic source(s), and very low for the Samothraki, Limnos, Paterma and Iasmos samples (on the average between 0.4 and 0.6 for (La + Ce)/Cr and 0.6–1 for (La + Ce)/Ni), suggesting a mainly basic source(s). The mineralogical composition coupled with the A‐CN‐K and A‐N‐K plots suggest a complex evolution. The clay mineral data (illite percentage in I/S and the stacking order R and the Kübler Index) coupled to vitrinite reflectance analysis indicate a high to intermediate diagenetic grade for the Middle to Upper Eocene samples (from Iasmos, Gratini, Organi, Paterma, Esimi and Samotraki sections) and a low diagenetic grade for the Upper Eocene to Oligocene samples (from Limnos and Avdira sections). These data helped in interpreting the geodynamic evolution of the studied basins where the magmatic activity plays an important role. In particular, Middle to Upper Eocene sediments show high to intermediate diagenetic grade since they are located in a portion of the basin dominated by Eocene to Oligocene magmatic activity and intrusion of granitoids, whereas, the Upper Eocene to Oligocene sediments are not involved in important magmatic activity and intrusion of granitoids and, thus, show low diagenetic grade. Furthermore, Middle to Upper Eocene sediments experienced deeper burial processes caused by lithostatic load, rather than the uppermost Eocene and Oligocene sediments, in relation of their position along the stratigraphic succession. These data suggest a burial depth of at least 3–4 km with a tectonic exhumation mainly related to the extensional phases of the Miocene age.     

3D numerical modelling of marine organic matter distribution: example of the early Jurassic sequences of the Lusitanian Basin (Portugal)

Due to the multiple controlling factors involved, it is a challenging task to identify and quantify the processes influencing the distribution and heterogeneity of marine organic‐rich rocks. To improve our understanding of these deposits, we model their burial history and stratigraphic evolution as well as processes linked to marine organic matter history throughout the Lower Jurassic in of the Northern Lusitanian Basin (Western Iberian Margin). This 15‐Ma‐long interval is modelled using 100‐kyr time steps to simulate lithologies and organic matter heterogeneity as layers with a thickness of 2–5 m, depending on the sedimentation rate in the basin. The model is calibrated by well and outcrop data which provide structural and biostratigraphic constraints, as well as information on the depositional facies and geochemistry of the sediments. The results show that the presence of organic‐rich intervals is linked to first‐order variations in the basin geometry and sedimentation rates. Without considering any variation of primary productivity or oxygen content in surface sea waters, the parameters of basin geometry and sedimentation rate are sufficient to predict the main characteristics of source rocks, i.e. their occurrence, thickness and mineralogy at the basin scale. However, to fit the measured organic carbon contents, we need to take account of other parameters such as variations of primary productivity or changes in dissolved oxygen concentration.     

Tectonic heritage in drainage pattern and dynamics: the case of the French South Alpine Foreland Basin (ca. 45–20 Ma)

The spatial and temporal organization of depositional environments in drainage networks of foreland basins reflect the tectonic and erosional dynamics associated with the development of mountain belts. We provide field evidences for the initiation and evolution of a complex drainage system in the French South Alpine Foreland Basin related to Western Alps exhumation. Sedimentological and structural analyses of the Eocene–Early Miocene succession were investigated in the (1) Argens/Peyresq, (2) Barrême/Blieux/Taulanne and (3) Montmaur/St‐Disdier sectors. Combined with the existing structural data set, we propose a new model that integrates the regional tectonic activity, the palaeovalley orientation and their dynamics through time. The Eocene–Miocene deposits clearly show the existence of N–S‐oriented palaeovalleys. The systematic presence of early NE–SW‐ to N–S‐oriented strike‐slip and extensional faults in the palaeovalleys suggests that these tectonic structures were responsible for the formation of the initial N–S‐oriented basin‐floor topographies. The vertical offset of the strike‐slip faults induced sufficient accommodation space for the Cenozoic sedimentation since the Middle Eocene. It implies the creation of N–S‐oriented palaeovalleys during the northward Pyrenean‐Provençal phase, pre‐dating westward Alpine compression. Later, the Oligocene Alpine tectonic phase induced drainage expansion toward the orogenic wedge and the erosion of the exhumed internal massifs by transverse streams. The establishment of new connections between the old topographic lows formed a longitudinal drainage pattern that remains the locus of deposition in a regional sedimentary routing system. In this model, former strike‐slip faults correspond to weakness zones overprinted by the westward Alpine shortening that allowed the formation of the modern piggyback basin structure of the foreland and the long‐time preservation of the palaeovalley geometry.     

Quantifying the influence of sill intrusion on the thermal evolution of organic‐rich sedimentary rocks in nonvolcanic passive margins: an example from ODP 210‐1276, offshore Newfoundland,Canada

Intrusive magmatism is an integral and understudied component in both volcanic and nonvolcanic passive margins. Here, we investigate the thermal effects of widespread (ca. 20 000 km²) intrusive magmatism on the thermal evolution of organic‐rich sedimentary rocks on the nonvolcanic Newfoundland passive margin. ODP 210‐1276 (45.41°N, 44.79°W) intersects two sills: an older, upper sill and a younger, lower sill that are believed to correspond to the high amplitude ‘U‐reflector’ observed across the Newfoundland Basin. A compilation of previous work collectively provides; (1) emplacement depth constraints, (2) vitrinite reflectance data and (3) ⁴⁰Ar/³⁹Ar dates. Collectively, these data sets provide a unique opportunity to model the conductive cooling of the sills and how they affect thermal maturity of the sedimentary sequence. A finite differences method was used to model the cooling of the sills, with the model outputs then being entered into the EASY%R_o vitrinite reflectance model. The modelled maturation profile for ODP 210‐1276 shows a significant but localized effect on sediment maturity as a result of the intrusions. Our results suggest that even on nonvolcanic margins, intrusive magmatism can significantly influence the thermal evolution in the vicinity of igneous intrusions. In addition, the presence of widespread sills on nonvolcanic passive margins such as offshore Newfoundland may be indicative of regional‐scale thermal perturbations that should be considered in source rock maturation studies.     

Comment on Non‐unique stratal geometries: implications for sequence stratigraphic interpretations,by: P.M. Burgess and G.D. Prince,Basin Research (2015) 27, 351–365

The non‐unique variability highlighted by Burgess & Prince (Basin Res. 2015, 27 , 351) (i.e. the origin and timing of maximum flooding surfaces, maximum regressive surfaces and subaerial unconformities; the process of topset aggradation in relation with the various types of shoreline trajectory; and the multiple controls that may affect the progradation and retrogradation of a shoreline) is irrelevant to the workflow of sequence stratigraphy. What is relevant is the observation of the unique stratal geometries that are diagnostic to the definition of all units and surfaces of sequence stratigraphy. In downstream‐controlled settings, these unique stratal stacking patterns relate to the forced regressive, normal regressive and transgressive shoreline trajectories. Multiple controls interplay during the formation of each type of stacking pattern, including accommodation, sediment supply and the energy of the sediment‐transport agents. This interplay explains the non‐unique variability, but does not change the unique criteria that afford a consistent application of sequence stratigraphy. Failure to rationalize the non‐unique variability within the context of unique stratal geometries is counterproductive, and obscures the simple workflow of sequence stratigraphy.     

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.     

Using basin thermal history to evaluate the role of Miocene–Pliocene flat‐slab subduction in the southern Central Andes (27° S–30° S)

Studies in both modern and ancient Cordilleran‐type orogenic systems suggest that processes associated with flat‐slab subduction control the geological and thermal history of the upper plate; however, these effects prove difficult to deconvolve from processes associated with normal subduction in an active orogenic system. We present new geochronological and thermochronological data from four depositional areas in the western Sierras Pampeanas above the Central Andean flat‐slab subduction zone between 27° S and 30° S evaluating the spatial and temporal thermal conditions of the Miocene–Pliocene foreland basin. Our results show that a relatively high late Miocene–early Pliocene geothermal gradient of 25–35 °C km⁻¹ was typical of this region. The absence of along‐strike geothermal heterogeneities, as would be expected in the case of migrating flat‐slab subduction, suggests that either the response of the upper plate to refrigeration may be delayed by several millions of years or that subduction occurred normally throughout this region through the late Miocene. Exhumation of the foreland basin occurred nearly synchronously along strike from 27 to 30° S between ca. 7 Ma and 4 Ma. We propose that coincident flat‐slab subduction facilitated this wide‐spread exhumation event. Flexural modelling coupled with geohistory analysis show that dynamic subsidence and/or uplift associated with flat‐slab subduction is not required to explain the unique deep and narrow geometry of the foreland basin in the region implying that dynamic processes were a minor component in the creation of accommodation space during Miocene–Pliocene deposition.