An artificial neural network‐hydrodynamic coupled modeling approach to assess the impacts of floods under changing climate in the East Rapti Watershed, Nepal

Recurring floods have devastating consequences on the East Rapti Watershed (ERW), but effective mitigation/adaptation measures are lacking. This article aims at establishing a rainfall‐runoff (RR) relationship; estimating depth and extent of inundation under climate change scenarios; assessing impacts on the socio‐economy; and identifying and evaluating adaptation strategies in the ERW. Artificial Neural Network (ANN) was used to generate peak flows which were then entered into a hydraulic model to simulate inundation. Results were validated with field survey. The calibrated and validated RR and hydraulic models were fed with projected future climate (2021–2050) derived from multiple regional‐climate‐models to assess the changes in inundation. Results showed the peak discharge likely exceeds 10,500 m3/s at the ERW outlet in the extreme future flood scenario with corresponding inundation of 80 km2 and up to a depth of 11 m sweeping away over 1000 houses and 19 km2 of agricultural land in the critical areas. Constructing a 17 km long embankment in the critical areas along the right bank of the East Rapti River could reduce the flood spread by 35%, safeguarding 78% of the houses and saving 51% agricultural land compared with the scenarios without the embankment.

damages. Moreover, they provide evidence-based information for decision-making on flood risk reduction and devising adaptation strategies (Schwarz et al., 2018). Onedimensional (1D) and/or two-dimensional (2D) models are used to simulate flood propagation depending upon the requirement, data availability and site conditions. 1D models are generally used for longer simulation periods in rivers where the flow is concentrated along the main river whereas 2D models find their application in areas where flow is likely to change direction and spread across a large area (Hutanu et al., 2020). Additionally, if only peak flow data is available, 1D steady flow modeling is preferred while 2D unsteady modelling is performed if continuous flood flow hydrographs are available as input to the model (HEC-RAS, 2020). Some popular hydrodynamic modelling software are HEC-RAS (Song et al., 2019), MIKE 11 (Chowdhury et al., 2020), FLO-2D (Mitra et al., 2021), LISFLOOD-FP (Rajib et al., 2020), etc. Among the many available models, HEC-RAS developed by U.S. Army Corps of Engineers has been widely used for flood hazard studies due to its flow modelling capability in 1D, 2D, and combined 1D-2D domain (Aryal et al., 2020;Song et al., 2019). In Nepal, several flood studies have applied 1D HEC-RAS in rivers such as Babai (Bhattarai et al., 2019), Kankai (Khadka & Bhaukajee, 2018), Daraudi (Banstola & Sapkota, 2019) and West Rapti (Devkota et al., 2020), to list a few.
Hydrodynamic models require flood hydrographs as the basic input. Hydrological models are generally used to generate these flood hydrographs at points of interest along the stream. Physically based distributed and conceptual models can predict the rainfall-runoff (RR) process at the catchment scale fairly accurately. However, they are not applicable to basins where detailed data and parameters are unavailable . In such cases, data driven models like Artificial Neural Network (ANN) can be used to optimally model the RR process in catchments even with insufficient input data. Such models have been found to be handy where runoff is the only output needed by developing empirical relations between the input (for example, precipitation and temperature) and output (water level and discharge) variables without considering the physics behind the RR process Wu et al., 2021). Due to such an ability, application of ANNs have increased over the years for RR modeling (Wu et al., 2021).
Climate change is expected to impact the hydrology by increasing the intensity of extreme precipitation events, leading to higher risk of flooding (Tabari, 2020). A study by MoFE (2019) projected increase in mean temperature of Nepal by 0.9 to 1.1 C and 1.3 to 1.8 C in the 2030s and 2050s, respectively, as a result of CC. Additionally, precipitation is projected to increase by 11% to 23% by 2100. Various studies have modeled RR responses to CC (Hakala et al., 2018;Skoulikaris et al., 2019). Outputs from Global Climate Models (GCMs) or Regional Climate Models (RCMs) are corrected for biases while projecting future climate of an area. Among the various methods available for bias correction, Quantile Mapping (QM) is considered as a relatively better technique for improving the performance of climate models (Tegegne et al., 2017). Studies like Dhaubanjar et al. (2019) and Marahatta et al. (2021) have applied this technique in CC studies in different basins of Nepal. However, climate projections are approximations which vary considerably with location and physiographic regions, thus, warranting careful interpretation of their results.
Climate-induced flood risks can be minimized by suitable adaptation strategies that are cost-effective and environment-friendly. Adaptation options can be structural and nonstructural (Martinez et al., 2021). Structural adaptation measures include flood embankments, culverts and improved drainage, among others; embankments are the most popular (Ward et al., 2017). Nonstructural measures like national adaptation plans, awareness raising, hazard and vulnerability mapping and zoning, community-based adaptation plans and land use planning have been used in different studies (Devkota et al., 2016;Devkota & Bhattarai, 2015;Gonzalez et al., 2016;Martinez et al., 2021;Thapa-Parajuli et al., 2018). In Bangladesh and Nepal, awareness raising, early warning systems and communitybased flood management have been seen to be quite successful (Dewan, 2015). Though there are various potential adaptation strategies documented in literature, suitability to a particular area is largely governed by its cost-effectiveness and preferences of the local community (Nofal & Lindt, 2020).
Despite abundant recent studies on future hydroclimatic projections due to CC (e.g., Bhattarai et al., 2022;Dhaubanjar et al., 2019;Marahatta et al., 2021), their impacts on flood inundation at the watershed scale are still limited in Nepal. Consequently, very few studies (Dewan, 2015) have quantified the flood impacts on the socio-economic sectors. Studies testing and recommending both structural and nonstructural flood adaptation strategies are scarce in Nepal. Furthermore, sparse meteorological stations pose an additional challenge in flood modelling. Majority of the past flood studies have used conventional hydrological models (e.g., HEC-HMS, SWAT, etc.) to simulate the RR process in ungauged river basins of Nepal despite difficulties in validation because of lack of sufficient observed data. However, recent global trend of RR modeling is the application of data-driven techniques such as ANN as a convenient alternate compared to the conventional approach (Meresa, 2019;Saez et al., 2018). Thus, this study tries to fill this research gap in Nepal by utilizing data-driven technique for RR modeling of ungauged/sparsely gauged watershed under changing climate, assessing socio-economic impacts of inundation, and proposing both structural and nonstructural measures for the flood-impacted communities. Therefore, the overarching objective of this study is to assess the impact of CC on inundation and its implication on the socio-economy of one of the data scarce watersheds-East Rapti Watershed (ERW)located in south-central Nepal. This study aims to answer the following research questions: (i) What will be the extent and depth of inundation in ERW under current and future climatic scenarios? (ii) What are the potential impacts of inundation on the socio-economic system of ERW?, and (iii) What are the possible adaptation strategies to minimize such impacts? The results will definitely be beneficial to the ERW in particular while the applied methodological framework consisting of the ANN-hydrodynamic modeling will find its applicability in other partly gauged flood-prone basins of the region.

