Estimating Reservoir Release Using Multi-Source Satellite Datasets and Hydrological Modeling Techniques

Reservoir release is an essential variable as it affects hydrological processes and water availability downstream. This study aims to estimate reservoir release using a satellite-based approach, specially focusing on the impacts of inflow simulations and reservoir water storage change (RWSC) on release estimates. Ten inflow simulations based on hydrological models and blending schemes are used in combination with three RWSC estimates based on two satellite-based approaches. A case study is performed at the Ankang reservoir, China. The results demonstrate that release estimates show high skill, with normalized root-mean-square error (NRMSE) less than 0.12 and Kling-Gupta Efficiency (KGE) over 0.65. The performance of release estimates is varying with and influenced by inflow simulations and RWSC estimates, with NRMSE ranging from 0.09–0.12 and KGE from 0.65–0.74. Based on time-varying Bayesian Model Averaging (BMA) approaches and synthetic aperture radar (SAR) satellite datasets, more accurate inflow and RWSC estimates can be obtained, thus facilitating substantially release estimates. With multi-source satellite datasets, temporal scale of reservoir estimates is increased (monthly and bi-weekly), acting as a key supplement to in situ records. Overall, this study explores the possibility to reconstruct and facilitate reservoir release estimates in poorly gauged dammed basins using hydrological modeling techniques and multi-source


Introduction
Reservoirs are essential tools in water resources management by providing immense societal and economic benefits in forms of hydroelectricity generation, flood control and water resources diversion [1,2]. The boom of reservoir impoundment has been underway since last century and will continue for the following decades, driven by climate change adaptation and population growth [3,4]. However, dam construction and regulation are disastrous for ecosystems and their negative impacts are well documented [5,6]. Great concerns have been raised about these impacts, particularly in transboundary river basins [7,8]. Therefore, reservoir monitoring is crucial for the well-being of regions with massive riparian population.
Reservoir release is an essential variable as it affects hydrological processes and water availability downstream [9,10]. Estimation of reservoir release is of great significance for In light of the above, this study aims to reconstruct reservoir release by combining multi-source satellite datasets and hydrological modeling techniques, focusing on the impacts of satellite dataset and inflow uncertainty on release simulations. To address such concerns, ten inflow simulations based on hydrological models and blending schemes are used in combination with three RWSC estimates based on two satellite-based approaches. The experiments are carried out at Ankang reservoir, located in the Han River of Central China. It is noteworthy that our case study is not only to improve the accuracies of individual estimations (i.e., inflow and RWSC), but more importantly to understand the impacts of propagation of associated uncertainties on reservoir release. The following two research questions are to be addressed: (1) What is the impact of different hydrological modeling and satellite-based approaches on the simulation of inflow and RWSC combined with multi-source satellite datasets?
(2) How do inflow simulations and RWSC estimates affect the accuracy of release estimates?

Study Area
The Ankang reservoir is located on the upstream of the Han River in China ( Figure 1). The basin is dominated by a subtropical monsoon climate with a majority of rainfall from June to October. The upstream has a drainage area of 38,640 km 2 with mean annual inflow of 568.78 m 3 /s [40]. We selected reservoir as a case study considering two aspects. First, concerning the availability of both in situ and satellite datasets, such a medium reservoir (2.58 km 3 in water capacity and 70 km 2 in water extent at the normal water level; see details in Table 1) [41] is suitable to investigate the effectiveness of our satellite-based approach, which has potential of being transferred to estimate release in ungauged basins with the same or larger reservoirs. Second, from a regional perspective, the release directly contributes to the inflow of Danjiangkou reservoir, thus, the case study can serve a basis for water resources management and provide a cooperative environment for all stakeholders of the Han River basin, which serves as the major water source of the South-to-North Water Diversion [42]. with a single hydrological model [29,37], while calculating RWSC using one approach and failing to exploit the contributions of satellite datasets (WSE and SWE) on RWSC [33]. In light of the above, this study aims to reconstruct reservoir release by combining multi-source satellite datasets and hydrological modeling techniques, focusing on the impacts of satellite dataset and inflow uncertainty on release simulations. To address such concerns, ten inflow simulations based on hydrological models and blending schemes are used in combination with three RWSC estimates based on two satellite-based approaches. The experiments are carried out at Ankang reservoir, located in the Han River of Central China. It is noteworthy that our case study is not only to improve the accuracies of individual estimations (i.e., inflow and RWSC), but more importantly to understand the impacts of propagation of associated uncertainties on reservoir release. The following two research questions are to be addressed: (1) What is the impact of different hydrological modeling and satellite-based approaches on the simulation of inflow and RWSC combined with multi-source satellite datasets?
(2) How do inflow simulations and RWSC estimates affect the accuracy of release estimates?

