4.1 Sinuosity Index
Natural River channel usually shows unpredictable behavior thus to identify their patterns within the floodplain sinuosity has been identified. The Sinuosity Index (SI) was determined by digitizing the river channel's centerline to establish a consistent river pathway (Stream Length), which differed between 53.67 km and 52.82 km, as illustrated in Fig. 6. Additionally, an imaginary representation of the river channel's pathway (Valley Length) was considered. The resulting numerical outputs for the years 1990, 1995, 2000, 2005, 2010, 2015, 2020, and 2023 are presented in Fig. 6 and Fig. 7 below.
As shown in Fig. 7, reach wise river pattern of flow was assessed under respective years where the value ranges from 1.02 to 1.38 i.e. <1.5 where river found multi-channel determined as braided nature on its planform. A higher sinuosity value observed in the year 2015 which was 1.38 at Reach C (CS-18 to CS-27) and lower value at Reach A was 1.02 respectively. Fig. 6 graph trend looks matching the Fig. 7 graph trend where river sinuosity and mid channel line length observed higher at the same year 2015 i.e., 54.79 km and SI- 1.02 to 1.38.
The pattern of the river channel-braiding index determined from GIS in the form of map as shown in Fig. 8 and Fig 12. The river channel 2015 found highly braided from the GIS mapping and analysis the braided trend started in 2015 where the bifurcation observed higher at CS-13 to CS-16 and CS-18 to CS-21. After the construction of embankment along both side of the river in year 2020, we found the regular channel was brought under control upto CS-16 but channel found evolved beyond CS-18.
4.2 Channel Migration (Mid-line Migration)
West Rapti River (WRR) shifting trend was analyzed using Remote Sensing data and GIS techniques within its planform under 27 cross-sections from east to west over 53 km range. Spatio-temporal maps of river channel generated by overlapping the mid-channel line of the years 1990, 1995, 2000, 2005, 2010, 2015, 2020 and 2023 respectively as shown in Fig. 9. The Fig. 9 shows those cross-sections, MCL overlay, and locations where the channel flows. The river channel shift was observed in 6 directions only whereas the channel flows from East to West.
When the channel shift was evaluated, the year 1990 was used as the base year. We observed that the channel shift was higher in the following years: 1995 under CS-19; 2000 under CS-14 and CS-23; 2005 under CS-11, CS-19, and CS-23; 2010 under CS-23; 2015 not found above 800 m but above 700 m on CS-11 and CS-14; 2020 under CS-23 and CS-25; as well as 2023 under CS-11 and CS-23 separately. In altogether the channel shifting was found highest over CS-11 and CS-23 from 1990 to 2023 i.e., 808.47 m and 1426.97 m respectively. While the channel movement was also recorded in 6 directions where we dig out that within CS-5, CS-8, CS-9, and CS-12 in Southwest direction since 1990 then within CS-11, CS-19, and CS-23 in Northeast direction and within CS-25 and CS-27 in South and Southeast direction respectively as shown in Table 2 and Fig. 11.
From Fig. 10 and table 2, We have figured out that there was highest shift of the river channel at CS-23 in Northeast direction, also at CS-19 channel shift was found in both direction North and South since year 1990.
Table 2: Tabular descriptions of WRR mid-channel line shifting distance and direction from 1990 to 2023 where year 1990 riverline taken as base year where N: North, NE:Northeast, NW: Northwest, SE: Southeast, SW:Southwest, and S:South. (Year_dist: shifting Distance and Year_dir: shifting direction)
Fig. 11 represents the channel pattern at different cross-sections from 1 to 27 where we can observe the shift of river channel from 1990 to 2023. The place where the channel was in 1990 and 2023 can be individually observed as provided from top to bottom CS-1 to CS-27 respectively with the flow, shift distance and direction. In CS-1 the river has shifted 304.98 m southeast whereas at CS-27 the river had shifted 386.73 m Southeast direction.
4.3 Braiding Index Calculation
Braiding pattern of the West Rapti River found decreasing in trend might be the impact of river embankments development by Nepal government. The trend was found decreasing from 1990 to 1995 then increasing until 2010 at last the trend decreased up to year 2023 i.e. 1.64 to 1.03. As provided in Fig. 12 below, maximum mean braiding index (BI) depicts that the channel has highest number of braid bars i.e. 1.76 in year 2010 and minimum mean BI depicts the lowest number i.e. 0.76 in year 1995 respectively throughout the river channel.
4.4 Geospatial determination of driver for channel migration
This study within the river area explains that the channel found braided according to SI and BI values and had undergone successive phases of erosion and deposition processes from 1990 to 2023. This clarifies that the channel had a higher deposition, about 15.35 Sq. km, at these cross-section areas as shown in Fig. 13, that forced river to shift higher and portions of river found colonized by the vegetation within the area.
Maximum historical deposition of sand bars found highly at CS-9 to 13 and CS-18 to 21 where the channel shifting found higher due the expansion and contraction of riverbed bars throughout the channel. Heavy deposition found halt at CS-13 to 16 after year 2020 due to engineering structures like embankment found on both sides of the channel at these cross-sections and others too.
Vegetation: As given in Fig. 14, 2023 google earth image showing vegetation colonized at CS-18 to CS-21 region that is comparable to Fig. 13(C) due to regular sediment deposition and the channel with island bar formed within the region since 2015.
4.5 Discussions
Every Monsoon brings a huge flood within the river channel. The dynamic nature of WRR from 1990 to 2023 within its planform enlightened through the GIS maps as shown above. We found irregular patterns of changes throughout the river channel under respective timeline of 33 years. These changes brought several shifting due to sandbars, point bars, braid bars, mid-channel bars, etc., formation as well as torrent rivers in the south which are the major drivers of west rapti river channel shifting. Confluence points were found to be the main widening part of this river channel that appeared differential patterns i.e., Arjun River, Ransing River and other torrent rivers within West Rapti River. Talchabhadel and Sharma (2014) found in his research that torrent river caused river channel shifting.
The channel shifting was maximum on CS-23 in the year 2023 i.e., 1426.97 m Northeast with respect to base year 1990. Also, within cross-section 25 there was instant shift from 54.39 m to 827.87 m from year 2005 to 2020. In average channel shifting was highly observed in between cross-section 18-21. Whereas CS-2 and CS-24 depict that the channel has been repositioned as it was in the year 1990.
The riverbanks were really conserved in the places where mitigation measures were implemented throughout the channel. Typical features like point bars, cut banks, meanders, anabranching, and braiding were found due to river channel migration. Baniya et. al., (2023) have also found a similar river pattern within Koshi river that shifted about 2 km from 1999 to 2019. The digitised imaginary line from the corresponding year follows the regular channel course, but since 2015, the channel has shifted primarily towards the northeast. We discovered that the trend of channel shifting was increasing, which led to continuous shifting that expanded the river floodplain structure within CS-18 to CS-21. Karnali river channel shifting study by Rakhal et. al., (2021) also found fang development i.e. heavy expansion of the river due to regular erosion and deposition.
With the observation of the field, we identified significant and transformative multichannel developments and movements within cross-sections 18 to 21, as well as in CS-9 to 13. Notably, cross-sections 13 to 16 exhibit a controlled channel behavior, attributed to the strategic implementation of embankments on both sides. This proactive measure, initiated in the year 2020, aims to effectively mitigate the impact of floods on human settlements. The channel's geomorphic characteristics reveal a distinctive meandering behavior, particularly evident from cross-section 22 onward. Noteworthy is the transition of the channel from a multichannel flow with braid bars to forming a single flow. Prior to this transition, the channel displayed a complex multichannel structure, characterized by the presence of braid bars throughout its course.
An additional observation of interest pertains to the colonization of the former riverbed area by vegetation. This phenomenon is indicative of dynamic changes and irregular shifts in the river channel. Such ecological alterations underscore the intricate interplay between geomorphological processes and ecological dynamics in the studied region. Further investigation reveals a discernible southwestward shift in the channel within CS-22 since the year 2000, albeit of minimal magnitude. This directional shift is a noteworthy aspect of the channel's trajectory and contributes to our understanding of its long-term spatial dynamics.
Bifurcation of the river channel first started at CS-14 where an expansion started to begin from year 2000 thus was brought into single channel from 2021. Similar changes in the pattern were found within CS-18 to CS-21 where it started to bifurcate from year 2015 expansion can be visualized in recent map i.e. year 2023.