Next Article in Journal
Challenges of Comparing Marine Microbiome Community Composition Data Provided by Different Commercial Laboratories and Classification Databases
Previous Article in Journal
Water Quality Predictions Based on Grey Relation Analysis Enhanced LSTM Algorithms
Previous Article in Special Issue
The Morphological Evolution of a Step–Pool Stream after an Exceptional Flood and Subsequent Ordinary Flow Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of River Morphology on the Remote Sensing Based Discharge Estimation: Implications for Satellite Virtual Gauge Establishment

1
Shaanxi Key Laboratory of Earth Surface System and Environmental Carrying Capacity, Northwest University, Xi’an 710127, China
2
China Water Resources Pearl River Planning Surveying & Designing Co., Ltd., Guangzhou 510630, China
3
College of Urban and Environmental Sciences, Northwest University, Xi’an 710127, China
4
Institute of Earth Surface System and Hazards, Northwest University, Xi’an 710127, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(23), 3854; https://doi.org/10.3390/w14233854
Submission received: 20 September 2022 / Revised: 29 October 2022 / Accepted: 24 November 2022 / Published: 26 November 2022

Abstract

:
Monitoring of river discharge is a key process for water resources management, soil and water conservation, climate change, water cycling, flood or drought warning, agriculture and transportation, especially for the sustainable development of rivers and their surrounding ecological environment. Continuous and comprehensive discharge monitoring was usually impossible before, due to sparse gauges and gauge deactivation. Satellite remote sensing provides an advanced approach for estimating and monitoring river discharge at regional or even global scales. River morphology is generally considered to be a direct factor that affects the accuracy of remote sensing estimation, but the specific indicators and the extent to which it affects the estimation accuracy have not yet been explored, especially for medium to small rivers (width < 100 m). In this paper, six sites with hydrological gauges in the upper Heihe River Basin (HRB) of northwestern China and the Murray Darling Basin (MDB) of southeastern Australia were selected as the study cases. River discharge was estimated from Landsat imagery using the C/M method accordingly. River gradient, sinuosity, and width were obtained from Digital Elevation Model data for each site. Global Surface Water Dataset (GSWD) was also employed for indicating the dynamic status of river morphology. A series of methods were applied to analyze the influence of river morphology on estimation accuracy qualitatively and quantitatively, based on which we established inference about the theory of selecting satellite virtual gauges (SVGs). The results confirm the feasibility of the C/M method for discharge estimation, with the accuracy affected by multiple river morphological indicators. Among them, river width was found to be the most significant one. Moreover, water occurrence and water extent extracted from GSWD also have impact on the discharge estimation accuracy. Another independent river section in MDB was set as an example to demonstrate the reasonability of the established theory. It is anticipated that this study would promote the application of remote sensing for discharge estimation by providing practical guidance for establishing appropriate SVGs.

Graphical Abstract

1. Introduction

River discharge, or called streamflow, is a vital component of the hydrological cycle. Therefore, accurately monitoring river discharge is crucial for a variety of hydrological studies and applications, such as water resources management, soil and water conservation, flood/drought monitoring under a changing climate [1,2]. River discharge monitoring via traditional stream gauging networks is usually the most appropriate approach of estimating discharge, which however is limited to point locations along the river in accessible and developed reaches [3]. On one hand, gauges are distributed unevenly and sparsely worldwide. One the other hand, the number of gauges is further decreasing in recent years due to various reasons, such as extreme weather, lack of maintenance, etc. [4]. This makes it even more difficult for natural rivers to be monitored continuously and comprehensively [5].
Remote sensing technology provides an alternative and advanced discharge estimation approach due to its unique advantages in obtaining river water extent, width, water level, and other river characteristics timely and stably [6,7], which greatly extends the source of the traditional hydrological data [8]. Gleason et al. [9] suggested that for rivers with limited observations, remote sensing could act as important supplementation and expansion. While for rivers with completely no observation, it would be the core data source. Therefore, river discharge monitoring with satellite remote sensing is a burgeoning field attracting increasing attentions [10]. Since the Prediction in Ungauged Basins (PUB) was proposed in 2003 [11], satellite remote sensing has been more and more widely used for river discharge estimation through establishing virtual gauges [12]. So far, a series of remote sensing-based discharge estimation methods have been developed. Gleason et al. [9] made a comprehensive review of these methods regarding the availability of gauge data (gauged, semi-gauged, regionally ungauged, and totally ungauged basins).
The rating curve method is a traditional approach that estimates river discharge by establishing a relationship between available observed discharge and one of the remotely sensed river characteristics, such as water level [13] and water area [14]. This approach is suitable for rivers that have in situ observations, but it cannot be applied to ungauged rivers. Other studies used remotely sensed data as major input to estimate discharge based on hydraulic geometry [7,15,16,17,18,19], or flow wave propagation models [20,21]. However, some key information required by these methods, such as river bathymetry and riverbed morphology, is usually unavailable. Although it may be obtained from remote sensing techniques [22,23], extensive and various remote sensing data are needed, which may limit their applicability in some rivers.
The Calibration/Measurement (C/M) method, first proposed by Brakenridge et al. [24], can be considered as an improved rating curve approach. It estimates discharge based on the close correlation between observed discharge and the ratio of a land pixel for calibration (C) and a water pixel for measurement (M). It was first tested using Advanced Microwave Scanning Radiometer for EOS (AMSRE), and then extended to MODIS imagery for estimating river discharge for different river basins [25,26,27,28]. Due to the coarse spatial resolution of AMSRE and MODIS data, these studies only focused on large rivers. Other studies, therefore, tried to extend the C/M method to medium-to-high resolution images, such as Landsat and Sentinel-2 [29,30], for estimating discharge of medium-to-small rivers (width < 100 m). These studies demonstrate the feasibility of utilizing various remote sensing data for estimating river discharge for those rivers that have limited observations due to various reasons.
As more and more remotely sensed data are becoming available, every single river is likely to be fully covered by various data. In that case, how to locate an ideal reach for establishing satellite virtual gauge (SVG) to properly turn remote sensing signals to discharge information? Existing studies suggested that river morphology is likely to be an important factor that affects discharge estimation accuracy. For example, Durand et al. [31] compared five different methods for estimating discharge in 19 selected rivers and found obvious accuracy difference due to different river morphology. Brinkerhoff et al. [32] demonstrated the important influence of river morphology on discharge estimation using a hydraulic method. Hou et al. [12] concluded that wide channels with strong temporal variations, broad floodplains and multiple braided or anastomosing channels provided the best conditions for SVG. They tested a series of different discharge estimation methods at different river reaches globally. However, as the C/M method is not a hydraulic method, would the river morphology also affect its accuracy? How much uncertainty would be induced by river morphology in the C/M modelling process? Answers to these questions are still not clear.
This study aims to investigate the influence of river morphology on estimating discharge using remote sensing. We focus on medium to small rivers (width < 100 m) through implementing C/M method with medium-resolution Landsat imagery, considering that smaller rivers are more extensively existing and ungauged. We try to investigate how river morphological indicators, including river gradient, sinuosity, and width, affect discharge estimation accuracy through statistical analysis, based on which we try to provide guidance for selecting appropriate river reaches to establish SVGs through a case demonstration. As most of the current remote sensing-based discharge estimation studies overlooked the site selection, this study is expected to help reexamine how to carry out C/M discharge estimation more effectively, especially when remote sensing data are usually available for the whole river. As the C/M method requires a time series of discharge data as the input, we select a series of reaches located near several documented gauges as test sites. Our objectives include: (1) qualitatively and quantitatively revealing the impact of river morphology on the accuracy of C/M method, and (2) providing guidance for selecting appropriate sites for establishing SVGs based on finding of the above analysis.

