Confinement width and inflow-to-sediment discharge ratio control the morphology and braiding intensity of submarine channels: Insights from physical experiments and reduced-complexity models

. Submarine channels conveying sediment gravity flows are often topographically confined, but the effect of confinement width on channel morphodynamics is incompletely understood. We use physical experiments and a reduced-complexity model to investigate the effects of confinement width ( 𝐵𝐵 ) and inflow-to-sediment discharge ratio ( 𝑄𝑄 𝑖𝑖𝑖𝑖 / 𝑄𝑄 𝑠𝑠 ) on the evolution of submarine braided channels. The results show that a larger confinement width results in increased active braiding 15 intensity ( 𝐵𝐵𝐵𝐵 𝐴𝐴 ), and that 𝐵𝐵𝐵𝐵 𝐴𝐴 takes longer to stabilize (i.e., a longer critical time, 𝑡𝑡 𝑐𝑐 ). At a fixed confinement width, a higher 𝑄𝑄 𝑖𝑖𝑖𝑖 / 𝑄𝑄 𝑠𝑠 slightly decreases the 𝐵𝐵𝐵𝐵 𝐴𝐴 . Digital Elevation Models of Difference (DoD) of the experiments allow measurement of the morphological active width ( 𝑊𝑊 𝑎𝑎 ) of the submarine channels and the bulk morphological change ( 𝑉𝑉 𝑏𝑏𝑏𝑏𝑏𝑏𝑏𝑏 ) within an experiment, defined as the sum of total erosion and deposition. We find that 𝑊𝑊 𝑎𝑎 and 𝑉𝑉 𝑏𝑏𝑏𝑏𝑏𝑏𝑏𝑏 are proportional to 𝐵𝐵 . We further confirm that 𝐵𝐵𝐵𝐵 𝐴𝐴 is proportional to both dimensionless sediment-stream power ( 𝜔𝜔 ∗∗ ) and dimensionless stream power ( 𝜔𝜔 ∗ ). These trends 20 are consistent for submarine braided channels both with and without confinement width effects. Furthermore, we built a reduced-complexity model (RCM) that can simulate flow bifurcation and confluence of submarine braided channels. The simulated flow distribution provides reliable predictions of flow depth and

For instance, experiments of laterally unconfined submarine braided channels were first confirmed at two laboratory facilities (Foreman et al., 2015), and demonstrated that density currents readily produce submarine braided channels for flow aspect ratios (depth-to-width) similar to those of fluvial braided rivers.The authors confirmed that the stability theory for river planform morphology (Parker, 1976) successfully describes submarine braided channels in both experiments and the field.Subsequent experiments (Lai et al., 2017) further demonstrated that laterally unconfined submarine braided channels (Fig. 1a) have similar responses to fluvial braided rivers (Egozi and Ashmore, 2009), for which both the active braiding intensity (  ) and total braiding intensity (  ) increase with increasing inflow discharge and bed slope.  is defined as the total number of channels in cross-section, while   is a subset of cross-sectional channel numbers, which reflects channels that can transport bed load, influence channel morphology and relate to channel bifurcations and avulsions (Ashmore, 2009;Egozi and Ashmore, 2009).Past fluvial and submarine braiding experiments (Ashmore 1991;Egozi and Ashmore, 2009;Bertoldi et al. 2009;Lai et al., 2017) show that   scales linearly with dimensionless stream power ( * ).While the recently proposed entropic braiding index (Tejedor et al., 2022) accounts for variation in the width of individual channel threads, in this study we use the traditional braiding index for ease of comparison to these previous studies.
In natural rivers, lateral confinement reduces braiding intensity (Garcia Lugo et al., 2015;Carbonari et al., 2020).In submarine settings, canyon walls, erodible terraces, and basement ridges could also limit the development of submarine braided channels.
For instance, Orinoco deep-sea fan is the first reported field evidence of submarine braided channels and have a width-to-depth ratio around 60 to 70 (Belderson et al., 1984;Ercilla et al., 1998).Further downstream of these submarine braided channels, confinement by the Barbados deformation front and a basement ridge force the channels to transform into a single channel (Belderson et al., 1984;Callec et al., 2010).Along the Sicilian margin, submarine braided channels occur on the Stromboli slope valley (south eastern Tyrrhenian Sea) that is confined by a terrace and valley wall (Gamberi et al., 2011).In the same region, the Baia di Levante Fan (Vulcano Island, Italy), a modern submarine volcaniclastic fan, was formed by repetitive stacking of gravity flows (Romagnoli et al., 2012), with clear erosive banks confining the observed submarine braided channels (Fig. 1b).and (4) to establish scaling relationships across fluvial and submarine braided channel systems.

