Flood Hazard Mapping of Rivers in Snow- and Glacier-Fed Basins of Different Hydrological Regimes Using a Hydrodynamic Model under RCP Scenarios

: The global warming trends have accelerated snow and glacier melt in mountainous river basins, which has increased the probability of glacial outburst ﬂooding. Recurrent ﬂood events are a challenge for the developing economy of Pakistan in terms of damage to infrastructure and loss of lives. Flood hazard maps can be used for future ﬂood damage assessment, preparedness, and mitigation. The current study focused on the assessment and mapping of ﬂood-prone areas in small settlements of the major snow- and glacier-fed river basins situated in Hindukush–Karakoram– Himalaya (HKH) under future climate scenarios. The Hydrologic Engineering Center-River Analysis System (HEC-RAS) model was used for ﬂood simulation and mapping. The ALOS 12.5 m Digital Elevation Model (DEM) was used to extract river geometry, and the ﬂows generated in these river basins using RCP scenarios were used as the inﬂow boundary condition. Severe ﬂooding would inundate an area of ~66%, ~86%, ~37% (under mid-21st century), and an area of ~72%, ~93%, ~59% (under late 21st century RCP 8.5 scenario) in the Chitral, Hunza, and Astore river basins, respectively. There is an urgent need to develop a robust ﬂood mitigation plan for the frequent ﬂoods occurring in northern Pakistan. maps of Chitral, Hunza, and Astore basins using the HEC-RAS model under RCP climate scenarios. Flood maps may be used for the implementation of a ﬂood management program, which can signiﬁcantly reduce damages in the study area. The results of this study can potentially assist policymakers and planners to establish and implement effective ﬂood mitigation and management policies in the study area.


Introduction
Floods are a recurrent and frequent natural disaster in Pakistan [1][2][3][4]. In Pakistan, over 60% of the land area is vulnerable to floods, 40% to earthquakes, and 6% to cyclones [1]. Indus, the longest river in Pakistan, originates from the Hindukush-Karakoram-Himalaya (HKH) region and, along with its tributaries Chenab, Ravi, Sutlej, and Jhelum, runs north to south through the country [5]. Due to snowmelt and heavy monsoon rainfalls, significant flows are generated during summer [6]. As a result, the rivers overflow and enter their floodplains. The flood peaks of the rivers usually do not overlap but when they do, massive flooding results.
The floods in the Upper Indus Basin (UIB) are different from those in the Lower Indus Basin. The geology, climate of the northern areas of Pakistan, weak geological formations, intense rainfall, erosion, steepness followed by erosion, sedimentation, and meandering rivers make the region prone to natural disasters [1]. Although Pakistan was hit by severe floods in 1973, 1992, and 2006, the flood of 2010 was significantly impactful. In the possible alternative to observe and evaluate flood risk. Remote sensing data (landcover, elevation, river width and depth, flood extent, water level, etc.) used as an input in flood simulation models have a great potential to improve the understanding and forecasting of flood events [32].
Various computational models (such as HEC-RAS, HEC-HMS, MIKE FLOOD, MIKE SWMM, MIKE11, and InfoWORKS) exist to simulate velocity, depth, and other flood characteristics [22,33,34]. The Hydrologic Engineering Center River Analysis System (HEC-RAS) is among the best hydrodynamic models for flood simulation studies [2,3,21,23,[35][36][37][38][39][40][41][42][43][44][45]. HEC-RAS has been used for modeling flood inundation, management, and mitigation across the world [2,22,46,47]. The flood event of 2010 involving the Swat River was modeled [23] using the 2D application of the HEC-RAS model and reported an increase in the extent of the inundation area with increasing flood return periods. The HEC-RAS 1D-2D integrated model was reported to accurately display susceptible points in the representation of flooding in the coastal urban area of the UK for sea-level changes up to the year 2100. The study recommended considering the use of 1D-2D integrated flood modeling before decision making [48]. Similarly, ref [49] assessed the capability of the HEC-RAS model to simulate flood as a combined 1D and 2D model and compared it with Gerris and FLUMEN models. Modeling was performed for the 2011 Baeksan River flood in South Korea. It was concluded that HEC-RAS 1D-2D modeling is a robust method for the accurate simulation of flood events as compared to the other models. HEC-RAS was also applied in combination with the Watershed Modelling System and GIS to analyze the flood event in Ain Sefra city (Algeria) with or without built concrete walls for return periods of 10, 100, and 1000 years. It was concluded that retaining walls decrease the inundation extent but are not efficient for all return periods. Therefore, flood-retaining infrastructures were recommended so to ensure the safety during floods in the city in future [50]. Similarly, another study [51] used the HEC-RAS model to analyze the streamflow and develop a floodplain model for Carson River and reported an increase in the flood level 100 years into the future under climate change.
The Fifth Assessment Report (AR5) [52] of the Intergovernmental Panel on Climate Change (IPCC) established greenhouse gas concentration climate change projections recognized as Representative Concentration Pathways (RCPs). The projections' near term is fixed in 2035, while the long term is established for 2100. Four scenarios have been developed, named after the level of radiative forcing and measured in Watts per square meter (W/m 2 ) by the year 2100 [53]. These scenarios include a mitigation scenario (RCP 2.6), two mediumstabilization scenarios (RCP 4.5 and RCP 6), and a high-emission scenario (RCP 8.5). The RCP 8.5 appears to lessen greenhouse gases and to limit warming by 2100 through a little endeavor. Many researchers have used Representative Concentration Pathways (RCP) scenarios to evaluate the impact of climate change on river floods using HEC-RAS [45,54,55]. A study mapped flood hazard in the two municipalities of the Petite-Nation River basin in Quebec, Canada, under projected RCP 4.5 and RCP 8.5 climate scenarios (2050 and 2080) [56]. Flood hazard maps were generated using the GARI (Gestion et Analyse du Risque d'Inondation) tool to determine flood risk due to the changing climate. Flow hydrographs were derived for 25, 50, and 100 years from the flood event of 2007. Similarly, another work [45] employed HEC-RAS to simulate floods under future climate scenarios driven from three General Circulation Models (GCMs) in the 2020s, 2050s, and 2080s in the Yang River Basin of Thailand.
Therefore, owing to the need of predicting flood inundation extents in the rivers of snow-and glacier-fed catchments of HKH region Pakistan, this study aimed to develop flood maps of Chitral, Hunza, and Astore basins using the HEC-RAS model under RCP climate scenarios. Flood maps may be used for the implementation of a flood management program, which can significantly reduce damages in the study area. The results of this study can potentially assist policymakers and planners to establish and implement effective flood mitigation and management policies in the study area.