2. Materials and Methods

2.1. Materials

2.1.1. Study Area

As shown in Figure 1, two study areas, the Heihe River Basin (HRB) (Figure 1b), China, and the Murray Darling Basin (MDB) (Figure 1c), Australia, both in semi-arid region, were selected. For each study area, three gauges with observed hydrological records were selected. For the HRB, Gauge Qilian at the Babao River, Gauge Zhamashike at the Heihe River, and Gauge Sunan at the Liyuan River were selected. For each gauge, four river sections (called study sites hereafter) were selected nearby. They were marked as 1-1, 1-2, 1-3, 1-4, 2-1, 2-2, 2-3, 2-4, 3-1, 3-2, 3-3 and 3-4 (Figure 1d–f). For the MDB, Edwards River at Leiwah (Gauge ID: 409035), downstream of Balranald Weir (Gauge ID: 410130), and Murray River at Below Wakool Juction (Gauge ID: 414200) were selected, along with their corresponding nearby study sites marked as 4-1, 4-2, 4-3, 4-4, 5-1, 5-2, 5-3, 5-4, 6-1, 6-2, 6-3 and 6-4 (Figure 1g–i). These sites were selected where there are no tributaries to ensure that the streamflow in these sites is consistent with that recorded at the gauges. The size of site (21 × 21 pixels) was determined considering the computation efficiency of C/M modelling. Darling River at Wilcannia Main Channel (Gauge ID: 425008) (Figure 1c) was selected to demonstrate the implications for SVG establishment.

2.1.2. Data

Observed discharge in the HRB were obtained from Li et al. [29]. Observed discharge in the MDB were obtained from New South Wales Office of Water and Victorian Department of Environment and Primary Industries, after being converted from Megaliters/day (ML/d) to m3/s (1ML/d = 0.011574 m3/s). Their hydrographs are displayed in Figure 2. The modelling period and validation period were listed in Table 1. As the MDB experienced a 10-yr dry period with very low streamflow observed (2000–2009), which was then reversed later (2010–2011) with high streamflow [33], we selected two periods for modelling and two periods for validation for the MDB sites to ensure that both high flow and low flow periods have been included for modelling and validation. For each site, we used Level 2 atmospherically corrected surface reflectance of Landsat data, which were acquired from Google Earth Engine (GEE) [34]. We checked the data for each site, and manually picked out those cloud contaminated ones. Only the NIR band of Landsat was employed as the input of C/M method because of its sensitivity to surface water.
MERIT HYDRO dataset [35,36] derived from the global hydrologically corrected DEM product (MERIT DEM), was used in this study to generate static river morphological factors, including river gradient, sinuosity, and river width. River width in the MERIT HYDRO dataset is derived through integrating flow direction and SRTM Water Body Database (SWBD) [37], which was later further improved to handle subpixel water fraction [36], enabling width estimation for rivers smaller than 90 m. It is noted that this value is a proxy of river width that does not take river channel and water level variation into consideration. In addition, we adopted water dynamic layers from Global Surface Water Dataset (GSWD) developed by European Joint Research Centre (JRC) [38] as indicators for the dynamic river morphology, which were used for qualitative analysis. Specifically, Water Occurrence and Maximum Water Extent of this dataset were used as indicators for overall river water dynamics. Water Occurrence stands for the frequency with which water was present. It shows where surface water occurred between 1984 and 2019 and provides information concerning overall water dynamics. Maximum Water Extent provides information on all the locations ever detected as water over the 36-year period.

2.2. Methodology

2.2.1. Extraction of River Morphology

In this study, river morphology indicators including river gradient, sinuosity, and width. River network is first needed before these indicators can be derived. We employed the river network data in the MERIT HYDRO dataset directly and processed the missing river network in some of the sites from MERIT DEM using the hydrology tools in ArcGIS. Based on the river centerline of the river network, river gradient (G) in this study refers to the ratio of elevation difference (H) to the horizontal distance (L). Sinuosity (S) was calculated as the ratio of the linear distance (D) to the actual river length (l). Extracting river width from MERIT DEM was following Yamazaki et al. [36]. River centerline pixels were first identified, whose corresponding river width were calculated accordingly. For each site, there are a series of river centerline pixels. Among them, the one that has the largest river width was selected, whose river width was taken as that of site. A detailed demonstration is referred to Figure A1 in the Appendix A.

2.2.2. C/M Modelling and Validation

The basic principle of C/M method is to calculate the ratio between a stable land pixel for calibration (C) and a pixel within river for measurement (M), and then apply a linear regression between C/M series and observed discharge series to give discharge information at a site [39]. In this study, time series of Landsat NIR band were used as the input of the C/M model, together with the corresponding gauge observations. As land cover is generally changing over time, to obtain a more stable C pixel, a coefficient of variation (CV) was computed and applied to each pixel, considering that CV can indicate the variability of surface reflectance [40]. Lower CV represents higher stability of land cover. Therefore, we considered those pixels that have low CV values (lower than the 5th percentile) to be stable pixels, and took their average as a surrogate of C pixel [39]. For each Region of Interest (ROI), all pixels were considered as M candidates, and a time series of C/M was calculated for each candidate. The correlation between its C/M value and the observed discharge was quantified by Pearson’s correlation coefficient (r) and validated with the corresponding p-value. Finally, the candidate pixel with the highest r value and a p-value < 0.05 (correlation is significant at the 95% confidence level) was chosen to be the optimal location of M. This ensures the C/M model built for this ROI have the best performance. In this study, each selected site near the gauge was taken as an ROI.
Once the C and M pixels are determined, a linear relationship would be established between the C/M ratio and discharge, which could then be employed to estimate river discharge from remote sensing data. Corresponding observed discharge of the validating period were used as the reference to evaluate the accuracy of C/M method in each site. Relative Root Mean Square Error (RRMSE) (Equation (1)) and Nash-Sutcliffe Efficiency coefficient (NSE) (Equation (2)) were used for accuracy assessment. In Equation (1), n stands for the number of validation data of each site, Qrs and Qg stand for the estimated discharge and observed discharge, respectively. Q g ¯ is the mean value of the observed discharge. In Equation (2), the T is the total observation times, Q g t and Q r s t represent the observed value and estimated value at time t. RRMSE represents the normalized RMSE value. Lower RRMSE indicates higher accuracy. NSE is another widely used indicators for hydrological studies. Its value ranges from negative infinity to 1, with larger value representing more stable and reliable modelling results.
RRMSE = t = 1 n Q rs   Q g 2 n Q g ¯ ,
NSE = 1 t = 1 T ( Q g t Q rs t ) 2 t = 1 T ( Q g t Q g ¯ ) 2 ,