Experiments
The experimental basin (1.8 m long, 0.55 m wide, and 0.05 m deep) was submerged into a water tank (2.3 m long, 1 m wide, and 0.65 m deep) to investigate the effect of different confinement widths on the evolution of submarine braided channels (Fig. 2).The water tank provided still, ambient water to simulate a deep-sea environment.Three drains were set at the downstream end of the tank to remove turbid water.The initial bed slope ( 0 ) of the basin is adjustable.Prior to experiments, plastic sand (specific gravity = 1.5, median grain size  50 = 0.34 mm, uniformity coefficient  =  60 / 10 = 1.64) was formed into a planar bed within the basin.Then a straight channel was excavated in the middle with a fixed depth of 0.015 m and a prescribed confinement width.At the upstream end, the same dry plastic sand was provided steadily by a motor-controlled conveyor belt to simulate bed load sediment input.Dyed and saturated saline inflows (  = 1200 kg/m 3 ) were supplied upstream to simulate long-lived, unconfined hyperpycnal flows or mud-rich turbidity currents (Métivier et al., 2005;Spinewine et al., 2009;Sequeiros et al., 2010;Weill et al., 2014;Foreman et al., 2015;Lai et al., 2016;2017).Check-dam partitions at the inflow exits allowed the saline discharges and dry plastic sands to be fully mixed before entering to the basin.
Six runs of experiments were conducted to investigate the effects of different confinement width () and inflow-to-sediment discharge ratio (    ⁄ ) on the evolution of the submarine braided channels.Combinations of confinement  (at 0.12 m, 0.24 m, and 0.48 m) and     ⁄ (at ~ 60 and ~90) encompass Series A and Series B (Table 1).The confinement widths were designed to distinguish from previous experiments (Foreman et al., 2015 andLai et al., 2017)  = 60 and 90 are the two reasonable values of this ratio.The initial bed slope was kept at 0.114 (6.5 degrees, the most effective bed slope for forming submarine braided channels reported by Lai et al., 2017).Each experimental run lasted 100 min and was divided into 10 successive stages (10 min for each stage).During the experiments, time-lapse images were obtained using a Nikon D7100 at intervals of 5 sec.At the end of each stage, the inflow and sediment supply were stopped, and topography temporarily stabilized.Without draining the ambient water, the terrain surface was scanned.After that, a high-resolution (1 mm x 1 mm) digital elevation model (DEM), orthorectified image, and gradient map of each stage were constructed through a topographic imaging system (Lai et al., 2016;Lai et al., 2017).

Automated measurement of channel position
Braiding intensity has been widely used in fluvial braided rivers (Egozi and Ashmore, 2008;Bertoldi et al., 2009;Egozi and Ashmore, 2009) and applied to submarine braided channels (Lai et al., 2017;Limaye et al., 2018).In past analyses of submarine braided channels (Lai et al., 2017), the braiding intensity required manual interpretation of cross-sectional channel numbers.
This approach is time-consuming and can be subjective.In this study, we developed an RGB color separation algorithm to identify channel positions, calculate cross-sectional channel numbers, and obtain an objective, repeatable, and averaged braiding intensity for a selected color image (Fig. 3).For instance, a total of 1080 color orthophotos were analyzed using an algorithm based on Eq. ( 1) that obtains enhanced gray-scale images: where ,  and  are weighting coefficients to be calibrated by each run (e.g.,  = 1,  = 0.8,  = 0.35 for Run A2).R, G, B represent the red, green and blue bands of a color image, respectively.The weighting coefficients used to generate the grayscale images varied slightly between experimental runs due to different light settings and dye concentration.Then, the grayscale images were converted into binarized images by an adaptive threshold method (Bradley and Roth, 2007) (Fig. 3b).After that the cross-sectional illumination values (from x = 300 to 1500 mm, extracted every 1 mm) were smoothed by a Gaussian filter to locate the channel positions (Fig. 3c to Fig. 3e, showing only three example cross sections).After the above procedure, the channel positions and cross-sectional channel numbers could be automatically obtained (Fig. 3f).Then the average value would be taken as the representative active braiding intensity (  ) of that image.The results of braiding intensity are presented in Section 3.2.

