On the Contribution of Satellite Altimetry-Derived Water Surface Elevation to Hydrodynamic Model Calibration in the Han River

Satellite altimetry can fill the spatial gaps of in-situ gauging networks especially in poorly gauged regions. Although at a generally low temporal resolution, satellite altimetry has been successfully used for water surface elevation (WSE) estimation and hydrodynamic modeling. This study aims to investigate the contribution of WSE from both short-repeat and geodetic altimetry to hydrodynamic model calibration, and also explore the contribution of the new Sentinel-3 mission. Two types of data sources (i.e., in-situ and satellite altimetry) are investigated together with two roughness cases (i.e., spatially variable and uniform roughness) for calibration of a hydrodynamic model (DHI MIKE 11) with available bathymetry. A 150 km long reach of Han River in China with rich altimetry and in-situ gauging data is selected as a case study. Results show that the performances of the model calibrated by satellite altimetry-derived datasets are acceptable in terms of Root Mean Square Error (RMSE) of simulated WSE. Sentinel-3A can support hydrodynamic model calibration even though it has a relatively low temporal resolution (27-day repeat cycle). The CryoSat-2 data with a higher spatial resolution (7.5 km at the Equator) are proved to be more valuable than the Sentinel-3A altimetry data with a low spatial resolution (104 km at the Equator) for hydrodynamic model calibration in terms of RMSE values of 0.16 and 0.18 m, respectively. Moreover, the spatially variable roughness can also improve the model performance compared to the uniform roughness case, with decreasing RMSE values by 2–14%. Our finding shows the value of satellite altimetry-derived datasets for hydrodynamic model calibration and therefore supports flood risk assessment and water resources management.