2.2.3. Qualitative and Quantitative Analysis

Qualitative analysis is carried out by checking how C/M modelling accuracy varies with surface water dynamics, and river morphological indicators. Quantitative analysis is conducted using the Multivariate Analysis of Variance (MANOVA), a member of the General Linear Model. It is one of the most common multivariate statistical methods used for investigating whether two or more variables would significantly affect the observed variables [41]. It calculates an F value for each variable, which is taken as the ratio of its independent effect and the random effect. The independent effect is measured using Mean of Squares for factor A (MSA) which is calculated through averaging the inter-group sum of squares. The random effect is measured using Mean of Squares for Error (MSE) which is calculated through averaging the intra-group sum of squares. A significant test sig can then be derived according to the distribution of F values. A large F generally causes small sig value. If the sig is less than the significant level (usually set as 0.05), it suggests the pass of the significance test [42]. In this study, all 24 selected study sites were considered as independent samples with a series of river morphological indices (i.e., river width, sinuosity, and gradient) and accuracy indices (i.e., RRMSE and NSE), which were used as the input of MANOVA. Through comparing the F value of each river morphological indices, MANOVA is able to identify the factor that affects the discharge estimation most, with the largest F value and the least sig.

2.2.4. Demonstration of Satellite Virtual Gauge Selection

The above qualitative and quantitative analysis gives inference about how river morphology affects C/M discharge estimation. We employ another independent gauge in the MDB (Gauge 425008) to conduct a practice of applying the implications on the SVG establishment. 19 ROIs that have the same size with previously defined ROIs near the gauge are selected as candidate SVGs. We then try to find the optimal one based on the findings of the above analysis. C/M modelling and validation are carried out for all selected ROIs and the results are cross compared to confirm whether the one inferred by the river morphology is the optimal one.

3. Results

3.1. Descriptions of River Morphology

As river network in the MERIT HYDRO dataset was generated from DEM data, which sometimes embraces errors. We used river water extent visually inspected from Landsat imagery as the assistant to correct the river network. Figure 3 shows the corrected river networks in the HRB and the MDB. Accordingly, the river gradient, sinuosity, and river width of each study site were calculated and listed in Table 2. The overall gradient of the HRB is larger than that of the MDB, leaving its sites generally having higher gradients that those of the MDB. For example, site 2-2 near Gauge Zhamashike has a gradient of 39.07 m/km. Meanwhile, the difference of gradient among four sites of the same gauge can be large (e.g., gradients for Gauge Sunan). As the rivers in the HRB are relatively straight, sinuosity of the MDB sites is obviously higher than that of the HRB sites. River widths of the study sites are around 60~70 m. Sunan sites have a width less than 30 m, and one of the MDB sites have a width slightly greater than 100 m.

3.2. Discharge Estimation Using C/M Method

Pearson correlation coefficient (Pearson’s r) between C/M value for each candidate M pixel and observed discharge of the 6 hydrological gauges (24 sites in total) selected in the HRB and MDB is shown in Figure 4, with red dots representing the finally selected location of the M pixel. It is observed that the M pixels tend to locate at the edge of the river channel rather than the center, probably because the pixels in the center are always wet, and their reflectance thus changes relatively small as the streamflow increases. On the contrary, pixels on the river edge would experience more intensive reflectance variation as water level changes, making them more sensitive to discharge variation. This is consistent with the findings of Tarpanelli et al. [26].
Figure 5 shows the scatter plots of the optimized C/M value and the observed discharge. The established models and their corresponding coefficient of determination (R2) are also displayed. The observed discharge and C/M series are overall positively correlated for all sites as indicated by Pearson’s r and coefficient of determination (R2), particularly for sites near Gauge Qilian and Gauge Zhamashike in the Heihe River Basin, Gauge 409035 and Gauge 414200 in the Murray Darling Basin (R2 up to 0.7). However, for Gauge 410130 and Gauge Sunan, the correlation is noted to be relatively lower, which suggests considerable uncertainties may exist in discharge estimation for these sites. Moreover, four sites selected for each gauge may demonstrate quite different accuracy, with Pearson’s r differing up to 0.4 (e.g., Pearson’s r from 0.33 to 0.77 for Gauge Sunan).
Figure 6 shows the hydrograph of the estimated discharge and the observed discharge for the validation period. Estimated discharge from four selected sites and observed discharge from corresponding gauge exhibit general consistent pattern. However, when the observed discharge is high, the C/M method tends to underestimate the discharge. Table 3 shows the modelling accuracy in forms of RRMSE and NSE. Overall, river discharge of Gauge Qilian, Gauge Zhamashike in the HRB and Gauge 409035 and Gauge 414200 in the MDB has been estimated relatively accurately using the C/M method, while the results for the other two gauges are less accurate. It is noted that for each gauge, its four selected sites exhibit different accuracy.

3.3. Influence of River Morphology on Discharge Estimation

3.3.1. Qualitative Analysis