Reduced-complexity model
Reduced-complexity numerical models are often used to simulate macro-scale, long-term fluvial braided rivers (Murray and Paola, 1994;1997;Thomas and Nicholas, 2002).The model proposed by Murray andPaola (1994, 1997) is the most influential work on reduced-complexity models.The model predicts the generic dynamics of stream braiding, including channel shifting, avulsion and migration.Later research modified their flow routing schemes to make reduced-complexity models comparable to shallow water equations models (Thomas and Nicholas, 2002;Thomas et al., 2007;Nicholas, 2009;Williams et al., 2016).
Recently, the CAESAR-LISFLOOD model (Coulthard et al., 2013) was applied to map experimental submarine braided channels and bar geometry (Limaye et al., 2018).This analysis showed that compared to a subaerial braided channel with the same water and sediment fluxes, submarine braided channels developed fewer and deeper channels near the inlet, with declining topographic variations further downstream.In this study, we followed the numerical scheme of Thomas and Nicholas (2002) and developed a Matlab based reduced-complexity model, with modified hydraulic conditions for density currents.Our modifications assume large aspect ratios (width-to-depth ratio), reduced gravity and underflows entrain negligible ambient water, and that the flows of submarine braided channels may resemble fluvial braided rivers (Foreman et al., 2015;Lai et al., 2017;Limaye et al., 2018).
The model starts calculations when the discharge is provided in the cell at the upstream boundary.Then the discharge is transferred downstream according to the calculation path of each row, and the model step ends when the calculation reaches the downstream boundary.The discharge is distributed from each upstream cell to the five immediate downstream cells (Thomas and Nicholas, 2002) rather than three cells used by Murray andPaola (1994, 1997).The discharge could be transferred laterally up to 60 degrees to the downstream direction.The allocated discharge in the five downstream cells is given by: where  0 is the outflow discharge of the upstream cell;   is the routing potential, which determines the proportion of discharge that can be allocated to the five downstream cells.In the five downstream cells, the   of each cell was calculated from the upstream cell, starting from its water depth ℎ 0 .Assuming the density current is uniform, the flow depth at steady-state can be approximated by (Ippen and Harleman, 1952;Lofquist, 1960): where C is a coefficient to be calibrated by each experiment (C = 1.2 for this study);   = 1.68 × 10 −3 kgm -1 s -1 is the dynamic viscosity of saline water at 20 °;  is bed slope;   is inflow density;   is the density of ambient water;  is a dimensionless parameter, which is proportional to the ratio of interface velocity (  ) to maximum velocity (  ) of the density current, given by: For a uniform density current at steady state,     ⁄ ≅ 0.59 (suggested by Ippen and Harleman, 1952).Therefore, the calculated  value is 0.138 for this study.Flow depths (ℎ  ) at each of the five downstream cells were approximated from: where  0 is the bed elevation of the upstream cell;   is the bed elevation of each of the five downstream cells;  is the distance from the upstream cell to each of the five downstream cells.If the ℎ  calculated by Eq. ( 5) is a positive value, the flow can be routed; if ℎ  is a negative value, the flow cannot be routed.Therefore, the values of ℎ  affect the calculation method of   , which could be divided into the following two cases: (1) When at least one of the 5 downstream cells have ℎ  > 0, then   can be calculated.In the remaining cells, if its ℎ  ≤ 0, its corresponding routing potential will be set to 0.   is calculated as follows: where   is the local bed slope between the cells that flow is being routed from and to; the slope exponent 0.25 was suggested by Thomas and Nicholas (2002).
(2) If all the five downstream cells have ℎ  ≤ 0, the flow is assumed to be critical.This behavior often happens when density currents flood over sandbars, which is more likely to occur in submarine braided channels than in fluvial braided rivers (Sequeiros, 2012).Flow depths (ℎ  ) at each of the five downstream cells are then approximated by: where ℎ  is the maximum depth, ℎ  = � 0 � ′ ⁄ � 0.67 ;  ′ = (  −   )   ⁄ is reduced gravity;   is the lowest bed elevation in the downstream cell.The corresponding   is: The above rules together constitute our reduced-complexity model for simulating the flows of submarine braided channels.In our simulations, the initial bed condition is imported from the experimental DEM.The simulated inflow discharge is uniformly distributed at the upstream cells based on the experimental total inflow discharge.Notice that our reduced-complexity model only simulates the underflows at the stable phase (Stage 6 to Stage 10), i.e., the morphological changes are excluded.
Additionally, extra simulations are performed with double inflow discharge based on the experimental DEM at the final stage.
These simulations are used for testing whether the linear relationship still holds between dimensionless stream power and active braiding intensity under extreme events.These results will be addressed in Section 3.4 and Section 4.

