Morphometric analysis of submarine fan systems, the largest sedimentary deposits on Earth, demonstrates scaling relationships between genetically related channels and lobe-shaped bodies (LBs) deposited beyond the channel terminus, providing insight into the architectural development of these systems. Compiling dimensional data from depositional systems that cover a range of sediment supply characteristics, tectonic settings, and geographic locations enables investigation into global trends in depositional morphology. LBs have a consistent, scale-independent length-to-width ratio of ∼2:1. The thickness-to-area ratios for LBs show multiple morphologic trends, likely driven by topographic confinement, with LBs getting proportionally thicker in relation to increasing confinement. Morphometric analysis of genetically related channel dimensions (width, relief, cross-section area) and LB dimensions (length, width, thickness, area, volume) reveals robust scaling relationships; most notably, channel width and cross-sectional area can be used to predict the volume and depositional area of related LBs. These relationships demonstrate that LBs proportionally scale to their concomitant channels, and thus to the volume of sediment supplied prior to an avulsion. While the dimensions of submarine fans scale to associated terrestrial catchments, the building blocks of submarine fans (i.e., channels and LBs) do not, suggesting a down-system decoupling (or lack of scaling) at LB deposition time scales. Applying these morphometric trends and scaling relationships as input parameters for source-to-sink and reservoir models can improve predictions of stratigraphic architecture, sediment partitioning, and sediment/carbon flux in modern and ancient submarine fan systems.
Scaling relationships of sediment-routing systems provide insights into the intrinsic properties and processes of those systems by identifying linkages between morphometric and other quantitative characteristics (e.g., hydraulic properties and sediment flux). For example, scaling relationships have been identified between segments within source-to-sink systems (e.g., Sømme et al., 2009), as well as in finer-scale studies within specific depositional environments (rivers, estuaries, deltas; e.g., Leopold and Langbein, 1962). However, relatively few studies have focused on scaling relationships in submarine fans, which are net-depositional environments within continental-margin sediment-routing systems that are diverse in size and morphology (Fig. DR1 in the GSA Data Repository1). These composite depositional features consist of channels and multiple scales of lobe-shaped bodies (LBs), formed by channelized and unconfined sediment gravity flows, respectively (Piper and Normark, 1983). While previous studies have focused on the scaling relationships of submarine channels (e.g., Covault et al., 2012; Reimchen et al., 2016) or LBs (e.g., Jegou et al., 2008; Prélat et al., 2009, 2010), scaling between submarine channels and LBs remains unexplored. In this study of Eocene to modern fan systems, we (1) document the morphometrics of, and scaling relationships between, submarine channels and concomitant (i.e., genetically related) LBs from element to complex scale (i.e., larger than bed scale); (2) evaluate the utility of these relationships for predicting depositional body dimensions in data-poor areas; and (3) investigate potential relationships between terrestrial portions of sediment-routing systems and channel-LB morphology.
DATA AND METHODS
We compiled dimensional data for LBs (n = 271) and, when possible, concomitant channels (n = 52, 160 LBs have paired channel data) from 35 submarine-fan systems (Fig. 1A). The data set of channel and LB measurements made from sonar and seismic-reflection surveys (Fig. 1B) includes submarine-fan systems that cover a range of geographic locations, tectonic settings, source-to-sink configurations, and sediment-supply characteristics (e.g., catchment size; see Table DR1 in the Data Repository), which enables investigation into how system characteristics influence depositional morphology and scaling. Dimensional data from basin floor, slope, and ponded fan settings are included, but we excluded LBs showing evidence of significant sediment bypass (i.e., thoroughgoing channels). Confinement (sensuPrélat et al., 2010) at the system scale is documented as ‘confined’ when basin topography limits the runout of the LBs, or ‘unconfined’ when LB runout is not significantly impacted by topographic variability or basin margins. We acknowledge that the morphology of individual channel-LB systems changes as they grow (Deptuck et al., 2008); this variability is endemic to our data set, which provides the most comprehensive collection of channel-LB dimensions published to date.
These equations do not account for the range of LB morphologies, which causes a mean absolute percent error of calculated versus measured ALB and VLB of ∼25%. Compensational stacking occurs at multiple scales within LBs, and various hierarchical schemes are used in the literature to describe the resultant deposits (e.g., Deptuck et al., 2008; Prélat et al., 2009; Straub and Pyles, 2012). For consistency, we applied standardized hierarchical terminology to all LBs by qualitatively matching their original documentation to the hierarchy defined by Prélat et al. (2009). This hierarchy consists of four levels (Fig. 1C) in increasing size and complexity: ‘beds’ (H1) deposited by an individual event/turbidity current; ‘lobe elements’ (H2) composed of stacked beds/bed sets; ‘lobes’ (H3) formed by one or more stacked lobe elements fed by a single channel; and ‘lobe complexes’ (H4) that develop when avulsions or channel migrations result in development of multiple lobes (Prélat et al., 2009). Individual beds (H1) are not included in our analysis due to limited data availability.
Morphometrics of Submarine Channels
Submarine channel aspect ratios (WCH/HCH) in this study (median = 44) are consistent with those documented by Konsoer et al. (2013; Fig. 2A; Table DR1). Channel dimensions from this study are generally small compared to those of Konsoer et al., likely because submarine channels typically decrease in size downstream and, unlike Konsoer et al., we exclusively measured distal channel reaches (Fig. 2A).
Morphometrics of Lobe Bodies
A robust power-law relationship exists for the planform morphology (LLB and WLB) for all LBs (LLB = 1.74 WLB1.02, r2 = 0.86, n = 271) and similar trends exist within each hierarchical level (r2 values >0.72; Fig. 2B; Table DR1). The planform aspect ratios (LLB/WLB) of LBs demonstrate a relatively narrow range of planform geometries (circular versus elongate), with a median of 2:1 (Fig. 2B; P10 = 1.2, P50 = 2.0, P90 = 4.4). Planform dimensions generally increase from H2 to H4 within individual systems, but considerable overlap between hierarchical levels indicates that their dimensions are not globally consistent (Fig. 2B), likely due to differences in the formative properties (e.g., flow volumes and grain size). There is more than an order-of-magnitude difference between the average volumes of each hierarchical level of LB (H2 = 0.07 km3, H3 = 2.9 km3, H4 = 62.0 km3), similar to the findings of Prélat et al. (2010). However, significant overlap between the dimensions of hierarchical levels indicates that hierarchy cannot be inferred solely from LB dimensions (cf. Prélat et al., 2009).
In contrast to their planform morphology, the three-dimensional (3-D) morphologies of LBs do not show simple scaling. Figure 2C documents the significant variability in the thickness-to-area ratio (HLB/ALB), revealing two trends. Trend 1 illustrates the distribution of confined LBs, which are proportionally thicker than unconfined LBs illustrated in Trend 2. These trends agree with observations of ‘lobe’ (H3) morphologies made by Prélat et al. (2010), but our data incorporate more hierarchical levels and a wider variety and number of systems. The two trends are distinct at lower ALB values, but converge as ALB increases; accordingly, the trends are identified in H2 and H3, but are indistinguishable for H4. Trend 1 shows a shift to lower HLB/ALB (m/km2) ratios with increasing area and hierarchy level (from ∼10:1 to <1:1), while Trend 2 typically remains between 1:10 and 1:100 (Fig. 2C).
Scaling Relationships between Channels and Lobe Bodies
Our morphometric analysis demonstrates several strong scaling relationships between concomitant channels and LBs. Statistical analysis of channel (WCH, HCH, ACH) and LB (ALB, HLB, VLB) dimensions for hierarchical levels H2–H4 shows that WCH has a strong, positive power-law scaling with both ALB and VLB (r2 ≈ 0.7; Fig. 3). ACH also shows a positive, albeit weaker, power-law scaling with both ALB and VLB (r2 ≈ 0.6; Fig. 3). HCH does not correlate with any LB dimensions (r2 < 0.27), and no channel dimension correlates with HLB (r2 < 0.26; Fig. 3). Additionally, our investigation of possible relationships between channel-LB dimensions and sediment-routing–system parameters (fluvial catchment area, water discharge, and suspended sediment load) found no significant correlations (r2 < 0.35; Fig. 4).
Controls on Morphology and Volume of Lobe Bodies
Planform aspect ratios of LBs display a high degree of consistency across all hierarchical levels (Fig. 2B). Having a strong LLB-WLB relationship in spite of the wide range of system characteristics represented by the data set highlights a potentially universal trend in LB planform geometry, and indicates that turbidity flow properties (rather than external factors) impart primary control on LB planform geometries. This is supported by flume experiments, which show that flow properties (density, volume, sediment type) primarily impact deposit geometry (Baas et al., 2004), while factors that modify the flow processes (e.g., topographic confinement) seem to be secondary (Al Ja’Aidi et al., 2004). The dimensions of LBs change throughout the time span of their deposition, which is dictated by avulsions (Deptuck et al., 2008); accordingly, our measurements of modern, active LBs will underrepresent their final deposit dimensions.
The presence of multiple trends between HLB and ALB (Fig. 2C) indicates that the relative thickness of LBs is more sensitive to external factors. If the impacts of external factors (e.g., confinement) on LB geometry are known, it is possible to predict HLB-ALB relationships using the trends in Figure 2C. While flow properties such as sediment grain size and concentration impact sediment dispersal, and thus LB thickness (Baas et al., 2004), our data demonstrate that topographic confinement strongly influences the 3-D morphology of LBs (Fig. 2C; also see Covault and Romans 2009; Prélat et al., 2010). When unconfined, flows spread out to produce thin LBs with relatively large depositional areas, but when confined, the depositional areas of flows are limited, resulting in LBs with greater relative thickness (Al Ja’Aidi et al., 2004).
Additionally, the two trends in Figure 2C have different implications for the dependence of compensational stacking on VLB and/or the hierarchical level of LBs. Straub and Pyles (2012) suggest that compensational stacking should cause reduced HLB/ALB ratios with increased hierarchical level. Confined systems (Trend 1) in Figure 2C align with Straub and Pyles’s (2012) prediction; however, unconfined systems (Trend 2) have a relatively consistent HLB/ALB ratio regardless of LB volume or hierarchical level. If confinement is indeed driving the trends in Figure 2C, the impact of compensational stacking on deposit morphology is variable across LB hierarchies in confined systems, but consistent in unconfined systems.
Controls on Channel-Lobe Body Scaling
WCH and ACH have strong scaling relationships with ALB and VLB (Fig. 3). The lack of correlations between the vertical measurements (HCH, HLB) and other parameters of channel-LB systems may stem from basin confinement (as discussed above) as well as the lower relative accuracy of vertical (thickness, relief) measurements (Fig. 3). The data in this study support a simple association between turbidite flow properties, channel morphology, and lobe dimensions. Because the hydraulic geometries of submarine channels are proportional to their formative sediment-gravity-flow events (Konsoer et al., 2013) and VLB is proportional to the dimensions of the associated channel (Fig. 3), accordingly, VLB is proportional to the properties of the LB-forming sediment gravity flows (e.g., volume and number of flows).
Submarine-fan dimensions scale with associated catchment parameters (Sømme et al., 2009), but this study indicates that the dimensions of channels and LBs that comprise fans do not (Fig. 4). We did not find any significant correlations between terrestrial catchment parameters (area, water discharge, and suspended sediment load) and submarine channel or LB dimensions (Fig. 4). While these relationships exist when comparing catchments to overall fan dimensions (Sømme et al., 2009), our data suggest that, at smaller spatiotemporal scales (i.e., LB deposition over 102–104 yr; Deptuck et al., 2008; Jegou et al., 2008), the terrestrial and submarine portions of these systems are decoupled due to autogenic processes (e.g., avulsion) and incomplete transfer of sediment (Romans et al., 2016).
Predicting Channel and Lobe Body Dimensions
Scaling relationships between submarine channels and LBs (Figs. 3A, 3B, 3D, and 3E) provide a tool for predicting the dimensions of these sedimentary features, if the scale of one parameter is known (e.g., from subsurface seismic data). Figure 5 contains predicted ranges of ALB and VLB for channels with widths of 0.5 km and 1.0 km, using Figures 3A and 3B. Data limitations do not allow direct prediction of hierarchical levels of LB dimensions; however, because LB dimensions generally increase with hierarchy, H2 through H4 LBs should successively fall between P10 to P90 for a given channel dimension in Figure 3. This information is valuable because hierarchical level influences the architectural complexity of LBs (Pyrcz et al., 2005). Channel-LB scaling provides insight into the sedimentation dynamics that generate the depositional architecture of submarine fans; however, when applying these relationships to ancient systems, differences between active channel conduits and channel-fill deposits must be addressed. The lack of correlations involving HCH or HLB in Figure 3 may stem from their natural variability due to knickpoint migration in channels and antecedent topography for LBs, or the relatively low precision of vertical versus lateral measurements.
This study documents morphometric trends for lobe bodies in submarine fans, scaling relationships between concomitant submarine channels and LBs, and decoupling of sedimentary processes and fluxes between key parameters of submarine channel-LB deposits and associated terrestrial drainages. LBs have consistent planform length-width relationships (LLB = 1.74 WLB1.02; r2 = 0.86, n = 271) and a restricted range of planform aspect ratios (length/width): P10 = 1.2, P50 = 2.0, P90 = 4.4. There are two trends in the 3-D morphology of LBs: (1) LBs that form in topographically confined settings tend to have higher maximum thickness/area ratios that decrease with increased scale, and (2) LBs that form in topographically unconfined settings tend to have lower and hierarchy-independent maximum thickness/area ratios. While distinct for lower hierarchical levels, these trends converge and become indistinguishable at higher hierarchical levels. We interpret that topographic confinement is the primary driver of the thickness/area trends of LBs, but additional factors (e.g., grain size) may also influence these trends.
Submarine channel width and cross-sectional area show power-law scaling with both the depositional area and volume of concomitant LBs. The relationship between channel width and LB volume (VLB = 10–10.78·WCH3.95; r2 = 0.77) indicates that (1) channel dimensions scale to the total volume of sediment that passes through the channel before an avulsion or shift in the location of deposition causes a new LB to be established, and (2) LBs scale to the volumes of turbidite flows passing through channels. Applications of these dimensional ranges and scaling relationships include providing realistic input parameters for source-to-sink modeling of sediment budgets, and reservoir models for natural-resource characterization. The lack of correlations between terrestrial catchment-scale parameters and channel-LB dimensions suggests decoupling of sediment fluxes/processes between the terrestrial and deep-water portions of sediment routing systems.
Financial support for this work was provided by Chevron Corporation to the Chevron Center of Research Excellence at the Colorado School of Mines (USA). This paper benefited from discussions with Ashley Harris, Fabien Laugier, Jeremiah Moody, and Morgan Sullivan. This manuscript benefited substantially from reviews by Ian Kane, Amandine Prélat, Tor Sømme, an anonymous reviewer, and editor James Schmitt.