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.     

Late Cretaceous‐Palaeogene stratigraphic and basin evolution in the Zhepure Mountain of southern Tibet: implications for the timing of India‐Asia initial collision

This article presents combined stratigraphic, sedimentological, subsidence and provenance data for the Cretaceous–Palaeogene succession from the Zhepure Mountain of southern Tibet. This region records the northernmost sedimentation of the Tethyan passive margin of India, and this time interval represents the transition into continental collision with Asia. The uppermost Cretaceous Zhepure Shanpo and Jidula formations record the transition from pelagic into upper slope to delta‐plain environments. The Palaeocene–lower Eocene Zongpu Formation records a carbonate ramp that is overlain by the deep‐water Enba Formation (lower Eocene). The upper part of the Enba Formation records shallowing into a storm‐influenced, outer shelf environment. Detrital zircon U–Pb and Hf isotopic data indicate that the terrigenous strata of the Enba Formation were sourced from the Lhasa terrane. Unconformably overlying the Enba Formation is the Zhaguo Formation comprising fluvial deposits with evidence of recycling from the underlying successions. Backstripped subsidence analysis indicates shallowing during latest Cretaceous‐earliest Palaeocene time (Zhepure Shanpo and Jidula formations) driven by basement uplift, followed by stability (Zongpu Formation) until early Eocene time (Enba Formation) when accelerated subsidence occurred. The provenance, subsidence and stratigraphy suggest that the Enba and Zhaguo formations record foredeep and wedge‐top sedimentation respectively within the early Himalayan foreland basin. The underlying Zongpu Formation is interpreted to record the accumulation of a carbonate ramp at the margin of a submarine forebulge. The precursor tectonic uplift during latest Cretaceous time could either record surface uplift over a mantle plume related to the Réunion hotspot, or an early signal of lithospheric flexure related to oceanic subduction, continental collision or ophiolite obduction. The results indicate that the collision of India with Asia occurred before late Danian (ca. 62 Ma) time.     

Application of seismic and sequence stratigraphic concepts to a lava‐fed delta system in the Faroe‐Shetland Basin,UK and Faroes

Detailed seismic stratigraphic analysis of 2D seismic data over the Faroe‐Shetland Escarpment has identified 13 seismic reflection units that record lava‐fed delta deposition during discrete periods of volcanism. Deposition was dominated by progradation, during which the time shoreline migrated a maximum distance of ～44 km in an ESE direction. Localised collapse of the delta front followed the end of progradation, as a decrease in volcanic activity left the delta unstable. Comparison with modern lava‐fed delta systems on Hawaii suggests that syn‐volcanic subsidence is a potential mechanism for apparent relative sea level rise and creation of new accommodation space during lava‐fed delta deposition. After the main phase of progradation, retrogradation of the delta occurred during a basinwide syn‐volcanic relative sea level rise where the shoreline migrated a maximum distance of ～75 km in a NNW direction. This rise in relative sea level was of the order of 175–200 m, and was followed by the progradation of smaller, perched lava‐fed deltas into the newly created accommodation space. Active delta deposition and the emplacement of lava flows feeding the delta front lasted ～2600 years, although the total duration of the lava‐fed delta system, including pauses between eruptions, may have been much longer.     

Propagation history and passive rotation of mesoscale normal faults: implications for synrift stratigraphic development

Field data from onshore exposures of the Oligo-Miocene Gulf of Suez Rift in the Sinai document the passive rotation of early formed mesoscale synthetic and antithetic faults and associated half-graben due to long-lived activity on large displacement (2–5 km) block-bounding faults. Early formed small-displacement (<350 m) mesoscale antithetic faults and half-graben within regional-scale fault blocks underwent progressive steepening due to footwall uplift, rotational faulting and footwall flexing on large-displacement, block-bounding faults. In contrast, mesoscale synthetic faults were progressively rotated to shallower angles. Analysis of palaeohorizontal surfaces within synrift sediments deposited in half-graben adjacent to the mesoscale faults indicate passive rotations of up to 25° about horizontal axes since deposition. Passive burial and in-filling of early formed mesoscale faults and half-graben by synrift sediments is consistent with extension being transferred from numerous mesoscale faults to few block-bounding macroscale faults as extension preceded. Furthermore, this transfer of extension appears to be associated with a marked change in basin configuration, synrift sediment dispersal patterns and facies development. Identification of early formed, passively rotated normal faults and half-graben is important for correctly reconstructing the early stages of basin palaeogeography and sediment dispersal, and for addressing models of rift basin evolution.     