Morphological evolution of submarine braided channels
Physical experiments document the initiation and evolution of submarine braided channels in response to set, known input parameters and conditions (Fig. 4).For instance, in Run A2, the initial confinement width was set at  = 0.24 m, and the inflow-to-sediment discharge ratio was     ⁄ = 58.7.By  = 1200 s, two small channels and a few sandbars appeared within the straight channel (Fig. 4b).At  = 2400 s, multiple channels and sandbars gradually developed downstream.At  = 3600 s, the flows showed frequent bifurcation and convergence.At  = 4800 s, slight bank erosion appeared on both sides and some channels were abandoned.At  = 6000 s, channels frequently change their routes and sandbars grew and elongated.The DEM hillshades clearly document the channel relief and bar thickness of submarine braided channels (Fig. 5).These high-quality DEMs provided the basis for subsequent morphometric analyses.Amongst Series A and Series B (Fig. 6), under the same     ⁄ condition, a higher confinement width produces more channels and higher degree of braiding.On the contrary, when confinement width is fixed, a larger     ⁄ generates fewer channels but wider main channels.For instance, in Series A (    ⁄ = 60), Run A3 (Fig. 6c) has more channels and higher degree of braiding.Similarly, in Series B (    ⁄ = 90), Run B3 (Fig. 6f) is more braided.However, when confinement width is fixed 220 (e.g.,  = 0.24 m, Fig. 6b and Fig. 6e), a larger     ⁄ would have less channels but wider main channels.These characteristics are also valid for confinement  = 0.12 m (Fig. 6a and Fig. 6d) and  = 0.48 m (Fig. 6c and Fig. 6f).Additionally, a larger     ⁄ also makes sand bars more connected and facilitates small sandbars to merge into longer and larger sandbars (Fig. S1 to Fig. S10 and Movie S1 to Movie S6 provide more detailed evolution processes).

Evolution of active braiding intensity
The evolution of active braiding intensity reflects the critical time for submarine braided channels to reach a stable phase and the influence of confinement width and inflow-to-sediment discharge ratio (Fig. 7).For instance, the evolution of the braiding index over time for each run included two phases: (1) a rising phase and (2) a stable phase.We interpret the time for development of the stable phase as a critical time (  ).Our results showed that   is proportional to , e.g., when  = 0.12 m,   = 2400 s; when  = 0.24 m,   = 3600 s; when  = 0.48 m,   = 4800 s.In addition, when     ⁄ is fixed,   increases as  increases, e.g., Series A (Fig. 7a to Fig. 7c) and Series B (Fig. 7d to Fig. 7f).On the contrary, when  is fixed,   at the stable phase decreases slightly as     ⁄ increases, e.g.,   of Series B at  = 6000 s are all slightly smaller than that of Series A. Therefore, a larger confinement width postpones the critical time and increases the active braiding intensity at the stable phase.However, under the same confinement width, a larger inflow-to-sediment discharge ratio slightly decreases the   .

