Predicting Channel Conveyance and Characterizing Planform Using River Bathymetry via Satellite Image Compilation (RiBaSIC) Algorithm for DEM-Based Hydrodynamic Modeling

: Digital Elevation Models (DEMs) are widely used as a proxy for bathymetric data and several studies have attempted to improve DEM accuracy for hydrodynamic (HD) modeling. Most of these studies attempted to quantitatively improve estimates of channel conveyance (assuming a non-braided morphology) rather than accounting for the actual channel planform. Accurate representation of river conveyance and planform in a DEM is critical to HD modeling and can be achieved with a combination of remote sensing (e.g., satellite image) and field data, such as water surface elevation (WSE). Therefore, the objectives of this study are (i) to develop an algorithm for predicting channel conveyance and characterizing planform via satellite images and in situ WSE and (ii) to estimate discharge using the predicted conveyance via an HD model. The algorithm is named River Bathymetry via Satellite Image Compilation (RiBaSIC) and uses Landsat satellite imagery, Shuttle Radar Topography Mission (SRTM) DEM, Multi-Error-Removed Improved-Terrain (MERIT) DEM, and observed WSE. The algorithm is tested on four study areas along the Willamette River, Kushiyara River, Jamuna River, and Solimoes River. Channel slope and predicted hydraulic radius are subsequently estimated for approximating Manning’s roughness factor. Two-dimensional HD models using DEMs modified by the RiBaSIC algorithm and corresponding Manning’s roughness factors are employed for discharge estimation. The proposed algorithm can represent river planform and conveyance in single-channeled, meandering, wandering, and braided river reaches. Additionally, the HD models estimated discharge within 14%–19% relative root mean squared error (RRMSE) in simulation of five years period.


