Abstract

Understanding the history, causes, and impact of sea-level changes is a challenge for our societies that face accelerated global sea-level rise. In this context, improvement of our knowledge of sea-level changes and shoreline migration at geological time scales is critical. The preserved, laterally correlative sedimentary record of continental erosion on passive margins has been used to reconstruct past sea level. However, the detailed nature of a basic clinothem progradational pattern observed on many of these margins is still poorly known. This paper describes the sedimentary facies and interprets the depositional environments and the architecture of the clinothems of the New Jersey shelf (offshore northeastern USA) to depict the origin and controls of the distribution of the sediment on the margin. We analyze 612 cores totaling 1311 m in length collected at three sites 60 km offshore Atlantic City, New Jersey, during International Ocean Discovery Program–International Continental Scientific Drilling Program (IODP-ICDP) Expedition 313. The three sites sampled the lower to middle Miocene passive margin sediments of the New Jersey shelf clinothems. We also collected wireline logs at the three sites and tied the sedimentary architecture to the geometry observed on seismic profiles. The observed sediment distribution in the clinoform complex differs from that of current models based on seismic data, which predict a progressive increase in mud and decrease in sand contents in a seaward direction. In contrast, we observe that the clinoforms are largely composed of muds, with sands and coarser material concentrated at the rollover, the bottomset, and the toe of the slope. The shelf clinothem topsets are storm-influenced mud whereas the foreset slope is composed of a mud wedge largely dominated by density current deposits (e.g., low-density turbidites and debrites). The architecture of the clinothem complex includes a composite stack of ∼30-m-thick clinothem units each made up of four systems tracts (Transgressive, Highstand, Forced-Regressive, and Lowstand Systems Tract) building individual transgressive-regressive sequences. The presence of mud-rich facies deposited during highstands on the topset of the clinoform, 40–60 km offshore from the sand-prone shoreface deposit (observed in the New Jersey onshore delta plain), and the lack of subaerial erosion (and continental depositional environments) point to a depositional model involving a subaerial delta (onshore) feeding a distant subaqueous delta. During forced regressions, shelf-edge deltas periodically overstep the stacks of flood-influenced, offshore-marine mud wedges of the New Jersey subaqueous delta, bringing sand to the rollover and building up the large-scale shelf-prism clinothems. The clinothem complex develops on a gently dipping platform with a ramp-like morphology (apparent dip of 0.75°–0.5°) below mean storm wave base, in 30–50 m of water depth, 40–60 km seaward of the coastal area. Its shape depends on the balance between accommodation and sedimentation rates. Subaqueous deltas show higher accumulation rates than their subaerial counterparts and prograde three times further and faster than their contemporaneous shoreline. The increase in the intensity of waves (height and recurrence intervals) favors the separation between subaqueous and subaerial deltas, and as a consequence, the formation of a flat topset geometry, a decrease in flood events and fluvial discharge, an overall progressive decrease in sediment grain size (from sequence m5.45, ca. 17.8–17.7 Ma, onwards), as well as an increase in sedimentation rates on the foresets of the clinoforms. All of these are recognized as preliminary signals that might characterize the entry into the Neogene icehouse world.

INTRODUCTION

Passive margin successions are commonly characterized as relatively simple suites of sedimentary strata truncated by unconformities that can be correlated regionally, and in some cases worldwide. They show a gradational evolution of depositional environments from continental to deep marine realms, including the shoreline, and a well-known subsidence history. As such, passive margins have for a long time been used to reconstruct global sea-level variations on geological time scales (e.g., Vail et al., 1977; Miller et al., 2005). Sequence stratigraphy proved to be an effective tool to decipher the passive margin sedimentary record from the earliest outcomes of seismic stratigraphy (Vail et al., 1977; Posamentier et al., 1988) to the most recent standardization (e.g., Catuneanu et al., 2009), applications (e.g., Embry, 2009), and simplifications (Neal and Abreu, 2009; Miller et al., 2018). Numerous hypotheses, however, have not been fully tested yet, e.g., the nature of the sedimentary facies that compose the prograding clinothems on the shelf—a classic end member of most passive margin sedimentary records (e.g., Mitchum et al., 1977; Berg, 1982; Alexander et al., 1991; Pirmez et al., 1998; Hubbard et al., 2010, Helland-Hansen et al., 2012). Similarly, the timing and phase relationship of these sedimentary facies with respect to relative sea-level changes (Reynolds et al., 1991) or the paleo–water depth of sediment deposited at the top and the toe of the clinoforms (Greenlee and Moore, 1988) have not been determined definitively. Clinoforms are distinct sigmoidal geometric features associated with topset, foreset, and bottomset deposits that generally prograde seaward (termed a clinothem). Clinothem topsets were originally termed as the shelf and the rollover point at the shelf break (Vail et al., 1977). This has created confusion because the modern continental shelf-slope break is typically in 120–200 m of water (e.g., Heezen et al., 1959), and it has been shown that the rollover features (also called “depositional shelf breaks”) associated with clinothems are shallower, showing different features than the shelf-slope break.

This paper seeks to ground-truth the vertical succession and lateral facies associations of clinothems in the relatively simple passive margin system of the New Jersey shelf (offshore northeastern USA) by using the data provided by International Ocean Discovery Project–International Continental Scientific Drilling Program (IODP-ICDP) Expedition 313 (Mountain et al., 2010a, 2010b). We examine successively the sedimentary facies found in the coreholes, and the seismic and sedimentary architecture of the clinothems; we propose a stratigraphic and depositional model for the New Jersey clinothems that combines subaerial-subaqueous delta with shelf prism–scale clinoforms and discuss some important factors controlling this complex architecture.

BACKGROUND WORK

The New Jersey shelf is part of the U.S. mid-Atlantic margin that extends from New Jersey through Delaware to Maryland (Fig. 1). The U.S. Atlantic margin is a classic passive margin, which showed a rifting phase in the Late Triassic (ca. 230–198 Ma; Sheridan and Grow, 1988; Withjack et al., 1998) and seafloor spreading commencing during the Early to Middle Jurassic (180–165 Ma), followed by thermal subsidence, sediment loading, and flexure (Watts and Steckler, 1979; Reynolds et al., 1991). The sedimentary succession consists of Upper Triassic alluvial, evaporitic, and restricted marine sediments injected by dikes, sills, and lava flows, overlain by a thick (8–12 km) Jurassic to mid-Cretaceous limestone and shale succession fringed by a barrier reef complex (Poag, 1985a). Sedimentation rates decreased from the Late Cretaceous to the Paleogene, building up a mixed siliciclastic carbonate ramp ending in a condensed and starved clay-rich carbonate ramp with the general cooling of temperatures, by the late middle Eocene (onshore) to earliest Oligocene (offshore) (Miller and Snyder, 1997; Steckler et al., 1999). Due to the increase of tectonically and/or climatically driven denudation in the hinterland (Poag and Sevon, 1989; Poag, 1992; Poag and Ward, 1993), sediment supply increased in the late Oligocene and Miocene, building out progressively a large set of prograding clinothems (Fulthorpe and Austin, 1998) that were capped by Pleistocene deposits (Davies et al., 1992; Austin et al., 1995, 1996).