Water Occurrence and Maximum Water Extent obtained from the JRC GSWD product (Figure 7) were employed for analyzing the impact of surface water dynamics on river discharge estimation. When water Occurrence within the river channel basically exceeds 90%, this means that the permanent water has been overall properly detected by the Landsat imagery. When the Maximum Water Extent is obviously larger than the channel, this means that the water dynamics have been adequately captured by historic Landsat data. Among the four river sections that have JRC data covered, study sites (3-1, 3-2, 3-3 & 3-4) near Gauge Sunan on Liyuan River have relatively lower accuracy, with RRMSE greater than 1 and NSE less than 0 (Table 3). It is observed that maximum water extent is small and water occurrence is relatively low here, with the maximum occurrence equals to 71%, and most values close to 50%. It is obvious that the river water is mostly constrained within the river channel. This suggests that surface water cannot be observed clearly by Landsat, and its dynamics is limited. This means that the water dynamics here could not be properly captured by Landsat, and thus the discharge cannot be precisely estimated, as proved in Table 3. Sites near Gauge 409035 and Gauge 410130 in the MDB are in a similar position but with a slightly higher accuracy as those near Gauge Sunan.
On the contrary, sites near Gauge 414200 have relatively higher accuracy, with some obtaining RRMSE less than 0.55 and NSE greater than 0.60 (Table 3). For this site, Water Occurrence within the river channel generally exceeds 90%, and Maximum Water Extent is obviously larger than the channel. This implies that these sites have relatively higher quality Landsat observations, which adequately reveal its extensive surface water dynamics. Therefore, we conclude that available high-quality Landsat observations and surface water dynamics are likely to be two important determinants for the accuracy of discharge estimation. Since these determinants could be inferred from Water Occurrence and Maximum Water Extent layers in the JRC GSWD, the GSWD product can serve as important support that provides prior knowledge about the suitability of SVGs.

3.3.2. Quantitative Analysis

Multivariate analysis of variance was conducted to analyze how the three river morphological indicators (i.e., gradient, sinuosity, and width) affect the accuracy of river discharge estimation. As shown in Table 4, river width was found to be the most important factor affecting the discharge estimation, in forms of either RRMSE or NSE. It is also the only factor that passed the significance test (Sig < 0.05), which means the other two factors cannot be proved to be significantly affecting the discharge estimation accuracy. Nevertheless, it seems that the sinuosity tends to affect the estimation accuracy RRMSE (high F and low Sig value for RRMSE) and the gradient tends to affect the stability of the model NSE (high F and low Sig value for NSE). Even though river gradient and sinuosity did not pass the significance test of multivariate analysis of variance, it seems that the two indicators also affect the estimation accuracy for some of the selected study sites. Table 4 also shows the correlation coefficient (r) between RRMSE/NSE and the three river morphological indicators. It is also found that river width is the factor that has the closest relationship with discharge estimation accuracy, in forms of either RRMSE or NSE.

3.3.3. Practice on Another Independent Gauge

The above results indicate that the river width is a major factor that influences the discharge estimation when using the C/M method, and the surface water dynamics would also affect the accuracy, with high Water Occurrence and large Maximum Water Extent benefiting discharge estimation. Overall, higher estimation accuracy tends to appear at the sites with larger river width (mostly larger than 50 m). We employed another independent gauge in the MDB (Gauge 425008) to conduct a practice of applying the implications on the SVG establishment. The observed discharge of Gauge 425008 can be found in Figure 2.
19 ROIs were selected near the Gauge 425008 as candidate SVGs, with the same size as previously defined ROIs (21 × 21 pixels) (Figure 8a). These ROIs were selected based on the overall river width (>50 m) and surface water dynamics (Figure 8b,c) to ensure that the river section in the ROI is wide enough and the river water dynamics have been properly captured by historical Landsat data. River morphological indicators of these ROIs were calculated and listed in Table 5. Among them, ROI 5 has the largest river width and the highest sinuosity, a relatively higher water occurrence and large maximum water extent. According to the implications we got from the precious experiments, we infer ROI 5 to be the most suitable location as the SVG for the C/M method.
In order to further verify our inference, we calculated the Pearson’s r between the C/M ratio and the observed discharge using each pixel of the whole study site as a candidate M (Figure 9a). It is observed that pixels with high correlation (Pearson’s r > 0.7) generally fall into the previously established 19 ROIs (Figure 9b). For each ROI, we selected the pixel that has the highest Pearson’s r as the optimized M pixel, and constructed a C/M model to estimate the discharge for Gauge 425008. The accuracy of each ROI is listed in Table 6, indicating that the highest accuracy occurs in ROI 5 (the estimated discharge using ROI 5 and observed discharge for Gauge 425008 is shown in Figure 9c). A chart showing the difference of the estimated discharge for different ROIs can be found in Figure A2 in the Appendix A. This suggests that ROI 5 is the best location for establishing the SVG for this river reach, which further verifies the conclusions we made from the study cases in China and Australia. Except for ROI 5, ROI 2, ROI 4, and ROI 19 showed relatively high estimation accuracy. They could be used as excellent substitutions to ROI 5, when it is continuously obscured by clouds during the modelling period. As the C/M method is an iterative optimization method whose computation cost increases drastically as the size of virtual gauge increases, quickly locate a ROI that is much smaller but generates decent accuracy would be helpful for establishing SVGs for river discharge monitoring, especially for large scale monitoring where there are massive SVGs to be established.

4. Discussion

Based on experiments conducted in a series of sites selected from China and Australia, we tried to investigate the influence of river morphology on discharge estimation using the popular C/M method. We select river gradient, river sinuosity and river width as three indicators of river morphology, note that all three indicators are static variables that do not link to river water dynamics. Even though we only selected six reaches from China and Australia considering the data availability and workload, we applied four sites to each reach and thus built 24 cases that can provide fair support to investigate the scientific questions we aim to answer. We identified that river width is the most import factor that controls the reliability of discharge estimation using the C/M method. It is also noted that the other river morphology indicators, such as river gradient and sinuosity, may also affect the accuracy of the C/M method. This finding is overall consistent with existing studies (e.g., [12,31,32]). Even though they focused on different discharge estimation methods, there is a consensus that river morphology is critical for discharge estimation, and therefore should be taken into consideration when selecting SVG reaches.
River water dynamics, as constrained by the river morphology, were also found to be an important factor that may affect the C/M accuracy. As the river water dynamics can be easily obtained from remote sensing imagery and even readily available in several surface water products (e.g., GSWD), they can serve as useful resource for preselecting potential suitable SVG locations. The finding of this study was further applied an independent river reach. It was proved that it can help quickly locate a smaller ROI that generates the best C/M model for discharge estimation. Considering that the C/M method is a iterative algorithm whose computation cost increases drastically with the size of ROI, this study provides useful practice guidance for establishing SVGs for river discharge monitoring, especially for large scale monitoring where many SVGs are needed.
There are several issues that may limit our conclusions or need further attentions. First, extracting river morphology from remote sensing images and DEM may introduce uncertainties, especially in arid and semi-arid areas, where the river is either meandering or relatively narrow. Second, gradient and sinuosity failed to pass the significance test of multivariate analysis of variance, probably due to the limited number of cases. Third, water and land color change will affect the discharge estimation, as the C/M method utilizes the ratio of NIR reflectance of a C pixel and a M pixel [43]. In order to weaken the impact of land cover change (in particular vegetation phenology), we adopted the CV value of each pixel along the time series, and considered the 5th percentile of low CV pixels as land pixels. These pixels were averaged to represent the C pixel. Through this process, it is believed that those land pixels whose NIR reflectance changes drastically will not be used for the C/M modelling. M pixel was determined through an iterative optimization process. Watercolor, in particular the turbidity, would certainly affect the NIR reflectance, with higher turbidity usually causing higher NIR reflectance. Therefore, the M pixel embraces different levels of uncertainty according to the turbidity level in the SVG. Higher turbidity makes the M value greater than it should be when water is clean, and thus may lead to model overestimation. One thing that further complicates the situation is the uncertain and unstable relationship between turbidity and streamflow. Therefore, in addition to river morphology, watercolor may be another important factor that affects the spectral-based discharge estimation method. The SVG site selection may need to consider the watercolor variation as well, in particular for those rivers whose turbidity varies fiercely across space and time.
We only adopted the C/M method in this study. There are several different discharge estimation methods that are based on quite different principles. Whether the river morphology also affect the accuracy of these methods is still not clear. In the future, the input of experimental river sections, observed discharge data, and remote sensing images can be increased, and multi-source remote sensing images may also be integrated into different discharge estimation models, to comprehensively investigate how river morphology affects different approaches of discharge estimation.