| Study watershed
The ERW has a drainage area of 3202 km 2 (Figure 1a). In a study of 25 river basins of Nepal, out of the 41 extreme flood events that have occurred in those basins in between 1991 and 2015, 40% were in ERW (ADPC, 2016). ERW has reported 143 flood events that affected 3460 people and 139 died in these events during the same period. Among these, the most devastating 10 events reported during 1991-2015 caused damage of over one million Nepalese Rupees (NRs.) in 2015 prices. These evidences indicate that the ERW is a flood-prone basin which requires appropriate flood management and adaptation strategies for risk minimization.
The elevation of the watershed varies from 136 to 2579 m above the mean sea level. Although the ERW lies immediately below the Nepal Himalayas, it is not a snowfed watershed. There are seven meteorological and three hydrological stations administered by the Department of Hydrology and Meteorology (DHM), Government of Nepal within the ERW. The average annual rainfall of the basin is about 2008 mm. Similarly, the average maximum and F I G U R E 1 Location map of East Rapti Watershed (ERW) (top); land use/cover map of ERW (bottom) (ICIMOD, 2013). DHM: Department of Hydrology and Meteorology minimum temperatures are 30 C and 17 C, respectively. Out of the three gauging stations, Rajaiya (DHM Index no. 460) is located in the uppermost part of the basin (Figure 1a). The observed instantaneous maximum discharge at this station is 1010 m 3 /s in the baseline period . The most dominant land use/cover (LULC) of ERW is forest (66%) while agricultural land covers about 28% of the catchment area ( Figure 1b).

