Spatiotemporal Characterization and Analysis of River Morphology Using Long-Term Landsat Imagery and Stream Power

: Meandering rivers are among the most dynamic Earth-surface systems, which generally appear in fertile valleys, the most valuable lands for agriculture and human settlement. Landsat time series and morphological parameters are complementary tools for exploring river dynamics. Karun River is the most efﬂuent and largest meandering river in Iran, which keeps the Karun’s basin economy, agriculture, and industrial sections alive; hence, investigating morphological changes in this river is essential. The morphological characteristics of Karun have undergone considerable changes over time due to several tectonic, hydrological, hydraulic, and anthropogenic factors. This study has identiﬁed and analyzed morphological changes in Karun River using a time series of Landsat imagery from 1985–2015. On that basis, morphological dynamics, including the river’s active channel width, meander’s neck length, water ﬂow length, sinuosity index, and Cornice central angle, were quantitatively investigated. Additionally, the correlation between the stream power and morphological factors was explored using the data adopted from the hydrometric stations. The results show that the dominant pattern of the Karun River, due to the sinuosity coefﬁcient, is meandering, and the majority of the river falls in the category of developed meander rivers. Moreover, the number of arteries reduced in an anabranch pattern, and the river has been migrating towards the downstream and eastern sides since 1985. This phenomenon disposes a change in the future that can be hazardous to the croplands and demands speciﬁc considerations for catchment management.


Introduction
Karun basin is one of the most critical basins in West Asia and the second most important river in the Persian Gulf and Oman Sea catchments.It is a river with greater discharge and the only navigable river in Iran with a critical role in the Karun basin's ecological, environmental, and socio-economic conditions.The river runs through several big cities such as Ahvaz, Shushtar, and Khorramshahr, with more than five million inhabitants.It is the main water supply for many agro-industries near the city of Ahvaz and for producing over 5.5 million tons/year of sugarcane.Therefore, monitoring the morphology and spatiotemporal dynamics of the Karun River is essential as a prerequisite for effective river management.
Meanders are the most common river pattern in populated areas [1], which are a response of the river to decrease its energy, whereas it leads to an increase in length and reduction in the channel's slope.Meandering rivers time and again occur in fertile valleys, the most valuable lands for agriculture and human settlement.Given that floodplains have been gradually occupied by growing urbanization and industrial development [2], this process increases the importance of flood control, bank erosion control, and meander migration.Therefore, investigating and revealing river path alteration trends are essen-tial for assessing river dynamics, catchment management, risk reduction, and floodplain improvement [3,4].
Investigating the river path trend is feasible by studying the changes in morphological parameters of the river (e.g., width, wavelength, river length, sinuosity, and central angle).There is an extensive body of literature that has explored the linkage between the morphological characteristics and variables such as sediment rate, discharge fluctuation, river bed slope, climate change, gravel mining, dam construction, human interventions, and land use and land cover changes [5][6][7][8][9][10][11][12][13][14][15].Riverbank erosion can lead to landscape degradation as well as environmental and socio-economic impacts, which have been observed in various countries at different scales [16].Bank erosion and displacement of the river's path can result in river instability and meandering formation and development.In addition, human intervention in natural conditions and the riparian zone intensifies bank erosion [17][18][19][20][21].As a result of human intervention, dam construction affects the downstream, which can seriously change the hydro-morphological characteristics such as sediment transport pattern and flow discharge, which cause riverbed run narrowing and bank erosion [22][23][24][25][26][27].
Stream power can be used to investigate the influence of flow discharge changes on river morphology [28].Stream power, the product of discharge and the river's slope, is a helpful tool for describing the river's morphological changes [29,30].Stream power has been less explored in the literature and was mainly considered concerning the bedload distribution [31], bank strength, prediction of channel dynamic [32], and sinuosity [33].
The bodies of water and river planform can be extracted through satellite images and using relevant spectral indices.Time series of satellite images paved the way for depicting the channel layout over time driven by natural events such as floods or anthropogenic activities [28].Freely available satellite imagery with regular global coverage is a valuable resource for river monitoring.Although satellite imagery has a long-term archive, the spatial resolution of images has been increased over time, which is an essential factor, especially for monitoring narrow rivers.In many river studies, satellite images have been used for evaluating and monitoring river morphology [34][35][36][37][38][39][40][41][42][43].
Given the importance of the Karun River above, it is necessary to monitor the river's dynamics, for example, by examining and quantifying the geometry and pattern of the meanders over time.This paper investigates the trend of morphological changes of the Karun River and detects the underlying drivers over the last 30 years using Landsat imagery.In previous studies, the number of morphological characteristics was limited [44]; in this study, the sinuosity coefficient, central angle (Cornice angle), active channel width, wavelength, and river length, as well as arc properties, are analyzed.In addition, the river reach is categorized with three different methods (Schumm, Cornice, and Leopold), which have not been investigated previously.Finally, the relation between unit stream power (USP) and morphological parameters, essential for basin monitoring and management, will be quantified.These objectives, in conjunction, will help to understand how the changes in river flow discharge affect the river's morphological changes.