An integrated model of clastic injectites and basin floor lobe complexes: implications for stratigraphic trap plays

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.     

Factors controlling stratal pattern and facies distribution of fluvio‐lacustrine sedimentation in the Sivas mini‐basins,Oligocene (Turkey)

The Sivas Basin, located in the Central Anatolian Plateau of Turkey, is a foreland basin that records a complex interaction between sedimentation, salt tectonics and regional shortening during the Oligo‐Miocene leading to the formation of numerous mini‐basins. The Oligocene sedimentary infill of the mini‐basins consists of a thick continental succession, the Karayün Formation, comprising a vertical succession of three main sub‐environments: (i) playa‐lake, (ii) fluvial braided, and (iii) saline lacustrine. These sub‐environments are seen as forming a large Distributive Fluvial System (DFS) modified through time as a function of sediment supply and accommodation related to regional changes in climate and tectonic regime. Within neighbouring mini‐basins and despite a similar vertical stratigraphic succession, subtle variations in facies assemblages and thickness are observed in stratigraphic units of equivalent age, thus demonstrating the local control exerted by halokinesis. Stratigraphic and stratal patterns reveal in great detail the complex interaction between salt tectonics and sedimentation including different types of halokinetic structures such as hooks, wedges and halokinetic folds. The regional variations of accommodation/sediment supply led to coeval changes in the architectural patterns recorded in the mini‐basins. The type of accommodation regime produces several changes in the sedimentary record: (i) a regime dominated by regional accommodation limits the impact of halokinesis, which is recorded as very small variations in stratigraphic thickness and facies distribution within and between mini‐basins; (ii) a regime dominated by localized salt‐induced accommodation linked to the subsidence of each individual mini‐basin enhances the facies heterogeneity within the DFS, causing sharp changes in stratigraphic thickness and facies assemblages within and between mini‐basins.     

Low‐accommodation and backwater effects on sequence stratigraphic surfaces and depositional architecture of fluvio‐deltaic settings (Cretaceous Mesa Rica Sandstone,Dakota Group,USA)

The adequate documentation and interpretation of regional‐scale stratigraphic surfaces is paramount to establish correlations between continental and shallow marine strata. However, this is often challenged by the amalgamated nature of low‐accommodation settings and control of backwater hydraulics on fluvio‐deltaic stratigraphy. Exhumed examples of full‐transect depositional profiles across river‐to‐delta systems are key to improve our understanding about interacting controlling factors and resultant stratigraphy. This study utilizes the ~400 km transect of the Cenomanian Mesa Rica Sandstone (Dakota Group, USA), which allows mapping of down‐dip changes in facies, thickness distribution, fluvial architecture and spatial extent of stratigraphic surfaces. The two sandstone units of the Mesa Rica Sandstone represent contemporaneous fluvio‐deltaic deposition in the Tucumcari sub‐basin (Western Interior Basin) during two regressive phases. Multivalley deposits pass down‐dip into single‐story channel sandstones and eventually into contemporaneous distributary channels and delta‐front strata. Down‐dip changes reflect accommodation decrease towards the paleoshoreline at the Tucumcari basin rim, and subsequent expansion into the basin. Additionally, multi‐storey channel deposits bound by erosional composite scours incise into underlying deltaic deposits. These represent incised‐valley fill deposits, based on their regional occurrence, estimated channel tops below the surrounding topographic surface and coeval downstepping delta‐front geometries. This opposes criteria offered to differentiate incised valleys from flood‐induced backwater scours. As the incised valleys evidence relative sea‐level fall and flood‐induced backwater scours do not, the interpretation of incised valleys impacts sequence stratigraphic interpretations. The erosional composite surface below fluvial strata in the continental realm represents a sequence boundary/regional composite scour (RCS). The RCS’ diachronous nature demonstrates that its down‐dip equivalent disperses into several surfaces in the marine part of the depositional system, which challenges the idea of a single, correlatable surface. Formation of a regional composite scour in the fluvial realm throughout a relative sea‐level cycle highlights that erosion and deposition occur virtually contemporaneously at any point along the depositional profile. This contradicts stratigraphic models that interpret low‐accommodation settings to dominantly promote bypass, especially during forced regressions. Source‐to‐sink analyses should account for this in order to adequately resolve timing and volume of sediment storage in the system throughout a complete relative sea‐level cycle.     