| Methodology
The methodological framework used in this study is shown in Figure 2. The study started by developing a rainfall-runoff (RR) model applying ANN to estimate daily monsoon (June-September) discharge using precipitation and temperature as the input variables. The daily discharge was converted into instantaneous maximum discharge by using regression equation developed from the plot of yearly 24-h maximum discharge versus yearly instantaneous discharge available at Rajaiya for the baseline period. The next step included flood modeling using HEC-RAS to calculate the flood extent and depths along the river for the baseline period followed by inundation mapping. The study watershed was divided into four sections ( Figure 3a). Then, future climate data was derived from climate models and fed into the flood model considering five flood scenarios (Figure 3b) to assess the impact of CC on the ERW.
Details of the data used in this study are shown in Table 1. The historical hydro-meteorological data was obtained from DHM. Future climate data was derived from CORDEX-SA. The Land-use/cover (LULC) map was obtained from ICIMOD (2013) and the SRTM 90 m digital elevation model (DEM) was obtained from USGS.

| Rainfall-Runoff (RR) modelling
Unavailability of a gauging station at the basin outlet made the application of conventional hydrological models such as HEC-HMS and SWAT, among others, difficult. Hence, data driven ANN model was used to establish the RR process of the study basin. ANN model was trained, tested and validated using daily rainfall and temperature as input data and discharge as the output for the monsoon period. The same model was used to estimate the flows at downstream confluences of the East Rapti River with the other smaller tributaries assuming a similar RR pattern throughout the study watershed.
Mathematical representation of the ANN is given in Equation 1.
where Q, discharge; t, t À 1, t À 2 = day considered as t, the day before t as t À 1, the day before t À 1 as t À 2. R t , R tÀ1 , R tÀ2 = Weighted average rainfall (calculated using Thiessen polygon method) from seven meteorological stations at day t, t À 1, and t -2, respectively. T max = Average of maximum temperature of meteorological stations 902 and 906 (Figure 1a for location). Multi-Layer Perceptron (MLP) with feed forward back propagation is a popular architecture in RR modelling for simulating base flow, dynamic RR process and river flow simulation (Dastjerdi et al., 2019). Among the many available learning rules for the neutral network, delta learning rule is faster and accurate in deep learning and therefore has been used in many recent studies (Shaikh, 2020). To introduce nonlinearity into the output of a neural, variants of activation function can be used. In recent literatures, Rectified Linear Unit (ReLU) are used as activation fuctions for a neuron in MLP (Hawamdeh & Kuisi, 2021).
MLP usually consists of three layers: the input layers, where data are introduced to the network; hidden layer, where the data are processed (can be one or multiple) and the output layer, where the results for given inputs are produced (Saez et al., 2018). Selection of appropriate number of hidden layers depends on the complexity of both input and output parameters and has significant role in the quality of learning. Previous studies have shown that two hidden layers is sufficient to model complex rainfall-runoff process (Saez et al., 2018). Therefore, in this study, MLP architecture ( Figure 4) with one input layer, two hidden layers, and an output layer trained with a back-propagation algorithm is used to predict the  discharge. The hidden layers were trained using delta learning rule and ReLU as transfer function in python.
The ANN model was trained using monsoon data (June to September; calibration from 2004 to 2009) at daily resolution and then validated for the monsoon season of 2010 to 2013. Six indicators, namely, Nash-Sutcliffe Efficiency (NSE), coefficient of determination (R 2 ), index of agreement (d), standard deviation ratio (RSR), percent error in peak flow (PEPF), and percent bias (PBIAS) have been used to access the performance of the ANN model which is shown in Table 2 (Moriasi et al., 2007). The calibrated ANN model was used to generate daily discharges at the points of interest along the river reach for flood modelling.

| Flood frequency analysis
A regression analysis between the 24-h mean discharge and the yearly instantaneous maximum discharge available at Rajaiya for the baseline period (1979-2010) was carried out. The same equation was used to estimate the instantaneous maximum discharge for the future period (2021-2050) assuming that flood characteristics of the basin remains constant in the future too. Furthermore, flood frequency analysis was performed assuming Gumbel distribution to determine the flood magnitudes of 2, 5, 10, 25, 50, and 100 years return periods.