Materials and Methods
The objective of this study was to develop flood maps of a study area in Pakistan. HEC-RAS version 5 was used as a 2D model to map the flood area.

Study Area
The study sites were chosen within Chitral (Hindukush range), Hunza (Karakoram range), and Astore (western Himalaya) in the UIB, as shown in Figure 1. These basins were selected because each one is influenced by a slightly different climate and is situated in different geographical portions of the UIB. Small regions, mostly containing settlements, within Chitral, Hunza, and Astore basins ( Figure 1) were selected for the application of the model. The areas of these small regions within Chitral, Hunza, and Astore basins are 0.50, 0.96, and 0.83 km 2 , respectively. The basin boundaries and reach of the study areas mapped for flood are shown in Figure 1. this study can potentially assist policymakers and planners to establish and implement effective flood mitigation and management policies in the study area.

Materials and Methods
The objective of this study was to develop flood maps of a study area in Pakistan. HEC-RAS version 5 was used as a 2D model to map the flood area.

Study Area
The study sites were chosen within Chitral (Hindukush range), Hunza (Karakoram range), and Astore (western Himalaya) in the UIB, as shown in Figure 1. These basins were selected because each one is influenced by a slightly different climate and is situated in different geographical portions of the UIB. Small regions, mostly containing settlements, within Chitral, Hunza, and Astore basins ( Figure 1) were selected for the application of the model. The areas of these small regions within Chitral, Hunza, and Astore basins are 0.50, 0.96, and 0.83 km 2 , respectively. The basin boundaries and reach of the study areas mapped for flood are shown in Figure 1. Chitral basin is located at latitude 35°50′ N and longitude 71°47′ E and covers an area of ~12390.48 km 2 with a ~1988 km 2 (16%) area above 5000 m A.S.L ( Figure 2). Approximately 12.6% (~1563.025 km 2 ) of this area is covered with glaciers. This region has a subtropical Mediterranean climate, influenced by the mid-latitudes westerly and the Southwest Asian summer monsoon. Its mean elevation is about ~4586 m ASL (estimated by DEM), and the highest gauging station is located at 1500 m A.S.L ( Figure 2). The mean annual temperature of Chitral is 15.6 °C , and the total annual precipitation is 418 mm [57]. The mean annual flow of Chitral river is 291 m 3 /s ( Figure 3). Chitral basin is located at latitude 35 • 50 N and longitude 71 • 47 E and covers an area of~12390.48 km 2 with a~1988 km 2 (16%) area above 5000 m A.S.L ( Figure 2). Approximately 12.6% (~1563.025 km 2 ) of this area is covered with glaciers. This region has a subtropical Mediterranean climate, influenced by the mid-latitudes westerly and the Southwest Asian summer monsoon. Its mean elevation is about~4586 m ASL (estimated by DEM), and the highest gauging station is located at 1500 m A.S.L ( Figure 2). The mean annual temperature of Chitral is 15.6 • C, and the total annual precipitation is 418 mm [57]. The mean annual flow of Chitral river is 291 m 3 /s ( Figure 3).
The area of the Hunza basin is about 13,715 km 2 , and nearly its 25% (3433 km 2 ) is covered with glaciers ( Figure 2). Both Hunza and Astore basins are influenced by the westerlies and Indian summer monsoon [58][59][60]. In the Hunza basin, maximum precipitation is brought by westerlies in winter in the form of snow, with mean annual precipitation between 180 mm and 690 mm in the valleys. The mean annual discharge of Hunza river is 307 m 3 /s (Figure 3), while Astore is influenced by winter rainfall forced by westerlies [61], with a mean annual discharge of 135 m 3 /s. The area of the Hunza basin is about 13,715 km 2 , and nearly its 25% (3433 km 2 ) is covered with glaciers ( Figure 2). Both Hunza and Astore basins are influenced by the westerlies and Indian summer monsoon [58][59][60]. In the Hunza basin, maximum precipitation is brought by westerlies in winter in the form of snow, with mean annual precipitation between 180 mm and 690 mm in the valleys. The mean annual discharge of Hunza river is 307 m 3 /s ( Figure 3), while Astore is influenced by winter rainfall forced by westerlies [61], with a mean annual discharge of 135 m 3 /s. The Astore basin is located at latitude 35°33′ N and longitude 74°42′ E and covers an area of 3990 km ( Figure 1). The mean elevation is about 4100 m A.S.L., with maximum elevation reaching above 8000 m A.S.L ( Figure 2). Approximately 15% (607 km 2 ) of this area is covered with snow and glaciers [61]. The mean snow cover area, mean total precipitation, mean monthly discharge, and mean monthly temperature of the study area are shown in Figure 3.   The area of the Hunza basin is about 13,715 km 2 , and nearly its 25% (3433 km 2 ) is covered with glaciers ( Figure 2). Both Hunza and Astore basins are influenced by the westerlies and Indian summer monsoon [58][59][60]. In the Hunza basin, maximum precipitation is brought by westerlies in winter in the form of snow, with mean annual precipitation between 180 mm and 690 mm in the valleys. The mean annual discharge of Hunza river is 307 m 3 /s (Figure 3), while Astore is influenced by winter rainfall forced by westerlies [61], with a mean annual discharge of 135 m 3 /s. The Astore basin is located at latitude 35°33′ N and longitude 74°42′ E and covers an area of 3990 km (Figure 1). The mean elevation is about 4100 m A.S.L., with maximum elevation reaching above 8000 m A.S.L ( Figure 2). Approximately 15% (607 km 2 ) of this area is covered with snow and glaciers [61]. The mean snow cover area, mean total precipitation, mean monthly discharge, and mean monthly temperature of the study area are shown in Figure 3.   Figure 2). Approximately 15% (607 km 2 ) of this area is covered with snow and glaciers [61]. The mean snow cover area, mean total precipitation, mean monthly discharge, and mean monthly temperature of the study area are shown in Figure 3.
The soils in this region are made of the accumulated colluvial material on lower parts of mountain slopes and alluvial deposits in valleys. The foothills of the mountains are used for growing wheat and maize under rainfed cultivation, while fruit orchards are present near streams where water for irrigation is available [62]. Hindukush, Karakoram, and Himalaya are young mountain ranges and are vulnerable to erosion and disintegration by chemical, physical, and biological processes. Continuous tectonic activity and high altitude combined with seasonal effects of the climatic conditions such as heavy snowfall, freezing and thawing, and action of glaciers produce high sediment in these mountain ranges [63,64]. Sediments are carried away and deposited in the major rivers with meltwater from snow and glacial ice, by debris flow, channel erosion during flash floods, and outburst floods from landslides and glacial lakes. The lower part of the mountain ranges receives heavy Water 2021, 13, 2806 6 of 15 rainfall, and rivers in these ranges carry large amounts of sediments in flash floods. Among the tributaries of the Indus River, the Hunza river has the highest, while Chitral and Astore rivers have comparatively lower sediment flows [63,65]. A study reported the predominance of physical weathering in the release and dispersion of sediments in the Chitral basin due to the scarcity of moisture [66].