Study Area
The Ankang reservoir is located on the upstream of the Han River in China ( Figure  1). The basin is dominated by a subtropical monsoon climate with a majority of rainfall from June to October. The upstream has a drainage area of 38,640 km 2 with mean annual inflow of 568.78 m 3 /s [40]. We selected reservoir as a case study considering two aspects. First, concerning the availability of both in situ and satellite datasets, such a medium reservoir (2.58 km 3 in water capacity and 70 km 2 in water extent at the normal water level; see details in Table 1) [41] is suitable to investigate the effectiveness of our satellite-based approach, which has potential of being transferred to estimate release in ungauged basins with the same or larger reservoirs. Second, from a regional perspective, the release directly contributes to the inflow of Danjiangkou reservoir, thus, the case study can serve a basis for water resources management and provide a cooperative environment for all stakeholders of the Han River basin, which serves as the major water source of the South-to-North Water Diversion [42].

Satellite Radar Altimetry
The Ankang satellite altimetry-derived water surface elevations are retrieved from Sentinel-3A. Sentinel-3A is one of the constellations of Sentinel-3 mission with a 27-day cycle. It has been operational since 2016 and thus, nearly five years of data are used in this study. Level 2 data are retrieved using the ESA GPOD (Grid Processing On Demand) platform (https://gpod.eo.esa.int/, last accessed on 8 February 2022). The level 2 product contains 20Hz measurements including the range, satellite altitude, geoid height, waveforms and corrections. The range is retracked with SAMOSA+ retracker [43].
The water surface elevations WSE, are constructed via the following equations: where H alt is the satellite altitude, N geo is the geoid height (i.e., the EGM2008 geoid model), R range is the distance from satellite to water surface, which can be expressed as: where R trac is the onboard tracker range that is 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 ionospheric correction and the wet and dry tropospheric corrections. R geo is the geophysical corrections including pole tide, solid earth tide, and ocean loading tide. Altimetry measurements over reservoir are firstly extracted based on the global surface water occurrence data set (https://global-surface-water.appspot.com/, last accessed on 8 February 2022). To ensure that a higher number of potential measurements are preselected, a low threshold of 10% occurrence is used. Secondly, a filtering process is performed for each pass to remove outliers using the median of absolute deviation [44]. Thirdly, the remaining measurements are used to construct WSE time series using "tyHydro" [45].

Synthetic Aperture Radar Imagery
The Ankang SAR images from Sentinel-1 mission are used to map SWE series over the case study considering two aspects. First, global products (e.g., Global Surface Water Dynamics and Global Surface Water Explorer) from optical satellite observations will tend to have a bias when used at the individual reservoirs. Monitoring SWE using optical data will underestimate water presence in cloudy regions, and is probably constrained by light availability in the darker months at higher latitudes. Second, the major advantage of SAR compared to optical and infrared imagery lies in its day/night observation capacity and its ability to penetrate clouds, hence providing a source for consistent monitoring irrespective of weather. Therefore, we applied a 4-year time series of Sentinel-1 images (early 2016 to 2020 May) to map the SWE dynamics over the reservoir.
Sentinel-1 is composed of a constellation of two near-polar satellites. Both satellites carry a C-band SAR instrument which collects data in all-weather conditions enabling continuous 12 day repeat monitoring. The workflow was adapted from Druce et al. [46] and consists of pre-processing, classification, and post-processing steps. In the pre-processing step, precise orbit vectors and range-Doppler terrain correction are applied to obtain a georeferenced SAR image. To exploit the difference in backscatter from water and non-water surfaces, a trained logistic regression model, trained using data from Global Surface Water Explorer, was used to convert backscatter values into water probability scores. Pixels with a probability score greater than 50% were subsequently classified as water and followed by a simple sieve filtering technique to reduce some inherent speckle noise. Finally, to mask out topography-induced radar shadows falsely classified as water, SRTM DEM was used to prepare a topographic mask that constrained the water detection to be within the high contour line (+10 m) of maximum reservoir extent.

Model Input Forcings and Other Validated Datasets
To facilitate modeling simulations and minimize uncertainties of input forcings, highquality satellite datasets are used to serve the meteorological inputs of hydrological models, i.e., precipitation and temperature. The GPM-Final precipitation products are used for their high quality (https://pmm.nasa.gov/GPM, last accessed on 8 February 2022). ERA5-Land temperature data at hourly and 0.1 • × 0.1 • resolution are acquired from the Copernicus Climate Data Store (https://cds.climate.copernicus.eu/, last accessed on 8 February 2022). Daily water level, storage change, inflow, and release are from the local watershed agency and National Hydrological Information Centre (http://xxfb.mwr.cn/index.html, last accessed 8 February 2022) ( Figure 2). Benchmark curve of elevation-storage is established using in situ measurements. Sentinel-1 is composed of a constellation of two near-polar satellites. Both satellites carry a C-band SAR instrument which collects data in all-weather conditions enabling continuous 12 day repeat monitoring. The workflow was adapted from Druce et al. [46] and consists of pre-processing, classification, and post-processing steps. In the pre-processing step, precise orbit vectors and range-Doppler terrain correction are applied to obtain a georeferenced SAR image. To exploit the difference in backscatter from water and non-water surfaces, a trained logistic regression model, trained using data from Global Surface Water Explorer, was used to convert backscatter values into water probability scores. Pixels with a probability score greater than 50% were subsequently classified as water and followed by a simple sieve filtering technique to reduce some inherent speckle noise. Finally, to mask out topography-induced radar shadows falsely classified as water, SRTM DEM was used to prepare a topographic mask that constrained the water detection to be within the high contour line (+10 m) of maximum reservoir extent.

Model Input Forcings and Other Validated Datasets
To facilitate modeling simulations and minimize uncertainties of input forcings, high-quality satellite datasets are used to serve the meteorological inputs of hydrological models, i.e., precipitation and temperature. The GPM-Final precipitation products are used for their high quality (https://pmm.nasa.gov/GPM, last accessed on 8 February 2022). ERA5-Land temperature data at hourly and 0.1° × 0.1° resolution are acquired from the Copernicus Climate Data Store (https://cds.climate.copernicus.eu/, last accessed on 8 February 2022). Daily water level, storage change, inflow, and release are from the local watershed agency and National Hydrological Information Centre (http://xxfb.mwr.cn/index.html, last accessed 8 February 2022) ( Figure 2). Benchmark curve of elevation-storage is established using in situ measurements.

Methodology
The flowchart of the entire workflow consists of three major steps and is shown in Figure 3. First, reservoir inflow forecasts using four hydrological models and six blending schemes. Second, RWSC estimates by two satellite-based approaches with processed satellite datasets. Third, calculation of release from a reservoir model and evaluation of the contributions and uncertainties of our satellite-based method in estimating reservoir release.

Methodology
The flowchart of the entire workflow consists of three major steps and is shown in Figure 3. First, reservoir inflow forecasts using four hydrological models and six blending schemes. Second, RWSC estimates by two satellite-based approaches with processed satellite datasets. Third, calculation of release from a reservoir model and evaluation of the contributions and uncertainties of our satellite-based method in estimating reservoir release.

Hydrological Modeling Techniques
The four hydrological models (i.e., GR4J, IHACRES, SIMHYD and XAJ) are used to represent the uncertainty of inflow simulations. They are chosen considering three aspects. First, all of these models have been successfully applied in a plethora of catchments for a broad set of use cases such as streamflow simulation and parameter regionalization [47,48]. Second, considering complexities of structures and parameterizations (e.g., 4-15 parameters), Remote Sens. 2022, 14, 815 6 of 18 these models are suitable for blending schemes to facilitate hydrological simulations. Third, as hydrological models are widely used, it is crucial to compare the contributions of different hydrological models to reservoir release estimates using our satellite-based approach.
Remote Sens. 2021, 13, x 6 of 19 Figure 3. Flowchart of the methodology used to assess the contributions of multi-source satellite datasets and hydrological modeling techniques in estimating reservoir release.

Hydrological Modeling Techniques
The four hydrological models (i.e., GR4J, IHACRES, SIMHYD and XAJ) are used to represent the uncertainty of inflow simulations. They are chosen considering three aspects. First, all of these models have been successfully applied in a plethora of catchments for a broad set of use cases such as streamflow simulation and parameter regionalization [47,48]. Second, considering complexities of structures and parameterizations (e.g., 4-15 parameters), these models are suitable for blending schemes to facilitate hydrological simulations. Third, as hydrological models are widely used, it is crucial to compare the contributions of different hydrological models to reservoir release estimates using our satellite-based approach.
The GR4J (Génie Rural à 4 paramètres Journalier) model was proposed by Perrin et al. [49]. It is one of the most widely-used hydrological models for its simplicity and good performance. Based on the unit hydrograph, two runoff components are routed to the outlet, contributing to the hydrograph. The IHACRES model has eight parameters and includes two major modules [50]. Effective rainfall is firstly calculated by a non-linear function. Then, the total runoff is partitioned into fast and slow flows. The SIMHYD (Simple Lumped Conceptual Daily Rainfall-Runoff) model is a simple model of the HYDROL-OGY model [51] with nine parameters. Three storages (interception loss, soil moisture, and groundwater and routing process) are included and different runoff generation mechanisms are incorporated in this model. The XAJ (Xinanjiang) model was developed by Zhao [52], containing four major modules (evapotranspiration, runoff generation, separation, and routing) with 15 parameters.
The Bayesian Model Averaging (BMA) is employed to blend an ensemble of streamflow forecasts by different hydrological models, producing a new forecast posterior The GR4J (Génie Rural à 4 paramètres Journalier) model was proposed by Perrin et al. [49]. It is one of the most widely-used hydrological models for its simplicity and good performance. Based on the unit hydrograph, two runoff components are routed to the outlet, contributing to the hydrograph. The IHACRES model has eight parameters and includes two major modules [50]. Effective rainfall is firstly calculated by a non-linear function. Then, the total runoff is partitioned into fast and slow flows. The SIMHYD (Simple Lumped Conceptual Daily Rainfall-Runoff) model is a simple model of the HYDROL-OGY model [51] with nine parameters. Three storages (interception loss, soil moisture, and groundwater and routing process) are included and different runoff generation mechanisms are incorporated in this model. The XAJ (Xinanjiang) model was developed by Zhao [52], containing four major modules (evapotranspiration, runoff generation, separation, and routing) with 15 parameters.
The Bayesian Model Averaging (BMA) is employed to blend an ensemble of streamflow forecasts by different hydrological models, producing a new forecast posterior probability density function (PDF). Based on the law of total probability, given the observations Y and the streamflow predictions of S models [f 1 , f 2 , . . . , f s ], the posterior PDF of a forecast in BMA can be expressed as: where p(y| f s , Y) is the posterior distribution of y given the model prediction f s and observations Y. p( f s |Y) represents the likelihood of ensemble predictions being correct, which reflects the performance of model s in capturing observations Y. More simply, p( f s |Y) can be rewritten as w s , the weight of each model, therefore, the general BMA approach in Equation (3) returns the weight averaged predicts of each model and can be rewritten as Therefore, the posterior PDF is the weighted average of model predictions using a static value of w for each model. However, the general BMA may not always improve the accuracy of posterior PDF because ensemble predictions are time variant. In this work, time-varying (dynamic) w values are used to facilitate streamflow simulations considering the seasonality or annual effects (Figure 3). To make it clearer, taking the monthly scales (12 w values for each model) as example, the general BMA is repeated for each month, taking the ensemble simulations of each month as inputs, and then obtaining the w of each model in the corresponding month. Finally, the entire series of posterior PDF are constructed with the dynamic w values and ensemble simulations. The parameters w are often iteratively optimized using the Expectation-Maximization (EM) algorithm by maximizing the log likelihood function [53]. EM is an efficient algorithm to iteratively optimize BMA parameters and contains three major modules. Firstly, the Box-Cox transformation is used to follow the Gaussian assumption of the conditional probability distribution p(y| f s , Y) . Secondly, the objective function is set up, i.e., the log likelihood function I(θ), which is preferred and expressed as follows: where g(·) denotes Gaussian distribution, θ is parameters (i.e., weights w s and variance σ 2 s of each model) to be optimized. Thirdly, the EM algorithms were run for optimizing parameters.
The basin-averaged precipitation and potential evaporation are required as inputs to drive models. Potential evaporation is calculated based on the Oudin method, requiring temperature as inputs [54]. We split the meteorological forcings into three periods: 2001-2008 for calibration, 2009-2015 for validation, 1 January 2016 to 31 May 2020 for inflow forecast. Hydrological models are calibrated using the Shuffled Complex Evolution (SCE-UA) with the objective function of KGE [55]. Dynamic BMA at four varying-time scales (i.e., daily, ten-day, monthly, and seasonal) are developed and used to blend ensemble simulations, and two other blending approaches (static BMA, mean ensemble simulations) for comparison ( Figure 3). The weights of BMA are fitted with the ensemble simulations of 2001-2015. To validate the performance of blending schemes, the leave-one-out approach is used [56], fitting BMA weights with the ensemble simulations for all years except 1 year for validation (15 combinations for calibration/validation). Then, we simulate the reservoir inflow from 1 January 2016 to 31 May 2020 for further release estimates.

Satellite-Based RWSC Estimates
The core of satellite-based estimation of RWSC is to construct hypsometry (i.e., elevation-area-storage curve). Based on the multi-source satellite datasets, we can fit the H-A curve with power or polynomial functions [57]. A power function with clear physical meanings is used to fit H-A curve in this study. To minimize the uncertainty, the data (WSE and SWE) acquired within 5 days interval are used to fit the curve by Least Squares Method. When H-A curve is determined, H-V or A-V curve is subsequently derived following the mathematic method.
Remote Sens. 2022, 14, 815 where Pre-defined parameter ranges are the same with Zhong et al. [2], which suggests that a narrower range is enough to cover a wide range of elevation-area curves. RWSC [L 3 ] can be determined by H-V curve and WSE or A-V curve and SWE data. Here, we refer to as WSE-derived method and SWE-derived method. In cases where no satellite imagery available, we extend the WSE-derived method by assuming constant and stageindependent reservoir surface area. Therefore, RWSC can be calculated by multiplying the difference of satellite altimetry observations with a mean area [58]. The three RWSC (V(h), V(A), and constant A) estimates are then used as inputs of reservoir model for release estimates and further validation against in situ records (Figure 3).

A Simple Reservoir Model
Reservoir release can be obtained by using a reservoir model, Equation (8) [29,57].
where release O [L 3 ] is the difference between reservoir inflow I [L 3 ] and RWSC [L 3 ] at a certain time step (e.g., two adjacent satellite altimetry observations with a 27-day cycle). The seepage is assumed negligible. Although a small fraction of evaporation is lost, it can be compensated by precipitation over the reservoir (i.e., evaporation and precipitation account for 0.5% and 0.4% of mean monthly inflow, respectively). Particularly in the case study with subtropical monsoon climate, evaporation rate is small. Thus, it is reasonable to neglect these water fluxes for their little impacts on reservoir release [33,58].

Inflow Simulations
The performances of four hydrological models are assessed by KGE, PBIAS, NSE, logNSE, CC, BR, and RV ( Figure 4 and Table 2). The values of KGE and PBIAS reveal that all hydrological models perform well with KGE > 0.6 and absolute PBIAS < 0.15. Nonetheless, there are some differences. For example, the GR4J exhibits the best simulation skills in terms of highest KGE, NSE, and CC values during calibration period, while the XAJ achieves the best scores during validation period. The XAJ exhibits best skills in reproducing low flows according to its highest logNSE values over 0.65, while the IHACRES is optimal in RV with nearly one. None of these models can fully replicate the variations of streamflow. High flows are often underestimated while low flows tend to be overestimated. Large discrepancies can be found in the evaluation metrics and streamflow hydrographs, indicating that it is necessary to generate a more accurate forecast from the ensemble simulations considering their uncertainties. Figure 5 shows the performance of blending approaches as well as the individuals (i.e., simulations of each hydrological model). The blending approaches perform better than the individuals in terms of evaluation metrics. For example, the KGE values of the blending approaches range from 0.77-0.85, increased by 6.9-18.1% compared with those of the individuals. The uncertainties are greatly reduced by using the blending approaches according to the evaluation metrics values during validation period. It is informative to notice that dynamic BMA schemes achieve better scores than the other two blending approaches, except for dynamic BMA at seasonal scale. Dynamic BMA schemes outperform the other two approaches in terms of NSE, logNSE, CC, PBIAS, and RV metrics. However, static BMA is the best in terms of BR metrics, and also outperforms dynamic BMA at longer time scales such as monthly or seasonal scales, with KGE value improved by 0.2. Figure 5, the best-performing model is dynamic BMA at daily scale. To consider the impacts of different times in dynamic BMA, the trends of three evaluation metrics (KGE, NSE, and logNSE) are analyzed ( Figure 6). The results indicate that the performance of dynamic BMA improves as times increase. The optimal weights of dynamic BMA are shown in Figure 7, supporting the dynamics of probability of the ensembles being correct at different times. Overall, dynamic BMA can facilitate hydrological simulations through considering model uncertainties, and daily time scale is recommended when using dynamic BMA.

Metric Abbreviation Expression Perfect Score
Kling-Gupta Efficiency Root-mean-square error RMSE proved by 0.2. As shown in Figure 5, the best-performing model is dynamic BMA at daily scale. To consider the impacts of different times in dynamic BMA, the trends of three evaluation metrics (KGE, NSE, and logNSE) are analyzed ( Figure 6). The results indicate that the performance of dynamic BMA improves as times increase. The optimal weights of dynamic BMA are shown in Figure 7, supporting the dynamics of probability of the ensembles being correct at different times. Overall, dynamic BMA can facilitate hydrological simulations through considering model uncertainties, and daily time scale is recommended when using dynamic BMA.   uation metrics (KGE, NSE, and logNSE) are analyzed ( Figure 6). The results indicate that the performance of dynamic BMA improves as times increase. The optimal weights of dynamic BMA are shown in Figure 7, supporting the dynamics of probability of the ensembles being correct at different times. Overall, dynamic BMA can facilitate hydrological simulations through considering model uncertainties, and daily time scale is recommended when using dynamic BMA.

RWSC Estimates
The key factors concerning the performance of RWSC estimates are the accuracy of satellite WSE and SWE datasets, which are firstly evaluated against in situ records. The performance of altimetry measurements is evaluated by three statistical metrics (RMSE, BIAS and CC) and considered moderate based on statistical assessment, visual inspection and previous publications [59]. The values of RMSE (bias-removed), BIAS, and CC are 0.36 m, −1.32 m, and nearly one, respectively (Table 3). It should be noted that the bias is likely not a measurement bias but a geoid issue. Datum conversion is not publicly available; thus, we compared the WSE anomalies by removing their mean values. Seasonal deviations are also examined and it was found that high water levels are more accurate in terms of a lower RMSE value of 0.16 m, while low water levels are less accurate with a

RWSC Estimates
The key factors concerning the performance of RWSC estimates are the accuracy of satellite WSE and SWE datasets, which are firstly evaluated against in situ records. The performance of altimetry measurements is evaluated by three statistical metrics (RMSE, BIAS and CC) and considered moderate based on statistical assessment, visual inspection and previous publications [59]. The values of RMSE (bias-removed), BIAS, and CC are 0.36 m, −1.32 m, and nearly one, respectively (Table 3). It should be noted that the bias is likely not a measurement bias but a geoid issue. Datum conversion is not publicly available; thus, we compared the WSE anomalies by removing their mean values. Seasonal deviations are also examined and it was found that high water levels are more accurate in terms of a lower RMSE value of 0.16 m, while low water levels are less accurate with a higher RMSE value of 0.34 m. To validate the performance of SWE datasets, we plotted the time series of in situ water level records and satellite SWE at the overlapping dates. The time series show good agreements, with CC exceeding 0.97. Overall, clear patterns can be found in WSE and SWE, illustrating reservoir behaviors of seasonal filling and emptying (Figure 8).  Figure 3).

RWSC Estimates
The key factors concerning the performance of RWSC estimates are the accuracy of satellite WSE and SWE datasets, which are firstly evaluated against in situ records. The performance of altimetry measurements is evaluated by three statistical metrics (RMSE, BIAS and CC) and considered moderate based on statistical assessment, visual inspection and previous publications [59]. The values of RMSE (bias-removed), BIAS, and CC are 0.36 m, −1.32 m, and nearly one, respectively (Table 3). It should be noted that the bias is likely not a measurement bias but a geoid issue. Datum conversion is not publicly available; thus, we compared the WSE anomalies by removing their mean values. Seasonal deviations are also examined and it was found that high water levels are more accurate in terms of a lower RMSE value of 0.16 m, while low water levels are less accurate with a higher RMSE value of 0.34 m. To validate the performance of SWE datasets, we plotted the time series of in situ water level records and satellite SWE at the overlapping dates. The time series show good agreements, with CC exceeding 0.97. Overall, clear patterns can be found in WSE and SWE, illustrating reservoir behaviors of seasonal filling and emptying ( Figure 8).   Figure 9 shows the error series of RWSC. The performance of RWSC is assessed by four metrics (RMSE, NRMSE, BIAS, and CC) and considered good based on visual inspection and statistical assessment (Table 3). WSE-derived method performs better than SWEderived method in RWSC. For example, the values of RMSE, NRMSE, BIAS, and CC for WSE-derived and SWE-derived RWSC are 0.026/0.099 km 3 , 0.019/0.064, 0.004/−0.003 km 3 , and 0.996/0.914, respectively. Two types of WSE-derived methods perform very similarly, indicating that main source of RWSC error comes from SWE and a combination of SWE and WSE. To figure out the impacts of the propagation of associated satellite datasets uncertainties on RWSC, the error series of WSE, water level change, RWSC estimates as well as elevation-storage curves are depicted in Figure 9. There are three important error sources in the simulated RWSC: elevation error, surface extent error, and error from their combinations. The combinations represent the propagation processes from satellite datasets uncertainties to RWSC (satellite-fitted elevation-storage curve). We established the area-elevation relationship and then derived elevation-storage curve using satellite SWE and WSE (Figure 9d,e). Apparent shifts are found in the satellite-fitted elevation-storage curve (Figure 9d). To eliminate systematic bias, both in situ and satellite-fitted curves are transformed to storage anomaly. We found that the slope of satellite-fitted curve is smaller than that of in situ curve at high water levels, and larger at low water levels (Figure 9e). There is room for improvement and increasing the accuracy of the elevation-storage curve using satellite datasets. Nevertheless, the close agreement between satellite-fitted and in situ curves established a good level of trustworthiness in our satellite-based approaches of calculating RWSC. Based on such satellite-fitted curves, we could obtain RWSC estimates using either SWE or WSE time series. Satellite SWE errors are the main error source of RWSC, supporting the work of Bonnema and Hossain [58]. It is worth noting that there is small difference between two types of SWE-derived RWSC estimates, thereby demonstrating the potential of estimating RWSC with satellite WSE datasets and constant area. The advantage of current satellite sensors cannot be overlooked. Although satellite altimetry measurements are carefully processed, large errors occur sometimes, e.g., at 43rd observations, contributing to the large error of RWSC (Figure 9b). When assuming the reservoir extent is constant, RWSC is calculated by using elevation change, thus, large error of RWSC occurs when two adjacent altimetry measurements are separately underestimated and overestimated. Taking the 15th observations as example, we examined the seasonal performance of satellite altimetry measurements and found there is no obvious pattern in elevation error during the transiting periods. As a result, it is better to reduce the satellite datasets errors before calculating RWSC, which also suggests that WSE-derived method is more suitable to estimate reservoir storage dynamics.

Reservoir Release Estimates
The hydrographs of reservoir release with different input combinations are depicted in Figure 10, and the results of evaluation metrics are shown in Figure 11. Almost all estimated releases show similar patterns to the observations, with KGE and NSE values exceeding 0.65 and NRMSE values below 0.12 km 3 , favoring the success of satellite-based approach. Through improving the accuracy of individual estimates of inflow and RWSC, more accurate release estimates are obtained ( Figure 11). Moreover, impacts of reservoir on natural streamflow are assessed in Figure 10, highlighting the value of satellite-based approach.

Reservoir Release Estimates
The hydrographs of reservoir release with different input combinations are depicted in Figure 10, and the results of evaluation metrics are shown in Figure 11. Almost all estimated releases show similar patterns to the observations, with KGE and NSE values exceeding 0. 65 and NRMSE values below 0.12 km 3 , favoring the success of satellite-based approach. Through improving the accuracy of individual estimates of inflow and RWSC, more accurate release estimates are obtained ( Figure 11). Moreover, impacts of reservoir on natural streamflow are assessed in Figure 10, highlighting the value of satellite-based approach.

Reservoir Release Estimates
The hydrographs of reservoir release with different input combinations are depicted in Figure 10, and the results of evaluation metrics are shown in Figure 11. Almost all estimated releases show similar patterns to the observations, with KGE and NSE values exceeding 0.65 and NRMSE values below 0.12 km 3 , favoring the success of satellite-based approach. Through improving the accuracy of individual estimates of inflow and RWSC, more accurate release estimates are obtained ( Figure 11). Moreover, impacts of reservoir on natural streamflow are assessed in Figure 10, highlighting the value of satellite-based approach.

Discussion
The satellite datasets can characterize the dynamics of reservoirs in otherwise ungauged basins. We select the most recent products and explore how new satellites may deliver useful measurements over reservoirs. To our knowledge, this study firstly exploits the contributions of satellite datasets from Sentinel missions on RWSC and release estimates. SAR expands the applicability of satellite missions in monitoring any inland water body. This is the case for the Ankang, i.e., a reservoir with long-narrow shape and medium size in area and storage. Sentinel missions as Sentinel-1/-3 can deliver useful data and produce valuable information such as RWSC and release. Our SAR-based classification for mapping surface water areas performs considerably better than Landsat-based maps [46], and this will make the estimation of reservoir area more robust and hence improve the accuracy of RWSC and release estimates. The findings indicate that the accuracy of RWSC estimates is more sensitive to satellite SWE while less sensitive to WSE. These discrepancies can be attributed to the reservoir characteristic. To incorporate this property, we plan to evaluate the SWE maps as a next step and may lead to further improvements in the results. To improve satellite-based estimations of RWSC and release, some aspects can be considered. The first is to examine the variations of reservoirs in area [28]. For reservoirs with small variation in area, a simple method (i.e., V (constant A)) would be efficient; otherwise, a use of multi-source satellite datasets is recommended. The second aligns with advances in constructing reservoir geometry using digital elevation model. Liu et al. [9] modeled the reservoir bathymetry with SRTM DEM and satellite imagery. The principle is to utilize geomorphological similarity and topographical continuity of reservoirs to model bathymetry. Last but not least, by ensuring that satellite datasets are accurate with filtering algorithms, performances of RWSC and release estimates can be potentially greatly increased.

Discussion
The satellite datasets can characterize the dynamics of reservoirs in otherwise ungauged basins. We select the most recent products and explore how new satellites may deliver useful measurements over reservoirs. To our knowledge, this study firstly exploits the contributions of satellite datasets from Sentinel missions on RWSC and release estimates. SAR expands the applicability of satellite missions in monitoring any inland water body. This is the case for the Ankang, i.e., a reservoir with long-narrow shape and medium size in area and storage. Sentinel missions as Sentinel-1/-3 can deliver useful data and produce valuable information such as RWSC and release. Our SAR-based classification for mapping surface water areas performs considerably better than Landsat-based maps [46], and this will make the estimation of reservoir area more robust and hence improve the accuracy of RWSC and release estimates. The findings indicate that the accuracy of RWSC estimates is more sensitive to satellite SWE while less sensitive to WSE. These discrepancies can be attributed to the reservoir characteristic. To incorporate this property, we plan to evaluate the SWE maps as a next step and may lead to further improvements in the results. To improve satellite-based estimations of RWSC and release, some aspects can be considered. The first is to examine the variations of reservoirs in area [28]. For reservoirs with small variation in area, a simple method (i.e., V (constant A)) would be efficient; otherwise, a use of multi-source satellite datasets is recommended. The second aligns with advances in constructing reservoir geometry using digital elevation model. Liu et al. [9] modeled the reservoir bathymetry with SRTM DEM and satellite imagery. The principle is to utilize geomorphological similarity and topographical continuity of reservoirs to model bathymetry. Last but not least, by ensuring that satellite datasets are accurate with filtering algorithms, performances of RWSC and release estimates can be potentially greatly increased.
The inflow is a key factor influencing the performance of release estimates, which is reflected by different performances of release using various hydrological modeling techniques (Figures 10 and 11). Besides, the streamflow errors (e.g., timing errors) are likely eliminated when aggregating daily inflow at monthly scale, thus reducing the contributions of varying-time BMA to release estimates. It should be noted that release estimates at shorter scales are also dependent on the satellite-based estimation of RWSC, where we also contributed by deriving bi-weekly RWSC estimation using continuous surface extents. In the case study, reservoir experiences seasonal filling and emptying, which are fully captured by satellite datasets (Figure 8). This is illustrated by Figure A1, which shows monthly average reservoir surface elevation, demonstrating the variations of reservoir operating patterns as well as the potential of estimating release at shorter time scales based on this altimetry-based rule. Although such a satellite-based rule is less applicable, it could be a reliable reference of daily release estimates. For example, Han et al. [33] simulated daily release based on a reservoir rule, which is represented by a storage curve interpolated by monthly storage. The performance of RWSC and release estimates is also expected to be improved as the temporal resolution of satellite datasets is increasing (upcoming SWOT and multiple altimetry missions), potentially producing accurate release simulations at shorter time scales [60].
The results can be valid for the Han River basin and regions with similar climate conditions and water reservoirs in size. In spite of good performance of our approach, further studies should be done for a bigger picture of the performance of reservoir release estimates after repeating our works in other regions. First, the evaporation is negligible compared to inflow in the case study with humid climate conditions, while it could be a potential key variable in semiarid regions [61,62]. Second, the performance of satellite datasets varies with reservoirs of varying sizes and functions. Third, bias is associated with local climates and geographies, bringing a challenge for a wide generalization of the findings.
While the modeling framework adopted in-situ streamflow to calibrate hydrological models and the loss of accuracy is expected when this data is not available, the findings still foster further research initiatives in understanding impacts of reservoirs on river systems in ungauged basins. Progress in remote sensing technologies creates more opportunities in ungauged basins. Several studies have informed hydrological models (i.e., model parameter calibration) with multi-source satellite datasets such as river water level [63], soil moisture [64], evapotranspiration [65], and total water storage [66]. It is expected that more accurate inflow simulations can be still produced by our varying-time BMA schemes, bringing benefits for more accurate estimation of release. The uncertainties of satellite datasets should be paid more attention, and evaluation of the SWE maps accuracy would be desirable and may lead to further improvements in the results. Without additional observations, satellite-based estimations of reservoir records can be the best reference in ungauged basins.

Conclusions
This study evaluates the ability of reconstructing reservoir release for poorly gauged dammed basins using multi-source satellite datasets and hydrological modeling techniques. The experiment is performed at the Ankang reservoir of the Han River basin in China, using ten inflow simulations based on hydrological models and blending schemes and three RWSC estimates based on two satellite-based approaches. The main findings are as follows.
The satellite-based approach achieves satisfactory results in estimating reservoir states such as RWSC and release, thereby supporting further research on dammed basins. The performance of release estimates is variable and greatly affected by the accuracies of inflow and RWSC estimates. Through dynamically correcting the posterior probabilities at temporal scales, more accurate inflow simulations can be produced from our varying-time BMA schemes, bringing benefits for more accurate estimation of release. We contribute to the current understanding of satellite-based estimation of RWSC, which are calculated based on most recent SAR multi-source satellite WSE and SWE datasets. With multi-source approaches, both temporal scale and accuracy of reservoir estimates are increased (i.e., monthly and bi-weekly), acting as a key supplement to in situ records. The results indicate that the accuracy of RWSC estimates is more sensitive to satellite SWE while less sensitive to satellite WSE. In cases where no satellite imagery are available, we suggest calculating RWSC by multiplying elevation changes with constant area as an effective alternative source of information.
Accurate estimation of release is of great significance to hydrologists and local managements in the background of booming reservoir impoundments. The satellite-based estimations of reservoir states (e.g., release and RWSC) are expected to foster further research initiatives in understanding impacts of reservoirs on river systems. Given the background of booming satellite missions, this study can also broaden our understanding of the advantages and limitations when integrating satellite datasets into hydrological modeling for poorly gauged dammed basins.