| Hydraulic modelling
One dimensional steady flow simulation was performed to calculate the surface profile under a mixed flow regime in HEC-RAS (HEC-RAS, 2020) and flood extent maps were prepared using HEC-GeoRAS. During model calibration in HEC-RAS, the values of Manning's n (0.05, F I G U R E 4 Structure of multilayer feed forward network of the ANN used in this study (adapted from Saez et al., 2018) T A B L E 2 Performance ratings adopted for the ANN model in this

<±10
±10 to ±15 ±15 to ±25 >±25 Note: *Performance rating criteria adopted from (Moriasi et al., 2007); O i = observed discharge at the ith time; O = mean observed discharge of the total simulation time period; S i = simulated discharge at the ith time; S = mean simulated discharge of the total simulation time period; n = number of timesteps in the total simulation period; Q 0 peak ð Þ = observed peak during the total simulation period; Q s peak ð Þ = simulated peak during the total simulation period; V 0 = Total observed flood water volume; V s = total simulated flood water volume. 0.04, and 0.035 for left over bank, main channel and right over bank, respectively) were assigned based on the recommendations of Gautam and Dulal (2013) for the study basin. In most steady 1D flow computations, downstream boundary condition is set as normal depth (i.e., uniform depth) because it allows the user to provide an energy slope and automatically calculate the depth using Manning's equation (Banstola & Sapkota, 2019;Bhattarai et al., 2019). Thus, downstream boundary condition was provided as normal depth in this study. Generally, 5-and 100-year return period floods are considered as frequent and rarer events, respectively in flood studies (Hawker et al., 2020). In this study, simulation was carried out in HEC-RAS for 5-and 100-year return period flood scenarios in both the baseline and future periods for comparision. As the ERW is only partly gauged and continuous instantaneous maximum discharge for past years were unavailable, recent available instantaneous maximum flood discharge data of the 12th July, 2019 flood was considered for validation. Initially, an attempt was made to validate the model by making pixel-to-pixel comparisons of the flood maps with LANSAT 8 (250 m resolution, scene ID = "LC81410412019195LGN00", accessed date: August 15, 2021) and MODIS (500 m resolution, MOD13Q1, accessed date: August 14, 2021) satellite data (Enea et al., 2018). But it could not be done because of lack of error-free satellite data for this region for the selected flood event. Also, Sentinel satellite imagery was explored but no data were available for download for the selected region and the flood event. Hence, a field survey method of validation was adopted in which the inundation map as a result of the July 12, 2019 flood (i.e., 2350 m 3 /s) in the East Rapti River was generated by the model and verified by visual inspection of the watermarks and based on consultation with the local communities. A similar validation approach has been used satisfactorily in other parts of Nepal especially when observed data are not available or very scarce (Devkota et al., 2020;Sarchani et al., 2020). Field survey was conducted by visiting accessible sites (outside of the national park area which is under high security) during the COVID 19 lockdown period in September 2020 by the study team led by the first author.

| Future climate projection
Bias corrected nine RCMs using QM under two RCP scenarios (8.5 and 4.5) recommended by (Ray et al., 2022) were used to project future climate for the study basin. The projected data was input to the RR model to generate future daily discharge which was used in the hydraulic model simulation to generate inundation maps. Two RCMs (out of the shortlisted nine) for which the flood peaks of RCP8.5 were larger than RCP4.5, were selected. From the prepared inundation maps, the flood value corresponding to the maximum flood extent was considered as the extreme flood condition.

| Impacts of floods on socio-economy
To assess the potential flood impacts of CC in the future, inundation depth was divided into four classes (low: <1 m depth; medium: between 1 and 2 m; high: between 2 and 3 m; danger: >3 m) based on the severity. Usually inundation of depth <1 m would be easy for adults to wade through, and thus there would not be any lifethreatening damage potential. Flood depths of more than 1 m pose a significant risk to human life and property. While, single-storied houses get inundated and could be washed away (which is the case for most rural houses) if flood depth exceeds 3 m. Therefore, it was found rational to consider these classes of inundation depth for our study. Similar values have been taken in several other studies of Nepal and also recommended by the Rainfall-Runoff manual of Australia (Ball et al., 2019;Bhattarai et al., 2019).
Additionally, impact of floods on the livelihood of the community was assessed by estimating the number of households submerged, inundated land and agricultural production of the study basin. These were derived for the baseline and future from the inundation map of ERW in the extreme flood conditions. Consequently, the direct economic loss due to inundation of summer crops (June-September) in the baseline and future periods were determined by multiplying the market rate of the crops per hectare and total agricultural land inundated. Critical areas were then marked based upon the maximum number of houses and largest agricultural land inundated in the extreme flood condition. The economic loss in those critical areas as a result of historical and future inundation scenarios were compared.