Data Sets
The Digital Elevation Model (DEM) is the most essential input data required for 2D models. Generally, the efficacy of 2D hydraulic models depends on the DEM used, because the correct topographical data influence the flood simulation results. Even though the basic input data ensures the accuracy of the model, such data are not freely available in high resolution. For accessing the topographical data required by the Digital Elevation Model, Advanced Land Observation Satellite (ALOS) 12.5 m data were used. Daily river flow data under RCP scenarios of the study area for 2000 (base year) were taken from our previous study on the simulation of river flow under the RCP scenarios published by IPCC for the HKH region [28,29]. Land use land cover (LULC) data were taken from the land cover atlas of Pakistan.

Flood Mapping Using the Hydrologic Engineering Centre-River Analysis System (HEC-RAS)
The HEC-RAS model was used for assessing and mapping the flood-prone areas within the study region under future river flows. HEC-RAS can perform 1D and 2D flood simulations using the St. Venant's equation of Conservation of Mass and Conservation of Momentum along one and two dimensions, respectively [67]. In the case of flood, water starts to overflow and becomes a 2D phenomenon; therefore, 2D modeling becomes more suitable. It is more convenient to use it in urban settings, dam breaks, and wide floodplains. A 2D HEC-RAS model lays off a mesh on the terrain data provided. Then, 2D computations and terrain representation compute where the water will move and to what extent. Therefore, unlike 1D hydraulic models, all 2D models are run in unsteady flow conditions [68].
This study used the 2D capacity of HEC-RAS for developing floodplain maps of the study area. A flowchart of the methodology used is presented in Figure 4. The Digital Elevation Model (DEM) of the study area is a topographical map that includes terrain properties. The 12.5 m ALOS DEM was used to delineate the hydrographic network in the HEC-RAS version 5.0 for each basin. The geometry of the river reaches and the flow data are important parameters for flood mapping in HEC-RAS. The next step was the representation of the terrain, using RAS Mapper to view maps and data spatially. To create the geometry of the basin, a 2D flow area tool was used to draw a polygon, near the outlet of the basin, which defined the area selected for the simulation. The first boundary condition was drawn upstream, whereas second one was drawn downstream of the river. After defining the boundary conditions upstream and downstream of the river, computational points were generated with a grid size of 20 for both horizontal and vertical spacing. While low resolution yield faster simulations, smaller grid sizes represent the terrain better, but require a longer calculation time. Normal Depth was used for the downstream boundary condition, whereas Manning's Roughness coefficients for the river reaches were assigned according to the HEC-RAS Hydraulic Reference Manual. A Flow Hydrograph was used for the upstream boundary. The Flow hydrograph window was filled using a fixed start time and introducing flow values for 2 months (June-July), as these months show the highest flow peaks. Unsteady flow was simulated using unsteady flow analysis. The results were accessed using RAS Mapper and visualized by activating velocity, WSE (Water Surface Elevation), geometry, and depth. The Google Earth image was added from Map layers of HEC-RAS to visualize the flood-prone and -affected areas around the river's mapped area.
Flow Hydrograph was used for the upstream boundary. The Flow hydrograph window was filled using a fixed start time and introducing flow values for 2 months (June-July), as these months show the highest flow peaks. Unsteady flow was simulated using unsteady flow analysis. The results were accessed using RAS Mapper and visualized by activating velocity, WSE (Water Surface Elevation), geometry, and depth. The Google Earth image was added from Map layers of HEC-RAS to visualize the flood-prone and -affected areas around the river's mapped area.