Volume change and morphological active width
DEM of Differences (DoDs) demonstrate the sediment distribution patterns of erosion and deposition for submarine braided channels (Fig. 8).The results showed that confinement width and bank erosion was not proportional for Series A or Series B.
For instance, when  = 0.24 m, bank erosion on both sides was the most significant, especially for Run B2 (Fig. 8e).On the contrary, when confinement width is fixed, a larger     ⁄ would cause more severe bank erosion.Moreover, in Series B (Fig. 8d to 8f), the areas of erosion and deposition became more continuous and contiguous.In Series A (  ≅ 55 mm 2 /s,     ⁄ ≅ 60), most of the submarine braided channels develop within the given confinement width without much widening to accommodate the inflow.However, in Series B (  ≅ 84 mm 2 /s,     ⁄ ≅ 90), the submarine braided channels would widen the given confinement width to accommodate the larger inflow and cause stronger bank erosion.Besides, the DoDs also reveal the morphological active widths (  ) and bulk changes (  ) for each run.The relationship between   and   for both fluvial and submarine braided channels will be discussed in Section 4.

Simulated flow discharge and sediment transport rate
The reduced-complexity model for simulating submarine braided channels captured the pattern of flow bifurcation and confluence (Fig. 10).Modeled flow occupied the main channels that developed in each run of the experiment.In Series A (Fig. 10a to Fig. 10c), a larger confinement width increased the number, width and discharge of the main channels.In Series B (Fig. 10d to Fig. 10f), the main channels became wider and more significant.The simulated flow characteristics are similar to those qualitatively observed in the experiments.However, validating the reduced-complexity model by directly comparing the measured flow depths is unfeasible because the averaged depth of experimental density currents is less than 3 mm.Therefore, we converted the simulated and experimental flow maps (Fig. S11 to Fig. S16) into binarized images and compared their pixel differences.We examined the cross-sectional differences between the experimental and simulated flows every 1 mm along the  direction.The results (Fig. S17) showed that the errors are less than 30 % for  = 0.48 m (Run A3 and Run B3) and even less than 6.5 % for  = 0.12 m (Run A1 and Run B1).This indicates that our reduced-complexity model is a reliable and applicable tool for simulating submarine braided channels.affected by different confinement widths (Fig. 11).For instance, for Series A the modality of main channels belongs to unimodal distribution (Fig. 11a to Fig. 11c); for Series B the modality of main channels includes bimodal and multimodal distributions (Fig. 11d to Fig. 11f).Combining the simulated flow discharges and the bed load equation of Schoklitsch (1950), the sediment flux (kg/m-s) for each five downstream cells could be roughly estimated by: where   is local bed slope;   is flow discharge (m 2 /s);   is critical discharge (m 2 /s), given by: where G is the specific gravity of sediment (G = 1.5 for plastic sand);   is medium grain size (  = 0.34 mm for this study), which is suitable for the grain size between 0.31 to 7.02 mm; Note that there is no sediment flux when   <   .Kernel density estimation (KDE) confirms the probability distribution of simulated flows are robust and consistent.Under extreme events, the KDE provides the predicted pattern of cross-sectional discharge and corresponding   (Fig. 12).For instance, we extracted the simulated flow discharge at the middle reach (x = 950 to 1050 mm) to perform probability density function (PDF) and KDE statistical analyses.We excluded the data for which its unit discharge is less than 50 mm 2 /s to avoid too many outliers.The KDE shows that the PDF of simulated discharges are all positively skewed and less sensitive to confinement width.On the contrary, when doubling the inflow discharge, the skewness of KDE decreases, allowing the given discharges be allocated over a broader unit discharge.Therefore, the number of channels and corresponding   all increased.
Our reduced-complexity model appears consistent with the behavior of submarine braided channels, although the model still requires calibration and validation with more field bathymetric surveys.