5. Conclusions

When establishing SVGs and implementing discharge monitoring in the basin, finding the optimized location that is able to provide accurate river hydrological parameters for estimating discharge is critical for improving the efficiency of models. Combined with surface water dynamics derived from Landsat imagery and static river morphological indicators derived from DEM data, this paper qualitatively and quantitatively investigated the influence of river morphology on the C/M discharge estimation model, especially for medium to small rivers (width < 100 m). Through experimenting at 24 sites in semi-arid regions of China and Australia, guidance for selecting appropriate river reaches to establish SVGs in another river section (near Gauge 425008) were provided. Conclusions of this study are as follows.
(1)
The observed discharge and the C/M series in all the selected river sections showed medium-to-high positive correlation. For each gauge, the estimated results of its four study sites exhibit a moderate but different accuracy. This suggests that the C/M method is overall reliable for all the selected study sites.
(2)
Estimation results are generally better in study sites with higher water occurrence and larger maximum water extent. This implies river sections with high-quality Landsat observations and extensive surface water dynamics are important prerequisites for reliable discharge estimation. Therefore, some surface water products (e.g., GSWD) could be useful resources for assisting SVG location.
(3)
It is found that river width has the most significant influence on the discharge estimation accuracy and model stability (sig < 0.05). Even though river gradient and sinuosity did not pass the significance test of multivariate analysis of variance, it seems that the two indicators also affect the estimation accuracy for some of the selected river sections. This suggests that river morphology is truly an important factor that influencing the C/M river discharge estimation method.
(4)
Through testing at another independent river reach, our findings about how river morphology affects discharge estimation and how can they be applied for establishing SVGs are verified, which further confirmed the significance of this study. It provides useful practice guidance for establishing SVGs for river discharge monitoring, especially for large scale monitoring where many SVGs are needed.

Author Contributions

Conceptualization, C.H. and Z.S.; methodology, Z.S.; software, Z.S.; validation, Z.S.; formal analysis, Z.S.; investigation, Q.C.; data curation, Q.C.; writing—original draft preparation, Z.S.; writing—review and editing, C.H.; visualization, Q.C.; supervision, C.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China [grant number U2243205], the National Key R&D Program of China [grant number 2017YFC0404302], and Shaanxi Natural Science Foundation [grant number 2021JM-314].

Acknowledgments