| Evaluation of adaptation strategies
To deal with the potential flood damages, various structural and nonstructural adaptation strategies were reviewed (please refer to Section 1: Introduction). The most commonly used structural measure for disaster risk reduction is embankment due to its cost-effectiveness in areas with high population and asset concentrations (Ward et al., 2017). Therefore, the effect of addition of embankments in the critical areas on the right bank of the mainstream was simulated for extreme flood condition in future using the HEC-RAS model. The total number of houses that can be safeguarded, agricultural land reclaimed and the corresponding economic value of the assets due to the construction of embankments in critical areas in the projected future climatic conditions were determined. Flood forecasting and warning, flood risk management, public participation, and other nonstructural approaches are frequently employed to mitigate flood risk across the world (Gonzalez et al., 2016;Martinez et al., 2021). Some viable nonstructural adaptation measures which might be suitable for the study area are recommended at the end derived from literature review.

| ANN model performance
Details of the ANN model calibration, simulation results and model performance are shown in Table 3. The statistical evaluation for calibration and validation phases showed satisfactory performance of the model for all the indicators (Table 2) as recommended by Moriasi et al. (2007). It can be seen that Nash-Sutcliffe Efficiency (NSE) ranges from 0.55 (0.7) to 0.75 (0.9) in calibration (validation) phase. Similarly, the variation of coefficient of determination (R 2 ) in calibration (validation) is from 0.56 (0.74) to 0.79 (0.91). In addition, index of agreement (d) in calibration (validation) ranges from 0.84 (0.93) to 0.92 (0.98). Similarly, standard deviation ratio (RSR) in calibration (validation) varies from 0.5 (0.31) to 0.67 (0.55). Likewise, the range of percentage error in peak flow (PEPF) in calibration (validation) is from À23.95% (6.31%) to +25.47% (+29.79%). Correspondingly, percentage bias (PBIAS) in calibration (validation) phase ranges from 0.36 (0.29) to 13.48 (8.15). Interestingly, the validation results were found to be better because of the good correlation between rainfall and runoff than for calibration ( Figure 5).

| Regression analysis of flood peaks
The regression equation derived from the plot of maximum of the yearly 24-h mean discharge versus instantaneous maximum discharge available at Rajaiya for the baseline period is shown in Figure 6. The logarithmic regression equation with R 2 of 0.88 shows a good prediction capability. Thus, the same equation was used to predict the instantaneous maximum flow in the future period as well.

| Validation of hydraulic model
During the field survey, high water marks of the July 12, 2019 flood event was observed and their depths were measured at six locations, this flood being the most F I G U R E 6 Yearly instantaneous maximum versus yearly maximum 24-h mean discharge recent in the ERW. In addition, local people were consulted for cross-verification of this information. The observed and simulated water depths were compared with the model simulated inundation map (Figure 7). The problem of inundation is significant in the downstream part of the ERW where major portion of the Chitwan National Park is located. Due to high security in this protected area and unavailability of access road in downstream part of the river, the observation points were limited to six locations. In the uppermost validation point (Rajaiya gauging station), it is seen that the simulated depth (4.39 m) and observed depth (4.3 m) are very close. As flood spreads in the downstream sections, the difference between the observed and simulated flood depth becomes significant with the maximum deviation of 0.24 m at Kumroj (immediate upstream location of riverside Sauraha). Thus, a good agreement is seen between the observed and simulated values at the considered locations indicating the model is well validated.

| Future climate and discharge
It can be seen from Table 4 that the future maximum temperature for RCP4.5 (RCP8.5) is projected to increase by 6.3% (7.5%) as compared with the baseline period. Also, the future minimum monsoon temperature is projected to decrease by 0.6% in both RCP scenarios. Likewise, maximum monsoon rainfall is projected to increase for RCP4.5 (RCP8.5) by 25.4% (29.7%). Two RCMs (MPI_REMO and NOAA_RegCM4) show increase in RCP8.5 flood peaks compared with the RCP4.5 case for most of the return period floods. Generally, high emission scenarios lead to higher values of mean daily (and corresponding instantaneous maximum) discharge and therefore, MPI_REMO and NOAA_RegCM4 were selected for hydraulic modelling. The peak discharge in RCP8.5 is likely to double from the RCP4.5 value in the case of RCMs MPI_REMO and increase by 60% for NOAA_RegCM4 for a 100-year return period flood.