Discussion
The purpose of this study is to understand the effects of confinement width and inflow-to-sediment discharge ratio on the evolution of submarine braided channels using physical experiments and reduced-complexity models.We propose that confinement width has a stronger influence than inflow-to-sediment discharge ratio on the morphological evolution of submarine braided channels.The scaling relationships established by experiments and simulations can help us better interpret the evolution of field-scale submarine braided channels.
First, we find that both active width (  ) and bulk change (  ) are proportional to confinement width ().Under a fixed reach length, a wider  would generate a larger active area.Then, the larger active area would give a larger active width.By definition (Ashmore et al., 2011;Peirce et al., 2018), bulk change is the sum of the absolute values of erosion volume (  ) and deposition volume (  ) in a DoD.A wider  increases the non-zero areas in a DoD, therefore increase the   .For example, the   = 0.12−0.14m for  = 120 mm (Fig. 13); the   = 0.24−0.28m for  = 240 mm; the   = 0.48−0.49m for  = 480 mm.Our results show that   is proportional to   , which agrees with the response of fluvial braided rivers (gray dots in Fig. 13).Even if the Reynolds number () and densimetric Froude number (  ) of our submarine braided channels indicate that flows in these channels are laminar and subcritical (Table 2), as long as the depth-to-width ratio (ℎ/) decreases and slope-to-densimetric Froude number ratio (/  ) increases, submarine braided channels will naturally form (Fig. 14).Our results agree with the morphological classification of fluvial rivers (Parker, 1976).In short, the active width and bulk change of both submarine and fluvial braided channels share the same trend.This implies that the two braiding systems may have similar physical mechanisms., where ℎ is the flow depth of the saline underflow, estimated from the experiments (ℎ ≈ 2 mm);  =   (ℎ) ⁄ is the cross-sectional averaged flow velocity;   is inflow total discharge;  is confinement width;  is the kinematic viscosity of water ( = 10 -6 m 2 s -1 at 20 ° C).
4 Reduced bed shear stress,   = (  −   )ℎ, where  is bed slope.Second, we demonstrate that active braiding intensity (  ) is proportional to both dimensionless sediment-stream power ( * * ) and dimensionless stream power ( * ), despite the dimensionless critical time (  * ) of submarine and fluvial braided channels are different (Fig. 16).The dimensionless critical time is defined as   * =       3 ⁄ (Lai et al., 2017).We found that the   * of lateral confined or unconfined submarine braided channels all fall around 10 8 to 10 9 (Table 2 and Fig. 16a); while the fluvial braided rivers fall around 10 11 .We think our smaller flume size, steeper bed slope and lighter bed material together make the submarine braided channels to reach   * earlier than that of fluvial braided rivers.In addition, the Shields parameter () corresponding to our sediment is about 0.46, which is much larger than the critical threshold for transport (~0.04 to 0.06) and larger than occurs in the fluvial braided rivers (Egozi and Ashmore, 2009) (Fig. 15).The dimensionless sediment-stream power, defined as the combination of width-to-depth ratio ( ℎ ⁄ ), sediment-to-inflow discharge ratio (    ⁄ ) and bed slope () (Fig. 16b, modified from Paola, 2001), is proportional to confinement width, sediment supply and bed slope.When   is fixed, a larger   reduces  * * , making the corresponding   become smaller.For instance, in Series B (  ≅ 84 mm 2 /s,     ⁄ ≅ 90), the submarine braided channels widen the imposed confinement width to accommodate the larger inflow, making the area of erosion and deposition in the DoDs more continuous and contiguous.This explains why increasing   slightly reduces the   of Series B at stable phase.The  * * is applicable to both unconfined and confined submarine braided channels (Fig. 16b).
On the other hand, the dimensionless stream power ( * ) considers density differences of the underflow, grain size (  ) and settling velocity (  ) (Fig. 16c), but does not include the effects of  and   .The  * works well for both unconfined fluvial braided rivers (Zanoni et al., 2007;Egozi and Ashmore, 2009;Bertoldi et al., 2009) and confined fluvial braided rivers (Garcia Lugo et al., 2015)(Fig.16c).In short, both  * * and  * are proportional to   .For an inverse problem, if   is known, these two relationships can be used to estimate the corresponding parameters (e.g., ℎ,   ,   , density differences), which provides a way to apply this work to paleoflow reconstruction.show strong linear relationships and support us to infer that the submarine braided channels on Baia di Levante Fan would continue to develop and prograde downstream.However, we think that the downstream lateral confinement may force the braided channels to transform into a single channel.The confinement width of these submarine braided channels gradually decreases from upstream to downstream and is well bounded by two erosive banks.The Baia di Levante Fan was formed by gradual stacking of gravity flows fed by repeated small-scale mass-wasting events rather than a single catastrophic collapse (Romagnoli et al., 2012).This modern field example agrees with the formation and morphological features observed in our laboratory experiments.Additionally, recent studies on the submarine flanks of La Reunion also confirm that long-term erosive processes combined with high sediment supply can form wide volcaniclastic deep-sea fans similar to the siliciclastic ones (Saint-Ange et al., 2011;Sisavath et al., 2011;Casalbore et al., 2021).Based on the reported bathymetry (Fig. 1b), when the water depth is 1000 m, the bed slope of Baia di Levante Fan is 4 degrees, the confinement width is 3500 m (Romagnoli et al., 2012), and the number of submarine braided channels is about 3. Assuming these three submarine channels are active, we obtain   = 3.Through the experimental scaling relationship (Fig. 16b), we obtain the corresponding dimensionless sedimentstream power is  * * = 2.5.Therefore, we can estimate the reasonable flow depth range of the turbidity current by assuming     ⁄ .For instance, when     ⁄ = 13, the estimated flow depth is 63.2 m.This depth agrees with the reported depth of the turbidity current (60 m) and the height of erosive walls (70−80 m) (Romagnoli et al., 2012).When     ⁄ = 9, the estimated flow depth is 91.4 m.It is possible to cause overspill and form sediment waves (or cyclic steps) outside the bank.Similar morphology appears on the left outer bank of Baia di Levante submarine fan.Overall, the morphologic estimates derived from our proposed scaling relations are consistent with field observations, using plausible values for the inflow-to-sediment discharge ratio.Constraining these inputs for field cases would be major step toward predicting the morphology of submarine braided channels.
Additionally, the experimental results are intriguing in that they further emphasize the similarities in the propensity for and properties of braiding in fluvial and submarine systems (Lai et al., 2017).In isolation, the experimental results would predict braided submarine channels are common in modern oceans, however, they are exceedingly rare (Foreman et al., 2015).Herein, we built upon previous experiments that documented the same relationship between dimensionless stream power and active braiding intensity among rivers and submarine channels (Lai et al., 2017) and expanded the experiments to show both the active width and bulk morphologic change increase as the confinement width () increases commensurate with experimental fluvial systems.Thus, braided submarine channels can be spontaneously generated both through self-channelization on an initially flat surface (Foreman et al., 2015;Lai et al., 2017;Limaye et al., 2018) and within a pre-existing confined channel (this study).Yet in nature both these bathymetric cases overwhelming result in sinuous, single-thread submarine channels.
This implies there may be a missing attribute within the experiments that in nature suppresses braiding.Lai et al. (2017) suggested the missing attribute was enhanced deposition of suspended sediment loads within inactive channel threads during low-flow conditions.In nature higher sediment concentrations within turbidity currents as compared to river fluid flows likely facilitate greater aggradation and in-filling of inactive channel threads.This process would limit segmentation of the flow by bars and promote development of a single thread over multiple threads.In the experiments there is no suspended sediment load.This phenomenon may be exacerbated in field submarine channels whose active width is confined by canyon or valley walls as compared to unconfined systems.In confined systems the fine-grained turbidity currents that overspill channel margins cannot be advected across the overbank area to zones more distal from the channel.Instead, the fine-grained material is likely retained and forced to deposit within the active width.Interestingly, the Baia di Levante Fan is comparatively coarse in grain size distribution (Ramagnoli et al., 2012), and in the absence of significant silt and clay components may represent a unique case in which a braided system may evolve.This suggests that provenance and grain size distribution of sediment supply may play a nontrivial role in the submarine channel planform morphology.Finally, the observation that as the confinement width, , increases the bulk morphologic change increases provides a testable hypothesis for stratigraphic studies of submarine channels.Namely, increased bulk morphologic change should produce more heterogeneous channel sand-body deposits, particularly related to the abundance bar reactivation surfaces.It may be productive for future stratigraphic studies to quantify this abundance scaled to the channel geometry.