Study Area and Data
Karun River is located in the Karun basin with a catchment area of 67,112 km 2 , which is characterized as the highest discharge and the longest river in Iran.The catchment is divided into three sub-basins and 42 plains (Table 1).With a total length of 950 km, Karun originates mainly from the Zardkuh Mountains of the Bakhtiari district and receives various tributaries such as Dez and the Kuhrang before passing through the city of Ahvaz and discharging into Arvandrud.In this study, a portion of the river with a length of 300 km (Figure 1) is investigated.Its upstream part is located close to one of the largest dams in Iran, Gotvand Dam, built in 2003, and the downstream region is located near the city of Shirinshahr.Several massive projects have been constructed since 70 years ago that show the socio-economic importance of this river [45].Significant changes were imposed on the river's morphology during these constructions; the most important outcomes were changing the river width and meanders.in Iran, Gotvand Dam, built in 2003, and the downstream region is located near the city of Shirinshahr.Several massive projects have been constructed since 70 years ago that show the socio-economic importance of this river [45].Significant changes were imposed on the river's morphology during these constructions; the most important outcomes were changing the river width and meanders.

Data
Multi-temporal Landsat images were used to detect the dynamic of Karun River in three time points, namely, 1985, 2000, and 2015.To account for the topography of the region, a digital elevation model (DEM) with a resolution of one arc-second (30 m) was acquired from the SRTM satellite (Table 2, Figure 2-left panel).A land use map was extracted from Landsat-8 in 2015 (Figure 2-right panel).Since the Karun River reaches its highest discharge in spring, the considered satellite data belonged to May.Three scenes of Landsat data have covered the study area.

Data
Multi-temporal Landsat images were used to detect the dynamic of Karun River in three time points, namely, 1985, 2000, and 2015.To account for the topography of the region, a digital elevation model (DEM) with a resolution of one arc-second (30 m) was acquired from the SRTM satellite (Table 2, Figure 2-left panel).A land use map was extracted from Landsat-8 in 2015 (Figure 2-right panel).Since the Karun River reaches its highest discharge in spring, the considered satellite data belonged to May.Three scenes of Landsat data have covered the study area.

Methods
Landsat data with 15-year intervals  were processed to delineate the river's bed extent.The modified normalized difference water index (MNDWI) was employed to reveal the water body extents.MNDWI is a powerful remote sensing index for water body extraction [46].MNDWI is computed using the near infrared and the short-wave infrared bands.Equation (1) shows the relation between NIR and SWIR for the calculation of MNDWI.MNDWI = (Green − SWIR)/(Green + SWIR) (1) A manual improvement followed the river's bed delineation; then, the centerline was extracted using the centerline function in ArcGIS.For a detailed exploration of the river, the study area was divided into eight segments based on Schumm's study in 1985, which categorizes the river's pattern according to three criteria: the number of courses, sinuosity coefficient (SC), and lateral stability (Figure 3).