Preserving coastal natural character: Court interpretations of a long‐standing New Zealand policy goal

New Zealand's policy goal to preserve coastal natural character was first incorporated into planning legislation in 1973 and protected area legislation in 1977. An evaluation of 100 relevant Resource Management Act Court decisions determined the Courts' interpretations of this policy goal. The Courts have repeatedly found that natural character is of nature (not culture), including natural elements, patterns, and processes. There is a continuum of naturalness from pristine indigenous ecosystems to the urban environment. A 2014 Supreme Court decision has upheld the use of environmental bottom lines for natural character and the importance of protection as part of sustainable management.     

Impact of sediment transport efficiency on large-scale sequence architecture: results from stratigraphic computer simulation

Our understanding of the formation of the wide range of sequence architectures we observe in the rock record is still somewhat limited. The sedimentary response to the complex interaction of various time-variable basin-scale processes such as subsidence, eustasy and sediment supply is difficult to understand without numerical models. The computer simulation model presented here, DEMOSTRAT, is a powerful tool to investigate the sequence development scenarios in a 2-D dip section. The model includes tectonic subsidence, eustasy, two-component (sand and mud) nonlinear diffusional sediment transport, compaction and isostasy. The transport coefficients in the diffusion equations express the system's ability to transport sand and mud, and are mainly dependent on climate and subaquatic processes. Keeping other model input parameters constant, the magnitude of transport coefficients seems to have an important impact on sequence development. With high transport coefficients, extensive erosion during sea-level fall and lack of sediment buildup above sea-level during rise may reduce the preservation potential for nonmarine sediments. In addition, the former slope break will be eroded during transgression, forming sand-rich slope or basin floor sediments that may be misinterpreted as lowstand fans. Moreover, the magnitude of transport coefficients has an impact on unconformity timing and development (shown in Wheeler plots).
     

Structural controls on the stratigraphic architecture of net‐transgressive shallow‐marine strata in a salt‐influenced rift basin: Middle‐to‐Upper Jurassic Egersund Basin,Norwegian North Sea

In this study, we integrate 3D seismic reflection, wireline log, biostratigraphic and core data from the Egersund Basin, Norwegian North Sea to determine the impact of syn‐depositional salt movement and associated growth faulting on the sedimentology and stratigraphic architecture of the Middle‐to‐Upper Jurassic, net‐transgressive, syn‐rift succession. Borehole data indicate that Middle‐to‐Upper Jurassic strata consist of low‐energy, wave‐dominated offshore and shoreface deposits and coal‐bearing coastal‐plain deposits. These deposits are arranged in four parasequences that are aggradationally to retrogradationally stacked to form a net‐transgressive succession that is up to 150‐m thick, at least 20 km in depositional strike (SW‐NE) extent, and >70 km in depositional dip (NW‐SE) extent. In this rift‐margin location, changes in thickness but not facies are noted across active salt structures. Abrupt facies changes, from shoreface sandstones to offshore mudstones, only occur across large displacement, basement‐involved normal faults. Comparisons to other tectonically active salt‐influenced basins suggest that facies changes across syn‐depositional salt structures are observed only where expansion indices are >2. Subsidence between salt walls resulted in local preservation of coastal‐plain deposits that cap shoreface parasequences, which were locally removed by transgressive erosion in adjacent areas of lower subsidence. The depositional dip that characterizes the Egersund Basin is unusual and likely resulted from its marginal location within the evolving North Sea rift and an extra‐basinal sediment supply from the Norwegian mainland.     