Conclusions
Physical experiments and a reduced-complexity model allow us to better explain the responses of submarine braided channels to confinement widths and inflow-to-sediment discharge ratios.The rules of stream power are applicable to field-scale submarine braided channels.Our main findings are: 1. We confirm that a larger confinement width () postpones the critical time (  ) and increases the active braiding intensity (  ) at the stable phase.At a fixed confinement width, a larger inflow-to-sediment discharge ratio (    ⁄ ) slightly decreases the   .
2. The active width (  ) and bulk change (  ) are proportional to confinement widths for submarine braided channels, in agreement with fluvial braided rivers.
3. The experimental and modeling results agree that active braiding intensity is proportional to both dimensionless sediment-stream power ( * * ) and dimensionless stream power ( * ).
4. Kernel density estimation (KDE) helps to visualize and predict the pattern of cross-sectional discharge distribution.
5. We infer that the submarine braided channels on Baia di Levante Fan may continue to develop and prograde downstream, but that the downstream lateral confinement will eventually force the submarine braided channels to transform into a single channel.
6. Overall the relations we develop from the experimental and modeling results are able to predict the first-order general morphological and sedimentological patterns for field scale submarine fans.

Figure 2 .
Figure 2. (a) Experimental setup for studying the evolution of submarine braided channels with different confinement widths and inflow-to-sediment discharge ratios.;(b) inflow width, confinement width and valley width of Runs A1 and B1; (c) inflow width, confinement width and valley width of Runs A2 and B2; (d) inflow width, confinement width and valley width of Runs A3 and B3.115