We acknowledge the New South Wales Office of Water (http://realtimedata.water.nsw.gov.au/ (accessed on 10 December 2021)) and Victorian Department of Environment and Primary Industries (http://data.water.vic.gov.au (accessed on 10 December 2021)) for providing gauge discharge data. MERIT DEM data were obtained from Yamazaki (http://hydro.iis.u-tokyo.ac.jp/~yamadai/MERIT_Hydro/ (accessed on 10 December 2021)). We would like to thank Tarpanelli for helping us with the code of C/M algorithm. We are grateful to two anonymous reviewers for their helpful comments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Selection of representative width value based on width from MERIT HYDRO in this study. River centerline pixels were first identified, and each had corresponding river width information. As for the selected river reach, there are a series of river centerline pixels. Among them, the one that has the largest river width was selected, and this largest river width was taken as the river width of the river reach (site identified by the pink ROI 2-2 for example, the white background stands for 0 value.).
Figure A1. Selection of representative width value based on width from MERIT HYDRO in this study. River centerline pixels were first identified, and each had corresponding river width information. As for the selected river reach, there are a series of river centerline pixels. Among them, the one that has the largest river width was selected, and this largest river width was taken as the river width of the river reach (site identified by the pink ROI 2-2 for example, the white background stands for 0 value.).
Water 14 03854 g0a1
Figure A2. Estimated discharge using all 19 ROIs and observed discharge in the site near Gauge 425008.
Figure A2. Estimated discharge using all 19 ROIs and observed discharge in the site near Gauge 425008.
Water 14 03854 g0a2

References

  1. Sheffield, J.; Wood, E.F.; Pan, M.; Beck, H.; Coccia, G.; Serrat-Capdevila, A.; Verbist, K. Satellite remote sensing for water resources management: Potential for supporting sustainable development in data-poor regions. Water Resour. Res. 2018, 54, 9724–9758. [Google Scholar] [CrossRef] [Green Version]
  2. Xu, M.; Kang, S.; Chen, X.; Wu, H.; Wang, X.; Su, Z. Detection of hydrological variations and their impacts on vegetation from multiple satellite observations in the Three-River Source Region of the Tibetan Plateau. Sci. Total Environ. 2018, 639, 1220–1232. [Google Scholar] [CrossRef] [PubMed]
  3. Kebede, M.G.; Wang, L.; Yang, K.; Chen, D.; Li, X.; Zeng, T.; Hu, Z. Discharge estimates for ungauged rivers flowing over complex high-mountainous regions based solely on remote sensing-derived datasets. Remote Sens. 2020, 12, 1064. [Google Scholar] [CrossRef] [Green Version]
  4. Biancamaria, S.; Hossain, F.; Lettenmaier, D.P. Forecasting transboundary river water elevations from space. Geophys. Res. Lett. 2011, 38, L11401. [Google Scholar] [CrossRef] [Green Version]
  5. Davids, J.C.; Rutten, M.M.; Pandey, A.; Devkota, N.; van Oyen, W.D.; Prajapati, R.; van de Giesen, N. Citizen science flow—An assessment of simple streamflow measurement methods. Hydrol. Earth Syst. Sci. 2019, 23, 1045–1065. [Google Scholar] [CrossRef] [Green Version]
  6. Döll, P.; Douville, H.; Güntner, A.; Müller Schmied, H.; Wada, Y. Modelling freshwater resources at the global scale: Challenges and prospects. Surv. Geophys. 2016, 37, 195–221. [Google Scholar] [CrossRef] [Green Version]
  7. Huang, Q.; Long, D.; Du, M.; Zeng, C.; Qiao, G.; Li, X.; Hou, A.; Hong, Y. Discharge estimation in high-mountain regions with improved methods using multisource remote sensing: A case study of the Upper Brahmaputra River. Remote Sens. Environ. 2018, 219, 115–134. [Google Scholar] [CrossRef]
  8. Lettenmaier, D.P.; Famiglietti, J.S. Hydrology: Water from on high. Nature 2006, 444, 562–563. [Google Scholar] [CrossRef]
  9. Gleason, C.J.; Durand, M.T. Remote sensing of river discharge: A review and a framing for the discipline. Remote Sens. 2020, 12, 1107. [Google Scholar] [CrossRef] [Green Version]
  10. Wang, L.; Sichangi, A.W.; Zeng, T.; Li, X.; Hu, Z.; Genanu, M. New methods designed to estimate the daily discharges of rivers in the Tibetan Plateau. Sci. Bull. 2019, 64, 418–421. [Google Scholar] [CrossRef]
  11. Sivapalan, M.; Takeuchi, K.; Franks, S.W.; Gupta, V.K.; Karambiri, H.; Lakshmi, V.; Liang, X.; McDonnell, J.J.; Mendiondo, E.M.; O’Connell, P.E.; et al. Decade on Predictions in Ungauged Basins (PUB), 2003–2012: Shaping an exciting future for the hydrological sciences. Hydrol. Sci. J. 2003, 48, 857–880. [Google Scholar] [CrossRef] [Green Version]
  12. Hou, J.; van Dijk, A.I.J.M.; Beck, H.E. Global satellite-based river gauging and the influence of river morphology on its application. Remote Sens. Environ. 2020, 239, 111629. [Google Scholar] [CrossRef]
  13. Coe, M.T.; Birkett, C.M. Calculation of river discharge and prediction of lake height from satellite radar altimetry: Ex-ample for the Lake Chad basin. Water Resour. Res. 2004, 40, W10205. [Google Scholar] [CrossRef]
  14. Papa, F.; Prigent, C.; Rossow, W.B. Monitoring flood and discharge variations in the large Siberian Rivers from a multi-satellite technique. Surv. Geophys. 2008, 29, 297–317. [Google Scholar] [CrossRef]
  15. Bjerklie, D.M.; Birkett, C.M.; Jones, J.W.; Carabajal, C.; Rover, J.A.; Fulton, J.W.; Garambois, P.-A. Satellite remote sensing estimation of river discharge: Application to the Yukon River Alaska. J. Hydrol. 2018, 561, 1000–1018. [Google Scholar] [CrossRef] [Green Version]
  16. Sichangi, A.W.; Wang, L.; Yang, K.; Chen, D.; Wang, Z.; Li, X.; Zhou, J.; Liu, W.; Kuria, D. Estimating continental river basin discharges using multiple remote sensing data sets. Remote Sens. Environ. 2016, 179, 36–53. [Google Scholar] [CrossRef] [Green Version]
  17. Gleason, C.J.; Smith, L.C. Toward global mapping of river discharge using satellite images and at-many-stations hy-draulic geometry. Proc. Natl. Acad. Sci. USA 2014, 111, 4788–4791. [Google Scholar] [CrossRef] [Green Version]
  18. Yang, S.; Wang, P.; Lou, H.; Wang, J.; Zhao, C.; Gong, T. Estimating river discharges in ungauged catchments using the slope–area method and unmanned aerial vehicle. Water 2019, 11, 2361. [Google Scholar] [CrossRef] [Green Version]
  19. Brombacher, J.; Reiche, J.; Dijksma, R.; Teuling, A.J. Near-daily discharge estimation in high latitudes from Sentinel-1 and 2: A case study for the Icelandic Þjórsá river. Remote Sens. Environ. 2020, 241, 111684. [Google Scholar] [CrossRef]
  20. Hirpa, F.A.; Hopson, T.M.; De Groeve, T.; Brakenridge, G.R.; Gebremichael, M.; Restrepo, P.J. Upstream satellite remote sensing for river discharge forecasting: Application to major rivers in South Asia. Remote Sens. Environ. 2013, 131, 140–151. [Google Scholar] [CrossRef]
  21. Sichangi, A.W.; Wang, L.; Hu, Z. Estimation of River Discharge Solely from Remote-Sensing Derived Data: An Initial Study Over the Yangtze River. Remote Sens. 2018, 10, 1385. [Google Scholar] [CrossRef] [Green Version]
  22. Ovakoglou, G.; Alexandridis, T.K.; Crisman, T.L.; Skoulikaris, C.; Vergos, G.S. Use of MODIS satellite images for detailed lake morphometry: Application to basins with large water level fluctuations. Int. J. Appl. Earth Obs. Geoinf. 2016, 51, 37–46. [Google Scholar] [CrossRef]
  23. Reil, A.; Skoulikaris, C.; Alexandridis, T.K.; Roub, R. Evaluation of riverbed representation methods for one-dimensional flood hydraulics model. J. Flood Risk Manag. 2018, 11, 169–179. [Google Scholar] [CrossRef] [Green Version]
  24. Brakenridge, G.R.; Nghiem, S.V.; Anderson, E.; Mic, R. Orbital microwave measurement of river discharge and ice status. Water Resour. Res. 2007, 43, W04405. [Google Scholar] [CrossRef]
  25. Robert Brakenridge, G.; Cohen, S.; Kettner, A.J.; De Groeve, T.; Nghiem, S.V.; Syvitski, J.P.M.; Fekete, B.M. Calibration of satellite measurements of river discharge using a global hydrology model. J. Hydrol. 2012, 475, 123–136. [Google Scholar] [CrossRef]
  26. Tarpanelli, A.; Brocca, L.; Lacava, T.; Melone, F.; Moramarco, T.; Faruolo, M.; Pergola, N.; Tramutoli, V. Toward the estimation of river discharge variations using MODIS data in ungauged basins. Remote Sens. Environ. 2013, 136, 47–55. [Google Scholar] [CrossRef]
  27. Revilla-Romero, B.; Beck, H.E.; Burek, P.; Salamon, P.; de Roo, A.; Thielen, J. Filling the gaps: Calibrating a rainfall-runoff model using satellite-derived surface water extent. Remote Sens. Environ. 2015, 171, 118–131. [Google Scholar] [CrossRef]
  28. Hou, J.; van Dijk, A.I.J.M.; Renzullo, L.J.; Vertessy, R.A. Using modelled discharge to develop satellite-based river gauging: A case study for the Amazon Basin. Hydrol. Earth Syst. Sci. 2018, 22, 6435–6448. [Google Scholar] [CrossRef] [Green Version]
  29. Li, H.; Li, H.; Wang, J.; Hao, X. Extending the ability of near-infrared images to monitor small river discharge on the Northeastern Tibetan Plateau. Water Resour. Res. 2019, 55, 8404–8421. [Google Scholar] [CrossRef]
  30. Shi, Z.; Chen, Y.; Liu, Q.; Huang, C. Discharge estimation using harmonized Landsat and Sentinel-2 product: Case studies in the Murray Darling Basin. Remote Sens. 2020, 12, 2810. [Google Scholar] [CrossRef]
  31. Durand, M.; Gleason, C.J.; Garambois, P.A.; Bjerklie, D.; Smith, L.C.; Roux, H.; Rodriguez, E.; Bates, P.D.; Pavelsky, T.M.; Monnier, J.; et al. An intercomparison of remote sensing river discharge estimation algorithms from measurements of river height, width, and slope. Water Resour. Res. 2016, 52, 4527–4549. [Google Scholar] [CrossRef] [Green Version]
  32. Brinkerhoff, C.B.; Gleason, C.J.; Feng, D.; Lin, P. Constraining remote river discharge estimation using reach-scale geomophology. Water Resour. Res. 2020, 56, e2020WR027949. [Google Scholar] [CrossRef]
  33. Huang, C.; Chen, Y.; Wu, J. Mapping spatio-temporal flood inundation dynamics at large river basin scale using time-series flow data and MODIS imagery. Int. J. Appl. Earth Obs. Geoinf. 2014, 26, 350–362. [Google Scholar] [CrossRef]
  34. Gorelick, N.; Hancher, M.; Dixon, M.; Ilyushchenko, S.; Thau, D.; Moore, R. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 2017, 202, 18–27. [Google Scholar] [CrossRef]
  35. Yamazaki, D.; Ikeshima, D.; Tawatari, R.; Yamaguchi, T.; O’Loughlin, F.; Neal, J.C.; Sampson, C.C.; Kanae, S.; Bates, P.D. A high-accuracy map of global terrain elevations. Geophys. Res. Lett. 2017, 44, 5844–5853. [Google Scholar] [CrossRef] [Green Version]
  36. Yamazaki, D.; Ikeshima, D.; Sosa, J.; Bates, P.D.; Allen, G.H.; Pavelsky, T.M. MERIT Hydro: A high-resolution global hydrography map based on latest topography dataset. Water Resour. Res. 2019, 55, 5053–5073. [Google Scholar] [CrossRef] [Green Version]
  37. Yamazaki, D.; O’Loughlin, F.; Trigg, M.A.; Miller, Z.F.; Pavelsky, T.M.; Bates, P.D. Development of the global width database for large rivers. Water Resour. Res. 2014, 50, 3467–3480. [Google Scholar] [CrossRef]
  38. Pekel, J.-F.; Cottam, A.; Gorelick, N.; Belward, A.S. High-resolution mapping of global surface water and its long-term changes. Nature 2016, 540, 418–422. [Google Scholar] [CrossRef]
  39. Tarpanelli, A.; Amarnath, G.; Brocca, L.; Massari, C.; Moramarco, T. Discharge estimation and forecasting by MODIS and altimetry data in Niger-Benue River. Remote Sens. Environ. 2017, 195, 96–106. [Google Scholar] [CrossRef]
  40. Brown, C.E. Coefficient of variation. In Applied Multivariate Statistics in Geohydrology and Related Sciences; Brown, C.E., Ed.; Springer: Berlin/Heidelberg, Germany, 1998; pp. 155–157. [Google Scholar]
  41. Warne, T.R. A Primer on multivariate analysis of variance (MANOVA) for behavioral scientists. Pract. Assess. Res. Eval. 2014, 19, 17. [Google Scholar] [CrossRef]
  42. Haase, R.F.; Ellis, M.V. Multivariate analysis of variance. J. Couns. Psychol. 1987, 34, 404. [Google Scholar] [CrossRef]
  43. Guo, K.; Zou, T.; Jiang, D.; Tang, C.; Zhang, H. Variability of Yellow River turbid plume detected with satellite remote sensing during water-sediment regulation. Cont. Shelf Res. 2017, 135, 74–85. [Google Scholar] [CrossRef]
Figure 1. (a) Locations of two study areas: (b) the Heihe River Basin and its selected sites (df); and (c) the Murray Darling Basin and its selected sites (gi).
Figure 1. (a) Locations of two study areas: (b) the Heihe River Basin and its selected sites (df); and (c) the Murray Darling Basin and its selected sites (gi).
Water 14 03854 g001
Figure 2. Observed discharge of the selected 7 hydrological gauges in the Heihe River Basin and the Murray Darling Basin.
Figure 2. Observed discharge of the selected 7 hydrological gauges in the Heihe River Basin and the Murray Darling Basin.
Water 14 03854 g002
Figure 3. River networks extracted from DEM. (a) The Heihe River Basin; (b) The Murray Darling Basin.
Figure 3. River networks extracted from DEM. (a) The Heihe River Basin; (b) The Murray Darling Basin.
Water 14 03854 g003
Figure 4. Pearson’s r correlation coefficient images of the C/M value of each pixel and observed discharge in time series. Red dots represent the finally selected location of the M pixel. (a) The Heihe River Basin; (b) The Murray Darling Basin.
Figure 4. Pearson’s r correlation coefficient images of the C/M value of each pixel and observed discharge in time series. Red dots represent the finally selected location of the M pixel. (a) The Heihe River Basin; (b) The Murray Darling Basin.
Water 14 03854 g004
Figure 5. The scatter plots of C/M series and observed discharge. (a) Sites of the Heihe River Basin; (b) Sites of the Murray Darling Basin.
Figure 5. The scatter plots of C/M series and observed discharge. (a) Sites of the Heihe River Basin; (b) Sites of the Murray Darling Basin.
Water 14 03854 g005
Figure 6. Estimated discharge and the observed discharge of each hydrological gauge for the validation period.
Figure 6. Estimated discharge and the observed discharge of each hydrological gauge for the validation period.
Water 14 03854 g006
Figure 7. Surface water dynamics of river sections in study regions. (a,c,e,g) are the water occurrence (the darker the color of water occurrence map, the easier it is to be extracted stably). (b,d,f,h) are the maximum water extent. (Gauge Qilian and Gauge Zhamashike are not shown due to data availability).
Figure 7. Surface water dynamics of river sections in study regions. (a,c,e,g) are the water occurrence (the darker the color of water occurrence map, the easier it is to be extracted stably). (b,d,f,h) are the maximum water extent. (Gauge Qilian and Gauge Zhamashike are not shown due to data availability).
Water 14 03854 g007
Figure 8. Selected ROIs near the Gauge 425008 as candidate virtual gauges, shown with (a) river width, (b) water occurrence, and (c) maximum water extent. The number labelled in the ROI stands for its ID.
Figure 8. Selected ROIs near the Gauge 425008 as candidate virtual gauges, shown with (a) river width, (b) water occurrence, and (c) maximum water extent. The number labelled in the ROI stands for its ID.
Water 14 03854 g008
Figure 9. (a) The selected 19 ROIs and the calculated Pearson’s r image within the whole study site. (b) Pearson’s r larger than 0.7. (c) Estimated discharge using the optimized ROI (ROI 5) and observed discharge for Gauge 425008.
Figure 9. (a) The selected 19 ROIs and the calculated Pearson’s r image within the whole study site. (b) Pearson’s r larger than 0.7. (c) Estimated discharge using the optimized ROI (ROI 5) and observed discharge for Gauge 425008.
Water 14 03854 g009
Table 1. Modelling period and validation period of each gauge, and the corresponding amount of Landsat images.
Table 1. Modelling period and validation period of each gauge, and the corresponding amount of Landsat images.
GaugeModeling PeriodAmountValidation PeriodAmount
Qilian7 January 2000–22 February 201112626 March 2011–2 December 201531
Zhamashike7 January 2000–26 March 201112711 April 2011–2 December 201532
Sunan25 January 2010–12 January 20174013 February 2017–30 December 201711
40903514 January 2000–20 July 2007 &
18 February 2010–26 April 2011
9721 August 2007–16 December 2009 &
1 September 2011–4 November 2011
24
41013014 January 2000–20 July 2007 &
18 February 2010–26 April 2011
9721 August 2007–16 December 2009 &
1 September 2011–4 November 2011
24
41420014 January 2000–20 July 2007 &
18 February 2010–26 April 2011
9721 August 2007–16 December 2009 &
1 September 2011–4 November 2011
24
4250085 January 2000–26 September 2006 &
24 January 2010–26 December 2010
8028 August 2007–20 October 2009 &
27 January 2011–11 November 2011
19
Table 2. The extraction results of gradient, sinuosity, and river width of each river section.
Table 2. The extraction results of gradient, sinuosity, and river width of each river section.
GaugeStudy SiteGradient (m·km−1)SinuosityRiver Width (m)
Qilian1-115.891.0864.09
1-216.711.2676.6
1-314.311.1054.99
1-48.171.0754.98
Zhamashike2-116.21.0668.52
2-239.071.1071.26
2-37.601.0155.82
2-420.461.0868.52
Sunan3-116.871.0023.4
3-230.091.1128.66
3-322.181.0063.67
3-418.871.0332.30
4090354-10.151.4162.27
4-20.481.3671.79
4-33.571.1160.92
4-41.401.0855.34
4101305-12.121.2841.15
5-20.321.6655.74
5-39.941.4753.02
5-43.871.2654.66
4142006-10.061.2278.65
6-20.071.60101.20
6-34.481.1151.78
6-41.151.1193.53
Table 3. The validation results of different hydrological gauges.
Table 3. The validation results of different hydrological gauges.
GaugeStudy SiteRRMSENSE
Qilian1-10.440.60
1-20.540.39
1-30.630.18
1-40.560.35
Zhamashike2-10.560.50
2-20.560.51
2-30.550.52
2-40.550.53
Sunan3-12.68−10.79
3-21.78−4.21
3-31.82−4.43
3-40.760.06
4090354-10.88−0.14
4-20.600.47
4-30.610.46
4-40.96−0.34
4101305-11.69−0.21
5-21.59−0.07
5-32.24−1.10
5-42.00−0.68
4142006-10.650.40
6-20.530.61
6-31.00−0.43
6-41.02−0.48
Table 4. F statistics and Sig values of the multivariate analysis of variance, and correlation coefficient (r) between RRMSE/NSE and river morphological indicators.
Table 4. F statistics and Sig values of the multivariate analysis of variance, and correlation coefficient (r) between RRMSE/NSE and river morphological indicators.
RRMSENSE
FSigrFSigr
gradient0.060.810.042.340.14−0.28
sinuosity0.890.360.120.580.450.26
width13.270.00−0.566.320.020.54
Table 5. Gradient, sinuosity, and width of the 19 ROIs.
Table 5. Gradient, sinuosity, and width of the 19 ROIs.
ROI No.Gradient (m·km−1)SinuosityRiver Width (m)
10.381.1174.69
21.121.0875.37
31.161.0666.38
40.901.0371.03
50.412.7489.28
60.661.0961.68
70.991.0884.96
80.581.1683.13
90.591.2781.62
100.601.1181.54
110.701.4768.05
120.661.1472.95
130.291.1978.99
140.471.8287.11
150.362.4765.44
160.981.0959.84
170.471.0757.70
180.381.0354.11
190.821.1351.65
Table 6. The accuracy of discharge estimation of 19 selected ROIs.
Table 6. The accuracy of discharge estimation of 19 selected ROIs.
ROI No.RRMSENSE
11.190.48
20.840.74
31.020.62
40.910.70
50.680.83
61.240.44
71.440.25
81.75−0.12
91.88−0.28
101.130.54
111.460.22
121.300.38
131.090.56
141.080.57
151.130.53
161.060.59
171.78−0.16
181.040.60
190.740.80
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shi, Z.; Chen, Q.; Huang, C. The Influence of River Morphology on the Remote Sensing Based Discharge Estimation: Implications for Satellite Virtual Gauge Establishment. Water 2022, 14, 3854. https://doi.org/10.3390/w14233854

AMA Style

Shi Z, Chen Q, Huang C. The Influence of River Morphology on the Remote Sensing Based Discharge Estimation: Implications for Satellite Virtual Gauge Establishment. Water. 2022; 14(23):3854. https://doi.org/10.3390/w14233854

Chicago/Turabian Style

Shi, Zhuolin, Qianqian Chen, and Chang Huang. 2022. "The Influence of River Morphology on the Remote Sensing Based Discharge Estimation: Implications for Satellite Virtual Gauge Establishment" Water 14, no. 23: 3854. https://doi.org/10.3390/w14233854

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop