Flood Mitigation in the Transboundary Chenab River Basin: A Basin-Wise Approach from Flood Forecasting to Management

Rapid and reliable flood information is crucial for minimizing post-event catastrophes in the complex river basins of the world. The Chenab River basin is one of the complex river basins of the world, facing adverse hydrometeorological conditions with unpredictable hydrologic response. Resultantly, many vicinities along the river undergo destructive inundation, resulting in huge life and economic losses. In this study, Hydrologic Engineering Centre–Hydrologic Modeling System (HEC-HMS) and HEC–River Analysis System (HEC-RAS) models were used for flood forecasting and inundation modeling of the Chenab River basin. The HEC-HMS model was used for peak flow simulation of 2014 flood event using Global Precipitation Mission (GMP) Integrated Multisatellite Retrievals-Final (IMERG-F), Tropical Rainfall Measuring Mission_Real Time (TRMM_3B42RT), and Global Satellite Mapping of Precipitation_Near Real Time (GSMaP_NRT) precipitation products. The calibration and validation of the HEC-RAS model were carried out for flood events of 1992 and 2014, respectively. The comparison of observed and simulated flow at the outlet indicated that IMERG-F has good peak flow simulation results. The simulated inundation extent revealed an overall accuracy of more than 90% when compared with satellite imagery. The HEC-RAS model performed well at Manning’s n of 0.06 for the river and the floodplain. From the results, it can be concluded that remote sensing integrated with HEC-HMS and HEC-RAS models could be one of the workable solutions for flood forecasting, inundation modeling, and early warning. The concept of integrated flood management (IFM) has also been translated into practical implementation for joint Indo-Pak management for flood mitigation in the transboundary Chenab River basin.

The study area is divided in two components: (1) the upper catchment, which includes 20 subbasins contributing the flows to the Marala barrage (outlet), and (2) the lower part, which includes the river floodplain from the Marala barrage to the Chenab-Ravi Rivers' confluence, as show in Figure 1. The geographic extent of the Chenab River catchment lies 74-77.85 • E and 32-34.3 • N, while elevation ranges from 240 to 7085 m. The geographic extent of the floodplain lies 71.6-74.47 • E and 30.4-32.8 • N, while elevation ranges from 131 to 235 m. The transboundary Chenab River basin has a total area of 67,430.34 km 2 , of which 28,480.97 km 2 of catchment area is lying above the Marala barrage. The catchment area above the Marala barrage is divided among Pakistan (402 km 2 ), Jammu and Kashmir (20,139.7 km 2 ), and India (7939.27 km 2 ). It clearly reflects that 27.87%, 70.71%, and 1.42% area of the catchment lies in India, Jammu and Kashmir, and Pakistan, respectively. The transboundary Chenab River in the Indian territory extends to about 230 km, and then the river enters the Jammu and Kashmir region, where it flows across a length of about 276 km. Then it enters the Marala barrage at Sialkot (Pakistan) after covering a total distance of about 506 km. The Chenab River then flows from the Marala barrage in the southwest direction to the Khanki barrage after covering a 59 km distance. The river further flows in the same direction for a downstream distance of 29 km to the Qadirabad barrage. It then proceeds toward the Trimmu barrage located at a distance of 247 km, where the Jhelum River diverts its water into the Chenab River. The last point of the study area is the confluence of the Chenab and Ravi Rivers, which is about 75 and 410 km downstream of the Trimmu and Marala barrages, respectively. The slope of the river in the upper catchment is about 25 m/km, and in the floodplain, it is about 0.4 m/km [101]. The Chenab River catchment is a snow-dominant area, and about 65% of precipitation is received in monsoon and pre-monsoon, while 26% precipitation is received in winter in the form of snow [102]. Therefore, the river receives a high magnitude of flow due to heavy precipitation in monsoon, which is augmented by the overlapped snowmelt due to the summer season. Conclusively, the peak flow from the Chenab River catchment occurs from June to September [99,103]. the HEC-RAS model and satellite imagery. The study also aims to propose a way forward for flood mitigation in the transboundary Chenab River basin through implementation of IFM.

Study Area
The study area is divided in two components: (1) the upper catchment, which includes 20 subbasins contributing the flows to the Marala barrage (outlet), and (2) the lower part, which includes the river floodplain from the Marala barrage to the Chenab-Ravi Rivers' confluence, as show in Figure 1. The geographic extent of the Chenab River catchment lies 74°-77.85° E and 32°-34.3° N, while elevation ranges from 240 to 7085 m. The geographic extent of the floodplain lies 71.6°-74.47° E and 30.4°-32.8° N, while elevation ranges from 131 to 235 m. The transboundary Chenab River basin has a total area of 67,430.34 km 2 , of which 28,480.97 km 2 of catchment area is lying above the Marala barrage. The catchment area above the Marala barrage is divided among Pakistan (402 km 2 ), Jammu and Kashmir (20,139.7 km 2 ), and India (7,939.27 km 2 ). It clearly reflects that 27.87%, 70.71%, and 1.42% area of the catchment lies in India, Jammu and Kashmir, and Pakistan, respectively. The transboundary Chenab River in the Indian territory extends to about 230 km, and then the river enters the Jammu and Kashmir region, where it flows across a length of about 276 km. Then it enters the Marala barrage at Sialkot (Pakistan) after covering a total distance of about 506 km. The Chenab River then flows from the Marala barrage in the southwest direction to the Khanki barrage after covering a 59 km distance. The river further flows in the same direction for a downstream distance of 29 km to the Qadirabad barrage. It then proceeds toward the Trimmu barrage located at a distance of 247 km, where the Jhelum River diverts its water into the Chenab River. The last point of the study area is the confluence of the Chenab and Ravi Rivers, which is about 75 and 410 km downstream of the Trimmu and Marala barrages, respectively. The slope of the river in the upper catchment is about 25 m/km, and in the floodplain, it is about 0.4 m/km [101]. The Chenab River catchment is a snow-dominant area, and about 65% of precipitation is received in monsoon and pre-monsoon, while 26% precipitation is received in winter in the form of snow [102]. Therefore, the river receives a high magnitude of flow due to heavy precipitation in monsoon, which is augmented by the overlapped snowmelt due to the summer season. Conclusively, the peak flow from the Chenab River catchment occurs from June to September [99,103].

Hydrometeorological, Soil, and Landuse Datasets Used in the Study
Flood forecasting and early warning through hydrological and hydrodynamic modeling in HEC-HMS and HEC-RAS models require high spatiotemporal and reliable landuse, soil types, and all the related hydrometeorological information about the river basin [104][105][106][107][108]. Most of the Chenab River catchment is under Indian control; therefore, timely and reliable in situ hydrometeorological information is a serious issue for hydrological studies and flood forecasting. For this purpose, several satellite and global-gridded datasets were used for landuse, soil types, precipitation, temperature, and SWE information. IMERG-F for GPM Level 3 Version 6 (GPM_3IMERGDF) was used in this research, which provides daily multi-satellite precipitation at a spatial resolution of 0.1 • [109]. The daily IMERG-F data depend upon the availability of half-hourly data, and the latency period for the daily data can be from 3 to 4 months. The data are processed using the Goddard Profiling Algorithm (GPROF2017) [110,111] and then gridded and intercalibrated to the GPM Combined Ku Radar-Radiometer Algorithm (CORRA) and merged into half-hourly data [112]. The daily Tropical Rainfall Measuring Mission (3B42 near Real Time (RT)) Version 7 was also used in this study at a spatial resolution of 0.25 • , which is generated from 3-hourly TRMM Multi-Satellite Precipitation Analysis (TMPA) [113,114]. The TRMM_3B42RT is available after 2 months of observation due to data incorporation from multi-satellites, gauge calibration, sensor calibration, and post-processing using algorithms [115]. The GSMaP estimates have various input data sources of high spatiotemporal resolution, comprising extensive satellite information from infrared and passive microwave radiometers [116][117][118][119]. The GSMaP_NRT was also used in this research with a spatial resolution of 0.1 • with a low latency period (4 h) [120,121], which is the core cause for its use in monitoring the hydrological response in the river basins. The Glob-Snow SWE Level 3A product, developed by the Finnish Meteorological Institute (FMI) and Environment Canada, was used in this study at a spatial resolution of 0.25 • (https: //www.globsnow.info/swe/GlobSnow2_SE_SWE_Product_User_Guide_v1_r1.pdf, accessed on 21 August 2021).
Land surface temperature (LST) also plays an important role in predicting the hydrological response of the snow-dominant river catchments in the form of snowmelt [122]. For the 1992 flood event, the daily mean temperature AphroClim_V1808 (over Monsoon Asia APHRO_MA_TAVE_CLM_V1808) data was used at a spatial resolution of 0.05 • (http://aphrodite.st.hirosaki-u.ac.jp/products.html, accessed on 21 August 2021). The said temperature data was developed by the APHRODITE project [123], including the Himalayas, South and Southeast Asia, and mountainous areas in the Middle East [124]. For the 2014 flood event, the global near-real-time (NRT) hourly LST data at a spatial resolution of 0.05 • of geostationary satellites were incorporated in the analysis. The said NRT hourly LST data were produced after less than 3h of observation of top-of-atmosphere brightness temperatures, from the infrared spectral channels of a constellation of geostationary satellites (Meteosat Second Generation, Geostationary Operational Environmental Satellite, and Multifunction Transport Satellite/Himawari) [125].
The landuse change detection of the river catchments is a major challenge [126] and helps to understand hydrological processes and associated systems, which define the water yield from the river catchments [38,[127][128][129]. Due to lack of in situ data, the global land cover map was used for 2014 for the Chenab River catchment, which was produced at a spatial resolution of 300 m by the European Space Agency (ESA) Climate Change Initiative (CCI; https://maps.elie.ucl.ac.be/CCI/viewer/, accessed on 21 August 2021) [130]. The landuse for the Chenab River floodplain was developed using Landsat 5 Thematic Mapper (TM) and Landsat 8 Operational Land Imager (OLI) for the 1992 and 2014 flood events, respectively. The supervised classification was performed using the training data (groundtruth data collected from the farming community) for each earth feature in the Chenab River floodplain for both flood events. The Normalized Difference Vegetation Index (NDVI) [131], as presented in Equation (1), was used for landuse classification at a spatial resolution of Remote Sens. 2021, 13, 3916 6 of 30 30 m for the Chenab River floodplain, using green and NIR bands of Landsat 5 TM and Landsat 8 OLI, as presented in Table 1.
The main purpose of the detailed landuse classification was to explore the earth features, especially the vegetation types, in the Chenab River floodplain. The Manning's roughness coefficient (n) for each earth feature was selected from the literature [132,133] and incorporated in the numerical simulation of the HEC-RAS model.
The soil types and texture are also important, along with landuse, to understand the hydrological response of a complex river catchment [134,135]. Harmonized World Soil Database (HWSD) version 1.2 was used, which is produced at a spatial resolution of 1 km by a collaboration between the UNESCO and the Food and Agriculture Organization (FAO) with the International Institute for Applied Systems Analysis, the International Soil Reference and Information Centre (ISRIC)-World Soil Information, the Institute of Soil Science of the Chinese Academy of Sciences, and the Joint Research Centre of the European Commission (http://www.fao.org/soils-portal/data-hub/soil-maps-anddatabases/harmonized-world-soil-database-v12/en/, accessed on 21 August 2021) [136].

Peak Flow Simulation Using the HEC-HMS Model
The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Digital Elevation Model (DEM) was used at a spatial resolution of 30 m for delineation of watershed, river, and river catchment using the HEC-GeoHMS model. The resultant delineated catchment and river network were imported into the HEC-HMS model for hydrological modeling. The HEC-HMS model is designed to simulate the hydrological processes of complex and dendritic river basins [137]. The HEC-HMS model can be applied to analyze the flood magnitude, flood frequency, flood warning, and spillway and reservoir planning, and can also be used for event-based and continuous hydrological modeling [46,138,139]. The soil moisture accounting (SMA) method was used for loss estimation for calibration, validation, and application periods of the HEC-HMS model, and it has been applied for event and continuous modeling worldwide [140][141][142]. The SMA method accounts for retention by canopy cover and simulates the water through soil layers to the groundwater [143]. The SMA requires the initial soil moisture condition for each subbasin at the beginning of the simulation. As a physical survey is not possible due to transboundary conflicts, the daily TMI/TRMM surface soil moisture information was retrieved from Land Parameter Retrieval Model (LPRM) Level 3 at a spatial resolution of 0.25 • [144].
The snowmelt model is important for quantification of the streamflow from snowdominant watersheds. The temperature index (TI) method was used for snowmelt runoff modeling in the HEC-HMS model due to its worldwide applications [145][146][147][148]. The TI method represents cold energy stored in the snowpack (SWE) with past conditions and all the other factors for estimation of melting for each degree above freezing temperature. The melt rate of the model is estimated by considering the present and past situations of the snowpack in the river catchment. The PX temperature (threshold temperature) differentiates snow from liquid rain, and the base temperature separates melt and nonmelt periods of snow. The Soil Conservation Services Unit Hydrograph was used as the transform method, the constant monthly method was used for baseflow estimation, and the lag method was used for river routing. Automatic calibration was performed in the HEC-HMS model, which is used to minimize the objective function by fitting the best-possible agreement between the observed and the simulated flow at the gauging station [137]. The observed and simulated flows were compared at the outlet Marala barrage, as presented in Figure 1. The R 2 , Pearson's r, Nash and Sutcliffe efficiency (NSE), and root-mean-square error (RMSE) were used to compare the simulated flow with the observed flow, and equations are provided in Appendix A.

Numerical Simulation in the HEC-RAS Model
Two-dimensional floodplain inundation modeling requires fine-resolution topographic information about the river channel and floodplain areas [149] to get better control over the flow distribution [72,150]. For this purpose, the RTC elevation dataset (DEM) of the Phased Array type L-band SAR (PALSAR), mounted on the Advanced Land Observing Satellite-1 (ALOS), was used in the study. There were 96 cross sections, and their extraction was performed by converting the elevation model to a triangulated irregular network using the Delaunay triangulation interpolation method, which is a simple alternative of a regular raster and continuous surface with triangular facets. The extracted cross sections were corrected using Google Earth imagery to give the river a suitable and realistic shape.
The HEC-RAS v5 is a 2D numerical simulation model that be used to solve the 2D Saint Venant equation or 2D diffusive-wave equations [73,151,152]: where h is the depth of water (m), g is the acceleration due to gravity (m s −2 ), p and q are specific flow (m 2 s −1 ) in cartesian directions, n is the Manning resistance, ζ is the surface elevation (m), ρ is the density of water (kg m −3 ), f is the Coriolis term (s −1 ), and τ xx , τ yy , and τ xy are the components of effective shear stress [73]. The Coriolis term of the equation was disregarded in agreement with shallow-water assumptions and gravity flow conditions [153,154]. The inertial terms of the last 2 equations are neglected for diffusive-wave equation-based simulation. Both the Saint Venant and diffusivewave equations were applied in the HEC-RAS model. The computing cells were created within a closed-boundary polygon, with staggered grids of rectangular cells of similar size (12.5 m × 12.5 m) and polygons of different shapes at the corners. The implicit finite volume method was applied for solving the equations, which approximates an average integral on a reference volume, and hydraulic property tables were computed before solving the equations. Elevation-hydraulic property relationships for each computational cell face and elevation-volume relationship for each cell were computed. Finally, both the equations were solved for the 2D model, and 30 iterations were preferred as compared to the coupled 1D-2D model due to limited information about overflow locations. The inflow hydrograph at the Marala barrage was applied as a boundary condition, and flow was defined in the cells covering the river width at the first cross section. The normal depth boundary conditions were applied at the borders of the model.

Satellite-Based Flood Extent Mapping
The Landsat 5 TM, Landsat 7 ETM+, Landsat 8 OLI, MOD09GA, and MOD09A1 datasets were used to map the flood extent for calibration and validation of the HEC-RAS model. The digital numbers of the Landsat images were converted to reflectance using the Fast Line-of-Sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) module used for atmospheric correction [99,[155][156][157]. Landsat 5 (TM), Landsat 7 ETM+, and Landsat 8 OLI have 30 m spatial resolution, while MOD09GA and MOD09A1 have 500 m spatial resolution. All the satellite datasets were then reprojected to the Universal Transverse Mercator projection and resampled at 12.5 m spatial resolution (similar to the cell size of the DEM) using the nearest-neighbor method [97,98].
By considering the higher concentration of sediments in the floodwater, different NDWI equations were applied for water body delineation for all the satellite products. The NDWI equation [93] is between the green and NIR bands: The modified NDWI (MNDWI) can also be used as an indicator for water body delineation using green and shortwave-infrared (SWIR) bands [95,158,159], with a minimum MNDWI threshold of 0.0 [94]: The bands used for calculation of the NDWI, MNDWI 1, and MNDWI 2 and their thresholds are presented in Table 1. In general, the NDWI, MNDWI 1, and MNDWI 2 range from −1 to +1, and theoretical threshold values for all the indices greater than 0 indicate water bodies in the above equations [94,155]. However, thresholds greater than 0 are not possible for all the methods [160] and all the water bodies due to sensor quality, atmospheric interference, and characteristics of water bodies, and slight adjustment of the threshold is necessary for better results [95,161]. For this purpose, all the three approaches were used with threshold values indicated in the literature by different scientists, as presented in Table 1.  [168,169] The NDWI, MNDWI 1, and MNDWI 2 methods were used to delineate the blended high spatiotemporal flood extent using Landsat and MODIS products. All the three algorithms (NDWI, MNDWI 1, and MNDWI 2) were applied, owing to the high turbidity Remote Sens. 2021, 13, 3916 9 of 30 of floodwater, using thresholds recommended in the literature for flood extent delineation, and a blended flood map was developed by combining outputs of all the algorithms. The NDWI and MNDWI 1 were used for flood extent mapping using Landsat 5 (TM) for the calibration year 1992, and the blended flood extent was prepared for analysis. The NDWI, MNDWI 1, and MNDWI 2 were used to map the flood extent using Landsat 7 ETM+, Landsat 8 OLI, MOD09GA, and MOD09A1 products, and a final blended flood map was developed using outputs of all the products for the 2014 flood.

Calibration and Validation of the HEC-RAS Model
The HEC-RAS model was calibrated and validated for the flood events of 1992 and 2014, respectively, using different combinations of Manning's n values for the main channel and floodplain. Only Landsat 5 (TM) was used for the 1992 flood event, while different satellite images were blended to produce high-resolution spatiotemporal flood extent for the 2014 flood, as presented in Table 2. Manning's n values were selected for rivers and the floodplain using land cover data of the European Space Agency for 1992 and 2014. The roughness coefficients were further parametrized based on sensitivity analysis for the calibration year 1992. The model was calibrated and validated by comparing the inundation extent with the satellite-based inundation extent, and by comparing the simulated flow against the observed flow at barrages within the study area. Information about flow at barrages and satellite imagery used in the analysis for both the flood events are presented in Table 2. The simulated flood extent, produced by HEC-RAS against different values of Manning's n for the river and floodplain, was compared with satellite-based inundation extent for accuracy assessment, and to fine-tune Manning's n for the river and floodplain.

Accuracy Assessment of the HEC-RAS Model
The overall accuracy was estimated by comparing the correctly predicted flood area only (simulated cells) with the observed (satellite-based) flooded cells. It is the simplest way and provides accurate measures: where P 0 is the overall accuracy, n is the total number of flooded cells observed through satellite imagery, and C p are correctly predicted (simulated) flood cells.
The performance evaluation of the model was also tested using F1 and F2 calculations recommended by [170,171]: where A is the correctly predicted flood area, which is also observed in satellite imagery; B is the predicted flood area that is not actually observed in satellite imagery; and C is the area that is not simulated as flooded area by the model but is actually observed as flooded area in the satellite imagery. F1 and F2 values ranged from 0 to +1 and −1 to +1, respectively, where higher values (near to +1) indicated a good relationship between predicted and observed flooded cells.

Performance Evaluation of the HEC-HMS Model
The landuse data of the ESA (2014) and HWSD are presented in Figure 2. There is a major contribution of loam and clay soil in the upper and middle catchment of the Chenab River, and a high percentage cover of a glacier is also present. The higher percentage of these mentioned classes represents major flow contribution from any high-precipitation event that can augment the surface flow significantly. The landuse map represents major landuse, such as forest, grassland, mosaic cropland, bare land and snow cover. The grassland and glacier are in the upper-to-middle part of the catchment, while forest and mosaic cropland are from the middle-to-lower parts of the catchment. The landuse of the river catchment reveals that the catchment is mostly covered by grassland and mosaic cropland and the combined effect of these two landuse features can be helpful to minimize the flood intensity. The landuse of the ESA has been successfully used by researchers for hydrological modeling of river basins [172,173] where is the overall accuracy, n is the total number of flooded cells observed through satellite imagery, and are correctly predicted (simulated) flood cells. The performance evaluation of the model was also tested using F1 and F2 calculations recommended by [170,171]: where is the correctly predicted flood area, which is also observed in satellite imagery; B is the predicted flood area that is not actually observed in satellite imagery; and C is the area that is not simulated as flooded area by the model but is actually observed as flooded area in the satellite imagery. F1 and F2 values ranged from 0 to +1 and -1 to +1, respectively, where higher values (near to +1) indicated a good relationship between predicted and observed flooded cells.

Performance Evaluation of the HEC-HMS Model
The landuse data of the ESA (2014) and HWSD are presented in Figure 2. There is a major contribution of loam and clay soil in the upper and middle catchment of the Chenab River, and a high percentage cover of a glacier is also present. The higher percentage of these mentioned classes represents major flow contribution from any high-precipitation event that can augment the surface flow significantly. The landuse map represents major landuse, such as forest, grassland, mosaic cropland, bare land and snow cover. The grassland and glacier are in the upper-to-middle part of the catchment, while forest and mosaic cropland are from the middle-to-lower parts of the catchment. The landuse of the river catchment reveals that the catchment is mostly covered by grassland and mosaic cropland and the combined effect of these two landuse features can be helpful to minimize the flood intensity. The landuse of the ESA has been successfully used by researchers for hydrological modeling of river basins [172,173]. The comparison of simulated and observed flows using all the three precipitation products for calibration, validation, and application periods are presented in Figure 3a-c. The IMERG-F-based simulated flow was in good agreement with the observed flow in validation and application periods, with higher values of R 2 , Pearson's r, and NSE. The IMERG-F-based simulated flow indicated slightly lower values of Pearson's r and the NSE in the calibration period, while IMERG-F showed good performance in the peak flow simulation, which is obvious in the flood event of 2014 (Figure 3a). The TRMM_3B42RTbased simulated flow was also in reliable agreement with the observed flow in validation and application periods but with slightly lower values of R 2 , Pearson's r, and NSE as compared to the IMERG-F-based simulation. The TRMM_3B42RT-based simulated flow indicated a slightly fair agreement with the observed flow in the calibration period, and it did not perform well for simulation of the peak flow, which is clear in Figure 3b. The GSMaP_NRT-based simulated flow was in reliable agreement with the observed flow in calibration and application periods, with reliable values of R 2 , Pearson's r, and NSE, while it revealed slightly fair results in the validation period, and it did not perform well for peak flow simulation, as can be seen in Figure 3c. The scatter plots of all the simulated flows compared to the observed flow are presented in Figure 4a-c. The scatter plots covered all the calibration, validation, and application periods of the simulations. It is clear from the scatter plots that there were few deviations of points from the 1:1 line in IMERG-F-based simulation results as compared to the other simulation results. The scatter plots reveal that GSMaP_NRT-based simulation results indicate slightly lower deviations as compared to the TRMM_3B42RT-based simulation results, with slightly higher values of R 2 , Pearson's r, and the NSE.
A research on the same study area also revealed similar flow simulation results of IMERG-F (R 2 = 0.66, NSE = 0.61) and TRMM_3B42 v7 (R 2 = 0.64, NSE = 0.54) on the daily timescale [50]. The authors also concluded that GPM IMERG-F can be a good replacement of TRMM_3B42 v7. A hydrological study conducted by the researchers in the Beijiang River basin also revealed that IMERG-F presents satisfactory results, with a higher correlation coefficient (0.63) and NSE (0.742), with a higher probability of detection as compared to TRMM_3B42 v7, IMERG-Early, and IMERG-Late [174]. A study conducted by researchers in the Upper Huaihe River basin revealed that IMERG-F has better performance as compared to IMERG-E and IMERG-L [175]. A study conducted by researchers in Yunnan Province of China revealed that IMERG-F can capture more flood events and can be effectively used to improve the accuracy of flash flood warnings [176]. A study conducted by researchers in 300 river catchments across mainland China revealed that TRMM_3B42 has poor flow simulation performance, while IMERG-F has good results and it can be used for hydrological studies in ungauged river basins [177]. A study conducted by researchers in Myanmar watershed revealed that both IMERG-F and TRMM_3B42RT performed better in hydrological modeling, with higher NSE values of 0.84 and 0.868, respectively [178]. A study conducted by researchers in the Jhelum and Chenab River catchments revealed that GSMaP_NRT-based simulated flow underestimates peak flow events and it cannot be used for flood forecasting [179]. Reliable flow simulation performance of the GSMaP_NRT has also been observed by researchers in the Jhelum River basin, the Fuji River basin, and the sparsely gauged watershed in the Peruvian Andes [180][181][182]. From the results and discussions, it can be concluded that IMERG-F has better flow simulation performance than the other products used in the study, and IMERG-F can also be used for flood forecasting, while GSMaP_NRT and TRMM_3B42RT exhibit almost similar flow simulation results but both these products cannot be used for flood forecasting.

Performance Evaluation of the HEC-RAS Model
The landuse of the floodplain of the Chenab River catchment is presented in Figure 5   The maximum roughness coefficient is offered by forests, orchards, and sugarcane, but these classes reveal an insignificant percentage cover in the floodplain. The major contribution of landuse in the floodplain was by rice and cotton in 1992, in addition to fodders in 2014, and all these classes offered Manning's roughness coefficient of about 0.05.
The performance assessment of the calibration event is presented in Table 3, with an actual observed (satellite-based) flooded area of 3806 km 2 . Manning's n values for the channel and floodplain were selected, ranging from 0.030 to 0.060 with an interval of 0.005 and 0.01 for the channel and the floodplain, respectively. It was observed that the simulated flood extent increased by increasing the n value for the channel and floodplain. The increment in the simulated flood extent was observed more by a small increase in channel n, while the increase in the simulated flood extent was observed as low with a high increase in the n value of the floodplain. The lowest inundated area was recorded as 2914 km 2 , and the highest inundation extent was recorded as about 4155 km 2 . The A represents overlapped inundation of both simulated and observed data, and the maximum correctly simulated area was observed as 3480 km 2 at channel n of 0.060 and floodplain n of 0.04 and 0.05. The mentioned combinations of n values had a maximum overall accuracy of about 91.43% and higher values of F1 and F2, representing good and reliable results that have also been observed by researchers [170,171]. It was also clear that the values of F1 and F2 were higher at channel n of 0.04, but the overall accuracy was less at this value of n. Keeping in view all the approaches, a floodplain n of 0.05 and a channel n of 0.06 indicate the best fit for the calibration period, and the overall accuracy for all the tested combinations was more than 75%, representing reliable results. A research on the same study area also revealed similar flow simulation results of IMERG-F (R 2 = 0.66, NSE = 0.61) and TRMM_3B42 v7 (R 2 = 0.64, NSE = 0.54) on the daily timescale [50]. The authors also concluded that GPM IMERG-F can be a good replacement of TRMM_3B42 v7. A hydrological study conducted by the researchers in the Beijiang River basin also revealed that IMERG-F presents satisfactory results, with a higher correlation coefficient (0.63) and NSE (0.742), with a higher probability of detection as compared to TRMM_3B42 v7, IMERG-Early, and IMERG-Late [174]. A study conducted by researchers in the Upper Huaihe River basin revealed that IMERG-F has better performance as compared to IMERG-E and IMERG-L [175]. A study conducted by researchers in Yunnan Province of China revealed that IMERG-F can capture more flood events and can be effectively used to improve the accuracy of flash flood warnings [176]. A study conducted by researchers in 300 river catchments across mainland China revealed that TRMM_3B42 has poor flow simulation performance, while IMERG-F has good results and it can be used for hydrological studies in ungauged river basins [177]. A study conducted by researchers in Myanmar watershed revealed that both IMERG-F and TRMM_3B42RT performed better in hydrological modeling, with higher NSE values of 0.84 and 0.868, respectively [178]. A study conducted by researchers in the Jhelum and Chenab River catchments revealed that GSMaP_NRT-based simulated flow underestimates peak flow events and it cannot be used for flood forecasting [179]. Reliable flow simulation performance of the GSMaP_NRT has also been observed by researchers in the Jhelum River basin, the Fuji

Performance Evaluation of the HEC-RAS Model
The landuse of the floodplain of the Chenab River catchment is presented in Figure  5, which was used in the HEC-RAS model for inundation modeling. The landuse of the 1992 flood event was developed using Landsat 5 TM, while the landuse of 2014 was developed using Landsat 8 OLI. The landuse plays a key role in flood propagation in the floodplain of the river basins due to variation in the roughness offered by the earth features. Manning's runoff coefficients are also presented, along with the earth features in the floodplain of the study area, in Figure 5. The maximum roughness coefficient is offered by forests, orchards, and sugarcane, but these classes reveal an insignificant percentage cover in the floodplain. The major contribution of landuse in the floodplain was by rice and cotton in 1992, in addition to fodders  The results of numerical simulation of 2D diffusive-wave equations were compared with the Landsat 5 (TM)-based inundation extent for the 1992 flood, as presented in Figure 6a,b. All the approaches (NDWI and MNDWI 1) were used, and outputs were blended to produce a single better inundation extent. It is clear from Figure 6a,b that simulated and observed inundated extents represented almost similar patterns in each district. However, the simulated inundation extent overestimated the observed flood extent in the Gujranwala, Hafizabad, and Chiniot districts. A flooded area was also observed on the northern side of the Ravi River due to high magnitude of the water. The area between the Chenab and the Jhelum River on the upper side of the Jhang district was inundated due to the high flood in the Jhelum River, which further augmented the flow of the Chenab River at the Trimmu barrage. literature. The final observed area was recorded as 3489 km 2 after blending outputs obtained using different approaches for each satellite product. The overall accuracy ranged from 76.6% to 94.2%. The n values of 0.06 for the floodplain and 0.06 for the main channel were the best combination with higher values of F1 and F2.   The validation of the model was performed by selecting four higher values of the main channel n against the floodplain n values, as used in the calibration, as presented in Table 4. Landsat 7 ETM+, Landsat 8 OLI, and MOD09A1 and MOD09GA were used for the NDWI, MNDWI 1, and MNDWI 2, with upper threshold values recorded from the literature. The final observed area was recorded as 3489 km 2 after blending outputs obtained using different approaches for each satellite product. The overall accuracy ranged from 76.6% to 94.2%. The n values of 0.06 for the floodplain and 0.06 for the main channel were the best combination with higher values of F1 and F2.
The comparison between simulated and observed inundation extents for 2014 is presented in Figure 7a,b, indicating good agreement between inundation extents. There was good agreement in almost all the districts, but the model simulated more inundated area in the lower part of the Hafizabad district and upper parts of the Chiniot district lying near the river. The comparison between simulated and observed inundation extents for 2014 is presented in Figure 7a,b, indicating good agreement between inundation extents. There was good agreement in almost all the districts, but the model simulated more inundated area in the lower part of the Hafizabad district and upper parts of the Chiniot district lying near the river. Flood hazard maps indicate flood risk areas and provide the most advanced flood probability zones, which help people to decide about their investment in the floodplains.
The NDWI/MNDWI has been applied successfully using Landsat 8 OLI, TM, and ETM+ by researchers for water body delineation in different areas of China, the Evros River, wetlands of the Macquarie Marshes, and Varanasi (Banaras), India [183][184][185][186]. The spatiotemporal inundation mapping of the Murray-Darling Basin, Australia, was performed by researchers using the Open Water Likelihood algorithm for MODIS, and MNDWI for Landsat TM and ETM+ [187].
In this study, both the Saint Venant equation and the 2D diffusive-wave equation were applied. The results were similar, but the 2D diffusive-wave equation was much faster than the other, and similar findings have been observed by researchers [73]. A flood inundation study conducted by researchers for the 2002 Baeksan flood event in Korea re- Flood hazard maps indicate flood risk areas and provide the most advanced flood probability zones, which help people to decide about their investment in the floodplains.
The NDWI/MNDWI has been applied successfully using Landsat 8 OLI, TM, and ETM+ by researchers for water body delineation in different areas of China, the Evros River, wetlands of the Macquarie Marshes, and Varanasi (Banaras), India [183][184][185][186]. The spatiotemporal inundation mapping of the Murray-Darling Basin, Australia, was performed by researchers using the Open Water Likelihood algorithm for MODIS, and MNDWI for Landsat TM and ETM+ [187]. In this study, both the Saint Venant equation and the 2D diffusive-wave equation were applied. The results were similar, but the 2D diffusive-wave equation was much faster than the other, and similar findings have been observed by researchers [73]. A flood inundation study conducted by researchers for the 2002 Baeksan flood event in Korea revealed good agreement between the simulated inundation and observed data [70]. A flood inundation study by researchers in the Yesil (Ishim) River, Kazakhstan, revealed that there was good agreement between the observed (satellite-based) and simulated flood cells [188]. The HEC-RAS model has been used with reliable results by researchers for flood inundation modeling in the lower region of the Brazos River watershed, Prayagraj, Kilicozu Creek, and Secchia River [69,[189][190][191].
The hydrographs at barrages downstream of the Marala barrage were also compared for accuracy assessment in calibration and validation periods, as shown in Figure 8a There existed a strong relationship between simulated and observed discharges for both the time periods at barrages, with a minimum R 2 of 0.98. The simulated hydrograph was lower as compared to the observed hydrograph at all the barrages for both the time periods, because of the overestimated inundation extent by the model for both the time periods. It was also observed that at the Trimmu barrage, the reduction in simulated flow was even more as compared to the first two barrages, which indicates a continuous overestimation of the riverbanks' overflow by the model throughout the river section. The researchers conducted a study on the Grand River of Ohio, USA, which revealed that HEC-RAS is a suitable model for flood inundation modeling and river flow routing [192].

Severe Flood Risk Areas and River Routing for Early Flood Warning
The flood hazard areas were declared for the 2014 flood using criteria of the Ministry of Land Infrastructure and Transport (MLIT), Japan, as presented in Table 5. In H1 and H2 categories, on-foot evacuation is possible, but slight complication is present in the latter case for adults, which creates problems during evacuation. On-foot evacuation is not possible in H3 category, and people should stay in their homes. The areas lying under H4 and H5 categories should be monitored on top priority during a flood as these categories are extremely dangerous for people, even inside their homes. It is also clear from Figure 7a,b that the areas in the Gujranwala, Hafizabad, Chiniot, and Jhang districts came under the H4 category in the 2014 flood. In such a situation, people are insecure inside their homes and have to wait on the roofs of their houses for emergency rescue. The areas lying nearby the river are at high risk (H5 category), and in this situation, early evacuation and emergency rescue should be treated as top priority.
The HEC-RAS model also provides travel time and estimates the point velocity along with the depth at any specific location in the study area [193]. The observed (actual) travel time of water from Marala to Khanki barrage is 12 h, Khanki to Qadirabad barrage is 6 h, and Qadirabad to Trimmu barrage is 48 h [194], while the simulated time of water from Marala to Khanki barrage is 10.8 h, Khanki to Qadirabad barrage is 6.2 h, and Qadirabad to Trimmu barrage is 47 h. The simulated and observed travel times of water in the river are almost similar, thereby indicating reliability of river routing performed by the HEC- There existed a strong relationship between simulated and observed discharges for both the time periods at barrages, with a minimum R 2 of 0.98. The simulated hydrograph was lower as compared to the observed hydrograph at all the barrages for both the time periods, because of the overestimated inundation extent by the model for both the time periods. It was also observed that at the Trimmu barrage, the reduction in simulated flow was even more as compared to the first two barrages, which indicates a continuous overestimation of the riverbanks' overflow by the model throughout the river section. The researchers conducted a study on the Grand River of Ohio, USA, which revealed that HEC-RAS is a suitable model for flood inundation modeling and river flow routing [192].

Severe Flood Risk Areas and River Routing for Early Flood Warning
The flood hazard areas were declared for the 2014 flood using criteria of the Ministry of Land Infrastructure and Transport (MLIT), Japan, as presented in Table 5. In H1 and H2 categories, on-foot evacuation is possible, but slight complication is present in the latter case for adults, which creates problems during evacuation. On-foot evacuation is not possible in H3 category, and people should stay in their homes. The areas lying under H4 and H5 categories should be monitored on top priority during a flood as these categories are extremely dangerous for people, even inside their homes. It is also clear from Figure 7a,b that the areas in the Gujranwala, Hafizabad, Chiniot, and Jhang districts came under the H4 category in the 2014 flood. In such a situation, people are insecure inside their homes and have to wait on the roofs of their houses for emergency rescue. The areas lying nearby the river are at high risk (H5 category), and in this situation, early evacuation and emergency rescue should be treated as top priority.
The HEC-RAS model also provides travel time and estimates the point velocity along with the depth at any specific location in the study area [193]. The observed (actual) travel time of water from Marala to Khanki barrage is 12 h, Khanki to Qadirabad barrage is 6 h, and Qadirabad to Trimmu barrage is 48 h [194], while the simulated time of water from Marala to Khanki barrage is 10.8 h, Khanki to Qadirabad barrage is 6.2 h, and Qadirabad to Trimmu barrage is 47 h. The simulated and observed travel times of water in the river are almost similar, thereby indicating reliability of river routing performed by the HEC-RAS model. The simulated travel time enables administrative authorities to declare an emergency in severe-risk areas, and start emergency rescue and evacuation.

Implementation of Integrated Flood Management for Flood Mitigation in the Transboundary Chenab River Basin
Flood mitigation is a long-term and ongoing process before the occurrence of flood disasters, and it teaches us how to live rationally with floods. The World Meteorological Organization (WMO) and the Global Water Partnership established an Associated Programme on Flood Management (APFM) to develop a concept of IFM [195]. IFM aims at optimizing the use of floodplains and minimizing loss of life from flooding by integrating land and water resource development in a river basin [196]. IFM offers an underway paradigm shift from flood control to flood management by managing risk in living with floods [195], as there is a greater probability of floods in the future due to climate variability and change [197]. In Pakistan, the National Flood Protection Plan (NFPP)-IV (from 2005 to 2015) is being carried out using both structural and non-structural measures, with an amount of Rs. 290,919 and 41,327 million, respectively [198].

Clear and Objective Policies Supported with Legislation and Regulation
IFM suggests sound policies for planning, allocation, and management of available resources in river basins through a legislative framework by clearly defining powers, responsibilities, and obligations of the concerned departments and floodplain residents. There is a need to implement IFM in the Chenab River basin, with proper legislation and regulations to put IFM from concept into practice. The responsibilities and powers of all the departments and individuals should be fixed for pre-flood, flood, and post-flood conditions, for appropriate flood management in a coordinated way [199]. It is the responsibility of the WAPDA, FFC, FFD, NDMA, PDMA, and PMD to address the regulatory framework of IFM, which comprises water resources development, floodplain zoning, flood forecasting, flood early warning, and quick disaster response. The Army can be called in a flood emergency in H5 and H4 categories to help all the institutions responsible for the regulatory framework of IFM. NGOs can be permitted to use their available resources for the evacuation of people under H3 and H4 categories to save lives, which is the main aim of IFM. The situation is more complex in the transboundary Chenab River basin, and IFM should be implemented across the river basin if both the riperians (Pakistan and India) agree to it.

Adopting Basin Approach to Flood Management
Environmental degradation, urbanization, over-grazing, deforestation, and changes in farming systems have increased the probability of floods, and there is no effective landuse planning in the Chenab River basin to control floods in Pakistan [200]. IFM focuses on taking possible benefits from floods, such as retaining part of the floodwater in dams at upper catchments or at barrages in the lower reaches of the Chenab River. The purpose can be attained by strong coordination between India and Pakistan for storage and release of water from the dams. The purpose can also be attained in the lower reaches of the Chenab River in Pakistan, by compromising on agricultural lands just upstream of the barrages to declare storage areas for floodwater. Such storage areas at the barrages will not only reduce the flood intensity by providing temporary water storage, but will also augment the groundwater recharge due to alluvial floodplains of the Chenab River basin. Low-lying depressions nearby the banks of the Chenab River can also play a key role in flood attenuation, but dumping of waste in those depressions not only reduced the water storage facilities but also causes serious groundwater quality issues due to recharge. In the 2014 flood, areas in the Jhang district were badly affected, so IFM also considers alternate plans for extreme flood events to identify areas to be sacrificed for flooding away from the main cities without any sociopolitical influence.

A Multidisciplinary Approach
IFM focuses on multi-dimensional aspects of flood management in the Chenab River basin instead of traditional one-dimensional economic considerations. IFM considers the socioeconomic, environmental, and institutional aspects to determine the best-possible options to flood management by getting public input in the decision-making process. All the above-mentioned institutions, departments, and NGOs should work in good collaboration and under central coordination (DA), Government of Pakistan, in order to ensure that flood management activities are being carried out in a true manner, rather than converting a hazard into a disaster.

Manage Risk and Reduce Vulnerability
The livelihood of the majority of people living near the Chenab River (floodplain) is directly linked to farming, agriculture, and pastoral activities. IFM focuses on policy formulation for considering the flood risk for individuals living in poverty, with suitable socioeconomic implications for the future. Unfortunately, Pakistan has both extremes, high population growth and weaker economic conditions, which increases the flood risk. In these conditions, the main goal of IFM implementation is to prevent a hazard from becoming a disaster [201] in the Chenab River basin by strengthening livelihoods, providing easy access to resources, maintaining the ecosystem, enhancing flood preparedness, ensuring quick response, and pursuing sustainable recovery.

Enabling Community Participation
The community participation in flood management activities in the Chenab River basin can play a key role not only in building resilience in the communities but also in obtaining a better understanding of the sociopolitical barriers to flood management. The Punjab Province of Pakistan holds complex regional politics, so community participation can also bring input from all the cultural, religious, and regional communities to a single table, for a detailed discussion on IFM without any sociopolitical interference. Moreover, community participation decentralizes the flood management process and ensures equity in resources and development [201]. Gender and religion-based community participation can further bring the women, children, weaker persons of society, and representatives of different religions to the decision-making process.

Adopting the Best Mix of Strategies
Climate, basin characteristics, and socioeconomic conditions determine the bestpossible strategies under the framework of IFM. As climate is unpredictable, catchment is under Indian control, and socioeconomic conditions of Pakistan are weaker, so IFM compares the available resources and suggests a combination of workable strategies that are resilient, flexible, and adaptable to changing conditions. IFM also suggests the implementation of structural and non-structural measures appropriate to the region.

Structural Measures Proposed for the Chenab River
There are three commissioned dams under Indian control in Jammu and Kashmir (Salal, Baglihar, and Dul Hasti dams), but in Pakistan, there is no suitable site for construction of a major water storage facility. Therefore, Pakistan has no dam option to control/manage floods, while mutual effective communication can be the only possible way to control floods downstream by managing the storage and release of water at upstream. There are four barrages in the floodplain areas of Pakistan, which are being used to store water temporarily in order to reduce the intensity of floods in the river. New inundation canals from Qadirabad and Trimmu barrages can be constructed in order to extract the excess water from the river, store it in several small ponds in the fields, and also divert water to the Cholistan Desert of Pakistan. River training works are of prime importance to control potential floods in this contemporary time [202]. Straightening of the river either by relocation or by cutoffs has been widely used to increase the river capacity, thereby minimizing overbank flooding [203]. The assessment and approval of new flood control schemes and flood protection works are under the domain of the PID and FFC. At present, the Chenab River from the Marala barrage to the Chenab-Ravi Rivers' confluence has been subjected to meandering and following a zigzag manner, as presented in Figure 9a-h.
water temporarily in order to reduce the intensity of floods in the river. New inundation canals from Qadirabad and Trimmu barrages can be constructed in order to extract the excess water from the river, store it in several small ponds in the fields, and also divert water to the Cholistan Desert of Pakistan. River training works are of prime importance to control potential floods in this contemporary time [202]. Straightening of the river either by relocation or by cutoffs has been widely used to increase the river capacity, thereby minimizing overbank flooding [203]. The assessment and approval of new flood control schemes and flood protection works are under the domain of the PID and FFC. At present, the Chenab River from the Marala barrage to the Chenab-Ravi Rivers' confluence has been subjected to meandering and following a zigzag manner, as presented in Figure 9ah.
Apart from floods in this condition, the Chenab River is continuously expanding its breadth at off-structure sites and taking fertile land in it each year in Hafizabad, Chiniot, and Jhang districts. Artificial cutoffs should be preferred wherever needed, which will reduce overbank flooding from the Chenab River not only by shortening the river length but also by decreasing the resistance to flow [204]. The resulting higher velocity will further reduce the probability of overbank flooding not only by increasing the discharge but also by augmenting the river bed and bank erosion [205,206]. The confinement of flood flows over floodplains with the formation of levees and river training works, in the inner channel to transport high sediment load by guiding the flow away from the levees, save the river from scouring near the levees' foundation [207], as proposed in Figure 9a-d,g. At the curved portions of the river throughout the river length, the river flow velocity remains maximum at the outer periphery, which ultimately loosens the nearby fertile soil, and the river expands its breadth by cutting that soil into Apart from floods in this condition, the Chenab River is continuously expanding its breadth at off-structure sites and taking fertile land in it each year in Hafizabad, Chiniot, and Jhang districts. Artificial cutoffs should be preferred wherever needed, which will reduce overbank flooding from the Chenab River not only by shortening the river length but also by decreasing the resistance to flow [204]. The resulting higher velocity will further reduce the probability of overbank flooding not only by increasing the discharge but also by augmenting the river bed and bank erosion [205,206].
The confinement of flood flows over floodplains with the formation of levees and river training works, in the inner channel to transport high sediment load by guiding the flow away from the levees, save the river from scouring near the levees' foundation [207], as proposed in Figure 9a-d,g. At the curved portions of the river throughout the river length, the river flow velocity remains maximum at the outer periphery, which ultimately loosens the nearby fertile soil, and the river expands its breadth by cutting that soil into itself. Short spurs should be constructed where there are steep curves in the river channel, as proposed in Figure 9h. Deflecting and attracting groynes in Figure 9d-g can also be used as a tool to divert water away from the weakened banks of the river to avoid land loss, as proposed. The hockey and T-shaped spurs can also be constructed at the banks to train the river flow in a specified path to minimize overbank flooding, as presented in Figure 9e,f.

Non-Structural Measures Proposed for the Chenab River Basin
IFM suggests that structural measures and non-structural measures including afforestation offer partial safety to the people, and these strategies require a huge economic budget. Some of the non-structural measures constitute cost-effective/inexpensive multiple sub-plans for flood management, so it is crucial to focus on these measures. The IHP of the UNESCO recommended some non-structural measures for flood mitigation [208], as discussed for the Chenab River basin.

Flood Preparedness Measures
The public awareness regarding flood information, training, expected response, and community responsibility, should be raised as a top priority by the NDMA, PDMA, DDMA, DA, and NGOs through news channels, radio communication, social media, cell-phone messages and alerts, in situ information transfer at public places, announcements from speakers installed on towers/buildings, announcements in educational institutions, etc. It will be more beneficial and of utmost importance to seek out and declare honest and non-political representatives from the community for a small and specific area of risky floodplains. All these representatives should be well trained in advance and should have strong coordination not only with the district administration but also with every person from the community of their domain without any socioeconomic, political, and religious conflict of interests. The NGOs and local/social representatives from each district should arrange flood rescue and evacuation rehearsals in the floodplains in order to ensure peaceful and smooth execution of flood management plans during an emergency. The DA and social representatives from the community should ensure good condition of roads parallel to the floodplains, and the availability of vehicles, boats, sufficient food, security, and other material-related necessary items in each district should be ensured.

Emergency Response Measures
Loss of life and property can be minimized with effective response measures with a reliable flood forecasting and early warning system. There should be effective planning and implementation of disaster management plans in the Chenab River basin, which is the responsibility of the NDMA and PDMA in collaboration with other provincial authorities. One of the main issues is the poverty of people living in the floodplains, who ignore the flood warnings and try to save their precious belongings, get stuck in the flooding water, and often die or wait on the rooftops for rescue by the Army. The DA and social representatives should be vigilant regarding rescue operations in the affected areas, and the Army should be called for rescue in H4 and H5 categories, wherever needed. The NDMA, PDMA, DDMA, and DA, along with social representatives and NGOs, should assess and plan financial management and donations from federal/provincial ministries and foreign institutions for early rehabilitation of the affected communities.

Preserving Ecosystems
IFM considers the entire river basin ecosystem as a single unit and focuses on the sustainability of the ecosystem, which is the main principle of the United Nations Convention on Biological Diversity (1998). It is observed that structural measures for flood control and management are threats to the sustainability of the environment due to alteration in the river's natural environment. [209]. The proposed structural measures on the Chenab River may alter the flow regimes, fix the shape and geometry of the river, and separate the river from the natural floodplain, and can cause loss of habitats, threats to biological diversity, and decrease in productivity. The structural and some of the non-structural measures may adversely affect the riverine ecosystem, thereby disturbing a large number of flora and fauna by reducing the frequency of floods. Therefore, a trade-off is needed between competing interests to appraise the required flow regimes in a river basin for the sustainability of the riverine ecosystem necessary for food security and livelihoods in the Chenab River basin [210].

Conclusions
The hydrological modeling in HEC-HMS indicated that the IMERG-F precipitation product performed well in calibration, validation, and application periods as compared to GSMaP_NRT and TRMM_3B42RT, with higher values of R 2 , Pearson's r, and NSE. The flow simulation using low-latency precipitation products (GSMaP_NRT and TRMM_3B42RT) indicated reliable results with suitable values of R 2 , Pearson's r, and NSE. IMERG-F can be used for flood forecasting in data-scarce river basins, and it is also suitable for other fields of studies, such as agriculture water management, drought monitoring, groundwater recharge, etc. The GSMaP_NRT-based simulations are slightly good as compared to TRMM_3B42RT, and both these products can be used for other fields of studies, but cannot be used for predicting the extreme hydrological response of complex river basins. Both the finer resolution landuse developed using in situ data and detailed in situ soil texture information can be helpful in predicting the hydrologic response in a more reliable way. The in situ snow and precipitation measurements can be used to correlate with satellite precipitation data before hydrologic simulations, in order to decide which precipitation product can be used effectively. The HEC-HMS model in this study is used as a semidistributed model, but a fully distributed model can be helpful to ensure better results. The calibration and validation results of HEC-RAS indicated that Manning's n value of 0.06 is optimum for inundation modeling for the Chenab River and its floodplain. The overall accuracy, at Manning's n value of 0.06 for the river and floodplain, was observed as 91.43% and 94.2% for the 1992 and 2014 floods, respectively. The detailed subsurface soil texture information of the Chenab River basin can be helpful to increase the accuracy of the model. The current landuse in the study area includes almost all the major crops, including some fodders, which are not suitable for the floodplain areas. Flood propagation can be minimized in the floodplain, but it requires massive changes in the landuse/agricultural measures, such as mulching, direct seeding, cover crops, perennial crops, and conversion of arable fields into grasslands. The detailed in situ information about river banks' vegetation can be helpful in assigning exact roughness values to the banks in the model, which will increase the efficiency of the HEC-RAS model. The comparison of simulated and observed hydrographs at the barrages indicated good relationship, which forecasts the traveling time across the river length, thus provides guideline to administrative authorities for better planning of evacuation and rescue. River training works are also proposed for improving the river and minimizing the probability of overbank flooding. The multi-disciplinary approaches of IFM proposed in the Chenab River basin can be helpful for effective flood management. Funding: This research has not received any external funding.

Acknowledgments:
We acknowledge financial support from the Deutsche Forschungsgemeinschaft and Universität Rostock within the funding program Open Access Publishing. We are thankful to the Water and Power Development Authority of Pakistan for sharing the Chenab River flow data. We are also thankful to the Federal Flood Commission of Pakistan for sharing the flood information.

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

Appendix A
where, F o = observed flow. F s = simulated flow.