Hydrodynamic Modelling of Floods and Estimating Socio-economic Impacts of Floods in Ugandan River Malaba Sub-catchment

River Malaba sub-catchment tends to experience dramatic flooding events, with several socio-economic impacts to the nearby communities, such as loss of lives and destructions of physical infrastructure. Analysis of spatiotemporal extents to which settlements, crops and physical infrastructures tend to be inundated are vital for predictive planning of risk-based adaptation measures. This paper presents a case study on flood risk assessment for Ugandan River Malaba sub-catchment. We applied the two-dimensional Hydraulic Engineering Center’s River Analysis System (2D HEC-RAS) for modelling of flooding extents. We considered extreme flow quantiles, lower and upper quantiles corresponding to the 95% confidence interval limits aimed at determining uncertainties in the flooding extents. Spatial extents of inundation on human settlement, land cover and infrastructure were analysed with respect to return periods of extreme flow quantiles. Finally, we estimated economic loss on infrastructure due to flooding. Results from the 2D HEC-RAS model were satisfactorily comparable with the results of observations. Amongst the land use types, cropland exhibited the highest vulnerability with at least 10,234.8 hectare (ha) susceptible to flooding event of 100-year return period (YRP). Inundated built-up land-use exhibited the highest vulnerability percentage increase (90%) between 2- and 100-YRP. In US Dollar, about US$ 33 million and US$ 39 million losses are estimated at 2- and 100-YRP, respectively, due to inundated rice gardens and these indicate a looming high risk of household food insecurity and poverty. Several infrastructure including 15 academic institutions, 12 health facilities, 32 worshiping places remain annually vulnerable to flooding. At least 6 km and 7 km of road network are also susceptible to flooding under extreme flows of return periods 2 and 100 years, respectively. Churches exhibited the highest economic losses of US$ 855,065 and US$ 1,623,832 at 2-YRP and 100-YRP, respectively. This study findings are relevant for planning the development of sustainable flood risk adaptation pathways given the established destructions within the sub-catchment due to flooding.