| Climate change impacts on inundation
The change in inundation extents in the future due to CC corresponding to the two selected RCMs with respect to the baseline period for the considered five flood conditions are summarized in Table 5. Flood condition I indicates flood in East Rapti mainstream only while flood condition V indicates simultaneous floods in the mainstream as well as all its tributaries. It is evident that a 5-years (100-years) return period flood is expected to cause a decrease (decrease) in the inundation extent by 1% (1%) for flood condition IV of NOAA_RegCM4 under RCP4.5 and increase by 17% (13%) for flood condition V of MPI_REMO under RCP8.5. The maximum inundation extent (79.85 km 2 ) is expected to occur in the future for a 100-years flood, RCP8.5, flood condition V of MPI_R-EMO. Considering the 100-years flood in RCP8.5 for condition V of MPI_REMO, the inundation extent is expected to increase by 17% for the "high" inundation level (depth <3 m) and by 3% for the "danger" inundation level (depth >3 m) compared with the baseline period (Table 6). Additionally, Figure 8a shows that inundation level for the ">3 m" depth class is mostly confined within the channel boundary in the upstream section and spreads towards the banks in the downstream section. From the inundation map prepared for baseline and future extreme flood conditions, the total number of inundated houses were estimated to be 679 and 1173 and agricultural land covering 29 and 34 km 2 , respectively, with the maximum inundation depth of over 11 m in the river channel.

| Socio-economic implications of inundation
Based on the number of inundated houses and the extent of agricultural land, four impacted areas were marked and classified as uppermost, upper, lower and lowermost critical areas. In the baseline and future periods, 677 and 1073 houses while 17 and 13 km 2 of agricultural land respectively fall under critical areas. Considering the average number of people per household to be 4.5 (CBS, 2016), total number people affected by inundation in the baseline and future periods are respectively 3056 and 5278. Similarly, if the economic gain from agricultural land is considered as NRs. 55,000/ha (MoALD, 2014), direct economic loss of NRs.145 million and NRs. 170 million can be expected, respectively, under extreme flood conditions. This is a huge loss for a developing country like Nepal.

| Evaluation of structural adaptation strategies
To simulate how potential flood damage in the future could be minimized, embankments were added to the right bank of the East Rapti River in the hydraulic model in critical areas where settlements are dense ( Figure 8).
As the length and height of embankments varied in each critical area, maximum and minimum heights of  (Table 7). Thus, about 17 km long embankment with maximum height ranging from 4 to 11 m were provided in different critical areas. The hydraulic model was then run for the extreme flood condition scenario. The flood inundation maps focused on the critical areas are presented in Figure 8. Addition of embankments resulted in reduction of flood extent with increased flood depths as shown in Figure 9. The future flood extent reduced by 35% because of the embankment.
To estimate the average cost of the embankment, average of maximum heights of embankment, that is, 7 m (average of 4 and 11 m) in all critical areas were considered. The average cost of an Reinforced Cement Concrete (RCC) embankment is Nepalese Rupees (NRs 1 ) 20,000 per cubic meter in the study basin (DOR, 2018). If a cantilever RCC embankment 17 km long having top width 0.5 m and base width 3.6 m (0.5 x average height of embankment) is considered, then the total cost of construction will be about NRs. 5 billion. The addition of embankment can reclaim a total of 1579 ha land, which is worth about NRs 9.3 billion (DOLMA, 2021). In addition, 832 houses with 3744 affected people and 6 km 2 of agricultural land equivalent to NRs. 33 million worth of direct agricultural produce can be saved in the future with the construction of embankments. Details of the flood extent, numbers of houses and agricultural land inundated in the baseline period and that can be saved due to the construction of embankment in the critical areas in the future are shown in Table 8. Compared with the baseline, flood extent, number of houses likely to be inundated and agricultural land are expected to increase in the critical areas by 28% 58%, and 44%, respectively, in the future due to CC. If embankment is constructed during the baseline period as well, significant cultivable land and houses can be saved as in the future.

