Using the Digital Elevation Model (DEM) to Improve the Spatial Coverage of the MODIS Based Reservoir Monitoring Network in South Asia

: Satellite remote sensing of near real-time reservoir storage variations has important implications for ﬂood monitoring and water resources management. However, satellite altimetry data, which are essential for estimating storage variations, are only available for a limited number of reservoirs. This lack of high-density spatial coverage directly hinders the potential use of remotely sensed reservoir information for improving the skills of hydrological modeling over highly regulated river basins. To solve this problem, a reservoir storage dataset with high-density spatial coverage was developed by combining the water surface area estimated from Moderate Resolution Imaging Spectroradiometer (MODIS) imageries with the Digital Elevation Model (DEM) data collected by the Shuttle Radar Topography Mission (SRTM). By including more reservoirs, this reservoir dataset represents 46.6% of the overall storage capacity in South Asia. The results were validated over ﬁve reservoirs where gauge observations are accessible. The storage estimates agree well with observations, with coe ﬃ cients of determination ranging from 0.47 to 0.91 and normalized root mean square errors (NRMSE) ranging from 15.46% to 37.69%. Given the general availability of MODIS and SRTM data, this algorithm can be potentially applied for monitoring global reservoirs at a high density.


Introduction
Human-made reservoirs, which are managed by storing and releasing water under predetermined operation rules, play an important role in mitigating floods and improving the efficiency of the water supply for municipal, industrial, and agricultural demands [1][2][3][4]. Although most (if not all) human operated reservoirs are monitored in real-time, reservoir storage information is not commonly available to the public. Indeed, this directly limits the effectiveness of reservoir flow regulation with regard to flood control, water supply, and other purposes-especially for those reservoirs located within transboundary river basins. For instance, the lack of reservoir information for the Mekong River delta has created challenges with regard to flood forecasting in this region [5,6]. In addition, when assessing and predicting the impacts of droughts, the lack of reservoir storage information reduces the reliability of drought analysis systems [7,8].
Due to the limited availability of gauge observations-especially with regard to remote locations, restricted locations, and/or observations over large geographical areas-remote sensing technology provides a promising alternative by monitoring reservoirs from space [4,[9][10][11][12]. With remotely sensed water surface area and elevation data, reservoir storage information can be inferred. Reservoir surface area is commonly estimated by classifying optical satellite imageries [13,14] and surface elevation values are typically obtained from satellite radar altimetry [15,16]. The Geoscience Laser Altimeter System (GLAS) onboard the Ice, Cloud, and Land Elevation Satellite (ICESat) and the Advanced Topographic Laser Altimeter System (ATLAS) onboard ICESat-2 were used to measure the elevation values of relatively small lakes and reservoirs [4,[17][18][19][20].
Even though a variety of remote sensing approaches were developed to monitor reservoir storage from space [21][22][23], they are still insufficient in terms of spatial and temporal coverage-which hinders their applications when high-density reservoir network information is required. For radar altimetry, the restrictions are mainly due to the coarse spatial resolution. With about 3-20 km footprints, it is difficult to capture water surface level values using radar altimetry over reservoirs that are either not large enough or do not overlap with the satellite tracks [24]. Even for lakes that are detectable by radar altimeters, the data may not be accurate enough for applications if the surrounding topography is complex. Consequently, as of 2015 less than 200 large lakes and reservoirs have been observed using the past and current set of radar altimeters [24]. Compared with radar altimeters, the ICESat/GLAS instrument has a distinct advantage with its small footprint (70 m)-but this comes at the cost of a very long return period (91 days). By combining ICESat elevation values and Moderate Resolution Imaging Spectroradiometer (MODIS) area estimations, Zhang et al. [25] developed an algorithm which is partially capable of monitoring South Asian reservoirs at 16-day intervals, with 28% of the total capacity of in the region covered. Despite such progress, the reservoir observation network is still too sparse due to the large spaces between satellites tracks. Water surface area from Landsat and the area-elevation relationship provided by the Shuttle Radar Topography Mission (SRTM) were combined to infer the water level and reservoir storage variations [26][27][28]. Landsat can be used to estimate water surface area for smaller reservoirs and lakes due to its high spatial resolution (30 m). However, its repeat period of 16 days limits its ability to monitor reservoir storage at high temporal resolution-especially when cloud coverage is too thick. Therefore, the lack of dense spatial and temporal representation from satellite altimeters remains a major challenge for collecting reservoir storage information on a large scale.
South Asia, which contains one of the largest and densest populations, suffers the most from the dearth of reservoir storage data sharing. The deficient communications with regard to reservoir storage (and management decisions) further exacerbate the casualties and economic losses from flood events. According to past statistical records, South Asia experiences one of the highest fatality rates in the world caused by floods [29]. The available remotely sensed reservoir storage datasets only sparsely cover the region. For instance, radar altimetry data are only available for six reservoirs in this region, which accounts for 10.70% of the total capacity in South Asia (according to Hydrology by altimetry data from Laboratoire d'Etudes en Géophysique et Océanographie Spatiales (LEGOS) [30] and the Global Reservoir Lake Monitor [31]). Although the use of ICESat elevation data improved the coverage to around 28% of South Asian reservoirs [25], it still does not meet the strong societal need. Therefore, acquiring reservoir storage information with large spatial coverage is critical for minimizing the vulnerabilities and maximizing the benefits to communities in this region through good reservoir management practices.
To extend the spatial coverage where remote sensing reservoir storage data are available, a reservoir storage dataset was developed by leveraging the global coverage capability of the Digital Elevation Model (DEM) collected by SRTM. Although DEMs have been most commonly used for generating river routing networks [32,33], they have also been adopted in studies to estimate glacier variations [34,35] and surface water storage change [36]. Due to its high consistency, accuracy, and global coverage [35,37], the SRTM DEM was used to extract the area-elevation (A-H) relationship for calculating reservoir storage in this study.
Our overarching goal was to improve the spatial coverage of the remotely sensed reservoir storage dataset in the South Asia region. To this end, the A-H relationship of a given reservoir was first derived from MODIS water surface area values and SRTM DEM surface heights, and then combined with the area time series to estimate storage variations. The results were validated with gauge observations. The performance of the generated reservoir dataset was compared with the ICESat based algorithm reported by Zhang et al. [25]. In addition to the data analysis and the results validation, storage estimation uncertainties due to reservoir surface area retrieval algorithm parameterization and elevation measurement errors were also quantified.

Remote Sensing Data
In this study, the two primary remote sensing datasets were the STRM DEM and the MODIS imageries. The DEM was used for inferring the A-H relationship. The DEM data were collected by SRTM during an 11-day mission in February 2000, covering a near-global domain from 56 • S to 60 • N [38]. The relative vertical accuracy was~6m, and the absolute accuracy was~16 m [37]. The NASA SRTM V3.0 dataset provides land surface elevation values at a 30-m spatial resolution globally. Here, the global SRTM DEM dataset was obtained from the U.S. Geological Survey's Long Term Archive [39].
For each given reservoir, MODIS imageries were used to derive surface area estimations, which were then applied to the A-H relationship to generate a long-term time series of reservoir storage. The reservoir surface area was calculated from the MODIS/Terra 16-day, 250-m resolution vegetation indices product (MOD13Q1). Specifically, an image classification algorithm (Section 3.2.1) was applied to the Normalized Difference Vegetation Index (NDVI) imageries to extract the reservoir area. From 2000 to 2015, a total of 365 imageries were processed for each reservoir.

Data for Validations
Gauge observations released by the Indian Central Electricity Authority (CEA, [40]) were used to validate the remotely sensed reservoir storage dataset. This gauge data contained daily reservoir water level and storage information for 30 hydropower reservoirs. We downloaded the record from 2008 to 2011 and from 2013 to 2016 in May 2016.
Additionally, the reservoir storage results derived from MODIS and SRTM were compared against the previous results from MODIS and ICESat [25]. Because the Zhang dataset contains results from 21 South Asian reservoirs, this cross-validation helped us to better understand the overall performance of this new dataset on a regional scale.

Reservoir Selection
Two criteria were used to identify the reservoirs included in this study: First, the reservoir maximum area at capacity needed to be larger than 55 km 2 . The threshold of 55 km 2 was based on a comprehensive consideration of both estimation accuracy and spatial coverage. This would guarantee that the surface area could be estimated with high accuracy using medium-resolution MODIS imageries. Reservoirs larger than 55 km 2 account for~46.6% of the total South Asian reservoir capacity. Second, the surface area according to the SRTM DEM for a reservoir of interest should not reach its maximum surface area (estimated from MODIS). Otherwise, the respective ranges of area and elevation detected by SRTM DEM would have been too small to infer the A-H relationship accurately. Following the above criteria, a total of 28 reservoirs were chosen from the Global Reservoir and Dam (GRanD) database [41]. Figure 1 shows the locations of these reservoirs, and compares the reservoirs from this study with those in Zhang et al. [25], with details shown in Table 1. Yellow dots represent reservoirs that can only be monitored by the Moderate Resolution Imaging Spectroradiometer-Ice, Cloud, and Land Elevation Satellite (MODIS-ICESat). Green dots are reservoirs that can only be monitored through the MODIS-Shuttle Radar Topography Mission (STRM). Red points are reservoirs that can be monitored by both approaches. For each reservoir, detailed information is provided in Table 1.

Methodology for Reservoir Storage Estimation
The MODIS-SRTM-based reservoir storage estimation algorithm-referred to as the "MODIS-SRTM algorithm" hereafter-is illustrated using the flowchart in Figure 2. It mainly contains three step: First, the water surface area was estimated from MODIS NDVI imageries via an enhanced classification procedure; second, the A-H relationship was generated from the DEM information by regressing the cumulative area values against their corresponding elevation values (within the delineated reservoir maximum domain; and third, by applying the water surface area estimations to the A-H relationship, the reservoir storage variations were calculated. Further details of these steps are provided as below.

Surface Area Estimation
For each given reservoir, the water surface area was estimated using the enhanced K-means classification approach developed by Zhang et al. [25]. First, a threshold of 0.1 was applied to each 16-day MODIS NDVI image from 2000 to 2015, where pixels with NDVI values less than 0.1 were considered water. Based on these simplified classifications, a mask image was created to represent the water coverage percentile and to delineate the domain of the reservoir. Then, the K-means clustering algorithm [42] was used to identify all water pixels within the masked area of the MODIS NDVI images. Finally, a classification enhancement procedure was applied to finetune the results from the K-means clustering. The main purpose of the enhancement was to use the water occurrence map as a reference to correct misclassified pixels and/or to assign an appropriate class to the unclassified pixels [25].

Area-Elevation (A-H) Relationship Development
The SRTM DEM data were used to develop the A-H relationship for each reservoir. As a valid approximation, the relationships for all reservoirs were assumed to be linear (H = kA + b, where k is the slope of the A-H relationship, and b is the intercept) [43]. To capture the relationship, we first delineated the water surface area from the DEM for each reservoir of interest. For a given reservoir, the water surface area during the SRTM acquisition time was expanded to include its surrounding pixels by gradually increasing the surface elevation threshold, with the water surface elevation corresponding to the DEM area as the initial value. During this process, all pixels that were not directly connected to the increasing water area were discarded as noise. This expansion continued until the new area on this DEM reached the maximum reservoir area estimated from the MODIS images (from 2000 to 2015). This maximum reservoir area was then delineated from the SRTM DEM. A simplified example of a delineated reservoir is shown in Figure 3a. After delineating the maximum coverage of the reservoir from the DEM, the cumulative area (e.g., A 3 ) at any given elevation value (e g., H 3 ) could be estimated by counting the number of pixels with elevations equal to or smaller than that value (i.e., H 3 ). By regressing the cumulative area values against the elevation values, the A-H relationship for the reservoir of interest was established (Figure 3b). A real example of the A-H relationship development for the Pong reservoir is provided in Figure 3c,d.

Storage Estimation
Reservoir storage can be estimated based on the remotely sensed water surface area and elevation using Equation (1): where V C , A C , and H C represent the storage, area, and water elevation at capacity, respectively. V RS , A RS , and H RS are the remotely sensed storage, area, and water height at the monitoring time.
In this MODIS-SRTM algorithm, since H RS can be calculated by applying the A-H relationship to the MODIS area estimation (i.e., A RS ), the reservoir storage is calculated through Equation (2) (which was transformed from Equation (1)).
Using the methods explained in this section, the reservoir storage was calculated for the 28 selected reservoirs in South Asia from 2000 to 2015. Using the Hirakud reservoir as an example, Figure 4b compares the time series of reservoir storage from this MODIS-SRTM algorithm with that from the MODIS-ICESat algorithm by Zhang et al. [25]. Results suggest that these two sets of storage estimations are in good agreement. However, compared with the MODIS-ICESat-based algorithm, the storage values from this study tend to be underestimated (due to the different A-H relationships). To better understand the error statistics of these two approaches, validations using gauge data were conducted and are reported on in Section 4.1.

Results
The MODIS-SRTM-based reservoir storage dataset was examined comprehensively from three perspectives: First, the reliability of the dataset was tested by validating the MODIS-SRTM based reservoir storage results with both in situ gauge data and the MODIS-ICESat based results. Second, the enhanced spatial coverage from this new dataset was compared with the existing reservoir storage dataset in South Asia. Third, the uncertainties associated with the algorithm and dataset were analyzed.

Validation Results
The MODIS-SRTM-based reservoir storage was validated over 11 reservoirs (Table 2) where gauge observation data were available. The performance of the results was evaluated using Equations (3)-(5), which represent three statistical criteria: the coefficient of determination (R 2 ), the relative bias (B), and the normalized root mean square error (NRMSE): where RS represents the remotely sensed results, Obs is the gauge data, i denotes the i th record, n is the total number of data points, and RS and Obs are the average values of the remote sensing results and the gauge data, respectively. As shown in Table 2, most of these results were highly correlated with CEA gauge observations. The R 2 values ranged from 0.47 to 0.91, with a mean of 0.8. The lowest R 2 was found over the Ranjit Sagar reservoir. This reservoir has a relatively small area (56 km 2 at capacity) and is very meandering with a high shoreline to area ratio, complicating the accurate estimation of the surface area from the medium spatial resolution MODIS data [9]. This multicriteria evaluation provided a comprehensive understanding of the results. Using the Srisailam reservoir as an example, its R 2 value was the second highest among all of the validated reservoirs, but its NRMSE was relatively large. Because the slope of the A-H relationship (k, in Equation (2)) is constant, a high R 2 value suggests that the area estimations are accurate. Thus, the large NRMSE was mainly caused by errors associated with the slope of the A-H relationship. Because the area error was proven to be small as indicated by the large R 2 , the SRTM DEM was thus the primary error source for the storage results for this reservoir. Another example is the Ranjit Sagar reservoir. Although it had an extremely low R 2 value due to the large amount of error in the surface area estimations, the storage bias was close to those of the Pong and Rihand reservoirs, which indicates a relatively more accurate A-H relationship over this reservoir.
The performance of this algorithm was also compared with the MODIS-ICESat algorithm by Zhang et al. [25] ( Table 3). The remotely sensed reservoir storage data from these two algorithms were validated over five reservoirs (Hirakud, N. J. Sagar, Pong, Rengali, and R. P. Sagar) where gauge observations and A-H values were available (from both MODIS-ICESat and SRTM). As shown in Figure 5a, both the MODIS-SRTM and MODIS-ICESat-based approaches performed well overall. The time series from these two algorithms closely matched the gauge values for reservoir storage. To highlight the differences between the DEM and ICESat based algorithms, Figure 5b compares the storage errors against the gauge observations from these two datasets. The error statistics are provided in Table 3. Among each of the five reservoirs, the NRMSE of the MODIS-SRTM algorithm ranged from 18.14% to 27.95%, with a mean value of 21.94%. The relative bias values ranged from −11.07% to 19.25%. The NRMSE of the MODIS-ICESat algorithm ranged from 14.58% to 26.50%, with a mean value of 18.83%. The bias values ranged from −8.97% to 4.41%. In terms of accuracy, the two approaches performed relatively similarly, with the MODIS-ICESat algorithm slightly better than the DEM based algorithm. For the N. J. Sagar reservoir, the NRMSE was 27.95% for the MODIS-SRTM algorithm and 26.50% for the ICESat-based algorithm. For this reservoir, the DEM results were more accurate than the ICESat results. The NRMSE was 15.00%, which was 3.18% better than the ICESat based algorithm. For the Hirakud, Pong, and R. P. Sagar reservoirs, the MODIS-ICESat algorithm showed a superior accuracy when validated against the gauge data. The higher accuracy of the MODIS-ICESat algorithm at these three reservoirs may be attributed to the higher vertical accuracy of the ICESat elevation values, and/or the longer observation period of ICESat (than the DEM, which results in a more representative A-H relationship). Because the ICESat and SRTM approaches both use the same MODIS water area values, the larger bias of storage from the SRTM DEM implies that the lower accuracy of SRTM could reduce the quality of the reservoir storage product. As stated by the authors of [44], the components of the SRTM error include baseline roll error, phase error, beam differential errors, and timing and position errors. However, the SRTM DEM errors related to terrain, timing, and position-along with the low vertical resolution (1-m intervals)-still influenced the accuracy of the A-H relationship, which led to a higher bias of the storage calculation. Overall, the MODIS-SRTM algorithm performed reasonably well.

Spatial Coverage of the Reservoir Storage Dataset
With full-coverage of two-dimensional elevation data at a fine spatial resolution (30 m), the MODIS-DEM algorithm generated storage time series for the 28 reservoirs in South Asia from 2000 to 2015 ( Figure 6). These reservoirs had an integrated capacity of 124.17 km 3 (46.6% of the region's total capacity). Compared with the MODIS-ICESat algorithm, the MODIS-SRTM algorithm enabled the monitoring of eight additional reservoirs (Figure 1), which represented a 18.6% increase of the overall storage capacity. Sriram Sagar, which was almost at its maximum level during the SRTM flight time, was the only reservoir for which the A-H relationship could be generated by MODIS-ICESat but not by the DEM.  The new dataset contains the storage variation information over multiple reservoirs at the basin scale, which is essential for regional water management purposes. For instance, with two additional reservoirs included in the dataset, the total storage of the monitored reservoirs in the Krishna river basin (KRB) increased from 33.4% to 67.0% (i.e., from 9.70 to 19.44 km 3 ) of the basin's capacity (29 km 3 ). The Krishna River is the fourth largest river in India, with its basin extending over an area of 259,948 km 2 (about 8% of India). Most of the KRB is relatively flat, with about 76% of the basin covered by agricultural land. Many hydroelectric power stations are distributed along the Krishna River, providing clean energy to a large area of India. Therefore, the improved spatial representativeness of reservoirs in this river basin is essential for hydrologic modeling and water management. The Ukai Dam across the Tapti River was constructed for the purposes of irrigation, hydropower generation, and flood control. The Tapti River basin accounts for nearly two percent of the total area of India. However, before this study, no reservoir in this basin had remotely sensed elevation or storage data from space. In 2000, a severe drought occurred in the Tapti basin, causing drinking water scarcity in some villages [45]. In 2009, many districts in this basin were declared to be under drought conditions due to the deficiency of rainfall from June to September [46]. The low storage values around 2000 and 2009 ( Figure 6) reflect this water scarcity. Figure 6 also shows an increase of maximum storage in the Mangla Reservoir after 2012. This is attributed to the enhanced storage capacity, that was used to increase the reservoir's irrigation capability [47]. Another example is the Yeldari reservoir. According to media reports, two severe drought events occurred in the region in 2004 and from 2012 to 2015-and, in both cases, the Yeldari reservoir almost dried up [48,49].

Uncertainty Analysis
The storage uncertainty associated with the A-H relationship is primarily attributed to the use of partial bathymetry information to represent the A-H relationship for the entire reservoir. Because the DEM dataset only represents the part of the bathymetry that was above the water surface when the SRTM measurements were collected, it assumes that the unmeasured part below the water surface shares the same A-H relationship. To quantify the uncertainty associated with this assumption, we compared the storage estimations from three different scenarios (Figure 7). In each case, a simplified cross-sectional view of the reservoir was used-with the water surface area collected by the SRTM (in 2000) indicated as A 1 , and the area of the reservoir bottom indicated as A 2 . Under all scenarios, the storage volume below the DEM water surface was preserved. The first scenario (Figure 7a) follows the algorithm used in this study, which assumes that the A-H relationship remains the same across the entire profile. The second scenario (Figure 7b) assumes that the area of the reservoir bottom is zero, and thus the A-H relationship of the unknown part below the water surface has the smallest possible slope of k min . The third scenario (Figure 7c) assumes that the minimum area from the MODIS estimations is the area of the reservoir bottom, and thus the A-H relationship of the unknown part below the water surface has the largest possible slope of k max . Using the slope of the upper portion (i.e., k) as estimated from the DEM, the reservoir storage value when the DEM was constructed (i.e., V 2 ) can be calculated using Equation (6): As shown in Figure 7b, the minimum value of A 2 -which is 0-can be used to estimate k min via Equation (7): Similarly, the maximum value of A 2 -which is equal to the minimum water surface area from MODIS during the research period (Figure 7c)-can be used to estimate k max after Equation (8): Thus, for any MODIS remotely sensed area (A RS ) that is less than A 1 , the storage can range between a minimum possible value of V RS min (Equation (9)) and a maximum value of V RS max (Equation (10)): Therefore, the uncertainties associated with the constant slope assumption can be represented by the difference between the two storage estimates described below using Equation (11): The uncertainties associated with this source are illustrated in Figure 8. For all 28 reservoirs, the absolute uncertainty due to the unmeasured A-H relationship ranged from 0 km 3 to 0.54 km 3 , with an average of 0.23 km 3 . Among these reservoirs, the Rihand reservoir had the largest absolute uncertainty (0.54 km 3 ), primarily because this large reservoir was at a relatively high level when the DEM data were collected. The surface area of the Rihand reservoir-as measured by DEM-was 388.96 km 2 , whereas its surface area at full capacity is 485 km 2 . Considering all of the reservoirs, we found a significant increasing trend of the absolute uncertainty as the reservoir capacity increased. For every 1 km 3 increase in reservoir capacity, the uncertainty increased by 0.034 km 3 (p < 0.01). The averaged relative uncertainty caused by the unmeasured A-H relationship was 4.68%. However, we observed no significant relationship between the relative uncertainty and the capacity.
The uncertainties from the area estimation algorithm were quantified thoroughly by Zhang et al. [2014]. From this source, the absolute uncertainties were also found to be highly correlated with the storage at capacity, where the absolute uncertainties had an average value of 3.90%. This is a similar uncertainty range but lightly larger than the unmeasured A-H relationship.
The vertical error of SRTM DEM could be another source of uncertainty. This was not analyzed in this study because the storage calculation (in this study) was based on the slope of the A-H relationship and area estimations, rather than using the absolute elevation value from the SRTM DEM directly. Since the slope of the A-H relationship is determined by the elevation difference of reservoir pixels, the absolute vertical DEM error can be offset during the process, reducing its influence on the storage estimation.

Conclusions
In this study, an algorithm that leverages the SRTM DEM data was developed to improve the spatial coverage of the reservoir monitoring network in South Asia. By combining water surface area from MODIS for reservoir storage estimations, we were able to take the advantage of high temporal resolution of MODIS and large spatial coverage of SRTM. Furthermore, validation results against gauge observations over 11 reservoirs in South Asia suggested that the storage estimations had a good level of accuracy (with R 2 values ranging from 0.47 to 0.91). The integrated storage capacity of these reservoirs was 118.76 km 3 , which represents 46.6% of the overall storage in the region.
This algorithm still has some limitations that need to be noted. First, the accuracy of the proposed algorithm depends on the water level at the time the DEM data were collected. For certain reservoirs that were almost full during the SRTM acquisition time, this approach did not work. Due to the assumption that the A-H relationship derived from the DEM above the water surface represented the full bathymetry, uncertainties in storage estimations were introduced in addition to those from the area retrieval algorithm. Second, the low vertical resolution of SRTM DEM and the errors from different sources may reduce the accuracy of the storage estimation [44]. Therefore, examining the DEM errors with respect to the terrain of the reservoirs could help us to better understand the error characteristics of the storage estimation bias. Third, due to the medium resolution of MODIS, the accuracy of reservoir storage estimation decreased for the reservoir with the smallest surface area (56 km 2 ). Nonetheless, the benefits of the extended number of reservoirs outweigh the constraints.
The algorithm proposed in this study can provide reservoir storage products that support water management on a large scale. For instance, given the long-term availability of high spatial resolution sensors, this approach could be used to monitor much smaller sized reservoirs than possible using existing techniques. This algorithm may also contribute to future satellite missions such as the Surface Water Ocean Topography (SWOT) mission, which will provide a direct water surface measurement for about two-thirds of global lakes and reservoirs, including those with an individual water area > 0.06 km 2 . Funding: This work was supported by the NASA Science of Terra, Aqua, and Suomi NPP (TASNPP) Program (80NSSC18K0939) provided to Texas A&M University. It has benefitted from the computing resources of the Texas A&M Supercomputing Facility (http://sc.tamu.edu). The authors would also like to thank the Central Electricity Authority of India for providing the reservoir gauge observations.

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