Methods
Landsat data with 15-year intervals  were processed to delineate the river's bed extent.The modified normalized difference water index (MNDWI) was employed to reveal the water body extents.MNDWI is a powerful remote sensing index for water body extraction [46].MNDWI is computed using the near infrared and the shortwave infrared bands.Equation (1) shows the relation between NIR and SWIR for the calculation of MNDWI.MNDWI = (Green − SWIR)/(Green + SWIR) (1) A manual improvement followed the river's bed delineation; then, the centerline was extracted using the centerline function in ArcGIS.For a detailed exploration of the river, the study area was divided into eight segments based on Schumm's study in 1985, which categorizes the river's pattern according to three criteria: the number of courses, sinuosity coefficient (SC), and lateral stability (Figure 3).The SC and central angles are used to characterize meandering rivers [33,34,47,48].The value of SC ranges between 1.5 and 4 [49][50][51].Channels with an SC value of less than 1.05, from 1.06-1.3,and greater than 1.3 are called straight, sinuous, and meandering rivers [52].The central angle was proposed by Cornice (1985).It works according to the central angle to quantify the extent of meandering development and distinguishing them in The SC and central angles are used to characterize meandering rivers [33,34,47,48].The value of SC ranges between 1.5 and 4 [49][50][51].Channels with an SC value of less than 1.05, from 1.06-1.3,and greater than 1.3 are called straight, sinuous, and meandering rivers [52].The central angle was proposed by Cornice (1985).It works according to the central angle to quantify the extent of meandering development and distinguishing them in alluvial rivers (1980).Using the river route, morphological parameters were extracted for each time point as below: Width: To obtain the river's width for each arc, four or five sections, including the beginning, end, and middle of the arc with different widths, were measured, and their average was introduced as a representative for each arc (Equation ( 2)).
Wavelength: For this purpose, inflection points (i) for each arc were determined, and the direct distance between each pair of points was measured as half of the wavelength because each wave consists of two arcs in the form of a sine function.
Central angle: The Cornice central angle indicates the river's form development, which can be deduced by enclosing a circle in an arc.In doing so, a circle is surrounded by an arc, then the first and endpoints of their connections are adjoined to the center of the circle; the angle created in the center is called the Cornice central angle.
Sinuosity: Sinuosity or the coefficient of curvature for each arc is the ratio of the river's length to half of the wavelength.For each arc, the length of the river's centerline was measured between the inflection points and half of the wavelength, which is the direct length between these points (Equation ( 3)).
Slope: A digital elevation map (Figure 2-left panel) was used to calculate the slope for each arc.The points closest to the river (right and left bank of the arc) were identified, and their slope was calculated.According to Figure 4, the slope between three points (on each side of the arc) was calculated (Equations ( 4) and ( 5), and after averaging for each side, one slope was obtained, and an equivalent slope (S) was extracted from the mean of the two values obtained for each side (Equation ( 6)).
The hydraulic geometric flow discharge relation for the studied sections was obtained using the historical data.The relations between river discharge (Q) and some geometric (H, B, and A) and hydraulic parameters (V and τ) are shown in Figure 5. Flow parameters including H (average height), A (bed surface area), V (flow velocity), and τ (bed shear stress) were calculated based on observation in Ahvaz hydrometric station.Then, stream power per unit length, width, area, and weight were calculated using Equations ( 7)- (10).
where γ is water-specific weight, Q denotes the river flow discharge, and S shows the longitudinal river's bed slope obtained from the DEM layer.The hydraulic geometric flow discharge relation for the studied sections was obtained using the historical data.The relations between river discharge (Q) and some geometric (H, B, and A) and hydraulic parameters (V and τ) are shown in Figure 5. Flow parameters including H (average height), A (bed surface area), V (flow velocity), and τ (bed shear stress) were calculated based on observation in Ahvaz hydrometric station.Then, stream power per unit length, width, area, and weight were calculated using Equations ( 7)- (10).Morphological parameters and stream power have been calculated cumulatively across the study area.Since the integral for each parameter was converted to sigma, the distance between the arcs was assumed to be equal.As shown in Equations ( 11) and ( 12), X is the cumulative number of parameters that include the USP, morphological, hydraulic, and geometric parameters, x is the number of parameters in each arc, ds is equal, and n is the number of arcs along the river (Figure 4).
We use the normalized mass method to investigate the relationship between the UPS and morphological parameters of the river.Their graphs were drawn, the trend line with the best regression coefficient was plotted on each curve, and its equation was determined (1985, 2000, and 2015).Finally, the general relation was obtained.In the last step, the centerline of the river for each time point (1985, 2000, and 2015) was superimposed in AutoCAD, and common arcs were examined to determine the arc's migrations.Displacement and its clockwise angle relative to the north were calculated for each time interval (every 15 years and at the end of the last interval).Finally, the changes in the river's routes and the arcs with the highest migration were identified and analyzed.

Schumm Approach for Channel Pattern
According to Schumm's channel classification scheme (1985), the Karun River demonstrated three primary forms: straight, meandering, and braided in the eight identified intervals.Table 3 summarizes the river's morphologic properties, whereas Tables 4 and 5 provide the temporal changing meandering and braiding development along the river path.As shown in Table 4, rotation and conversion are the dominant properties in the meandering segments, where segment 4 has the most active part.On the other hand, sinuous side channels, split channels, sub-parallel anabranches, and cutoff loops were the main properties of the braided patterns (Table 5).

Sinuosity Coefficient
Figure 6 demonstrates the changes in arc sinuosity along the river from Gotvand to Shirinshahr.As shown, after 55 km from the beginning of the reach, the river altitude decreases to 40 m, while the altitude variation within 300 km of the river is only 60 m, so the river proceeds with a gentle slope at the beginning of the reach.Based on Figures 7  and 2, the conclusion can be drawn that the slope has an inverse relation with the SC.In the first 55 km, the highest SC was 3.18 (1985), 1.38 (2000), and 1.28 (2015), while throughout the reach (within 150-160 km from the beginning), it was 4.67 with an altitude of about 20 m in 2015.As seen in Table 6, the highest SCs in 1985 and 2000 were 3.44 and 3.99, respectively, which were identified within the range of 205 to 212 km from the beginning at 11 m.

Sinuosity Coefficient
Figure 6 demonstrates the changes in arc sinuosity along the river from Gotvand to Shirinshahr.As shown, after 55 km from the beginning of the reach, the river altitude decreases to 40 m, while the altitude variation within 300 km of the river is only 60 m, so the river proceeds with a gentle slope at the beginning of the reach.Based on Figures 2 and 7, the conclusion can be drawn that the slope has an inverse relation with the SC.In the first 55 km, the highest SC was 3.18 (1985), 1.38 (2000), and 1.28 (2015), while throughout the reach (within 150-160 km from the beginning), it was 4.67 with an altitude of about 20 m in 2015.As seen in Table 6, the highest SCs in 1985 and 2000 were 3.44 and 3.99, respectively, which were identified within the range of 205 to 212 km from the beginning at 11 m. and 2, the conclusion can be drawn that the slope has an inverse relation with the SC.In the first 55 km, the highest SC was 3.18 (1985), 1.38 (2000), and 1.28 (2015), while throughout the reach (within 150-160 km from the beginning), it was 4.67 with an altitude of about 20 m in 2015.As seen in Table 6, the highest SCs in 1985 and 2000 were 3.44 and 3.99, respectively, which were identified within the range of 205 to 212 km from the beginning at 11 m.As shown in Table 7, the meandering pattern appeared as the dominant morphological pattern at each time interval, whereas the arcs with the straight SC disappeared gradually.Additionally, the sinusoidal was and the meandering pattern was reduced to less than 2%.Although the upper limit of SC equals approximately three, the maximum SC of Karun (Table 6) in 1985, 2000, and 2015 was 3.44, 3.99, and 4.67, respectively.Our results are consistent with other studies such as Hoseynzade and Ismaili (2016) on Babol (3.62) and Talar River (3.04) in Iran and Kiss (2008) on Tisza River (more than 4) in Hungary.According to Figure A1, the SC increases as the altitude decreases.In addition, the river's migration capacity was minimized at high altitudes due to the mountainous texture and steep slopes; furthermore, as the river moves toward lower altitudes, the land's texture changes from the mountainous to the plain form and river migration increases.

Cornice Central Angle
Cornice's central angle was calculated for all arcs (Figure A2 (right panel) in Appendix A), and the percentage of meandering arcs of the Karun River was calculated for the studied time points (Table 8).As shown in Table 8, the central angles in the study area most frequently occurred within 85-158 degrees, related to the developed meandering category (70%, 73.5%, and 75% of arcs for 1985, 2000, and 2015, respectively).The second frequency was related to the non-developed meandering in all time intervals.Additionally, the lowest frequency in 1985 was related to 0-41 degrees (similar to meandering), which falls in the more developed meandering category (158-296) for 2000 and 2015.Additionally, according to Table 8, there is no case (central angle of arc) with more than 296 degrees and no more arcs from 2000 to 2015 in 0-41 degrees.

River Width
The results show that the river width reduced by 32% during the studied periods (Figure A2 (left panel) in Appendix A).The overall trend of the river's width showed widening towards the lower altitudes at the river's downstream, where the landscape changes from high altitudes to plain.

Stream Power
USP (Ω) was calculated for each arc along the river (Equation ( 7)), and its relation with the morphological factors (η i ) was explored.K i and R 2 are coefficients of USP and determination (Equation ( 13)), respectively, presented in Table 9 and the last column is the general coefficient for each morphological parameter.The stream power per unit length has the most prediction ability for the morphological parameters.Since all the linear equations show a more than 88 percent correlation, it could be concluded that these equations are suitable for understanding the dynamic of morphological parameters using the power stream.

River Migration
The following results were obtained by comparing the river pattern in the study years.8 shows the spatiotemporal morphological evolution of the Karun River over space and time.The red parts are the dried and disappeared sections; the green color shows the emerging sections, and the blue part depicts the unaltered locations.Our results show that the river migrates gradually toward the east and southeast (Figure A3 in Appendix A), indicating that the river's east bank is unstable.Figure 8 shows the arcs with the most displacement along the Karun River reach (Gotvand to Shirinshahr).These changes indicate that the riverbanks in these arcs are unstable, so they need to be considered more than other parts.According to Figure 2 (right panel), agriculture and the built-up areas are the dominant land use.Therefore, a sudden event such as heavy rain or flood could have many socio-economic consequences in the floodplain.

Influence of Gotvand Dam on the Karun River Morphological Pattern
To investigate the influence of Gotvand Dam on the morphological changes in Karun River, the river path was examined in 1985, 2000, and 2015 and the common arcs (including 57 arcs for each year) were extracted (Table 10).Our exploration shows that the width of Karun River has decreased after the dam construction, which was expected due to the reduction and changes in the water and sediment discharge.However, sinuosity and Cornice central angle of more than half of the arcs in both periods have shown a positive change.For this specific case study, the percentage and type of changes for the curvature coefficient were almost the same both pre-and post-dam construction, whereas an increasing trend was observed for the changes in central angle of a positive type.As a result, the decrease in the river sediment and the increase in the river's channel erosion are consistent with the changes in sinuosity and Cornice central angle.

Influence of Gotvand Dam on the Karun River Morphological Pattern
To investigate the influence of Gotvand Dam on the morphological changes in Karun River, the river path was examined in 1985, 2000, and 2015 and the common arcs (including 57 arcs for each year) were extracted (Table 10).Our exploration shows that the width of Karun River has decreased after the dam construction, which was expected due to the reduction and changes in the water and sediment discharge.However, sinuosity and Cornice central angle of more than half of the arcs in both periods have shown a positive change.For this specific case study, the percentage and type of changes for the curvature coefficient were almost the same both pre-and post-dam construction, whereas an increasing trend was observed for the changes in central angle of a positive type.As a result, the decrease in the river sediment and the increase in the river's channel erosion are consistent with the changes in sinuosity and Cornice central angle.This study investigated long-term morphological changes of the Karun River from downstream of Gotvand Dam to Shirinshahr (~300 Km).Gotvand Dam is one of the largest dams in Iran and last downstream dam on the Karun River, which has led to changes in water and sediment volumes [44].The average daily discharge of the river was 534 m 3 s −1 whereas, after the dam construction, it decreased to 490 m 3 s −1 .The average daily suspended sediment loads before and after the dam construction were 16.36 and 9.12 gL −1 , respectively.Based on Equation ( 14), the river width directly relates to flow discharge and is inversely related to the sediment discharge.Therefore, decreasing the flow discharge of the Karun River can be a reason for the shrinkage in the river width, which is similar to the findings of Mengen et al. (2020) [54] on the Mekong River and Byishimo et al. (2014) [55].
Morphological changes happen gradually in uncontrolled rivers, but in rivers under human control, these changes occur in a shorter time.Dam construction, land use change, water transmission networks, urban development, sand mining, and bank protection are influential factors that accelerate river changes [56].Therefore, humans play a crucial role in the river system, hence considering that they are an effective measure for a better realization of human-river interaction.

Morphological Pattern
According to Schumm's approach, the Karun River's reaches were divided into eight intervals; four sections have shown the meandering form, two of them belonged to the braided category, and the other two had straight form.The majority of the morphological changes happened in the first (braided pattern), third (braided pattern), and fourth intervals (meandering pattern), where cropland and residential area are the dominant land use in the floodplain (Figure 2-right panel) and river reach.Therefore, the sudden morphological changes in the river caused by extreme floods will affect human lives and property.On the other hand, the pattern of morphological changes indicated that the slope was inversely related to the river sinuosity, which conforms to the Timár flume experiments [57] and Petrovzszki and Frasson's research along natural rivers at any discharge value [48,58].
The number of river arteries decreased over time, where the flow and sediment discharge reduced after the construction of Gotvand Dam.Furthermore, adjacent arcs were interconnected over time and evolved into a single arc.Shrinkage of the river width is another observation within the 30 years from 1985 to 2015.The studies of Nelson et al. (2013) [59] on the Jackson Dam and its morphological effects on the Snake River have shown that some of the river's arteries were dried up by decreasing flow and sediment discharge.Most of the arc's movements occurred in Gotvand and Mollasani, indicating the greater slope and remarkable migration in the reach.It shows that the river tends to deplete its high energy through transverse displacement floodplain.Three ranges are highlighted in Figure 8, which had the most changes in the river reach during the study period .They are also disposed to change in the future, which can be hazardous to cropland and residential areas.
The river's arrival into the plain provides a condition for morphological changes.Reduction in the slope and river velocity and changes in the texture of river reach are some potential reasons for the creation of the meandering form in the floodplain.Sinuosity coefficient, the river width, and central angle had an inverse relation with the slope, and by gentling the slope, their value increases, which is approved by Frasson et al [48].
Morphological parameters were calculated based on USP, making it possible to predict changes' feasibility before natural events such as floods.In this research, the prediction of morphological parameter changes based on USP has been studied and evaluated for the first time, which had an acceptable result.Human activities in the river reach and floodplains lead to changes in river morphology.When there is a good correlation between USP and morphological parameters, it indicates that the river eventually becomes a USP function.

Conclusions
This paper explored the long-term trend of the Karun River as the most important river of Iran and revealed that it evolved into a developed meandering system.Due to the decreasing amount of water discharge, some arteries were lost, the width of the river has decreased, and it turned into a single stream in line with the majority of previous studies.If the current trend continues, cropland and urban areas on the fragile parts of the river will be at high risk.In addition, since the river migrates toward the east and southeast, the east coastal lands of the river will be more susceptible to change.Additionally, the relevance of morphology parameters and stream power has been surveyed in this study, which helps to anticipate the possible alteration of morphological characteristics.
We suggest considering the influence of hydrological events such as extreme floods in recent years and other relevant variables such as soil texture and the riparian zone on the Karun River in future studies.Moreover, in the coming years, Sentinel-2 images can be considered as an alternative source due to their higher spatial resolution, which shows more details of spatial variations in river dynamics.

Figure 1 .
Figure 1.Study area in the Karun River basin.(a) Iran's second-grade catchment, (b) Great Karun catchment, (c) a river with a 3 km buffer zone (flow direction is north to south).Figure 1. Study area in the Karun River basin.(a) Iran's second-grade catchment, (b) Great Karun catchment, (c) a river with a 3 km buffer zone (flow direction is north to south).

Figure 1 .
Figure 1.Study area in the Karun River basin.(a) Iran's second-grade catchment, (b) Great Karun catchment, (c) a river with a 3 km buffer zone (flow direction is north to south).Figure 1. Study area in the Karun River basin.(a) Iran's second-grade catchment, (b) Great Karun catchment, (c) a river with a 3 km buffer zone (flow direction is north to south).

Figure 2 .
Figure 2. (Left panel) Digital elevation map of Karun Basin (SRTM).The 55 km segment is mentioned in the results section.(Right panel) Land use map of the study area extracted from Landsat data in 2015.

Water 2022 , 21 Figure 2 .
Figure 2. (Left panel) Digital elevation map of Karun Basin (SRTM).The 55 km segment is mentioned in the results section.(Right panel) Land use map of the study area extracted from Landsat data in 2015.

Figure 3 .
Figure 3. Segmentation of the study area for meandering characterization.

Figure 3 .
Figure 3. Segmentation of the study area for meandering characterization.

Figure 4 .
Figure 4. Calculation of the morphological parameters from satellite images.

Table 4 .Table 5 .Table 5 .
Observed changes are shown by (*) from 1985-2015 in the meandering pattern along Karun River based on Schumm's studies [53].Observed changes are shown by (*) in the braided pattern along the Karun River in 30 years (1985-2015) in the braided pattern based on Schumm's studies [53].Observed changes are shown by (*) in the braided pattern along the Karun River in 30 years (1985-2015) in the braided pattern based on Schumm's studies [53].
the study reach of river in three years(1985-2000-2015)

Water 2022 , 21 Figure 7 .
Figure 7. Plan morphological changes in Karun over 30 years: (a) Changes in the study area during three intervals shown at the bottom.(b) Zoomed-in view of 9 parts that are marked in (a).

Figure 7 .
Figure 7. Plan morphological changes in Karun over 30 years: (a) Changes in the study area during three intervals shown at the bottom.(b) Zoomed-in view of 9 parts that are marked in (a).

Figure 7
illustrates the trend of the river's dynamic.Part (a) shows the river changes from 1985-2000 and 2000-2015 and the changes that occurred from 1985-2015.Part (b) represents the trend of changes from 1985-2015.

Figure A2 .
Figure A2.(Right) Cornice central angle graph of the study region in 1985, 2000, and 2015, and (Left) the process of changing the river width in the study region (1985, 2000, and 2015).

Figure A2 .
Figure A2.(Right) Cornice central angle graph of the study region in 1985, 2000, and 2015, and (Left) the process of changing the river width in the study region (1985, 2000, and 2015).

Table 1 .
Details of Great Karun and its three sub-basins (Iran's Ministry of Energy, 2012).

Table 1 .
Details of Great Karun and its three sub-basins (Iran's Ministry of Energy, 2012).

Table 2 .
DEM and satellite data * to monitor the dynamic of Karun River.

Table 2 .
DEM and satellite data * to monitor the dynamic of Karun River.

Table 3 .
Properties of river pattern and changes in each interval.

Table 4 .
Observed changes are shown by (*) from 1985-2015 in the meandering pattern along Karun River based on Schumm's studies [53].

Table 8 .
Rate of the meandering pattern in the Karun River based on the Cornice index in 1985, 2000, and 2015.

Table 9 .
Relationship between morphological parameters and unit stream power (k i ) and coefficient determination (r 2 ) in 1985, 2000, and 2015.

Table 10 .
The percentage (frequency) and type of morphological changes pre-and post-Gotvand Dam construction.