Introduction
Floods are natural disasters with the highest global occurrence and broader physical distribution producing severe damages in many countries (Najibi and Devineni 2018). The occurrence of floods could be attributed to human factors. For instance, conversation of land use types from one state to another for various purposes, such as settlement, agriculture, urbanisation, and other recreational activities especially near the river banks (Alawamy et al. 2020). This can result in extensive unexpected implications, such as alteration of the hydrological processes, hence exacerbating the occurrence of weather disasters (Apollonio et al. 2016;Kilama Luwa et al. 2020;Szwagrzyk et al. 2018;Zope et al. 2016). In addition, climate variability could be regarded as one key cause of the disasters. The most updated climate projections indicate that floods in several parts of the world will likely increase in intensity and frequency in the future (Hirabayashi et al. 2021).
Several flooding events (Table 6 of Appendix A) have impacted Uganda. These events result in considerable socioeconomic impacts including damage of crops, destruction of infrastructures including roads, bridges, schools, and loss of livestock and human lives. River Malaba sub-catchment is an example of the floodplains which have suffered dramatic floods, especially in its low-lying areas (Mayega et al. 2015;Ministry of Water and Environment 2018). The nearly flat topography facilitates the emplacement of houses or infrastructures. However, human lives and their property are at high risk of inundation. In River Malaba sub-catchment and the neighbouring areas, rainfall shows an increasing trend and multi-decadal variability (Mubialiwo et al. 2020;Onyutha et al. 2020). Furthermore, the study by Onyutha et al. (2021a) reported a likely general increase in precipitation in the 2050s, 2070s and 2090s over the entire Uganda, where the study area is located. The positive trends in rainfall have resulted in more recurrent floods over the period 2010-2019 than from 1960 to 2009 . However, the positive trends in rainfall may be affected by rising temperatures due to likely increase in evapotranspiration (Xu et al. 2017). Besides, flooding events may increase at a slower rate than rainfall which is an indication of climate change impact (Bucchignani et al. 2018;Padulano et al. 2021). The sub-catchment hydrological regimes have also been greatly influenced by land use changes (Jiang et al. 2014). The possibility of experiencing more floods of varying magnitudes in future in the study area presents a worrying situation for water resources management. To effectively manage the associated impacts of flooding, there is need for predictive planning and operation of climate change adaptation measures (Xu et al. 2017). Such planning can be supported by provision of comprehensive information or maps of flood inundation and extents of infrastructure and property that can be affected by flooding as well as vulnerability of the local population. In addition, to reduce the effects of natural disasters, there is need to evaluate their socio-economic impacts and to consider management strategies (structural and non-structural).
Flood inundation mapping can be achieved by mimicking the physics of floodwaters in the river systems through hydrodynamic modelling using a number of numerical tools. Numerical (hydrodynamic) models can be categorised as one-dimensional (1D), two-dimensional (2D) or threedimensional (3D) (Teng et al. 2017). 1D models consider one directional flow along the centreline of the river channel and are suitable when detailed analysis is not essential (Alkema 2007). 2D models make it possible to evaluate the landscape effect on the magnitude of flooding especially in near-flat topography (Alkema 2007), such as the present study area. 3D models are rarely used, since 2D shallow water approximations are considered adequate especially in areas with limited high resolution data, such as the current study area (Teng et al. 2017). 2D models have exhibited commendable performance. For instance, the study by Dhungel et al. (2019) revealed superior performance of 2D model in simulating flow depth and velocity in the lower Provo river, USA. Ghimire (2019) established that with fine resolution survey data and availability of local data, 2D models can present sound performance.
There are various 2D hydrodynamic models, such as the Hydraulic Engineering Center's River Analysis System (2D HEC-RAS) (Brunner 2016), SOBEK Suite DELFT3D (Deltares systems 2019), MIKE 21 (DHI 2017), TUFLOW HPC (Heavily Parallelised Compute) (BMT-WBM 2018). Whereas most hydrodynamic models are commercial, 2D HEC-RAS is freely available for use. HEC-RAS is a popular modelling program that computes water surface profiles in natural rivers and other channels applicable in floodplain management.
Recently, 2D HEC-RAS model has attracted extensive applications in flood risk assessment. For instance, Pinos and Timbe (2019) compared the performance of four twodimensional hydraulic models 2D HEC-RAS included, in regard to estimation of flood inundation maps. The study revealed reliable performance of 2D HEC-RAS (Pinos and Timbe 2019). Ongdas et al. (2020) presented commendable performance of 2D HEC-RAS when applied for flood hazard maps generation for Yesil (Ishim) river in Kazakhstan. Other studies that applied 2D HEC-RAS include (Garcia et al. 2020;Kumar et al. 2019;Rangari et al. 2019).
2D models have been identified to provide suitable information to policymakers as they can clearly indicate which areas are to be inundated faster to make proper evacuation and/or resettlement plans (Alkema 2007). To the best of the authors' knowledge, by the time of conducting the present Published in partnership with CECCR at King Abdulaziz University study, there were no previous studies that performed hydrodynamic modelling of floods in the study area which could contribute to establishment of sustainable risk-based water resources management strategies. Therefore, the present study applied the 2D HEC-RAS hydrodynamic model for the purpose of identifying the inundation extents across the study area during flooding events of varying return periods. The study further assessed the socio-economic impact of floods on the society.

The Study Area
River Malaba sub-catchment ( Fig. 1) is transboundary with a total drainage area of about 3500 km 2 , shared between Uganda (approximately 69% or 2395 km 2 ) and Kenya (around 31% or 1100 km 2 ). The sub-catchment stretches between latitudes 0º 19′ N and 1º 07′ N and longitudes 33º 37′ E and 34º 37′ E. It comprises of River Malaba formed by two tributaries of Lwakhakha and Makalisi and later joined by Lumbaka/Kibimba, Osia and Mahanga tributaries. The sub-catchment is part of the Mpologoma catchment within Kyoga Water Management Zone (KWMZ). River Malaba originates from Mount Elgon at about 4315 m above sea level criss-crossing through districts of Bududa, Manafwa, Tororo, Butaleja reaching the shores of Lake Kyoga on the western side at about 1055 m above sea level, where it discharges into Mpologoma River.
Because of the steep topographic gradient between the highlands of Mount Elgon and the lowlands, the upstream areas of Manafwa and Bududa experience landslides, while the downstream regions suffer from floods (Mayega et al. 2015;Ministry of Water and Environment 2018). The subcatchment main land use is rain-fed subsistence agriculture engaging almost 85% of the population. While dense forests are found in the highlands of Mount Elgon, the remaining areas comprise of agricultural and grassland, fallow land, and isolated woodlots. Petric Plinthosols and Gleysols form the major soils types in the sub-catchment but with other Published in partnership with CECCR at King Abdulaziz University categories including Lixic ferrasols, Acric ferrasls and Nitisols (Barasa et al. 2017;Kitutu et al. 2009).
At the national level, Uganda's climate is mostly influenced by the Monsoons and Inter-Tropical Convergence Zone (ITCZ). This results in a bimodal rainfall distribution as they move from the North to South and viceversa (Ministry of Water and Environment 2014). However, locally, the climate of the present study area is somewhat influenced by the presence of large water bodies (such as Lake Victoria and Kyoga) and the Mount Elgon slope breezes with a tendency of disturbing the afternoon convection (Camberlin 2009). The average annual rainfall in the area reaches 1375 mm, but Bududa and Manafwa districts experience slightly much more rains (of up to 1800 mm). The March to May and June to October happen to be the wettest and driest seasons, respectively, with an average temperature range of 15.8-30.6 °C (Barasa et al. 2013).