Introduction
Climate change and human activities have altered the river flow regimes and led to more frequent and severe extreme natural disasters, such as floods and droughts [1,2]. Monitoring hydraulic variables and relevant hydrological-hydraulic processes of rivers usually play an essential role in flood risk assessment and water resources management [3,4]. One of the key hydraulic variables is the water 2 of 18 surface elevation (WSE), which refers to the height of surface water above a given datum (i.e., 1985 National Height Datum of China for in-situ records, and EGM2008 datum for satellite WSE) [5]. Reliable and timely estimation of WSE is important for flood and water resources management tasks [6,7]. Traditionally, WSE measurements relied on in-situ gauging stations and often suffered from insufficient spatial coverage and data access constraints [8,9]. Hydraulic infrastructure and engineering projects are often designed in data sparsely or ungauged river reaches, which hampers the stakeholders to make decisions about river system management [10]. Hence, it is difficult for engineers to find an efficient way for WSE estimation along the entire river reach. Currently, a practical way to achieve this goal is the comprehensive estimation of WSE and other hydraulic variables on a whole river reach by using hydrodynamic models [4,11].
A hydrodynamic model is an efficient tool for flood forecasting and water resources management, as it can simulate the spatio-temporal dynamics of hydraulic variables given the boundary and initial conditions [12,13]. Observed data including WSE is often employed to calibrate the parameters of the hydrodynamic model and thus influences the model performance of simulated hydraulic variables [14]. However, observations from in-situ gauging stations are often sparse, and model parameters (e.g., spatially variable roughness parameters) probably cannot be estimated well from in-situ data through a calibration procedure [15]. The hydrodynamic model application is constrained in flood and water resources management, especially in ungauged or poorly gauged basins [16][17][18].
With the rapid advances in remote sensing sciences and technologies, space-borne sensors have provided valuable estimates of hydrological-hydraulic variables such as river extents and WSE [19][20][21][22]. WSE can be freely obtained from satellite altimetry missions (e.g., Envisat and Jason-2/-3), while river extents can be derived from optical or SAR (synthetic aperture radar) imagery satellite missions (e.g., Landsat and Sentinel-1/-2) [23,24]. These data from a satellite mission called satellite observations have been taken as a key supplement to in-situ data for hydrodynamic modeling. There have been many studies focusing on integrating satellite observations in facilitating hydrodynamic modeling. For example, satellite imagery observations have been adopted to calibrate two-dimensional (2-D) hydrodynamic models [25,26], to validate the efficiency of hydrodynamic models [27,28], and even integrated into model structures for data assimilation [29]. Although promising results have been shown by integrating satellite imagery observations in hydrodynamic modeling, there are also a lot of challenges. For example, optical imagery data are easily affected by cloud conditions while SAR imagery data always suffer the low spatial resolution [30], inundation extents are hard to be interpreted due to the noise of the backscatter and are also unavailable during the flood period.
WSE characterizes the surface water variations [31]. Satellite altimetry-derived WSE has been widely used to calibrate hydrodynamic models as an alternative to in-situ records, especially in ungauged areas [32,33]. For example, short-repeat missions such as Envisat and ERS-2 have been adopted in calibrating a quasi-two-dimensional hydrodynamic model in Italy [30]. Recently, CryoSat-2 was used to calibrate a one-dimensional hydrodynamic model in the Brahmaputra River [34].
Hydrodynamic models can benefit from satellite altimetry data in model calibration, improving the skills of simulating hydraulic variables such as WSE and discharge. However, the contributions of various altimetry-derived WSE datasets to hydrodynamic model calibration are different. The data-quality and sampling density are two potential pivotal factors influencing the contribution of satellite WSE on hydrodynamic model calibration. Previous works have revealed the importance of altimetry spatial sampling density to hydrodynamic model calibration through providing a high number of virtual river stage gauging stations from satellite altimetry [35,36], exploring the value of CryoSat-2 altimetry with dense sampling patterns [34], or applying various altimetry-derived datasets to hydrodynamic model calibration [37]. However, few studies have simultaneously evaluated the contributions of WSE from both short-repeat and geodetic altimetry missions on hydrodynamic model calibration, and also little is known about the contribution of the new Sentinel-3 mission [38].
Short-repeat missions such as Sentinel-3 have a repeat cycle of 10 to 35 days with a relatively sparse ground track pattern, while geodetic missions such as CryoSat-2 have dense point observations of the river longitudinal profile due to their drifting ground track pattern [39,40]. As a new generation satellite altimetry mission with a 27-day repeat cycle, Sentinel-3 is likely to deliver high-quality WSE observation for its new open-loop tracking and SAR modes [38]. The CryoSat-2 was launched in July 2010 to an orbit with a 369-day full repeat cycle and dense ground tracks [41]. Sentinel-3 is operating in a constellation of two-satellites (Sentinel-3A and Sentinel-3B). The Sentinel-3A (S3A) and Sentinel-3B (S3B) were launched in February 2016 and April 2018, respectively. Detailed information on Sentinel-3 and CryoSat-2 is available at https://earth.esa.int/web/guest/missions/esa-operational-eo-missions/.
This study aims to investigate the contribution of altimetry-derived WSE from short-repeat and geodetic missions (Sentinel-3A and CryoSat-2) to hydrodynamic model calibration. The 150 km long reach of the Han River with rich data in China was selected as a case study area. Model calibration and validation are carried out under fourteen scenarios by combining different data sources (i.e., in-situ and satellite altimetry) and two roughness cases (i.e., spatially variable and uniform roughness). Therefore, the results of our study will help to implement satellite altimetry-derived WSE into hydrodynamic model calibration and therefore support flood forecasting and water resources management.

Study Area
The Han River offers a large amount of biological and hydropower resources and serves as one of the major water sources for the South-to-North Water Diversion in China [42]. Concerning the availability of in-situ hydrometric and satellite altimetry data, a 150 km long reach of Han River in China ( Figure 1) is selected to explore the contributions of satellite altimetry-derived datasets on hydrodynamic model calibration. The case area belongs to the middle stream of the Han River between Huangjiagang and Yicheng hydrometric stations. Three main tributaries are located in the study area, and two hydraulic projects, the Wangbuzhou and Cuijiaying reservoirs, are also located with the functions of flood control and hydropower production ( Figure 1).

Hydrometric Data
Daily discharge of four hydrometric stations and daily water level of five hydrometric stations from 01/01/2016 to 12/31/2018 are collected from the Bureau of Hydrology, Changjiang Water Resources Commission, affiliated to the Ministry of Water Resources of China (Table 1). The main-channel geometry of the reach is described by the 94 surveyed cross-section profiles, which are surveyed by the Bureau of Hydrology in 2016. The information on the two hydraulic projects can be found in Table 1 and Figure 2. The operational rules applied for the two projects are dependent on their water levels and water inflows. Specifically, the outflow equals the inflow when the inflow is lower than the maximum allowable release under the normal water level, otherwise, the outflow equals the maximum allowable release under its corresponding water level (Figure 2c,d).