Introduction
Digital Elevation Models (DEMs) are commonly used for representing terrain in many geospatial, hydrologic, and hydraulic applications. Some of the commonly used DEM datasets have global coverage, including the Shuttle Radar Topography Mission (SRTM), Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), Advanced Land Observing Satellite Global Digital Surface Model (AW3D30), Multi-Error-Removed Improved-Terrain (MERIT), TerraSAR-X add-on for Digital Elevation Measurement (TanDEM-X), and Hydrological data and maps based on SHuttle Elevation Derivatives at multiple Scales (HydroSHEDS). These global DEMs are often the best available terrain data in many data-poor regions. Global DEMs are usually produced via remote sensing techniques and/or by merging multiple datasets. Thus, global DEMs are prone to errors in localized applications [1]. Studies have shown that fluvial features (e.g., river bathymetry, planform, and floodplain) are often not accurately represented in DEMs [2,3] and, thus, can introduce significant errors into hydrodynamic (HD) modeling. The errors may originate from the data source, data acquisition, signal processing algorithm, and/or spatial and temporal resolution of the DEM used in HD modeling [4][5][6][7][8][9].
Inaccurate representation of channel conveyance in a DEM can be averted by merging surveyed bathymetric data. However, a bathymetric survey can be both time and resource-intensive [10]. Therefore, DEM-based HD modeling in a data-poor region is often undertaken without modifying the DEM hence the changes in river bathymetry and planform are not accurately represented. Nevertheless, several studies have investigated on improving DEM accuracy over floodplains for hydrodynamic modeling [11,12]. For instance, MERIT DEM was introduced by minimizing elevation errors [13]. However, a unified DEM correction mechanism functional for both floodplain and channel is rare. It is particularly important in dynamic river systems where the river bathymetry, planform, floodplain connection, etc. are continuously evolving. This mutual interaction is eventually linked to the flow regime [14]. Therefore, it is important to search for a solution that can incorporate the three-dimensional configuration of a river into a DEM for an accurate representation of conveyance and planform into HD modeling. Such a solution would be especially useful for data-poor regions if it can perform well, despite limited field data. Certain gaps in field data can be addressed by remotely acquired data such as satellite imagery, which provides near real-time observation of a large part of the Earth and its waterbodies.
Satellite and aerial images have long been used to study river hydraulics [8,15]. Landsat-based river widths over long period have been correlated with altimetry-based water surface elevation (WSE) for producing equivalent cross-sections for an ungauged watershed [16]. However, this approach cannot predict riverbed elevation below the low water level and is limited to trapezoidal channel assumption. Their approach would be especially useful if depth prediction below the low water level could be made. A curve-fitting method is also developed to correlate remotely-sensed WSE and river width to predict river bathymetry [17]. However, this approach is only applicable to single-channeled rivers with predefined cross-sections. Several other methods based on open channel flow equations have been proposed to predict the depth of water (below low water level or a given stage) to improve remote sensing-derived DEMs [2,6,18,19]. These methods are commonly based on conceptual channel shapes (e.g., rectangular, prismatic, parabolic) for deriving equivalent cross-sections using known reach-averaged hydraulic properties (e.g., rating curve, channel slope, roughness factors, etc.). Two studies used conceptual cross-sections to calibrate channel dimension and Manning's roughness factors simultaneously and found these two parameters can be interdependent [20,21]. Similarly, Manning's roughness factor is found to be sensitive to DEM resolution for a DEM raster-based (i.e., DEM) steady HD model [22]. The dependency arises from the fact that a coarse DEM may show a higher river bed elevation due to the applied data interpolation method. Such simplification of river dimensions in large and morphologically complex rivers can lead to misrepresentation in spatial flood analysis. Additionally, DEMs only represent a snapshot in time, therefore, it becomes outdated if the river planform is changed. Nevertheless, remote sensing data such as time series of satellite imagery has the capability of capturing spatial and temporal dynamics of large and complex rivers.
A few remote sensing-based methods have been developed to characterize river dimensions [23][24][25][26][27]. Some of these methods utilize passive optical imaging to map river depths. The fundamental principle in retrieving river depth from a passive optical image is that depth in shallow clear-flowing water is linearly (or monotonically) correlated to an image-derived quantity, which is a function of reflectance [28]. These methods have shown promising results in clear-flowing water, but their performance over turbid water and unsurveyed reaches are limited. Two methods named Hydraulically Assisted Bathymetry (HAB), and Flow Resistance Equation-Based Imaging of River Depths (FREEBIRD) are capable of mapping bathymetry in the absence of field-based depth measurements [24,29]. However, they require an initial estimate of the channel aspect (width-to-depth) ratio or knowledge of discharge [30]. Another method using Image-to-Depth Quantile Transformation (IDQT) was proposed for inferring river bathymetry using images [30]. This technique can be applied in the absence of depth data for depth-to-image quantity calibration and does not necessarily require a linear correlation. IDQT can replace field surveyed depth measurements with a probabilistic distribution of the ratio between spot depths and reach-averaged depth. However, obtaining reach-averaged depth might be challenging in a river where no prior knowledge of the river bathymetry is available. Another method was presented based on entropy theory in which Moderate Resolution Imaging Spectroradiometer (MODIS) data and a predefined horizontal velocity distribution along the channel cross-section were used to predict bathymetry [31]. This study noted the method's applicability in narrow rivers is uncertain due to the spatial resolution of MODIS data (250 m). Additionally, the elliptical velocity distribution approximation might not be valid for braided or meandering river sections where maximum velocity is dictated by river morphology. Therefore, this approach is only applicable to selected river sections.
To summarize, existing indirect depth measurement techniques are either dependent upon prior knowledge and/or predefined cross-section shapes. A study based on an assessment of different bathymetry prediction methods noted that the predefined river cross-section assumption might be acceptable for one-dimensional (1D) modeling but its usefulness for two-dimensional (2D) models is uncertain [32]. This study also recommended that a more robust and spatially accurate bathymetry prediction method is needed for large rivers and streams. The existing bathymetry prediction methods usually approximate river cross-sections at defined transects and interpolate in between those transects. This is a critical obstacle to spatial flood analysis via 2D modeling because river width and planform is not the same along the river length. Additionally, reliance on the dimension, discharge, and/or velocity data make them inapplicable in data-scarce areas. Therefore, being able to utilize minimal hydrologic data for prediction of river bathymetry and hence conveyance could be beneficial. Linking WSE with earth-observing satellite images can be particularly useful because ground-based WSE measurements are more frequently available than flow measurements and current satellite imaging missions scan a large part of the earth at a convenient time interval (e.g., approximately 16 days for Landsat). Additionally, space-based WSE measurement is already being performed from current missions such as Jason-3, Satellite with ARgos and ALtiKa (SARAL), and Ice, Cloud and land Elevation Satellite-2 (ICWSat-2) [33]. Future missions, such as National Aeronautics and Space Administration's (NASA) Surface Water and Ocean Topography (SWOT) are also scheduled to estimate WSE, water mask (water extent), and slope [34] of rivers wider than 100 m. Thus, a WSE and satellite image-based conveyance-prediction algorithm can potentially be used in remote, inaccessible, and truly data-poor regions in the future. Hence, the objectives of this study are: 1. to develop an algorithm for predicting channel conveyance and characterizing planform via satellite images and observed (in situ) WSE as input; 2. to estimate river discharge using the predicted conveyance via an HD model after applying the algorithm.
In the following sections, a novel algorithm named the River Bathymetry via Satellite Image Compilation (RiBaSIC) is proposed and validated. The purpose of this research is conveyance prediction and planform characterization of a river reach using publicly available remote sensing data and observed WSE that would be suitable for 2D HD modeling. In this paper, the terms "observed" refers to any data that is measured in the field or surveyed or processed using measured data (e.g., streamflow, rating curve).

Materials and Method
The methodology of this research follows a two-step approach that includes: (i) developing a river conveyance prediction and planform characterization algorithm via satellite images and observed WSE (ii) estimating river discharge using a moderate resolution (i.e., 30 m) global DEM-based HD model after applying the conveyance-prediction algorithm. The input data needed for this study are satellite images, observed WSE, and DEM. A brief description of the methodology is provided in the following subsections.

Developing a River Conveyance Prediction and Planform Characterization Algorithm via Satellite Image and Observed WSE
Satellite images are useful in analyzing the change of river extent over time. An algorithm to predict river conveyance via linking satellite images and observed WSE is presented in the following paragraphs. A schematic of the method is shown in Figure 1, a qualitative example is provided in Appendix A, and a list of acronyms is presented in Abbreviations. The only observed data needed for proposed RiBaSIC algorithm is time series of WSE at one location near the study area; thus, it can be applied to a region where hydrologic (i.e., flow, rating curve) and hydraulic data (e.g., roughness factors, channel aspect, conceptual channel shape, etc.) are limited. However, the length of the study area should be limited such that it does not include any major inflow or outflow and seasonal WSE variation is uniform within the reach length. To execute the algorithm, a programming language (MATLAB) and a geospatial analysis tool (ArcGIS) is used. The study area for implementing the RiBaSIC algorithm is hereafter called "RiBaSIC extent". The assumptions for executing the RiBaSIC algorithm is also listed below: 1. the WSE of the study area is not affected by the tide; 2. the seasonal variation of WSE within the study area is uniform; 3. the river reach maintains flow continuity meaning no major lateral inflow or outflow; 4. the reach for which the channel profile is determined does not have any abrupt changes in elevation such as a cascade or rapid; 5. the river reach is free from any natural and/or artificial storage areas such as lakes; 6. no significant natural and/or artificial changes happened in the river planform or conveyance within the study period; 7. the width of the actively flowing channel increases with increasing WSE; 8. the satellite image has sufficient bands to cover the visible and Infra-Red (IR) spectrum; 9. the river bankfull width is at least three times the spatial resolution of the satellite image for identifying the channel and side slopes; 10. there are a sufficient number of satellite images available to capture a minimum of one hydrologic cycle.

Time Window Selection
The first task of the RiBaSIC algorithm is to select a time window within which planform shifting (i.e., change in riverbank lines, islands, confluence, etc.) in the RiBaSIC extent is nominal. This can be done by comparing satellite images of the RiBaSIC extent immediately before and after the selected study period. Additionally, the time window should not include any major natural (such as avulsion, cutoff, meander migration) or manmade (such as dredging, channelization) alteration of the river conveyance. However, it is recommended that the time window should not be less than one year so that the seasonal variation of the river stage is captured.

Input Data Processing
The second task for the algorithm is processing the input data for implementing in RiBaSIC algorithm. Available satellite images within the time window for the RiBaSIC extent is collected. For ease of processing, it is recommended that all the collected images have the same source and same spatial resolution. Historic WSE for a nearby streamflow gauge is also collected from which the lowest water level (LWL) is calculated. It should be noted that low water level is the minimum WSE for any given period (e.g., one hydrologic year) but LWL is the minimum WSE from the entire record of historic WSE.
An approximate river polygon (covering the main channel, overbank areas, and adjacent floodplains) for RiBaSIC extent is delineated and converted to a raster with the same spatial resolution and cell alignment as the satellite image. The raster is then converted to an Environmental Systems Research Institute (ESRI) point shapefile (hereafter called as "river points") where each point represents the center of a corresponding raster cell and has a unique identification number (hereafter called as "cell id").
A DEM for the RiBaSIC extent is also collected. This is used to determine a longitudinal channel profile of the river. The profile is estimated by fitting DEM elevations and distance along the center of the main channel of the river. The second-order polynomial channel profile is used because a river profile may not be linear in longer (for example in watershed-scale) reach lengths [35]. The elevation (z) at a distance x is expressed as Equation (1) by this channel profile. The coefficients a, b, and c are derived by fitting DEM elevations and distance along the main channel of the river. Differentiation of the channel profile provides Equation (2) that represents the tangential or linear slope of the channel profile at distance x.

Channel profile
Linear slope: where, x = distance