| DISCUSSION
In this study, the ANN model configuration is based on five inputs (R t , R tÀ1 , R tÀ2 , T max , T min ) for a single output (Q) which showed satisfactory results in calibration and validation phases. Due to limited hydro-meteorological data within the study area, variation in the input (e.g., with the addition of humidity, wind, sunshine) were not possible. However, the ANN model could be set for multiple set of inputs to enhance the performance of the model. Rajurkar et al. (2002) compared the ANN configurations for three cases (single input, two inputs and three inputs) and concluded that the performance capacity of ANN model increases as the number of inputs increases. Also, ANN model performance in this study shows that it can be used in watershed with no gauging stations at the outlet when the hydrological response does not change much throughout the watershed. This is comparable to the studies by Meresa (2019) and Kassem et al. (2020) which revealed that ANN approach is suitable to predict runoff in ungauged catchment with a reasonable accuracy. We have assumed that the value of peak discharge at any day t is a function of the maximum rainfall for consecutive 3 days (R t , R tÀ1 , R tÀ2 ) and temperature (T max and T min ) of the flood day t. As a result, high total monsoon rainfall and high average monsoon temperature for any RCM may not always necessarily result in higher peak discharge.
Validation of the hydraulic model in ungauged river basins is quite challenging. When water level and/or discharge gauging stations are available at the basin outlet or confluence with another river downstream, direct validation of the model is done using the observed data. However, in the case of ungauged or partly gauged basins, the model is indirectly validated by either satellite image or field survey. A recent study by Tamiru and Dinka (2021) evaluated HEC-RAS model by comparing inundation map generated in HEC-RAS and water body delineated from remotely sensed LANDSAT 8 imagery in lower Baro Akobo River basin in Ethiopia. Similarly, several other studies evaluated flood inundation maps from MODIS and ALOS satellite images in flood affected areas in Myanmar and Pakistan (Amarnath & Rajah, 2016;Khaing et al., 2019). In our case, while an attempt was made to perform satellite image validation by comparing LANDSAT 8 and MODIS data with the model simulated results, due to high cloud cover (above 80%) validation was not possible. As a result, this study resorted to field-based indirect validation. Similar postflood field measurement-based validation have been previously applied in poorly gauged river basins in Greece (Sarchani Results of our study show that the future monsoon climate is projected to be warmer and wetter due to increased average monsoon temperature and rainfall. On the other hand, the maximum monsoon rainfall is projected to decrease. Additionally, the frequency of occurrence of rainfall is expected to increase in the monsoon. These can be attributed to the increase in number of rainy days with rainfall of low intensity. Such results are highly dependent on the choice of the RCMs and their underlying uncertainties. Similar findings have been reported by recent flood studies in Thuli Bheri (Aryal et al., 2022), and Karnali basins of Nepal (Dhaubanjar et al., 2019) underscoring the high uncertainties related to the choice of climate data, hydrological model and the time horizon of analysis.
The instantaneous maximum peak flood is projected to increase in future resulting in increased flood extent in different climate scenarios compared with the baseline. Also, flood extent of low inundation depth decreases and high depth increases as flood magnitude increases in the future. Similar results are depicted in southern Italy (Apollonio et al., 2020), South China (Li et al., 2019), West Rapti (Devkota et al., 2020) and Karnali River basins (Aryal et al., 2020) in western Nepal. Moreover, we can infer from the results that plain areas with higher flow velocity despite low flood depths can be severely damaged due to possible scouring, bank cutting and washing away of physical infrastructure and belongings. F I G U R E 9 Water surface elevation before and after addition of embankment in extreme flood condition at five typical cross-sections (a) X 1 ; (b) X 2 ; (c) X 3 ; (d) X 4 ; and (e) X 5 . Refer to Figure 8 for the spatial locations of the cross-sections in the basin.
Our results show that settlements at high elevation (upper part of ERW) were not affected by floods as much as the lower reach settlements. The narrow gorges through which the river flows in the upstream area are generally bounded by cliffs on the sides and the settlements are usually located on highlands. As flood emerge into low plains from these high elevations, they spread out with abrupt gradient decrease resulting in altered river courses and sweeping of settlements and agricultural land in the floodplains. Because of this very reason, most of the critical areas are in lower reach of ERW. This study is confined to three aspects of socio-economic impact of flooding: impacted population, number of houses and agricultural land considering the direct economic losses. However, there are various other impacted sectors as well such as infrastructure, type of houses (public and private), roads, electricity, among others that could be areas of future research for the ERW as in the study of River Benue of Nigeria by Abubakar et al. (2020).
The effectiveness of addition of embankment on right bank of the East Rapti River as a structural adaptation measure was analyzed based upon the current data. However, population growth, rate of urbanization, change in agricultural practices, crop production rates and their economic values, among many other socioeconomic and biophysical dimensions are highly dynamic in nature which evolve continuously over time. Assumption of constant landuse for the future is a limitation of this study. Modelling all these aforementioned aspects to analyze their impacts on the flood induced damage either in silo or comprehensively could be future research arenas.
Community Based Disaster Management Project piloted by United Nations Development Program in during 2006-2008 implementing nonstructural flood adaptation measures in the Chitwan and Makwanpur districts of ERW has been considered a successful project (UNDP, 2008). More recently, ERW was one of the core areas where Flood Risk Management Project was implemented by the Asian Development Bank (ADB, 2020). The project recommended construction/improvement of flood shelters be designed and implemented through community-based consultation. Moreover, based on literature review (Devkota et al., 2016;Devkota et al., 2020;Devkota & Bhattarai, 2015;Thapa-Parajuli et al., 2018) nonstructural adaptation strategies such as communicating, community-based flood management planning taking care of affected people, exchanging help, and coordinating with government, among others are predicted to be suitable for the ERW.
Overall, we can infer that the aggregation of a datadriven models like ANN and hydrodynamic flood models is conveniently applicable for flood studies with fair accuracy in data-scarce watersheds. Similar modelling approach including analysis of output of various climate models, modelling of structural adaptation measures, and socio-economic analysis can be helpful for CC impact assessments due to inundation in other ungauged watersheds of Nepal and other regions with similar RR characteristics.

| CONCLUSION
This study carried out flood modeling in the ERW to assess the socio-economic impacts under current and future climate scenarios. Rainfall-runoff modelling was done using ANN while the hydraulic modelling was performed using HEC-RAS 1D. Further, suitable adaptation strategies were identified and evaluated to minimize the potential flood induced loss. Based on our findings, we conclude the following: • The mean monsoon temperature and rainfall of the ERW were 10 C and 1263 mm, respectively, for the baseline period. However, an increase of 0.4 C in the temperature and 98 mm in the precipitation is projected in the future under RCP 8.5 based on MPI_R-EMO data. Additionally, maximum rainfall is projected to increase by 40.5% and the magnitude of a 100-year return period flood is expected to increase to 3134 m 3 /s with flood extent to 80 km 2 and a maximum inundation depth of over 11 m corresponding to the RCP 8.5 scenario. • The extreme flood in baseline (future) period is likely to inundate a total of 712 (1173) houses with 3056 (5278) inhabitants and 29 km 2 (34 km 2 ) agricultural land, among which 9 km 2 (13 km 2 ) and 679 (1073) houses were seen to fall under critical areas. The economic loss in baseline (future) period due to inundation of agricultural land in ERW and the critical areas were estimated to be NRs 145 million (170 million) and NRs 45 million (65 million), respectively. • The use of embankments as structural adaptation measure on the right bank of critical zones of the East Rapti River could be effective in saving 3744 people living in the floodplain areas and NRs 33 million equivalent of agricultural land in the future by decreasing the flood extent by 35%. In addition, effectiveness of the structural measures can be increased by implementing community-based nonstructural adaptation measures.
This study assessed CC impacts using an ANN-RR model and two forcing variables (rainfall and temperature), neglecting other near-surface atmospheric factors, which might affect runoff generation. Evaluating impacts of using multiple variables on the results could be an area of future research. This study assesses inundation impact on limited socio-economic sectors. Incorporating other multidisciplinary dimensions may be an extension of this study. Overall findings of this study suggest that a combination of data-driven hydrological and hydrodynamic models can be a promising approach for flood impact assessments under changing climate in data-scarce watersheds of Nepal and beyond.

ACKNOWLEDGMENTS
We would like to acknowledge Mr. Bharat Bhatta for his support in the ANN work; Ms. Riyanka Pokhrel for her help in the field survey; the Department of Hydrology and Meteorology (DHM), Government of Nepal and Ms. Anupama Ray for the hydro-climatic data and Mr. Sunil Adhikari for his constant support throughout the research.

DATA AVAILABILITY STATEMENT
Data available on request from the authors.