Hydrometric Data
Daily discharge of four hydrometric stations and daily water level of five hydrometric stations from 01/01/2016 to 12/31/2018 are collected from the Bureau of Hydrology, Changjiang Water Resources Commission, affiliated to the Ministry of Water Resources of China (Table 1). The mainchannel geometry of the reach is described by the 94 surveyed cross-section profiles, which are surveyed by the Bureau of Hydrology in 2016. The information on the two hydraulic projects can be found in Table 1 and Figure 2. The operational rules applied for the two projects are dependent on their water levels and water inflows. Specifically, the outflow equals the inflow when the inflow is lower than the maximum allowable release under the normal water level, otherwise, the outflow equals the maximum allowable release under its corresponding water level (Figure 2c,d). Note: Chainage refers to the distance along the river centerline between points of interest. River width indicates the bank's full river width from surveyed cross-section river profiles.

Satellite Altimetry Data
The altimetry-derived WSE are retrieved from Sentinel-3A and CryoSat-2. Sentinel-3A data are collected from the ESA level 2 product "Enhanced measurements", which are downloaded from the Copernicus Open Access Hub (https://scihub.copernicus.eu/dhus/). The ESA level 2 product contains 20Hz measurements including the waveform, tracker range, satellite altitude, and corrections. The waveforms are retracked with a traditional offset center of gravity retracker [7]. For CryoSat-2, we use the ESA level 1b product from baseline C (https://science-pds.cryosat.esa.int/). The ESA level 1b product contains 20Hz measurements including the waveform, tracker range, satellite altitude, and corrections. The CryoSat-2 waveforms were retracked with a narrow primary peak threshold retracker [43].
The water surface elevations WSE, are constructed via the following equation where h is the altitude of the satellite, and N is the geoid height, which is the EGM2008 geoid model [44]. R is the range that is the distance from the satellite to the surface, which can be expressed as where trac R is the distance measured by the onboard tracker, hence the distance to the nominal bin in the waveform and retrac R is the re-tracking correction. atm R is the atmospheric corrections

Satellite Altimetry Data
The altimetry-derived WSE are retrieved from Sentinel-3A and CryoSat-2. Sentinel-3A data are collected from the ESA level 2 product "Enhanced measurements", which are downloaded from the Copernicus Open Access Hub (https://scihub.copernicus.eu/dhus/). The ESA level 2 product contains 20Hz measurements including the waveform, tracker range, satellite altitude, and corrections. The waveforms are retracked with a traditional offset center of gravity retracker [7]. For CryoSat-2, we use the ESA level 1b product from baseline C (https://science-pds.cryosat.esa.int/). The ESA level 1b product contains 20Hz measurements including the waveform, tracker range, satellite altitude, and corrections. The CryoSat-2 waveforms were retracked with a narrow primary peak threshold retracker [43].
The water surface elevations WSE, are constructed via the following equation where h is the altitude of the satellite, and N is the geoid height, which is the EGM2008 geoid model [44]. R is the range that is the distance from the satellite to the surface, which can be expressed as where R trac is the distance measured by the onboard tracker, hence the distance to the nominal bin in the waveform and R retrac is the re-tracking correction. R atm is the atmospheric corrections including the wet and dry tropospheric corrections and the ionospheric correction. R geo is the geophysical corrections including solid earth tide, pole tide, and ocean loading tide. All corrections are taken from the altimetry products.

WSE Data Processing
The altimetry-derived WSE are extracted by a water mask from Global Surface Water Explorer (https://global-surface-water.appspot.com/). Specifically, the maximum extent is used to pre-select as many observations as possible [45]. The identical satellite ground tracks with a repeat cycle of 27 days allow to derive the WSE time series at two virtual stations (VSs) from Sentinel-3A over the study reach ( Figure 1). Firstly, the outliers are discarded by comparing the altimetry-derived WSE with ASTER GDEM v.3 (Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model) [46] with an absolute bias threshold of 20 m. The DEM datasets have been reprojected to the EGM2008 geoid at the platform of VDatum software [47]. Secondly, a filtering process is performed for each track to discard outliers based on the work by Jiang et al. [40] and finally, the median value of all WSE observations for a given track is used to derive WSE time series. The data processing for CryoSat-2 altimetry is slightly different from that of Sentinel-3A due to its drifting ground track pattern. A series of virtual stations were selected and all observations were relocated to the nearest virtual stations based on their local slopes [43]. The virtual stations are determined based on the coverages of CryoSat-2 observations along the river reach and placed on the model river line (i.e., the centerline of the river from the extracted water mask). The procedure of relocating CryoSat-2 observations onto virtual stations has been proved to be an efficient way of constructing WSE time series at VSs and strengthening the hydrologic and hydrodynamic applications of CryoSat-2 [43]. Balancing the maximization of the VS numbers and the minimization of error from projecting observations to virtual stations, 16 VSs have been figured out in the study reach (Table 2). An outlier filtering algorithm is performed to obtain the median values of each track: a measurement is discarded if |h-u|>3δ (h is the point observation, u is the mean value and δ is the standard deviation) [40]. All the median observations are relocated to the nearest virtual stations through local slope correction. It should be noted that there are often systematic biases between satellite altimetry-derived observations and in-situ data due to the unknown local vertical datums and other factors. If their mean values of every series are removed, both in-situ data and satellite altimetry-derived datasets are transformed to WSE anomalies.

Hydrodynamic Model in the Study Area
A one-dimension (1D) hydrodynamic model is used to simulate flow regimes and water levels of the study reach. The model was running at the DHI MIKE 11 software [48]. The software solves the Saint-Venant equations with an implicit 6-point finite-difference scheme. The model is directly forced by the observed daily discharge data from the main upstream and three tributaries, while the daily water level of Yicheng hydrometric station is served as the downstream boundary condition (shown in Figure 1). The study reach is sketched by 94 cross-section profiles, and the profiles are derived from the detailed bathymetric map with an interval of around 50 m and thus can capture the variations of riverbed geometry. The operational rules of Wangbuzhou and Cuijiaying reservoirs are also simulated in the model as it is described in Section 2.2.

Hydrodynamic Model Calibration
The hydrodynamic model calibration focuses on the estimation of the channel roughness, i.e., the Strickler coefficient Ks (the reciprocal of the Manning coefficient, see in Manning's Equation (3)).
With given bathymetry data, the cross-section area (A), the hydraulic radius (r), the slope of the hydraulic gradient line (S) at the discharge (Q) can be described as a function of water depth relating to WSE. Only Ks is a parameter that varies along the river and is often calibrated through maximizing the degree of goodness of fit between observed and simulated WSE. To figure out the effect of heterogeneity of roughness distribution on model performance, two roughness cases are assumed in the calibration process. One is that only a single Ks coefficient is assigned to the whole reach (uniform Ks). The other one is that the river reach is partitioned into sub-reaches according to the hydro-morphological characteristics of the river and the distribution of hydrometric stations (variable Ks). The variable Ks with sub-reaches has also been proved to be a feasible approach to improve the hydrodynamic model performance [34]. Specifically, four sub-reaches are used in our case study for the variable Ks. The first 29 km long reach is very gentle (0.035 m/km) and the following 44 km long reach is narrow and steep with a complex river morphology. In the following 38 km long reach, the river channel is steeper with a gradient of 0.125 m/km. The last reach is a little wider and deeper.
To investigate the efficiency and reliability of the various satellite altimetry-derived datasets, three configurations are constructed for the hydrodynamic model calibration. The first one is from the in-situ data (called Configuration A), the second one is from the altimetry-derived data from Sentinel-3A and CryoSat-2 (called Configuration B) and the last one is from the combination of in-situ and satellite altimetry-derived data (called Configuration C). Combining the three configurations and two roughness cases, there are fourteen scenarios for calibrating the hydrodynamic model as shown in Table 3. Table 3. Summary of hydrodynamic model calibration.

Scenarios Design
Using In-Situ Using Sentinel-3A

CryoSat-2
Configuration A: using in-situ data (two schemes) BV3 -Y Y Configuration C: using both in-situ and satellite altimetry data (six schemes) The 1D hydrodynamic model is calibrated by the fourteen scenarios from 01/01/2016 to 12/31/2017 through minimizing the Root Mean Square Error (RMSE) between simulated and observed WSE. The model is also validated by the data from 01/01/2018 to 12/31/2018. The manual trial-and-error method is adopted for testing the variable Ks case and the scenario-based method is used (running models with a range of values of Ks at an interval of 0.5/0.1 m 1/3 s −1 ) for the uniform Ks case [33,49].
where O i andp i are the observed and simulated WSE at the ith point, respectively. n is the total number of WSE observations.

Calibrated Strickler Coefficient Ks
The calibrated uniform Ks plays a vital role in hydrodynamic simulation, and the optimal Ks value can be easily determined from the clear patterns shown in Figure 3a (e.g., 34.6 m 1/3 /s as the optimal value for the scheme AU1). The model performance is sensitive to uniform Ks value (RMSE ranging from 0.105-0.450 m) and the optimal RMSE values are assigned near the high bound of the parameter range. The calibrated Ks values of four sub-reaches reflect the characteristics of each sub-reach (Figure 3b). The Ks values of the second and the final sub-reaches range from 27-33 m 1/3 /s. The Ks is relatively high at the first upstream 29-km reach (approximately 37 m 1/3 /s) due to the flat and straight river channel. On the contrary, the Ks value has decreased significantly in the third reach due to the increasing water depth and highly vegetated channel. Roughness is closely related to flow resistance, which is positively influenced by river depth and vegetable conditions. More specifically, flow resistance increases with a highly vegetated channel and river depth, and leads to higher roughness (smaller Ks) [49,50]. It is interesting to find that the RMSE values jump near the Strickler coefficient close to 30 in Figure 3a. As there are discharge bursts in 2017 that can challenge the calibration of parameters (Figure 2a), the RMSE value is sensitive to the bias of peak simulation. The peaks of discharge should be processed separately and the model parameters (Ks) should be chosen in a special way. There is not enough available data (in-situ records and satellite observations) to separately calibrate the parameters. In addition, all the jumping curves are from the uniform Ks scenarios, which proves that the uniform parameter can bring more uncertainties of discharge simulation and cut down the contribution of the satellite altimetry-derived WSE to hydrodynamic model calibration at the whole river reach.

Model Calibration with the in-situ Observations
Regarding configuration A, hydrodynamic models are calibrated by the in-situ observations using the uniform and variable Ks (i.e., the schemes AU1 and AV1). The results during the calibration and validation periods reveal that the model can reproduce WSE well, where the RMSE ranging from 0.100-0.216 m (Table 4). However, there are few obvious differences in the accuracies of WSE simulation from the two roughness cases. The averaged performances from the variable Ks are slightly better than those from the uniform Ks during the calibration periods while the results are not as good as those from the uniform Ks during the validation periods. Nevertheless, the improved results from the variable Ks can be found with lower RMSE values at the specific station (Figures 4 and 5). For example, the RMSE from the variable Ks is 0.074 m at Xiangyang station, decreased by 11.5% compared to that from the uniform Ks ( Figure 5). Overall, the improvements with the variable Ks highlight the potential superiority of the variable Ks in hydrodynamic model calibration.

Model Calibration with the Satellite Altimetry-derived Observations
Concerning configuration B, the model is calibrated by the satellite altimetry-derived WSE. There are six model scenarios, i.e., calibration by Sentinel-3A, CryoSat-2, and the combinations of these two datasets taking the uniform and variable Ks cases (see the details in Table 3). The calibrated models show the comparable performance in comparison to those by using in-situ data in terms of RMSE values during both calibration and validation periods where the RMSE ranges 0.153-0.179 m

Model Calibration with the Satellite Altimetry-derived Observations
Concerning configuration B, the model is calibrated by the satellite altimetry-derived WSE. There are six model scenarios, i.e., calibration by Sentinel-3A, CryoSat-2, and the combinations of these two datasets taking the uniform and variable Ks cases (see the details in Table 3). The calibrated models show the comparable performance in comparison to those by using in-situ data in terms of RMSE values during both calibration and validation periods where the RMSE ranges 0.153-0.179 m in Table 4. Thus, the altimetry-derived datasets from Sentinel-3A and CryoSat-2 can be applied in the hydrodynamic model calibration. However, the contributions of WSE from different altimetry to hydrodynamic model calibration are different. The schemes BU(V)2 (against CryoSat-2) show the lowest values of RMSE while the schemes BU(V)1 (against Sentinel-3A) perform worst during the calibration periods. Interestingly, the model calibrated by the schemes BU(V)2 shows the best performance in the calibration period while it performs worst in the validation period probably due to their different sample size during the calibration and validation periods. The simulated WSE at virtual stations from the variable Ks consistently shows a better performance than those from the uniform Ks. The distributions of difference between model simulations and satellite altimetry-derived WSE have also been shown in Figures 6 and 7, which further proves the advantages of the variable Ks.
performance in the calibration period while it performs worst in the validation period probably due to their different sample size during the calibration and validation periods. The simulated WSE at virtual stations from the variable Ks consistently shows a better performance than those from the uniform Ks. The distributions of difference between model simulations and satellite altimetryderived WSE have also been shown in Figures 6 and 7, which further proves the advantages of the variable Ks.   virtual stations from the variable Ks consistently shows a better performance than those from the uniform Ks. The distributions of difference between model simulations and satellite altimetryderived WSE have also been shown in Figures 6 and 7, which further proves the advantages of the variable Ks.   To show the performances of hydrodynamic model calibrated by altimetry-derived datasets at three in-situ stations, their statistical metrics (RMSE) are shown in Figure 5 and their hydrographs are depicted in Figure 8. The models calibrated by satellite altimetry-derived datasets can also provide reliable WSE estimations at these in-situ stations. Moreover, the statistical metrics are different for every station due to their flow regimes and channel geometry. For example, WSE estimations at Miaogang station are not as good as those at other two in-situ stations with higher RMSE values ( Figure 5). The model performance can also be dependent on the schemes with different satellite altimetry-derived datasets or roughness cases. The schemes BU(V)2 (against CryoSat-2) show better skills than the schemes BU(V)1 (against Sentinel-3A) or BU(V)3 (against all satellite altimetry-derived datasets) in simulating WSE at in-situ stations ( Figure 5). For example, the RMSE values from the scheme BU2 at Miaogang station are 0.158/0.316 m for the calibration/validation periods, while the RMSE values from the schemes BU1 and BU3 are 0.162/0.420 and 0.162/0.306 m, respectively. RMSE values ( Figure 5). The model performance can also be dependent on the schemes with different satellite altimetry-derived datasets or roughness cases. The schemes BU(V)2 (against CryoSat-2) show better skills than the schemes BU(V)1 (against Sentinel-3A) or BU(V)3 (against all satellite altimetryderived datasets) in simulating WSE at in-situ stations ( Figure 5). For example, the RMSE values from the scheme BU2 at Miaogang station are 0.158/0.316 m for the calibration/validation periods, while the RMSE values from the schemes BU1 and BU3 are 0.162/0.420 and 0.162/0.306 m, respectively.

Model Calibration with Both in-situ and Satellite Altimetry-derived Observations
In the case of configuration C, the hydrodynamic model was calibrated by all available observations at in-situ and virtual stations. The contributions of additional satellite altimetry-derived WSE to hydrodynamic model calibration are relatively small at the evaluated in-situ and virtual stations (Table 4). This can be attributed to the fact that the temporal resolution of the available insitu observations (daily frequency) is higher than that of satellite altimetry-derived datasets (e.g., 27day frequency for S3A observations). Despite the little contribution, the advantages of additional satellite altimetry-derived datasets to the hydrodynamic model can still be found through the

Model Calibration with Both In-Situ and Satellite Altimetry-Derived Observations
In the case of configuration C, the hydrodynamic model was calibrated by all available observations at in-situ and virtual stations. The contributions of additional satellite altimetry-derived WSE to hydrodynamic model calibration are relatively small at the evaluated in-situ and virtual stations (Table 4). This can be attributed to the fact that the temporal resolution of the available in-situ observations (daily frequency) is higher than that of satellite altimetry-derived datasets (e.g., 27-day frequency for S3A observations). Despite the little contribution, the advantages of additional satellite altimetry-derived datasets to the hydrodynamic model can still be found through the reproduction of the flood event (i.e., the flow processes with the maximum discharge at the upper boundary) (Figure 9). The WSE obtained from the model by calibration using all available observations is different from that obtained by calibration using only in-situ observations in the region where the satellite altimetry-derived observations can be found (shown in Figure 9d). With more spatially-distributed WSE observations, the hydrodynamic model provides more reliable WSE estimations on the whole river reach. The integration of altimetry-derived datasets and in-situ observations has led to the best performance for the hydrodynamic model calibration. Therefore, the satellite altimetry-derived datasets can be taken as a significant complementary to in-situ gauging networks, especially for improving the reliable simulation of WSE at the whole river reach scale. region where the satellite altimetry-derived observations can be found (shown in Figure 9d). With more spatially-distributed WSE observations, the hydrodynamic model provides more reliable WSE estimations on the whole river reach. The integration of altimetry-derived datasets and in-situ observations has led to the best performance for the hydrodynamic model calibration. Therefore, the satellite altimetry-derived datasets can be taken as a significant complementary to in-situ gauging networks, especially for improving the reliable simulation of WSE at the whole river reach scale.

Effects of Satellite Altimetry-Derived WSE on Hydrodynamic Model Calibration
Two pivotal factors are influencing the performance of the hydrodynamic model calibrated by satellite altimetry-derived WSE. One is the accuracy of satellite data itself and the other one is the spatio-temporal resolution of satellite altimetry missions. The two VSs of Sentinel-3A (i.e., VS1 and VS2) are close to the Huangjiagang and Xiangyang hydrometric stations (eight km from VS1 to Huangjiagang, four km from VS2 to Xiangyang). The slopes between the two VSs and their corresponding hydrometric stations are too flat to impact the WSE. The WSE at the two VSs can replace the data observed at the Huangjiagang and Xiangyang stations. The RMSE values are 0.28 m and 0.16 m for the VS1 and VS2, respectively. In the same way, the accuracies of CryoSat-2 at VS3 and VS7 are also evaluated by the observation data from Lianghekou (one km from VS3) and Xiangyang (eight km from VS7). The RMSE values are 0.14 m and 0.71 m for the VS3 and VS7, respectively. According to the previous studies [7,51], the RMSE values for satellite altimetry below 0.3 m are taken as good, and below 0.6 m are moderate. Sentinel-3A and CryoSat-2 data are not only found to be good and moderate in the literature [7,34,43] but also have been proved to be good and moderate in our case study river, a medium river (around 400 m) with braided and single-threaded reaches. By discarding outliers based on DEM approaches, a good or moderate accuracy of altimetry observations can be obtained in the case study, indicating that the land relief issue can be well

Effects of Satellite Altimetry-Derived WSE on Hydrodynamic Model Calibration
Two pivotal factors are influencing the performance of the hydrodynamic model calibrated by satellite altimetry-derived WSE. One is the accuracy of satellite data itself and the other one is the spatio-temporal resolution of satellite altimetry missions. The two VSs of Sentinel-3A (i.e., VS1 and VS2) are close to the Huangjiagang and Xiangyang hydrometric stations (eight km from VS1 to Huangjiagang, four km from VS2 to Xiangyang). The slopes between the two VSs and their corresponding hydrometric stations are too flat to impact the WSE. The WSE at the two VSs can replace the data observed at the Huangjiagang and Xiangyang stations. The RMSE values are 0.28 m and 0.16 m for the VS1 and VS2, respectively. In the same way, the accuracies of CryoSat-2 at VS3 and VS7 are also evaluated by the observation data from Lianghekou (one km from VS3) and Xiangyang (eight km from VS7). The RMSE values are 0.14 m and 0.71 m for the VS3 and VS7, respectively. According to the previous studies [7,51], the RMSE values for satellite altimetry below 0.3 m are taken as good, and below 0.6 m are moderate. Sentinel-3A and CryoSat-2 data are not only found to be good and moderate in the literature [7,34,43] but also have been proved to be good and moderate in our case study river, a medium river (around 400 m) with braided and single-threaded reaches. By discarding outliers based on DEM approaches, a good or moderate accuracy of altimetry observations can be obtained in the case study, indicating that the land relief issue can be well addressed by discarding invalid measurements (e.g., a careful analysis of satellite waveforms identification). Our findings suggest that a good or moderate accuracy of altimetry observations can provide satisfactory results for hydrodynamic model calibration.
The spatio-temporal resolution of satellite altimetry missions directly influences the performance of hydrodynamic model calibration. Our findings further revealed that CryoSat-2 geodetic altimetry with high spatial resolution shows better performance in hydrodynamic model calibration than Sentinel-3A altimetry with a low spatial resolution (Table 4 and Figures 3-5). It should be noted that CryoSat-2 geodetic altimetry has a full cycle of 369 days while Sentinel-3A has a 27-day cycle at the specified virtual stations. The values of the new Sentinel-3A altimetry to hydrodynamic model calibration are highlighted due to its high data-quality with a new on-board tracking mode and easier work in deriving WSE time series at virtual stations. But it is challenging to process the spatially-distributed CryoSat-2 observations and use them for hydrodynamic model calibration [34,39]. In our case study, the CryoSat-2 data are relocated onto a series of virtual stations and then transformed into anomalies for hydrodynamic model calibration [43]. The systematic bias still remains when processing such few CryoSat-2 observations into anomalies, and can be propagated into WSE simulation of the model and reduces its accuracy. This could be attributed to the fact that the results are not the best among the schemes of configuration B when pooling the Sentinel-3A and CryoSat-2 datasets for model calibration (shown in Table 4 and Figures 3-5). Further works can focus on producing blending products with more densified data from different satellite altimetry missions through the statistical and hydraulic blending approaches and use them for hydrodynamic modeling [52,53].

Comparison of Model Performance from Uniform and Variable Roughness Parameters
The effects of the two parameterizations on hydrodynamic model calibration have been reported in the literature [30,34], and are further explored on the use of various altimetry-derived datasets in our case study. The uniform Ks is widely used in hydrodynamic model calibration for its simplicity and promising efficiency at the evaluated stations. However, this common approach fails to characterize the variations of the riverbed and is probably unable to reproduce accurate WSE at other locations of the river under high heterogeneity [4,54]. The hydrodynamic model with uniform Ks calibrated by Sentinel-3A derived-WSE shows poor skills with high RMSE values at in-situ stations ( Figure 5). Moreover, the uniform Ks can also be attributed to the equifinality issue, i.e., a model with different roughness sets performs equally well at the evaluated gauging stations (shown in Figure 3a). To find a more efficient approach for more accurate WSE estimation along the river under the heterogeneity of the channel, the variable Ks is recommended for hydrodynamic model calibration as it takes the spatial variability of flow resistance into account [34]. The model with the variable Ks can simulate WSE more accurately with lower RMSE values compared to that with the uniform Ks.
Two pivotal factors are influencing the hydrodynamic model calibration with the variable Ks. One is the availability of data for model calibration and the other is the density of resolved variable roughness parameters. With satellite altimetry-derived data available, the variable Ks could be an efficient approach to facilitate hydrodynamic modeling. The model performance is also closely dependent on the satellite altimetry datasets with different spatio-temporal resolution. The CryoSat-2 geodetic mission with dense observations is more helpful than those with a low spatial resolution (Sentinel-3A) in hydrodynamic model calibration (Table 4 and Figures 3-5). The contributions of the densified roughness parameters on model calibration have been highlighted in the literature [4,34], while the improvements of the densified roughness parameters come at the expense of uncertainties and computation time. The more parameters commonly required more data for calibration [4]. Schneider et al. [34] used the 10-km variable roughness parameters to facilitate hydrodynamic modeling and found that the performances of the model were not refined well with large uncertainties. This issue is further revealed in our case study that the accuracies of WSE from the variable Ks are slightly better than those from the uniform Ks during the calibration periods with reducing RMSE values by 2-14%, while the results are not as good as those from the uniform Ks during the validation periods. However, the advantages of the variable Ks are not only to improve the accuracy of simulated WSE at the evaluated gauging stations but also to improve the hydrodynamic simulation on the entire river reach via a physically-based strategy. With more satellite altimetry-derived WSE available, the variable Ks could be estimated well and therefore increases the performance of the hydrodynamic model [55]. As the calibrated Ks values are reasonable compared to previous works [34,56], the variable Ks is recommended for hydrodynamic model calibration.

Conclusions
The contributions of current satellite altimetry missions (i.e., short-repeat and geodetic altimetry) on hydrodynamic model calibration are reviewed. The contribution of the new Sentinel-3A altimetry is also assessed. Hydrodynamic model (DHI MIKE 11) calibrations with available bathymetry data are carried out by combining seven calibration datasets (i.e., in-situ, satellite altimetry, and both of them) and two roughness cases (i.e., spatially variable roughness and uniform roughness). In the case study river reach with a river width of around 400 m, both Sentinel-3A and CryoSat-2 can deliver useful observations with the averaged RMSE values of 0.22 and 0.49 m, respectively. Our findings suggest that a good or moderate accuracy of altimetry observations (RMSE below 0.6 m) can contribute to the hydrodynamic model calibration. Specifically, calibration against altimetry-derived WSE achieves promising results in terms of RMSE as 0.100-0.216 m during the calibration periods, which are in line with previous studies [30,34]. Moreover, model performance varies with different altimetry-derived datasets due to their different characteristics such as spatio-temporal sampling patterns.
Although the satellite altimetry has a low temporal resolution, both Sentinel-3A and CryoSat-2 can provide a comparable result as those from the daily in-situ observations. CryoSat-2 with a higher spatial resolution (7.5 km at the Equator) brings better performances than short-repeat Sentinel-3A altimetry with a low resolution (104 km at the Equator) in the hydrodynamic model calibration according to their RMSE values of 0.16 and 0.18 m, respectively. Moreover, the variable Ks improves the performance of hydrodynamic simulation compared to the uniform Ks reducing the RMSE values by 2-14%, and is recommended in hydrodynamic model calibration for its physically-based strategy. Our study not only enriches the existing knowledge on integrating satellite altimetry-derived datasets into hydrodynamic model calibration but also supports the future hydrology-related satellite altimetry design.

Conflicts of Interest:
The authors declare no conflict of interest.