Shoreline Delineation
The third task is to determine the edges of the water extent at different WSE following an approach similar to the waterline method [12,[36][37][38][39][40]. To determine the edges of the water, all the images within the time window are processed using the Equation (3) to calculate a Normalized Difference Water Index (NDWI) [41]. During this process, cloudy pixels in an image are filtered out via a cloud-masking tool called "Fmask" [42,43]. A positive NDWI indicates water whereas a zero or negative value indicates a non-water surface. However, to avoid noisy pixels (e.g., saturated soil), small NDWI values (NDWI < 10) are excluded and, hence, the threshold is set as 10.
The wet pixels identified in each image are assigned WSE corresponding to the image acquisition date. After compiling NDWI for all the images, the lowest WSE for each cell is calculated. The lowest WSE per cell estimated from NDWI is, hereafter, called "shoreline elevation" for respective river points (or cells).

Bathymetry Prediction
The fourth task of the RiBaSIC algorithm is bathymetry prediction. It has multiple steps that are presented below. a.
Image Processing: typically, water reflects 10% of the incoming radiation that is limited to visible to NIR spectrum but soil reflects throughout the IR bands [44,45]. The highest reflectance is provided by turbid water [45] which is commonly found in river channels. Additionally, pixels in a river channel are mainly covered by soil that is periodically inundated. Therefore, short wavelength infra-red (SWIR) bands are useful for separating inundated and non-inundated river pixels in an image. So, the ratio of visible and IR (NIR and SWIR) bands are used in this step for predicting bathymetry. Similar to Section 2.1.3, cloud-covered pixels are removed first via the "Fmask" tool. Once the cloudy pixels are removed, Equation (4) b. Generating Library of ADWI and WSE: a library of ADWI and corresponding WSE for each cell within the RiBaSIC extent (also referred to as river points) is generated. The library is used to determine the lowest and highest WSE at which a cell is found to be within a channel (LWSE and HWSE, respectively) along with the highest WSE at which a cell is found to be dry or non-channel (DWSE). If any cell is found to be in the channel at only one instance, then that WSE would be categorized as LWSE for that cell. The ADWI value corresponding to a cell's LWSE is stored as that cell's "Low Water Index" (LWI). c. ADWI versus Depth Correlation: the historic LWL of the nearby streamflow gauge is required to establish the ADWI versus depth correlation. If the riverbed is static (i.e., no significant bed level aggradation-degradation), then the thalweg elevation (lowest elevation of the riverbed) should be less than the LWL. Riverbed elevation for each cell or river point is approximated by two trials using the equations provided below.
1st Trial In the 1st trial, thalweg elevation is assumed the same as LWL. Hence, thalweg elevation, Zt = LWL. Where, LWL = lowest water level from historic data maximum WSE in image acquisition dates Depth factor, (maximum-minimum)WSE in image acqusition dates Then, for any cell that is wet more than once,  Depth during satellite image scene , WSE After analyzing each cell in all the images, ADWIs are sorted and corresponding depths are stored. For each integer ADWI value, the average depth (from all the cells in every image that has the same ADWI) is calculated to create an ADWI versus depth correlation.
2nd Trial Ideally, the thalweg elevation of a river should be lower than LWL. Therefore, the depth values obtained in the 1st trial is used to make an approximation for the depth of thalweg at LWL. This depth value is called delta (Δd) which is approximated as the first one percentile of all the average depth values calculated in the 1st trial. Thus, thalweg elevation in the second trial is estimated as the following: After updating the thalweg elevation (Zt), the same approach as the 1st trial is repeated. By the end of the 2nd trial, the final ADWI versus depth correlation is established. This correlation is used as a lookup table for estimating ADWI bathymetry.
ADWI bathymetry: the ADWI versus average depth correlation is used to derive the river bathymetry. For each cell, the average depth corresponding to LWI (from ADWI versus average depth correlation) is subtracted from LWSE to get the bed elevation. This is hereafter called "ADWI bathymetry".

Post-Processing
The last task of the RiBaSIC algorithm is post-processing the shoreline and bathymetry elevations for the cells obtained in Sections 2.1.3 and 2.1.4. Post-processing is needed for accurately assigning bathymetric, overbank and floodplain elevations of the river reach. It has three main steps that are presented below with a schematic in Figure 2. b. RiBaSIC Bathymetry: the shoreline delineation method using NDWI (please refer to Section 2.1.3) cannot predict the bathymetry in the main channel that is under water during the dry season. Therefore, the shorelines delineated from a dry season cloud-free image is assumed as the main channel. The cells inside the main channel are assigned ADWI bathymetry estimated in Section 2.1.4. The cells outside of the main channel but found wet in at least one image in Section 2.1.3 are considered as overbank areas. The shoreline elevations are assigned for cells in overbank areas. This merged dataset is, hereafter, called "RiBaSIC bathymetry". c. RiBaSIC DEM: the cells that are neither in the main channel nor in overbank areas are considered as floodplain areas. The elevations of the floodplain cells are obtained from a DEM. Finally, RiBaSIC bathymetry and the DEM are merged to create a raster surface that is, hereafter, called "RiBaSIC DEM".

Discharge Estimation from HD Model using a Moderate Quality Global DEM
Discharge estimation using observed WSE and moderate resolution DEM is performed via Hydrologic Engineering Center's River Analysis System (HEC-RAS) model version 5.0.6 that allows 2D HD simulation. The study area for HD modeling is hereafter called a "model extent". Observed WSE is used as the upstream boundary condition, Manning's n values are used to define surface roughness and friction slope at normal depth is used as the downstream boundary condition. The friction slope at downstream is determined by differentiating the second-order channel profile [please refer to Equation (2) in Section 2.1.2] at the downstream location. The methodology for approximating Manning's roughness factors is presented in the next sub-section.

Manning's Roughness Factor Approximation
To perform HD modeling in this study, appropriate Manning's n factors are needed. Since the RiBaSIC algorithm is intended to be applied in regions without prior knowledge of geomorphology and flow, the authors demonstrate an approach to approximate Manning's roughness factors using RiBaSIC bathymetry and channel profile. Based on a review of the literature, several equations are available for predicting Manning's roughness factors from the morphological properties of rivers. A few equations are selected that predict Manning's roughness factors from hydraulic radius and slope (refer to Table 1). These equations are chosen because slope can be obtained from channel profile [assuming channel slope (S) ≈ water surface slope (Sw) ≈ energy gradient (Sf)] and hydraulic radius can be obtained from RiBaSIC bathymetry. Each of these equations is developed for particular types of rivers and they may not be the same as the given study areas. Hence, a Manning's roughness factor approximation methodology is followed for the model extents of this study. It has five steps that are presented below with a schematic in Figure 3.  [48,49] Step 1: hydraulic radius of the cross-sections at upstream of the model extent is calculated.
Step 2: linear slope at the upstream location is estimated. For more details on estimating linear slopes at upstream and downstream locations, the readers are referred to Section 2.1.2 and Appendix A.
Step 3: Manning's n values for the upstream location is calculated using the equations listed in Table  1.
Step 4: any Manning's n value less than 0.010 is considered as "unreasonable" for natural rivers (for the flow within the main channel) and hence discarded. This filtering is important because Manning's n values can be as small as 0.025 for minor streams (width at flood stage < 30 m) although for major streams, the value can be smaller than that of a minor stream of similar characteristics [51]. Moreover, alluvial sand bed rivers could have Manning's n value ranging from 0.018 to 0.035 [52]. The smallest prescribed Manning's n value is 0.016, which is applicable for straight and uniform earthen channels [51]. Therefore, Manning's n value as small as 0.010 is likely to be an erroneous one and needs to be excluded.
Step 5: the average of Manning's n estimates (from step 4) is calculated.
The Manning's roughness factor downstream of the model extent is also estimated following the above-mentioned steps. The upstream and downstream Manning's roughness factors are spatially interpolated along the river length to use in HD modeling.

Study Areas
Four study areas are selected where the RiBaSIC algorithm is tested and 2D HD modeling is performed. The four study areas are Willamette River in Oregon, USA (SA1), Kushiyara River in Bangladesh (SA2), Jamuna River in Bangladesh (SA3), and Solimoes River (SA4) in Brazil. The study areas are chosen to cover a wide range of rivers spread across different geographic locations. The RiBaSIC extent and model extents are the same for SA3 and SA4 but the model extents in SA1 and SA2 are shorter than their RiBaSIC extents to save the computational time of HD simulation. The HD model for each study area is simulated for five consecutive years to determine the efficiency of the proposed methodology for estimating long-term flow patterns. The study areas are shown in Figure  4 and a brief description is presented in Table 2.

Data
Satellite imagery, WSE, and DEM are the main data needed for executing the RiBaSIC algorithm. Landsat satellite image is collected from the United States Geologic Survey (USGS) EarthExplorer database. Observed WSE is collected from USGS, Bangladesh Water Development Board (BWDB), and the SO-HYBAM (previously known as Environmental Research Observatory). The one-arc second (~30 m) SRTM DEM is also collected from USGS EarthExplorer for obtaining channel profile. Additionally, MERIT DEM is collected from the Institute of Industrial Sciences at The University of Tokyo. MERIT DEM is used for representing floodplain elevation because it is derived from SRTM version 2.1 and AW3DDEM (ALOS: Advanced Land Observing Satellite, World 3D-DEM) that has an improved representation of floodplains compared to SRTM [13]. Since the RiBaSIC bathymetry has a spatial resolution of one-arc-second, MERIT DEM is resampled to one-arc second from three-arc second (~90 m) using bilinear interpolation for producing RiBaSIC DEM. Surveyed bathymetry data and observed discharge are also collected from nearby streamflow gauges for validating RiBaSIC algorithm and accuracy assessment of estimated discharge. The available literature is also studied to fill up the gaps in field data.   [53]. 3 Daily discharge at Bahadurabad is calculated using a published rating curve [54]. 4 The date of the bathymetric survey is not mentioned in the navigation chart. 5 Bathymetry survey conducted by USGS for NASA SWOT mission [55] is merged with LiDAR DEM from the Oregon Department of Geology and Mineral Industries (DOGAMI). Note: Please refer to Abbreviations for a list of acronyms.

Results
The RiBaSIC algorithm is developed (see Section 2.1) and applied to produce RiBaSIC DEM for the study areas. The validation of the algorithm and performance of HD simulations using The RiBaSIC DEM are presented in this section.

Validation of the RiBaSIC Algorithm
RiBaSIC algorithm is implemented in the study areas following the methodology presented in Section 2.1 and data listed in Table 2. The spatial resolution of the RiBaSIC bathymetry raster is the same as the Landsat image (30 m). Therefore, the spatial resolution of RiBaSIC DEM (RiBaSIC bathymetry merged with resampled MERIT DEM) is also one arc-second or 30 m. A planform view of the RiBaSIC bathymetries over some reaches in the study areas is presented in Figure 5. It can be noticed that RiBaSIC bathymetry can capture morphological features such as planform shifting, confluence, bifurcation, anabranches, islands, sand bars, etc. The deviation of river planform in RiBaSIC bathymetry and SRTM waterbody data (representing the year 2000, the year of SRTM mission) is also noticeable. Therefore, it can be useful for spatial flood hazard analysis in large, complex, and morphologically active rivers. Additionally, the cross-sectional comparison of RiBaSIC bathymetry at a few locations in SA1, SA2, SA3, and SA4 is presented in Figure 6. It is noticeable that the bed level elevation in SRTM and MERIT DEM in all the cross-sections is significantly higher than the actual bed level. The RiBaSIC estimated bed level elevations although not exact but shows a better approximation than that of SRTM and MERIT DEMs.

Discharge Estimation by a Moderate Resolution Global DEM-Based HD Modeling
The first task for discharge estimation is Manning's roughness factor approximation. For SA1 and SA2, hydraulic radii at upstream and downstream locations are calculated. In SA3 and SA4, considering their reach length to width ratio, a reach-averaged R is calculated. The upstream and downstream linear slopes for each study area is combined with a corresponding hydraulic radius to calculate Manning's roughness factors using the equations listed in Table 1. However, Jarrett's equation resulted in very small Manning's roughness factors (i.e., <0.010) in low gradient rivers (SA2, SA3, and SA4) hence are excluded from roughness factor estimation for those rivers. In fact, Jarrett's formula is developed for steep gradient rivers and the filtering process has identified its unsuitability for these three locations. The upstream and downstream Manning's roughness factors are first calculated and spatially interpolated along the river length to be used in the HD modeling. The result of Manning's roughness factor approximation is presented in Table 3.
The modeled discharge estimates in the study areas are compared with observed or rated discharge obtained from nearby streamflow gauges. Observed and rated discharge are hereafter used synonymously. Simulation of five hydrologic years for all study areas allowed determining the annual and seasonal error of the models. The modeled flow hydrographs along with the observed discharge are shown in Figure 7 and their performance statistics are presented in Table 3.

Discussion
The current study introduces the RiBaSIC algorithm and demonstrates that it can predict the conveyance and planform of a river without predefined geometric properties (e.g., approximate cross-sectional dimension, aspect ratio, etc.). The only observed data needed is WSE at a nearby location, which indicates RiBaSIC can be applied to areas where no other data (e.g., discharge, roughness factors, channel dimensions, bank slope, etc.) is available. However, the performance of the algorithm is limited to the spatial and temporal coverage of the satellite image being used. The length of the time window and the number of images is also important because the larger the range of WSE (high and low stages during the image acquisition) the better the river bathymetry prediction should be. In this study, red, green, blue, NIR, SWIR 1, and SWIR 2 bands of Landsat image are used. The aforementioned Landsat bands are available at 30 m spatial resolution. However, any other satellite image with reasonable radiometric, spatial, and temporal resolution can also be used for bathymetric prediction by the RiBaSIC algorithm. Therefore, further research should be directed towards evaluating the performance of RiBaSIC algorithm with multi mission satellite images such as Sentinel 2 and upcoming Landsat 9. The width of the river (relative to image spatial resolution), the hydrological cycle, and the availability of cloud-free images should also be considered before selecting the image source.
The experiments in this study showed the RiBaSIC algorithm can predict the hydraulic radius of the actively flowing channel using DEM derived slope, satellite image, and WSE. This indicates that bathymetry predicted by RiBaSIC can provide an equivalent conveyance for HD modeling. The Manning's n approximation method followed in this research only requires slope and hydraulic radius, hence, it is well suited for RiBaSIC bathymetry-based HD modeling. Unlike an existing method proposed by Bjerklie, it needs no calibration and is effective in both braided and non-braided rivers [56]. River planform features (i.e., bank elevation, side slope, islands) are also better represented by this algorithm compared to the global DEMs that are often decades old (refer to Figures 5 and 6). Therefore, it is not limited to simplified river assumptions where a river is often assumed single-channel with a predefined shape. This would be especially helpful to produce synthetic bathymetry for rivers with significant morphological activity and to fill data gaps where bathymetry data for a part or entire study area is missing. However, RiBaSIC is not able to predict local scour and pools because the algorithm is expected to provide a general shape of the active channel and is limited to the spatial resolution of the satellite image.
The simulated discharge using the RiBaSIC bathymetry-based model for the study areas showed an annual RRMSE ranging from 14% to 19%. It is noticeable that the performance of long-term simulation on these widely varying rivers yielded similar modeling performance. However, separating the errors by season showed the error is greater in the dry season compared to the wet season. This is primarily because of the small percentage of flow available in the dry season in the study areas. The high error in the dry season could also be due to uncertainty in the RiBaSIC algorithm to predict the bathymetry elevation of the deepest part of the river. However, seasonal variation of the slope and flow roughness in the study areas could be contributing factors as well.
As evident from the results, the RiBaSIC algorithm is effective in predicting conveyance, planform, and estimating discharge using observed WSE and approximated Manning's roughness factor. For absolutely data-poor regions, altimeter-based WSE is a suitable replacement for observed WSE and several methods have been developed for discharge estimation from altimeter-based WSE [31,[57][58][59][60][61][62][63]. The existing methods of discharge estimation from altimeter-based WSE are based upon the conceptual correlation of simplified river geometry and/or prior knowledge of the study area (e.g., discharge, velocity). Based on an intercomparison of five established remote sensing river discharge estimation algorithms, one study reported that the applicability of their studied algorithms varies in truly ungauged reaches [64]. Additionally, the effectiveness of remote sensing and observed data-based bathymetry and discharge prediction algorithms on braided rivers are largely unexplored. The RiBaSIC algorithm presented in this research may, therefore, be useful in truly ungauged and complex rivers via altimeter-based WSE. One major challenge, in that case, would be accommodating the temporal resolution of altimeter-based WSE. Several researchers have overcome that issue by linear interpolation and data synthesis [57,59,62,65]. The expected improvements in WSE measurements from the proposed SWOT mission make it encouraging to extend this research in the future to truly ungauged settings using altimeter-based WSE measurements. SWOT is expected to provide WSE, width, and slope estimates for rivers larger than 100 m [34]. Therefore, the SWOT estimations may provide more accurate and seasonally varied slope approximation that will be useful for the RiBaSIC algorithm. The SWOT data will be useful in directly assigning WSE to wet cells thus enabling spatial and seasonal slope implementation in the current RiBaSIC algorithm. Additionally, SWOT's capability of identifying waterbodies might improve the RiBaSIC algorithm's image processing efficiency from satellite images.

Conclusions
The first objective of this research is to develop an algorithm to predict river conveyance and characterize planform via satellite images and observed WSE. The proposed RiBaSIC algorithm can predict river shape, depth, and characterize planform features. Therefore, this could be useful for areas where limited streamflow gauges are available, but river bathymetric survey data are missing or sparse. It predicted the bed level below LWL, thus, making it applicable for modeling annual events where both dry season and wet season flows are of interest. Additionally, the algorithm is based upon satellite image observation; therefore, updated bathymetric data can be produced (annually) in response to river dynamics such as planform shifting via erosion deposition, cut-off, avulsion, etc. This study focused on utilizing publicly available remote sensing data such as SRTM DEM, MERIT DEM, and Landsat images so that data-poor regions could benefit at no cost. The outcomes of this study also indicate this algorithm may be implemented over absolutely data-poor areas by replacing observed WSE with altimeter-based WSE. Therefore, the effectiveness of the altimeter-based RiBaSIC algorithm for ungauged and complex river systems would be future research avenue.
The second objective of this research is to estimate discharge from observed WSE using a moderate resolution global DEM-based HD model after applying the RiBaSIC algorithm. Global DEMs are not suitable for HD modeling due to lack of conveyance over rivers [66]. On the other hand, HD models need periodic "re-initialization" with observation for the best performance [67]. Hence, it is challenging to rely upon the aging global DEMs for HD modeling. Nevertheless, this study showed that the RiBaSIC algorithm can produce spatially and temporally relevant DEM (using satellite image as an observation) that can be effectively used in HD modeling. The 14%-19% RRMSE in discharge estimation over four non-tidal rivers where the maximum discharge varies from 1700 m 3 /s to 160,000 m 3 /s indicates the proposed RiBaSIC algorithm and Manning's roughness factor approximation method is robust enough to be applied on other non-tidal rivers.
The integration of the RiBaSIC algorithm and 2D modeling shows the application of the presented approach is not limited to discharge assessment only. It can be useful for applications such as inundation mapping, flow distribution in braided channels, flow patterns as a response to channel shifting, etc. The ability of the RiBaSIC algorithm to combine planform and conveyance is also useful for areas where river courses have been altered and wetlands have been converted to croplands. For example, floodplains in Mississippi River alluvial valleys in West Tennessee, Mississippi, and Kentucky show traces of abandoned channels and drained wetlands that are not well represented in moderate resolution SRTM DEM. RiBaSIC can be applied in such areas for correcting DEM and subsequently improving 2D HD modeling. Nevertheless, the RiBaSIC algorithm at its current form employs only one channel profile (slope) but it can also be tested if using seasonally varied slope (if available) has any impact on the results of the RiBaSIC algorithm and corresponding HD modeling performance.

Appendix A.2. Input Data Processing
Three images within the selected time window (i.e., Year 10) is found from a source that has 10 m spatial resolution. The WSE corresponding image acquisition dates, historic hydrograph, and LWL (16 m) are presented in Figure A2. The river polygon is converted to a raster (spatial resolution 10 m) and then to river point shapefile. The river polygon, raster grids, river points, and the cell ID for the river points are presented in Figure A3. Figure A3. (a) River polygon, (b) river raster, (c) river points, and (d) unique cell ID for the river points.
A longitudinal channel profile of the river is also determined using the DEM and shown in Figure A4.

Appendix A.3. Shoreline Delineation
As mentioned in Section 2.1.3, the wet pixels in three images are determined using NDWI and the shoreline elevations are calculated. For example, cell ID 7 (third column from the left in the second row) is wet in all three images; hence, the shoreline elevation is the smallest WSE at which it is wet which is 18 m. NDWI for wet pixels in the selected images and resulting shoreline elevation for cells from three images are shown in Figure A5.

Appendix A.4. Bathymetry Prediction
The steps of river bathymetry prediction using ADWI and WSE is presented here. a. Image Processing: the ADWI values for channel cells in three images are shown in Figure A6. b. Generating Library of ADWI and WSE: the HWSE, LWSE, DWSE, and LWI is calculated in this step are presented in Figure A7. For cell ID number 7, from image processing:  the highest WSE at which it is in the channel is 20 m, therefore, HWSE = 20 m;  the lowest WSE at which it is in the channel is 18 m, therefore, LWSE = 18 m;  the cell is not found dry in any image, therefore, DWSE = null;  ADWI at 18 m WSE is 60, therefore, LWI = 60. Figure A7. Lowest water surface elevation (WSE) at which a cell is found to be within a channel (LWSE), highest WSE at which a cell is found to be within a channel (HWSE), highest WSE at which a cell is found to be dry (DWSE), and low water index (LWI) for the cells calculated from ADWI.
c. ADWI versus Depth Correlation: 1st Trial Zt = LWL = 16 (see Figure A1   The ADWI versus average depth values calculated from all the river points (cells) from three images is calculated and presented in Figure A10.  b. RiBaSIC Bathymetry: the main channel of the Fusion River is obtained from a dry season satellite image as shown in Figure A13. The cells inside the main channel are assigned. Slope transferred ADWI bathymetry and shoreline elevations are assigned for the cell outside of the main channel to get the RiBaSIC bathymetry as shown in Figure A13. Cell ID 7 falls inside the main channel; therefore, ADWI bathymetry has the same elevation as slope transferred ADWI bathymetry (5 m). c. RiBaSIC DEM: the empty cells in RiBaSIC bathymetry are filled by elevations extracted from the DEM to produce RiBaSIC DEM as shown in Figure A14. It can be seen that for cell ID 7, the DEM elevation is 9m but that has been modified to 5 m in RiBaSIC DEM. The RiBaSIC algorithm is implemented on the Fusion River for Year 10. A cross-sectional comparison at Section 2 (along row 2; shown in the dotted box) is presented in Figure A15. It can be noticed that the RiBaSIC bathymetry predicted the bed level below the LWL whereas the shoreline elevations predicted overbank areas. Therefore, RiBaSIC DEM could represent an improved channel conveyance that also reflects the planform Year 10.