A three‐dimensional stratigraphic model for the Messinian crisis in the Pannonian Basin,eastern Hungary

A three‐dimensional quantitative stratigraphic forward model is employed to investigate the controls leading to the Messinian events in the lacustrine Pannonian Basin of Central Paratethys, and the link between the Messinian salinity crisis in the Mediterranean and the late Miocene‐Pliocene stratigraphy of the Pannonian Basin. Subsurface geological data show that a prominent unconformity surface formed during Messinian time in the Pannonian Basin associated with a sudden forced regression, abrupt basinward shift of facies and a subsequent, prolonged lowstand normal regression. The lowstand prograding series filled up the shallow basin fast, while, at the same time, the marginal areas of the basin were subject to tectonic inversion. The Dionisos program used in this research is built on a nonlinear water‐driven sediment diffusion process, and it employs multiple sediment classes, basin flexure and compaction. Four different scenarios were built in the experiments to test possible basin histories with different rates and timing of tectonic inversion. Each scenario was modelled in two versions: including and not including a lake‐level fall in the Messinian. The results confirm that the Pannonian Basin in the study area has undergone a tectonic inversion since the Messinian, although the exact rates of uplift at different locations remain uncertain. The unconformity and the observed stratigraphic architecture and facies pattern could be modelled adequately only in the versions that applied a Messinian lake‐level fall. Our research concludes that the Messinian unconformity in the Pannonian Basin was caused by an absolute lake‐level drop, likely linked to the desiccation of the Mediterranean, followed by subsidence and normal regression in the basin centre and concomitant tectonic inversion and uplift along the basin margins.     

Forward stratigraphic modelling of sediment pathways and depocentres in salt‐influenced passive‐margin basins: Lower Cretaceous,central Scotian Basin

Source‐to‐sink studies and numerical modelling software are increasingly used to better understand sedimentary basins, and to predict sediment distributions. However, predictive modelling remains problematic in basins dominated by salt tectonics. The Lower Cretaceous delta system of the Scotian Basin is well suited for source‐to‐sink studies and provides an opportunity to apply this approach to a region experiencing active salt tectonism. This study uses forward stratigraphic modelling software and statistical analysis software to produce predictive stratigraphic models of the central Scotian Basin, test their sensitivity to different input parameters, assess proposed provenance pathways, and determine the distribution of sand and factors that control sedimentation in the basin. Models have been calibrated against reference wells and seismic surfaces, and implement a multidisciplinary approach to define simulation parameters. Simulation results show that previously proposed provenance pathways for the Early Cretaceous can be used to generate predictive stratigraphic models, which simulate the overall sediment distribution for the central Scotian Basin. Modelling confirms that the shaly nature of the Naskapi Member is the result of tectonic diversion of the Sable and Banquereau rivers and suggests additional episodic diversion during the deposition of the Cree Member. Sand is dominantly trapped on the shelf in all units, with transport into the basin along salt corridors and as a result of turbidity current flows occurring in the Upper Missisauga Formation and Cree Member. This led to sand accumulation in minibasins with a large deposit seawards of the Tantallon M‐41 well. Sand also appears to bypass the basin via salt corridors which lead to the down‐slope edge of the study area. Sensitivity analysis suggests that the grain size of source sediments to the system is the controlling factor of sand distribution. The methodology applied to this basin has applications to other regions complicated by salt tectonics, and where sediment distribution and transport from source‐to‐sink remain unclear.     