Results
HEC-RAS was used for generating floodplain maps for Chitral, Hunza, and Astore basins under RCP scenarios.

Flood Assessment and Mapping under Mid-and Late-21st Century RCP Scenarios
The flow simulated under RCP scenarios from all three basins [28,29] was used as input in flow conditions of the HEC-RAS model. A summary of the results is presented in Table 1.

Results
HEC-RAS was used for generating floodplain maps for Chitral, Hunza, and Astore basins under RCP scenarios.

Flood Assessment and Mapping under Mid-and Late-21st Century RCP Scenarios
The flow simulated under RCP scenarios from all three basins [28,29] was used as input in flow conditions of the HEC-RAS model. A summary of the results is presented in Table 1. The RCP mid-and late-21st century scenarios were applied to the simulated flow to compare it with the flow of the base year and create an inundation map of flood-prone areas. The results showed some overflow, when compared to the base year. The yellow color in the flood maps represents the areas affected due to flood, with depth ranging from 0 to 15 ft; light blue represents the minimum depth, and dark blue the maximum depth.

Flood Assessment and Mapping in Chitral Basin
The maximum runoff in the Chitral river was observed on 14 July and was 1123.30, 1154.60, and 1194.70 m 3 /s for RCP mid-(2.5, 4.5, and 8.5) and 1181.80, 1180.60, and 1309.50 m 3 /s for RCP late-(2.5, 4.5, and 8.5) 21st century scenarios. The underlying area of the Chitral river simulation is shown in Figure 5a. The results of RCP late 21st century showed greater flow depth and extent as compared to the flow simulated under mid-21st century scenarios ( Figure 5).
The RCP mid-and late-21st century scenarios were applied to the simulated flow to compare it with the flow of the base year and create an inundation map of flood-prone areas. The results showed some overflow, when compared to the base year. The yellow color in the flood maps represents the areas affected due to flood, with depth ranging from 0 to 15 ft; light blue represents the minimum depth, and dark blue the maximum depth.