The New Jersey shelf prograding clinothems were first recognized with (low-resolution) multi-channel seismic profiles collected by Grow et al. (1979) and Schlee (1981) to image rift-stage sediments. Later studies outlined Paleogene and Neogene sequences (e.g., Poag and Schlee, 1984; Poag, 1985b; Poag and Ward, 1987; Greenlee and Moore, 1988) and, by using higher-resolution seismic data (from cruise Ew9009 [http://www-udc.ig.utexas.edu/sdc/cruise.php?cruiseIn=ew9009] on the R/V Ewing in 1990, cruise Oc270 [http://www-udc.ig.utexas.edu/sdc/cruise.php?cruiseIn=oc270] on the R/V Oceanus in 1995, and cruise CH0698 [http://www-udc.ig.utexas.edu/sdc/cruise.php?cruiseIn=ch0698] on the R/V Cape Hatteras in 1998; Fig. 1), the three-dimensional (3-D) geometry and lateral variability of the Miocene depocenters (Fulthorpe and Austin, 1998; Fulthorpe et al., 1999; Poulsen et al., 1998; Monteverde et al., 2008) tied to the available industry and Ocean Drilling Program (ODP) coreholes drilled on the coastal plain (Legs 174AX and 150X) and on the outer shelf and slope (Legs 174A and 150) (Fig. 2A).

These studies have provided the chronology of sedimentation on the New Jersey shelf for the last 100 m.y. (e.g., Poag, 1985b; Miller et al., 1998, 2005; Browning et al., 2006, 2013). The sequence boundaries recognized onshore from facies successions are correlated to the continental slope with a ±0.5 m.y. accuracy and tied to the glacio-eustatic lowerings associated with δ18O increases (Miller et al., 1998). The timing and the number of Miocene sequences fit with the Bahamas carbonate sequences (ODP Leg 166; Eberli et al., 1997) and the published Exxon Production Research global sea-level curve (Vail et al., 1977; Haq et al., 1987), although they significantly differ in amplitudes and rates (Van Sickel et al., 2004; Miller et al., 2005).

Despite the real progress presented by these studies, coastal plain cores missed most of the lowstand deposits due to changes in sea level of low to intermediate amplitudes that induce large unconformities and hiatuses, while the coreholes from the shelf edge and slope either failed to sample the sandy intervals (ODP Leg 174A) or provided a muddy homogeneous facies succession (ODP Leg 150), forming a poor basis for detailed inference of sea-level changes (Fig. 2). The IODP-ICDP Expedition 313 was specifically designed to fill this gap of knowledge and complement the New Jersey U.S. margin transect by drilling three holes in an intermediate position in the inner part of the New Jersey shelf in order to (1) recover most of the early to mid-Miocene interval in sea-level history including that inferred from the lowstand deposits and (2) provide a large and diversified set of sedimentary facies very sensitive to changes in sea level (Mountain et al., 2009) that can be easily interpreted in terms of paleo–water depths. Three holes were drilled in the inner shelf targeting the upper Oligocene to middle Miocene seismically imaged prograding clinothems (Fig. 2B). The three holes drilled in 34–36 m of water, 45–67 km offshore, sample a 22-km-long clinothem transect including the topset, foreset, and toeset of several clinothems from 180 to 750 m below seafloor (mbsf).

Building on seismic stratigraphy of Monteverde et al. (2008), Expedition 313 drilled through 25 regionally mapped Oligocene to Miocene seismic surfaces that correlate to facies changes in the coreholes (Mountain et al., 2010; Inwood et al., 2013; Miller et al., 2013a, 2013b; Browning et al., 2013). The lithostratigraphic description of split cores shows silt-rich supply systems that reveal a notable depletion in clays and a marked difference between the top and the toe of the clinoform bodies (e.g., Lofi et al., 2013). The topset facies succession shows well-sorted silts and sands deposited in offshore to shoreface, mixed wave- to river-dominated shelf environments. The toe of the clinothems show silts and silty clays deposited below wave base. These sediments are typically interbedded with poorly sorted silts and sands deposited by continuous downslope gravity transport processes such as debris flows and turbidity currents during periods of clinoform rollover and upper slope degradation (Hodgson et al., 2018). The open shelf experienced frequent cycles of dysoxia. In situ and reworked glauconite is a common component of topset and bottomset strata that also show sharp changes in pore water salinity. Strontium isotopic ages measured on molluscs and foraminifers, reliable biostratigraphic zonation of multiple fossil groups (foraminifers, dinocysts and nannofossils; Katz et al., 2013; Browning et al., 2013), and specific pollen markers (McCarthy et al., 2013; Kotthoff et al., 2014) verify a nearly continuous record of ∼1 m.y. sea-level cycles and climate variations that may explain facies changes along the slopes of the clinoforms. We found no evidence of sea-level drop below the clinoform inflection point—i.e., depositional shelf break or offlap break—but the occurrence of shoreface deposits along the slope of the clinoforms and of deep water facies on their topsets suggest large changes in amplitude of relative sea level in the range of 60 m (Mountain et al., 2010a, 2010b).

DATA AND METHOD

Despite the difficulties of coring the sandy material of the shallow New Jersey shelf, 612 cores were collected at three sites (Sites M0027, M0028, and M0029) with 80% recovery for a total of 1311 m (Mountain et al., 2010a, 2010b). The deepest hole (M0029A) reached 757 mbsf; the oldest sediment (uppermost Eocene) was recovered in Hole M0027A. Besides the cores, the expedition collected wireline logs at the three sites—gamma ray, resistivity, magnetic susceptibility, sonic, acoustic televiewer, and vertical seismic profiles—which, together with multisensor core logs on unsplit cores, provide precise ties between core logs and seismic profiles (Mountain et al., 2010a, 2010b; Miller et al., 2013b). This data set was complemented post-cruise by the shipboard party by visual core descriptions and smear slide analyses, biostratigraphy (calcareous nannoplankton, diatoms, dinoflagellate cysts), magnetostratigraphy, and pore water geochemistry; paleobathymetry and paleoenvironments inferred from benthic foraminifera assemblages, dinocysts, pollen analyses, and grain size analysis (Ando et al., 2014). Multi-sensor core logger (MSCL) and petrophysical data were directly measured on cores and core sections (Inwood et al., 2013); Sr analyses were performed as an aid to get additional age control on the sediment and to evaluate the length of time represented by key stratal surfaces (Miller et al., 2013b; Browning et al., 2013).

Here we present a detailed analysis of the lithofacies successions and their interpretation in terms of depositional environment based on visual core description. Genetically related facies successions are bounded by unconformities tied to the interpreted seismic line 529 (cruise Oc270) (Monteverde et al., 2008) crossing the three holes, by using MSCL and downhole (impedance) log data (Miller et al., 2013b). In the following, we use the updated version of the stratal ages of Mountain et al. (2010b) provided by Browning et al. (2013).

RESULTS

The most striking feature of the early to mid-Miocene sediments of the New Jersey passive margin, as seen on the seismic line 529 of cruise Oc270, is the clinoform shape of the reflectors (Fig. 2). Here we present a description of the main lithofacies of the clinothems sampled in Site M0027, M0028, and M0029 coreholes of Expedition 313 and an interpretation of their depositional environments (Fig. 3; Supplemental Files [Figs. S1–S3]1). Based on the established ties between the cores and the Oc529 seismic line (Mountain et al., 2010b; Miller et al., 2013a, 2013b), we reconstruct the lateral facies distribution along a two-dimensional (2-D) dip-oriented line across the margin, and interpret their depositional environments, systems tracts, and sequences.

Lithofacies and Depositional Environments

The lithofacies observed in cores are described below from the deep offshore marine to the shallow marine shoreface environments (Table 1). Numbers in the text refer to the hole and the lithofacies label (e.g., 27-1 refers to Hole M0027 and lithofacies 1; 28-4 refers to Hole M0028 and lithofacies 4; etc.) shown in Figure 3 and with more details in the Supplemental Files (Figs. S1–S3 [footnote 1]).

Offshore Marine Gravity-Flow Deposits—Lithofacies Association 1 (FA1)

Gravity flow deposits are abundant at the base of the three holes (M0027, 480–625 m; M0028, 500–610 m; M0029, 600–730 m; Fig. 3) and in the middle part of Hole M0029 (450–500 m; Fig. 3). Their lithological diversity is very large. They are here tentatively grouped in two broad categories, but a detailed description of their occurrence along the bottomsets of the clinothems can be found in Hodgson et al. (2018).

Amalgamated toe-of-slope debrites and turbidites—FA1.1. Lithofacies association FA1.1 occurs above seismic surfaces at the bottomset of clinoforms (Fig. 3). The lithofacies comprises a few-meters-thick, poorly sorted, glauconitic and quartzitic, medium to coarse sand with angular quartz granules and gravels, clay clasts, and, in places, a muddy matrix. Glauconite content can decrease upsection from 80%–90% (lithofacies 27-4, 28-6) to <40% (28-7), or to 5% with increasing mud content and decreasing quartz grain concentration (29-15). At the base, the facies is homogeneous and well bedded (27-4, 27-5, 28-6, 28-2), alternating with pale brown clay beds with nannofossils, locally cemented by carbonates and bioturbated by clay-lined horizontal and vertical burrows filled by glauconite (up to 87 cm below erosion surfaces). The gravel to medium sand beds show crude normal and reverse grading, with 1-m-high oblique stratification (15° apparent dip) underlined by muddy lamina cut by clay-filled Chondrites trace fossils (27-5, 28-6) and ripple and low-angle lamination. Upsection, the facies thins and fines upward (rarely coarsening), organized in meter-scale successions (28-7, 27-6, 28-9, 29-9, 29-13, 29-15, 29-18) capped by bioturbated horizons. Parallel and ripple cross-lamination give rise to finer-grained sediment with faint low-energy horizontal lamination, convolution, and laminae with a chaotic pattern. Thick- and thin-walled shell fragments are scattered throughout but show signs of severe dissolution. Locally, concretions and quartz granules have a patchy green coating or shattered fabric (e.g., 29-15). No plant debris or mica is observed.

The basal graded beds with angular quartz and rounded glauconite in inclined beds and the well-stratified succession with nannofossils in clay lamina suggest a sediment supply via sediment gravity flows (turbidity currents and debris flows) in a deep water environment (well below wave base). The apparent cross-bedding may indicate a channelized context where turbidity currents built dune-scale migrating bedforms on the toe-of-clinoform slope. Burrows and clay laminae show that the dune-scale bedforms migrated episodically, rather than from sustained flows. The sedimentary structures in the quartz-rich glauconite sands point to some winnowing of sediment at the seafloor, which may have favored the authigenesis of glauconite grains, their relative concentration, and then possibly their cementation with carbonate along irregular erosion surfaces. Sediment was likely sourced from a clinoform rollover location that may have been exposed, as the subangular nature of the quartz sand grains point to a very early burying of freshly eroded material. This basal part of FA1.1 is interpreted as a basal channel lag, recording bypass along erosion surfaces (Stevenson et al., 2015).

The upper part of FA1.1 is well organized in meter-scale graded successions, but it is still poorly sorted, with a mix of quartz, glauconite, and shell fragments with traces of dissolution in a muddy matrix. These are probable indications of downslope gravity sediment transport. The fining-upward trend together with the suite of sedimentary structures showing a progressive decrease in the energy of transport (horizontal parallel lamination, cross-bedding, and low-energy planar lamination) point to a probable channelized transport of a stack of debrites and turbidites at the toe of a clinoform slope apron. The lack of mica and plant debris is noticeable and might indicate a progressive overall retreat of the regional river outlets. The gradational thinning and fining of beds uphole indicate a possible sediment supply response to transgression in an overall marine onlap. Those facies likely were emplaced below the maximum storm wave base (i.e., 80–120 m water depth) as suggested by the lack of traces of hydraulic reworking by waves.

Slope apron deposits—FA1.2. Lithofacies association FA1.2 comprises a wide range of facies that are usually arranged as glauconite-rich decameter-thick fining-upward successions. A distinction is made between coarse-grained deposits (FA1.2a) and the overlying silty-clay deposits (FA1.2b). FA1.2a is composed of stacked, meter- to 10-m-thick coarsening- and/or fining-upward successions. The coarsening-upward succession (lithofacies 27-2, 28-4, 28-5, 29-5, 29-8, 29-17) is made up of homogeneous, poorly sorted, coarse-grained quartz and glauconitic (5%–10% to 30%) sand with dispersed quartz granules, concentrations of pebbles (up to 15%), and common weak normal grading; granule and pebble grains tend to be less rounded than sand-grade grains. The fining-upward successions (27-3, 28-5, 29-3, 29-5, 29-17) are coarse- to medium-grained quartz and glauconitic (up to 25%–60%; 27-3) sand lying on a coarse basal lag with granules along an erosive or bioturbated surface. The original subhorizontal or ripple lamination with normal grading are in places largely obscured by bioturbation (e.g., Planolites, Teichichnus, Diplocraterion, and Zoophycos; 27-3) and dish, contorted, and convolute lamination. These coarsening-upward–fining-upward successions are locally interbedded with mud-supported granule and pebble sands with centimeter-thick clay beds (e.g., 29-3). Glauconite grains are a mix of pale-green (immature) granules and smaller black grains (mature). FA1.2b is fining-upward brown silt to silty clay with glauconitic, fine- to medium-grained sand beds, floating coarse and very coarse quartz and glauconite sand grains (29-10), foraminifers and thin-walled shells, traces of mica, and plant debris. Granules and aggregates of pale-green (immature) and subangular glauconite are observed with fine-sand grains of dark green and subrounded glauconite (mature). The sediment is moderately bioturbated with identifiable Chondrites and Planolites (Fig. 4C) and at one place a 30-cm-long Thalassinoides burrow (27-1). Glauconite is commonly concentrated in burrows (vertical burrows and Teichichnus) along erosion surfaces together with rare shell fragments, benthic foraminifers, and disseminated pyrite but including quartz granules (29-11, 29-12) or not (29-7). The sand beds form weakly normally or reversely graded layers with rare ripple-scale cross-lamination (29-4, 29-6). Thick-walled shells occur but show signs of severe dissolution (29-12).

FA1.2a facies are dominated by moderately to poorly sorted, coarse sand-prone sediment with dispersed granules and pebbles and local concentrations of gravel. Weak normal grading supports the partial transformation of debris flows into high-density turbidity currents, but the general poor sorting can be attributed to mixing through bioturbation and/or cohesive debris flow deposition, as shown by the presence of muddy sand with floating granules and pebbles. Centimeter-thick clay and sand interbeds might represent waning-stage deposits. These poorly sorted coarse sediments are interpreted as high-concentration flows of coarse material eroded from updip positions settled at the toe of a clinoform slope apron. These facies may correspond to simple coalesced toe-of-slope fans. The presence of micas and plant debris in the finer FA1.2b is indicative of a distant but quasi-permanent river sediment input on the shelf. The abundance, in places exclusive (29-7), of glauconite grains—either (1) dark green, well-rounded and subangular, and broken, or (2) pale green, in situ aggregates—as well as thick shells with traces of dissolution (29-12), points to temporary local deep marine starved conditions. The common concentration of the glauconite grains, together with foraminifers or quartz sand grains, in burrows along erosion surfaces points to occasional current-swept floor conditions, which probably in turn favors the generation of the glauconite. The glauconite grains can be produced both in situ or reworked from the clinoform slope (e.g., Hesselbo and Huggett, 2001). The common occurrence of fining-upward, centimeter- to decimeter-thick sand beds (29-10, 29-12), or thick mud (29-7) with floating coarse (29-10) to granule-size (29-12) quartz and glauconite grains clearly indicates either a component of downslope transport of clastic sediment in dilute turbidity currents or a stack of poorly mature, cohesive(?) debris flows provided by the slope of a clinoform. Overall the decameter-thick vertical stack of FA1.2a and FA1.2b corresponds to the development of fans on the bottomsets of clinoforms that develop below the maximum storm wave base (i.e., 80–120 m water depth) as suggested by the lack of traces of hydraulic reworking by waves.

Storm- and River-Influenced Offshore Marine Deposits—FA2

Offshore marine deposits consisting of marine silts and silty clays are widely distributed in the three holes. They encompass a large number of facies deposited below storm wave base (lower offshore, FA2.1) and between fair-weather and storm wave base (upper offshore, FA2.2). In the literature, this “upper offshore” domain is often termed the “offshore-shoreface transition” (as used by Mountain et al. [2010b]). In cores, lower (FA2.1) and upper (FA2.2) offshore deposits are either river- (FA2.1a, 2.2a) or storm- (FA2.1b, 2.2b) influenced.

River flood–influenced deposits. The river-influenced deposits consist of coarsening- and thickening-upward meter- to 10-m-thick successions of silty clay to fine to medium quartz sand with abundant mica, plant debris, shell fragments (lithofacies 28-3, 28-12), sponge spicules, and diatoms. The silty clays (FA2.1a) are nannofossil bearing. They can be tan colored (29-2) showing micro-laminated clay-rich layers (28-12). Sandier beds are parallel or ripple laminated, and normally (28-3, 28-10, 28-12) to inversely graded (28-3, 28-10) with gradational boundaries (28-3). Bioturbation is absent in clay-rich layers but obscures laminations in sands. Identified frequent burrows include Chondrites and Planolites but are predominantly simple meniscate backfilled forms, e.g., Teichichnus (29-2) or Taenidium (Fig. 4A), filled by glauconite (28-3) or replaced by pyrite that may mimic fecal pellets (?Ophiomorpha). Also present are common precompactional carbonate concretions (29-2, 28-3, 27-7, 28-3). The silt to fine sands (FA2.2a) show abundant plant debris, scattered pebbles and granules (28-11), parallel and climbing ripple lamination, and centimeter-scale reverse and normal grading, with low-angle lamination above a swaley scoured surface.

FA2.1a lacks abundant colloidal organic matter, but the abundance of plant debris and mica indicates a fluvial source for the sediment. The well-laminated silty clays suggest that the depositional environment was probably very quiet and periodically subject to dysoxic bottom waters, as shown by the presence of concretions that excluded burrowing fauna. At times of better oxygenation, the infauna was dominated by horizontally mining deposit-feeding organisms. The gradational upper and lower boundaries and local basal reverse grading of the very fine, parallel- and ripple-laminated sand beds, as well as the millimeter-scale clay-silt laminations in the clay, point to a river flood origin of deposition with a possible seasonality. The depositional environment is interpreted as river-dominated lower offshore (prodelta) settings. FA2.2a shows indications of hydraulic sediment transport on the seafloor (parallel, climbing ripple lamination) and also the possible distal influence of storms (low-angle lamination on swaley scoured surfaces) indicative of a progressive shallowing from the lower to the upper offshore environments.

Low-energy, storm-influenced deposits. The storm-influenced deposits consist of coarsening-upward, meter- to 10-m-thick successions of silty clay to quartz and glauconitic sands with gastropods (e.g., Turritella), thin-walled bivalve shell fragments, and benthic foraminifers. The clayey silt to silty clay (e.g., 27-16, 28-8, 28-21, 28-23) (FA2.1b) contains abundant dispersed glauconite sand grains and macroscopic plant debris and mica grains, landward of clinoform rollover points (e.g., 27-16, 27-17). In the upper part of the three holes is a dark brown massive clay with abundant benthic foraminifers, rare shell fragments, faint lamination defined by weak concentrations of plant debris, bioturbation (Planolites), and scattered calcilutite beds (27-17), and a pale and gray-yellow-brown color-banded clay with Chondrites, pyritic and silty laminae, and common soft-sediment contorted bedding at different scales (28-31, 29-26). Sporadic millimeter- to centimeter-scale silts and very fine sand beds with faint laminae and a fining-upward trend are locally present within the silty clay (e.g., 27-19, 29-16). Bioturbation is usually pervasive (29-14, 29-16, 29-21), and burrow infills locally contain coarse glauconite and quartz sand or granules (e.g., 29-16). The thoroughly bioturbated silty clays are interbedded with coarsening- (28-25) or fining-upward (28-22, 29-22, 29-27) successions of centimeter-scale, silty to very fine sand beds (FA2.2b). Bioturbation is pervasive in silty clay beds and also underlines erosion surfaces (Chondrites, Planolites, Cylindrichnus, Teichichnus burrows; Fig. 4B) (27-15, 29-23) or crosscuts the lamination (28-22). Sand beds are normally graded (29-22, 29-23) with low-angle shell layers (28-25, 29-20), hummocky cross-stratification (27-15, 28-22), and ripple and low-angle cross-laminations (28-29). FA2.2b typically overlies FA2.1b; together they form meter- to decimeter-thick coarsening-upward stacks of facies.

The silty clay material of facies FA2.1b shows indications of deep lower offshore marine conditions with an abundant autochthonous fauna and active bioturbation, and episodic starvation (in situ glauconite), but no trace of hydraulic reworking by wave action. In places, the banded character of the clays with Chondrites and pyrite and the presence of calcilutite beds confirm periodic sediment-starved bottom conditions. Color banding in clays underlines discrete depositional events that are interpreted as low-density turbidity currents possibly triggered by storm events. Occasional downslope sediment transport is suggested by the presence of sporadic fining-upward centimeter-thick sand beds. Deformation of cohesive sediment indicates some instability of the seafloor, and a clastic sill near the base of the clays may relate to overpressure in underlying sands. The mixing of mud with sand-size glauconite grains might imply some reworking from upslope of poorly mature, cohesive(?) debris flow deposits with burrowed erosive surfaces indicative of exposure of the seafloor between two successive sedimentation events. The heterolithic nature of FA2.2b points to upper offshore conditions with alternating quiet environment and reworking or downslope transport processes. The occurrence of well-sorted, fining-upward sand beds and shell layers with low-angle and ripple laminations are indicative of hydraulic reworking of the seafloor by the distal influence of storm waves (i.e., storm-graded layers) in slightly shallower waters than FA2.1b in upper offshore environment. The meter-scale coarsening-upward trends point to probably recurrent periods of shoaling from the lower to the upper offshore environments. Paleoenvironmental reconstructions suggest that FA2.2b facies were deposited below maximum fair-weather wave base (i.e., 20–30 m water depth) and above maximum storm wave base (i.e., 80–120 m water depth).

Lower Shoreface Deposits—FA3

The lower shoreface facies form regular coarsening-upward progradational successions at the toe of the shoreface. They are of two types, depending on their content in plant debris and mica (FA3.1) and in glauconite (FA3.2).

FA3.1 is rich in coarse sands, micas, and macroscopic plant debris. It is located in the middle part of Hole M0027 (332–337 m, 375–415 m) and the lower part of Hole M0028 (400–440 m). It shows coarsening- and thickening-upward successions of moderately bioturbated silts to well-sorted medium quartz sand with granules (lithofacies 27-10, 28-13, 28-14, 28-27). There are many distinctive two-part sand beds in the lower part showing a slight grain-size change at the lamina scale (28-13). These beds have sharp bases overlain by clean quartz sand (no silt or mica) and silty sand with subparallel low-angle lamination (Fig. 4D), current ripple laminations, and symmetrical-aggradational ripples at the top. Shells and plant debris are concentrated near the tops of beds (Fig. 4D), which are typically bioturbated (28-13). Shell hash layers increase upward along low-angle slightly erosive soles (28-14) or thick-thin beds (27-10). In the latter, thick beds with low-angle oblique lamination or hummocky cross-stratification, capped by mud laminae, alternate with intensely bioturbated horizons. Less-bioturbated places show some distinctive burrows such as Cylindrichnus, Thalassinoides, and Teichichnus. FA3.2 is rich in glauconite grains. It is found in the upper part of Holes M0028 and M0029. It comprises coarsening- and thickening-upward, meter- to decimeter-scale successions from muddy, glauconitic silt to medium quartz sand (28-17, 28-19). Hints of laminations are delineated by shell lags (28-17–28-19).

The abundance of mica and plant debris in FA3.1 suggests that the supply had a strong fluvial influence. The two-part sand beds are interpreted as combined-flow storm deposits in a high-sediment-supply environment. The lower clean and slightly coarser quartz sand is interpreted as having been deposited from geostrophic (unidirectional) currents that transported sediment entrained from the coeval shoreface to deeper water, as is commonly observed in subaqueous deltas (Mitchell et al., 2012; Patruno et al., 2015b). The upper sediment was locally agitated, and deposits from (orbital) currents set up by storm waves formed low-amplitude hummocks, current ripples, and late-stage symmetrical (wave) ripples before organisms could start to mix the sediment again. The succession overlies upper offshore FA2.1a and FA2.2a stacks and is interpreted as having been deposited at the toe of the shoreface in a river-influenced setting at water depth ranging from 10 to 20–30 m. The presence of glauconite grains and the overall absence of plant debris in FA3.2 suggest reduced river sediment supply as compared with FA3.1.

Upper Shoreface Sand Deposits—FA4

Upper shoreface facies comprises quartz-rich sand successions (FA4.1) overlying facies FA3.1, and glauconite-rich sand (FA4.2) overlying FA3.2.

FA4.1 is thick-bedded, coarsening-upward, medium- to coarse-grained quartz sands with subangular granules. The sediment is poorly sorted. It includes mica, wood and plant fragments (lithofacies 27-11, 27-12, 29-29), subangular granules (28-15, 28-28), black non-organic grains, and dispersed shell debris. The sand usually appears homogeneous and loose with remnants of deeply scoured (28-15, 29-29) or cemented channelized erosion surfaces (27-12). Sedimentary structures are faint and, when preserved, show high-angle trough cross-bedding and low-angle laminations (27-11, 27-12). Bioturbation is represented by distinctive Cylindrichnus backfilled burrows (27-11, 27-12). In the upper part of Hole M0029 (29-29, 29-30, 29-31), the stack of coarsening-upward trends is capped by interbedded, centimeter-scale, sharp-based, bioturbated blue-gray silty clay and well-sorted, medium to fine sand showing in places ripple cross-lamination (Fig. 4E). FA4.2 is thick bedded, poorly sorted, coarsening-upward, medium- to very coarse-grained glauconite-rich (up to 20%) quartz sand with granules (27-13) and dispersed shell debris (28-18). Glauconite content tends to increase upward from 1%–3% to as much as 20%. Sedimentary structures are usually poorly preserved and the sand appears homogeneous and loose (28-18), although sparse high-angle cross-bedding and scoured surfaces are locally underlined by coarse quartz sand laminae (27-13). Sparse vertical burrows and local carbonate-cemented horizons are recorded along internal discontinuities.

The coarse-grained, thick-bedded coarse sand with subangular granules and high-energy sedimentary structures (high- and low-angle cross-bedding, scouring) of both facies associations are typical of prograding wave-influenced shoreface sediment bodies. The abundance of mica flakes, angular granules, and wood fragments (FA4.1) points to the proximity of a river feeder at the top of the shoreface or delta lobe sands as shown by the presence of channelized scours (e.g., 27-12) in Hole M0027. The more massive and loose character of sand of FA4.2 together with high glauconite content might indicate the deposition of sediment rapidly transported across (bypassing?) the shelf from outer estuarine environments to the rollover. These sands deposited below fair-weather wave base at the rollover recall the sand-prone subaqueous clinoform facies of Mitchell (2012), Mitchell et al. (2012), and Patruno et al. (2015b).

Gully and Clinoform-Edge Fan Deposits—FA5

Facies association FA5 underlines seismic unconformities m5.2 (ca. 15–15.6 Ma; Browning et al., 2013) and m5 (ca. 13.7 Ma) (Fig. 2) landward of clinoform rollover points and overlies FA4.2 facies in Hole M0028. It comprises crudely coarsening-upward, coarse glauconite sand (up to 40%) intercalated with clay beds (lithofacies 28-26) or lamina (28-20). FA5 is poorly sorted and dominantly coarse grained, with a mix of mud, sand, and gravel. Generally, sediment is low in mica and organic matter and deficient in silt to medium sand grain sizes compared to the surrounding stratigraphy. Overall, the succession becomes more glauconitic (up to 40%) and increasingly stratified uphole, with highly bioturbated mud-rich units with gravel layers or weakly graded beds. FA5 shows numerous internal erosion surfaces overlain by (1) weakly fining-upward and bioturbated gravel-rich sand and sandstone with complete benthic foraminifera (28-26), and (2) 0.1–0.2-m-thick matrix-supported, muddy glauconitic (up to 40%) poorly sorted coarse sand with up to 20% quartz completed by lithic and mud granules and pebbles, and scattered pyrite grains, lacking shells but containing molds. Granules and pebbles are well rounded to angular. Bioturbation is pervasive, with mud- and glauconite-filled burrows (28-20). Mud pebbles exhibit bioturbation identical to that of the host bed (28-20).

The poorly sorted sand with mud and gravel, with shell molds and rip-up clasts, indicates erosion and reworking in a marine setting. Rounded clay clasts that have a bioturbated fabric are interpreted as intraclasts entrained from the substrate. The weak grading of some deposits supports waning sediment gravity flows. Possible interpretations of these deposits at the clinoform rollover include gully fills for the fining-upward successions and/or outbuilding of small deltas with steep fronts for the coarsening-upward trends that formed during periods of lowering relative sea level. The lack of plant debris and mica and the low proportion of very fine to medium sand compared to the surrounding stratigraphy suggest that these components may have been transported downdip, with the coarser material and muds from eroded surrounding substrate preferentially retained. The extent of the drainage basin that fed the deltas and/or gullies is not clear. However, the source for the coarse extrabasinal material need not to be tied directly to the hinterland and may have been from older deposits.

Transgressive Lag Deposits—FA6

Facies association FA6 lies on erosion surfaces at the topset of clinoforms and underlines seismic surfaces m5 (ca. 13.7–14.8 Ma; Browning et al., 2013), m5.2 (ca. 15–15.6 Ma), m5.3 (ca. 16.3–17.3 Ma), m5.32 (ca. 16.7–17.4 Ma), m5.4 (ca. 17–17.7 Ma), and m5.45 (ca. 18 Ma) in Hole M0027, and m4.1 (ca. 12.9 Ma) in M0028 (Figs. 2 and 3). FA6 is composed of massive, sharp-based, decimeter- to meter-scale, usually fining-upward successions from (1) clean, coarse to very coarse glauconitic and quartz sand with subangular quartz and glauconite granules to (2) bioturbated fine sand, containing moderately sorted medium sand, with large pieces of wood and plant debris, large thick-walled shells, mica (2%), benthic foraminifers, and clay-filled burrows, and (3) clayey silts with dispersed organic matter, isolated sand grains and granules, common shell layers (bivalves and gastropods), foraminifers, diatoms, and sand-filled burrows (lithofacies 27-18). Bored cemented sand nodules with glauconitized rims and Thalassinoides burrows going down to 25 cm filled by coarse sand and granules are observed, in places, along erosion surfaces (27-14, 28-16, 28-30).

The coarse-grained sediment with large, thick-walled shell and wood debris originates from a very proximal and high-energy shoreface environment close to a river outlet. The poor sorting with a mixing of granules with fines, arranged in both fining- and coarsening- upward trends, point to a transport mechanism by a mix of sediment traction and gravity flows along the erosional surfaces. The lack of clayey material suggests some reworking by wave energy in a depositional environment close to the toe of the shoreface. Coarse to fine sand with bored nodules and glauconite grains are indicative of condensed deposition and subsequent erosion and exhumation. These sediments underlining erosion surfaces are interpreted as relics of ravinement deposits associated with sea-level rise and transgression.

Seismic Stratigraphic Architecture

The interpretation of the New Jersey shelf seismic data showed that the early and middle Miocene sediment pile can be divided into 18 seismic units or clinothems bounded by unconformities of regional lateral extent designated by the name of their underlying unconformity (m6–m4.1; Monteverde et al., 2008) (Figs. 2 and 5A). In their regional study, Monteverde et al. (2008) distinguished two packages of reflectors on the New Jersey shelf, interpreted as lowstand and highstand systems tracts. The detailed observation of the 2-D seismic line 529 of cruise Oc270 led us to individualize four sets of reflection packages in the seismic units (Figs. 5, 6, 7, and 8) based on their shape and the internal reflector terminations and configurations. From base to top, the sets are the following:

  • Set 1—A lens-shaped unit bounded at the base and top by high-amplitude reflectors onlapping the toe of the clinoform front. The set is transparent to reflection free.

  • Set 2—A wedge- to lens-shaped unit largely sitting on the shelf with high-amplitude reflectors at the base and top, and more locally (distally) sitting on the bottomset or foreset of clinoforms. Internal reflectors show onlap and downlap basal terminations and toplap terminations at the top. Reflectors show high amplitude, good continuity, and low frequency on the shelf. The configuration is subparallel on the shelf, and low-angle oblique tangential seaward, in the clinothems.

  • Set 3—A slope-front fill unit with a bank to lens shape. Reflectors show moderate-angle downlap terminations and local onlaps, or are conformable at base, and show erosional or apparent truncations at top. Internal reflectors show a mix of high- and low-amplitude reflections, low to moderate continuity, and low frequency. The configuration is complex oblique to chaotic.

  • Set 4—A slope-front fill unit with a lens shape. Reflectors exhibit onlap-downlap terminations, locally conformable at the base and with apparent truncation at the top. Internal reflectors are of moderate to high amplitude and largely discontinuous, with a moderate frequency. The configuration is tangential oblique to complex oblique sigmoidal.

Sedimentary Facies Architecture

Preliminary ties were established between the cores, the holes (through downhole logs), and the Oc270–529 seismic line (Fig. 2) during Expedition 313 (Mountain et al., 2010b) and refined by Miller et al. (2013b). A reappraisal of these ties according to sedimentological observations made in this study provide the opportunity to correlate corehole facies interpretations from one hole to another (Fig. 3) and along the regional 2-D dip section following geometries observed on the seismic line (Fig. 7).

Clinothem Units

The individual seismic units appear basically composed of a deepening (transgressive) and a shallowing (regressive) suite of facies made up of coarse-grained debrites overlain by a progradational stack of lower to upper offshore, and shoreface sediments (Fig. 7; see a sketch of a typical clinothem units in Fig. 8). In detail, however, each set of seismic reflectors in the seismic unit shows a specific suite of facies or systems tract in cores (Figs. 7 and 8). Seismic set 1 comprises amalgamated debrite deposits (lithofacies association FA1.1) at the toe of the slope and transgressive lag deposits (FA6) on the topset of clinoforms. Seismic set 2 is well represented on the shelf but rarely recognized at the toe of the slope. It shows lower offshore sediments (FA2.1) giving rise upsection to upper offshore deposits (FA2.2). Seismic set 3 forms the bulk of the clinothems. It shows a large suite of facies from slope apron fans (FA1.2), to lower (FA2.1) and upper offshore deposits (FA2.2), giving rise to lower shoreface sediments (FA3). Gully-fill deposits (FA5) are observed twice capping set 3 deposits landward of the rollover. Seismic set 4 exhibits a large range of facies and the shallowest depositional environments. It shapes the rollover and drapes the slope. It comprises lower offshore to upper shoreface deposits (FA2.1 to FA4).

Key Bounding Surfaces

The seismic sequence boundaries (e.g., m5.8–4.1; Fig. 7) are the most prominent surfaces that bound the clinothem units in the coreholes (see Miller et al. [2013b] for a detailed description). They are sharp, erosional surfaces separating shallow marine facies below from deeper marine facies above (see a sketch of a typical clinothem unit in Fig. 8 and details in Fig. 3 and the Supplemental Files [Figs. S1, S2, and S3; footnote 1]). On the shelf, these surfaces are underlain by poorly sorted, immature sand material (transgressive ravinement deposits, lithofacies association FA6). At the toe of the slope, they are overlain by amalgamated sediment gravity flows (FA1.1), which progressively wedge out further seaward where the surfaces become conformable. The seismic set boundaries within the seismic units can also be delineated in cores. The surfaces at the boundary between seismic sets 1 and 2 correspond to a change from fining-upward to coarsening-upward trends at the deepest paleoenvironments, the finest-grain-size horizon, and gamma-log peaks. The surfaces at the boundary between seismic sets 2 and 3 and sets 3 and 4 are less straightforward in cores than on seismic lines. A seaward shift of the facies (i.e., systems) tracts and, in places, a sharp increase in grain size explains the high amplitude of the reflectors and characterizes the changes in the angle of the reflectors between the seismic sets on seismic lines.

Stacks of Clinothem Units and Clinoform Trajectory Analysis

Despite the lack of 3-D control on geometries, the 2-D seismic section displays the gross sediment architecture of the stack of New Jersey margin clinothem units (Monteverde et al., 2008). As shown on Figure 7A, and following Helland-Hansen and Hampson (2009), the trajectory of the clinoform rollovers shows repeating features that characterize different stacks of clinothem units (stacks A to D, Fig. 7A). Each stack exhibits a landward shift of the clinoforms’ break in slope followed by a seaward shift, which correspond respectively to a retrogradational and a progradational stacking pattern of the clinothem units into a larger-scale transgressive-regressive (T-R) sequence. Each T-R sequence is underlain on the topset by an erosional surface at the maximum regression that corresponds to the transgressive surface of the classical stratigraphic sequences. This surface is overlain by transgressive ravinement deposits (lithofacies association FA6) on the shelf, glauconitic shoreface sands (FA3, FA4) at the rollover, and submarine fan debrites (FA1.1) of the healing phase at the toe of the clinoform slope. In places, the surface can go through the debrites facies (Hodgson et al., 2018). The early stage of the progradational part of the T-R sequence comprises offshore marine silts and clays enriched with mica and wood fragments sourced by rivers (FA2) while the late stage of the progradation shows relatively clean and sandier lower (FA3) to upper shoreface (FA4) material influenced by permanent wave action giving rise downslope to slope apron sediment gravity flows (FA1.2).

DISCUSSION

Clinoforms are a common structure of continental margins in a number of settings and a key component of sedimentary successions used to interpret the interplay between changes in relative sea level and sediment supply (Vail et al., 1977; Posamentier et al., 1988; Van Wagoner et al., 1990). A review of the literature shows that different scales of progressively large-scale clinoforms can prograde simultaneously along a transect from the shoreline to the abyssal plain, merging to form compound clinoforms, depending on changes in sediment supply, sea level, and subsidence (Helland-Hansen and Hampson, 2009; Henriksen et al., 2009; Swenson et al., 2005; Patruno et al., 2015a). Expedition 313 focused on the intermediate-scale clinoforms located on the shelf (Hodgson et al., 2018; Mountain et al., 2010a, 2010b), between (1) the (subaerial) shoreline delta and its subaqueous or prodelta counterpart, of a few meters to tens of meters high, formed by T-R processes (i.e., that can migrate landward during transgressions; e.g., Helland-Hansen et al., 2012; Jones et al., 2015); and (2) the shelf-slope structures, at continental-margins scale, of hundreds to thousands of meters high formed by long-term subsidence (Helland-Hansen et al., 2012; Jones et al., 2015; see Patruno et al. [2015c] for a recent review). Although the New Jersey “intermediate-scale” clinoforms could be referred to as “mid-shelf delta” (Porębski and Steel, 2006) or “shelf-edge” or “shelf-prism” clinoforms (Helland-Hansen and Gjelberg, 2012 ; Patruno et al., 2015c), they show a range of scales from the small “delta” scale of a few tens of meters to the large “shelf-prism” scale >100 m thick during their progressive growth (see below and, e.g., Figs. 5 and 7). Such different scales can be preserved on present-day margins with sufficient accommodation space, such as in the Mediterranean Sea (e.g., Rabineau et al., 2006, 2014). The “shelf-prism” clinoforms prograde and accrete basinward or remain fixed through time (Helland-Hansen and Hampson, 2009) at 104 to 106 yr time scales (Vanney and Stanley, 1983), but can hardly backstep.

The New Jersey clinothems sit on low-angle (<1°) Eocene carbonate ramp deposits controlled by differential subsidence that were starved of siliciclastic input until the late Oligocene (Steckler, 1999). Due to global climate cooling and the correlative increase in sediment supply, the clinoforms started to nucleate (by delta progradation) during the late Oligocene to early Miocene, progressively building up the present-day shelf margin. These clinoforms built by the progressive growth of small-scale prograding subaerial-subaqueous deltas merged into large-scale shelf-prism (sensuPatruno et al., 2015a) or shelf-edge (sensuHelland-Hansen and Hampson, 2009) clinoforms, building a compound clinoform (sensuJones et al., 2015; Patruno et al., 2015a). This large-scale geometry is controlled by a number of parameters, including long-term subsidence, global sea-level changes, sediment supply, and oceanic currents, among others (e.g., Patruno et al., 2015a). The New Jersey clinoform complex sits on an initial topography that was three to five times steeper than the present-day shelf profile (1:300–1:500 slope versus 1:1000 at present) (Steckler et al., 1999). The height of the clinoform rose progressively in a concave-up ascending trajectory, from 26 m in seismic unit m6 to a maximum of 250 m in m4 (46 m in m5.6, 116 m in m5.4, 131 m in m5) (Steckler et al., 1999), with a quite constant flooding water depth of 60–120 m at the toe of the clinoform (Katz et al., 2013) except at m4, the discussion of which is beyond the scope of this paper.

Characteristically, this intrashelf clinoform complex is largely composed of muds with <∼30% clean sands and coarse material transported at the rollover and beyond at the toe of the slope (see example and discussion in Poyatos-Moré et al. [2016]). The shelf (i.e., the clinothem topset) is storm-influenced mud where the slope comprises a mud wedge largely dominated by density current deposits (flood, low-density turbidites). Its architecture includes a composite stack of clinothem units (Fig. 7).

Sequence Stratigraphic Interpretation of the Clinothem Deposits

Individual Clinothem Units

The seismic (or clinothem) units, which correspond in cores to T-R depositional sequences, are composed of seismic sets, which correspond in cores to the following systems (facies) tracts: (1) seismic set 1 with amalgamated toe-of-slope debrites that drape the toe of the slope (lithofacies association FA 1.2); (2) seismic set 2 with low-angle progradational muds deposited at the clinothems topset (storm-influenced FA2.1b, FA2.2b); (3) seismic set 3 with flood- to storm-influenced offshore muds (FA2.1 and FA2.2) to shoreface silty sands (FA3) building the foreset where debrites and slope apron deposits (FA1.2) drape the toe of slope; and (4) seismic set 4, which is largely shoreface sand (FA3, FA4), pebbles, and gravels and gully fills (FA5) expressed as high-angle progradational reflectors at, or just beyond, the clinothem rollover (Figs. 5, 6, 7, and 8).

Set 3 shows a composite set of progradational reflectors with a marked increase in reflector angle and progradational rate with respect to seismic set 2. Internal unconformities of limited lateral extent produced by lobe switching during avulsion correlate to parasequence boundaries in coreholes. Seismic set 4 is bounded above by a maximum regressive surface capped by transgressive ravinement deposits on the shelf. This maximum regressive surface has been correlated in several instances with seismic sequence boundaries recognized by reflector terminations (Monteverde et al., 2008; Mountain et al., 2010; Miller et al., 2013a, 2013b).

According to Catuneanu’s (2006) sequence stratigraphy nomenclature, the suite of systems tracts could be interpreted in one of two ways: hypothesis 1 (forced regressive delta), in which a healing-phase transgressive systems tract (in set 1) is overlain by a composite package comprising a highstand systems tract (low-angle mud, set 2), a forced regressive wedge (high-angle, upper offshore to shoreface sand downlapping onto coarse slope apron, set 3), and a lowstand systems tract (coarse shoreface sands and gully fills, set 4); or alternatively, hypothesis 2 (highstand delta), in which a lowstand systems tract (toe-of-slope debrites in set 1) is overlain by a thick and composite highstand systems tract (offshore muds to shoreface sands in sets 2, 3, 4). In both interpretations, ravinement deposits at the top of the shelf and glauconitic condensed horizons at the toe of the slope correspond to the transgressive systems tract.

The first hypothesis (“forced regressive” delta) is supported by (1) the lateral change in reflector configurations from low-angle sigmoidal (progradational-aggradational) in highstand systems tract (set 2) due to the predominance of fine-grained lithologies (Patruno et al., 2015c) to sandier oblique (progradational) in forced regressive systems tract (set 3) and high-angle oblique in lowstand systems tract (set 4) building up the upper slope; (2) the lack of clear aggradational topsets at the clinothem unit scale; (3) the sharp angular surfaces at base of seismic sets 3 and 4; (4) the presence of a zone of separation and bypass (tens of kilometers large) between the shoreface sands at the rollover and subaerial delta at the shoreline; (5) the increasing average grain size in regressive deposits from proximal to distal settings, which implies an increasing winnowing of earlier highstand deposits; and (6) coarse-grained sedimentation on the slope apron at the base of set 3 (sequence boundary sensuPosamentier and Morris [2000]). The progressive increase in the angle of progradation together with the height of the foreset correspond to the lateral change from low angle subaqueous delta sedimentation during highstand time to the outbuilding of shelf-prism scale clinoforms during forced regressive and lowstand times.

The second hypothesis (“highstand” delta) is supported by (1) the overall progradational, i.e., normal regressional, seismic reflector configuration; (2) the lack of seaward downstepping subaerial erosion and incision at the rollover, indicative of a net sea-level fall; (3) the presence of offshore muds on the shelf instead of relics of shoreface-delta sands, which in the first hypothesis should bring the sand to the rollover; and (4) the lack of clear erosion surfaces in cores at the base of sets 3 (force regressive) and 4 (lowstand). In this “highstand hypothesis”, the dip of the clinoform foresets precludes the preservation of lowstand shoreface sands, which are reworked in mass flows downslope beyond the rollover in the slope apron.

The New Jersey clinoforms follow a progradational, horizontal trajectory, with clinoforms as high as the mid-shelf water depth, with no paralic or coastal plain deposits, and a shelfal foreset underlined by turbidite deposits, which correspond to the main basic characters of mid-shelf deltas (Porębski and Steel, 2006). In such a context, the normal regressional shape of the clinoforms in seismic data does not necessarily mean that the bulk of the clinoforms was deposited during highstand to stillstand times of relative sea level. Indeed, shelf-prism clinoforms can only accrete basinward and then can remain fixed through time (Helland-Hansen and Hampson, 2009) over periods as long as 104 to 106 yr (Vanney and Stanley, 1983). They bathe in a maximum water depth of about the same magnitude as the height of the clinoforms before starvation. Sigmoidal morphologies in set 2 (in hypothesis 2, “highstand delta”) are due to the predominance of fine-grained lithologies, the influence of basinal processes, the relatively slow progradation, and the high-angle trajectories (Porębski and Steel, 2006), as one can find in subaqueous deltas. With relative sea-level fall and the correlative increase in wave action on the shelf, sediments bypassing the shelf augment the sand content and sorting in the distal seaward-prograding increments (Sydow and Roberts, 1994; Anderson et al., 2004), producing basinward detached delta bodies generating the typical oblique profiles in sets 3 and 4 (in hypothesis 1, “forced regressive and lowstand deltas”), but no exposure of the shelf. The latter, at the difference of subaerial deltas, do not show important downcuttings during relative sea-level falls (Helland-Hansen and Hampson, 2009), and the unconformities can be quite remote (Porębski and Steel, 2006), but significant amounts of sand are delivered to the slope and the basin floor. After a major flooding, the shelf-prism clinoforms are separated from the highstand shoreline by a muddy shelf, but both can coincide when the shoreline delta transits across the shelf, giving rise to the so-called shelf-edge deltas (Porębski and Steel, 2003; Burgess et al., 2008; Helland-Hansen and Gjelberg, 2012). Transgressive as well as regressive wave ravinement can sharply truncate the shelf-prism clinoforms in wave-dominated coasts (15–40 m of erosion; Rodriguez et al., 2001; Bhattacharya and Willis, 2001; Cattaneo and Steel, 2002) as can geostrophic currents during early highstand (Rine et al., 1991).

The debate to determine whether the bulk of those shelf-prism clinoforms corresponds to highstand or forced regressive sediment bodies compares to the distinction made between supply-driven deltas (or “highstand deltas”) and accommodation-driven deltas (or “forced regressive deltas”) (Porębski and Steel, 2006). Supply-driven deltas (or “highstand deltas”) show thick, stacked, sandy parasequences sets and extensive muddy low-energy delta front deposits damped in increasing water depths. Relative sea-level falls produce no incision at the shelf edge, and shelf-attached sandy delta fronts rather than a basin-floor fan. Supply-driven deltas are very stable features where sequence boundaries (type 2 sensuPosamentier et al., 1988) can be unidentifiable at the base of shelf margin systems tracts. Accommodation-driven deltas (or “forced regressive deltas”) show horizontal progradational rather than progradational-aggradational trajectories, no coastal plain deposits, and a proximal-distal increase in wave action transporting sands at the rollover. They are more unstable features with internal unconformities. Large-amplitude sea-level falls produce shelf-edge detached delta bodies that steepen the delta front when those shelf-edge deltas reach the deeper mid-shelf. These sea-level falls do not expose the shelf because the sediments that bypass the shelf-enhance the sand content in the more distal portion of the seaward progradational elements (Sydow and Roberts, 1994; Anderson et al., 2004). Both supply- and accommodation-driven deltas are truncated in the landward direction by wave and geostrophic-current ravinement, which control the preservation of topset-foreset deposits, and are overlain by open-shelf mudstones. New Jersey clinothems share more characteristics with accommodation-driven (forced regressive) deltas, but some clinoform progradational-aggradational trajectories, low-energy muddy clinoform foresets, and the lack of incision at the shelf edge are indicative of the influence of sediment supply in the overall architecture.

Stacks of Clinothem Units

The clinothem units are stacked in four large T-R sequences (stacks A–D, Fig. 7). Each T-R sequence comprises three to five clinothem units showing a period of backstepping followed by a progressive but marked seaward migration of the offlap break. This evolution is interpreted as the result of a rise followed by a fall in relative sea level. The landward-stepping clinothem units are sand rich, enriched with glauconite and thinner than seaward-stepping units. During the rise in sea level, the rollovers build out and then backstep, accumulating clean sands in wave-dominated shoreface deposits during late lowstand and transgressive times. On the shelf, the clinothem units are truncated by wave ravinement surfaces overlain by storm-influenced clayey silts that drape the flat shelf. This surface, which marks the transition from low-angle to flat-topped shelf style, is formed by wave-induced current erosion possibly linked to the successive seaward and landward shifts of the shoreline, attested by the presence of relics of wave-dominated sharp-based shoreface sand deposits at the rollover. The seaward-stepping clinothem units are thicker and show an overall increase in the sand-mud ratio through time and thinning out by erosion at the top. Each of them comprises a suite of well-developed sediment packages including, from the oldest to the youngest: (1) low-angle, sigmoidal to oblique progradation of muddy clinothems of subaqueous deltas (highstand time), bypassed by (2) high-angle, shelf-edge deltas with switching lobes, high terrestrial inputs, and flood-induced hyperpycnal flows (forced regression time) and (3) disconnected delta front sands that accumulate at the shelf edge (lowstand time). In contrast with the landward-stepping units, the last two packages of sediments are enriched with terrestrial material like micas and organic detritus. The large T-R sequences show thin backstepping patterns that appear contradictory with Helland-Hansen and Hampson’s (2009) observations that shelf-scale clinoforms cannot backstep as the smaller-scale subaqueous deltas can do.

The New Jersey cross-section shows a progressive construction of the clinoform slope through time. Stack A is largely progradational. Individual sequences are strongly eroded, with just the lower part of their slopes preserved possibly due to higher-order forced regression. Stack B is progradational to aggradational, with a well-preserved sand-rich clinoform rollover. From this stack upward, the clinothems are taller (>250 m) and steeper (dip >2.5°), still migrating seaward drastically (>33 km laterally in stack D; Steckler et al., 1999). Large volumes of fine-grained material drape the clinoform front whose sediment differential compaction produced a depression at its toe moving in a seaward direction (rolling syncline effect of Steckler et al. [1999]).

A Depositional Model for the New Jersey Miocene Clinoforms

Core analysis and the reinterpretation of seismic line 529 show that a specific facies architecture, with mud-dominated deposits on the topset, foreset, and bottomset and sand-dominated facies at the rollover and the toeset, characterize the clinothems (Fig. 5). This spatial organization differs from that of predictive clinothem models (e.g., Catuneanu, 2006; Mountain et al., 2010b) and of previous studies based on seismic data alone (e.g., Poulsen et al., 1998; Monteverde et al., 2008). These models described a progressive increase in mud content in a seaward direction and a net accumulation of sand in the clinothem at the depositional shelf break.

ODP coreholes drilled on the New Jersey coastal plain record a stack of ∼30-m-thick depositional sequences (Browning et al., 2008). The sequences display a basal muddy sandstone layer overlain by prodelta offshore marine muds giving rise upsection to delta front wave-dominated shoreface sands (Browning et al., 2008; Fig. 7). These successions are interpreted as amalgamated transgressive and highstand deposits of Miocene depositional sequences coeval with the clinothems drilled by IODP Expedition 313 coreholes 40–60 km in the offshore (Browning et al., 2008, 2013; Kominz et al., 2016). The falling stage and the lowstand deposits interpreted in the offshore are both missing in the onshore sequences.

Sediment Transport

The presence of mud-rich facies during highstand times on the topset of the clinoform on the shelf, between shoreface sand deposits at the rollover and sand-prone shoreline delta deposits 40–60 km landward, implies two possibilities: either (1) an important and rapid erosion that removed the entire uppermost sand-prone material on the topset of the clinoform after a 60+ km shoreface progradation and a shift of all of the facies belts to the depositional shelf break (rollover) (e.g., Porębski and Steel, 2006), or (2) a compound depositional system that involves a subaerial delta that fed a distant subaqueous delta (e.g., Helland-Hansen and Hampson, 2009; Patruno et al., 2015) by current transportation. The first hypothesis requires the transit of deltas and shoreline bodies across the shelf during forced regression and lowstand times, providing sands to the rollover (Posamentier and Morris, 2000), with deltas being subsequently erased by transgressive ravinement (Porębski and Steel, 2006). The shelf is flattened twice by wave action during forced regression and transgression-early highstand. The second hypothesis requires significant near-bed shear stresses in shallow marine areas marked by high-energy waves, tides, or currents causing shelf flattening and bypass through lateral advection and resuspension of sediment (Pirmez et al., 1998; Driscoll and Karner, 1999; Swenson et al., 2005; Cattaneo et al., 2007; Patruno et al., 2015a; Poyatos-More et al., 2016). Deposition occurs preferentially seaward of the rollover where near-bed agitation decreases below the threshold of sediment motion (Mitchell, 2012; Mitchell et al., 2012). The latter hypothesis is well documented for subaqueous deltas, which are the prodelta counterparts of subaerial deltas. But, these compound bodies are small, at delta scale, a few meters to tens of meters high, and can step back and forth during transgression-regression processes, which is different from what we observe on the very stable New Jersey shelf-prism clinoforms. Although many aspects of the dynamic and geometries of the shelf-prism and delta-scale clinoforms are reputedly scale invariant (Patruno et al., 2015a), no clear sediment transport process is documented for sand transportation across the shelf by suspension and advection at the large-scale shelf-prism clinoforms. Moreover, the accommodation-driven cross-shelf transportation of sand postulated in the first hypothesis is attested by the relics of shoreface sands in the transgressive lags, the flatness of the shelf enhanced by the back-and-forth passage of the “wave base razor”, a highstand water depth over 60 m beyond the fair-weather and storm wave base, and the presence of sharp surfaces at the base of the interpreted forced regressive and lowstand systems tracts indicative of drastic downward shifts of sediment. Nevertheless, it is very likely that both processes interact. A net fall in sea level, with a seaward migration of the shoreline delta and its subaqueous counterpart, would enhance near-bed shear stresses on the shelf, transporting sand to the rollover, without the strict necessity for the subaerial delta to prograde to the shelf edge. This mechanism would easily explain the lack of subaerial erosion (or deposition) in the corehole at the rollover and on the topset of clinoforms.

The Abundance of Mud

The New Jersey clinothems are made up of 70% muddy (clayey silt to silty clay) sediment mainly concentrated in highstand systems tracts, with sands being mostly concentrated at the rollover and at the toe of the slope. Shelfal muds result from the dispersal of river supply by dip-oriented current in high-sediment-supply conditions, or from alongshore currents settling contourite drifts in shallow waters (Hanebuth et al., 2015; Patruno et al., 2015a, 2015b, 2015c). Depending on sediment supply and tide and wave energy, mud dispersal on the shelf can form muddy clinothems detached from their subaerial counterparts (“subaqueous delta clinoforms”, e.g., Fly River [Papua New Guinea] or Amazon River, in Walsh and Nittrouer, 2009) or shelf mud wedges (“marine dispersal dominated”, e.g., Eel River [northern California, USA], in Walsh and Nittrouer, 2009). In the Holocene, shelf mud wedges seems mostly recorded in enclosed seas (e.g., Adriatic Sea: Trincardi, et al., 1994; Cattaneo et al., 2003; Yellow Sea and East China Sea: Liu et al., 2004, 2006), while deep water muddy clinothems are reported in open marine conditions, such as in the western and eastern Atlantic Ocean (Nittrouer et al., 1996; Hanebuth and Lantzsch, 2008; Hanebuth et al., 2015). Shelf mud wedges require a bathymetric low on the shelf to accumulate. An inner-shelf 6–12-m-scale low, deduced from backstripping, is thought to be the result of a dynamic topography effect (Kominz et al., 2016), but no wedging out of the seismic reflectors in a seaward direction on seismic sections substantiates this calculation, although it may be too small to be discriminated. The New Jersey clinoforms develop on a gently dipping shelf with a ramp-like morphology (apparent dip of 0.75°–0.5°; Miller et al., 2013a) bathed below mean storm wave base, in 30–50 m of water depth (e.g., Katz et al., 2013), 40–60 km seaward of the shoreline. This distance from the shoreline seems very likely because subaqueous delta clinoforms show higher accumulation rates than their subaerial counterparts and prograde three times further and faster than their contemporaneous shoreline (Patruno et al., 2015a). Such a separation is observed in the Holocene subaqueous deltas of the Amazon, with the rollover in ∼80 m water depth, 10–100 km from the shoreline (Nittrouer et al., 1996; Walsh and Nittrouer, 2009), and in the ancient Blackhawk subaqueous delta (Santonian–Campanian [Utah]) sited in 50–80 m of water depth, 25–70 km from the paleocoastline (e.g., Hampson, 2010).

A Composite Scenario of Deposition

A composite scenario, comprising the interbedding of subaerial and subaqueous deltas and shelf-edge deltas, probably best represents the depositional model of the New Jersey clinoforms (Fig. 9). The subaerial delta and its muddy flood-influenced subaqueous delta prograde at low angle during normal regression of highstand time. During a net sea-level fall, the subaqueous delta is progressively overstepped by the subaerial delta. The increase of wave and storm currents and near-bed stresses on the shelf brings sand to the rollover, building shelf edge deltas (Dixon et al., 2012) in forced regressive and lowstand systems tracts of shelf-prism clinothems. During the subsequent rise in sea level, the rollover builds out and then backsteps, accumulating clean sands in wave-dominated shoreface deposits. The overlying transgressive deposits confirm that the shelf-edge delta deposits represent the lowstand system tract (shelf margin wedge type, above a type 2 sequence boundary sensuPosamentier et al. [1988]) not preserved in ODP onshore coreholes (Browning et al., 2008). On the shelf, the subaqueous deltas are truncated by wave ravinement surfaces overlain by storm-influenced clayey silts that drape a flattened shelf.

Controls on Clinoform Architecture

Clinoforms are the basic building blocks of passive margin sedimentary architectures (e.g., Walsh et al., 2004; Helland-Hansen and Hampson, 2009; Patruno et al., 2015a). Their architecture is controlled by a number of forcing parameters that includes autocyclic (e.g., waves, currents, topography) and allocyclic (sea level, subsidence, sediment supply, and climate) controls.

The Clinoforms Record Changes in Accommodation Space

The New Jersey clinoforms as a whole show an overall concave-up ascending trajectory typical of a long-term relative sea-level rise compensated by a steady sediment supply filling space made available by eustasy and subsidence. The heights of the clinoforms increased progressively (26 m in sequence m6, to 46 m in sequence m5.6, 116 m in sequence m5.4, 131 m in m5, to a maximum of 250 m in sequence m4; Steckler et al., 1999) with a quite constant flooding water depth of 10 m to a few tens of meters. In detail, however, this architecture is composed of four stacks of clinothems, which are the consequences of marked differences in the response of the sedimentary system to changes in relative sea level (Figs. 7 and 9). Stack A is made up of low-gradient, sigmoidal, progradational, early Miocene clinoforms built by the wandering of subaqueous and shelf-edge delta lobes sharply truncated by an erosion surface during a net fall in relative sea level. Stack B is composed of oblique, progradational-aggradational, early middle Miocene clinoforms with well-preserved rollovers and clinoform slopes built during an overall rise in relative sea level. The clinothems comprise flat-top, low-gradient subaqueous delta lobes bypassed by wave-dominated, shoreface delta sands during forced regression due to the reduced accommodation space on the shelf (Fig. 9). The wave-dominated, shoreface delta sands accumulated with high-angle, oblique geometries, where space is available for sediment to accumulate, i.e., beyond the rollover. The combined action of (1) storm waves during highstand and fair-weather waves during forced regression and (2) the seaward shift of the shoreface toward the rollover (so-called regressive surface of marine erosion; e.g., Proust et al., 2001) shaped the flat surface at the top of the subaqueous delta. The shoreface sand and the subaqueous delta may be separated by tens of kilometers laterally and tens of meters vertically in water depth, depending on the intensity of wave and tide activity like that observed in the Ganges, Brahmaputra, Amazon, and Huang He (Yellow) river examples (Nittrouer et al. 1996; Kuehl et al., 1997; Michels et al., 1998; Hernández-Molina et al., 2000; Walsh et al., 2004; Patruno et al., 2015a, 2015c). From stack C upward (Fig. 7), the clinothems became definitively taller (>250 m) and steeper (>2.5°) and started migrating seaward at a faster pace. Large volumes of fine-grained sediment draped the slope during highstand times.

The Importance of a Flat Shelf (Clinoform Topset)

The progressive flattening of the New Jersey shelf during the Miocene increased the sedimentation rate along the clinoform foreset (cf. stacks B and C) feeding slope-apron and toe-of-slope fans (cf. stack A). This increased (1) the progradation rate of the shelf (3–5 km/m.y. in the early Miocene to 16–17 km/m.y. in the middle Miocene; Steckler et al., 1999; Carvajal et al., 2009), (2) the height of the slope (22–46 m in the early Miocene to 175–256 in the middle Miocene) and (3) the dip of the slope (1°–2.7°) of the clinoforms with almost no change in the highstand shelf water depth (30–75 m; Katz et al., 2013). Surprisingly, the New Jersey shelf-prism clinoforms do not show any indication of significant autoretreat (Helland Hansen et al., 2012) due to overextension and insufficient sediment supply, nor oversteepening and increasing water depth during the early Miocene to early mid-Miocene. These long-term rates of progradation are similar to those retained globally (10–500 km/m.y.; Patruno et al., 2015a), but the rate of progradation seems to have slowed down a little bit later (m4, 12.5 Ma; Steckler et al., 1999) probably due to oversteepening when the dip of the clinoforms reached 2.5°. With such a dip, 100 m of drop in sea level, which would have been unlikely, would have only produced 3 km of progradation (Burgess et al., 2008).

The Importance of Climate Change

The increase in the intensity of waves (height, recurrence intervals) favors the separation between subaqueous and subaerial deltas (Kuehl et al., 1997; Swenson et al., 2005), and as a consequence, the formation of a flat topset, the decrease in flood frequency and fluvial discharge, the overall progressive decrease in sediment grain size (from clinothem m5.45 onwards), as well as the increase in sedimentation rate on the slope of the subaqueous delta clinoforms. All of these are recognized as preliminary signals that might indicate the entry into the Neogene icehouse world (Helland Hansen et al., 2012). Large continental ice sheets grew and decayed in Antarctica in the early Miocene to early middle Miocene (Barrett et al., 1987; Miller et al., 1991; Zachos et al., 2001) (reduced ice sheets in the Northern Hemisphere at that time; Wolf-Welling et al., 1996; Wright and Miller, 1996) but stabilized with the installation of a permanent large continental ice sheet in East Antarctica at ca. 14.6–13.0 Ma (deep-sea Miocene isotope events Mi2a, Mi3, Mi4; Miller et al., 1991) during deposition of clinothems m5.2 to m4.1 (Browning et al., 2013). The later m4.1 corresponds to the onset of large-amplitude and low-frequency sea-level changes that resulted in the increasing rate of clinoform progradation.

CONCLUSION

Passive margins are considered as the places of simple, undeformed sedimentation that can be easily correlated over large distances, which theoretically offers the possibility of quantifying changes in global sea level. The classical depositional model is composed of prograding clinoforms well imaged at a large scale by seismic techniques. Their internal structure is rarely exposed or cored, although detailed knowledge of their architecture can definitively improve our understanding of their origin and the distribution pattern of geological reservoirs.

The New Jersey clinothem complex comprises mud with ∼30% sand and coarse material concentrated at the rollover, the foreset, and the toe of the clinoform and beyond. Clinothems nucleate by delta progradation on topography three to five times steeper than the present-day shelf. Their compound depositional system involves a subaerial delta feeding a shelf-edge subaqueous delta lying on a gently deepening ramp-like morphology, typically set below fair-weather wave base. The low-angle foreset muds prograde during highstand times when subaerial shoreface sand and coarse material transit to the rollover and the toe of the clinoform during forced regression and lowstand times, where it accumulates in large shelf prism–scale clinothems. Sea level has never dropped below the rollover, but the back-and-forth displacements of the wave-base razor as well as the geostrophic current flatten the shelf and control the decoupling of the subaerial and subaqueous shelf-edge deltas.

The New Jersey clinothems show a curious example of the dynamic equilibrium between the geometry of the sediment body and control parameters (eustasy, subsidence, and sediment supply): the height of the clinoforms, the dip of the slope, and the progradation rate increase with no change in the water depth on the shelf, strictly controlled by the wave-base razor. The latter appears to be a first-order parameter controlling sediment progradation responsible for sediment transportation and bypass from the subaerial delta shoreline to the toe of the clinoforms, and, depending on sea-level and climate changes—climate cooling intensifies the wave regime—controls the overall geometry of the clinothems.

ACKNOWLEDGMENTS

This paper used samples and data provided by the Integrated Ocean Drilling Program (IODP) and the International Continental Scientific Drilling Program (ICDP). The research was supported by the CNRS-Institut National des Sciences de l’Univers 2010 Post-Campagne Research Program and a CNRS-IODP France postdoctoral grant A64622 to Hugo Pouderoux. We acknowledge funding from NERC (NE/F001428/1) to Hesselbo, and NERC (NE/H014306/1) to Hodgson. We are grateful to two anonymous reviewers for the careful review and editing of a previous version of this manuscript. This work has benefited from their valuable expertise and ideas. This paper is part of Geosphere themed issue Results of IODP Exp313: The History and Impact of Sea-Level Change Offshore New Jersey with guest associate editors G. Mountain, J.-N. Proust, and D.B. McInroy.

1Supplemental Files. Figure S1 is detailed lithological sections of hole M0027 drawn from the visual core descriptions. Figure S2 is detailed lithological sections of hole M0028 drawn from the visual core descriptions. Figure S3 is detailed lithological sections of hole M0029 drawn from the visual core descriptions. Please visit https://doi.org/10.1130/GES01545.S1 or the full-text article on www.gsapubs.org to view the Supplemental Files.
Science Editor: Shanaka de Silva
Guest Associate Editor: David Burns McInroy
Gold Open Access: This paper is published under the terms of the CC-BY-NC license.