Hydrologic Model
The hydrologic model of the upstream Imus River Basin is composed of 10 subbasins, 7 reaches, and 7 junctions (Figure 1). It was calibrated using the rainfall data from Storm Jolina (Figure 2) that occurred on August 25, 2017. The simulated discharge from August 25, 2017, 15:00 to August 26, 2017, 15:00 was then generated using the estimated baseline parameters and compared with the actual hydrograph observed at the Daang Hari Bridge (Figure 3).
The baseline parameters were then adjusted both manually and by autocalibration using the Optimization Manager of HEC-HMS. The new simulation run of the calibrated model (Figure 4) was tested with NSE, PBIAS, and RSR values of 0.903, -0.0376, and 0.30, respectively, indicating a very good fit between the observed data and simulated results but with a slight overestimation bias [25].
The calibrated model was then subjected to validation using the rainfall data from Storm Falcon (Figure 5) that occurred on July 17, 2019. The simulated hydrograph was compared to the observed flow at the Daang Hari Bridge from July 17, 2019, 00:00 to July 18, 2019, 00:00 as shown in Figure 6. The model obtained NSE, PBIAS, and RSR scores of 0.818 (very good fit), -0.044 (near-optimal, overestimation biased), and 0.4 (better model simulation), respectively. Upon validation, the hydrologic model of the upstream Imus River Basin performed well and can be used for accurate simulation of streamflow in the river basin.
Generated flood maps for Storm Jolina and Storm Falcon were compared to satellite images of the actual flooding obtained from Sentinel-2: MultiSpectral Instrument, Level 2A to verify the accuracy of the model as shown in Figures 7 and 8. Both storms brought only low-intensity rainfall causing low runoff from the upstream basin. Storm waters traveled down the river and filled the wetlands near the coast and did not cause much flooding on the areas located adjacent to the river.
Parametric Sensitivity Results
The initial abstraction (Ia), time of concentration (Tc), storage coefficient (R), recession constant (Rc), and initial discharge (Qi) were the parameters tested for their sensitivity to the results of the model. Based on the sensitivity test, changes in Ia did not affect both the peak discharge (Qp) and the time to peak (Tp) and caused a very small variation to the runoff volume (V). Changes in Rc and Qi did not also affect Qp, Tp and V. Changes in the time of concentration did not affect peak discharges but affected both the time to peak and total runoff volume. The Tc showed a direct relationship with Tp and inverse relation to V. On average, a ±10% change on Tc caused a 30-min advance or delay on Tp, and 0.24 mm (12,204 m3) increase or decrease in V. Storage coefficient is inversely proportional to Qp and V, and directly proportional to Tp. On average, a ±10% change in R caused a 1.2 m3/s increase or decrease on Qp, a 40-min advance or delay on Tp, and a 2.90 mm (147,637 m3) increase or decrease in V. Figure 9 shows generated hydrographs indicating the model sensitivity to each parameter.
Summary of the relative sensitivity coefficient of each parameter to the hydrograph properties of the model output was shown in Table 1. The relative sensitivity coefficient was used to facilitate the comparison between the parameters. A positive value indicates a direct relationship between the parameter and model output, and vice versa. The higher the absolute value, the greater is the effect of the parameter on the model output. Based on the analysis, the storage coefficient (R) had the greatest effect for all the three model outputs.
Table 1. Computed relative sensitivity coefficient for each parameter in respect to the hydrograph properties.
Parameter
|
Relative sensitivity (sr)
|
Qp
|
Tp
|
V
|
Ia
|
0.000
|
0.000
|
-0.006
|
Tc
|
-0.009
|
0.219
|
-0.016
|
R
|
-0.690
|
0.306
|
-0.198
|
Rc
|
0.000
|
0.000
|
0.000
|
Bi
|
0.000
|
0.000
|
0.000
|
Model Response to Land Cover Change
The projected change in land cover was demonstrated in terms of an increase in built-up areas only. Based on the 2010 and 2015 survey, the rate of built-up expansion was 2.08% annually. This modest rate of increase was assumed to continue for the next 10 years from the 2015 survey indicating a little over 20% increase from 2015 to 2025.
Results showed that the peak discharge may increase by 2.33% (Table 2) and runoff volume by 1.86% (Table 3) on average in the river basin. This increase was brought by the urbanization in the area that resulted in more impervious areas in the basin such as roofing, roads, sidewalks, pavements, and other artificial structures that caused less infiltration and hence greater flood peak and runoff [33,34].
Table 2. Change in peak discharge (Qp) in m3/s due to the land cover change.
Return period
|
2015
|
2025
|
Percent change
|
5-year
|
91.00
|
93.80
|
3.08
|
10-year
|
111.70
|
114.70
|
2.69
|
25-year
|
138.50
|
141.50
|
2.17
|
50-year
|
157.70
|
160.80
|
1.97
|
100-year
|
177.10
|
180.20
|
1.75
|
Mean
|
2.33
|
Table 3. Change in runoff volume (V) in mm due to the land cover change.
Return period
|
2015
|
2025
|
Percent change
|
5-Year
|
154.17
|
157.98
|
2.47
|
10-Year
|
193.85
|
197.92
|
2.10
|
25-Year
|
246.09
|
250.42
|
1.76
|
50-Year
|
283.98
|
288.47
|
1.58
|
100-Year
|
322.61
|
327.33
|
1.46
|
Mean
|
1.87
|
Simulation of Hydrograph for Normal and Climate-affected Conditions
Based on the climate projections of PAGASA, there will be a 12.67% and 15.67% increase from the normal annual rainfall amount in the province for RCP 4.5 (moderate-level risk) and RCP 8.5 (high-level risk), respectively. Using the normal and projected climate change hyetographs, synthetic hydrograph for each storm return period were generated at the Daang Hari Bridge. The simulation period was three days with a 10-minute time interval. Figure 10 shows the synthetic hydrograph for each return period under normal conditions while Figure 11 and Figure 12 show the synthetic hydrographs at different return periods under RCP 4.5 and RCP 8.5, respectively.
Hydraulic Model and Flood Inundation Maps
Areas affected by the flood were those situated in the cities of Imus and Bacoor as well as some portion of the municipality of Kawit. The majority of the higher portion of the floodplain may experience maximum flood depth ranging from 0.21 m to 1.0 m. The lowest portions of the floodplain which extends to the coast are considered as the most susceptible with possible maximum flood depth ranging from 1.01 to 3.5 m. Intensified monsoons, tropical cyclones, and rising sea levels will enhance the destructive effect of floods in areas near the coast. As the stormwater increases due to the intensified rains and climate change, more flooded areas were observed in Imus City with depths of 0.2 to 1 m. As the return period increases, the amount of floodwaters coming upstream increases resulting in potentially wider and deeper flooding downstream of the basin.
Due to the low elevation, the floodwaters concentrated along the coastal areas. Maximum flood depths on the lower reaches were almost the same for the normal conditions, RCP 4.5 and RCP 8.5 under the same return period. The overflow of floodwater increases as the return period increases and due to the increase of total rainfall projected under RCP 4.5 and 8.5. The generated flood maps in HEC-RAS for each hypothetical rainfall scenario are shown in Figures 13 to 17.