Flood Assessment and Mapping in Chitral Basin
The maximum runoff in the Chitral river was observed on 14 July and was 1123.30, 1154.60, and 1194.70 m 3 /s for RCP mid-(2.5, 4.5, and 8.5) and 1181.80, 1180.60, and 1309.50 m 3 /s for RCP late-(2.5, 4.5, and 8.5) 21st century scenarios. The underlying area of the Chitral river simulation is shown in Figure 5a. The results of RCP late 21st century showed greater flow depth and extent as compared to the flow simulated under mid-21st century scenarios ( Figure 5).  The late RCP 8.5 generated a flood covering an area of 0.36 km 2 (72%), which was largest compared to the inundated area under all other mid-and late-21st century scenarios; almost all nearby populated areas appeared vulnerable to it, as shown in Figure 5g. An approximate area of 0.04 (8%), 0.30 (60%), and 0.33 km 2 (66%) was flooded under RCP mid-21st century scenarios (2.5, 4.5, and 8.5), while an area of 0.05 (10%), 0.32 (64%), and 0.36 km 2 (72%) was flooded under RCP late-(2.5, 4.5 and 8.5) 21st century scenarios.

Flood Assessment and Mapping in Hunza Basin
In Hunza river, the maximum runoff was observed on 3 August and was 1606.30, 1681.3, 1771.80 m 3 /s for RCP mid-(2.5, 4.5, and 8.5) and 1585.60, 1745.90, 2058.20 m 3 /s for RCP late-(2.5, 4.5, and 8.5) 21st century scenarios. The affected area of the Hunza river simulation is shown in Figure 6a. The late RCP 8.5 generated a flood covering an area of 3.67 km 2 (93%), which was the largest compared to inundated areas under all other midand late-21st century scenarios ( Figure 6).

Flood Assessment and Mapping in Astore Basin
Maximum runoff was observed on 7 August in the Astore river and was 513.20, 538.8, and 1099.50, m 3 /s for RCP mid-(2.5, 4.5, and 8.5) and 509.30, 558.80, and 1288.20 m 3 /s for RCP late-(2.5, 4.5, and 8.5) 21st century scenarios. The affected area of the Astore river simulation is shown in Figure 7a. The highest flow was recorded under the late 8.5 scenario with respect to all other scenarios (Figure 7).  With the increase in temperature and precipitation in the different scenarios, the depth of the flood also increased. An area of 0.49 km 2 (59%) appeared to be vulnerable to flood for late 8.5, as compared to other scenarios. In Astore, an approximate area of 0.05 (6%), 0.12 (14%), and 0.31 km 2 (37%) was flooded under RCP mid-21st century scenarios (2.5, 4.5, and 8.5), and area of approximately 0.12 (14%), 0.18 (21%), and 0.49 km 2 (59%) were flooded under RCP late-(2.5, 4.5, and 8.5) 21st century scenarios.

