Sediment transport at the river confluences: few observations from a sub-tropical plateau fringe river of eastern India

ABSTRACT Confluences are often attributed by abrupt change in the load transport and encompass significant regional differences. Tributary watersheds of the Chotonagpur plateau fringe fans are widely reported to be the source of premier sediment yields causing amplification of sediment load in the trunk river. In this direction this article scrutinizes the role of tributaries on sediment transport at 11 confluence locations along 86.4 km stretched river Bakreshwar. The work revolves around two fundamental questions- How do tributaries contribute sediment load to the main river respective to their catchment size, length, runoff and sediment yield volume? and how do tributaries' interaction with the trunk stream influence the longitudinal continuum of sediment load in the receiving river? The study affirms that the junctions represent locations of abrupt change in the longitudinal gradient of sediment load. Size and length of the tributaries relative to the main stream have explicit controls on the post-confluence change in sediment load. Tributary watershed sediment yield is the foremost factor for amplifications of sediment load and morphological complexities at and downstream of the junctions. This work is a preliminary attempt and can hopefully stimulate ideas for fruitful and coherent scientific inquiry on river confluences particularly for Indian circumstances.


Introduction
Within the fluvial networks, confluences are that critical nodes where two or more channels having independent flow and sediment regimes converge in a single channel downstream (Mosley, 1976). During the 1950s and 1960s Leopold and his colleagues (Leopold & Maddock, 1953;Leopold et al., 1964) theorized that a river from its headwaters to the mouth exhibits a continuous gradient of channel form variables and concave longitudinal profile in response to the downstream increase in discharge of water and sediment. According to Vannote et al. (1980), this continuous gradient in physical conditions along the river should result in a continuum of ecology along the river (River Continuum Concept, RCC). McManus (1985) revealed that in a relatively stable stream the volume of sediment flux generally maintains a balanced condition. However, the physical structures of the main stream are not "smooth monotonic function of distance downstream"  but exhibit step changes wherever significant tributaries join a channel (Rhoads, 1987;Richards, 1980) and many physical scientists (e.g., Montgomery, 1999) have challenged the "linear perspective" of riverine ecosystems (Benda et al., 2004a;Bruns et al., 1984). At the tributary junctions there are dramatic changes in discharge and sediment loads (Richards, 1980) that force adjustment to the interrelated channel form variables (Leopold et al., 1964).
Actually at the confluences of a river, tributary influx of water and sediment bring out sudden change/adjustment in the hydro-sedimentological condition of the recipient channel (Biron & Lane, 2008;Ferguson et al., 2006;Rhoads, 1987;Richards, 1980) leading to changes in physical conditions (Best, 1986;Miller, 1958). Therefore, there is a need to evaluate the confluences as discontinuous steps in the longitudinal continuum of the flow and sediment discharge regimes of a river.
The step changes in the physical gradient of a river at the confluences reflect dynamic interaction between tributary and main stream's water and sediment. As the tributaries have their typical catchment character and their get-together with the main stream lead to the change in the sediment load and consequent hydromorphological change of the receiving river downstream. These changes have a causal link to the tributary catchment dimensions likebasin size and shape, length of the tributary, catchment runoff, sediment yield, drainage density and pattern, spatial distribution, local geometry and so on (Benda et al., 2004b). Remarks such as "streams are the manifestation of the landscapes that they drain" (Hynes, 1975); "A river is a report card for its watershed"; etc. affirm this causal link. Despite this ubiquity of catchment influence on water and sediment transport within the confluence hydrodynamic zone it is rather surprising that, only a handful of studies have paid consistent attention to this issue. Study of Mueller (1972), Mesa and Mifflin (1986), Robinson et al. (1995), Benda et al. (2004b), Macnab et al. (2006), and Sklar et al. (2006), Ghosh (2016aGhosh ( , 2019a are some of the important exceptions that have paid attention to the confluence hydro-sediment dynamics viewed from the tributary catchment context. Besides, study of Richards (1980), Bruns et al. (1984), Rhoads (1987), Ward and Stanford (1995), Montgomery (1999), Ferguson et al. (2006), , and Duncan et al. (2009) have highlighted confluences as the point of discontinuity in the hydrosediment continuum within a fluvial network. Apart from these, study of Best (1985Best ( , 1986Best ( , 1987Best ( , 1988, Biron et al. (1993), Boyer et al. (2006), Ghobadian and Bejestan (2007), Szupiany et al. (2009), Ribeiro et al. (2012), Martin-Vide et al. (2015), and Nazari- Giglou et al. (2016) are also important in this regard to explore the nature of hydro-sediment dynamics and related issues and explanations at the confluence hydro-dynamic zone. Works on the related ground is very few in Indian circumstances.
Small to medium watersheds (10-100 km 2 area) of the Chotonagpur plateau fringe fans (eastern counterpart), breed a large amount of runoff during June and September . The high intensity seasonally skewed rainfall during the monsoon (Ghosh, 2018) is primarily responsible for the sharp increase in seasonal discharge of the rivers (contributing as much as 90% to the annual discharge) unlike other rivers of the gangetic flood plain (Bandyopadhyay et al., 2014). Moreover, the undulating lateritic terrain conditions helps to generates some of the premier sediment yields of Eastern India which added to the main river from the tributaries. Such area-specific high sediment yields are primarily due to erosion propensity of the secondary laterite formation in this counterpart (Chakraborty, 1970;Jha & Kapat, 2003). Subsoil sediments produced by gully incision in Lateric patches of this region provide most of the sediment load to the headwater sub-catchments (Ghosh & Bhattacharya, 2012) and supply it to the main stream. For example, the river bed of the lower-middle and lower Bakreshwar River is being aggrading (Ghosh, 2015;Ghosh, 2016c) owing to the high rates of sediment supply particularly from the tributaries. Such high rates of sediment influx from the tributaries have both on-site effects as well as offsite impacts. Therefore, a comprehensive account concerning addition of sediment from the tributaries, sediment load dynamics of the main river at the confluence hydrodynamic zone and its causal link with the catchment scales river network functions is needful for scaling-up our insight so as to understand the regional spirit.
In view of the above, the present paper concentrates on a few tributary junction sites along 86.4 Km. stretched river Bakreshwar having western extension to the Chotonagpur plateau fringe and eastern extension towards the Gangetic flood plain of eastern India. The work revolves around two fundamental research questions -(1) How do tributaries its interaction with main stream influence the longitudinal gradient in sediment discharge of the main river within riverine settings? and (2) How do tributaries contribute sediment load of the main river respective to their physical settings namely catchment size, basin order, main stream length, catchment runoff, and sediment yield volume?
The region bears the imprints of Late Pleistocene ferruginous formations with the remnants of tropical deciduous forest. Hard massive Basalt of Jurassic to Cretaceous age, soft and medium hard laterites of Cenozoic age and China Clay of Late Pleistocene to Early Eocene age are also found in different parts of the basin in different depths where deep percolation is negligible (Saha, 1961). Spatially, 35% of the basin area in the upper catchment is composed with granitic gneissic rock of Pleistocene age (50 lakh years old) overlain by coarse grain lateritic soil; 27% and 38% area at the middle and lower catchment respectively are made up with older and newer alluvium of Holocene period over granitic basement (Figure 1). The catchment area covers mostly with the reddish, loose and friable sterile soil with ferruginous concretion called laterite soil (local name "kankara"). Due to strong infertility of land it is named as "Rarh" or "sterile land" (O'Malley, 1914). This secondary laterite formation (Chakraborty, 1970) is primarily due to the deposition by the small tributaries coming from Chotonagpur Plateau (Jha & Kapat, 2003, 2009). The soil catenas consist of plateau fringe and high plain with Laterite Soil (Ultisols), adjacent slopes with sandy and loamy soils and small valley floors with older alluvial soils (Ulfisols). No important lineaments are there but few locally confined aquifer structures are found in this region. The thickness of alluvium (ranges from 12-20 m) increases towards the eastern part of the catchment and groundwater yield potentialities of the aquifers also increases eastward with a rate about 5-15 m 3 /h (Ray & Shekhar, 2009). The water table is moderately deep (5-10 mbgl) with moderate to high seasonal fluctuation (Mukherjee et al., 2007). Flood incidence is common particularly in the lower catchment during peak monsoon period. Discharge from Nilnirjon dam and additional influx of water through canals from other river systems encourage flood incidences (Ghosh, 2015) and bank erosion (Ghosh et al., 2016;K.G. Ghosh, 2016c) particularly in the lower part of the river basin.
Agricultural land use (>53%) dominates the region. Permanent forest cover lies exclusively in the lateritic upper catchment area constitutes only about 16% of the basin area and is decaying both qualitatively and quantitatively at a very fast rate. Moderate-maximum chemical weathering, moderate physical weathering, moderate mass wasting, moderate-maximum fluvial processes (sheet wash, rain splash, rill and gully erosion), laterisation, etc. are some major pedo-geomorphic processes. Rolling lateritic uplands, wide planation surface, deep red weathering zone, duricrust of Feoxides (Ferrucrete), badlands in the upper part and low lying flat land in the lower part of the basin characterize the morphological features. The region The selected tributaries have delineated by their respective IDs (t1, y2 . . .) and the surveyed sample confluence segments have indicated by red square and labelled correspondingly as S1, S2, . . . The spatial geology shown in the map is based on the Geological Quadrangle map (73 M) of the Geological Survey of India (GSI). experiences sub-humid and sub-tropical monsoon type climate with alternative wet and dry period where rainfall (varies with the elevation between 925 and 2178 mm during 1901-2017) shows declining trend and the drought condition is reported to be intensified (Ghosh, 2018(Ghosh, , 2019b. Average potential evaporation of this area from 1901 to 2015 is 73.45 mm/year. Presence of laterite indicates the former existence of tropical wet dry periods (Tertiary). During South-West Monsoon (June to September) advective rainfall yields more than 80% out of the total annual average rainfall (1383.84 mm) and all the rivers (tributary and main) receive copious amounts of rainfall during these months and thereby intensify the problem of soil erosion. It is found that, typical heavy monsoon rainfall intensity of around 23.51 mm/h (WBPCB: West Bengal Pollution Control Board, 2009) and median drop diameter of almost 2.2 mm with terminal velocity of about 6 m/s (Hudson, 1965) in this area enhancing rainfall erosivity (Ghosh, 2016b). In the dry period the western highlands experience excessive dry condition causing physical weathering. Average estimated soil loss of this basin area is about 7 tons/ha/year . During the monsoonal months particularly in the undulating western extreme heavy runoff (runoff coefficient in monsoon is as high as 0.705)  is the prime cause for adding the eroded sediment to the drainage networks. Owing to these high rates of sediment supply the river bed of the lowermiddle and lower Bakreshwar river is aggrading (Ghosh, 2015). Hence, attention should be given upon tributary contribution of sediment load to the main river Bakreshwar and associated causal effects.

Materials and methods
The present study dealt with the fifth order Bakreshwar river (Area: 714.56 sq. km.; Length: 86.4 km.) and its major tributaries (Figures 1 and 2) located in the eastern India. The catchment area is primarily a part of the farthest eastern extension of the Chotonagpur Plateau in the West and Moribund Delta in the East and often referred to Rarh Bengal (Bagchi & Mukerjee, 1983;Biswas, 1987). The whole study is based on 2 fundamental data types (Table 1): Channel dimensions related data which is collected from field survey and the catchment dimensions related data which is secondary in nature and mostly calculated in GIS environment. In order to scale-up our understanding of the catchment level fluvial dimensions, important characters of the tributary catchments ( Figure 2 and Table 2) are examined scientifically.
Data pertaining to primarily the discharge and sediment transport within the confluence hydrodynamic zone is obtained from subsequent field survey (   The other catchment scale variables that are supposed to influence the sediment transport within the confluence hydrodynamic zone are soil erosion and runoff. The Revised Universal Soil Loss Equation, RSULE (Renard et al., 1997(Renard et al., , 1991 and SCS-CN (NEH-4, 1954;USDA, 1972) methods are used to show the spatial pattern of soil loss and runoff, respectively. The rational method (Goyal, 2016, p. 132) and Garde and Kothyari's equation (Garde & Kothyari, 1987) are followed for the calculation of Catchment yield (i.e., runoff) and the volume of sediment yield for the respective tributary catchments and the main river catchment area up to the confluence junctions, respectively.
Catchment yield (runoff) is obtained by assuming a suitable runoff co-efficient.
where Crunoff coefficient Aarea of the catchment (in sq. m), P -Precipitation (in m.). The values of runoff coefficients -C fluctuate depending upon the soil type, vegetation, geology, etc. As the present catchment area is primarily a rural tract (Pal, 2016) therefore a uniform C value of 0.4 have taken for the calculation (vide Stephenson, 1981, pp. 39;Goyal, 2016, p. 132).
Amount of sediment yield is estimated using Garde and Kothyari (1987) where A is catchment area (km 2 ), H is catchment slope, Dd is drainage density, P is annual mean precipitation (mm) and Fc is erosion factor (vegetative cover factor). Again FA is the area of arable land in the catchment, FG is the area occupied by grass and scrub, FF is the forested area while Fw is the area of waste land.
The basic characteristics of the tributaries and their respective catchments are forwarded in the following table (Table 2). Strahler, 1957) Where, Bifurcation ratio, Rb = Nu/(Nu+1) and Nu is the Number of stream segment of a given order. Nu+1 is that of the next higher order. Strahler, 1952) where, Ls -Proportion of uneroded part; Lus -Proportion of eroded part Smith, 1935) Average slope; tan ϕ ¼ N Â I 636:6 ½9 (Wentworth, 1930) Where, Nnumber of contour cuttings per miles or km, Icontour interval, 636.6constant. 3.3. Summarizing the channel dimensions related data collection methods To have a better insight into the flow and sediment influx at the tributary junctions focus has been given on the stream reach meaning thereby segment of a river with fairly uniform size and shape, water surface slope, channel materials and flow characteristics (Dingman, 2009). Apart from these, variables such as underlying geology, perenniality of the river segment, stream confinement (or valley width), stage of channel evolution, dominant channel materials, major human intervention, etc. are (Harmar & Clifford, 2007;Montgomery & Buffington, 1997;NEH, 1988;Paz & Collischonn, 2007;Rosgen, 1998) considered for reach identification. The entire river stretch from source to outlet is divided into seven total reaches ( Figure 1) and the flow line diagram of the river Bakreshwar and its tributaries are shown in Figure 4. Ten major tributaries are identified with their corresponding 10 junctions with the main river Bakreshwar. At least one sample tributary junction (Table 3) from each reach is preferred for study. Three cross-sections from each sample confluence segment (two on main river: pre and post-confluence and one on the tributary river) have taken and their corresponding hydraulic (flow), sediment transport and morphological (form) attributes occurring in the confluences are also being investigated. Additionally, we have decided to investigate the confluence segment of the river Bakreshwar with the main river Kuya (up to the confluence tip area: 509.75 sq. km. and length: 110.09 km). Altogether, 7 reaches, 11 confluence segments, and 33 cross-sections are investigated. Rigorous field survey has been carried out with the help of four local field personnel to unveil the changes in the sediment load character at each confluence hydrodynamic zone. As stated earlier cross sections across the river are taken at the immediate upstream and downstream of the confluence points during the post-monsoon period. The sediment influxes at the confluences are measured for both under steady flow conditions in the monsoon (during June-September) and dying flow conditions during post-monsoon (October-January).
From each confluence hydrodynamic zone two water samples (one from upstream, one from downstream of the confluence junction) of handful amount (1 L sample) are taken in each time of surveying (during September 2011 and September 2017) and sediment load concentration (Both suspended i.e., fine sediments that are transported in suspension and dissolved i.e., sediments that are carried in solution) excluding bed load are computed in laboratory setup following standard manual. In this case, precise weight of the collected samples, using a digital weighing machine are measured and compared with that of the equal volume of fresh water. Thereby, mean sediment concentration and standard deviation of the samples are calculated. The velocity (v) is measured in the field using Digital Current Meter. Volumetric flow rates or Discharges (Q, in m 3 /s) are calculated using formula, Q = A(Cross section area)*v(Velocity). For calculation of suspended sediment discharge the following formula (Gray & Simoes, 2008) is used: where, Q s -Suspended sediment discharge in metric tones; Q w -Water discharge in m 3 /s; C s -Mean concentration of suspended sediment in the cross-section in mg/L; Ka coefficient based on the unit of measurement of water discharge which is 0.0864 when discharge is in m 3 /s. After combining the water and sediment discharges data collected form subsequent field survey for each tributary junction the power function model of the form Yc = aX b that relates the sediment concentrations to water discharge is estimated. The equation is- where, Q L -Suspended sediment discharge in (mg/L), Qw -Stream discharge in (m 3 /s), a and b are the coefficients that provide the best relationship between discharge and the sediment load. In this case, a and b are accepted as constant throughout the calculations. After log-transformation to the arithmetic domain and application of the Ferguson (1986) correction factor, the sediment load occurring at a specific discharge is estimated using the following expression: where, C f is the log-transformation bias correction factor- where, exp is the exponential function and s is the standard error of the regression equation.

Catchment scale characteristics of the sample tributaries
From Table 4 and Figure 5(a) it is quite evident that, only in three cases the drainage area of the tributary river is >60% of the drainage area of the main river. The first is at the S1 (Bansbuni) when the river Bakreshwar is receiving the first significant tributary, second is at the S5 (Kanspai) where tributary Chandrabhaga (areal share is maximum in compare to other tributary catchments) and third is at location S11 (Milonpur) where the river Bakreshwar itself became tributary of the master river Kuya. It is interesting to note that the drainage area of the tributary Bakreshwar (714.56 km) is more than that of the river Kuya (Kopai) up to the confluence tip (509.75 sq km) at S11. In case of the relative drainage length ( Figure 5 (b)) the picture remains almost the same for other sample tributary junctions except S11 (Milonpur). At this confluence location (S11) the length of the main river Kuya (110.09 km) is more than that of the tributary river Bakreshwar (86.36 km) although the catchment area of the river Kuya is less than its tributary Bakreshwar.
From the analysis of the drainage areas and length of the tributaries in compare to that of the main river up to the confluence location, it can be seen that the size of the tributary basins and length of the tributaries in compared to that of the main river at up to the confluence nodes decreases downstream, as illustrated in Figure 5(a,b). This behaviour indicates that downstream tributary river influence on channel hydrolosedimentological modification of the main river decreases downstream.
On the basis of previous studies (Rhoads, 1987;Roy & Woldenberg, 1986) it is postulated that sediment transport within the confluence hydrodynamic zone of main stem channels depends on the size of the (3) (8) (1) Reach number/ID successively from upstream to downstream, (2) Length of a reach, (3) Cumulated drainage area at the end of a reach, (4) Sample segment name is alike the nearby settlement, (5) Id/name of the selected tributary, (6) Typical form/processes occur in a reach, (7) Dominant channel material, (8) Longitudinal slope is the ratio of unit drop and distance, (10) Average stream bank top to top width, (11) dominant underlying Geology and (12) annotations tributary relative to the size of the main stream. For instance, larger basins typically produce larger quantities of sediment, so that larger tributaries generally create larger fans (Bull, 1964;May & Gresswell, 2003). In addition, the length of the tributary channel is supposed to be associated with increase of sediment load and energy downstream that would affect the ability to transport tributary inputs of sediment downstream of the confluences. Since main stream channels increase in size, length and energy downstream, increasingly larger tributaries should also be required to create significant change of the hydro-sedimentological panorama of the main river. Another two forces that are thought to have significant impact in the confluence hydro-dynamic zone are runoff volume and sediment yield volume of the tributary basins (Table 5). The relative contribution of the runoff volume and sediment yield volume for tributary and main river drainage areas up to the confluence tip are shown in the Figure 6(a,b), respectively.
The spatial patterns of runoff and soil erosion as well as relative strength of runoff and sediment yield volume for tributaries and main river catchment areas up to the confluence tip are shown in the Figures 7 and 8.
Two interesting factors to note down both in case of the runoff volume and sediment yield First: the upper and upper middle catchment tributaries contribute greater share of runoff and sediment yield than that of their areal/length share and the   circumstances get reversed in case of the lower and lower middle catchment. It means the tributaries located in the lower middle and lower catchments, particularly the sediment yield volume is relatively less than their areal contribution. This signifies that from upstream to downstream the possibility of addition of water and sediment to the main river from the tributary river respective to their size and length slow down.
Second: The left hand side tributary river and its catchment constitutes relatively more areal share (59%); length share (51%) as well as runoff volume (60%) and sediment yield (58%) ( Table 6). This signifies that the probability of the addition of water and sediment will be more in case of the left hand side tributaries than that of the right hand side tributaries. The field results (Table 7) show that the velocity is higher in the immediate downstream of the confluence in comparison to the immediate upstream in all the cases  (ranges between 02% and 68% increase). The increase of velocity and discharge is maximum in the Kopai-Bakreshwar confluence followed by Chandrabhaga-Bakreshwar confluence and minimum in Khagradanga where an almost dry channel eventually meets the river Bakreshwar (Figures 9 and 10). However, it is interesting to note that, the increase of discharge (ranges between 10.95% and 58.65%) is relatively more to that of the increase in velocity (ranges between 02% and 68% increases) in most of the sample locations (except Milonpur and Tekadda) along the main channel in post-confluence conditions. This anomaly may be due  to 15.78% to 23.39% increase in sediment concentration and cross sectional area as well as hydraulic radius after the confluence point. The eventual increase of velocity in comparison to the discharge in Kopai-Bakreshwar confluence and Tekadda may be attributed to the additional supply of water from river Kuya and irrigation canal network (Ghosh, 2015), respectively.

Downstream pattern of sediment load and sediment discharge: Pre-post confluence dynamism
Due to successive addition of eroded material from several non point sources, sediment concentration generally increases downstream. However, discontinuities in the fining of sediments are often observed at tributary junctions (Rice, 1998). Moreover, due to the anthropogenic restriction in form of dam (e.g., Nilnirjon dam in present case) and bridge culvert (e. g., Hatikra Bridge in present case), some amount of sediment settles down and thereby lowers the load. In some cases large boulders in the bed and minimum slope of the bed eventually helps to settle down the sediment load ( Figure 14). On an average sediment load grossly increases downstream and the tributary confluences show some anomalous change in the sediment load character. This also leads to a sudden change in channel bed morphology.
The results (Table 8) indicate that all the tributaries have their influence in increasing the suspended sediment discharge of the main river by −14% to 48%.
Besides, the mean grain size of bed load varies from fine sand to medium sand in each upstream and downstream of tributary junctions, respectively.
From Table 8 and Figures 11-13 it is quite evident that the sediment concentration and sediment discharge of the river Bakreshwar has changed abruptly at most of the tributary junctions. The standard deviation of the sediment load for individual locations at the downstream of the confluence is much higher than that of the upstream. This affirms that tributaries have their impact to break the longitudinal rhythm of the hydro-morphological panorama of the main river. Another two points should be noted: One: where percentage increase of water discharge is high the % increase of sediment load is quite low; Two: in all the cases % increase of sediment discharge is more than the % increase of sediment load.

Downstream pattern of sediment load: some locations with irregularities
There are few confluence zones in the present case where the nature of change in the sediment concentration has shown some abnormality. Out of the 11 sample confluence segment three of such location deserves worth mentioning: First: In most of the cases velocity is found significantly higher in the immediate downstream of the confluence in comparison to the immediate upstream (Table 7 and Figure 9) and the percentage increase in  velocity ranges from 2.8% (Khagradanga, S9) to 68.6% (Milonpur, S11). However, the insignificant increase in velocity (8.6% increase) and substantial amount of decrease in sediment load (−13.7%) at Bakreshwar (S2) may be attributed to the dissected channel pattern ( Figure 14) as a result of positioning of exposed hard rock (Granite and Quartz type) at different alignment in the river bed which virtually helps to settled down the sediment load.
Second: In Hatikra (S6) both the sediment concentration (−14.3% decrease) and suspended sediment discharge (+3.3% increase) is abnormally low in the post confluence section then that of the velocity (+11.9% increase) as well as water discharge (+20.6% increase). The field records show that there is an elevated bridge culvert for the Bolpur-Suri State Highway (SH) which has exerted considerable impact on the free flow of water and some amount of sediments settled down. The nearby brick kiln industry has also caused destruction of riparian areas, sand excavation and associated channel course shifting (Figure 15) which eventually helps to increase the sediment concentration in the pre-confluence section.
Third: In Tekadda (Figure 16) velocity and water discharge increases to about 15%. However, in comparison to that sediment load counts an abnormal decrease (−6.6%). This is because of the fact that, substantial amount of relatively freshwater is added to this point (S8) through the tributary via Kopai-Bakreshwar irrigation canal network (Ghosh, 2015). Note.: a and b implies pre-confluence and post-confluence results, respectively. At S11 Bakreshwar is the tributary river and Kopai is the main river and the combined flow is known as the river Kuya (Kopai). % increase in Sediment Concentration and discharge has been obtained by calculating the net change in post-confluence segment [simple subtraction of post-pre confluence condition) and thereby converted the change into % in respect of the pre-tributary condition. Simply (Pre-post channel values)/preconfluence condition*100] Figure 11. Sediment load (suspended and dissolved) and Sediment discharge pattern in the pre-confluence and post-confluence segments of different sample sites from source to mouth of the river Bakreshwar.
4.2.4. Downstream trend of sediment concentration: break in the longitudinal continuum From Table 8 and Figures 11 and 12 it is evident that the sediment load is gradually increasing downstream. Active liquefaction process has made the water more turbid downstream. There is wide variation in the sediment concentration up and downstream of the Nilnirjon Dam. The average sediment load during mid-monsoon at the upstream sites of the dam is 2.27 mg/L. But, it is 3.501 mg/L at the downstream sites. At the upstream up about 5-7 Km. due to the arrest of water sediment concentration is quite more. In Kanspai (S5) and Milonpur (Kuya) (S11) sediment concentration is high due to mixing of water from   tributary river Chandrabhaga and Bakreshwar, respectively. At Tekadda (S8) influx of water from Mayurakshi Canal has raised the sediment concentration slightly.   In Table 9 downstream trends of sediment loads and sediment discharges are synchronized with the theoretical trend.
4.3. Influence of catchment scale characteristics of the sample tributaries on reach scale flow character at the confluences Quantitative statistical relation diagram (Figures 17-21) may help to determine the role of tributaries, particularly their relative area, length, runoff, and sediment yield volume to the increase/decrease of the sediment load and sediment discharge in the immediate downstream of the confluence segments.
The linear regression relationship between the tributary drainage area relative to the main stream drainage area up to the confluence tip (in %) and percentage increase in sediment load and sediment discharge in the post confluence segments ( Figure 17) shows that the greater the drainage area of the tributary more the possibility of increase sediment load and sediment discharge.
Marginal positive relation between the tributary river length relative to main stream length up to the confluence tip and sediment load/sediment discharge (Figure 17(a,b)) signifies that greater the length of the tributary higher the amount of sediment load it receives from different non-point sources.
Tributary river runoff volume relative to the main river runoff volume up to the confluence tip positively and significantly controls the sediment load and sediment discharge (Figure 19(a,b)).
The control of tributary basin sediment yield volume (Figures 8 and 20) relative to the main river catchment sediment yield volume is found to be the most influential factor (R 2 = 0.683 and 0.681) to regulate post-confluence increase in sediment load and sediment discharge.
In comparison to the immediate upstream sediment load and discharge to that of the confluence downstream sediment load and discharge there is reasonable differences. The net difference and percentage change in the confluence downstream in compare to the upstream have a good relation with the catchment size, its location, length of the tributary, individual catchment runoff, and sediment yield. However, their functional relations are more prominent in case of the sediment discharge than that of the sediment loading. Their statistical relations are forwarded in the following correlation matrix Table 10.

Some other points of concern
Combining all the discharges measured along the main river in pre and post confluence segments during the entire field visit and thereby calculating the sediment discharge of the river in pre-and post-confluence segments and by this means plotting (Figure 21) with the tributary catchment sediment yield volume in respect of their annual expedience probability it is found that, tributary catchment sediment yield volume has a significant good relation with the net increase of sediment discharges in the post confluence segments. Therefore, the underlying fact behind (Table 11) the tributary catchment sediment yield volume to become the most influential factor to regulate the increase of sediment discharge at the post confluence segments are needed to be explored. Light textured lateritic soils of this region having shallow depths are very prone to erosion in case of high intensity rainfall under the undulating terrain conditions. It has been reported that, typical heavy monsoon rainfall intensity is around 23.51 mm/h in this area (Water Resource and its Quality in West Bengal, A State of Environmental Report, WBPCB, 2009) and median drop diameter of monsoon rainfall is almost 2.2 mm with terminal velocity of almost 6 m/ s (Hudson, 1965).
The frequency distribution of different intensity ranges of erosive events show an increasing trend as  the number of short durated heavy rainfall events has increased, particularly from the mid-1990s. As soils of the region are shallow and lands are undulating, of considerable concern is that so much surface erosion propensity at the catchment scale and sediment transfer through the tributary can be affected by so small a quantity of water, typically, 6-8 mm. This is an intense landscape response to weather events that reflects the sensitivity and fragility of the sub-tropical plateau fringe environment. Given the degree of erodibility and the efficiency of sediment transfer, hydrologists need to take heed that impounded water resources are at risk from high rates of sediment generation in this region.

Conclusion and future research prospects
At a river confluence, an abrupt meeting of two channels each having independent flow and sediment discharge regimes creates unique hydro-sediment and morphological conditions. In the present study, based on field data obtained from 11 confluence locations along 5 order of magnitude river Bakreshwar with the catchment area of about 714.56 sq. km located over the Chotonagpur plateau fringe fans of eastern India there has been an attempts to explore the sediment load dynamicity of those critical nodes with the objectives to anticipate the effect of tributaries to create discontinuities in the longitudinal rhythm of sediment load as a Figure 21. Plotting of tributary watershed estimated sediment yield volume and net increase of sediment discharge in post confluence segment against the Probability [p = r/(n + 1), where r is rank order and n is the total number of events of exceedance. Numerical values of the exponents are shown to be informed of the power function equation ax b . Where, values of the exponents "a" and "b" are given individually, and x represents water discharge in cubic metre per second (Cumec); However in case of specific energy and roughness factor linear regression model has been used due to the −ve value of the dependent variables function of spatial scale tributary catchment personality. The analysis of 11 confluence junctions from the 86.4 km stretched river Bakreshwar, constituting 7 channel reaches and 33 cross-sections in total, show that the tributary junctions represent areas of point specific change in the longitudinal gradient of sediment load. The study affirms that the tributaries have played a decisive role to locally alter the main stem sediment load respective with their catchment size relative to the main stem. Tributary attributes such as size and length, relative to the main stream up to the confluence tip found to have specific controls on the post-confluence change in sediment load. Tributary catchment sediment yield and runoff volume is found to be the most influential factor to regulate sediment transport within the confluence hydrodynamic zone. However there are several other vectors in addition to those discussed in this article that could control confluence effects as well. All of these vectors undoubtedly contribute to the variation exhibited in the data. Assessment of the role of catchment scale character of the tributary river on post confluence change in sediment concentration is thought to be the keystones for scaling-up our understanding of riverine processes to the drainage network scale. The analytical tools and approaches described in this article constitute a preliminary attempt towards understanding the said link using empirical data sets. Several other issues are yet to be explored. For instance, recent modelling work on sediment routing at tributary junctions illustrates that longitudinal discontinuities in physical and biological diversity can be predicted at network scales Rice et al., 2006). Such predictive models can be tested against empirical patterns revealed in this work. Moreover, spatial analysis of river networks is currently a major focus of environmental monitoring and assessment at local and regional scales and can provide useful tools for understanding and modelling tributary impacts in a network context (Peterson & Urban, 2006;Theobald et al., 2006;Ver Hoef et al., 2006). Additionally some other issues like changes in hydrological and morphological condition at the confluences are yet to be explored. The author is working in this aspect and will be reported in the subsequent study. The step-like changes in discharge and sediment load as outlined in the present piece of writing may be accompanied by discontinuities in longitudinal profile (Hanks & Webb, 2006;Rice & Church, 2001), changes in planform style (e.g., Galay et al., 1998;Mosley & Schumm, 2001;Russell, 1954), the growth of tributary fans (e.g., May & Gresswell, 2004), shifts in fish and insect community composition (e.g., Fernandes et al., 2004;Kiffney et al., 2006; and other geomorphological and ecological impacts (Ferguson & Hoey, 2008;Rice et al., 2008) are also needed to be examined for this particular geographical region. Present analysis of the effects of spatial scale tributary catchment character could also be applied to other riverine landscapes of the Gangetic floodplain to unveil the true regional character.
Furthermore, viewing confluences as the discontinuous steps in the longitudinal continuum of flow and sediment discharge regimes of a river and critical nodes of complex hydraulics, routing of water, sediments, pollutants as well as ecology have apparent implications for the management and restoration of rivers. For instance, sediment fluxes what Rice and Church (1998) called "sedimentary links" from tributaries which become particularly dramatic during significant flooding episodes causes abrupt change of the main-channel form (river bed aggradation/degradation, bank erosion, etc.) and therefore has obvious implications for the assessment of bank erosion and channel stability above and below significant tributary junctions. These main stream river problems are possible to manage by the careful manipulation of  Fournier (1960) index of concentration of rainfall (p. 2/P) is a measure of intensity of rainfall, where prainfall in the month with greatest precipitation (mm) and Pannual rainfall (mm).
tributary characteristics. Moreover, identification and mapping of geomorphically significant confluences based on relative strength of the tributaries and significant flow-related changes could contribute to identify probable ecological "hot spots" by identifying particularly vulnerable aquatic resources which in turn could guide for restoration efforts. The work reported here is thus will hopefully bring forward valuable knowledge base and stimulating ideas for fruitful and coherent future work on river confluences particularly in Indian circumstances and thence future efforts could be directed for management and restoration of the confluences.

Acknowledgments
An earlier version of this article has been presented at the 30th National Conference of 'The Indian Institute of Geomorphologists (IGI)' on 'Geomorphology, Environment and Society', held at the Department of Geography, Jamia Millia Islamia, New Delhi, India during October 2018 and the author likes to express his sincere thanks to Prof Sutapa Mukhopadhyay of the Visva-Bharati University and Dr Swades Pal of the University of Gour Banga, West Bengal, India for their reasonable suggestion during initial draft of the manuscript. In addition, the author is also thankful to all the dedicated field crews, especially Dr Abhishek Bhattacharya of Kalna BEd. College; Dr Surojit Let of Hetampur College; Dr Priyanka Das of Wada Women's College; Dr Sunil Saha of the University of Gour Banga and Dr Kaustuv Mukherjee of Chandidas Mahavidyalaya who spent long hours walking kilometres to collect the data that were necessary to develop the ideas, concepts and analyses in this article. Thanks should also go to the local field personnel, for their fair and priceless field assistance. Lastly my heartfelt thanks go to the anonymous reviewers for their constructive comments and suggestions to improve the manuscript.

Disclosure statement
No potential conflict of interest was reported by the author.