Evaluating alluvial stratigraphic response to cyclic and non‐cyclic upstream forcing through process‐based alluvial architecture modelling

Formation of alluvial stratigraphy is controlled by autogenic processes that mix their imprints with allogenic forcing. In some alluvial successions, sedimentary cycles have been linked to astronomically‐driven, cyclic climate changes. However, it remains challenging to define how such cyclic allogenic forcing leads to sedimentary cycles when it continuously occurs in concert with autogenic forcing. Accordingly, we evaluate the impact of cyclic and non‐cyclic upstream forcing on alluvial stratigraphy through a process‐based alluvial architecture model, the Karssenberg and Bridge (2008) model (KB08). The KB08 model depicts diffusion‐based sediment transport, erosion and deposition within a network of channel belts and associated floodplains, with river avulsion dependent on lateral floodplain gradient, flood magnitude and frequency, and stochastic components. We find cyclic alluvial stratigraphic patterns to occur when there is cyclicity in the ratio of sediment supply over water discharge (Q_s/Q_w ratio), in the precondition that the allogenic forcing has sufficiently large amplitudes and long, but not very long, wavelengths, depending on inherent properties of the modelled basin (e.g. basin subsidence, size, and slope). Each alluvial stratigraphic cycle consists of two phases: an aggradation phase characterized by rapid sedimentation due to frequent channel shifting and a non‐deposition phase characterized by channel belt stability and, depending on Q_s/Q_w amplitudes, incision. Larger Q_s/Q_w ratio amplitudes contribute to weaker downstream signal shredding by stochastic components in the model. Floodplain topographic differences are found to be compensated by autogenic dynamics at certain compensational timescales in fully autogenic runs, while the presence of allogenic forcing clearly impacts the compensational stacking patterns.     

Central‐Auckland rainfall, 1853–2020: Sites histories and implications for developing a long‐term rainfall record

Previous work related to central‐Auckland site Albert Park showed that sparse site histories for inner‐city weather stations can be usefully supplemented by mining searchable online historical archives (Papers Past, DigitalNZ), at least as far back as 1909. Assessment of potential is extended here back to 1853, within the context of developing a long‐term rainfall record for Auckland. Although late 19th century sources are thinner and of poorer quality, much useful information could be gleaned about the locations and exposures of early sites. Hypothesised implications for the rainfall record suggest that all early site data need to be adjusted upwards.     

Middle Eocene‐Oligocene broken‐foreland evolution in the Andean Calchaqui Valley,NW Argentina: insights from stratigraphic,structural and provenance studies

Two end‐member models have been proposed for the Paleogene Andean foreland: a simple W‐E migrating foreland model and a broken‐foreland model. We present new stratigraphic, sedimentological and structural data from the Paleogene Quebrada de los Colorados (QLC) Formation, in the Eastern Cordillera, with which to test these two different models. Basin‐wide unconformities, growthstrata and changes in provenance indicate deposition of the QLC Formation in a tectonically active basin. Both west‐ and east‐vergent structures, rooted in the basement, controlled the deposition and distribution of the QLC Formation from the Middle Eocene to the Early Miocene. The provenance analysis indicates that the main source areas were basement blocks, like the Paleozoic Oire Eruptive Complex, uplifted during Paleogene shortening, and that delimits the eastern boundary of the present‐day intraorogenic Puna plateau. A comparison of the QLC sedimentary basin‐fill pattern with those of adjacent Paleogene basins in the Puna plateau and in the Santa Bárbara System highlights the presence of discrete depozones. These reflect the early compartmentalization of the foreland, rather than a stepwise advance of the deformation front of a thrust belt. The early Tertiary foreland of the southern central Andes is represented by a ca. 250‐km‐wide area comprising several deformation zones (Arizaro, Macón, Copalayo and Calchaquí) in which doubly vergent or asymmetric structures, rooted in the basement, were generated. Hence, classical foreland model is difficult to apply in this Paleogene basin; and our data and interpretation agree with a broken‐foreland model.