Figure 3 .
Figure 3. Procedure of image processing for identifying channel positions for each experiment.(a) orthorectified color image; (b) converted binarized image; (c)-(e) extracted examples of cross-sectional illumination values and captured channel positions; (f) verification for the automatically captured channel positions and cross-sectional channel numbers.

Figure 4 .
Figure 4. Evolution of submarine braided channels of Run A2 from t = 0 to 6000 s.The confinement width is  = 240 mm.The inflow unit width discharge is   = 55.68 mm 2 /s.The inflow-to-sediment discharge ratio is     ⁄ = 59.

Figure 9 .
Figure 9. Normalized cross-sectional volume changes for all runs.

Figure 10 .
Figure 10.Submarine braided channels simulated by the reduced-complexity model for each run at  = 6000 s.

Figure 11 .
Figure 11.Simulated discharge, flow depth and sediment flux extracted at the midstream ( = 900-1000 mm) for each 300 run.The initial bed surface is based on the experimental DEM at the final stage.

Figure 12 .
Figure 12.Kernel density estimation (KDE) for the probability density function (PDF) of simulated flow distribution for each run.Color thick lines represent the KDEs of different runs.Black dash lines are the KDEs of double inflow discharge.

Figure 13 .
Figure 13.The trend of active width (  ) to bulk change (  ).Color dots represent data of submarine braided channels (this study).Gray dots are experimental data of fluvial braided rivers (Peirce et al., 2018).340 Table 2. Physical parameters and dimensionless numbers related to the submarine braided channels.

Figure 16 .
Figure 16.(a) Dimensionless critical time (  * ) of submarine and fluvial braided channels; (b) Relation of dimensionless sediment-stream power ( * * ) to active braiding intensity (  ) for two sets of submarine braiding experiments; (c) Relation of dimensionless stream power ( * ) to active braiding intensity (  ) for two sets of submarine braiding experiments, four sets of fluvial braiding experiments and one set of our reduced-complexity model.
with full valley width (), and the  in this study is based on the proportion of 21.8%, 43.6% and 87.3% of .The inflow unit width discharge is calibrated by   =    ⁄ , where b is inflow width (not confinement width ), because the submarine braided channels do not occupy the entire confinement width.The inflow-to-sediment discharge ratio ( Lai et al. (2017) (2015)h discharge (  ) were controlled at the same time.Based on previous reported successful runsof Foreman et al. (2015)andLai et al. (2017),     ⁄