Discussion
Based on the results obtained, it was observed that an increase in river discharge under mid-and late-21st century RCPs would result from an increase in temperature and precipitation. Increasing precipitation augments snow and glacial ice, while increasing temperature increases snow and glacial melt, resulting in an increase in river discharge [45,69]. The flood generated by the RCP 2.6 was of less depth and extent as compared to that generated by RCP 4.5, which was less than that of RCP 8.5 for both mid as well as late 21st century. More structures (houses, hotels, institutions, and fields) appeared vulnerable to flood in the case of RCP 8.5, as compared to RCP 4.5, which involved more vulnerable structures than RCP 2.6. These findings can help policymakers to plan and prevent damages to the society and economy of Pakistan through effective natural disaster management practices. The Intergovernmental Panel on Climate Change (IPCC)'s report of 2013, reported a 0.85 • C increase in the global temperature from 1880 and projected an increase of 3.7 • C at the end of the 21st century [52]. A study [70] investigated the impact of climate change on precipitation and flood using 24 daily precipitation simulations from the Coupled Model Intercomparison Project Phase 5 (CMIP533) GCMs regarding the period of 1971-2000 and the future period of 2070-2099 under the Representative Concentration Pathway (RCP) 8.5 scenario. The results showed an increase in precipitation and frequency of floods in all climatic regions. RCP climate scenarios presented by IPCC for the Indus River Basin (IRB) are useful to assess the river flow [51,71,72].
In a similar study [73], the authors performed flood hydraulic modeling and floodplain mapping for the Chashma-Taunsa (252 km) reach of the Indus basin. HEC-RAS, ArcGIS, and its extension Hec-GeoRas were used for hydraulic modeling and mapping of floodplains. Flood maps show that the depth of the flooded area has increased from 7.1 m to 9.1 m over 25 years, with 500 floods over the return period. A study used HEC-RAS to develop floodplain maps for a certain reach of the Kabul River in Pakistan. SRTM DEM was used for preparing geometrical data of the river [2]. The highest flows from frequency analysis were used as input data into HEC-RAS to estimate the corresponding flood levels expected along river reaches. The interpretation of floodplain maps suggested that more than 400% of the area, predominantly agricultural, is susceptible floods, contrary to normal river discharge.
The China-Pakistan Economic Corridor (CPEC), a joint project of China and Pakistan, passes through the Hunza basin in Pakistan. Hunza has many surging glaciers which create glacial lakes [74]. Such glacial lakes are used for freshwater but are extremely dangerous when their dams are breached, resulting in glacier lake outburst floods and flash floods. Since the route passes through a flood-and other hazard-prone area, GLOFs and flash floods should be considered in the CPEC policy formulation. This would be helpful, as the Government may take concrete steps in the mitigation and prevention of floods. Another study stressed the usefulness of inundation maps for flood preparedness, immediate response to future floods, assessment of damage, mitigation, and flood risk analysis [75]. Tools like HEC-RAS along with experience and skills can provide important guidelines for policymakers, helping them to elaborate plans to respond, mitigate, and manage future floods.

Conclusions
Flood hazard mapping of snow and glacier melt mountainous basins is important for flood mitigation and preparedness to floods as it allows identifying the vulnerability of a place and its community. The HEC-RAS model was used in order to develop a flood map of the study area and assess flood risk under mid-and late-21st century RCP scenarios. The results of this study indicated that:

•
Floods of mid-as well as late-21st century RCP 8.5 are predicted to have greatest depth and velocity, so will lead to more devastating inundations as compared to floods of mid-and late-21st century RCP 2.6 and 4.5.

•
Under the RCP 8.5 late-21st century scenario, a peak flood was simulated covering an area of~72,~93, and~59% of the studied reaches within Chitral, Hunza, and Astore river basins, respectively. • Flood simulation and GIS can be used effectively to develop flood management strategies. The digitization of floods under RCP scenarios can help develop plans to minimize and prevent damage from annual flooding.
High-resolution satellite images incorporated with terrain information should be used effectively for continuous monitoring of flood-prone areas to mitigate the risk of floods in the future. Given the frequent nature of floods in Pakistan, there is a pressing demand to develop a robust flood mitigation plan. Other than structural flood mitigation measures, nonstructural measures such as awareness campaigns, natural disaster management, early warning systems, and flood forecasting are likely to be very productive for the mitigation of the direct and indirect impacts of floods.