Research Data
In this study, both primary and secondary data sets were considered as detailed below. Primary data included the river geometry and socio-economic information. Secondary data comprised hydrometric, physiographic, physical infrastructure. Other secondary information included the household size which was adopted from the previous study (Mubialiwo et al. 2021a).

Digital Elevation Model and River Geometry
The digital terrain model (DEM) of the study area with a resolution of 12.5 m was obtained from the ALOS PAL-SAR database, managed by the Alaska Satellite Facility (via https:// asf. alaska. edu/ data-sets/ sar-data-sets/ alos-palsar/; Alaska Satellite Facility 2021) (accessed: 10 June, 2021). The river geometry including cross sections with an interval of 500 m were obtained using the Engineer's Geographic Positioning System Real-Time Kinematic (GPS-RTK). The floodplain extents of the October 2007 and March 2010 events at various locations were obtained by established the flood marks (mainly flood depth) following well-documented approaches (Koenig et al. 2016). This information was used to validate the flood inundation polygons. Furthermore, information about structures along the river such as bridges within the floodplains was obtained. In this study (especially in low-lying flood prone areas), most of the bridges were simple structures. Based on this, culverts were applied as an estimate, representing the simple bridge openings.

Land Use Map
Land use land cover (LULC) maps with a resolution of 30 m for three different years ( , 2015 were obtained from the National Forest Authority (NFA) of Uganda. The obtained maps were already classified by NFA by unsupervised classification. This approach has recently been recommended by Ma et al. (2020) due to the lower field data sample requirements as compared to the supervised method. In the present study, LULC changes over the study area from 2010 to 2017 were quantified prior to hydrodynamic modelling (Fig. 3). Bare land, built-up areas and cropland increased from 2010 to 2017, while water surface decreased (Fig. 3). A decrease in water surface is an indication of reduced available water. Despite a decrease in forestland from 2010 and 2015, a slight increase occurred from 2015 to 2017. Forestland exhibited the highest decrease from 2010 to 2015, while crop land had the highest increase from 2010 to 2017. Grassland increased from 2010 to 2015 and later decreased.

Flow Data and Processing
Simulated daily rainfall-runoff was adopted from the earlier study . The series covered the period 1999-2016 totalling to 18 years. Daily rainfall-runoff was simulated using seven lumped conceptual rainfall-runoff models (Table 7 of Appendix B). Evaluation of the mismatches between observed and simulated flows were based on a total of nine "goodness-of-fit" metrics ). Rainfall-runoff discharge was adopted instead of observed river flow to capture the dynamics of modifying land use types and their corresponding effect on hydrology. Besides, observed river flow cannot reflect the amount of flow in the floodplain as it is only measured at one location. It is vital to differentiate between rainfall-runoff discharge and river flow. Rainfall-runoff discharge refers to the inflowing discharge into the river, and was considered (in a lumped way) upstream of a given location along the river; or along a given river stretch. In other words, rainfall over the catchment generates runoff which flows into the river. The amount of runoff depends on the intensity of the rainfall. In turn, the magnitude of flooding event depends on the amount of runoff generated. The river flow (at a given point) comprises the transformation of the rainfall-runoff discharges (upstream of that point, and at previous time steps), after being routed through the river network.
To guide further use of the rainfall-runoff discharge data in hydrodynamic modelling, trend statistics were computed to assess the presence of either monotonic trends or seasonal components (Table 1). Analyses adopted the same approaches previously applied by Mubialiwo et al. (2020Mubialiwo et al. ( , 2021c, and Onyutha et al. (2021b) as follows: The linear 1 3 Published in partnership with CECCR at King Abdulaziz University trend analyses on rainfall-runoff discharge followed the approach by Sen (1968) and Theil (1950), while significance of the trend was done using the cumulative sum of the difference between exceedance and non-exceedance counts of data (Onyutha 2021). Details of the adopted approaches for trend analyses are provided in Appendix C.
Rainfall-runoff discharge exhibited positive trends in both seasonal and annual time scales.
Except for MAM, the JJAS, OND and annual rainfall-runoff discharges had H o (no trend) rejected (Z > 1.96, p < 0.05) . This could be attributed to the increasing precipitation amount and decreasing potential evapotranspiration. Significant positive sub-trends occurred from around 2007. An important step in the hydrodynamic modelling is selection of return periods for hydrographs to drive the model. From the fitted extreme value distribution (EVD) , the extreme flow quantiles (return levels) at six different return periods (2, 5, 10, 25, 50, and 100 years) were extracted (Table 2). Apart from using the extreme flow quantiles (Q), the lower (LQ) and upper (UQ) quantiles corresponding to the 95% confidence interval limits were computed to determine uncertainties in the flooding extents (Table 2) (Ialongo 2019). Synthetic hydrographs were developed for each return period and used as

Physical Infrastructures
Distribution of settlement clusters, and location of physical infrastructures (such as churches, schools, health facilities) were obtained from the Uganda Bureau of Statistics (Uganda Bureau of Statistics 2020, 2018) ( Fig. 4).
Ground measurements were also taken using the GPS-RTK to confirm the location of some infrastructure especially the airport, train station and some roads.

Description of the Used Model
As stated before, the Hydrologic Engineering Center's River Analysis System (HEC-RAS) was applied in this study. The model can simulate both steady and unsteady flow components. The steady flow element has the competence of modelling the subcritical, supercritical as well as mixed flow regime, water surface profiles. On the other hand, the unsteady flow component of HEC-RAS is capable of modelling independent 1D or 2D unsteady flow or combined 1D and 2D unsteady flows (Brunner 2016). Herein, 2D HEC-RAS model 5.0.7 was used to simulate floods. The study involved computation of flood depth, and assessment of inundation extent on population, crops and infrastructure. The model simulated floods according to the following expressions (Brunner 2016): where h is the water depth (m), p and q denote the specific flow in the directions of x and y measured in m 2 s −1 . r is the net rainfall in m . H is the surface elevation in m, while g stands for acceleration due to gravity measured in ms −2 . n is the manning's coefficient sm −1∕3 , is the density in kgm −3 . Published in partnership with CECCR at King Abdulaziz University xx , xy , yy represent the effective shear stress components, while k is the Coriolis measured in s −1 . For the case when the analysis considered diffusive wave equation, then the inertial components in Eqs. (2) and (3) is left out.
HEC-RAS has the capability to analyse either a fully 2D Saint Venant equation or the 2D diffusive wave equation following an implicit finite-volume approach (Brunner 2016). Both 2D diffusive wave and Saint Venant equations can exhibit close performance (Dhungel et al. 2019;Ongdas et al. 2020;Shustikova et al. 2019). However, in the present study, the 2D diffusive wave equation was adopted, because it supports longer time steps, is faster, and produces stable and accurate solutions than the Saint Venant equations (Ongdas et al. 2020;Sarchani et al. 2020).

Uncertainty Analysis and Calibration of Manning's Roughness Coefficients
In this study, we assessed the possible uncertainties induced in the flow/flood simulation by changes in land use from 2010 to 2017. This was done by evaluating the difference between the three LULC maps with respect to Manning's roughness coefficients. First, calibration was performed on the six major land use types (Fig. 2) by applying the Manning override regions function in HEC-RAS (Brunner 2016). Subsequently, the calibrated parameters were then used to simulate the floods in HEC-RAS. Because of the small sample size n of 6, the Student's t− distribution was adopted (Brereton 2015). with where s p is an estimator of the collective standard deviation of the two samples, s 1 and s 2 denotes the standard deviation, x 1 and x 2 are the mean values, n 1 and n 2 are the sample size of the manning's roughness coefficients from say 2010 and 2015, respectively. The p value can be computed using Ms Excel function TDIST (value, df , tail). Considering the period of flow data (1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016), the 2015 LULC map (Fig. 2b) was adopted in defining the manning's coefficient in the HEC-RAS model.

Hydrodynamic Modelling of Floods
Modelling was done using the inputs as DEM, 2015 land cover map (for manning's coefficients), and synthetic hydrographs for different return periods (see Sects. 2.2.1 to 2.2.3 for details). The 2D flow area covered the low-lying region which is prone to flooding (Fig. 5). Two different boundary conditions (hydrograph and normal depth) were defined. Composite hydrograph boundary conditions were defined upstream of every reach representing the runoff discharges. The normal depth boundary conditions were defined at the downstream of the river and at the edges of the model, where the runoff is presumed to flow. 2D HEC-RAS calculated the outflow, flood depth, water surface elevation, and velocity for each grid as the flood-water propagated throughout the domain to the normal boundary conditions. The 2D computational mesh was set up with a cell size of 12.5 × 12.5 m. This selection of grid was to ensure minimal or no deviation from the DEM resolution (12.5 × 12.5 m). Selection of particular time steps in HEC-RAS may influence the results. Small time steps can result in very long computation times, while too large time steps may as well cause numerical diffusion hence model instability (Brunner 2016). In the present study, stability of the model was ensured by estimating the time step based on the Courant-Friedrichs-Lewy approach as follows (Brunner 2016): where C stands for the Courant number (dimensionless), V is the flood wave velocity (ms −1 ), Δt is the computational time step (s),Δx denotes the grid size (12.5 m).
Evaluation of model results is an essential component to measure performance. Evaluation of the 2D HEC-RAS model necessitated historical record of flood events. This could be in form of delineated extent from observations or remote sensing images. Unfortunately, it is worth noting that there are no official historical records of observed flood events extents for the current study area despite the recurrent disasters. Even the available satellite images do not clearly capture the entire present study area. Nevertheless, the flood image for the 20 October 2007 flood by Red cross was used (accessed: 11 July, 2021). This event affected almost the entire Uganda (Reliefweb 2007). Prior to its use in model evaluation, the flood image (map) was corrected with guidance of the flood marks information obtained as per Sect. 2.2.1. The image was also compared with another one by Reliefweb (2020) (accessed: 11 July, 2021). Model performance was assessed based on several metrics including False Alarm Ratio ( FAR ), Hit Rate ( HR ), Critical Success Index ( CSI ) and Bias Error ( B e ) as shown in Eqn. (7) to (10) (Sampson et al. 2015). These performance criteria have been applied in recent studies (Costabile et al. 2020; Ongdas Published in partnership with CECCR at King Abdulaziz University et al. 2020; Sampson et al. 2015). Consider A sim ,A obs as the simulated and observed flooded areas, respectively. Assume A cor , A over , A under denote the correctly, over-and under-predicted flood areas, respectively, all measured in km 2 : FAR is a measure of model overprediction and it varies between 0 and 1. A value of 1 implies "all false alarm", while 0 denotes "no false alarm". HR shows how well a model mimics the benchmark inundation without necessarily chastising for overprediction. HR ranges between 0 and 1, with 1 indicating that all wet areas in the benchmark inundation map are as well wet in the simulated output. CSI gives the combined effect of FAR and HR catering for both overprediction and underprediction. A value of 1 means perfect match between simulated and benchmark, while 0 indicates no match between the two. Finally, B e has two ranges (0-1 and 1-∞ ). The range of 0-1 indicates underprediction of the model, and 1-∞ shows propensity of a model to overpredict (Sampson et al. 2015).

Flood Damage Estimation
Determination of economic losses due to flooding necessitated the flood depth and flood depth-damage functions for each infrastructure. The flood depths for individual infrastructure (churches, mosques, academic institutions and health facilities) spatially distributed within the inundated area were obtained from the generated flood depth map for each return period. We estimated the flood economic Published in partnership with CECCR at King Abdulaziz University damage values for the extreme flow (Q) conditions. There are no site-specific damage functions for the present study area neither for the neighbouring areas. Therefore, depthdamage functions from other areas of similar topographical representation were applied to estimate the economic damage value for each of the selected infrastructure. The depth-damage function were adopted from (Chen 2007;Pinos et al. 2020;Scorzini and Frank 2017). These were applied in conjunction with the flood damage functions and damage values provided by Huizinga et al. (2017) at national level. The depth-functions for a particular infrastructure at each depth location (e.g., church 1, church 2) were obtained by interpolating from the known depth-functions. The applied depth-damage functions for each infrastructure at different return period can be found in Table 8

Uncertainty Analysis and Calibration of Manning's Roughness Coefficients
The land use types and their calibrated manning's roughness coefficients are shown in Table 3. Calibration was intended to improve the HEC-RAS model performance.
The p values of 0.061, 0.064 and 0.379 for the combinations of 2010 and 2015, 2010 and 2017, and 2015 and 2017, respectively, were greater than (0.05) implying that H o was not rejected. This implies that changes in LULC from 2010 to 2017 was insignificant. Therefore, the river flow variation and or damage from floods within River Malaba sub-catchment could be attributed to other factors beyond LULC changes, such as climate variability. Previous study by Onyutha et al. (2021d) revealed that changes in LULC contributed only 8% to the River Mpanga flow variation. This contribution was smaller compared to other factors, such as climate variability (70%) (Onyutha et al. 2021d). Furthermore, uncertainties in flood inundation simulations may result from the adopted manning's roughens coefficient values. In most cases, uniform instead of spatially varying values of roughness coefficient are assumed which could be a source of uncertainties (Huang and Qin 2014). However, characteristic of flooding (extent and depth) may possess a low sensitivity to changes in manning's roughness coefficient (Höffken et al. 2020).

Model Performance in Mimicking the Historical Flood Event
The performance of the HEC-RAS model in capturing the historical flood image are presented in Table 4. The value of CSI indicates that the simulated flooded areas generally concur with the observed data. The value of HR above 0.5 and close to 1 further confirms an acceptable fit of wet areas in both simulated and observed data set. However, values of B e above 1 indicate an overestimation of the model.

Spatial Inundation Extent Analyses
Vulnerability analyses were done to identify the extent of flooding on settlement, land cover, roads and other infrastructures. Figure 6 shows the inundation extent of settlements considering the LQ, Q and UQ for varying return periods. Flooding was highly noticeable in the low-lying areas. The spatial comparison established that inundation extent of settlements in the study area varies significantly. Settlement clusters are much concentrated in the lowest part of the flood-prone area. This could be attributed to the societal stalemates such as less onerous way of setting up buildings with always profitable investments. This statement conforms with the findings by Kiyengo et al. (2019). The study focused on establishing the spatio-temporal distribution of flash floods and evaluate why people preferred  Published in partnership with CECCR at King Abdulaziz University settling in flood-prone area in Lubiji micro-catchment, a Kampala suburb (Kiyengo et al. 2019). The escalating settlements in flood-prone areas by low-income communities worsens their vulnerability and habits. The poor land use management further contributes to increasing floods in such areas (Ramiaramanana and Teller 2021). It is noticeable that the difference between results of 50and 100-year return periods in terms of the spatial coverage of flooded area was small even though it was expected to be large. This can be explained in three ways. First, the difference between quantiles of two return periods depends on the scale parameter (or the slope of the weighted linear regression line in quantile-quantile plot) especially for normal tailed EVD (Beirlant et al. 1996;Onyutha and Willems 2015;Willems et al. 2007). A large-scale parameter indicates high extreme flow variations (Onyutha and Willems 2015). The second explanation of the small difference between 50-and 100-year quantiles is due to the uncertainty in the statistical estimation of quantiles. In extreme value analysis, uncertainty in estimated flow quantiles increases as the considered return period goes higher than the length of the observed data series (Onyutha and Willems 2013). In this line, extrapolation of quantiles is recommended to be limited for return periods not exceeding two or three times the length of the observed data (Kangieser and Blackadar 1994). If, say, a series longer than 50 years was used, the 100-year return level could have probably been larger than the one used in this study. Third, the topography/bathymetry which can influence LULC of the study area cannot change very much in these few years, unless the region as suffered a natural disaster. This justified by the insignificant LULC changes from 2000 to 2014 over two catchments within Uganda (Onyutha et al. 2021c, d). Thus, although the 100year return period was considered for relevance, corresponding results should be cautiously interpreted given the uncertainty involved in the estimation of the extreme quantiles.
Field activities revealed that on average, a settlement cluster had 10 households, with an average household size of 9.7 people in the present study area (Mubialiwo et al. 2021a). These values were adopted in estimating the population exposed to floods (Fig. 7). It should be noted that this number excludes those affected by landslides in the highland areas of Bududa and Manafwa districts. At 2-year return period, at least 7760 people are exposed to floods considering the lower limit of extreme flow at 95% confidence interval. The exposed population would increase by 12% and 16% for extreme and upper limit of extreme flow at 95% confidence interval, respectively, considering the 2-year return period (Fig. 7). To avert future loss of lives, authorities could make arrangement for relocation of population from risky to safer places.
There was an increase of about 13.0%, 13.3% and 15.9% for LQ, Q and UQ, respectively, of the exposed population from 2-to 100-year return period (Fig. 7). The government through ministry of relief, disaster preparedness and refugees has always resettled people from high risk areas in the gazetted refuge settlements (Reliefweb 2019). However, the population relocated to safe places usually returns to flood-prone areas. The key underlying reasons could include (1) population pressure in refugee camps, (2) inappropriate information on disaster preparedness, (3) cultural beliefs influencing people's ability to cope with the new environment, (4) poor service delivery in the camps, and (5) infertile land in refugee settlements (Osuret et al. 2016). Successful resettlement necessitates sensitisation of the community about its rationale as a way of building trust in the government (Neema et al. 2018). Figure 8 shows the variation of flooding depth across the study area simulated for different return periods (2-100 years). The maximum depth varied from 3.79 to 4.53 m at LQ for 2-year return period and 4.62 to 5.36 m at UQ for 100-year return period. As earlier observed with  The flooding extent on various land cover types is shown in Fig. 9 for the extreme flows. Results on inundation of land cover types under LQ and UQ can be found in Fig. 13 of Appendix E. Flood vulnerability analyses results conform with the October to November 2020 field survey (Mubialiwo et al. 2021a) indicating highest percentage of vulnerability in the low-lying areas. The highest vulnerability percentage occurred in the cropland (cultivated land) with 10,234.8 hectare (ha) flooded by 100-YRP, followed by 3793.9 ha of bare land. Inundation of cropland increased by 15% from 8685 ha at 2 years to 10,235 ha at 100 years. Despite the relatively small areas of inundated built-up, these places have high population density of around 9.7 people per household (Mubialiwo et al. 2021a) mainly staying along the river banks. Besides, the inundated built-up land-use type exhibited the highest vulnerability percentage increase (90%) between 2-and 100-YRP followed by grassland (80%) (Fig. 9). Therefore, in the low-lying areas, effective land-use conservation and management should be emphasised.
For the inundated cropland, field survey analyses done in October to November 2020 revealed that rice is the most cultivated crop taking approximately 40%, followed by maize (20%), millet/sorghum (15%). The other crops grown include beans, ground nuts, sweet potatoes, cassava sharing the remaining percentage. In the present study, rice which mainly doubles as both food and cash crop has been taken as an example to quantity the would-be loss if the land is inundated. Focus group discussions conducted with Namunasa stream rice farmers cooperative society limited revealed that on a good season (without inundation), one hectare of land can on average, yield 2.74 tons. This value is not much different from the country-wide rice yields of say 2.84 ton/ ha in 2019 (Knoema 2021). The average cost of a ton is approximately UGX 2,400,000 (US$ 675). 1 By adopting these figures, the estimated loss on rice at varying return periods was computed, as shown in Fig. 10.
Destruction of plantations signals a high risk of food insecurity and poverty in the area. In 2019, gardens were destroyed in Tororo and Butaleja district leaving the community in looming fears of famine (Daily Monitor 2019). A similar event happened in 2020 (Uganda Radio Network 2020). With the above, communities may need to develop best management practices such as restoration of wetlands and embracing of cover crops (Antolini et al. 2020). Alternatively, farmers may need to opt cultivating low yielding traditional crops that are rather resistant to flooding instead of experiencing total losses. Besides, observation of flood trends may help famers adjust their agricultural activities in the flood prone areas. Figure 11a, b shows the inundated infrastructures (churches, mosques, academic institutions and health facilities). There was no noticeable difference in the number of inundated infrastructures at the six different return periods (2, 5, 10, 25, 50 and 100 years). Therefore, in Fig. 11a, b, only results obtained at 100-YRP are presented considering the extreme flow (Q). Considering the 100-YRP, fifteen academic institutions, twelve health facilities were inundated and remain susceptible to future inundation (Fig. 11a). Similarly, twelve mosques and twenty churches are at a risk of flooding (Fig. 11b).
Based on the available information, economic losses due to flooding of infrastructures were estimated following the procedure in Sect. 2.3.4. The estimated economic losses for  (Table 5). For higher return periods beyond 100 years, the risk needs to be assessed in future when adequate flow data is available. Improvement of catchment management can minimise the inundation of such key infrastructures. The major transport roads including highway are at a risk of inundation as presented in Fig. 12.
A total of about 4.6 km and 6.7 km of road network are prone to inundation considering the LQ and UQ at 2-year return period. At 100-year return period, the affected roads will increase to approximately 6.54 km and 8.43 km for LQ and UQ, respectively (Fig. 12). The affected sections of the roads have bridges and/or culverts. Inundation of the road infrastructure indicates probable inadequate capacity of the existing bridges and/or culverts. This puts the roads in danger of collapsing hence disrupting transport and livelihoods in the area. The authorities may use these findings as baseline information in development of more all-inclusive approaches for improving safety of these infrastructures.

Conclusion
The study area tends to experience dramatic flooding events, with several socio-economic impacts, such as loss of lives and destructions of physical infrastructure. Performing flood risk assessment is vital for predictive planning of risk-based adaptation measures. In this context, the 2D HEC-RAS model was employed on Ugandan River Malaba sub-catchment. Spatial extent of inundation on human settlement, land cover and several infrastructures were assessed with respect to return periods of extreme flow quantiles. HEC-RAS model results generally confirmed an acceptable fit with the observed data. The changes in LULC from 2010 to 2017 had no significant influence on the flow/flood simulations ( p > 0.05). Flooding was highly noticeable in the low-lying areas of Butaleja, Tororo and some parts of Namutumba and Bugiri. Cropland was and may continue experiencing the highest vulnerability followed by bare land. Large amounts of money losses are estimated due to inundated rice gardens. This signals a high risk of food insecurity and poverty in the study area. Therefore, fast growing crops and developing best agricultural management practices such as restoration of wetlands and embracing of cover crops could be adopted by farmers. Human settlement and several infrastructures (e.g., churches, mosques, academic institutions, health facilities, roads, etc.) are at a risk flooding. Churches displayed the highest economic loss overall. Therefore, there is need to sensitise the community on the dangers of setting-up structures in flood-prone areas. Besides, authorities should provide land in safe places for setting-up some of the communal infrastructures. In addition, several structural and non-structural sustainable management strategies including embankment/river training structures, flood forecasting and early warning, modification of farming practices, watershed management, relocation of personal property can help in reducing the flood-related loss of human life and property (Mubialiwo et al. 2021a).
A few limitations to this study which form recommendations for future research studies are worth mentioning. Future research studies should consider finer spatial resolution LULC data. Similarly, the use of high-resolution terrain data, for instance, by processing aerial images acquired by unmanned aerial vehicle (UAV) flights over the flood risk area is recommended for future research studies. To eliminate the bias from using assumed damage functions, future research studies are recommended when site-specific damage functions are developed. Besides, future research studies should explore applying flood loss models to estimate crop losses (Pinos et al. 2020;Tanir et al. 2021). Future research studies are recommended to advance the present study findings by performing hybrid 1D-2D model when the overflow locations information become available. This study provides information of the inundation extent and level of vulnerability across the study area which is vital for development of risk-based adaptation pathways. Government authorities can make use of this study findings to identify the modelled floodplains and vulnerability, which can be circulated to the local establishments for effective land-use planning and management.

B. Lumped Conceptual Rainfall-Runoff Models
See Appendix Table 7.

C. Trend Analyses
Trend magnitude or slope (m) in seasonal and annual rainfall-runoff discharge was calculated using (Theil 1950) and (Sen 1968) as per the following expression: where x j and x i are, respectively, the jth and ith value from the sample of size n such that 1 < i ≤ (n − 1) and 1 < j ≤ n. The null hypothesis H o (no trend) was tested using the method developed by (Onyutha 2021) such that for a given data set Y of sample size n, we can re-scale Y into series d y using where t y,i denotes the number of times the ith observation exceeds other data points in Y . In the same line, w y,i refers to number of times the ith data point appears within Y.
The trend statistic T was given by (Onyutha 2021) where The mean of T is zero and for large n the distribution of T is approximately normal with the variance of T or given by (Onyutha 2021) (C.5) V(T) = n n 2 − 1 12 .  Hydrological model focusing on sub-flows' variation (HMSV) (Onyutha 2019) 1 3 Published in partnership with CECCR at King Abdulaziz University The standardized test statistic Z which follows the standard normal distribution with mean (variance) of zero (one) is given by After application of Eq. (C.6) following correction of V(T) from the effect of scaling according to the procedure in Onyutha (2021), the H o (no trend) was rejected for , where Z ∕2 denotes the standard normal variate . at the selected ; otherwise the H o was not rejected at . In this study, was taken as 0.05 thereby corresponding with the Z value of 1.96.

D. Depth-Damage Functions
See Appendix Table 8.
Published in partnership with CECCR at King Abdulaziz University Table 8 Depth-damage function for each infrastructure at different return periods obtained for each location