LiDAR Surveys and Flood Mapping of Bauang River PDF

Hazard Mapping of the Philippines Using LIDAR (Phil-LIDAR 1)
LiDAR Surveys and Flood Mapping of Bauang River
© University of the Philippines and the University of the Philippines Baguio 2017
Published by the UP Training Center for Applied Geodesy and Photogrammetry (TCAGP)
College of Engineering
University of the Philippines – Diliman
Quezon City
1101 PHILIPPINES
This research project is supported by the Department of Science and Technology (DOST) as part of its Grants-in-Aid
(GIA) Program and is to be cited as:
E.C. Paringit and C. Pascua, (Eds.). (2017), LiDAR Surveys and Flood Mapping of Bauang River. Quezon City: University
of the Philippines Training Center on Geodesy and Photogrammetry-302pp.
The text of this information may be copied and distributed for research and educational purposes with proper
acknowledgement. While every care is taken to ensure the accuracy of this publication, the UP TCAGP disclaims all
responsibility and all liability (including without limitation, liability in negligence) and costs which might incur as a
result of the materials in this publication being inaccurate or incomplete in any way and for any reason.
For questions/queries regarding this report, contact:
Dr. Chelo Pascua

Project Leader, Phil-LiDAR 1 Program
University of the Philippines Baguio
Baguio City, Philippines 2600
E-mail: [email protected]
Enrico C. Paringit, Dr. Eng.

Program Leader, PHIL-LiDAR 1 Program
University of the Philippines Diliman
Quezon City, Philippines 1101
E-mail: [email protected]
National Library of the Philippines

ISBN: 978-621-430-088-4
i
ii
TABLE OF CONTENTS
LIST OF TABLES................................................................................................................................... v
LIST OF FIGURES............................................................................................................................... ix
LIST OF ACRONYMS AND ABBREVIATIONS.....................................................................................xiii
CHAPTER 1: OVERVIEW OF THE PROGRAM AND BAUANG RIVER......................................................... 1
1.1 Background of the Phil-LiDAR 1 Program..................................................................................... 1
1.2 Overview of the Bauang River Basin............................................................................................ 1
CHAPTER 2: LIDAR DATA ACQUISITION OF THE BAUANG FLOODPLAIN................................................. 3
2.1 Flight Plans................................................................................................................................... 3
2.2 Ground Base Stations.................................................................................................................. 5
2.3 Flight Missions........................................................................................................................... 12
2.4 Survey Coverage........................................................................................................................ 14
CHAPTER 3: LIDAR DATA PROCESSING OF THE BAUANG FLOODPLAIN................................................ 17
3.1 Overview of the LIDAR Data Pre-Processing.............................................................................. 17
3.2 Transmittal of Acquired LiDAR Data........................................................................................... 18
3.3 Trajectory Computation............................................................................................................. 18
3.4 LiDAR Point Cloud Computation................................................................................................ 21
3.5 LiDAR Data Quality Checking..................................................................................................... 21
3.6 LiDAR Point Cloud Classification and Rasterization.................................................................... 27
3.7 LiDAR Image Processing and Orthophotograph Rectification.................................................... 30
3.8 DEM Editing and Hydro-Correction........................................................................................... 31
3.9 Mosaicking of Blocks................................................................................................................. 34
3.10 Calibration and Validation of Mosaicked LiDAR DEM.............................................................. 36
3.11 Integration of Bathymetric Data into the LiDAR Digital Terrain Model.................................... 40
3.12 Feature Extraction................................................................................................................... 42
3.12.1 Quality Checking of Digitized Features’ Boundary.......................................................... 42
3.12.2 Height Extraction............................................................................................................. 42
3.12.3 Feature Attribution.......................................................................................................... 43
3.12.4 Final Quality Checking of Extracted Features.................................................................. 44
CHAPTER 4: LIDAR VALIDATION SURVEY AND MEASUREMENTS OF THE BAUANG RIVER BASIN.......... 45
4.1 Summary of Activities................................................................................................................ 45
4.2 Control Survey........................................................................................................................... 47
4.3 Baseline Processing................................................................................................................... 52
4.4 Network Adjustment................................................................................................................. 53
4.5 Cross-section and Bridge As-Built Survey and Water Level Marking......................................... 57
4.6 Validation Points Acquisition Survey.......................................................................................... 65
4.7 Bathymetric Survey.................................................................................................................... 67
CHAPTER 5: FLOOD MODELING AND MAPPING................................................................................. 71
5.1 Data Used for Hydrologic Modeling........................................................................................... 71
5.1.1 Hydrometry and Rating Curves......................................................................................... 71
5.1.2 Precipitation...................................................................................................................... 71
5.1.3 Rating Curves and River Outflow....................................................................................... 72
5.2 RIDF Station............................................................................................................................... 74
5.3 HMS Model................................................................................................................................ 76
5.4 Cross-section Data..................................................................................................................... 81
5.5 Flo 2D Model............................................................................................................................. 82
5.6 Results of HMS Calibration........................................................................................................ 83
5.7 Calculated outflow hydrographys and Discharge values for different rainfall return periods... 87
5.7.1 Hydrograph using the Rainfall Runoff Model.................................................................... 87
5.8 River Analysis (RAS) Model Simulation...................................................................................... 88
5.9 Flow Depth and Flood Hazard.................................................................................................... 89
5.10 Inventory of Areas Exposed to Flooding of Affected Areas...................................................... 96
5.11 Flood Validation..................................................................................................................... 227
REFERENCES ................................................................................................................................. 230
ANNEXES ................................................................................................................................. 231
Annex 1. Technical Specifications of the Pegasus LiDAR Sensor used in the Bauang
Floodplain Survey.......................................................................................................................... 231
Annex 2. NAMRIA Certification of Reference Points used in the LiDAR Survey............................. 232
Annex 3. Baseline Processing Reports of Control Points used in the LiDAR Survey....................... 238
Annex 4. The LiDAR Survey Team Composition............................................................................. 239
iii
Annex 5. Data Transfer Sheets for the Bauang Floodplain Flights................................................. 240
Annex 6. Flight Logs for the Flight Missions................................................................................... 241
Annex 7. Flight Status Reports....................................................................................................... 250
Annex 8. Mission Summary Reports.............................................................................................. 260
Annex 9. Bauang Model Basin Parameters.................................................................................... 280
Annex 10. Bauang Model Reach Parameters................................................................................. 283
Annex 11. Bauang Field Validation Points...................................................................................... 285
Annex 12. Educational Institutions Affected by Flooding in Bauang Floodplain............................ 290
Annex 13. Medical Institutions Affected by Flooding in Bauang Floodplain.................................. 292
iv
LIST OF TABLES
Table 1. Flight planning parameters for the Pegasus LiDAR system............................................................ 3
Table 2. Details of the recovered NAMRIA horizontal control point LUN-62, used as a base station
for the LiDAR acquisition................................................................................................................... 6
Table 5. Details of the established horizontal control point LUN-3047, used as a base station for the
LiDAR acquisition............................................................................................................................... 9
for the LiDAR acquisition................................................................................................................. 10
Table 7. Details of the recovered NAMRIA vertical control point LU-94, used as a base station for the
LiDAR acquisition........................................................................................................................ 11
Table 8. Ground control points used during the LiDAR data acquisition................................................... 12
Table 9. Flight missions for the LiDAR data acquisition in the Bauang floodplain..................................... 13
Table 10. Actual parameters used during the LiDAR data acquisition....................................................... 14
Table 11. List of municipalities and cities surveyed during the Bauang floodplain LiDAR survey............. 15
Table 12. Self-calibration results values for the Bauang flights................................................................. 21
Table 13. List of LiDAR blocks for the Bauang floodplain........................................................................... 23
Table 14. Bauang classification results in TerraScan.................................................................................. 27
Table 15. LiDAR blocks with their corresponding areas............................................................................. 32
Table 16. Shift values of each LiDAR block of the Bauang floodplain........................................................ 34
Table 17. Calibration statistical measures................................................................................................. 38
Table 18. Validation statistical measures................................................................................................... 39
Table 19. Quality checking ratings for the Bauang building features........................................................ 42
Table 20. Building features extracted for the Bauang floodplain.............................................................. 43
Table 21. Total length of extracted roads for the Bauang floodplain......................................................... 43
Table 22. Number of extracted water bodies for the Bauang floodplain.................................................. 43
Table 23. List of reference and control points occupied in the Bauang River Survey (Source: NAMRIA;
UP-TCAGP)....................................................................................................................................... 47
Table 24. Baseline processing report for the Bauang River survey............................................................ 53
Table 25. Constraints applied to the adjustments of the control points................................................... 54
Table 26. Adjusted grid coordinates for the control points used in the Bauang floodplain survey........... 54
Table 27. Adjusted geodetic coordinates for control points used in the Bauang River floodplain
validation......................................................................................................................................... 56
Table 28. Reference and control points used in the Bauang River Static Survey, with their
corresponding locations (Source: NAMRIA, UP-TCAGP).................................................................. 56
Table 29. RIDF values for the Baguio Rain Gauge, computed by PAGASA................................................. 75
Table 30. Range of calibrated values for the Bauang River Basin model................................................... 85
Table 31. Efficiency Test of the Bauang HMS Model.................................................................................. 85
Table 32. Peak values of the Bauang HEC-HMS Model outflow, using the Baguio RIDF............................ 88
Table 33. Municipalities affected in the Bauang floodplain....................................................................... 89
Table 34. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - A...................... 97
Table 35. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - B...................... 97
Table 36. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - C...................... 98
Table 37. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - D..................... 98
Table 38. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - E...................... 99
Table 39. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - F.................... 100
Table 40. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - G................... 100
Table 41. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - H................... 101
Table 42. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - I..................... 101
Table 43. Affected areas in Kapangan, Benguet during a 5-year rainfall return period........................... 106
Table 44. Affected areas in La Trinidad, Benguet during a 5-year rainfall return period - A.................... 108
Table 45. Affected areas in La Trinidad, Benguet during a 5-year rainfall return period - B.................... 108
Table 46. Affected areas in Sablan, Benguet during a 5-year rainfall return period................................ 110
Table 47. Affected areas in Tuba, Benguet during a 5-year rainfall return period................................... 112
Table 48. Affected areas in Tublay, Benguet during a 5-year rainfall return period................................ 114
Table 49. Affected areas in Agoo, La Union during a 5-year rainfall return period.................................. 115
Table 50. Affected areas in Aringay, La Union during a 5-year rainfall return period - A......................... 117
Table 51. Affected areas in Aringay, La Union during a 5-year rainfall return period- B.......................... 117
v
Table 52. Affected areas in Bagulin, La Union during a 5-year rainfall return period.............................. 119
Table 53. Affected areas in Bauang, La Union during a 5-year rainfall return period - A......................... 121
Table 54. Affected areas in Bauang, La Union during a 5-year rainfall return period - B......................... 121
Table 55. Affected areas in Bauang, La Union during a 5-year rainfall return period - C......................... 122
Table 56. Affected areas in Bauang, La Union during a 5-year rainfall return period - D........................ 122
Table 57. Affected areas in Burgos, La Union during a 5-year rainfall return period............................... 125
Table 58. Affected areas in Caba, La Union during a 5-year rainfall return period - A............................. 127
Table 59. Affected areas in Caba, La Union during a 5-year rainfall return period - B............................. 127
Table 60. Affected areas in Naguilian, La Union during a 5-year rainfall return period - A..................... 129
Table 61. Affected areas in Naguilian, La Union during a 5-year rainfall return period - B...................... 129
Table 62. Affected areas in Naguilian, La Union during a 5-year rainfall return period - C...................... 130
Table 63. Affected areas in Naguilian, La Union during a 5-year rainfall return period - D..................... 130
Table 64. Affected areas in Pugo, La Union during a 5-year rainfall return period.................................. 133
Table 65. Affected areas in Rosario, La Union during a 5-year rainfall return period.............................. 135
Table 66. Affected areas in San Fernando City, La Union during a 5-year rainfall return period - A........ 136
Table 67. Affected areas in San Fernando City, La Union during a 5-year rainfall return period - B........ 136
Table 68. Affected areas in San Gabriel, La Union during a 5-year rainfall return period....................... 138
Table 69. Affected areas in Santo Tomas, La Union during a 5-year rainfall return period...................... 139
Table 70. Affected areas in Tubao, La Union during a 5-year rainfall return period - A........................... 140
Table 71. Affected areas in Tubao, La Union during a 5-year rainfall return period - B........................... 140
Table 72. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - A.................. 142
Table 73. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - B.................. 142
Table 74. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - C.................. 143
Table 75. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - D................. 143
Table 76. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - E.................. 144
Table 77. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - F.................. 144
Table 78. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - G................. 145
Table 79. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - H................. 145
Table 80. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - I................... 146
Table 81. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - J................... 146
Table 82. Affected areas in Kapangan, Benguet during a 25-year rainfall return period........................ 152
Table 83. Affected areas in La Trinidad, Benguet during a 25-year rainfall return period - A.................. 153
Table 84. Affected areas in La Trinidad, Benguet during a 25-year rainfall return period - B.................. 153
Table 85. Affected areas in Sablan, Benguet during a 25-year rainfall return period.............................. 155
Table 86. Affected areas in Tuba, Benguet during a 25-year rainfall return period................................. 156
Table 87. Affected areas in Tublay, Benguet during a 25-year rainfall return period.............................. 157
Table 88. Affected areas in Agoo, La Union during a 25-year rainfall return period................................ 158
Table 89. Affected areas in Aringay, La Union during a 25-year rainfall return period - A....................... 159
Table 90. Affected areas in Aringay, La Union during a 25-year rainfall return period - B....................... 159
Table 91. Affected areas in Bagulin, La Union during a 25-year rainfall return period............................ 161
Table 92. Affected areas in Bauang, La Union during a 25-year rainfall return period - A....................... 162
Table 93. Affected areas in Bauang, La Union during a 25-year rainfall return period - B....................... 162
Table 94. Affected areas in Bauang, La Union during a 25-year rainfall return period - C....................... 163
Table 95. Affected areas in Bauang, La Union during a 25-year rainfall return period - D...................... 163
Table 96. Affected areas in Burgos, La Union during a 25-year rainfall return period............................. 166
Table 97. Affected areas in Caba, La Union during a 25-year rainfall return period - A........................... 168
Table 98. Affected areas in Caba, La Union during a 25-year rainfall return period - B........................... 168
Table 99. Affected areas in Naguilian, La Union during a 25-year rainfall return period - A................... 170
Table 100. Affected areas in Naguilian, La Union during a 25-year rainfall return period - B.................. 170
Table 101. Affected areas in Naguilian, La Union during a 25-year rainfall return period - C.................. 171
Table 102. Affected areas in Naguilian, La Union during a 25-year rainfall return period - D................. 171
Table 103. Affected areas in Pugo, La Union during a 25-year rainfall return period.............................. 174
Table 104. Affected areas in Rosario, La Union during a 25-year rainfall return period.......................... 176
Table 105. Affected areas in San Fernando City, La Union during a 25-year rainfall return period - A.... 177
Table 106. Affected areas in San Fernando City, La Union during a 25-year rainfall return period - B.... 177
Table 107. Affected areas in San Gabriel, La Union during a 25-year rainfall return period................... 179
Table 108. Affected areas in Santo Tomas, La Union during a 25-year rainfall return period.................. 180
Table 109. Affected areas in Tubao, La Union during a 25-year rainfall return period - A....................... 181
Table 110. Affected areas in Tubao, La Union during a 25-year rainfall return period - B....................... 181
Table 111. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - A............. 183
Table 112. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - B.............. 183
Table 113. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - C.............. 184
Table 114. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - D............. 184
vi
Table 115. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - E.............. 185
Table 116. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - F.............. 185
Table 117. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - G............. 186
Table 118. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - H............. 186
Table 119. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - I............... 187
Table 120. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - J............... 187
Table 121. Affected areas in Kapangan, Benguet during a 100-year rainfall return period..................... 193
Table 122. Affected areas in La Trinidad, Benguet during a 100-year rainfall return period - A.............. 194
Table 123. Affected areas in La Trinidad, Benguet during a 100-year rainfall return period - B.............. 194
Table 124. Affected areas in Sablan, Benguet during a 100-year rainfall return period.......................... 196
Table 125. Affected areas in Tuba, Benguet during a 100-year rainfall return period............................. 197
Table 126. Affected areas in Tublay, Benguet during a 100-year rainfall return period.......................... 198
Table 127. Affected areas in Agoo, La Union during a 100-year rainfall return period............................ 199
Table 128. Affected areas in Aringay, La Union during a 100-year rainfall return period - A................... 200
Table 129. Affected areas in Aringay, La Union during a 100-year rainfall return period - B................... 200
Table 130. Affected areas in Bagulin, La Union during a 100-year rainfall return period........................ 202
Table 131. Affected areas in Bauang, La Union during a 100-year rainfall return period - A................... 204
Table 132. Affected areas in Bauang, La Union during a 100-year rainfall return period - B................... 204
Table 133. Affected areas in Bauang, La Union during a 100-year rainfall return period - C................... 205
Table 134. Affected areas in Bauang, La Union during a 100-year rainfall return period - D.................. 205
Table 135. Affected areas in Burgos, La Union during a 100-year rainfall return period......................... 208
Table 136. Affected areas in Caba, La Union during a 100-year rainfall return period - A....................... 210
Table 137. Affected areas in Caba, La Union during a 100-year rainfall return period - B....................... 210
Table 138. Affected areas in Naguilian, La Union during a 100-year rainfall return period - A............... 212
Table 139. Affected areas in Naguilian, La Union during a 100-year rainfall return period - B................ 212
Table 140. Affected areas in Naguilian, La Union during a 100-year rainfall return period - C................ 213
Table 141. Affected areas in Naguilian, La Union during a 100-year rainfall return period - D............... 213
Table 142. Affected areas in Pugo, La Union during a 100-year rainfall return period............................ 216
Table 143. Affected areas in Rosario, La Union during a 100-year rainfall return period........................ 218
Table 144. Affected areas in San Fernando City, La Union during a 100-year rainfall return period - A.. 219
Table 145. Affected areas in San Fernando City, La Union during a 100-year rainfall return period - B.. 219
Table 146. Affected areas in San Gabriel, La Union during a 100-year rainfall return period................. 221
Table 147. Affected areas in Santo Tomas, La Union during a 100-Year rainfall return period................ 222
Table 148. Affected areas in Tubao, La Union during a 100-year rainfall return period - A..................... 223
Table 149. Affected areas in Tubao, La Union during a 100-year rainfall return period - B..................... 223
Table 150. Area covered by each warning level, with respect to the rainfall scenario............................ 226
Table 151. Actual flood depth vs. simulated flood depth in the Bauang River Basin.............................. 229
Table 152. Summary of Accuracy Assessment in the Bauang River Basin............................................... 229
vii
viii
LIST OF FIGURES
Figure 1. Location map of the Bauang River Basin (in brown)..................................................................... 2
Figure 2. Flight plans and base stations used to cover the Bauang floodplain survey................................ 4
Figure 3. (a) GPS set-up over LUN-62, as recovered 15 m south from the first access ladder of the river
control and about 100 m north from the end. It is also situated 300 m south of a hanging
bridge. (b) NAMRIA reference point LUN-62, as recovered by the field team........................... 6
Figure 4. (a) GPS set-up over LUN-176, recovered near a corner of a farm dike, about 15 m southeast
of the well, and about 20 m southwest of the nearest house. (b) NAMRIA reference point
LUN-176, as recovered by the field team.................................................................................. 7
Figure 5. (a) GPS set-up over LUN-3062, recovered at the top of a dike, approximately 100 m north
of Philippine Central College of Arts Sciences & Technology, and 80 m north of Naguilian
emission testing center in Barangay Natividad, Naguilian, La Union. (b) NAMRIA reference
point LUN-3062, as recovered by the field team....................................................................... 8
Figure 6. LUN-3047, as recovered in Barangay Nazareno, Agoo, La Union.................................................. 9
Figure 7. (a) GPS set-up over LUN-3129, recovered beside the National Road about 50 meters
northeast of the nearest house. (b) NAMRIA reference point LUN-3129, as recovered
by the field team...................................................................................................................... 10
Figure 8. NAMRIA reference point LU-94, as recovered by the field team................................................ 11
Figure 9. Actual LiDAR survey coverage for the Bauang floodplain........................................................... 16
Figure 10. Schematic diagram for the Data Pre-Processing Component................................................... 18
Figure 11. Smoothed Performance Metric Parameters of Bauang Flight 1153P...................................... 19
Figure 12. Solution Status Parameters of Bauang Flight 1153P................................................................. 20
Figure 13. The best estimated trajectory conducted over the Bauang floodplain.................................... 21
Figure 14. Boundaries of the processed LiDAR data over the Bauang floodplain..................................... 22
Figure 15. Image of data overlap for Bauang floodplain........................................................................... 24
Figure 16. Density map of merged LiDAR data for the Bauang floodplain................................................ 25
Figure 17. Elevation difference map between flight lines for the Bauang floodplain................................ 26
Figure 18. Quality checking for Bauang Flight 1153P, using the Profile Tool of QT Modeler..................... 27
Figure 19. (a) Tiles for Bauang floodplain, and (b) classification results in TerraScan............................... 28
Figure 20. Point cloud (a) before and (b) after classification..................................................................... 29
Figure 21. The production of (a) last return DSM and (b) DTM, (c) first return DSM and (d) secondary
DTM in some portion of Bauang floodplain............................................................................ 30
Figure 22. The Bauang floodplain, with available orthophotographs....................................................... 31
Figure 23. Sample orthophotograph tiles for the Bauang floodplain........................................................ 31
Figure 24. Portions in the DTM of the Bauang floodplain – a bridge (a) before and (b) after manual
editing; an embankment (c) before and (d) after data retrieval; and a building (e) before
and (f) after manual editing.................................................................................................... 33
Figure 25. Map of processed LiDAR data for the Bauang floodplain......................................................... 35
Figure 26. Map of the Bauang floodplain, with the validation survey points in green.............................. 37
Figure 27. Correlation plot between the calibration survey points and the LiDAR data........................... 38
Figure 28. Correlation plot between the validation survey points and the LiDAR data............................. 39
Figure 29. Map of the Bauang floodplain, with bathymetric survey points shown in blue....................... 41
Figure 30. Blocks (in blue) of Bauang building features that were subjected to QC................................. 42
Figure 31. Extracted features for the Bauang floodplain........................................................................... 44
Figure 32. Extent of the bathymetric survey (in blue line) in Bauang River and the LiDAR data
validation survey (in red)......................................................................................................... 46
Figure 33. GNSS network of the Bauang River field survey....................................................................... 48
Figure 34. GNSS receiver Trimble® SPS 985 at LU-67, located at the Bauang Bridge approach in
Barangay Nagrebcan, Municipality of Bauang, La Union......................................................... 49
Figure 35. GNSS receiver Trimble® SPS 985 at LUN-61, located near the BPPC Housing in Barangay
Guerrero, Municipality of Bauang, La Union........................................................................... 49
Figure 36. GNSS receiver Trimble® SPS 882 at UP-AGY, located at the Aringat Bridge approach in
Barangay Poblacion, Municipality of Aringay, La Union........................................................... 50
Figure 37. GNSS receiver Trimble® SPS 882 at UP-CBN, located at the Cabaroan Bridge approach in
Barangay Cabaroan, Municipality of San Juan, La Union......................................................... 50
Figure 38. GNSS receiver Trimble® SPS 985 at UP-GST, located at the Guesset Bridge approach in
Barangay Ambaracao, Municipality of Naguilian, La Union..................................................... 51
Figure 39. GNSS receiver Trimble® SPS 985 at UP-NG, located at the Naguilian Bridge approach in
Barangay Cabaritan, Municipality of Naguilian, La Union........................................................ 51
Figure 40. Cross-section survey at the Guesset Bridge............................................................................. 57
Figure 41. Cross-section survey at the Mamat-ing Bridge......................................................................... 58
Figure 42. Guesset Bridge cross-section location map.............................................................................. 59
Figure 43. Bridge cross-section diagram of the Guesset Bridge................................................................ 60
Figure 44. Mamat-ing Bridge cross-section loaction map......................................................................... 61
ix
Figure 45. Bridge cross-section diagram of the Mamat-ing Bridge........................................................... 62

Figure 46. Bridge data form for the Guesset Bridge.................................................................................. 63
Figure 47. Bridge data form for the Mamat-ing Bridge............................................................................. 64
Figure 48. Set-up for the validation points acquisition in the Bauang River Basin.................................... 65
Figure 49. Extent of the LiDAR ground validation survey of the Bauang River Basin................................ 66
Figure 50. Bathymetry survey set-up for Bauang River............................................................................. 67
Figure 51. Extent of the bathymetric survey of the Bauang River Basin................................................... 69
Figure 52. Bauang centerline riverbed profile........................................................................................... 70
Figure 53. The location map of the Bauang HEC-HMS model, which was used for calibration................. 72
Figure 54. Cross-section plot of the Mamat-Ing Bridge............................................................................. 73
Figure 55. Rating curve at the Mamat-Ing Bridge, San Fernando City, La Union....................................... 73
Figure 56. Rainfall and outflow data at the Mamat-Ing Bridge, which were used for modeling............... 74
Figure 57. Baguio RIDF location relative to the Bauang River Basin.......................................................... 75
Figure 58. Synthetic storm generated from a 24-hour period rainfall, for various return periods............ 76
Figure 59. Soil map of the Bauang River Basin (Source: DA)..................................................................... 77
Figure 60. Land cover map of the Bauang River Basin (Source: NAMRIA)................................................ 78
Figure 61. Slope map of the Bauang River Basin....................................................................................... 79
Figure 62. Stream delineation map of the Bauang River Basin................................................................. 80
Figure 63. The Bauang River Basin model, generated using HEC-HMS..................................................... 81
Figure 64. River cross-section of the Bauang River, generated through the ArcMap HEC GeoRAS tool... 82
Figure 65. A screenshot of the river sub-catchment with the computational area to be modeled in
FLO-2D Grid Developer System Pro (FLO-2D GDS Pro)............................................................ 83
Figure 66. Outflow hydrograph of Bauang produced by the HEC-HMS model, compared with
observed outflow..................................................................................................................... 84
Figure 67. Outflow hydrograph at the Bauang Station, generated using the Baguio RIDF simulated
in HEC-HMS.............................................................................................................................. 87
Figure 68. Sample output map of the Bauang RAS Model........................................................................ 88
Figure 69. 100-year flood hazard map for the Bauang floodplain, overlaid on Google Earth imagery..... 90
Figure 70. 100-year flow depth map for the Bauang floodplain, overlaid on Google Earth imagery........ 91
Figure 71. 25-year flood hazard map for the Bauang floodplain, overlaid on Google Earth imagery....... 92
Figure 72. 25-year flow depth map for the Bauang floodplain, overlaid on Google Earth imagery.......... 93
Figure 73. 5-year flood hazard map for the Bauang floodplain, overlaid on Google Earth imagery......... 94
Figure 74. 5-year flow depth map for the Bauang floodplain, overlaid on Google Earth imagery............ 95
Figure 75. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - A.................. 102
Figure 76. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - B.................. 102
Figure 77. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - C.................. 103
Figure 78. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - D.................. 103
Figure 79. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - E.................. 104
Figure 80. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - F................... 104
Figure 81. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - G.................. 105
Figure 82. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - H.................. 105
Figure 83. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - I................... 106
Figure 84. Affected areas in Kapangan, Benguet during a 5-year rainfall return period......................... 107
Figure 85. Affected areas in La Trinidad, Benguet during a 5-year rainfall return period - A.................. 109
Figure 86. Affected areas in La Trinidad, Benguet during a 5-year rainfall return period - B.................. 109
Figure 87. Affected areas in Sablan, Benguet during a 5-year rainfall return period.............................. 111
Figure 88. Affected areas in Tuba, Benguet during a 5-year rainfall return period................................. 113
Figure 89. Affected areas in Tublay, Benguet during a 5-year rainfall return period............................... 115
Figure 90. Affected areas in Agoo, La Union during a 5-year rainfall return period................................ 116
Figure 91. Affected areas in Aringay, La Union during a 5-year rainfall return period - A....................... 118
Figure 92. Affected areas in Aringay, La Union during a 5-year rainfall return period - B....................... 118
Figure 93. Affected areas in Bagulin, La Union during a 5-year rainfall return period............................ 120
Figure 94. Affected areas in Bauang, La Union during a 5-year rainfall return period - A....................... 123
Figure 95. Affected areas in Bauang, La Union during a 5-year rainfall return period - B....................... 123
Figure 96. Affected areas in Bauang, La Union during a 5-year rainfall return period - C....................... 124
Figure 97. Affected areas in Bauang, La Union during a 5-year rainfall return period - D....................... 124
Figure 98. Affected areas in Burgos, La Union during a 5-year rainfall return period............................. 126
Figure 99. Affected areas in Caba, La Union during a 5-year rainfall return period - A........................... 128
Figure 100. Affected areas in Caba, La Union during a 5-year rainfall return period - B......................... 128
Figure 101. Affected areas in Naguilian, La Union during a 5-year rainfall return period - A.................. 131
Figure 102. Affected areas in Naguilian, La Union during a 5-year rainfall return period - B.................. 131
Figure 103. Affected areas in Naguilian, La Union during a 5-year rainfall return period - C.................. 132
Figure 104. Affected areas in Naguilian, La Union during a 5-year rainfall return period - D.................. 132
Figure 105. Affected areas in Pugo, La Union during a 5-year rainfall return period.............................. 134
Figure 106. Affected areas in Rosario, La Union during a 5-year rainfall return period.......................... 135
Figure 107. Affected areas in San Fernando City, La Union during a 5-year rainfall return period - A.... 137
x
Figure 108. Affected areas in San Fernando City, La Union during a 5-year rainfall return period - B.... 137
Figure 109. Affected areas in San Gabriel, La Union during a 5-year rainfall return period.................... 138
Figure 110. Affected areas in Santo Tomas, La Union during a 5-year rainfall return period.................. 139
Figure 111. Affected areas in Tubao, La Union during a 5-year rainfall return period - A....................... 141
Figure 112. Affected areas in Tubao, La Union during a 5-year rainfall return period - B....................... 141
Figure 113. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - A.............. 147
Figure 114. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - B.............. 147
Figure 115. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - C.............. 148
Figure 116. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - D.............. 148
Figure 117. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - E.............. 149
Figure 118. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - F.............. 149
Figure 119. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - G.............. 150
Figure 120. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - H.............. 150
Figure 121. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - I............... 151
Figure 122. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - J............... 151
Figure 123. Affected areas in Kapangan, Benguet during a 25-year rainfall return period.................... 152
Figure 124. Affected areas in La Trinidad, Benguet during a 25-year rainfall return period - A.............. 154
Figure 125. Affected areas in La Trinidad, Benguet during a 25-year rainfall return period - B.............. 154
Figure 126. Affected areas in Sablan, Benguet during a 25-year rainfall return period.......................... 155
Figure 127. Affected areas in Tuba, Benguet during a 25-year rainfall return period............................. 156
Figure 128. Affected areas in Tublay, Benguet during a 25-year rainfall return period........................... 157
Figure 129. Affected areas in Agoo, La Union during a 25-year rainfall return period............................ 158
Figure 130. Affected areas in Aringay, La Union during a 25-year rainfall return period - A................... 160
Figure 131. Affected areas in Aringay, La Union during a 25-year rainfall return period - B................... 160
Figure 132. Affected areas in Bagulin, La Union during a 25-year rainfall return period........................ 161
Figure 133. Affected areas in Bauang, La Union during a 25-year rainfall return period - A................... 164
Figure 134. Affected areas in Bauang, La Union during a 25-year rainfall return period - B................... 164
Figure 135. Affected areas in Bauang, La Union during a 25-year rainfall return period - C................... 165
Figure 136 - Affected areas in Bauang, La Union during a 25-year rainfall return period - D.................. 165
Figure 137. Affected areas in Burgos, La Union during a 25-year rainfall return period......................... 167
Figure 138. Affected areas in Caba, La Union during a 25-year rainfall return period - A....................... 169
Figure 139. Affected areas in Caba, La Union during a 25-year rainfall return period - B....................... 169
Figure 140. Affected areas in Naguilian, La Union during a 25-year rainfall return period - A................ 172
Figure 141. Affected areas in Naguilian, La Union during a 25-year rainfall return period - B................ 172
Figure 142. Affected areas in Naguilian, La Union during a 25-year rainfall return period - C................ 173
Figure 143. Affected areas in Naguilian, La Union during a 25-year rainfall return period - D................ 173
Figure 144. Affected areas in Pugo, La Union during a 25-year rainfall return period............................ 175
Figure 145. Affected areas in Rosario, La Union during a 25-year rainfall return period........................ 176
Figure 146. Affected areas in San Fernando City, La Union during a 25-year rainfall return period - A.. 178
Figure 147. Affected areas in San Fernando City, La Union during a 25-year rainfall return period - B.. 178
Figure 148. Affected areas in San Gabriel, La Union during a 25-year rainfall return period.................. 179
Figure 149. Affected areas in Santo Tomas, La Union during a 25-year rainfall return period................ 180
Figure 150. Affected areas in Tubao, La Union during a 25-year rainfall return period - A..................... 182
Figure 151. Affected areas in Tubao, La Union during a 25-year rainfall return period - B..................... 182
Figure 152. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - A............ 188
Figure 153. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - B............ 188
Figure 154. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - C............ 189
Figure 155. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - D............ 189
Figure 156. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - E............ 190
Figure 157. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - F............ 190
Figure 158. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - G............ 191
Figure 159. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - H............ 191
Figure 160. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - I............. 192
Figure 161. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - J............. 192
Figure 162. Affected areas in Kapangan, Benguet during a 100-year rainfall return period................... 193
Figure 163. Affected areas in La Trinidad, Benguet during a 100-year rainfall return period - A............ 195
Figure 164. Affected areas in La Trinidad, Benguet during a 100-year rainfall return period - B............ 195
Figure 165. Affected areas in Sablan, Benguet during a 100-year rainfall return period........................ 196
Figure 166. Affected areas in Tuba, Benguet during a 100-year rainfall return period........................... 197
Figure 167. Affected areas in Tublay, Benguet during a 100-year rainfall return period......................... 198
Figure 168. Affected areas in Agoo, La Union during a 100-year rainfall return period.......................... 199
Figure 169. Affected areas in Aringay, La Union during a 100-year rainfall return period - A................. 201
Figure 170. Affected areas in Aringay, La Union during a 100-year rainfall return period - B................. 201
Figure 171. Affected areas in Bagulin, La Union during a 100-year rainfall return period...................... 203
Figure 172. Affected areas in Bauang, La Union during a 100-year rainfall return period - A................. 206
Figure 173. Affected areas in Bauang, La Union during a 100-year rainfall return period - B................ 206
xi
Figure 174. Affected areas in Bauang, La Union during a 100-year rainfall return period - C................. 207
Figure 175. Affected areas in Bauang, La Union during a 100-year rainfall return period - D................. 207
Figure 176. Affected areas in Burgos, La Union during a 100-year rainfall return period....................... 209
Figure 177. Affected areas in Caba, La Union during a 100-year rainfall return period - A..................... 211
Figure 178. Affected areas in Caba, La Union during a 100-year rainfall return period - B..................... 211
Figure 179. Affected areas in Naguilian, La Union during a 100-year rainfall return period - A.............. 214
Figure 180. Affected areas in Naguilian, La Union during a 100-year rainfall return period - B.............. 214
Figure 181. Affected areas in Naguilian, La Union during a 100-year rainfall return period - C.............. 215
Figure 182. Affected areas in Naguilian, La Union during a 100-year rainfall return period - D.............. 215
Figure 183. Affected areas in Pugo, La Union during a 100-year rainfall return period.......................... 217
Figure 184. Affected areas in Rosario, La Union during a 100-year rainfall return period...................... 218
Figure 185. Affected areas in San Fernando City, La Union during a 100-year rainfall return period - A.220
Figure 186. Affected areas in San Fernando City, La Union during a 100-year rainfall return period - B.220
Figure 187. Affected areas in San Gabriel, La Union during a 100-year rainfall return period................ 221
Figure 188. Affected areas in Santo Tomas, La Union during a 100-Year rainfall return period.............. 222
Figure 189. Affected areas in Tubao, La Union during a 100-year rainfall return period - A................... 224
Figure 190. Affected areas in Tubao, La Union during a 100-year rainfall return period - B................... 224
Figure 191. Flood validation points of the Bauang River Basin............................................................... 228
Figure 192. Flood map depth vs. actual flood depth for Bauang............................................................ 228
xii
LIST OF ACRONYMS AND ABBREVIATIONS

AAC Asian Aerospace Corporation LC Low Chord
Ab abutment LGU local government unit
ALTM Airborne LiDAR Terrain Mapper LiDAR Light Detection and Ranging
ARG automatic rain gauge LMS LiDAR Mapping Suite
ASTI Advanced Science and Technology m AGL meters Above Ground Level
Institute MMS Mobile Mapping Suite
AWLS Automated Water Level Sensor MSL mean sea level
BA Bridge Approach NAMRIA National Mapping and Resource
BM benchmark Information Authority
BSWM Bureau of Soils and Water NSO National Statistics Office
Management NSTC Northern Subtropical Convergence
CAD Computer-Aided Design PAF Philippine Air Force
CN Curve Number PAGASA Philippine Atmospheric
CSRS Chief Science Research Specialist Geophysical and Astronomical
DA Department of Agriculture Services Administration
DAC Data Acquisition Component PDOP Positional Dilution of Precision
DEM Digital Elevation Model PPK Post-Processed Kinematic
[technique]
DENR Department of Environment and
Natural Resources PRF Pulse Repetition Frequency
DOST Department of Science and PTM Philippine Transverse Mercator
Technology QC Quality Check
DPPC Data Pre-Processing Component QT Quick Terrain [Modeler]
DREAM Disaster Risk and Exposure RA Research Associate
Assessment for Mitigation RBCO River Basin Control Office
[Program]
RIDF Rainfall-Intensity-Duration-
DRRM Disaster Risk Reduction and Frequency
Management
RMSE Root Mean Square Error
DSM Digital Surface Model
SAR Synthetic Aperture Radar
DTM Digital Terrain Model
SCS Soil Conservation Service
DVBC Data Validation and Bathymetry
Component SRTM Shuttle Radar Topography Mission
FMC Flood Modeling Component SRS Science Research Specialist
FOV Field of View SSG Special Service Group
GiA Grants-in-Aid TBC Thermal Barrier Coatings
GCP Ground Control Point UPB University of the Philippines
Baguio
GNSS Global Navigation Satellite System
UP-TCAGP University of the Philippines
GPS Global Positioning System – Training Center for Applied
HEC-HMS Hydrologic Engineering Center - Geodesy and Photogrammetry
Hydrologic Modeling System UTM Universal Transverse Mercator
HEC-RAS Hydrologic Engineering Center - WGS World Geodetic System
River Analysis System
HC High Chord
IDW Inverse Distance Weighted
[interpolation method]
IMU Inertial Measurement Unit
kts knots
LAS LiDAR Data Exchange File format
xiii
xiv
CHAPTER 1: OVERVIEW OF THE PROGRAM AND

BAUANG RIVER
Enrico C. Paringit, Dr. Eng., Dr. George Puno, and Eric Bruno
1.1 Background of the Phil-LiDAR 1 Program
The University of the Philippines Training Center for Applied Geodesy and Photogrammetry (UP-TCAGP)
launched a research program in 2014 entitled “Nationwide Hazard Mapping using LiDAR” or Phil-LiDAR 1,
supported by the Department of Science and Technology (DOST) Grants-in-Aid (GiA) Program. The program
was primarily aimed at acquiring a national elevation and resource dataset at a sufficient resolution to
produce information necessary to support the different phases of disaster management. Particularly, it
targeted to operationalize the development of flood hazard models that would produce updated and
detailed flood hazard maps for the major river systems in the country.
The program was also aimed at producing an up-to-date and detailed national elevation dataset suitable
for a 1:5,000 scale mapping, with 50 cm and 20 cm horizontal and vertical accuracies, respectively. These
accuracies were achieved through the use of the state-of-the-art Light Detection and Ranging (LiDAR)
airborne technology procured by the project through the DOST. The methods applied in this report are
thoroughly described in a separate publication entitled “Flood Mapping of Rivers in the Philippines Using
Airborne LiDAR: Methods” (Paringit, et. al., 2017), available separately.
The implementing partner university for the Phil-LiDAR 1 Program is the University of the Philippines Baguio
(UPB). UPB is in charge of processing LiDAR data and conducting data validation reconnaissance, cross
section, bathymetric survey, validation, river flow measurements, flood height and extent data gathering,
flood modeling, and flood map generation for the twelve (12) river basins in the Northern Luzon Region.
The university is located in Baguio City in the province of Benguet.
1.2 Overview of the Bauang River Basin
The Bauang River Basin covers five (5) municipalities and one (1) city in the province of La Union, and four
(4) municipalities and one (1) city in the province of Benguet. According to the Department of Environment
and Natural Resources - River Basin Control Office (DENR-RBCO), it has a drainage area of 353 square
kilometers, and an estimated 863 million cubic meters (MCM) in annual run-off (River Basin Control Office,
2017).
Its main stem, the Bauang River, is part of the river systems in the Northern Luzon Region. The river stream
network runs along the Municipalities of Bauang, Naguilian, and San Fernando City in the province of La
Union.
1
Figure 1. Location map of the Bauang River Basin (in brown)
According to the 2015 national census conducted by the National Statistics Office (NSO), the total of
population of residents within the immediate vicinity of the river is 42,169, distributed among the twenty-
four (24) barangays in the Municipalities of Bauang, Naguilian and San Fernando City (Philippine Statistics
Authority, 2016).
The primary sources of livelihood for the locals are fishing, farming, and quarrying. The main agricultural
products of Bauang include rice, corn, guapple, and native grapes. Bauang is also a tourist destination,
known for its fine sand beaches (Inigo, 2015).
On July 2, 2015, the region was placed under a state of calamity due to floods caused by Typhoon Egay.
2
CHAPTER 2: LIDAR DATA ACQUISITION OF THE

BAUANG FLOODPLAIN
Engr. Louie P. Balicanta, Engr. Christopher Cruz, Lovely Gracia Acuña, Engr. Gerome Hipolito, Jasmine T.
Alviar, and Engr. Brylle Adam G. De Castro
The methods applied in this Chapter were based on the DREAM methods manual (Sarmiento, et al., 2014)
and further enhanced and updated in Paringit, et al. (2017).
2.1 Flight Plans
To initiate the LiDAR acquisition survey of the Bauang floodplain, the Data Acquisition Component (DAC)
created flight plans within the delineated priority area for the floodplain in La Union. These missions
were planned for sixteen (16) lines and ran for at most four and a half (4.5) hours, including take-off,
landing, and turning time. The Pegasus LiDAR System was utilized for the missions (See Annex 1 for the
sensor specifications). The flight planning parameters for the LiDAR system are found in Table 1. Figure 2
illustrates the flight plans and base stations for the Bauang floodplain survey.
Table 1. Flight planning parameters for the Pegasus LiDAR system

Pulse
Block Field of Repetition Scan Fre- Average Average
Flying Height Overlap
Name View Frequen- quency Speed Turn Time
(m AGL) (%)
(θ) cy (PRF) (Hz) (kts) (Minutes)
(kHz)
10C 1200 / 1500 / 1800 30 50 200 30 130 5
10D 1200 / 1500 30 50 200 30 130 5
10F 1200 30 50 200 30 130 5
10G 1500 30 50 200 30 130 5
3
Figure 2. Flight plans and base stations used to cover the Bauang floodplain survey
4
2.2 Ground Base Stations
The field team for this undertaking was able to recover the following NAMRIA control stations: (i.) LUN-62
and (ii.) LUN-176), with second (2nd) order accuracy; (iii.) LUN-3062, (iv.) LUN-3047, and (v.) LUN-3129,
with fourth (4th) order accuracy. The field team also recovered benchmark BM LU-94, which is of first
(1st) order accuracy. The certifications for the NAMRIA reference points are found in Annex 2; while the
baseline processing reports for the established control points are found in Annex 3. These were used as
base stations during the flight operations for the entire duration of the survey, held on February 25 –
March 8, 2015. The base stations were observed using dual frequency GPS receivers, TRIMBLE SPS 852 and
SPS 985. The flight plans and the locations of base stations used during the aerial LiDAR acquisition in the
Bauang floodplain are presented in Figure 2. The composition of the full project team is shown in Annex 4.
Figure 3 to Figure 8 exhibit the recovered NAMRIA reference points within the area. Table 2 to Table 7
provide the details about the NAMRIA control stations and established points. Table 8 lists all of the ground
control points occupied during the acquisition, together with the dates of utilization.
5

(b)
(a)
Figure 3. (a) GPS set-up over LUN-62, as recovered 15 m south from the first access ladder of the river
control and about 100 m north from the end. It is also situated 300 m south of a hanging bridge. (b)
NAMRIA reference point LUN-62, as recovered by the field team
Table 2. Details of the recovered NAMRIA horizontal control point LUN-62, used as a base station for the
LiDAR acquisition
Station Name LUN-62
Order of Accuracy 2nd
Relative Error (horizontal positioning) 1 in 50,000
Latitude 16° 33’ 19.98115” 120° 23’
Geographic Coordinates, Philippine Refer- 28.76004”
Longitude
ence of 1992 Datum (PRS 92)
Ellipsoidal Height 33.18400 m
Grid Coordinates, Philippine Transverse Easting 435034.926 m
Mercator Zone 5 (PTM Zone 5 PRS 92) Northing 1831016.667 m
Geographic Coordinates, World Geodetic Latitude 16° 33’ 14.07106”
System 1984 Datum Longitude 120° 23’ 33.49149” 69.44500
(WGS 84) Ellipsoidal Height m
Grid Coordinates, Universal Transverse Mer-
Easting 221592.72 m
cator Zone 51 North
Northing 1832084.35 m
(UTM 51N WGS 1984)
6
(b)
(a)
Figure 4. (a) GPS set-up over LUN-176, recovered near a corner of a farm dike, about 15 m southeast
of the well, and about 20 m southwest of the nearest house. (b) NAMRIA reference point LUN-176, as
recovered by the field team
Table 3. Details of the recovered NAMRIA horizontal control point LUN-176, used as a base station for the
LiDAR acquisition
Station Name LUN-176
Latitude 16° 46’ 14.35394” 120° 24’
Longitude
Grid Coordinates, Philippine Transverse Mer- Easting 436193.115 m
cator Zone 3 (PTM Zone 5 PRS 92) Northing 1854816.574 m
Latitude 16° 46’ 8.39718”
Geographic Coordinates, World Geodetic
Longitude 120° 24’ 10.13252”
System 1984 Datum (WGS 84)
Easting 222990.04 m
cator Zone 51 North
(UTM 51N WGS 1984)
7
(b)
(a)
Figure 5. (a) GPS set-up over LUN-3062, recovered at the top of a dike, approximately 100 m north of
Philippine Central College of Arts Sciences & Technology, and 80 m north of Naguilian emission testing
center in Barangay Natividad, Naguilian, La Union. (b) NAMRIA reference point LUN-3062, as recovered
by the field team
Table 4. Details of the recovered NAMRIA horizontal control point LUN-3062, used as a base station for
the LiDAR acquisition
Order of Accuracy 4th
Latitude 16° 31’ 55.00993”
Geographic Coordinates, Philippine Refer-
Longitude 120° 23’ 12.50504”
Geographic Coordinates, World Geodetic Latitude 16° 31’ 49.10470” 120° 23’
System 1984 Datum Longitude 17.23850”
(WGS 84) Ellipsoidal Height 61.64400 m
Easting 221076.59 m
cator Zone 51 North
(UTM 51N WGS 1984)
8
Figure 6. LUN-3047, as recovered in Barangay Nazareno, Agoo, La Union
Table 5. Details of the established horizontal control point LUN-3047, used as a base station for the LiDAR
acquisition
Latitude 16° 20’ 55.96430” North
Longitude 120° 21’ 47.08672” East
Ellipsoidal Height 43.62100 meters
Grid Coordinates, Philippine Transverse Mer- Easting 431948.446 meters
cator Zone 3 (PTM Zone 5 PRS 92) Northing 1808156.256 meters
Latitude 16°20’ 50.09786” North
Geographic Coordinates, World Geodetic
Longitude 120° 21’ 51.83567” East
System 1984 Datum (WGS 84)
Ellipsoidal Height 80.44800 meters
Easting 218278.33 meters
cator Zone 51 North
Northing 1809242.68 meters
(UTM 51N PRS 1992)
9
(b)
(a)
Figure 7. (a) GPS set-up over LUN-3129, recovered beside the National Road about 50 meters northeast
of the nearest house. (b) NAMRIA reference point LUN-3129, as recovered by the field team
Table 6. Details of the recovered NAMRIA horizontal control point LUN-3129, used as a base station for
the LiDAR acquisition
Latitude 16° 31’ 55.00993” 120° 23’
Longitude
Geographic Coordinates, World Geodetic Latitude 16° 31’ 49.10470” 120° 23’
System 1984 Datum Longitude 17.23850”
Easting 221076.59 m
cator Zone 51 North
(UTM 51N WGS 1984)
10
Figure 8. NAMRIA reference point LU-94, as recovered by the field team
Table 7. Details of the recovered NAMRIA vertical control point LU-94, used as a base station for the
LiDAR acquisition
Station Name LU-94
Relative Error (horizontal positioning) 1:50,000
Latitude 16°42’ 38.64674”
Longitude 120°20’35.05091”
Geographic Coordinates, World Geodetic Latitude 16°42’ 38.64674”
System 1984 Datum Longitude 120°20’35.05091”
11
Table 8. Ground control points used during the LiDAR data acquisition
Date Surveyed Flight Number Mission Name Ground Control Points
February 25, 2014 1153P 1BLK10AS056B LUN-176 and LUN-3129
February 26, 2014 1155P 1BLK10C057A LUN-176 and LU-94
February 27, 2014 1159P 1BLK10GD058A LUN-62 and LUN-302
February 27, 2014 1161P 1BLK10D058B LUN-62 & LUN-3062
February 28, 2014 1163P 1BLK10F059A LUN-72 & LUN-3047
March 2, 2014 1171P 1BLK10CDS061A LUN-62 and LUN-3062
March 2, 2014 1173P 1BLK10DS061B LUN-62 and LUN-3062
March 3, 2014 1177P 1BLK10CS062B LUN-176 and LUN-3129
March 8, 2014 1197P 1BLK10GCS067B LUN-62 and LUN-3062
2.3 Flight Missions
A total of nine (9) flight missions were conducted to complete the LiDAR data acquisition in the Bauang
floodplain, for a total of twenty-six hours and seven minutes (26+7) of flying time for RP-C9022. All missions
were acquired using the Pegasus LiDAR system. The flight logs for the missions are provided in Annex 6.
Table 9 indicates the total area of actual coverage and the corresponding flying hours per mission; while
Table 10 presents the actual parameters used during the LiDAR data acquisition.
12
Table 9. Flight missions for the LiDAR data acquisition in the Bauang floodplain
Flying
Area Area Sur- Hours
Flight Surveyed veyed No. of
Date Sur- Flight Surveyed
Plan Area within the outside the Images
veyed Number Area (km2)
(km2) Floodplain Floodplain (Frames)
Min
Hr
(km2) (km2)
February
25, 2014 1153P 223.74 0.000 0.000 92.675 NA 1 45
February
26, 2014 1155P 298.492 18.526 18.526 233.071 574 3 31
February
27, 2014 1159P 389.031 67.904 67.904 214.796 75 3 31
February
1161P 253.115 33.684 33.684 152.408 369 2 43
27, 2014
February
1163P 275.809 48.444 48.444 292.968 819 3 37
28, 2014
March 2,
2014 1171P 551.607 30.883 30.883 161.940 334 3 15
March 2,
2014 1173P 253.115 56.898 56.898 117.251 334 2 13
March 3,
2014 1177P 298.492 47.130 47.130 228.987 501 2 43
March 8,
2014 1197P 434.408 5.468 5.468 103.405 257 2 49
TOTAL 2977.809 1906.437 308.936 1597.502 3263 26 7
13
Table 10. Actual parameters used during the LiDAR data acquisition
Flying Scan Fre- Average Average

Flight Overlap PRF
Height FOV (θ) quency Speed Turn Time
Number (%) (Hz)
(m AGL) (kHz) (kts) (Minutes)
1153P 1200 30 50 200 30 130 5
1155P 1200 30 50 200 30 130 5
1159P 1200 30 50 200 30 130 5
1161P 1200 50 30 200 50 130 5
1163P 1200 50 30 200 50 130 5
1171P 1500 30 50 200 30 130 5
1173P 1200 30 50 200 30 130 5
1177P 1800 30 50 200 30 130 5
1197P 1500 30 50 200 30 130 5
2.4 Survey Coverage
The Bauang floodplain is located in the province of La Union. The municipalities of Bauang, San Fernando
City, San Juan, Naguilian, Pugo, Santol, Caba, Burgos, and Tubao in La Union were mostly covered by
the survey. The municipalities and cities surveyed with at least one (1) square kilometer coverage are
enumerated in Table 11. The actual coverage of the LiDAR acquisition for the Bauang floodplain is presented
in Figure 9. See Annex 7 for the flight status reports.
14
Table 11. List of municipalities and cities surveyed during the Bauang floodplain LiDAR survey
Area of Municipality Total Area Sur- Percentage of Area
Province Municipality/City
(km2) veyed (km2) Surveyed (%)
Bauang 85.260 85.260 100%
San Fernando City 121.053 121.033 100%
San Juan 53.442 53.308 100%
Naguilian 86.388 85.801 99%
Pugo 60.537 56.057 93%
Santol 97.972 85.920 88%
Caba 56.191 46.611 83%
Burgos 51.920 38.164 74%
Tubao 53.873 37.759 70%
Rosario 64.328 44.480 69%
La Union
San Gabriel 154.189 88.684 58%
Aringay 95.648 48.624 51%
Bacnotan 80.671 38.888 48%
Bagulin 77.967 33.921 44%
Sudipen 75.745 28.585 38%
Balaoan 60.960 15.085 25%
Bangar 45.120 10.915 24%
Agoo 33.707 5.938 18%
Luna 50.658 4.962 10%
Santo Tomas 58.531 3.714 6%
Sugpon 180.280 49.474 27%
Tagudin 54.348 11.711 22%
Ilocos Sur
Alilem 132.176 17.845 14%
Suyo 148.521 16.742 11%
Pangasinan Sison 151.962 50.992 34%
Sablan 90.217 21.162 23%
Benguet
Tuba 322.020 7.266 2%
Total 2543.68 1108.9 43.59%
15
Figure 9. Actual LiDAR survey coverage for the Bauang floodplain
16
CHAPTER 3: LIDAR DATA PROCESSING OF THE BAUANG

FLOODPLAIN
Engr. Ma. Rosario Concepcion O. Ang, Engr. John Louie D. Fabila, Engr. Sarah Jane D. Samalburo, Engr.
Harmond F. Santos, Jovy Anne S. Narisma, Engr. Ma. Ailyn L. Olanda, Engr. Merven Mattew D. Natino,
Engr. Kenneth A. Solidum, Engr. Jommer M. Medina, and Carl Joshua S. Lacsina
The methods applied in this Chapter were based on the DREAM methods manual (Ang, et al., 2014) and
further enhanced and updated in Paringit, et al. (2017).
3.1 Overview of the LIDAR Data Pre-Processing
The data transmitted by the DAC were checked for completeness based on the list of raw files required
to proceed with the pre-processing of the LiDAR data. Upon acceptance of the LiDAR field data,
georeferencing of the flight trajectory was done in order to obtain the exact location of the LiDAR sensor
when the laser was shot. Point cloud georectification was performed to incorporate the correct position
and orientation for each point acquired. The georectified LiDAR point clouds were subjected to quality
checking to ensure that the required accuracies of the program, which are the minimum point density,
and the vertical and horizontal accuracies, were met. The point clouds were then categorized into various
classes before generating Digital Elevation Models (DEMs), such as the Digital Terrain Model (DTM) and
the Digital Surface Model (DSM).
Using the elevation of points gathered from the field, the LiDAR-derived digital models were calibrated.
Portions of the river that were barely penetrated by the LiDAR system were replaced by the actual river
geometry, measured from the field by the Data Validation and Bathymetry Component (DVBC). LiDAR
acquired temporally were then mosaicked to completely cover the target river systems in the Philippines.
Orthorectification of images acquired simultaneously with the LiDAR data was accomplished through the
help of the georectified point clouds, and the metadata containing the time the image was captured.
These processes are summarized in the diagram in Figure 10.
17
Figure 10. Schematic diagram for the Data Pre-Processing Component
3.2 Transmittal of Acquired LiDAR Data
The data transfer sheets for all the LiDAR missions for the Bauang floodplain can be found in Annex 5.
Missions flown over Bauang, La Union during the first survey conducted in February 2014 and the second
survey conducted in March 2014 both used the Airborne LiDAR Terrain Mapper (ALTM™ Optech Inc.)
Pegasus system. The DAC transferred a total of 142.68 Gigabytes of Range data, 1.23 Gigabytes of POS
data, 121.72 Megabytes of GPS base station data, and 160 Gigabytes of raw image data to the data server
on February 25, 2014 for the first survey, and on March 1, 2014 for the second survey. The Data Pre-
processing Component (DPPC) verified the completeness of the transferred data. The whole dataset for
the Bauang survey was fully transferred on May, 17 2014, as indicated on the data transfer sheets for the
Bauang floodplain.
3.3 Trajectory Computation
The Smoothed Performance Metric parameters of the computed trajectory for Flight 1153P, one of the
Bauang flights, which are the North, East, and Down position RMSE values, are illustrated in Figure 11. The
x-axis corresponds to the time of flight, which is measured by the number of seconds from the midnight of
the start of the GPS week, which fell on February 25, 2014 at 00:00 hrs. on that week. The y-axis represents
the RMSE value for that particular position.
18
Figure 11. Smoothed Performance Metric Parameters of Bauang Flight 1153P
The time of flight was from 266,000 seconds to 275,000 seconds, which corresponds to the afternoon of
February, 2014. The initial spike reflected on the data corresponds to the time that the aircraft was getting
into position to start the acquisition, and the POS system was starting to compute for the position and
orientation of the aircraft. Redundant measurements from the POS system quickly minimized the RMSE
values of the positions. The periodic increase in RMSE values from an otherwise smoothly curving set
of RMSE values indicates the turn-around period of the aircraft, when the aircraft makes a turn to start
a new flight line. Figure 11 demonstrates that the North position RMSE peaked at 1.90 centimeters, the
East position RMSE peaked at 2.20 centimeters, and the Down position RMSE peaked at 3.60 centimeters,
which are within the prescribed accuracies described in the methodology.
19
Figure 12. Solution Status Parameters of Bauang Flight 1153P
The Solution Status parameters of Flight 1153P, one of the Bauang flights, which are the number of GPS
satellites, Positional Dilution of Precision (PDOP), and the GPS processing mode used, are depicted in Figure
12. The graphs indicate that the number of satellites during the acquisition did go down to 6. Majority of
the time, the number of satellites tracked was between 5 and 8. The PDOP value did not go above the
value of 3, which indicates optimal GPS geometry. The processing mode remained at the value of 0 for
majority of the survey, with some peaks to up to 1, attributed to the turns performed by the aircraft.
The value of 0 signifies a Fixed, Narrow-Lane mode, which is the optimum carrier-cycle integer ambiguity
resolution technique available for POSPAC MMS. All of the parameters satisfied the accuracy requirements
for optimal trajectory solutions, as indicated in the methodology. The computed best estimated trajectory
for all Bauang flights is exhibited in Figure 13.
20
Figure 13. The best estimated trajectory conducted over the Bauang floodplain
3.4 LiDAR Point Cloud Computation
The produced LAS data contains two hundred and sixty (260) flight lines, with each flight line containing
two (2) channels, since the Pegasus system contains two (2) channels. The summary of the self-calibration
results for all flights over the Bauang floodplain, obtained through LiDAR processing in the LiDAR Mapping
Suite (LMS) software, is given in in Table 12.
Table 12. Self-calibration results values for the Bauang flights

Parameter Absolute Value Computed Value
Boresight Correction stdev (<0.001degrees) 0.000370
IMU Attitude Correction Roll and Pitch Corrections stdev (<0.001degrees) 0.001263
GPS Position Z-correction stdev (<0.01meters) 0.0076
Optimum accuracy was obtained for all Bauang flights, based on the computed standard deviations of the
corrections of the orientation parameters. The standard deviation values for individual blocks are available
in Annex 8: Mission Summary Reports.
3.5 LiDAR Data Quality Checking
The boundaries of the processed LiDAR data on top of a SAR Elevation Data over the Bauang Floodplain are
represented in Figure 14. The map shows gaps in the LiDAR coverage that are attributed to cloud coverage.
21
Figure 14. Boundaries of the processed LiDAR data over the Bauang floodplain
The total area covered by the Bauang missions is 1171.69 square kilometers, comprised of eight (8) flight
acquisitions that are grouped and merged into four (4) blocks, as outlined in Table 13.

22
Table 13. List of LiDAR blocks for the Bauang floodplain

LiDAR Blocks Flight Numbers Area (sq. km)
1153P
1169P
LaUnion_Blk10C 399.93
1175P
1197P
1157P
LaUnion_Blk10D 1159P 326.71
1171P
LaUnion_Blk10F 1161P 324.81
1157P
LaUnion_Blk10G 120.23
1197P
TOTAL 1171.69 sq.km
The overlap data for the merged LiDAR blocks, showing the number of channels that pass through a
particular location, is presented in Figure 15. Since the Pegasus system employs two (2) channels, it is
expected to have an average value of 2 (blue) for areas where there is limited overlap, and a value of 3
(yellow) or more (red) for areas with three or more overlapping flight lines.
23
Figure 15. Image of data overlap for Bauang floodplain
The overlap statistics per block for the Bauang floodplain can be found in Annex 8. One (1) pixel corresponds
to 25.0 square meters on the ground. For this area, the minimum and maximum percent overlaps were
35.05% and 68.41%, respectively, which passed the 25% requirement.
The pulse density map for the merged LiDAR data, with the red parts showing the portions of the data that
satisfy the two (2) points per square meter criterion, is exhibited in Figure 16. It was determined that all
LiDAR data for the Bauang floodplain satisfy the point density requirement, and that the average density
for the entire survey area is 3.12 points per square meter.
24
Figure 16. Density map of merged LiDAR data for the Bauang floodplain
The elevation difference between overlaps of adjacent flight lines is shown in Figure 17. The default color
range is from blue to red. Bright blue areas correspond to portions where elevations of a previous flight
line, identified by its acquisition time, are higher by more than 0.20 meters relative to the elevations of its
adjacent flight line. Bright red areas indicate portions where elevations of a previous flight line are lower
by more than 0.20 meters relative to the elevations of its adjacent flight line. Areas with bright red or
bright blue colors were investigated further using the Quick Terrain (QT) Modeler software.
25
Figure 17. Elevation difference map between flight lines for the Bauang floodplain
A screen capture of the processed LAS data from Bauang Flight 1153P loaded in QT Modeler is provided
in Figure 18. The upper left image shows the elevations of the points from two (2) overlapping flight strips
traversed by the profile, illustrated by a dashed red line. The x-axis represents the length of the profile. It
is evident that there were differences in elevation, but the differences did not exceed the 20-centimeter
mark. This profiling was repeated until the quality of the LiDAR data became satisfactory. No reprocessing
was done for this LiDAR dataset.
26
Figure 18. Quality checking for Bauang Flight 1153P, using the Profile Tool of QT Modeler
3.6 LiDAR Point Cloud Classification and Rasterization
Table 14. Bauang classification results in TerraScan

Pertinent Class Total Number of Points
Ground 178,475,094
Low Vegetation 187,594,122
Medium Vegetation 317,240,334
High Vegetation 638,889,274
Building 23,295,587
The tile system that TerraScan employed for the LiDAR data, as well as the final classification image for
a block in the Bauang floodplain, are presented in Figure 19. A total of 1,477 1km by 1km tiles were
produced. The number of points classified according to the pertinent categories is illustrated in Table 14.
The point cloud had a maximum and minimum height of 1038.37 meters and 45.5 meters, respectively.
27
Figure 19. (a) Tiles for Bauang floodplain, and (b) classification results in TerraScan
An isometric view of an area before and after running the classification routines is shown in Figure 20. The
ground points are in orange, the vegetation is in different shades of green, and the buildings are in cyan.
It is visible that the residential structures adjacent or even below canopy were classified correctly, due to
the density of the LiDAR data.
28
Figure 20. Point cloud (a) before and (b) after classification
The production of last return (V_ASCII) and the secondary (T_ ASCII) DTM, and the first (S_ ASCII) and last
(D_ ASCII) return DSM of the area are illustrated in Figure 21, in top view display. The images convey that
the DTMs are a representation of the bare earth; while the DSMs reflect all features that are present, such
as buildings and vegetation.
29
Figure 21. The production of (a) last return DSM and (b) DTM, (c) first return DSM and (d) secondary DTM
in some portion of Bauang floodplain
3.7 LiDAR Image Processing and Orthophotograph Rectification
The 1,278 1km by 1km tiles area covered by the Bauang floodplain is displayed in Figure 22. After
employing tie point selection to fix photo misalignments, color points were added to smoothen out visual
inconsistencies along the seamlines where photos overlap. The Bauang floodplain survey attained a
total of 971.99 square kilometers in orthophotographic coverage, comprised of 2,345 images. Zoomed-in
versions of sample orthophotographs, identified by their tile numbers, are provided in Figure 23.
30
Figure 22. The Bauang floodplain, with available orthophotographs
Figure 23. Sample orthophotograph tiles for the Bauang floodplain
3.8 DEM Editing and Hydro-Correction
Four (4) mission blocks were processed for the Bauang floodplain. These blocks are composed of La Union
blocks, with a total area of 1173.92 square kilometers. Table 15 summarizes the names and corresponding
areas of the blocks, in square kilometers.
31
Table 15. LiDAR blocks with their corresponding areas

LiDAR Blocks Area (sq.km)
LaUnion_Blk10C 402.25
LaUnion_Blk10D 326.70
LaUnion_Blk10F 325.22
LaUnion_Blk10G 119.75
TOTAL 1173.92 sq.km
Portions of the DTM before and after manual editing are exhibited in Figure 24. The bridge (Figure 24a) was
considered to be an obstruction to the flow of water along the river, and had to be removed (Figure 24b)
in order to hydrologically correct the river. The embankment (Figure 24c) was misclassified and removed
during the classification process, and had to be retrieved to complete the surface (Figure 24d) to allow for
the correct flow of water. Another case was a building that was still present in the DTM after classification
(Figure 24e), and had to be removed through manual editing (Figure 24f).
32
Figure 24. Portions in the DTM of the Bauang floodplain – a bridge (a) before and (b) after manual
editing; an embankment (c) before and (d) after data retrieval; and a building (e) before and (f) after
manual editing
33
3.9 Mosaicking of Blocks
The LaUnion_Blk10A block was used as the reference block at the start of mosaicking, because it was
referred to a base station with an acceptable order of accuracy. Table 16 specifies the shift values applied
to each LiDAR block during the mosaicking process.

The mosaicked LiDAR DTM for the Bauang floodplain is displayed in Figure 25. It is demonstrated that the
entire Bauang floodplain is 99.69% covered by LiDAR data.
Table 16. Shift values of each LiDAR block of the Bauang floodplain
Shift Values (meters)
Mission Blocks
x y z
LaUnion_Blk10C 0.00 0.00 0.00
LaUnion_Blk10D 0.00 0.00 0.00
LaUnion_Blk10F 0.00 0.00 0.00
LaUnion_Blk10G 0.00 0.00 0.00
34
Figure 25. Map of processed LiDAR data for the Bauang floodplain
.
35
3.10 Calibration and Validation of Mosaicked LiDAR DEM
To undertake the data validation of the mosaicked LiDAR DEMs, the DVBC conducted a validation survey
along the Bauang floodplain. The extent of the validation survey done in Bauang to collect points with which
the LiDAR dataset was validated is illustrated in Figure 26, with the validation survey points highlighted in
green. A total of 4,794 survey points were used for the calibration and validation of the Bauang LiDAR data.
Random selection of 80% of the survey points resulted in 3,835 points, which were used for calibration.
A good correlation between the uncalibrated mosaicked LiDAR elevation values and the ground survey
elevation values is reflected in Figure 27. Statistical values were computed from the extracted LiDAR values
using the selected points to assess the quality of the data and to obtain the values for vertical adjustment.
The computed height difference between the LiDAR DTM and the calibration elevation values is 3.64
meters, with a standard deviation of 0.17 meters. Calibration of the Bauang LiDAR data was performed by
subtracting the height difference value, 3.64 meters, from the Bauang mosaicked LiDAR data. Table 17 lists
the statistical measurements of the compared elevation values between the LiDAR data and the calibration
data.
36
Figure 26. Map of the Bauang floodplain, with the validation survey points in green

37
Figure 27. Correlation plot between the calibration survey points and the LiDAR data
Table 17. Calibration statistical measures
Calibration Statistical Measures Value (meters)

Height Difference 3.64
Standard Deviation 0.17
Average -3.64
Minimum -3.98
Maximum -3.29
The remaining 20% of the total survey points, resulting in 212 points, were used for the validation of
the calibrated Bauang DTM. A good correlation between the calibrated mosaicked LiDAR elevation values
and the ground survey elevation, which reflects the quality of the LiDAR DTM, is demonstrated in Figure
28. The computed RMSE between the calibrated LiDAR DTM and the validation elevation values is 0.14
meters, with a standard deviation of 0.12 meters, as indicated in Table 18.
38
Figure 28. Correlation plot between the validation survey points and the LiDAR data
Table 18. Validation statistical measures
Validation Statistical Measures Value (meters)

RMSE 0.14
Standard Deviation 0.12
Average 0.07
Minimum -0.17
Maximum 0.38
39
3.11 Integration of Bathymetric Data into the LiDAR Digital Terrain Model
For bathy integration, centerline and cross-section data were available for Bauang, with 4,351 bathymetric
survey points. The resulting raster surface produced was obtained through the Local Polynomial interpolation
method. After burning the bathymetric data to the calibrated DTM, assessment of the interpolated
surface is represented by the computed RMSE value of 0.50 meters. The extent of the bathymetric survey
conducted by the DVBC in Bauang, integrated with the processed LiDAR DEM, is illustrated in Figure 29.
40
Figure 29. Map of the Bauang floodplain, with bathymetric survey points shown in blue
41
3.12 Feature Extraction
The features salient in flood hazard exposure analysis include buildings, road networks, bridges, and water
bodies within the floodplain area, with a 200-meter buffer zone. Mosaicked LiDAR DEM with a 1-meter
resolution was used to delineate footprints of building features, which consist of residential buildings,
government offices, medical facilities, religious institutions, and commercial establishments, among
others. Road networks – comprised of main thoroughfares, such as highways, and municipal and barangay
roads – are essential for routing disaster response efforts. These features are represented by a network of
road centerlines.
3.12.1 Quality Checking of Digitized Features’ Boundary

The Bauang floodplain, including its 200-meter buffer zone, has a total area of 153.56 square kilometers.
Of this area, a total of 5.0 square kilometers, corresponding to a total of 913 building features, were
considered for quality checking (QC). Figure 30 displays the QC blocks for the Bauang floodplain.
Figure 30. Blocks (in blue) of Bauang building features that were subjected to QC
Quality checking of the Bauang building features resulted in the ratings given in Table 19.
Table 19. Quality checking ratings for the Bauang building features
FLOODPLAIN COMPLETENESS CORRECTNESS QUALITY REMARKS
Bauang 99.87 100.00 98.50 PASSED
3.12.2 Height Extraction

Height extraction was done for 29,837 building features in the Bauang floodplain. Of these building
features, 1,845 building features were filtered out after height extraction, resulting in 27,992 buildings
with height attributes. The lowest building height is at 2.00 meters, while the highest building is at 9.52
meters.
42
3.12.3 Feature Attribution

Data collected from various sources, including OpenStreetMap, Google Maps, and Google Earth, were
used in the attribution of building features. Areas without available data were subjected to field attribution
using ESRI’s Collector App. The application can be accessed offline, and the data collected can be synced
with ArcGIS Online when WiFi or mobile data is available.
Table 20 summarizes the number of building features per type. Table 21 indicates the total length of each
road type, and Table 22 specifies the number of water features extracted per type.
Table 20. Building features extracted for the Bauang floodplain

Facility Type No. of Features
Residential 26,955
School 408
Market 57
Agricultural/Agro-Industrial Facilities 27
Medical Institutions 13
Barangay Hall 17
Military Institution 0
Sports Center/Gymnasium/Covered Court 11
Telecommunication Facilities 9
Transport Terminal 0
Warehouse 3
Power Plant/Substation 53
NGO/CSO Offices 17
Police Station 0
Water Supply/Sewerage 2
Religious Institutions 65
Bank 3
Factory 15
Gas Station 17
Fire Station 2
Other Government Offices 21
Other Commercial Establishments 296
Total 27,991
Table 21. Total length of extracted roads for the Bauang floodplain
Road Network Length (km)
Floodplain Barangay City/Municipal Provincial National Total
Others
Road Road Road Road
Bauang 90.56 72.41 32.78 6.33 0.00 205.59
Table 22. Number of extracted water bodies for the Bauang floodplain
Water Body Type
Floodplain Total
Rivers/Streams Lakes/Ponds Sea Dam Fish Pen
Bauang 113 41 0 0 0 154
A total of sixty (60) bridges and culverts over small channels that are part of the river network were also
extracted for the floodplain.
43
3.12.4 Final Quality Checking of Extracted Features

All extracted ground features were completely given the required attributes. All these output features
comprise the flood hazard exposure database for the floodplain. This completes the feature extraction
phase of the project.
Figure 31 represents the Digital Surface Model (DSM) of the Bauang floodplain, overlaid with its ground
features.
Figure 31. Extracted features for the Bauang floodplain
44
CHAPTER 4: LIDAR VALIDATION SURVEY AND

MEASUREMENTS OF THE BAUANG RIVER BASIN
Engr. Louie P. Balicanta, Engr. Joemarie S. Caballero, Ms. Patrizcia Mae. P. dela Cruz, Engr. Dexter T.
Lozano, For. Dona Rina Patricia C. Tajora, Elaine Bennet Salvador, For. Rodel C. Alberto, Cybil Claire
Atacador, and Engr. Lorenz R. Taguse
The methods applied in this Chapter were based on the DREAM methods manual (Balicanta, et al., 2014)
4.1 Summary of Activities
The DVBC, together with the UPB Phil-LiDAR 1 team, conducted field surveys in three (3) river basins in
Northern Luzon, including the Bauang River Basin, on July 3-17, 2015, and on November 26, 2015. The
scope of work was comprised of: (i.) initial reconnaissance; (ii.) control point survey; (iii.) cross-section and
bridge as-built surveys of the Guesset and Mamat-ing Bridges; (iv.) validation point acquisition of about
50 kilometers covering the Bauang River Basin; and (v.) bathymetric survey from Barangay Cabaritan,
Municipality of Naguilian down to the mouth of the river in Barangay Pugo, Municipality of Bauang, with
an estimated length of 11.189 kilometers, using Trimble® SPS 882 GNSS in PPK and RTK survey technique
and an OHMEX™ single beam echo sounder. See Figure 32 for the illustration of the extent of the surveys.
45
Figure 32. Extent of the bathymetric survey (in blue line) in Bauang River and the LiDAR data validation
survey (in red)
46
4.2 Control Survey
The GNSS network used for the Bauang River Basin is composed of five (5) loops established on July 4, 11,
and 16, 2015, occupying the following reference points: (i.) LU-67, a first-order BM in Barangay Nagrebcan;
and (ii.) LUN-61, a second-order GCP in Barangay Guerrero. Both are located in the Municipality of Bauang.
Four (4) control points were established along the approach of bridges, namely: (i.) UP-AGY, at the Arigay
Bridge in Barangay Poblacion, Municipality of Aringay; (ii.) UP-CBN, at the Cabaroan Bridge in Barangay
Cabaroan, Municipality of San Juan; (iii.) UP-GST, at the Guesset Bridge in Barangay Ambaracao, Municipality
of Naguilian; (iv.) UP-MMG, at the Mamat-ing Bridge in Barangay Mamat-ing Norte, Municipality of
Naguilian; and (v.) UP-NG, at the Naguilian Bridge in Barangay Cabaritan Sur, Municipality of Naguilian.
The summary of reference and control points and their corresponding locations is given in Table 23; while
the established GNSS network is illustrated in Figure 33.
Table 23. List of reference and control points occupied in the Bauang River Survey (Source: NAMRIA; UP-
TCAGP)
Con- Geographic Coordinates (WGS 84)

Order of
trol
Accuracy Ellipsoid Elevation in Date Estab-
Point Latitude Longitude
Height (m) MSL (m) lished
LUN- 2nd order,
16°29’30.73239” 120°19’46.68075” 49.954 9.213 2007
61 GCP
1st order,
LU-67 - - 49.515 - 2007
BM
UP- UP Estab-
- - - - 7-4-2015
AGY lished
UP- UP Estab-
- - - - 7-11-2015
CBN lished
UP- UP Estab-
- - - - 7-4-2015
GST lished
UP- UP Estab-
- - - - 7-11-2015
NG lished
47
Figure 33. GNSS network of the Bauang River field survey
The GNSS set-ups of reference points and established control points in the Bauang survey are exhibited in
Figure 34 to Figure 39.
48
Figure 34. GNSS receiver Trimble® SPS 985 at LU-67, located at the Bauang Bridge approach in Barangay
Nagrebcan, Municipality of Bauang, La Union
Figure 35. GNSS receiver Trimble® SPS 985 at LUN-61, located near the BPPC Housing in Barangay
Guerrero, Municipality of Bauang, La Union
49
Figure 36. GNSS receiver Trimble® SPS 882 at UP-AGY, located at the Aringay Bridge approach in Barangay
Poblacion, Municipality of Aringay, La Union
Figure 37. GNSS receiver Trimble® SPS 882 at UP-CBN, located at the Cabaroan Bridge approach in
Barangay Cabaroan, Municipality of San Juan, La Union
50
Figure 38. GNSS receiver Trimble® SPS 985 at UP-GST, located at the Guesset Bridge approach in
Barangay Ambaracao, Municipality of Naguilian, La Union
Figure 39. GNSS receiver Trimble® SPS 985 at UP-NG, located at the Naguilian Bridge approach in
Barangay Cabaritan, Municipality of Naguilian, La Union
51
4.3 Baseline Processing
The GNSS baselines were processed simultaneously in TBC by observing that all baselines have fixed
solutions with horizontal and vertical precisions within the +/- 20-centimeter and +/- 10-centimeter
requirement, respectively. In cases where one or more baselines did not meet all of these criteria, masking
was performed. Masking is the removal of portions of baseline data using the same processing software. It
is repeatedly processed until all baseline requirements are met. If the reiteration yields out of the required
accuracy, a re-survey is initiated. The baseline processing results of the control points in the Bauang River
Basin, generated by the TBC software, are summarized in Table 24..
52
Table 24. Baseline processing report for the Bauang River survey
Date of Ob- Solution H. Prec. V. Prec. Geodetic Ellipsoid Dist. ΔHeight
Observation
servation Type (Meter) (Meter) Az. (Meter) (Meter)
UPAGY --- LUN61 7-4-2015 Fixed 0.004 0.027 166°05’07” 11694.520 0.436
LUN61 --- UPGST 7-4-2015 Fixed 0.003 0.018 82°19’22” 8640.260 13.306
LU67 --- LUN61 7-4-2015 Fixed 0.002 0.010 1°33’38” 2754.678 -0.469
UPAGY --- UPGST 7-4-2015 Fixed 0.005 0.025 24°42’08” 13764.436 12.818
UPAGY --- UPNG 7-11-2015 Fixed 0.004 0.019 12°25’55” 15950.207 7.998
UPAGY --- UPNG 7-11-2015 Fixed 0.007 0.030 12°25’56” 15950.202 8.024
UPNG --- LU67 7-11-2015 Fixed 0.002 0.011 76°34’58” 6342.089 8.805
UPCBN --- LU67 7-11-2015 Fixed 0.003 0.014 21°30’22” 17868.131 3.769
UPNG --- UPGST 7-11-2015 Fixed 0.003 0.013 142°57’42” 3848.102 4.774
UPNG --- LUN61 7-11-2015 Fixed 0.003 0.013 235°55’50” 7539.286 -8.335
UPNG --- LUN61 7-11-2015 Fixed 0.003 0.019 235°55’50” 7539.295 -8.524
UPNG --- UPCBN 7-11-2015 Fixed 0.004 0.017 1°27’31” 15157.305 -5.053
UPCBN --- LUN61 7-11-2015 Fixed 0.004 0.018 198°53’37” 20479.168 -3.321
A total of thirteen (13) baselines were processed for the static survey. All of the baselines satisfied the
required accuracy of 0.20 meters for horizontal, and 0.10 meters for vertical, for the GNSS surveys of the
project.
4.4 Network Adjustment
After the baseline processing procedure, network adjustment was performed using TBC. Looking at the
adjusted grid coordinates table of the TBC-generated Network Adjustment Report, it is observed that
the square root of the sum of the squares of x and y must be less than 20 centimeters, and z less than 10
centimeters, or in equation form:
<20cm and
Where:
xe is the Easting Error,
ye is the Northing Error, and
ze is the Elevation Error
The six (6) control points – LUN-61, LU-67, UP-AGY, UP-CBN, UP-GST, and UP-NG – were occupied and
observed simultaneously to form a GNSS loop. The coordinates of LUN-61 and the elevation values of LU-
67 were held fixed during the processing of the control points, as presented in Table 25. Through these
reference points, the coordinates and elevation of the unknown control points were computed.
53
Table 25. Constraints applied to the adjustments of the control points

East σ North σ Height σ Elevation σ
Point ID Type
(Meter) (Meter) (Meter) (Meter)
LU67 Grid Fixed
LUN61 Global Fixed Fixed
Fixed = 0.000001(Meter)
The list of adjusted grid coordinates; i.e., Northing, Easting, Elevation, and computed standard errors of
the control points in the network, is outlined in Table 26. All fixed control points did not yield values for
grid and elevation errors.
Table 26. Adjusted grid coordinates for the control points used in the Bauang floodplain survey
East- North- Eleva- Eleva-
Easting Northing
Point ID ing Error ing Error tion tion Error Constraint
(Meter) (Meter)
(Meter) (Meter) (Meter) (Meter)
LU67 215027.120 0.008 1827985.553 0.007 8.867 ? e
LUN61 214915.551 ? 1825231.473 ? 9.213 0.047 LL
UPAGY 217579.313 0.012 1813837.247 0.010 9.197 0.089
UPCBN 221802.212 0.011 1844530.550 0.009 12.830 0.064
The network was fixed at the reference points. The list of adjusted grid coordinates of the network is
shown in Table 26.Using the equation for horizontal accuracy, and for vertical accuracy, following are the
computations for accuracy:
54
a. LUN-61
Horizontal Accuracy = Fixed
Vertical Accuracy = 4.7 cm < 10 cm
b. LU-67
Horizontal Accuracy = √ ((0.8) ² + (0.7) ²
= √(0.64+ 0.49)
= 1.06 cm < 20 cm
Vertical Accuracy = Fixed
c. UP-AGY
= √(1.44 + 1.0)
= 1.56 cm < 20 cm
d. UP-CBN
= √(1.21 + 0.81)
= 2.02 cm < 20 cm
e. UP-GSN
= √(1.0 + 0.81)
= 1.34 cm < 20 cm
f. UP-NG
= √(0.49 + 0.49)
= 0.98 cm < 20 cm
Following the given formula, the horizontal and vertical accuracy results of the four (4) occupied control
points are within the required accuracy of the project.
55
Table 27. Adjusted geodetic coordinates for control points used in the Bauang River floodplain validation
Point Height Height Error
Latitude Longitude Constraint
ID (Meter) (Meter)
LU67 N16°31’00.31120” E120°19’49.21075” 49.515 ? e
LUN61 N16°29’30.73239” E120°19’46.68075” 49.954 0.047 LL
UPAGY N16°23’21.45560” E120°21’21.45371” 50.363 0.089
UPCBN N16°40’01.07082” E120°23’30.27604” 53.288 0.064
UPGSN N16°30’08.22705” E120°24’35.41408” 63.176 0.075
UPNG N16°31’48.15542” E120°23’17.25573” 58.355 0.050
The corresponding geodetic coordinates of the observed points are within the required accuracy, as
demonstrated in Table 27. Based on the results of the computation, the accuracy conditions are satisfied;
hence, the required accuracy for the program was met.
The computed coordinates of the reference and control points utilized in the Bauang River GNSS Static
Survey are indicated in Table 28.
Table 28. Reference and control points used in the Bauang River Static Survey, with their corresponding
locations (Source: NAMRIA, UP-TCAGP)
Geographic Coordinates (WGS 84) UTM Zone 51

Con- Order
Ellip- Eleva-
trol of Accu-
soid tion in
Point racy Latitude Longitude Northing Easting
Height MSL
(m) (m)
1st
LUN-
order, 16°29’30.73239” 120°19’46.68075” 49.954 1825231.473 214915.551 9.213
61
BM
2nd
LU-
order, 16°31’00.31120” 120°19’49.21075” 49.515 1827985.553 215027.12 8.867
67
GCP
UP
UP-
Estab- 16°23’21.45560” 120°21’21.45371” 50.363 1813837.247 217579.313 9.197
AGY
lished
UP
UP-
Estab- 16°40’01.07082” 120°23’30.27604” 53.288 1844530.559 221802.212 12.83
CBN
lished
UP
UP-
Estab- 16°30’08.22705” 120°24’35.41408” 63.176 1826272.861 223497.8 21.702
GST
lished
UP
UP-
Estab- 16°31’48.15542” 120°23’17.25573” 58.355 1829376.051 221218.673 17.207
NG
lished
56
4.5 Cross-section and Bridge As-Built Survey and Water Level Marking
The cross-section and bridge as-built surveys were conducted at the upstream side of two (2) separate
bridges along the Bauang River. The first was at the Guesset Bridge on July 9, 2015, with control point
UP-GST as the base station; and the second was at the Mamat-ing Bridge on July 15, 2015, with UP-MMG
as the base station. Both surveys were conducted using Trimble® SPS 882 GNSS in PPK survey technique
(Figure 40).
Figure 40. Cross-section survey at the Guesset Bridge
The length of the cross-sectional line surveyed in the Guesset Bridge was 306.47 meters, with a total of
fifty-six (56) points; while that of the Mamat-ing Bridge was 263.42 meters, forty-three (43) points. UP-
GST and UP-MMG were utilized as the GNSS base stations, respectively. The location maps, cross-section
diagrams, and the bridge data forms for both bridges are presented in Figure 42 to Figure 47.
57
Figure 41. Cross-section survey at the Mamat-ing Bridge
The water surface elevation in MSL of the Bauang River was determined on July 9, 2015 at 10:59 hrs.,
using Trimble® SPS 882 in PPK mode. A value of 5.12 meters in MSL was obtained, and was translated into
markings on the Guesset Bridge’s pier using a digital level scale (Figure 41). This was used by the UPB Phil-
LiDAR 1 team as reference for flow data gathering and depth gauge deployment for the Bauang River.
58
Figure 42. Guesset Bridge cross-section location map
59
Figure 43. Bridge cross-section diagram of the Guesset Bridge
60
Figure 44. Mamat-ing Bridge cross-section location map
61
Figure 45. Bridge cross-section diagram of the Mamat-ing Bridge
62
Figure 46. Bridge data form for the Guesset Bridge
63
Figure 47. Bridge data form for the Mamat-ing Bridge
64
4.6 Validation Points Acquisition Survey
The validation points acquisition survey was conducted on July 13-14, 2015 using a survey-grade GNSS
Rover receiver, Trimble® SPS 985. The receiver was mounted on a pole that was attached in front of a
vehicle, as depicted in Figure 48. It was secured with a nylon rope and tied to the side mirrors of the vehicle
to ensure that it was horizontally and vertically balanced. The PPK technique utilized for the conduct of
the survey was set to continuous top mode, with UP-ABY, UP-CBN, and UP-NG occupied as the GNSS base
stations.
Figure 48. Set-up for the validation points acquisition in the Bauang River Basin
The validation points acquisition survey for the Bauang River Basin traversed eight (8) municipalities and
San Fernando City in the province of La Union. The route of the survey aimed to perpendicularly traverse
the LiDAR flight strips for the basin. A total of 3,824 points was acquired for the validation points acquisition
survey, with an approximate length of 51.482 kilometers, as illustrated in the map in Figure 49.
65
Figure 49. Extent of the LiDAR ground validation survey of the Bauang River Basin
66
4.7 Bathymetric Survey
A bathymetric survey of the Bauang River was conducted on July 14, 2015 using an OHMEX™ and a Trimble®
SPS 882 GNSS rover receiver attached to a pole on the side of the boat, as demonstrated in Figure 50.
The survey began in the upstream part of the river in Barangay Acao, with coordinates 16°31’09.7466”N
120°20’55.5551”E; and ended at the mouth of the river in Barangay Pugo with coordinates 16°31’23.9111”N
120°18’36.9876”E. Both points are located in the Municipality of Bauang.
A manual bathymetric survey, on the other hand, was executed on November 26, 2015 using a Trimble®
SPS 882 in GNSS RTK survey technique. The survey began at the upstream portion of the river in Barangay
Cabaritan in the Municipality of Naguilian, with coordinates 16°31’52.88892” 120°23’18.53929”; and
ended at the starting point of the bathymetric survey by boat in Barangay Acao, Municipality of Bauang.
The UP-established control point UP-NG, was occupied as the base station all throughout the surveys.
Figure 50. Bathymetry survey set-up for Bauang River
67
The bathymetric survey for the Bauang River gathered a total of 10,646 points, covering 11.189 kilometers
of the river, traversing the six (6) barangays in Bauang and two (2) barangays in Naguilian, as displayed in
Figure 51.
A CAD drawing was also produced to illustrate the riverbed profile of the Bauang River, as presented in
Figure 52. The profile demonstrates that the elevation was steady, with a total elevation difference of
approximately 11 meters for the surveyed area, except for the sudden drop near Barangay Nagrebcan. The
highest elevation observed was 8 meters, located in Barangay Cabaritan, Municipality of Naguilian; while
the lowest was -2.874 meters, located in Barangay Pugo, Municipality of Bauang.
68
Figure 51. Extent of the bathymetric survey of the Bauang River Basin
69
Figure 52. Bauang centerline riverbed profile
70
CHAPTER 5: FLOOD MODELING AND MAPPING

Dr. Alfredo Mahar Lagmay, Christopher Uichanco, Sylvia Sueno, Marc Moises, Hale Ines, Miguel del
Rosario, Kenneth Punay, Neil Tingin, and Hannah Aventurado
The methods applied in this Chapter were based on the DREAM methods manual (Lagmay, et al., 2014)
5.1 Data Used for Hydrologic Modeling
5.1.1 Hydrometry and Rating Curves

Rainfall, water level, and flow in a certain period of time, which are components and data that may affect
the hydrologic cycle of the Bauang River Basin, were monitored, collected, and analyzed.
5.1.2 Precipitation
Precipitation data was taken from an automatic rain gauge (ARG) installed by the Department of Science
and Technology – Advanced Science and Technology Institute (DOST-ASTI). This said rain gauge is the
Sablan ARG (16°29’50.63” N, 120°29’13.21” E), located in Sablan, Benguet (Figure 53). The precipitation
data collection started on July 4, 2015 at 00:00 hrs. until July 8, 2015 at 23:45 hrs., with a fifteen-minute
recording interval.
The total precipitation for this event in the Sablan ARG was 286.2 millimeters. It had a peak rainfall of 8.6
millimeters on July 5, 2015 at 06:45 hrs. The lag time between the peak rainfall and the discharge was
three (3) hours and fifty-five (55) minutes.
71
Figure 53. The location map of the Bauang HEC-HMS model, which was used for calibration
5.1.3 Rating Curves and River Outflow

A rating curve was computed at the prevailing cross-section (Figure 54) at the Mamat-Ing Bridge, San
Fernando City, La Union (16°34’46.44” N, 120°24’39.42” E). It establishes the relationship between the
observed water levels (H) from the Mamat-Ing Bridge and outflow (Q) of the watershed at this location.
For the Mamat-Ing Bridge, the rating curve was expressed as Q = 1.0901E-06e0.54736x, as illustrated in Figure
55.
72
Figure 54. Cross-section plot of the Mamat-Ing Bridge
Figure 55. Rating curve at the Mamat-Ing Bridge, San Fernando City, La Union
73
The rating curve equation was used to compute for the river outflow at the Mamat-Ing Bridge, for the
calibration of the HEC-HMS model for Bauang presented in Figure 56. The total rainfall for this event was
286.2 millimeters, and the peak discharge was 119.53 m3/s on July 5, 2015 at 10:40 hrs.
Figure 56. Rainfall and outflow data at the Mamat-Ing Bridge, which were used for modeling
5.2 RIDF Station
The Philippine Atmospheric Geophysical and Astronomical Services Administration (PAGASA) computed for
the Rainfall Intensity Duration Frequency (RIDF) values for the Baguio Rain Gauge (Table 29). This station
was selected based on its proximity to the Bauang watershed (Figure 57). The RIDF rainfall amount for
twenty-four (24) hours was converted into a synthetic storm by interpolating and re-arranging the values
such that certain peak values were attained at a certain time. The extreme values for this watershed were
computed based on a 59-year record.
74
Table 29. RIDF values for the Baguio Rain Gauge, computed by PAGASA
COMPUTED EXTREME VALUES (in mm) OF PRECIPITATION
T (yrs) 10 mins 20 mins 30 mins 1 hr 2 hrs 3 hrs 6 hrs 12 hrs 24 hrs
2 27.4 41.6 51.9 72.5 108 136.3 199.5 258.6 355.1
5 49.3 75.2 94.2 127.7 189.7 235.7 334.8 436.9 563.6
10 63.8 97.5 122.1 164.3 243.8 301.5 424.3 555 701.7
15 72 110 137.9 184.9 274.4 338.6 474.8 621.6 779.6
20 77.7 118.8 149 199.3 295.7 364.6 510.2 668.2 834.1
25 82.1 125.6 157.5 210.5 312.2 384.6 537.5 704.1 876.1
50 95.8 146.4 183.7 244.7 362.9 446.3 621.4 814.8 1005.5
100 109.3 167.1 209.7 278.7 413.2 507.5 704.7 924.7 1134
Figure 57. Baguio RIDF location relative to the Bauang River Basin
75
Figure 58. Synthetic storm generated from a 24-hour period rainfall, for various return periods
5.3 HMS Model
The soil shapefile was taken in 2004 from the Bureau of Soils and Water Management (BSWM) under
the Department of Agriculture (DA). The land cover dataset is from the National Mapping and Resource
information Authority (NAMRIA). The soil and land cover maps of the Bauang River Basin are displayed in
Figures 59 and 60, respectively.
76
Figure 59. Soil map of the Bauang River Basin (Source: DA)
77
Figure 60. Land cover map of the Bauang River Basin (Source: NAMRIA)
Six (6) soil classes were identified in the Bauang River Basin. These are clay, clay loam, loam, sand, silt
loam, and undifferentiated soil. Moreover, nine (9) land cover classes were identified. These are brush
lands, built-up areas, closed canopies, cultivated areas, grasslands, inland water, open areas, open canopy
forests, and tree plantations.
78
Figure 61. Slope map of the Bauang River Basin
79
STREAM DELINE-ATION
MAP OF
BAUANG BASIN
Figure 62. Stream delineation map of the Bauang River Basin
Using the SAR-based DEM, the Bauang basin was delineated and further subdivided into sub-basins. The
model consists of sixty (60) sub basins, thirty-one (31) reaches, and thirty-one (31) junctions, as shown in
Figure 63. The main outlet is at 194. See Annex 10 for the Bauang Model Reach Parameters.
80
Figure 63. The Bauang River Basin model, generated using HEC-HMS
5.4 Cross-section Data
Riverbed cross-sections of the watershed were necessary in the HEC-RAS model set-up. The cross-section
data for the HEC-RAS model was derived from the LiDAR DEM data. It was defined using the Arc GeoRAS
tool, and was post-processed in ArcGIS (Figure 64).
81
Figure 64. River cross-section of the Bauang River, generated through the ArcMap HEC GeoRAS tool
5.5 Flo 2D Model
The automated modeling process allowed for the creation of a model with boundaries that are almost
exactly coincidental with that of the catchment area. As such, they have approximately the same land
area and location. The entire area was divided into square grid elements, 10 meters by 10 meters in size.
Each element was assigned a unique grid element number, which served as its identifier. The elements
were then attributed with the parameters required for modeling, such as x- and y-coordinates of centroid,
names of adjacent grid elements, Manning’s coefficient of roughness, infiltration, and elevation values.
The elements were arranged spatially to form the model, allowing the software to simulate the flow of
water across the grid elements in eight (8) directions (i.e., north, south, east, west, northeast, northwest,
southeast, and southwest).
Based on the elevation and flow direction, it was determined that the water will generally flow from the
east of the model to the west, following the main channel. As such, boundary elements in those particular
regions of the model were assigned as inflow and outflow elements, respectively.
82
Figure 65. A screenshot of the river sub-catchment with the computational area to be modeled in FLO-2D
Grid Developer System Pro (FLO-2D GDS Pro)
The simulation was then run through the FLO-2D GDS Pro. This particular model had a computer run time
of 94.73828 hours. After the simulation, the FLO-2D Mapper Pro was used to transform the simulation
results into spatial data that shows the flood hazard levels, as well as the extent and inundation. Assigning
the appropriate flood depths and velocity values for Low, Medium, and High generated the flood hazard
map. Most of the default values given by the FLO-2D Mapper Pro were used, except for those in the Low
hazard level. For this particular level, the minimum h (maximum depth) was set at 0.2 meters; while the
minimum vh (product of maximum velocity (v) and maximum depth (h)) was set at 0 m2/s. The generated
flood hazard maps for Bauang are in Figures 69, 71, and 73.
The creation of a flood hazard map from the model also automatically generated a flow depth map,
depicting the maximum amount of inundation for every grid element. The legend used by default in the
Flo-2D Mapper was not a good representation of the range of flood inundation values, so a different
legend was used for the layout. In this particular model, the inundated parts covered a maximum land area
of 45654200.00 m2. The generated flood depth maps for Bauang are in Figures 70, 72, and 74.
There was a total of 374159128.62 m3 of water that entered the model. Of this amount, 40277538.13 m3
was due to rainfall, while 333881590.49 m3 was inflow from areas outside the model. 8907278.00 m3 of
this water was lost to infiltration and interception, while 78375370.45 m3 was stored by the flood plain.
The rest, amounting to up to 286876474.75 m3, was outflow.
5.6 Results of HMS Calibration
After calibrating the Bauang HEC-HMS River Basin model, its accuracy was measured against the observed
values. Figure 66 depicts the comparison between the two discharge data. The Bauang Model Basin
Parameters are available in Annex 9.
83
Figure 66. Outflow hydrograph of Bauang produced by the HEC-HMS model, compared with observed
outflow
Enumerated in Table 30 are the adjusted ranges of values of the parameters used in calibrating the model.
84
Table 30. Range of calibrated values for the Bauang River Basin model
Hydrologic Calculation Range of Cali-

Method Parameter
Element Type brated Values
Initial Abstraction (mm) 8.84 - 454

Loss SCS Curve Number
Curve Number 35 - 66.58
Time of Concentration (hr) 0.22 - 4.61

Basin Transform Clark Unit Hydrograph
Storage Coefficient (hr) 0.25 – 8.076
Recession Constant 0.588 - 1

Baseflow Recession
Ratio to Peak 0.15 – 0.33
Reach Routing Muskingum-Cunge Manning’s Coefficient 0.0001 - 0.0042
The initial abstraction defines the amount of precipitation that must fall before surface runoff. The
magnitude of the outflow hydrograph increases as the initial abstraction decreases. A range of values from
8.84 to 454 millimeters means that the amount of infiltration or rainfall interception by vegetation all over
the basin varies greatly.
The curve number is the estimate of the precipitation excess of soil cover, land use, and antecedent
moisture. The magnitude of the outflow hydrograph increases as the curve number increases. The range
of 65 to 90 for the curve number is advisable for Philippine watersheds, depending on the soil and land
cover of the area (M. Horritt, personal communication, 2012). For Bauang, the basin consists mainly of
brush lands; and the soil consists of mostly undifferentiated soil and clay.
The time of concentration and storage coefficient are the travel time and the index of temporary storage
of runoff in a watershed. The range of calibrated values from 0.22 hours to 8.076 hours determines the
reaction time of the model, with respect to the rainfall. The peak magnitude of the hydrograph decreases
when these parameters are increased.
The recession constant is the rate at which the baseflow recedes between storm events; and ratio to peak
is the ratio of the baseflow discharge to the peak discharge. Recession constant values within the range
of 0.588 to 1 imply that the basin is unlikely to quickly revert to its original discharge, and will be higher
instead. Values of ratio to peak within the range of 0.15 to 0.33 indicate a steeper receding limb of the
outflow hydrograph.
The Manning’s roughness coefficients correspond to the common roughness of Philippine watersheds.
The Bauang River Basin reaches’ Manning’s coefficients range from 0.0001 to 0.0042, signifying that there
is a variety of surface roughness all over the catchment (Brunner, 2010).
Table 31. Efficiency Test of the Bauang HMS Model

RMSE 3.9
r2 0.9738
NSE 0.97
PBIAS 0.88
RSR 0.17
The Root Mean Square Error (RMSE) method aggregates the individual differences of these two
measurements. It was computed as 3.9 m3/s.
85
The Pearson correlation coefficient (r2) assesses the strength of the linear relationship between the
observations and the model. A coefficient value close to 1 represents an almost perfect match of the
observed discharge and the resulting discharge from the HEC-HMS model. Here, it measured was at 0.9738.
The Nash-Sutcliffe (E) method was also used to assess the predictive power of the model. Here, the optimal
value is 1. The model attained an efficiency coefficient of 0.97.
A positive Percent Bias (PBIAS) indicates a model’s propensity towards under-prediction. Negative values
indicate a bias towards over-prediction. The optimal value is 0. In the model, the PBIAS is 0.88.
The Observation Standard Deviation Ratio (RSR) is an error index. A perfect model attains a value of 0
when the error units of the values are quantified. The model attained an RSR value of 0.17.
86
5.7 Calculated outflow hydrographys and Discharge values for different rainfall
return periods
5.7.1 Hydrograph using the Rainfall Runoff Model

The summary graph (Figure 67) illustrates the Bauang outflow using the Baguio RIDF curves in five (5)
different return periods (i.e., 5-year, 10-year, 25-year, 50-year, and 100-year rainfall time series), based on
the data from PAGASA. The simulation results reveal a significant increase in outflow magnitude as the
rainfall intensity increases, for a range of durations and return periods.
Figure 67. Outflow hydrograph at the Bauang Station, generated using the Baguio RIDF simulated in HEC-
HMS
A summary of the total precipitation, peak rainfall, peak outflow, and time to peak of the Bauang discharge
using the Baguio RIDF curves in five (5) different return periods is provided in Table 32.
87
Table 32. Peak values of the Bauang HEC-HMS Model outflow, using the Baguio RIDF
RIDF Pe- Total Precipitation Peak rainfall Peak outflow

Time to Peak
riod (mm) (mm) (m3/s)
8 hours
5-Year 563.85 49.3 743.2
20 minutes
6 hours
10-Year 701.7 63.8 1448.7
40 minutes
5 hours
25-Year 876.1 82.1 2509.1
10 minutes
4 hours
50-Year 1005.5 95.8 3376.7
40 minutes
4 hours
100-Year 1134 109.3 4281.4
20 minutes
5.8 River Analysis (RAS) Model Simulation
The HEC-RAS flood model produced a simulated water level at every cross-section, for every time step, for
every flood simulation created. The resulting model will be used in determining the flooded areas within
the model. The simulated model will be an integral part in determining the real-time flood inundation
extent of the river, after it has been automated and uploaded on the DREAM website. For this publication,
only a sample output map river is presented. The sample generated map of the Bauang River using the
calibrated HMS base flow is shown in Figure 68.

Figure 68. Sample output map of the Bauang RAS Model
88
5.9 Flow Depth and Flood Hazard
The resulting hazard and flow depth maps have a 10-meter resolution. The 5-year, 25-year, and 100-year
rain return scenarios of the Bauang floodplain are shown in Figures 69 to 74. The floodplain, with an
area of 1,027.9 square kilometers, covers nineteen (19) municipalities from two (2) provinces. Table 33
summarizes the percentage of area affected by flooding per municipality.
Table 33. Municipalities affected in the Bauang floodplain

Area
Province Municipality Total Area % Flooded
Flooded
Benguet Baguio City 60.8979 33.175 54.48%
Benguet Kapangan 133.9 44.5919 33.30%
Benguet La Trinidad 74.908 63.2633 84.45%
Benguet Sablan 90.2168 90.2168 100.00%
Benguet Tuba 322.02 137.594 42.73%
Benguet Tublay 63.209 35.2287 55.73%
La Union Agoo 33.7074 2.6153 7.76%
La Union Aringay 95.6482 95.2867 99.62%
La Union Bagulin 77.9673 65.656 84.21%
La Union Bauang 85.26 84.1416 98.69%
La Union Burgos 51.9196 51.9196 100.00%
La Union Caba 56.1911 56.1165 99.87%
La Union Naguilian 86.3875 86.3875 100.00%
La Union Pugo 60.537 50.1692 82.87%
La Union Rosario 64.3275 4.53607 7.05%
La Union San Fernando City 121.053 56.8372 46.95%
La Union San Gabriel 154.189 16.9234 10.98%
La Union Santo Tomas 58.5311 5.21613 8.91%
La Union Tubao 53.8728 45.1865 83.88%
89
Figure 69. 100-year flood hazard map for the Bauang floodplain, overlaid on Google Earth imagery
90
Figure 70. 100-year flow depth map for the Bauang floodplain, overlaid on Google Earth imagery
91
92
93
94
95
5.10 Inventory of Areas Exposed to Flooding
The affected barangays in the Bauang River Basin, grouped by municipality, are listed below. For the said
basin, two (2) provinces with nineteen (19) municipalities, consisting of three hundred and sixty-six (366)
barangays are expected to experience flooding when subjected to 5-year, 25-year, and 100-year rainfall
return periods.
For the 5-year return period, 48.86% of the Municipality of Baguio City, with an area of 60.8979 square
kilometers, will experience flood levels of less than 0.20 meters. 2.08% of the area will experience flood
levels of 0.21 to 0.50 meters. Meanwhile, 1.18%, 1.01%, 0.72%, and 0.15% of the area will experience
flood depths of 0.51 to 1 meter, 1.01 to 2 meters, 2.01 to 5 meters, and more than 5 meters, respectively.
Listed in Tables 34-42 are the affected areas, in square kilometers, by flood depth per barangay.
96
Table 34. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - A
Affected Barangays in Baguio City, Benguet
Aban-
ARINGAY-BAU- A. Bonifa- Auro- Auro-
ao-Zandu- Aurora
ANG BASIN cio-Cagu- Alfonso Andres Asin ra Hill, ra Hill, Bakakeng Bakakeng Bal-Marco-
eta-Kay- Ambiong Hill Prop-
ioa-Riman- Tabora Bonifacio Road North South Central North ville
ong-Chu- er
do Central Central
gum-Otek
0-0.20 0.054 0.037 0.08 0.24 0.075 4.35 0.025 0.044 0.07 0.72 0.084 0.12
0.21-0.50 0.0031 0.006 0.0026 0.0091 0.003 0.18 0.00085 0.0021 0.00082 0.028 0.004 0.0034
0.51-1.00 0.00036 0.0007 0.0073 0.0073 0.0001 0.082 0.0013 0.0023 0 0.024 0 0.0065
km.)
1.01-2.00 0.00038 0 0.0082 0.0003 0 0.073 0.0024 0.002 0 0.013 0 0.0022
2.01-5.00 0.00021 0 0.0055 0 0 0.031 0.0007 0 0 0.0088 0 0.0001
Affected Area (sq.

> 5.00 0 0 0.00061 0 0 0.0013 0 0 0 0.0022 0 0
Table 35. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - B
97
ARINGAY-BAU- Bayan Cabinet Camdas
ANG BASIN Bayan Bgh Com- Brook- Camp Campo City Camp
Balsigan Park Brookside Hill-Teach- Subdivi- Camp 8
Park East pound spoint Allen Filipino Central
West er’s Camp sion
0-0.20 0.054 0.17 0.076 0.017 0.092 0.1 0.36 0.16 0.09 0.13 0.13 0.069
0.21-0.50 0.000047 0.0053 0.0012 0.00078 0.0035 0.0051 0.0045 0.0012 0.0027 0.0057 0.017 0.0024
0.51-1.00 0.00019 0.0083 0.0027 0.0021 0.0015 0.0055 0.0058 0.00075 0 0.0068 0.024 0.0001
km.)
1.01-2.00 0 0.0022 0.0007 0.00028 0 0.0043 0.014 0.00021 0 0.0023 0.042 0
2.01-5.00 0 0 0 0 0 0.0011 0.0046 0.0011 0 0.00012 0.0001 0
Affected Area (sq.

> 5.00 0 0 0 0 0 0 0 0 0 0 0 0
Table 36. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - C
ARINGAY-BAU- Domin- General
Country General
ANG BASIN City Camp Cresencia Dizon Sub- ican Engineers’ Fairview Luna,
Club Vil- Dontogan Dps Area Ferdinand Emilio F.
Proper Village division Hill-Mira- Hill Village Lower
lage Aguinaldo
dor
0-0.20 0.085 0.2 0.17 0.26 1.23 0.56 0.12 0.1 0.26 0.095 0.065 0.029
0.21-0.50 0.00089 0.012 0.0055 0.0025 0.042 0.009 0.0015 0.0012 0.02 0.008 0.0014 0.0023
0.51-1.00 0 0.0051 0.0076 0.001 0.022 0.0055 0.0027 0 0.019 0.0047 0.0013 0.0016
km.)
1.01-2.00 0 0.0052 0.013 0.0015 0.025 0.0057 0.0055 0 0.014 0.0012 0.0024 0.00099
2.01-5.00 0 0.0007 0.0014 0.0011 0.055 0.0087 0.021 0 0.0024 0 0.0034 0.0029
Affected Area (sq.

> 5.00 0 0 0 0 0.043 0.000098 0 0 0 0 0 0
Table 37. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - D
98
ARINGAY-BAU- General Holy Imelda

ANG BASIN Green-wa- Guisad Guisad Happy Happy Holy Ghost Hon- Imelda R.
Luna, Gibraltar Ghost Village
ter Village Central Sorong Hollow Homes Extension ey-moon Marcos
Upper Proper
0-0.20 0.034 0.2 0.13 0.077 0.3 0.78 0.19 0.3 0.047 0.073 0.025 0.08
0.21-0.50 0.0014 0.0022 0.0033 0.0004 0.012 0.039 0.002 0.009 0.0029 0.0043 0.0001 0.0064
0.51-1.00 0.000007 0.00024 0.00055 0.0002 0.0041 0.014 0.0017 0.0059 0.0012 0.0014 0 0.0048
km.)
1.01-2.00 0 0.000071 0.00013 0 0.0012 0.015 0.0035 0.013 0.0023 0.0007 0 0.0056
2.01-5.00 0 0 0 0 0.00019 0.0045 0.012 0.008 0.0091 0 0 0.0018
Affected Area (sq.

> 5.00 0 0 0 0 0 0.0001 0.001 0 0 0 0 0
Table 38. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - E
ARINGAY-BAU- Lourdes Lourdes Mag-
Legar- Lourdes
ANG BASIN Kabaya-ni- Kagi- Kayang Kayang-Hill- Lopez Sub-di- Sub-di- say-say
Irisan da-Burn- Subdivision Lualhati
han ti-ngan Extension top Jaena vision, vision, Private
ham-Kisad Extension
Lower Proper Road
0-0.20 5.77 0.062 0.033 0.091 0.01 0.37 0.047 0.028 0.032 0.081 0.33 0.027
0.21-0.50 0.23 0.0025 0.00081 0.004 0.00019 0.025 0.0023 0.005 0.0023 0.0017 0.016 0.0039
0.51-1.00 0.14 0.0047 0.0015 0.00065 0.000001 0.024 0 0 0 0 0.016 0.00077
km.)
1.01-2.00 0.11 0.0018 0.0024 0.000086 0 0.04 0 0 0 0 0.006 0.0031
2.01-5.00 0.086 0 0.0054 0 0 0.04 0 0 0 0 0 0.0066
Affected Area (sq.

> 5.00 0.026 0 0 0 0 0.0054 0 0 0 0 0 0
99
Table 39. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - F
Middle
ARINGAY-BAU- Mag- Malcolm Que- Modern
ANG BASIN Magsay-say, Manuel A. Military Mrr-Queen New Outlook Padre
say-say, Square-Per- zon Hill Site, Pacdal
Upper Roxas Cut-Off Of Peace Lucban Drive Burgos
Lower fecto Sub-divi- West
sion
0-0.20 0.097 0.0096 0.045 0.099 0.15 0.29 0.098 0.062 0.12 0.14 1.1 0.12
0.21-0.50 0.011 0.0007 0.00012 0.0038 0.0072 0.017 0.0047 0.0036 0.0087 0.0037 0.044 0.0045
0.51-1.00 0.0025 0.0028 0 0.0072 0.0008 0.0038 0.0002 0.0011 0.0076 0.0029 0.018 0.0082
km.)
1.01-2.00 0.0068 0.0031 0 0.011 0 0.00097 0 0.00051 0.0088 0.00073 0.003 0.0044
2.01-5.00 0.0087 0.00092 0 0.0022 0 0 0 0 0.0069 0 0.0006 0.01
Affected Area (sq.

> 5.00 0 0 0 0 0 0 0 0 0 0 0 0
Table 40. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - G
100

ARINGAY-BAU- Pinsao Quezon Quezon Quirino
ANG BASIN Padre Palma-Ur- Pinsao Quirino
Phil-Am Pinget Pilot Proj- Poliwes Pucsusan Hill Prop- Hill, Up- Hill, Low-
Zamora bano Proper Hill, East
ect er per er
0-0.20 0.031 0.077 0.034 0.36 0.094 2.18 0.063 0.15 0.17 0.32 0.055 0.034
0.21-0.50 0.0054 0.00029 0.0043 0.014 0.0019 0.092 0.0024 0.0091 0.019 0.021 0.0029 0.0028
0.51-1.00 0.0093 0 0.0026 0.0055 0.0009 0.038 0.0003 0.0028 0.015 0.007 0.0041 0.0021
km.)
1.01-2.00 0.025 0 0.0022 0.0054 0.0013 0.019 0 0.0001 0.007 0.0013 0.0045 0.0003
2.01-5.00 0.024 0 0 0.0059 0.001 0.0064 0 0 0 0.0005 0.0077 0.0008
Affected Area (sq.

> 5.00 0 0 0 0.0037 0.00014 0.0024 0 0 0 0 0 0
Table 41. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - H
ARINGAY-BAU- Quiri- Rizal
Quirino Rock Rock Saint San San
ANG BASIN no-Mag- Monu- Rock Quar- Salud San Luis San Vicen-
Hill, Mid- Quarry, Quarry, Joseph Antonio Roque
say-say, ment ry, Upper Mitra Village te
dle Lower Middle Village Village Village
Upper Area
0-0.20 0.058 0.1 0.074 0.037 0.059 0.065 0.65 0.15 0.13 0.62 0.057 0.2
0.21-0.50 0.0092 0.003 0.0027 0.0026 0.003 0.011 0.03 0.0048 0.0031 0.022 0.0019 0.0088
0.51-1.00 0.0029 0.0026 0.0011 0.0011 0.00033 0.00076 0.022 0.000093 0.0003 0.008 0 0.0041
km.)
1.01-2.00 0.0009 0.0056 0 0 0.000014 0 0.01 0 0 0.0054 0 0.0036
2.01-5.00 0.0002 0.019 0 0 0 0 0.0031 0 0 0.0002 0 0.0027
Affected Area (sq.

> 5.00 0 0.0048 0 0 0 0 0 0 0 0 0 0
Table 42. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - I
101
ARINGAY-BAU- Sanitary Sanitary Santo Slaughter Slu-Svp
ANG BASIN Santa Es- Santo Session South Teodora Tran- Victoria
Camp, Camp, Tomas House Housing
colas-tica Rosario Road Area Drive Alonzo co-ville Village
North South Proper Area Village
0-0.20 0.15 0.085 0.08 0.2 0.22 0.21 0.025 0.21 0.34 0.058 0.23 0.19
0.21-0.50 0.0029 0.006 0.0037 0.0083 0.0074 0.006 0.00038 0.0035 0.034 0.0053 0.011 0.012
0.51-1.00 0.00046 0.0013 0.0016 0.0011 0.0024 0.013 0.00063 0.0017 0.015 0.0026 0.0072 0.0013
km.)
1.01-2.00 0 0 0.0026 0.00029 0.0005 0.0081 0.0016 0.0032 0.0086 0.0071 0.0081 0.0008
2.01-5.00 0 0 0.0002 0 0 0 0.0041 0.0015 0.004 0.0032 0.0041 0.0001
Affected Area (sq.

> 5.00 0 0 0 0 0 0 0 0 0 0 0.00059 0
Figure 75. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - A
Figure 76. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - B
102
Figure 77. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - C
Figure 78. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - D
103
Figure 79. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - E
Figure 80. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - F
104
Figure 81. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - G
Figure 82. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - H
105
Figure 83. Affected areas in Baguio City, Benguet during a 5-year rainfall return period - I
For the 5-year return period, 31.09% of the Municipality of Kapangan, with an area of 133.9 square
Listed in Table 43 are the affected areas, in square kilometers, by flood depth per barangay.
Table 43. Affected areas in Kapangan, Benguet during a 5-year rainfall return period
ARIN- Affected Barangays in Kapangan, Benguet

GAY-BAU- Gas- Poblacion Pon-
ANG BASIN Datakan Gadang Labueg Paykek Sagubo
we-ling Central ga-yan
0-0.20 3.91 0.19 15.45 2.06 0.036 1.88 9.9 8.2
0.21-
Affected Area (sq. km.)
0.12 0.0006 0.32 0.1 0.0013 0.068 0.36 0.3

0.50
0.51-
0.061 0 0.15 0.082 0.000037 0.048 0.17 0.12
1.00
1.01-
0.033 0 0.12 0.062 0.000038 0.05 0.1 0.055
2.00
2.01-
0.053 0 0.12 0.067 0 0.041 0.064 0.027
5.00
> 5.00 0.094 0 0.15 0.013 0 0.00093 0.0011 0.0001
106
Figure 84. Affected areas in Kapangan, Benguet during a 5-year rainfall return period
For the 5-year return period, 73.97% of the Municipality of La Trinidad, with an area of 74.908 square
Listed in Tables 44-45 are the affected areas, in square kilometers, by flood depth per barangay.
107
Table 44. Affected areas in La Trinidad, Benguet during a 5-year rainfall return period - A
ARINGAY-BAU- Affected Barangays in La Trinidad, Benguet

ANG BASIN Alapang Alno Ambiong Bahong Balili Beckel Betag Bineng
0-0.20 4.11 10.09 1.59 4.54 0.34 1.48 0.21 11.2
0.21-0.50 0.19 0.33 0.055 0.28 0.11 0.055 0.11 0.43
0.51-1.00 0.097 0.15 0.039 0.13 0.2 0.019 0.17 0.13
km.)
1.01-2.00 0.096 0.088 0.041 0.072 0.21 0.011 0.47 0.056
2.01-5.00 0.072 0.14 0.01 0.04 0.11 0.0016 0.59 0.029
Affected Area (sq.

> 5.00 0.0025 0.53 0.0002 0.0021 0.0012 0 0.0044 0.1
Table 45. Affected areas in La Trinidad, Benguet during a 5-year rainfall return period - B
ARINGAY-BAU- Affected Barangays in La Trinidad, Benguet

ANG BASIN Cruz Lubas Pico Poblacion Puguis Shilan Tawang Wangal
0-0.20 0.47 1.43 1.11 1.85 3.58 6.69 2.29 4.41
108
0.21-0.50 0.0087 0.069 0.36 0.084 0.15 0.2 0.067 0.16

0.51-1.00 0.00089 0.033 0.31 0.058 0.055 0.12 0.014 0.11
km.)
1.01-2.00 0.0004 0.03 0.2 0.071 0.024 0.07 0.01 0.074
2.01-5.00 0 0.02 0.14 0.096 0.0065 0.041 0.0031 0.069
Affected Area (sq.

> 5.00 0 0.0005 0.0067 0 0 0.0071 0.0001 0.009
Figure 85. Affected areas in La Trinidad, Benguet during a 5-year rainfall return period - A
Figure 86. Affected areas in La Trinidad, Benguet during a 5-year rainfall return period - B
For the 5-year return period, 93.11% of the Municipality of Sablan, with an area of 90.22 square kilometers,
will experience flood levels of less than 0.20 meters. 2.68% of the area will experience flood levels of 0.21
to 0.50 meters. Meanwhile, 1.36%, 0.90%, 0.71%, and 1.24% of the area will experience flood depths of
0.51 to 1 meter, 1.01 to 2 meters, 2.01 to 5 meters, and more than 5 meters, respectively. Listed in Table
109
46 are the affected areas, in square kilometers, by flood depth per barangay.
Table 46. Affected areas in Sablan, Benguet during a 5-year rainfall return period
Affected Barangays in Sablan
ARINGAY-BAUANG
BASIN Baneng-
Bagong Balluay Banangan Bayabas Kamog Pappa Poblacion
beng
0-0.20 5 8.66 17.83 13.08 12.79 10.91 10.45 5.3
0.21-0.50 0.083 0.29 0.5 0.39 0.43 0.26 0.28 0.18
0.51-1.00 0.04 0.16 0.25 0.18 0.2 0.15 0.16 0.095
1.01-2.00 0.038 0.077 0.19 0.1 0.12 0.096 0.11 0.074
(sq. km.)
Affected Area
2.01-5.00 0.078 0.035 0.14 0.13 0.068 0.032 0.12 0.031
> 5.00 0.71 0.08 0.013 0.27 0.012 0.034 0.0048 0
110
Figure 87. Affected areas in Sablan, Benguet during a 5-year rainfall return period
111
For the 5-year return period, 39.63% of the Municipality of Tuba, with an area of 322.02 square kilometers, will experience flood levels of less than 0.20 meters. 1.39% of
the area will experience flood levels of 0.21 to 0.50 meters. Meanwhile, 0.60%, 0.36%, 0.43%, and 0.32% of the area will experience flood depths of 0.51 to 1 meter, 1.01
to 2 meters, 2.01 to 5 meters, and more than 5 meters, respectively. Listed in Table 47 are the affected areas, in square kilometers, by flood depth per barangay.
Table 47. Affected areas in Tuba, Benguet during a 5-year rainfall return period
Affected Barangays in Tuba
ARINGAY-BAUANG
BASIN Tabaan Tabaan Taloy
Nangalisan Poblacion San Pascual Tadiangan Taloy Sur Twin Peaks
Norte Sur Norte
0-0.20 10.95 20.69 6.44 16.95 23.47 12.08 14.8 19.58 2.66
0.21-0.50 0.41 0.85 0.2 0.55 0.83 0.45 0.48 0.6 0.11
0.51-1.00 0.21 0.36 0.097 0.18 0.37 0.21 0.21 0.26 0.027
1.01-2.00 0.17 0.18 0.085 0.098 0.19 0.14 0.12 0.17 0.0011
(sq. km.)
Affected Area
2.01-5.00 0.29 0.25 0.057 0.1 0.26 0.16 0.11 0.18 0.0003
> 5.00 0.5 0.032 0.007 0.13 0.15 0.11 0.045 0.062 0
112
Figure 88. Affected areas in Tuba, Benguet during a 5-year rainfall return period
113
For the 5-year return period, 51.58% of the Municipality of Tublay, with an area of 63.21 square kilometers, will experience flood levels of less than 0.20 meters. 1.90% of
Table 48. Affected areas in Tublay, Benguet during a 5-year rainfall return period
Affected Barangays in Tublay
ARINGAY-BAUANG
BASIN Ambassa- Ambong- Tublay
Ba-Ayan Basil Caponga Daclan Tuel
dor dolan Central
0-0.20 0.023 3.74 2.05 3.89 2.96 4.95 3.56 11.43
0.21-0.50 0.0001 0.13 0.065 0.13 0.17 0.16 0.15 0.39
0.51-1.00 0 0.08 0.019 0.05 0.064 0.072 0.11 0.17
1.01-2.00 0 0.076 0.0097 0.03 0.028 0.032 0.083 0.11
(sq. km.)
Affected Area
2.01-5.00 0 0.017 0.0004 0.05 0.0035 0.0022 0.083 0.13
> 5.00 0 0.014 0 0.089 0 0 0.012 0.096
114
Figure 89. Affected areas in Tublay, Benguet during a 5-year rainfall return period
For the 5-year return period, 6.53% of the Municipality of Agoo, with an area of 33.71 square kilometers,
to 0.50 meters. Meanwhile, 0.38%, 0.23%, and 0.15% of the area will experience flood depths of 0.51 to
1 meter, 1.01 to 2 meters, and 2.01 to 5 meters, respectively. Listed in Table 49 are the affected areas, in
square kilometers, by flood depth per barangay.
Table 49. Affected areas in Agoo, La Union during a 5-year rainfall return period
Affected Barangays in Agoo
ARINGAY-BAUANG
BASIN Macalva Santa Rita Santa Rita
Nazareno Santa Ana Santa Rita
Norte East Norte
0-0.20 1.55 0.016 0.54 0.031 0.012 0.058
0.21-0.50 0.071 0.0002 0.038 0.019 0.0027 0.029
Affected Area
(sq. km.)
0.51-1.00 0.039 0 0.021 0.016 0.0041 0.047

1.01-2.00 0.028 0 0.027 0.0025 0.0046 0.015
2.01-5.00 0.0067 0 0.033 0 0.0045 0.0062
> 5.00 0 0 0.0001 0 0 0
115
Figure 90. Affected areas in Agoo, La Union during a 5-year rainfall return period
For the 5-year return period, 68.77% of the Municipality of Aringay, with an area of 95.65 square kilometers,
0.51 to 1 meter, 1.01 to 2 meters, 2.01 to 5 meters, and more than 5 meters, respectively. Listed in Tables
50-51 are the affected areas, in square kilometers, by flood depth per barangay.
116
Table 50. Affected areas in Aringay, La Union during a 5-year rainfall return period - A
Affected Barangays in Aringay
ARINGAY-BAUANG Pan- Pan-
BASIN Maca-ba- Pobla-ci- San An- San Beni-
Alaska Basca Dulao Gallano Manga gao-Aoan gao-Aoan Samara
to on tonio to Norte
East West
0-0.20 0.45 5.9 0.92 22.74 5.68 7.58 2.11 1.56 0.65 0.063 0.78 1.77
0.21-0.50 0.089 0.24 0.56 0.84 0.2 0.26 0.1 0.14 0.025 0.12 0.12 0.14
0.51-1.00 0.11 0.15 0.75 0.42 0.12 0.16 0.066 0.082 0.015 0.63 0.17 0.031
1.01-2.00 0.84 0.1 1.22 0.25 0.11 0.11 0.094 0.11 0.03 2.05 0.25 0.028
(sq. km.)
Affected Area
2.01-5.00 0.5 0.11 0.019 0.37 0.14 0.053 0.47 0.053 0.079 1.21 0.93 0.014
> 5.00 0 0.12 0 1.49 0.0082 0.0005 1.37 0 0.44 0.04 0.85 0
Table 51. Affected areas in Aringay, La Union during a 5-year rainfall return period- B
ARINGAY-BAUANG San San Santo Santo
117
BASIN San Beni- San Eu- San Juan San Juan Santa Santa Santa Santa
Simon Simon Rosario Rosario
to Sur genio East West Cecilia Lucia Rita East Rita West
East West East West
0-0.20 0.6 5.27 2.19 1.93 1.28 1.68 0.25 0.25 1.63 0.29 0.21 0.0073
0.21-0.50 0.034 0.25 0.23 0.086 0.064 0.14 0.017 0.13 0.065 0.038 0.093 0.025
0.51-1.00 0.024 0.11 0.19 0.067 0.032 0.098 0.012 0.2 0.032 0.0071 0.027 0.038
1.01-2.00 0.039 0.079 0.28 0.1 0.027 0.083 0.02 1.18 0.023 0.0092 0.059 0.23
(sq. km.)
Affected Area
2.01-5.00 0.063 0.048 0.34 0.61 0.026 0.054 0.58 0.8 0.0055 0.0031 0.018 1.06
> 5.00 0.13 0 0.054 0.72 0.096 0.0025 1.52 0.26 0 0 0 0
Figure 91. Affected areas in Aringay, La Union during a 5-year rainfall return period - A
Figure 92. Affected areas in Aringay, La Union during a 5-year rainfall return period - B
118
For the 5-year return period, 75.19% of the Municipality of Bagulin, with an area of 77.97 square kilometers, will experience flood levels of less than 0.20 meters. 2.28% of
Table 52. Affected areas in Bagulin, La Union during a 5-year rainfall return period
Affected Barangays in Bagulin
ARINGAY-BAU-
ANG BASIN Ali- Tagud- Tio-An- Wal-
Baay Cambaly Cardiz Dagup Libbo Suyo
bang-say tud gan la-yan
0-0.20 11.32 1.88 3.15 4.6 5.18 4.83 7.71 2.82 13.64 3.51
0.21-0.50 0.27 0.056 0.1 0.13 0.2 0.12 0.31 0.11 0.35 0.12
0.51-1.00 0.15 0.0071 0.064 0.061 0.13 0.037 0.16 0.05 0.14 0.053
1.01-2.00 0.097 0.0054 0.035 0.022 0.094 0.016 0.094 0.025 0.12 0.032
(sq. km.)
Affected Area
2.01-5.00 0.092 0.014 0.016 0.022 0.086 0.028 0.061 0.019 0.28 0.045
> 5.00 0.85 0.32 0.052 0 0.5 0.082 0.39 0.14 0.3 0.56
119
Figure 93. Affected areas in Bagulin, La Union during a 5-year rainfall return period
120
For the 5-year return period, 64.67% of the Municipality of Bauang, with an area of 85.26 square kilometers, will experience flood levels of less than 0.20 meters. 6.93%
of the area will experience flood levels of 0.21 to 0.50 meters. Meanwhile, 5.64%, 6.71%, 8.26%, and 6.60% of the area will experience flood depths of 0.51 to 1 meter,
1.01 to 2 meters, 2.01 to 5 meters, and more than 5 meters, respectively. Listed in Tables 53-56 are the affected areas, in square kilometers, by flood depth per barangay.
Table 53. Affected areas in Bauang, La Union during a 5-year rainfall return period - A
Affected Barangays in Bauang
ARINGAY-BAU- Caba-
ANG BASIN Baccuit Baccuit Bawan- Boy- Cabi-
Acao Bagbag Ballay Bucayab la-yan-
Norte Sur ta Utan si-lan
gan
0-0.20 3.09 0.17 0.023 0.57 0.28 0.77 3.36 3.02 4.79 1.82
0.21-0.50 0.17 0.11 0.075 0.12 0.0058 0.034 0.29 0.14 0.27 0.15
0.51-1.00 0.15 0.27 0.13 0.084 0.012 0.016 0.15 0.1 0.21 0.058
1.01-2.00 0.29 0.34 0.17 0.086 0.046 0.0043 0.086 0.14 0.26 0.038
(sq. km.)
Affected Area
2.01-5.00 0.65 0.07 0.22 0.039 1.2 0.00056 0.13 0.078 0.53 0.09
> 5.00 2.93 0 0 0.0004 0.76 0 0.67 0.022 0.56 0.0081
121
Table 54. Affected areas in Bauang, La Union during a 5-year rainfall return period - B
ARINGAY-BAU- Lower
ANG BASIN Calum- Casi- Central Central Guerre- Nagreb-
Carmay Dili Disso-Or San
baya la-gan East West ro can
Agustin
0-0.20 1.03 0.027 5.97 0.31 0.0068 0.37 0.51 0.52 2.28 0.22
0.21-0.50 0.2 0.051 0.42 0.021 0.0018 0.052 0.094 0.1 0.11 0.039
0.51-1.00 0.11 0.096 0.29 0.01 0.0035 0.051 0.045 0.0098 0.059 0.054
1.01-2.00 0.09 0.045 0.29 0.023 0.015 0.031 0.032 0.0036 0.049 0.028
(sq. km.)
Affected Area
2.01-5.00 0.0055 0.084 0.2 0.34 0.23 0.0003 0.016 0.0001 0.11 0.15
> 5.00 0 0 0.0035 0.21 0.0087 0 0 0 0.011 0.41
Table 55. Affected areas in Bauang, La Union during a 5-year rainfall return period - C
ARINGAY-BAUANG Palug-
BASIN Payocpoc Payocpoc
Pagdalag- Palintu- si-Lim- Parian Parian Payocpoc
Paringao Norte Norte Pilar
an Sur cang mansan- Este Oeste Sur
Este Oeste
gan
0-0.20 0 0.49 5.16 0.42 0.13 0.023 1.66 0.66 0.47 0.034
0.21-0.50 0 0.061 0.31 0.46 0.21 0.042 0.17 0.48 0.15 0.03
0.51-1.00 0.000003 0.092 0.19 0.32 0.6 0.094 0.13 0.4 0.082 0.014
1.01-2.00 0 0.053 0.089 0.11 0.8 0.27 0.094 0.26 0.065 0.000094
(sq. km.)
Affected Area
2.01-5.00 0 0.044 0.022 0.28 0.13 0.072 0.089 0.0003 0.0021 0
> 5.00 0 0.0017 0.0002 0.0096 0 0 0 0 0 0
Table 56. Affected areas in Bauang, La Union during a 5-year rainfall return period - D
122

ARINGAY-BAUANG Upper
BASIN Santa
Pottot Pudoc Pugo Quinavite Santiago Taberna San Agus- Urayong
Monica
tin
0-0.20 0.55 0.066 0.0002 0.56 5.84 1.24 0.01 8.24 0.46
0.21-0.50 0.15 0.071 0.0033 0.1 0.26 0.2 0.012 0.58 0.16
0.51-1.00 0.089 0.076 0.043 0.033 0.17 0.083 0.044 0.35 0.071
1.01-2.00 0.01 0.41 0.6 0.031 0.13 0.041 0.28 0.39 0.025
(sq. km.)
Affected Area
2.01-5.00 0 0.38 0.78 0.051 0.048 0.0028 0.85 0.16 0.0001
> 5.00 0 0 0.01 0 0.0007 0 0 0.0024 0
Figure 94. Affected areas in Bauang, La Union during a 5-year rainfall return period - A
Figure 95. Affected areas in Bauang, La Union during a 5-year rainfall return period - B
123
Figure 96. Affected areas in Bauang, La Union during a 5-year rainfall return period - C
Figure 97. Affected areas in Bauang, La Union during a 5-year rainfall return period - D
124
For the 5-year return period, 91.67% of the Municipality of Burgos, with an area of 51.92 square kilometers, will experience flood levels of less than 0.20 meters. 2.85% of
Table 57. Affected areas in Burgos, La Union during a 5-year rainfall return period
Affected Barangays in Burgos
ARINGAY-BAUANG
BASIN Lower New Po- Old Po- Upper
Agpay Bilis Caoayan Dalacdac Delles Imelda Libtong Linuan
Tumapoc blacion blacion Tumapoc
0-0.20 3.07 4.75 0.38 2.34 11.75 3.95 4.77 4.12 1.09 3.4 1.36 6.61
0.21-0.50 0.1 0.13 0.011 0.064 0.36 0.14 0.17 0.087 0.035 0.13 0.052 0.2
0.51-1.00 0.049 0.049 0.011 0.035 0.17 0.074 0.065 0.052 0.016 0.031 0.018 0.078
1.01-2.00 0.027 0.029 0.006 0.03 0.11 0.056 0.043 0.033 0.021 0.027 0.022 0.073
(sq. km.)
Affected Area
2.01-5.00 0.028 0.051 0.01 0.12 0.032 0.055 0.05 0.041 0.045 0.013 0.046 0.11
> 5.00 0.077 0.022 0.061 0.28 0.0025 0.057 0.033 0.42 0.0044 0.0055 0.065 0.1
125
Figure 98. Affected areas in Burgos, La Union during a 5-year rainfall return period
126
For the 5-year return period, 74.14% of the Municipality of Caba, with an area of 56.19 square kilometers,
0.51 to 1 meter, 1.01 to 2 meters, 2.01 to 5 meters, and more than 5 meters, respectively. Listed in Tables
58-59 are the affected areas, in square kilometers, by flood depth per barangay.
Table 58. Affected areas in Caba, La Union during a 5-year rainfall return period - A
Affected Barangays in Caba
ARIN-
GAY-BAUANG Pobla- San
Juan Pobla- San
BASIN Bautista Gana Las-Ud Liquicia cion Corne-
Cartas cion Sur Carlos
Norte lio
0-0.20 2.68 1.21 2.44 0.13 13.46 0.44 0.38 0.37 5.23
0.21-
0.35 0.61 0.17 0.037 0.51 0.054 0.096 0.3 0.39
0.50
Affected Area
0.51-
0.13 0.16 0.079 0.042 0.38 0.098 0.027 0.32 0.24
(sq. km.)
1.00
1.01-
0.05 0.1 0.031 0.031 0.42 0.21 0.018 0.27 0.18
2.00
2.01-
0.026 0.055 0.013 0.0063 0.63 0.16 0.011 0.046 0.26
5.00
> 5.00 0.00013 0 0.0003 0 0.84 0.00047 0 0 0.12
Table 59. Affected areas in Caba, La Union during a 5-year rainfall return period - B
ARIN- Affected Barangays in Caba

GAY-BAU- San Fer- San Gre- Santiago Santiago Sobre- Wences-
ANG BASIN San Jose Urayong
min gorio Norte Sur dillo lao
0-0-
5 1.97 2.33 0.084 0.05 3.81 1.9 0.17
.20
0.21-
0.31 0.13 0.2 0.094 0.13 0.14 0.23 0.14
0.50
Affected Area
0.51-
(sq. km.)
0.19 0.11 0.11 0.23 0.29 0.084 0.06 0.083

1.00
1.01-
0.16 0.081 0.17 0.21 1.02 0.042 0.014 0.01
2.00
2.01-
0.27 0.12 0.5 0.044 1.4 0.027 0.0013 0
5.00
> 5.00 0.23 0.077 0.094 0 0 0.014 0 0
127
Figure 99. Affected areas in Caba, La Union during a 5-year rainfall return period - A
Figure 100. Affected areas in Caba, La Union during a 5-year rainfall return period - B
128
For the 5-year return period, 69.65% of the Municipality of Naguilian, with an area of 86.39 square kilometers, will experience flood levels of less than 0.20 meters. 3.41%
Table 60. Affected areas in Naguilian, La Union during a 5-year rainfall return period - A
Affected Barangays in Naguilian
ARINGAY-BAUANG
BASIN Al-Alinao Al-Alinao Ambara- Ambara- Baraoas Baraoas
Aguioas Angin Balecbec Bancagan
Norte Sur cao Norte cao Sur Norte Sur
0-0.20 0.53 2.13 3.12 2.4 1.72 0.24 3.66 0.43 0.082 0.058
0.21-0.50 0.016 0.2 0.11 0.27 0.15 0.011 0.14 0.012 0.083 0.0057
0.51-1.00 0.011 0.22 0.042 0.27 0.11 0.0024 0.091 0.04 0.1 0.044
1.01-2.00 0.025 0.23 0.029 0.21 0.035 0.00048 0.087 0.1 0.13 0.067
(sq. km.)
Affected Area
2.01-5.00 0.2 0.37 0.065 0.74 0.054 0.000085 0.17 0.31 0.26 0.3
> 5.00 0.49 0.66 0.22 1.26 0.026 0 0.12 1.09 0 0
129
Table 61. Affected areas in Naguilian, La Union during a 5-year rainfall return period - B
ARINGAY-BAUANG
BASIN Bimmoto- Cabaritan Cabaritan Dal-Li- Dara- Gusing
Bariquir Bato Casilagan Guesset
bot Norte Sur paoen muangan Norte
0-0.20 1.06 1.64 4.16 0.0012 0.04 5.55 0.0087 0.14 0.44 3.64
0.21-0.50 0.025 0.081 0.16 0.00016 0.0022 0.21 0 0.0062 0.025 0.11
0.51-1.00 0.012 0.073 0.092 0.00082 0.0049 0.11 0.0003 0.0052 0.013 0.065
1.01-2.00 0.0059 0.031 0.052 0.0025 0.022 0.11 0.0004 0.023 0.003 0.05
(sq. km.)
Affected Area
2.01-5.00 0.015 0.019 0.051 0.04 0.18 0.28 0.19 0.23 0 0.083
> 5.00 0.012 0.0028 0.002 0.25 0.092 0.14 0.28 0.49 0 0.67
Table 62. Affected areas in Naguilian, La Union during a 5-year rainfall return period - C
ARINGAY-BAUANG
BASIN Gusing Lioac Lioac Magungu- Mamat-Ing Mamat-Ing Nagsido-
Lioac Sur Natividad Ortiz
Sur Norte Norte nay Norte Sur risan
0-0.20 2.01 0.27 0.55 0.53 2.93 1.24 1.22 0.56 1.17 0.024
0.21-0.50 0.071 0.0081 0.012 0.022 0.18 0.028 0.039 0.046 0.13 0.013
0.51-1.00 0.06 0.0037 0.009 0.0098 0.17 0.016 0.025 0.018 0.052 0.0075
1.01-2.00 0.065 0.0033 0.018 0.015 0.17 0.023 0.029 0.017 0.048 0.017
(sq. km.)
Affected Area
2.01-5.00 0.2 0.0003 0.36 0.042 0.43 0.1 0.1 0.014 0.28 0.29
> 5.00 1.22 0 0.08 0.0076 0.18 0.69 1.7 0.00021 0.094 0.14
Table 63. Affected areas in Naguilian, La Union during a 5-year rainfall return period - D
ARINGAY-BAUANG
BASIN San Anto- Suguidan Suguidan
Ribsuan San Isidro Sili Tuddingan
130
nio Norte Sur

0-0.20 1.03 4.77 7.41 0.75 0.33 2.24 2.14
0.21-0.50 0.043 0.2 0.27 0.049 0.013 0.13 0.07
0.51-1.00 0.0052 0.095 0.15 0.028 0.0096 0.11 0.057
1.01-2.00 0.0049 0.067 0.1 0.007 0.0071 0.14 0.071
(sq. km.)
Affected Area
2.01-5.00 0.019 0.036 0.14 0.017 0.1 0.71 0.16
> 5.00 0.17 0.003 0.041 0.1 0.46 1.66 0.28
Figure 101. Affected areas in Naguilian, La Union during a 5-year rainfall return period - A
Figure 102. Affected areas in Naguilian, La Union during a 5-year rainfall return period - B
131
Figure 103. Affected areas in Naguilian, La Union during a 5-year rainfall return period - C
Figure 104. Affected areas in Naguilian, La Union during a 5-year rainfall return period - D
For the 5-year return period, 70.12% of the Municipality of Pugo, with an area of 60.54 square kilometers,
132
Table 64. Affected areas in Pugo, La Union during a 5-year rainfall return period
Affected Barangays in Pugo
ARINGAY-BAUANG Mao- Pobla-
BASIN Amban- Mao- Pobla- Tavora Tavora
Ambalite Cares Duplas asoas Palina cion San Luis
gonan asoas Sur cion East East Proper
Norte West
0-0.20 2.71 4.22 3.33 0.91 1.28 1.62 20.44 3.04 2.23 1.13 0.3 1.24
0.21-0.50 0.19 0.19 0.22 0.04 0.099 0.08 0.77 0.14 0.13 0.038 0.061 0.13
0.51-1.00 0.21 0.14 0.23 0.027 0.12 0.051 0.35 0.078 0.08 0.018 0.16 0.29
1.01-2.00 0.24 0.16 0.26 0.041 0.11 0.039 0.2 0.042 0.095 0.01 0.28 0.31
(sq. km.)
Affected Area
2.01-5.00 0.18 0.078 0.19 0.046 0.093 0.028 0.19 0.021 0.14 0.0012 0.51 0.26
> 5.00 0.0014 0 0.096 0.0001 0.084 0.051 0.047 0.031 0.041 0 0.024 0.0088
133
Figure 105. Affected areas in Pugo, La Union during a 5-year rainfall return period
134
For the 5-year return period, 6.50% of the Municipality of Rosario, with an area of 64.33 square kilometers,
Table 65. Affected areas in Rosario, La Union during a 5-year rainfall return period
Affected Barangays in Rosario
ARINGAY-BAUANG
BASIN Amban- Inabaan
Marcos Parasapas San Jose
gonan Norte
0-0.20 0.26 0.76 1.27 1.25 0.64
0.21-0.50 0.011 0.018 0.061 0.052 0.026
Affected Area
(sq. km.)
0.51-1.00 0.004 0.011 0.035 0.027 0.012

1.01-2.00 0.0033 0.0043 0.044 0.02 0.0059
2.01-5.00 0.0004 0.0006 0.0075 0.0054 0.003
> 5.00 0 0 0 0 0
Figure 106. Affected areas in Rosario, La Union during a 5-year rainfall return period
135
For the 5-year return period, 37.86% of the Municipality of San Fernando City, with an area of 121.05 square kilometers, will experience flood levels of less than 0.20
meters. 2.05% of the area will experience flood levels of 0.21 to 0.50 meters. Meanwhile, 1.19%, 1.04%, 1.15%, and 3.68% of the area will experience flood depths of 0.51
to 1 meter, 1.01 to 2 meters, 2.01 to 5 meters, and more than 5 meters, respectively. Listed in Tables 66-67 are the affected areas, in square kilometers, by flood depth per
barangay.
Table 66. Affected areas in San Fernando City, La Union during a 5-year rainfall return period - A
Affected Barangays in San Fernando City
ARINGAY-BAUANG
BASIN Bangban- Cabarsi- Cala- Nagyu- Namtut-
Apaleng Bacsil Birunget Bungro Cadaclan Langcuas Masicong
golan can bugao buyuban an
0-0.20 4.26 5.59 0.9 3.05 0.31 1.41 2.3 1.66 0.21 1.36 0.74 0.16
0.21-0.50 0.2 0.28 0.053 0.19 0.079 0.042 0.1 0.076 0.0098 0.041 0.029 0.0073
0.51-1.00 0.14 0.19 0.016 0.13 0.052 0.029 0.045 0.049 0.004 0.019 0.017 0.0022
1.01-2.00 0.14 0.14 0.003 0.061 0.064 0.047 0.026 0.023 0.0006 0.012 0.0069 0.0021
(sq. km.)
Affected Area
2.01-5.00 0.26 0.16 0.0005 0.0091 0.04 0.17 0.047 0.0081 0.0002 0.034 0.001 0
> 5.00 0.11 0.027 0 0 0 0.89 0.0001 0 0 1.56 0 0
136
Table 67. Affected areas in San Fernando City, La Union during a 5-year rainfall return period - B
ARINGAY-BAUANG
BASIN Narra Narra Pao Sibo- Tan-
Pacpaco Pao Sur Pias Puspus Sacyud Sagayad Sevilla
Este Oeste Norte an-Otong quigan
0-0.20 4 2.48 0.74 2.54 3.3 0.97 0.65 3.64 0.49 0.41 3.02 1.65
0.21-0.50 0.19 0.2 0.033 0.096 0.19 0.1 0.028 0.14 0.053 0.028 0.18 0.12
0.51-1.00 0.11 0.091 0.014 0.046 0.11 0.042 0.015 0.11 0.048 0.039 0.095 0.033
1.01-2.00 0.095 0.1 0.0068 0.014 0.065 0.0096 0.0063 0.098 0.055 0.18 0.082 0.015
(sq. km.)
Affected Area
2.01-5.00 0.05 0.11 0.0021 0.0022 0.1 0.0064 0.0006 0.22 0.048 0.08 0.039 0.0056
> 5.00 0 0.003 0 0 0.074 0.0001 0 1.78 0 0.0026 0 0
Figure 107. Affected areas in San Fernando City, La Union during a 5-year rainfall return period - A
Figure 108. Affected areas in San Fernando City, La Union during a 5-year rainfall return period - B
137
For the 5-year return period, 10.27% of the Municipality of San Gabriel, with an area of 154.19 square
Table 68. Affected areas in San Gabriel, La Union during a 5-year rainfall return period
ARINGAY-BAUANG Affected Barangays in San Gabriel

BASIN Amontoc Apayao Bayabas
0-0.20 3.35 12.07 0.42
0.21-0.50 0.072 0.49 0.000089
Affected Area
(sq. km.)
0.51-1.00 0.021 0.2 0.0001

1.01-2.00 0.011 0.12 0
2.01-5.00 0.01 0.13 0
> 5.00 0.0075 0.035 0
Figure 109. Affected areas in San Gabriel, La Union during a 5-year rainfall return period
138
For the 5-year return period, 8.19% of the Municipality of Santo Tomas, with an area of 58.53 square
levels of 0.21 to 0.50 meters. Meanwhile, 0.19%, 0.15%, and 0.05% of the area will experience flood
depths of 0.51 to 1 meter, 1.01 to 2 meters, and 2.01 to 5 meters, respectively. Listed in Table 69 are the
affected areas, in square kilometers, by flood depth per barangay.
Table 69. Affected areas in Santo Tomas, La Union during a 5-year rainfall return period
ARINGAY-BAUANG Affected Barangays in Santo Tomas

BASIN Ambitacay Bail Pongpong
0-0.20 0.36 0.14 4.29
0.21-0.50 0.011 0.008 0.17
Affected Area
(sq. km.)
0.51-1.00 0.0068 0.0029 0.099

1.01-2.00 0.0033 0.0013 0.085
2.01-5.00 0.0023 0 0.029
> 5.00 0 0 0.0015
Figure 110. Affected areas in Santo Tomas, La Union during a 5-year rainfall return period
139
For the 5-year return period, 52.15% of the Municipality of Tubao, with an area of 53.87 square kilometers, will experience flood levels of less than 0.20 meters. 2.74%
Table 70. Affected areas in Tubao, La Union during a 5-year rainfall return period - A
Affected Barangays in Tubao
ARINGAY-BAU-
ANG BASIN Francia Francia Halog Leones
Anduyan Caoigue Garcia Gonzales Halog East
Sur West West East
0-0.20 5.36 1.4 0.94 0.36 0.49 0.14 0.1 0.14 0
0.21-0.50 0.37 0.094 0.037 0.013 0.022 0.0042 0.0027 0.016 0
0.51-1.00 0.19 0.055 0.019 0.0062 0.012 0.0046 0.002 0.019 0
1.01-2.00 0.12 0.15 0.016 0.0061 0.012 0.0086 0.0037 0.048 0
(sq. km.)
Affected Area
2.01-5.00 0.95 0.13 0.052 0.015 0.052 0.028 0.021 0.088 0.000097
> 5.00 0.67 0 0.0072 0.015 0.55 0.89 0.74 0.58 0.93
140
Table 71. Affected areas in Tubao, La Union during a 5-year rainfall return period - B
ARINGAY-BAU-
ANG BASIN Leones Magsay- Santa
Linapew Lloren Pideg Poblacion Rizal
West say Teresa
0-0.20 1.02 2.01 2.74 0.38 2.49 0.28 4.68 5.59
0.21-0.50 0.069 0.065 0.23 0.012 0.11 0.016 0.21 0.21
0.51-1.00 0.054 0.022 0.12 0.0099 0.078 0.0085 0.13 0.12
1.01-2.00 0.084 0.0029 0.14 0.013 0.094 0.011 0.16 0.16
(sq. km.)
Affected Area
2.01-5.00 0.16 0.003 0.29 0.039 0.1 0.038 0.31 0.51
> 5.00 0.26 0.0001 0.23 0.82 0.0063 0.71 2.79 1.79
Figure 111. Affected areas in Tubao, La Union during a 5-year rainfall return period - A
Figure 112. Affected areas in Tubao, La Union during a 5-year rainfall return period - B
141
For the 25-year return period, 47.52% of the Municipality of Baguio City, with an area of 60.8979 square kilometers, will experience flood levels of less than 0.20 meters.
2.51% of the area will experience flood levels of 0.21 to 0.50 meters. Meanwhile, 1.39%, 1.19%, 1.18%, and 0.32% of the area will experience flood depths of 0.51 to 1
meter, 1.01 to 2 meters, 2.01 to 5 meters, and more than 5 meters, respectively. Listed in Tables 72-81 are the affected areas, in square kilometers, by flood depth per
barangay.
Affected Barangays in Baguio City
Aban-
ARINGAY-BAUANG A. Bonifa- ao-Zandu- Aurora Aurora
BASIN cio-Cagu- Alfonso Andres Aurora Hill Bakakeng Bakakeng
eta-Kay- Ambiong Asin Road Hill, North Hill, South
ioa-Riman- Tabora Bonifacio Proper Central North
ong-Chu- Central Central
do
gum-Otek
0-0.20 0.049 0.034 0.075 0.24 0.073 4.22 0.024 0.042 0.069 0.71 0.082
0.21-0.50 0.0067 0.0057 0.0028 0.009 0.005 0.24 0.0002 0.0026 0.0019 0.03 0.0048
0.51-1.00 0.00061 0.0041 0.0027 0.009 0.0004 0.11 0.00086 0.0018 0 0.025 0.00059
1.01-2.00 0.00053 0 0.0098 0.0027 0 0.083 0.0021 0.0038 0 0.02 0
(sq. km.)
142
Affected Area
2.01-5.00 0.00053 0 0.011 0 0 0.065 0.0032 0.00097 0 0.012 0
> 5.00 0 0 0.0038 0 0 0.0036 0 0 0 0.0055 0
ARINGAY-BAUANG Cabinet Camdas
BASIN Bal-Marco- Bayan Park Bayan Park Bgh Com- Brook- Camp
Balsigan Brookside Hill-Teach- Subdivi- Camp 8
ville East West pound spoint Allen
er’s Camp sion
0-0.20 0.12 0.054 0.17 0.075 0.017 0.091 0.098 0.35 0.16 0.087 0.13
0.21-0.50 0.0025 0.00018 0.0071 0.0014 0.00049 0.0042 0.0075 0.01 0.0022 0.0058 0.0055
0.51-1.00 0.0054 0.00023 0.0069 0.0025 0.0016 0.0023 0.0047 0.0053 0.0012 0 0.008
1.01-2.00 0.0051 0.000012 0.0052 0.0014 0.0011 0 0.0055 0.0092 0.00032 0 0.0046
(sq. km.)
Affected Area
2.01-5.00 0.0007 0 0.00035 0 0 0 0.0031 0.017 0.00091 0 0.0012
> 5.00 0 0 0 0 0 0 0 0.00021 0.00067 0 0
ARINGAY-BAUANG Country Dominican
BASIN Campo City Camp City Camp Cresencia Dizon Sub- Engineers’ Fairview
Club Vil- Hill-Mira- Dontogan Dps Area
Filipino Central Proper Village division Hill Village
lage dor
0-0.20 0.12 0.069 0.083 0.19 0.16 0.25 1.18 0.55 0.12 0.1 0.25
0.21-0.50 0.014 0.0024 0.0028 0.016 0.0047 0.0044 0.045 0.016 0.0015 0.0021 0.017
0.51-1.00 0.021 0.0007 0 0.0068 0.0072 0.0013 0.026 0.0077 0.0019 0.0001 0.021
1.01-2.00 0.043 0 0 0.007 0.013 0.0017 0.024 0.0077 0.006 0 0.02
(sq. km.)
Affected Area
2.01-5.00 0.019 0 0 0.0013 0.0079 0.0011 0.063 0.0093 0.022 0 0.0067
> 5.00 0 0 0 0 0 0 0.077 0.00041 0.0017 0 0
143
ARINGAY-BAUANG General Gener- General
BASIN Greenwa- Guisad Guisad Happy Happy Holy Ghost
Ferdinand Emilio F. al Luna, Luna, Up- Gibraltar
ter Village Central Sorong Hollow Homes Extension
Aguinaldo Lower per
0-0.20 0.091 0.062 0.027 0.034 0.2 0.13 0.077 0.29 0.74 0.19 0.29
0.21-0.50 0.0085 0.0029 0.0016 0.0016 0.0035 0.0048 0.0005 0.015 0.061 0.00093 0.0079
0.51-1.00 0.0057 0.0013 0.0026 0.0002 0.00034 0.0012 0.0002 0.0064 0.02 0.0014 0.0064
1.01-2.00 0.0041 0.0023 0.00094 0 0.000071 0.00029 0 0.002 0.02 0.0026 0.0081
(sq. km.)
Affected Area
2.01-5.00 0.00025 0.0049 0.0032 0 0 0.000026 0 0.00039 0.0097 0.011 0.021
> 5.00 0 0 0.0011 0 0 0 0 0 0.0015 0.007 0.000015
ARINGAY-BAUANG Legar-
BASIN Holy Ghost Honey- Imelda R. Imelda Kabayani- Kagitin- Kayang Kayang-Hill- Lopez
Irisan da-Burn-
Proper moon Marcos Village han gan Extension top Jaena
ham-Kisad
0-0.20 0.043 0.071 0.024 0.074 5.61 0.061 0.031 0.089 0.0099 0.35 0.047
0.21-0.50 0.0029 0.005 0.00064 0.0084 0.28 0.0021 0.00073 0.0045 0.00048 0.024 0.0016
0.51-1.00 0.0013 0.0026 0 0.0044 0.17 0.0041 0.0015 0.0013 0.000098 0.028 0.001
1.01-2.00 0.0025 0.0008 0 0.0069 0.13 0.0037 0.0022 0.00022 0 0.037 0
(sq. km.)
Affected Area
2.01-5.00 0.01 0.0002 0 0.0048 0.12 0 0.0083 0 0 0.059 0
> 5.00 0.0024 0 0 0 0.052 0 0.000055 0 0 0.0091 0
144
ARINGAY-BAUANG Lourdes
Lourdes Lourdes Magsaysay Malcolm Middle
BASIN Subdi- Magsaysay, Magsaysay, Manuel A. Military
Subdivision Subdivi- Lualhati Private Square-Per- Quezon Hill
vision, Lower Upper Roxas Cut-Off
Extension sion, Lower Road fecto Subdivision
Proper
0-0.20 0.027 0.03 0.078 0.32 0.022 0.089 0.0086 0.045 0.095 0.14 0.28
0.21-0.50 0.0039 0.0033 0.005 0.019 0.0055 0.014 0.00062 0.00016 0.0037 0.012 0.021
0.51-1.00 0.0028 0.0011 0 0.017 0.0011 0.0021 0.0008 0.000023 0.0036 0.0012 0.0062
1.01-2.00 0 0 0 0.013 0.002 0.0049 0.0036 0 0.011 0 0.0018
(sq. km.)
Affected Area
2.01-5.00 0 0 0 0 0.0076 0.015 0.0036 0 0.0095 0 0.00018
> 5.00 0 0 0 0 0.0026 0.00065 0 0 0 0 0
ARINGAY-BAUANG
BASIN Modern Mrr-Queen New Luc- Outlook Padre Bur- Padre Palma-Ur- Pinsao Pi-
Pacdal Phil-Am Pinget
Site, West Of Peace ban Drive gos Zamora bano lot Project
0-0.20 0.095 0.061 0.11 0.14 1.08 0.11 0.024 0.076 0.032 0.35 0.091
0.21-0.50 0.0074 0.0029 0.0062 0.0046 0.048 0.0048 0.0024 0.0011 0.0047 0.019 0.0043
0.51-1.00 0.00045 0.0026 0.0067 0.0027 0.028 0.0051 0.0044 0 0.0027 0.0071 0.00088
1.01-2.00 0 0.00072 0.0038 0.0023 0.0067 0.01 0.015 0 0.0034 0.0065 0.0013
(sq. km.)
Affected Area
2.01-5.00 0 0 0.016 0 0.001 0.013 0.046 0 0 0.0076 0.0015
> 5.00 0 0 0.0033 0 0 0.0036 0.0025 0 0 0.0053 0.00033
ARINGAY-BAUANG Quiri-
Rizal
145
BASIN Pinsao Quezon Hill Quezon Quirino Quirino Quirino Quirino no-Mag-
Poliwes Pucsusan Monument
Proper Proper Hill, Upper Hill, East Hill, Lower Hill, Middle Hill, West saysay,
Area
Upper
0-0.20 2.13 0.062 0.14 0.17 0.32 0.052 0.033 0.05353 0.054 0.099 0.072
0.21-0.50 0.11 0.0031 0.013 0.016 0.024 0.0022 0.0028 0.00999 0.01 0.0025 0.004
0.51-1.00 0.053 0.0006 0.0043 0.017 0.012 0.0035 0.0029 0.005461 0.0055 0.0039 0.0018
1.01-2.00 0.027 0 0.0003 0.014 0.0024 0.0027 0.0006 0.0015 0.0015 0.0054 0.0001
(sq. km.)
Affected Area
2.01-5.00 0.012 0 0 0.0008 0.00081 0.013 0.0009 0.0002 0.0002 0.021 0
> 5.00 0.0029 0 0 0 0 0.00069 0 0 0 0.0065 0
ARINGAY-BAUANG Saint Sanitary Sanitary
BASIN Rock Quar- Rock Quar- Rock Quar- San Anto- San Luis San Roque San Vicen-
Joseph Salud Mitra Camp, Camp,
ry, Lower ry, Middle ry, Upper nio Village Village Village te
Village North South
0-0.20 0.036 0.055 0.061 0.63 0.15 0.13 0.61 0.055 0.2 0.15 0.082
0.21-0.50 0.0017 0.0032 0.011 0.03 0.0061 0.0042 0.031 0.004 0.012 0.0053 0.0072
0.51-1.00 0.0018 0.0012 0.0044 0.027 0.0005 0.00067 0.011 0 0.0058 0.00082 0.003
1.01-2.00 0.00085 0.0018 0.00022 0.017 0 0 0.0081 0 0.0038 0 0.0001
(sq. km.)
Affected Area
2.01-5.00 0.00062 0.001 0 0.0061 0 0 0.0022 0 0.0037 0 0
> 5.00 0 0 0 0 0 0 0 0 0 0 0
Table 81. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - J
146
ARINGAY-BAUANG
Slu-Svp
BASIN Santa Esco- Santo Rosa- Santo Tomas Session Slaughter Teodora Victoria
Housing South Drive Trancoville
lastica rio Proper Road Area House Area Alonzo Village
Village
0-0.20 0.077 0.19 0.21 0.2 0.024 0.2 0.33 0.056 0.22 0.18
0.21-0.50 0.006 0.01 0.01 0.0054 0.0005 0.0052 0.038 0.0063 0.011 0.013
0.51-1.00 0.0012 0.0035 0.0025 0.01 0.0006 0.0033 0.017 0.0012 0.0069 0.0039
1.01-2.00 0.003 0.00069 0.0014 0.014 0.0011 0.003 0.014 0.0051 0.0089 0.0012
(sq. km.)
Affected Area
2.01-5.00 0.00067 0 0 0.0014 0.0051 0.0025 0.0062 0.0082 0.008 0.0001
> 5.00 0 0 0 0 0.00085 0 0 0 0.0025 0
147
148
149
150
Figure 122. Affected areas in Baguio City, Benguet during a 25-year rainfall return period - J
151
Affected Barangays in Kapangan
ARIN-
GAY-BAU- Pobla-
Gaswel- Pon-
ANG BASIN Datakan Gadang Labueg Paykek cion Sagubo
ing gayan
Central
0-0.20 3.81 0.19 15.19 1.98 0.036 1.83 9.7 8.03
0.21-
0.16 0.0024 0.45 0.11 0.0018 0.077 0.42 0.38

0.50
0.51-
0.071 0 0.17 0.081 0 0.05 0.21 0.17
1.00
1.01-
0.046 0 0.14 0.078 0.000037 0.058 0.14 0.068
2.00
2.01-
0.054 0 0.16 0.1 0.000038 0.066 0.11 0.056
5.00
> 5.00 0.13 0 0.21 0.037 0 0.0052 0.014 0.0013
152
For the 25-year return period, 71.44% of the Municipality of La Trinidad, with an area of 74.908 square kilometers, will experience flood levels of less than 0.20 meters.
barangay.
ARINGAY-BAUANG Affected Barangays in La Trinidad

BASIN Alapang Alno Ambiong Bahong Balili Beckel Betag Bineng
0-0.20 3.97 9.68 1.55 4.38 0.25 1.45 0.14 10.9
0.21-0.50 0.25 0.41 0.063 0.34 0.076 0.067 0.056 0.59
0.51-1.00 0.11 0.18 0.038 0.16 0.089 0.029 0.036 0.18
1.01-2.00 0.11 0.11 0.052 0.1 0.26 0.015 0.26 0.089
(sq. km.)
Affected Area
2.01-5.00 0.12 0.16 0.025 0.073 0.28 0.0041 0.93 0.047
> 5.00 0.015 0.79 0.0028 0.013 0.017 0 0.14 0.14
153

BASIN Cruz Lubas Pico Poblacion Puguis Shilan Tawang Wangal
0-0.20 0.46 1.39 0.96 1.79 3.5 6.55 2.24 4.31
0.21-0.50 0.018 0.084 0.27 0.089 0.19 0.26 0.11 0.17
0.51-1.00 0.0016 0.04 0.31 0.064 0.08 0.14 0.022 0.13
1.01-2.00 0.0006 0.035 0.3 0.051 0.034 0.097 0.012 0.095
(sq. km.)
Affected Area
2.01-5.00 0 0.038 0.25 0.17 0.015 0.069 0.0064 0.099
> 5.00 0 0.0009 0.032 0 0 0.018 0.0005 0.027
154
For the 25-year return period, 91.04% of the Municipality of Sablan, with an area of 90.22 square kilometers,
ARIN- Affected Barangays in Sablan

GAY-BAU- Banan- Baneng- Pobla-
ANG BASIN Bagong Balluay Bayabas Kamog Pappa
gan beng cion
0-0.20 4.82 8.48 17.48 12.72 12.49 10.72 10.25 5.18
0.21-
0.1 0.36 0.67 0.48 0.57 0.34 0.34 0.23

0.50
0.51-
0.05 0.18 0.29 0.24 0.23 0.17 0.18 0.11
1.00
1.01-
0.036 0.12 0.22 0.13 0.17 0.13 0.13 0.092
2.00
2.01-
0.077 0.059 0.22 0.15 0.13 0.077 0.19 0.068
5.00
> 5.00 0.87 0.1 0.052 0.42 0.024 0.056 0.039 0.004
155
For the 25-year return period, 38.63% of the Municipality of Tuba, with an area of 322.02 square kilometers,
ARIN- Affected Barangays in Tuba

GAY-BAU- Nan- Pobla- San Pas- Tabaan Tabaan Tadian- Taloy Taloy Twin
ANG BASIN galisan cion cual Norte Sur gan Norte Sur Peaks
0-0.20 10.49 20.11 6.3 16.59 22.91 11.74 14.48 19.16 2.6
0.21-
0.47 1.07 0.26 0.69 1.03 0.59 0.61 0.78 0.14

0.50
0.51-
0.23 0.51 0.11 0.28 0.46 0.26 0.26 0.31 0.057
1.00
1.01-
0.18 0.25 0.085 0.13 0.26 0.17 0.16 0.22 0.0025
2.00
2.01-
0.35 0.3 0.11 0.14 0.27 0.21 0.15 0.24 0.0005
5.00
> 5.00 0.81 0.12 0.021 0.18 0.32 0.17 0.094 0.14 0
156
For the 25-year return period, 50.22% of the Municipality of Tublay, with an area of 63.21 square kilometers,
ARIN- Affected Barangays in Tublay

GAY-BAU- Ambas- Ambong- Tublay
ANG BASIN sador Ba-Ayan Basil Caponga Daclan Tuel
dolan Central
0-0.20 0.023 3.65 2 3.77 2.86 4.86 3.45 11.13
0.21-
0.0001 0.16 0.097 0.16 0.23 0.21 0.17 0.5

0.50
0.51-
0 0.08 0.028 0.074 0.085 0.087 0.094 0.21
1.00
1.01-
0 0.095 0.012 0.052 0.04 0.058 0.13 0.14
2.00
2.01-
0 0.05 0.0049 0.044 0.015 0.0098 0.12 0.16
5.00
> 5.00 0 0.018 0 0.15 0 0 0.039 0.18
157
ARINGAY-BAUANG
Norte East Norte
0-0.20 1.5 0.016 0.5 0.022 0.0093 0.037
0.21-0.50 0.083 0.00046 0.05 0.022 0.002 0.028
Affected Area
(sq. km.)
0.51-1.00 0.049 0 0.027 0.019 0.0029 0.048

1.01-2.00 0.037 0 0.026 0.0065 0.0071 0.032
2.01-5.00 0.017 0 0.05 0 0.0065 0.0095
> 5.00 0 0 0.0017 0 0 0
158
For the 25-year return period, 64.69% of the Municipality of Aringay, with an area of 95.65 square kilometers, will experience flood levels of less than 0.20 meters. 4.57%
BASIN Macaba- Pobla- San An- San Beni-
to cion tonio to Norte
East West
0-0.20 0.35 5.68 0.64 21.76 5.54 7.43 1.91 1.19 0.57 0.032 0.62 1.65
0.21-0.50 0.14 0.28 0.46 1 0.24 0.3 0.097 0.16 0.031 0.052 0.063 0.22
0.51-1.00 0.082 0.17 0.75 0.53 0.14 0.18 0.043 0.09 0.017 0.29 0.051 0.044
1.01-2.00 0.4 0.11 1.42 0.34 0.12 0.13 0.045 0.12 0.025 1.7 0.15 0.031
(sq. km.)
Affected Area
2.01-5.00 1.01 0.12 0.21 0.38 0.19 0.11 0.19 0.3 0.073 1.95 0.63 0.028
> 5.00 0 0.26 0 2.1 0.048 0.0032 1.92 0.087 0.53 0.075 1.58 0.0003
159
Table 90. Affected areas in Aringay, La Union during a 25-year rainfall return period - B
East West East West
0-0.20 0.54 5.1 1.91 1.8 1.21 1.59 0.21 0.14 1.59 0.26 0.14 0.0012
0.21-0.50 0.041 0.34 0.21 0.096 0.068 0.14 0.013 0.15 0.073 0.068 0.13 0.013
0.51-1.00 0.029 0.14 0.16 0.087 0.049 0.11 0.0099 0.091 0.048 0.012 0.035 0.035
1.01-2.00 0.035 0.1 0.25 0.086 0.041 0.11 0.017 0.66 0.028 0.0096 0.054 0.11
(sq. km.)
Affected Area
2.01-5.00 0.083 0.084 0.57 0.46 0.039 0.11 0.059 1.47 0.017 0.0069 0.05 1.21
> 5.00 0.16 0.0003 0.18 1 0.11 0.0062 2.09 0.3 0.00088 0 0 0
160
For the 25-year return period, 72.84% of the Municipality of Bagulin, with an area of 77.97 square
ARIN- Affected Barangays in Bagulin

GAY-BAU- Alibang- Cam- Tagud- Tio-An- Wal-
ANG BASIN Baay Dagup Libbo Suyo
say baly tud gan layan
0-0.20 10.97 1.81 3.06 4.53 4.95 4.7 7.41 2.74 13.3 3.32
0.21-
0.35 0.079 0.13 0.15 0.22 0.17 0.39 0.12 0.46 0.14
0.50
0.51-
0.17 0.015 0.076 0.079 0.13 0.056 0.2 0.07 0.19 0.07
1.00
1.01-
0.12 0.006 0.051 0.03 0.097 0.03 0.14 0.035 0.13 0.039
2.00
2.01-
0.12 0.015 0.032 0.037 0.1 0.03 0.097 0.029 0.27 0.051
5.00
> 5.00 1.03 0.36 0.066 0.69 0.14 0.49 0.16 0.49 0.69
161
ARINGAY-BAUANG
BASIN Baccuit Cabalayan-
Acao Baccuit Sur Bagbag Ballay Bawanta Boy-Utan Bucayab Cabisilan
Norte gan
0-0.20 2.76 0.1 0.0062 0.44 0.25 0.75 3.08 2.9 4.55 1.73
0.21-0.50 0.16 0.1 0.023 0.17 0.0033 0.041 0.32 0.18 0.3 0.16
0.51-1.00 0.12 0.19 0.11 0.096 0.0065 0.021 0.19 0.1 0.2 0.086
1.01-2.00 0.21 0.44 0.22 0.13 0.015 0.0089 0.12 0.13 0.31 0.047
(sq. km.)
Affected Area
2.01-5.00 0.66 0.12 0.27 0.063 0.11 0.0022 0.15 0.17 0.32 0.095
> 5.00 3.37 0 0 0.0017 1.91 0 0.83 0.028 0.94 0.045
162
ARINGAY-BAUANG
BASIN Central Lower San
Calumbaya Carmay Casilagan Central East Dili Disso-Or Guerrero Nagrebcan
West Agustin
0-0.20 0.94 0.0037 5.66 0.28 0.0035 0.34 0.45 0.45 2.15 0.17
0.21-0.50 0.18 0.019 0.46 0.025 0.0011 0.056 0.12 0.16 0.12 0.027
0.51-1.00 0.17 0.073 0.32 0.013 0.0015 0.061 0.064 0.022 0.065 0.02
1.01-2.00 0.12 0.1 0.34 0.02 0.0061 0.046 0.037 0.0054 0.057 0.017
(sq. km.)
Affected Area
2.01-5.00 0.034 0.11 0.35 0.19 0.21 0.0005 0.027 0.0004 0.12 0.14
> 5.00 0.00038 0.0017 0.034 0.38 0.047 0 0 0 0.096 0.52
ARINGAY-BAUANG Palug- Payocpoc
BASIN Pagdalagan Parian Payocpoc Payocpoc
Palintucang si-Limman- Parian Este Paringao Norte Pilar
Sur Oeste Norte Este Sur
sangan Oeste
0-0.20 0 0.44 4.99 0.28 0.072 0.0029 1.51 0.4 0.37 0.028
0.21-0.50 0 0.061 0.35 0.22 0.078 0.02 0.2 0.47 0.17 0.025
0.51-1.00 0 0.066 0.23 0.6 0.29 0.062 0.18 0.55 0.12 0.025
1.01-2.00 0.000003 0.1 0.15 0.15 1.18 0.25 0.11 0.35 0.093 0.00037
(sq. km.)
Affected Area
2.01-5.00 0 0.068 0.05 0.28 0.24 0.17 0.13 0.021 0.0049 0
> 5.00 0 0.0043 0.0003 0.068 0 0 0.0097 0 0 0
ARINGAY-BAUANG
163
BASIN Santa Mon- Upper San
Pottot Pudoc Pugo Quinavite Santiago Taberna Urayong
ica Agustin
0-0.20 0.5 0.01 0 0.5 5.68 1.09 0.0056 7.84 0.39
0.21-0.50 0.12 0.027 0.0003 0.14 0.3 0.28 0.0013 0.63 0.18
0.51-1.00 0.13 0.11 0.009 0.05 0.2 0.12 0.02 0.41 0.12
1.01-2.00 0.055 0.22 0.15 0.025 0.16 0.07 0.13 0.36 0.037
(sq. km.)
Affected Area
2.01-5.00 0.0002 0.64 1.21 0.07 0.11 0.013 1.03 0.47 0.002
> 5.00 0 0 0.064 0 0.0053 0 0.0015 0.013 0
164
Figure 136 - Affected areas in Bauang, La Union during a 25-year rainfall return period - D
165
For the 25-year return period, 89.12% of the Municipality of Burgos, with an area of 51.92 square kilometers, will experience flood levels of less than 0.20 meters. 3.74% of
ARINGAY-BAUANG
0-0.20 2.98 4.64 0.36 2.24 11.51 3.82 4.64 3.97 1.06 3.3 1.31 6.44
0.21-0.50 0.14 0.18 0.0083 0.08 0.45 0.19 0.23 0.094 0.041 0.2 0.067 0.26
0.51-1.00 0.067 0.058 0.013 0.046 0.22 0.093 0.082 0.056 0.022 0.037 0.025 0.1
1.01-2.00 0.036 0.045 0.0084 0.029 0.13 0.066 0.054 0.037 0.013 0.032 0.019 0.071
(sq. km.)
Affected Area
2.01-5.00 0.036 0.046 0.01 0.067 0.09 0.067 0.057 0.041 0.054 0.028 0.044 0.13
> 5.00 0.1 0.053 0.076 0.41 0.015 0.099 0.064 0.55 0.022 0.0091 0.097 0.17
166
167
For the 25-year return period, 69.01% of the Municipality of Caba, with an area of 56.19 square kilometers, will experience flood levels of less than 0.20 meters. 7.78% of
to 2 meters, 2.01 to 5 meters, and more than 5 meters, respectively. Listed in Tables 97-98 are the affected areas, in square kilometers, by flood depth per barangay.
ARINGAY-BAUANG
BASIN Juan Poblacion Pobla- San Cor-
Bautista Gana Las-Ud Liquicia San Carlos
Cartas Norte cion Sur nelio
0-0.20 2.47 0.89 2.35 0.097 12.94 0.33 0.31 0.2 4.94
0.21-0.50 0.44 0.76 0.19 0.045 0.61 0.088 0.14 0.3 0.4
0.51-1.00 0.2 0.3 0.11 0.042 0.42 0.073 0.042 0.31 0.29
1.01-2.00 0.086 0.11 0.048 0.047 0.41 0.19 0.025 0.41 0.26
(sq. km.)
Affected Area
2.01-5.00 0.041 0.095 0.024 0.012 0.74 0.25 0.016 0.082 0.31
> 5.00 0.001 0 0.0011 0 1.13 0.02 0 0 0.23
168
ARINGAY-BAUANG
BASIN San Fer- San Gre- Santiago Santiago
San Jose Sobredillo Urayong Wenceslao
min gorio Norte Sur
0-0.20 4.7 1.8 2.11 0.015 0.02 3.71 1.79 0.11
0.21-0.50 0.34 0.12 0.24 0.073 0.048 0.18 0.26 0.13
0.51-1.00 0.24 0.11 0.15 0.17 0.24 0.1 0.12 0.14
1.01-2.00 0.19 0.12 0.14 0.29 0.84 0.067 0.022 0.025
(sq. km.)
Affected Area
2.01-5.00 0.29 0.18 0.56 0.11 1.74 0.035 0.0052 0
> 5.00 0.41 0.15 0.2 0 0 0.029 0 0
169
ARINGAY-BAUANG
BASIN Al-Alinao Al-Alinao Ambaracao Ambaracao Baraoas Baraoas
Norte Sur Norte Sur Norte Sur
0-0.20 0.47 1.92 3.01 2.19 1.64 0.23 3.52 0.38 0.000051 0.041
0.21-0.50 0.016 0.15 0.13 0.18 0.14 0.013 0.17 0.0057 0.000043 0.00091
0.51-1.00 0.012 0.14 0.054 0.24 0.14 0.0049 0.1 0.006 0.00016 0.0023
1.01-2.00 0.019 0.23 0.028 0.17 0.062 0.0022 0.092 0.01 0.00074 0.0036
(sq. km.)
Affected Area
2.01-5.00 0.062 0.32 0.047 0.41 0.079 0.00065 0.19 0.11 0.25 0.092
> 5.00 0.68 1.04 0.3 1.97 0.044 0.000085 0.21 1.45 0.4 0.34
170
ARINGAY-BAUANG
BASIN Bimmoto- Cabaritan Cabaritan Dal-Li- Daramuan- Gusing
bot Norte Sur paoen gan Norte
0-0.20 1.03 1.58 4.05 0.00051 0.03 5.35 0.0073 0.091 0.42 3.54
0.21-0.50 0.034 0.086 0.19 0 0.0008 0.23 0.0002 0.0055 0.027 0.14
0.51-1.00 0.015 0.074 0.11 0 0.0013 0.14 0.0001 0.0083 0.02 0.076
1.01-2.00 0.0072 0.058 0.072 0.00035 0.0044 0.11 0.00067 0.016 0.0043 0.064
(sq. km.)
Affected Area
2.01-5.00 0.01 0.03 0.069 0.0045 0.046 0.19 0.0013 0.033 0.0001 0.074
> 5.00 0.027 0.017 0.03 0.29 0.26 0.37 0.47 0.74 0 0.73
ARINGAY-BAUANG
BASIN Magungu- Mamat-Ing Mamat-Ing Nagsido-
Gusing Sur Imelda Lioac Norte Lioac Sur Natividad Ortiz
nay Norte Sur risan
0-0.20 1.9 0.26 0.5 0.48 2.74 1.15 1.15 0.53 1.07 0.00019
0.21-0.50 0.079 0.011 0.013 0.031 0.18 0.03 0.05 0.043 0.12 0.0016
0.51-1.00 0.047 0.0054 0.0073 0.013 0.13 0.015 0.022 0.028 0.065 0.0037
1.01-2.00 0.064 0.0043 0.0099 0.015 0.14 0.016 0.023 0.016 0.046 0.0088
(sq. km.)
Affected Area
2.01-5.00 0.1 0.0006 0.042 0.035 0.45 0.041 0.065 0.028 0.1 0.045
> 5.00 1.43 0 0.46 0.054 0.43 0.84 1.81 0.013 0.37 0.43
ARINGAY-BAUANG
171
nio Norte Sur
0-0.20 0.99 4.64 7.23 0.7 0.31 1.99 2.05
0.21-0.50 0.052 0.24 0.31 0.057 0.014 0.12 0.08
0.51-1.00 0.014 0.13 0.18 0.04 0.008 0.065 0.048
1.01-2.00 0.0052 0.083 0.12 0.012 0.0072 0.1 0.063
(sq. km.)
Affected Area
2.01-5.00 0.013 0.071 0.16 0.017 0.021 0.37 0.13
> 5.00 0.19 0.0055 0.093 0.13 0.55 2.34 0.4
172
173
For the 25-year return period, 67.28% of the Municipality of Pugo, with an area of 60.54 square kilometers, will experience flood levels of less than 0.20 meters. 3.67% of
Norte West
0-0.20 2.56 4.09 3.1 0.88 1.14 1.57 19.93 2.95 2.09 1.1 0.23 1.08
0.21-0.50 0.15 0.21 0.23 0.041 0.1 0.084 0.94 0.16 0.15 0.05 0.015 0.094
0.51-1.00 0.11 0.13 0.22 0.032 0.098 0.069 0.47 0.095 0.082 0.025 0.069 0.11
1.01-2.00 0.27 0.15 0.25 0.028 0.16 0.047 0.26 0.073 0.12 0.015 0.24 0.33
(sq. km.)
Affected Area
2.01-5.00 0.4 0.2 0.35 0.073 0.16 0.033 0.26 0.032 0.18 0.0031 0.67 0.56
> 5.00 0.03 0.011 0.16 0.0068 0.13 0.074 0.12 0.041 0.098 0 0.12 0.059
174
175
For the 25-year return period, 6.31% of the Municipality of Rosario, with an area of 64.33 square kilometers,
ARINGAY-BAUANG
BASIN Ambango- Inabaan
nan Norte
0-0.20 0.25 0.75 1.22 1.22 0.62
0.21-0.50 0.014 0.022 0.066 0.067 0.034
Affected Area
(sq. km.)
0.51-1.00 0.0059 0.014 0.051 0.035 0.016

1.01-2.00 0.0045 0.0083 0.041 0.023 0.01
2.01-5.00 0.0011 0.0018 0.036 0.018 0.0052
> 5.00 0 0 0 0 0
176
to 1 meter, 1.01 to 2 meters, 2.01 to 5 meters, and more than 5 meters, respectively. Listed in Tables 105-106 are the affected areas, in square kilometers, by flood depth
per barangay.
ARINGAY-BAU-
ANG BASIN Bangban- Cabarsi- Nagyu-
Apaleng Bacsil Birunget Bungro Cadaclan Calabugao Langcuas Masicong Namtutan
golan can buyuban
0-0.20 4.04 5.37 0.87 2.94 0.2 1.31 2.22 1.61 0.21 1.31 0.73 0.15
0.21-0.50 0.22 0.32 0.068 0.19 0.11 0.047 0.12 0.085 0.012 0.05 0.035 0.011
0.51-1.00 0.15 0.22 0.028 0.17 0.067 0.024 0.069 0.063 0.005 0.027 0.021 0.0029
1.01-2.00 0.16 0.19 0.0075 0.1 0.087 0.032 0.041 0.035 0.0013 0.014 0.011 0.0028
(sq. km.)
Affected Area
2.01-5.00 0.3 0.21 0.0011 0.04 0.074 0.11 0.058 0.019 0.0002 0.025 0.0022 0
> 5.00 0.25 0.089 0 0.0006 0.0003 1.07 0.011 0.0004 0 1.6 0 0
177
ARINGAY-BAUANG
BASIN Narra Sibo-
Narra Este Pacpaco Pao Norte Pao Sur Pias Puspus Sacyud Sagayad Sevilla Tanquigan
Oeste an-Otong
0-0.20 3.87 2.3 0.72 2.49 3.12 0.92 0.64 3.34 0.43 0.37 2.91 1.57
0.21-0.50 0.21 0.26 0.041 0.11 0.22 0.1 0.031 0.15 0.052 0.03 0.21 0.16
0.51-1.00 0.15 0.13 0.018 0.066 0.12 0.074 0.021 0.1 0.057 0.022 0.12 0.05
1.01-2.00 0.11 0.12 0.012 0.024 0.087 0.017 0.011 0.081 0.081 0.085 0.098 0.03
(sq. km.)
Affected Area
2.01-5.00 0.1 0.15 0.0044 0.0071 0.12 0.0098 0.0009 0.17 0.076 0.21 0.079 0.011
> 5.00 0.00014 0.02 0.0001 0 0.17 0.0013 0 2.13 0.0006 0.015 0 0
178

0-0.20 3.29 11.73 0.42
0.21-0.50 0.11 0.63 0.0002
Affected Area
(sq. km.)
0.51-1.00 0.03 0.27 0.0001

1.01-2.00 0.02 0.16 0
2.01-5.00 0.015 0.16 0
> 5.00 0.014 0.1 0
179
Table 108. Affected areas in Santo Tomas, La Union during a 25-year rainfall return period

0-0.20 0.36 0.14 4.17
0.21-0.50 0.011 0.0097 0.2
Affected Area
(sq. km.)
0.51-1.00 0.0081 0.0042 0.12

1.01-2.00 0.0055 0.0024 0.11
2.01-5.00 0.0032 0 0.065
> 5.00 0.0002 0 0.0062
Figure 149. Affected areas in Santo Tomas, La Union during a 25-year rainfall return period
180
ARINGAY-BAUANG
BASIN Francia Francia Halog Halog Leones
Anduyan Caoigue Garcia Gonzales
Sur West East West East
0-0.20 5.07 1.3 0.87 0.34 0.46 0.12 0.081 0.033 0
0.21-0.50 0.43 0.086 0.044 0.015 0.021 0.0033 0.0036 0.0052 0
0.51-1.00 0.24 0.057 0.023 0.0092 0.012 0.0033 0.0065 0.0083 0
1.01-2.00 0.15 0.076 0.023 0.007 0.0079 0.0052 0.0083 0.017 0
(sq. km.)
Affected Area
2.01-5.00 0.27 0.29 0.052 0.015 0.018 0.022 0.012 0.08 0
> 5.00 1.49 0.015 0.057 0.029 0.61 0.92 0.76 0.74 0.93
181
ARINGAY-BAUANG
BASIN Leones Magsay- Pobla- Santa
Linapew Lloren Pideg Rizal
West say cion Teresa
0-0.20 0.79 1.96 2.52 0.34 2.4 0.26 4.35 5.23
0.21-0.50 0.055 0.093 0.24 0.016 0.12 0.013 0.23 0.24
0.51-1.00 0.034 0.035 0.12 0.0085 0.093 0.01 0.14 0.12
1.01-2.00 0.051 0.009 0.11 0.012 0.085 0.0077 0.13 0.095
(sq. km.)
Affected Area
2.01-5.00 0.2 0.0027 0.27 0.033 0.15 0.031 0.3 0.19
> 5.00 0.51 0.0015 0.48 0.85 0.039 0.75 3.12 2.51
182
For the 100-year return period, 46.61% of the Municipality of Baguio City, with an area of 60.8979 square kilometers, will experience flood levels of less than 0.20 meters.
barangay.
Aban-
ARINGAY-BAU- A. Bonifa-
ao-Zandu- Aurora Aurora
ANG BASIN cio-Cagu- Alfonso Andres Aurora Hill Bakakeng Bakakeng
eta-Kay- Ambiong Asin Road Hill, North Hill, South
ioa-Riman- Tabora Bonifacio Proper Central North
ong-Chu- Central Central
do
gum-Otek
0-0.20 0.046 0.031 0.072 0.24 0.072 4.13 0.023 0.04 0.068 0.7 0.081
0.21-0.50 0.0093 0.0065 0.0027 0.01 0.0055 0.28 0.0008 0.0035 0.0025 0.033 0.0061
0.51-1.00 0.001 0.0059 0.002 0.008 0.0007 0.12 0.0004 0.0014 0.0001 0.025 0.00087
1.01-2.00 0.00052 0.00046 0.0077 0.0052 0 0.089 0.0017 0.0033 0 0.024 0
(sq. km.)
183
Affected Area
2.01-5.00 0.00094 0 0.014 0 0 0.087 0.0043 0.0026 0 0.015 0
> 5.00 0 0 0.0062 0 0 0.0094 0.0004 0 0 0.0076 0
ARINGAY-BAU- Cabinet Camdas
ANG BASIN Bal-Mar- Bayan Park Bayan Park Bayan Park Bgh Com- Brook- Camp
Balsigan Brookside Hill-Teach- Subdivi- Camp 8
coville East Village West pound spoint Allen
er’s Camp sion
0-0.20 0.12 0.054 0.17 0.075 0.016 0.089 0.095 0.34 0.16 0.084 0.12 0.124508
0.21-0.50 0.0021 0.00017 0.0075 0.0016 0.00042 0.0049 0.0089 0.016 0.0028 0.0085 0.0055 0.005502
0.51-1.00 0.0049 0.00018 0.0056 0.0023 0.00094 0.0029 0.0044 0.0057 0.0013 0.0005 0.0064 0.006351
1.01-2.00 0.0069 0.00011 0.0073 0.002 0.0022 0.00021 0.0062 0.0078 0.00073 0 0.0091 0.00911
(sq. km.)
Affected Area
2.01-5.00 0.00082 0 0.0013 0 0 0 0.0044 0.023 0.00098 0 0.0029 0.002882
> 5.00 0 0 0 0 0 0 0 0.00031 0.00091 0 0.000038 0.000038
ARINGAY-BAU- City Country
ANG BASIN Campo City Camp Cresencia Dizon Sub- Dominican Engineers’ Fairview
Camp Club Dontogan Dps Area
Filipino Proper Village division Hill-Mirador Hill Village
Central Village
0-0.20 0.11 0.068 0.082 0.19 0.16 0.25 1.15 0.54 0.12 0.1 0.24
0.21-0.50 0.015 0.0023 0.0033 0.017 0.0045 0.0056 0.055 0.022 0.0028 0.0033 0.017
0.51-1.00 0.018 0.0015 0 0.0088 0.0059 0.0015 0.029 0.0097 0.0016 0.0001 0.019
1.01-2.00 0.034 0 0 0.008 0.012 0.0019 0.023 0.0077 0.0049 0 0.023
(sq. km.)
Affected Area
2.01-5.00 0.038 0 0 0.0019 0.013 0.0011 0.058 0.01 0.021 0 0.011
> 5.00 0 0 0 0 0 0 0.1 0.00051 0.0055 0 0
184
ARINGAY-BAU-
General General
ANG BASIN General Greenwater Guisad Guisad Happy Hol- Happy Holy Ghost
Ferdinand Emilio F. Luna, Gibraltar
Luna, Lower Village Central Sorong low Homes Extension
Aguinaldo Upper
0-0.20 0.088 0.06 0.026 0.033 0.2 0.13 0.077 0.29 0.71 0.18 0.29
0.21-0.50 0.0079 0.0039 0.0013 0.0017 0.0038 0.0061 0.0006 0.017 0.076 0.0022 0.0081
0.51-1.00 0.0059 0.0013 0.0028 0.00024 0.00051 0.0013 0.0002 0.0087 0.028 0.0018 0.0061
1.01-2.00 0.0056 0.0027 0.0011 0 0.000071 0.00039 0 0.0029 0.02 0.0019 0.0074
(sq. km.)
Affected Area
2.01-5.00 0.0014 0.0058 0.003 0 0 0.000037 0 0.00049 0.014 0.0094 0.024
> 5.00 0 0 0.0023 0 0 0 0 0 0.0024 0.012 0.0017
ARINGAY-BAUANG Holy Legar-
BASIN Honey- Imelda R. Imelda Kabayani- Kayang Kayang-Hill- Lopez
Ghost Irisan Kagitingan da-Burn-
moon Marcos Village han Extension top Jaena
Proper ham-Kisad
0-0.20 0.04 0.07 0.023 0.072 5.5 0.06 0.029 0.089 0.0097 0.34 0.046
0.21-0.50 0.0034 0.0058 0.0012 0.0083 0.32 0.0017 0.00089 0.0046 0.0005 0.022 0.002
0.51-1.00 0.0019 0.0027 0 0.0041 0.19 0.0032 0.0014 0.0016 0.00028 0.027 0.0013
1.01-2.00 0.0015 0.0013 0 0.0067 0.15 0.0054 0.0022 0.00032 0.000001 0.037 0
(sq. km.)
Affected Area
2.01-5.00 0.0079 0.0002 0 0.0072 0.15 0.0003 0.009 0 0 0.073 0
> 5.00 0.0071 0 0 0 0.07 0 0.0012 0 0 0.012 0
185
ARINGAY-BAUANG Lourdes
Lourdes Lourdes Magsaysay Malcolm Middle
BASIN Subdi- Magsaysay, Magsaysay, Manuel A. Military
Subdivision, Subdivision, Lualhati Private Square-Per- Quezon Hill
vision Lower Upper Roxas Cut-Off
Lower Proper Road fecto Subdivision
Extension
0-0.20 0.026 0.029 0.076 0.31 0.02 0.084 0.0082 0.044 0.091 0.14 0.28
0.21-0.50 0.0029 0.0032 0.007 0.021 0.0062 0.014 0.0002 0.00013 0.0054 0.013 0.023
0.51-1.00 0.0043 0.002 0.0001 0.015 0.0009 0.0035 0.00082 0.00019 0.0028 0.0024 0.0078
1.01-2.00 0 0 0 0.017 0.0019 0.0032 0.0031 0 0.0098 0 0.0024
(sq. km.)
Affected Area
2.01-5.00 0 0 0 0.0018 0.0071 0.018 0.0048 0 0.014 0 0.00018
> 5.00 0 0 0 0 0.0053 0.0024 0 0 0 0 0
ARINGAY-BAU-
ANG BASIN Modern Mrr-Queen New Luc- Outlook Padre Bur- Padre Palma-Ur- Pinsao Pilot
Pacdal Phil-Am Pinget
Site, West Of Peace ban Drive gos Zamora bano Project
0-0.20 0.093 0.06 0.11 0.14 1.07 0.099 0.019 0.075 0.031 0.34 0.089
0.21-0.50 0.0086 0.0033 0.0059 0.0053 0.055 0.0058 0.0018 0.0019 0.0042 0.024 0.0057
0.51-1.00 0.0013 0.0027 0.0075 0.003 0.03 0.0034 0.003 0 0.0033 0.0088 0.00087
1.01-2.00 0.00015 0.00092 0.0037 0.0029 0.012 0.0082 0.008 0 0.0041 0.0061 0.00098
(sq. km.)
Affected Area
2.01-5.00 0 0.0001 0.017 0 0.0012 0.02 0.048 0 0.0001 0.0095 0.0019
> 5.00 0 0 0.0057 0 0 0.0071 0.015 0 0 0.006 0.00033
ARINGAY-BAUANG Quiri-
Pinsao Quezon Hill Quezon Hill, Quirino Hill, Quirino Hill, Quirino Hill, Quirino Hill, Rizal Monu-
186
BASIN
Poliwes Pucsusan no-Magsay-

Proper Proper Upper East Lower Middle West ment Area
say, Upper
0-0.20 2.11 0.061 0.14 0.16 0.31 0.05 0.032 0.05 0.095 0.07 0.034
0.21-0.50 0.12 0.0038 0.017 0.015 0.026 0.0019 0.003 0.011 0.0036 0.0046 0.0016
0.51-1.00 0.062 0.0009 0.005 0.018 0.013 0.0034 0.0024 0.0071 0.0036 0.0024 0.00077
1.01-2.00 0.033 0 0.0008 0.018 0.0043 0.0027 0.0013 0.0018 0.0064 0.0005 0.0014
(sq. km.)
Affected Area
2.01-5.00 0.015 0 0 0.0017 0.0011 0.012 0.0009 0.0003 0.022 0 0.0032
> 5.00 0.0034 0 0 0 0 0.0041 0 0 0.0079 0 0.000006
ARINGAY-BAUANG San Sanitary Sanitary
BASIN Rock Quar- Rock Quar- Rock Quar- Saint Jo- San Antonio San Luis San Vicen-
Salud Mitra Roque Camp, Camp,
ry, Lower ry, Middle ry, Upper seph Village Village Village te
Village North South
0-0.20 0.051 0.058 0.63 0.14 0.13 0.6 0.053 0.19 0.15 0.08 0.076
0.21-0.50 0.0029 0.0098 0.031 0.0075 0.0043 0.034 0.0061 0.015 0.0069 0.0083 0.0072
0.51-1.00 0.0015 0.007 0.028 0.00096 0.0011 0.014 0 0.006 0.00087 0.0039 0.0014
1.01-2.00 0.0025 0.0013 0.021 0 0 0.0086 0 0.0045 0.000052 0.0002 0.0029
(sq. km.)
Affected Area
2.01-5.00 0.0041 0.00029 0.0092 0 0 0.0043 0 0.004 0 0 0.0011
> 5.00 0.000014 0 0 0 0 0 0 0 0 0 0
Table 120. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - J
187
ARINGAY-BAUANG Slu-Svp
BASIN Santa Esco- Santo Ro- Santo To- Session Slaughter Teodora Victoria Santa Es-
Housing South Drive Trancoville
lastica sario mas Proper Road Area House Area Alonzo Village colastica
Village
0-0.20 0.076 0.19 0.21 0.2 0.023 0.2 0.32 0.054 0.22 0.18 1.673953
0.21-0.50 0.0072 0.012 0.011 0.0077 0.00031 0.006 0.034 0.006 0.012 0.012 0.10784
0.51-1.00 0.0014 0.0042 0.0038 0.0089 0.0004 0.0038 0.027 0.0023 0.0079 0.0067 0.065944
1.01-2.00 0.0029 0.0012 0.0021 0.016 0.001 0.0031 0.014 0.0033 0.0066 0.0011 0.05135
(sq. km.)
Affected Area
2.01-5.00 0.0011 0 0 0.0035 0.0049 0.003 0.0087 0.011 0.012 0.0003 0.044057
> 5.00 0 0 0 0 0.0024 0 0 0 0.0037 0 0.006074
188
189
190
191
Figure 161. Affected areas in Baguio City, Benguet during a 100-year rainfall return period - J
192
Listed in Tables 121 are the affected areas, in square kilometers, by flood depth per barangay.
Affected Barangays in Kapangan
ARIN-
GAY-BAUANG Pobla-
Gaswel-
BASIN Datakan Gadang Labueg Paykek cion Pongayan Sagubo
ing
Central
0-0.20 3.75 0.19 15.01 1.93 0.035 1.8 9.57 7.91
0.21-
0.18 0.0037 0.53 0.11 0.0023 0.086 0.45 0.43
0.50
Affected Area
0.51-
0.08 0 0.19 0.08 0.0002 0.056 0.25 0.19
(sq. km.)
1.00
1.01-
0.057 0 0.14 0.076 0 0.055 0.16 0.087
2.00
2.01-
0.056 0 0.18 0.12 0.000075 0.081 0.14 0.07
5.00
> 5.00 0.16 0 0.25 0.06 0 0.013 0.029 0.0038
193
For the 100-year return period, 69.90% of the Municipality of La Trinidad, with an area of 74.908 square kilometers, will experience flood levels of less than 0.20 meters.
barangay.

BASIN Alapang Alno Ambiong Bahong Balili Beckel Betag Bineng
0-0.20 3.87 9.44 1.53 4.26 0.22 1.43 0.1 10.7
0.21-0.50 0.29 0.45 0.07 0.38 0.05 0.076 0.048 0.69
0.51-1.00 0.12 0.21 0.04 0.19 0.07 0.035 0.027 0.23
1.01-2.00 0.1 0.13 0.048 0.11 0.18 0.015 0.06 0.1
(sq. km.)
Affected Area
2.01-5.00 0.16 0.16 0.04 0.1 0.41 0.0079 0.93 0.066
> 5.00 0.029 0.95 0.0058 0.022 0.041 0 0.39 0.17
194

BASIN Cruz Lubas Pico Poblacion Puguis Shilan Tawang Wangal
0-0.20 0.45 1.35 0.9 1.75 3.46 6.45 2.19 4.25
0.21-0.50 0.026 0.096 0.19 0.092 0.2 0.31 0.14 0.18
0.51-1.00 0.0027 0.049 0.28 0.068 0.096 0.14 0.026 0.13
1.01-2.00 0.0007 0.036 0.36 0.052 0.041 0.12 0.015 0.1
(sq. km.)
Affected Area
2.01-5.00 0 0.044 0.34 0.17 0.022 0.087 0.008 0.11
> 5.00 0 0.0054 0.053 0.026 0 0.031 0.0005 0.049
195
For the 100-year return period, 89.56% of the Municipality of Sablan, with an area of 90.22 square
ARIN- Affected Barangays in Sablan

GAY-BAUANG Baneng- Pobla-
BASIN Bagong Balluay Banangan Bayabas Kamog Pappa
beng cion
0-0.20 4.71 8.35 17.23 12.48 12.27 10.58 10.1 5.09
0.21-
0.11 0.41 0.8 0.55 0.68 0.39 0.4 0.25

0.50
0.51-
0.052 0.21 0.32 0.27 0.26 0.18 0.2 0.13
1.00
1.01-
0.041 0.13 0.24 0.15 0.2 0.14 0.14 0.11
2.00
2.01-
0.079 0.074 0.26 0.17 0.17 0.11 0.21 0.092
5.00
> 5.00 0.95 0.12 0.093 0.53 0.042 0.071 0.076 0.0091
196
For the 100-year return period, 37.96% of the Municipality of Tuba, with an area of 322.02 square
ARIN- Affected Barangays in Tuba

GAY-BAU- Nangal- Pobla- San Tabaan Tabaan Tadian- Taloy Twin
ANG BASIN Taloy Sur
isan cion Pascual Norte Sur gan Norte Peaks
0-0.20 10.24 19.73 6.21 16.34 22.53 11.5 14.24 18.88 2.56
0.21-
0.49 1.22 0.3 0.78 1.18 0.68 0.72 0.88 0.16

0.50
0.51-
0.25 0.59 0.12 0.36 0.52 0.3 0.29 0.36 0.069
1.00
1.01-
0.19 0.31 0.09 0.16 0.31 0.19 0.19 0.25 0.0071
2.00
2.01-
0.31 0.31 0.13 0.15 0.29 0.23 0.19 0.27 0.0006
5.00
> 5.00 1.05 0.2 0.036 0.23 0.43 0.25 0.13 0.21 0
197
For the 100-year return period, 49.35% of the Municipality of Tublay, with an area of 63.21 square
ARIN- Affected Barangays in Tublay
GAY-BAU-
ANG Ambas- Ambong- Tublay
Ba-Ayan Basil Caponga Daclan Tuel
BASIN sador dolan Central
0-0.20 0.023 3.59 1.97 3.69 2.79 4.8 3.38 10.94
0.21-
0.0001 0.18 0.12 0.18 0.26 0.23 0.19 0.58

0.50
0.51-
0 0.087 0.034 0.082 0.098 0.099 0.093 0.24
1.00
1.01-
0 0.098 0.014 0.047 0.051 0.07 0.11 0.15
2.00
2.01-
0 0.075 0.0079 0.07 0.022 0.019 0.16 0.18
5.00
> 5.00 0 0.021 0 0.18 0.0008 0.0001 0.066 0.23
198
ARINGAY-BAUANG
Norte East Norte
0-0.20 1.48 0.015 0.48 0.017 0.0081 0.023
0.21-0.50 0.092 0.00068 0.053 0.021 0.00084 0.03
Affected Area
(sq. km.)
0.51-1.00 0.055 0 0.033 0.021 0.0027 0.041

1.01-2.00 0.041 0 0.027 0.0096 0.0074 0.05
2.01-5.00 0.026 0 0.058 0 0.0087 0.011
> 5.00 0.0001 0 0.0045 0 0 0
199
For the 100-year return period, 62.16% of the Municipality of Aringay, with an area of 95.65 square kilometers, will experience flood levels of less than 0.20 meters. 4.80%
BASIN Macaba- Pobla- San An- San Beni-
to cion tonio to Norte
East West
0-0.20 0.28 5.52 0.51 21.16 5.44 7.33 1.79 0.91 0.55 0.021 0.55 1.58
0.21-0.50 0.17 0.3 0.39 1.08 0.26 0.34 0.1 0.14 0.035 0.024 0.064 0.27
0.51-1.00 0.088 0.18 0.69 0.58 0.14 0.19 0.046 0.1 0.017 0.17 0.035 0.055
1.01-2.00 0.3 0.13 1.39 0.37 0.13 0.14 0.042 0.12 0.027 1.4 0.063 0.032
(sq. km.)
Affected Area
2.01-5.00 1.14 0.12 0.49 0.4 0.19 0.15 0.1 0.3 0.071 2.38 0.55 0.037
> 5.00 0 0.36 0 2.51 0.09 0.011 2.12 0.37 0.55 0.11 1.84 0.0018
200
Table 129. Affected areas in Aringay, La Union during a 100-year rainfall return period - B
East West East West
0-0.20 0.51 4.98 1.73 1.73 1.16 1.54 0.17 0.095 1.57 0.22 0.096 0
0.21-0.50 0.042 0.38 0.2 0.094 0.069 0.14 0.013 0.15 0.076 0.094 0.13 0.0041
0.51-1.00 0.031 0.17 0.15 0.074 0.055 0.12 0.01 0.077 0.055 0.019 0.059 0.032
1.01-2.00 0.038 0.12 0.24 0.081 0.049 0.1 0.013 0.34 0.034 0.0086 0.051 0.082
(sq. km.)
Affected Area
2.01-5.00 0.092 0.11 0.57 0.21 0.064 0.15 0.045 1.82 0.023 0.0097 0.068 1.25
> 5.00 0.18 0.0031 0.38 1.34 0.12 0.012 2.14 0.33 0.0019 0 0 0
201
For the 100-year return period, 71.15% of the Municipality of Bagulin, with an area of 77.97 square kilometers, will experience flood levels of less than 0.20 meters. 3.26%
of the area will experience flood levels of 0.21 to 0.50 meters. Meanwhile, 1.54%, 0.99%, 1.14%, and 6.14% of the area will experience flood depths of 0.51 to 1 meter, 1.01
Affected Barangays in Bagulin
ARINGAY-BAU-
ANG BASIN Alibang- Cam- Tio-An- Wal-
Baay Cardiz Dagup Libbo Suyo Tagudtud
say baly gan layan
0-0.20 10.73 1.76 3.01 4.48 4.8 4.59 7.18 2.69 13.05 3.2
0.21-0.50 0.41 0.092 0.15 0.18 0.24 0.19 0.44 0.14 0.54 0.16
0.51-1.00 0.18 0.022 0.082 0.092 0.15 0.071 0.23 0.08 0.22 0.081
1.01-2.00 0.14 0.0068 0.06 0.04 0.11 0.036 0.16 0.039 0.14 0.044
(sq. km.)
Affected Area
2.01-5.00 0.15 0.015 0.045 0.044 0.11 0.039 0.13 0.037 0.27 0.056
> 5.00 1.17 0.39 0.077 0.0032 0.79 0.2 0.59 0.18 0.63 0.77
202
203
ARINGAY-BAUANG
BASIN Baccuit Cabalayan-
Acao Baccuit Sur Bagbag Ballay Bawanta Boy-Utan Bucayab Cabisilan
Norte gan
0-0.20 2.6 0.075 0.0028 0.37 0.24 0.74 2.94 2.83 4.44 1.67
0.21-0.50 0.16 0.098 0.0099 0.18 0.0019 0.045 0.33 0.19 0.29 0.16
0.51-1.00 0.12 0.15 0.078 0.11 0.0024 0.024 0.21 0.11 0.2 0.11
1.01-2.00 0.17 0.47 0.23 0.15 0.011 0.013 0.13 0.13 0.28 0.054
(sq. km.)
Affected Area
2.01-5.00 0.68 0.16 0.3 0.083 0.055 0.003 0.16 0.21 0.4 0.086
> 5.00 3.56 0 0 0.0028 1.99 0.000075 0.92 0.038 1.02 0.086
204
ARINGAY-BAUANG
BASIN Central Lower San
Calumbaya Carmay Casilagan Central East Dili Disso-Or Guerrero Nagrebcan
West Agustin
0-0.20 0.89 0.00059 5.5 0.26 0.0022 0.32 0.4 0.41 2.06 0.14
0.21-0.50 0.17 0.0051 0.47 0.029 0.00088 0.065 0.14 0.18 0.14 0.033
0.51-1.00 0.2 0.046 0.34 0.013 0.0012 0.06 0.08 0.036 0.069 0.029
1.01-2.00 0.11 0.13 0.35 0.021 0.004 0.058 0.042 0.0073 0.066 0.013
(sq. km.)
Affected Area
2.01-5.00 0.061 0.11 0.44 0.1 0.18 0.0018 0.035 0.0006 0.12 0.091
> 5.00 0.0098 0.0058 0.069 0.48 0.078 0 0 0 0.16 0.59
ARINGAY-BAUANG Palug- Payocpoc
BASIN Pagdalagan Parian Payocpoc Payocpoc
Palintucang si-Limman- Parian Este Paringao Norte Pilar
Sur Oeste Norte Este Sur
sangan Oeste
0-0.20 0 0.41 4.88 0.19 0.051 0.0004 1.37 0.26 0.33 0.021
0.21-0.50 0 0.068 0.38 0.18 0.049 0.007 0.26 0.48 0.17 0.025
0.51-1.00 0 0.05 0.24 0.61 0.18 0.047 0.19 0.55 0.15 0.031
1.01-2.00 0.000003 0.13 0.19 0.26 1.01 0.2 0.15 0.44 0.094 0.00078
(sq. km.)
Affected Area
2.01-5.00 0 0.076 0.075 0.23 0.58 0.25 0.15 0.072 0.017 0
> 5.00 0 0.008 0.0004 0.13 0 0 0.019 0 0 0
ARINGAY-BAUANG
Santa Mon- Upper San
205
BASIN Pottot Pudoc Pugo Quinavite Santiago Taberna Urayong
ica Agustin
0-0.20 0.48 0.0001 0 0.45 5.58 0.99 0.0043 7.58 0.35
0.21-0.50 0.1 0.0096 0 0.15 0.31 0.33 0.0006 0.65 0.17
0.51-1.00 0.13 0.063 0.0031 0.065 0.22 0.15 0.01 0.45 0.15
1.01-2.00 0.092 0.22 0.095 0.026 0.17 0.075 0.085 0.38 0.052
(sq. km.)
Affected Area
2.01-5.00 0.0004 0.71 1.22 0.079 0.15 0.03 1.09 0.57 0.0065
> 5.00 0 0.00035 0.12 0 0.019 0 0.0076 0.11 0
206
Figure 175. Affected areas in Bauang, La Union during a 100-year rainfall return period - D
207
For the 100-year return period, 87.39% of the Municipality of Burgos, with an area of 51.92 square kilometers, will experience flood levels of less than 0.20 meters. 4.44%
of the area will experience flood levels of 0.21 to 0.50 meters. Meanwhile, 1.80%, 1.19%, 1.43%, and 3.76% of the area will experience flood depths of 0.51 to 1 meter, 1.01
ARINGAY-BAUANG
0-0.20 2.92 4.57 0.35 2.17 11.34 3.74 4.54 3.9 1.04 3.21 1.28 6.32
0.21-0.50 0.16 0.22 0.007 0.093 0.52 0.22 0.28 0.097 0.048 0.25 0.078 0.31
0.51-1.00 0.074 0.071 0.011 0.053 0.25 0.1 0.1 0.059 0.025 0.052 0.03 0.11
1.01-2.00 0.046 0.047 0.012 0.034 0.16 0.078 0.062 0.04 0.015 0.033 0.017 0.08
(sq. km.)
Affected Area
2.01-5.00 0.041 0.052 0.012 0.07 0.12 0.078 0.064 0.038 0.047 0.037 0.045 0.13
> 5.00 0.12 0.069 0.086 0.45 0.037 0.11 0.086 0.61 0.041 0.013 0.12 0.21
208
209
For the 100-year return period, 66.27% of the Municipality of Caba, with an area of 56.19 square kilometers, will experience flood levels of less than 0.20 meters. 7.67% of
to 2 meters, 2.01 to 5 meters, and more than 5 meters, respectively. Listed in Tables 136-137 are the affected areas, in square kilometers, by flood depth per barangay.
ARINGAY-BAUANG
BASIN Juan Car- Poblacion Poblacion San Cor-
Bautista Gana Las-Ud Liquicia San Carlos
tas Norte Sur nelio
0-0.20 2.36 0.73 2.3 0.062 12.63 0.27 0.27 0.13 4.79
0.21-0.50 0.45 0.77 0.2 0.044 0.65 0.054 0.14 0.28 0.4
0.51-1.00 0.26 0.4 0.13 0.059 0.43 0.073 0.073 0.32 0.29
1.01-2.00 0.11 0.12 0.068 0.063 0.45 0.19 0.035 0.47 0.28
(sq. km.)
Affected Area
2.01-5.00 0.057 0.13 0.03 0.014 0.74 0.31 0.02 0.12 0.38
> 5.00 0.0028 0 0.002 0 1.37 0.046 0 0 0.29
210
ARINGAY-BAUANG
BASIN San Fer- San Gre- Santiago Santiago Wences-
San Jose Sobredillo Urayong
min gorio Norte Sur lao
0-0.20 4.54 1.73 1.96 0.0022 0.007 3.64 1.73 0.078
0.21-0.50 0.34 0.12 0.21 0.023 0.032 0.21 0.27 0.12
0.51-1.00 0.25 0.1 0.17 0.15 0.18 0.11 0.16 0.16
1.01-2.00 0.2 0.11 0.21 0.33 0.74 0.081 0.033 0.043
(sq. km.)
Affected Area
2.01-5.00 0.33 0.22 0.5 0.15 1.93 0.045 0.007 0.0002
> 5.00 0.51 0.2 0.35 0 0 0.038 0 0
PLACEHOLDER - Affected areas in Caba, La Union during a 100-year rainfall return period
PLACEHOLDER - Affected areas in Caba, La Union during a 100-year rainfall return period
211
ARINGAY-BAUANG
BASIN Al-Alinao Al-Alinao Ambaracao Ambaracao Baraoas Baraoas
Norte Sur Norte Sur Norte Sur
0-0.20 0.43 1.8 2.94 2.08 1.59 0.23 3.43 0.36 0 0.015
0.21-0.50 0.016 0.11 0.15 0.15 0.13 0.012 0.18 0.005 0 0.0089
0.51-1.00 0.013 0.11 0.066 0.21 0.13 0.0058 0.11 0.0035 0 0.0072
1.01-2.00 0.019 0.21 0.033 0.19 0.089 0.0032 0.089 0.014 0 0.007
(sq. km.)
Affected Area
2.01-5.00 0.054 0.31 0.038 0.26 0.092 0.0033 0.2 0.029 0.0019 0.012
> 5.00 0.73 1.27 0.34 2.27 0.065 0.00017 0.28 1.54 0.66 0.43
212
ARINGAY-BAUANG
BASIN Bimmoto- Cabaritan Cabaritan Dal-Li- Daramuan- Gusing
bot Norte Sur paoen gan Norte
0-0.20 1.02 1.55 3.96 0.00024 0.027 5.23 0.0067 0.02 0.42 3.48
0.21-0.50 0.039 0.095 0.21 0.000076 0.00052 0.25 0.0001 0.0072 0.026 0.15
0.51-1.00 0.015 0.065 0.12 0.0002 0.0015 0.16 0.0004 0.024 0.023 0.084
1.01-2.00 0.011 0.077 0.086 0 0.002 0.12 0.0002 0.028 0.0075 0.068
(sq. km.)
Affected Area
2.01-5.00 0.0099 0.034 0.075 0.002 0.016 0.17 0.0016 0.048 0.0035 0.078
> 5.00 0.034 0.028 0.06 0.29 0.29 0.47 0.47 0.77 0 0.76
ARINGAY-BAUANG
BASIN Magungu- Mamat-Ing Mamat-Ing Nagsido-
Gusing Sur Imelda Lioac Norte Lioac Sur Natividad Ortiz
nay Norte Sur risan
0-0.20 1.84 0.25 0.47 0.45 2.6 1.09 1.1 0.5 1 0
0.21-0.50 0.089 0.014 0.013 0.036 0.17 0.033 0.054 0.038 0.11 0
0.51-1.00 0.052 0.0059 0.0093 0.0091 0.14 0.016 0.025 0.026 0.076 0
1.01-2.00 0.054 0.004 0.0096 0.013 0.15 0.013 0.024 0.022 0.055 0.00019
(sq. km.)
Affected Area
2.01-5.00 0.1 0.0019 0.031 0.039 0.4 0.029 0.049 0.041 0.097 0.029
> 5.00 1.49 0 0.49 0.079 0.61 0.92 1.86 0.028 0.44 0.46
ARINGAY-BAUANG
nio Norte Sur
213
0-0.20 0.97 4.55 7.11 0.67 0.3 1.88 2
0.21-0.50 0.056 0.27 0.35 0.054 0.015 0.12 0.09
0.51-1.00 0.022 0.15 0.2 0.045 0.0094 0.057 0.048
1.01-2.00 0.0049 0.1 0.14 0.024 0.0067 0.077 0.062
(sq. km.)
Affected Area
2.01-5.00 0.013 0.092 0.17 0.018 0.019 0.3 0.12
> 5.00 0.2 0.0099 0.14 0.14 0.56 2.56 0.46
214
215
For the 100-year return period, 65.67% of the Municipality of Pugo, with an area of 60.54 square kilometers, will experience flood levels of less than 0.20 meters. 3.87% of
Norte West
0-0.20 2.49 4.03 2.99 0.86 1.07 1.53 19.61 2.89 1.97 1.08 0.2 1.02
0.21-0.50 0.17 0.22 0.2 0.044 0.11 0.084 1.05 0.17 0.16 0.058 0.0062 0.075
0.51-1.00 0.098 0.14 0.19 0.032 0.087 0.077 0.54 0.1 0.095 0.028 0.0075 0.055
1.01-2.00 0.17 0.13 0.28 0.031 0.12 0.054 0.31 0.08 0.11 0.019 0.013 0.15
(sq. km.)
Affected Area
2.01-5.00 0.54 0.25 0.44 0.077 0.25 0.043 0.31 0.049 0.24 0.0055 0.27 0.68
> 5.00 0.063 0.025 0.21 0.015 0.15 0.085 0.18 0.054 0.15 0 0.84 0.25
216
217
For the 100-year return period, 6.20% of the Municipality of Rosario, with an area of 64.33 square
levels of 0.21 to 0.50 meters. Meanwhile, 0.20%, 0.16%, and 0.13% of the area will experience flood
depths of 0.51 to 1 meter, 1.01 to 2 meters, and 2.01 to 5 meters, respectively. Listed in Table 143 are the
affected areas, in square kilometers, by flood depth per barangay.
ARINGAY-BAUANG
BASIN Ambango- Inabaan
nan Norte
0-0.20 0.25 0.74 1.2 1.19 0.61
0.21-0.50 0.017 0.027 0.071 0.076 0.036
Affected Area
(sq. km.)
0.51-1.00 0.0065 0.016 0.053 0.037 0.019

1.01-2.00 0.0052 0.01 0.048 0.028 0.012
2.01-5.00 0.0021 0.0028 0.05 0.024 0.0071
> 5.00 0 0 0 0 0
218
to 1 meter, 1.01 to 2 meters, 2.01 to 5 meters, and more than 5 meters, respectively. Listed in Tables 144-145 are the affected areas, in square kilometers, by flood depth
per barangay.
ARINGAY-BAU-
ANG BASIN Bangban- Cabarsi-
Apaleng Bacsil Birunget Bungro Cadaclan Calabugao Langcuas Masicong Nagyubuyuban Namtutan
golan can
0-0.20 3.91 5.2 0.85 2.87 0.17 1.25 2.17 1.58 0.2 1.28 0.72 0.15
0.21-0.50 0.24 0.35 0.073 0.19 0.084 0.051 0.13 0.093 0.012 0.057 0.038 0.012
0.51-1.00 0.15 0.22 0.037 0.17 0.089 0.027 0.076 0.069 0.0057 0.032 0.022 0.0045
1.01-2.00 0.16 0.21 0.011 0.12 0.09 0.027 0.056 0.046 0.0024 0.016 0.015 0.0031
(sq. km.)
Affected Area
2.01-5.00 0.3 0.26 0.0022 0.087 0.11 0.084 0.063 0.026 0.0002 0.023 0.0034 0
> 5.00 0.37 0.13 0 0.0054 0.002 1.15 0.024 0.0015 0 1.62 0.00018 0
219
ARINGAY-BAUANG
BASIN Narra Sibo-
Narra Este Pacpaco Pao Norte Pao Sur Pias Puspus Sacyud Sagayad Sevilla Tanquigan
Oeste an-Otong
0-0.20 3.8 2.2 0.7 2.45 3 0.9 0.62 3.08 0.39 0.35 2.83 1.52
0.21-0.50 0.22 0.29 0.047 0.13 0.24 0.097 0.034 0.15 0.056 0.028 0.24 0.19
0.51-1.00 0.16 0.15 0.022 0.075 0.13 0.087 0.023 0.094 0.053 0.021 0.14 0.063
1.01-2.00 0.13 0.13 0.014 0.034 0.078 0.027 0.015 0.079 0.092 0.048 0.1 0.035
(sq. km.)
Affected Area
2.01-5.00 0.13 0.18 0.0065 0.013 0.14 0.013 0.0015 0.16 0.1 0.26 0.11 0.021
> 5.00 0.0023 0.034 0.0002 0 0.26 0.0022 0 2.42 0.0029 0.026 0.0005 0
220

0-0.20 3.25 11.51 0.42
0.21-0.50 0.13 0.72 0.0003
Affected Area
(sq. km.)
0.51-1.00 0.034 0.31 0.0001

1.01-2.00 0.018 0.17 0
2.01-5.00 0.021 0.18 0
> 5.00 0.019 0.14 0
221
Table 147. Affected areas in Santo Tomas, La Union during a 100-Year rainfall return period

0-0.20 0.35 0.14 4.1
0.21-0.50 0.0098 0.011 0.22
Affected Area
(sq. km.)
0.51-1.00 0.01 0.0056 0.13

1.01-2.00 0.0062 0.0028 0.12
2.01-5.00 0.0042 0.0001 0.096
> 5.00 0.0003 0 0.0097
Figure 188. Affected areas in Santo Tomas, La Union during a 100-Year rainfall return period
222
ARINGAY-BAUANG
BASIN Francia Francia Halog Halog Leones
Anduyan Caoigue Garcia Gonzales
Sur West East West East
0-0.20 4.84 1.26 0.84 0.33 0.45 0.12 0.049 0.0005 0
0.21-0.50 0.46 0.083 0.048 0.016 0.022 0.004 0.0021 0.0033 0
0.51-1.00 0.29 0.062 0.024 0.0094 0.013 0.003 0.0039 0.0045 0
1.01-2.00 0.17 0.062 0.026 0.0074 0.0088 0.0034 0.01 0.011 0
(sq. km.)
Affected Area
2.01-5.00 0.14 0.28 0.057 0.017 0.014 0.018 0.039 0.036 0
> 5.00 1.76 0.084 0.076 0.036 0.63 0.93 0.77 0.83 0.93
223
ARINGAY-BAUANG
BASIN Leones Magsay- Pobla- Santa
Linapew Lloren Pideg Rizal
West say cion Teresa
0-0.20 0.67 1.93 2.41 0.33 2.35 0.24 4.18 5.03
0.21-0.50 0.044 0.11 0.25 0.017 0.12 0.016 0.25 0.25
0.51-1.00 0.029 0.042 0.15 0.0096 0.099 0.0098 0.14 0.13
1.01-2.00 0.044 0.015 0.11 0.012 0.083 0.009 0.14 0.088
(sq. km.)
Affected Area
2.01-5.00 0.13 0.0026 0.22 0.03 0.17 0.026 0.29 0.18
> 5.00 0.72 0.0023 0.62 0.87 0.064 0.76 3.27 2.7
224
Among the barangays in the Municipality of Baguio City in Benguet, Irisan is projected to have the highest
percentage of area that will experience flood levels, at 10.46%. Meanwhile, Asin Road posted the second
highest percentage of area that may be affected by flood depths, at 7.75%.
Among the barangays in the Municipality of Kapangan in Benguet, Gasweling is projected to have the
highest percentage of area that will experience flood levels, at 12.18%. Meanwhile, Pongayan posted the
second highest percentage of area that may be affected by flood depths, at 7.92%.
Among the barangays in the Municipality of La Trinidad in Benguet, Bineng is projected to have the highest
percentage of area that will experience flood levels, at 15.96%. Meanwhile, Alno posted the second highest
percentage of area that may be affected by flood depths, at 15.12%.
Among the barangays in the Municipality of Sablan in Benguet, Banangan is projected to have the highest
percentage of area that will experience flood levels, at 20.98%. Meanwhile, Banengbeng posted the second
Among the barangays in the Municipality of Tuba in Benguet, Tabaan Sur is projected to have the highest
percentage of area that will experience flood levels, at 7.84%. Meanwhile, Poblacion posted the second
Among the barangays in the Municipality of Tublay in Benguet, Tuel is projected to have the highest
percentage of area that will experience flood levels, at 19.50%. Meanwhile, Daclan posted the second
Among the barangays in the Municipality of Agoo in La Union, Macalva Norte is projected to have the
highest percentage of area that will experience flood levels, at 5.02%. Meanwhile, Santa Ana posted the
Among the barangays in the Municipality of Aringay in La Union, Gallano is projected to have the highest
percentage of area that will experience flood levels, at 27.30%. Meanwhile, Manga posted the second
Among the barangays in the Municipality of Bagulin in La Union, Tio-Angan is projected to have the highest
percentage of area that will experience flood levels, at 19.03%. Meanwhile, Alibangsay posted the second
Among the barangays in the Municipality of Bauang in La Union, Upper San Agustin is projected to have
the highest percentage of area that will experience flood levels, at 11.41%. Meanwhile, Acao posted the
Among the barangays in the Municipality of Burgos in La Union, Delles is projected to have the highest
percentage of area that will experience flood levels, at 23.94%. Meanwhile, Upper Tumapoc posted the
Among the barangays in the Municipality of Caba in La Union, Liquicia is projected to have the highest
percentage of area that will experience flood levels, at 28.94%. Meanwhile, San Cornelio posted the
Among the barangays in the Municipality of Naguilian in La Union, San Isidro is projected to have the
highest percentage of area that will experience flood levels, at 9.38%. Meanwhile, Casilagan posted the
Among the barangays in the Municipality of Pugo in La Union, Palina is projected to have the highest
225
percentage of area that will experience flood levels, at 36.33%. Meanwhile, Ambangonan posted the
Among the barangays in the Municipality of Rosario in La Union, Marcos is projected to have the highest
percentage of area that will experience flood levels, at 2.21%. Meanwhile, Parasapas posted the second
Among the barangays in the Municipality of San Fernando City in La Union, Bacsil is projected to have
the highest percentage of area that will experience flood levels, at 5.28%. Meanwhile, Sacyud posted the
Among the barangays in the Municipality of San Gabriel in La Union, Apayao is projected to have the
highest percentage of area that will experience flood levels, at 8.45%. Meanwhile, Amontoc posted the
Among the barangays in the Municipality of Santo Tomas in La Union, Pongpong is projected to have the
highest percentage of area that will experience flood levels, at 7.99%. Meanwhile, Ambitacay posted the
Among the barangays in the Municipality of Tubao in La Union, Santa Teresa is projected to have the
highest percentage of area that will experience flood levels, at 15.56%. Meanwhile, Rizal posted the second
The generated flood hazard maps for the Bauang floodplain were also used to assess the vulnerability
of the educational and medical institutions in the floodplain. Using the flood depth units of PAGASA for
the hazard maps – “Low”, “Medium”, and “High” – the affected institutions were given an individual
assessment for each flood hazard scenario (i.e., 5-year, 25-year, and 100-year).
Table 150. Area covered by each warning level, with respect to the rainfall scenario
Area Covered in sq. km.
Warning Level
5 year 25 year 100 year
Low 37.57 42.03 45.01
Medium 40.35 42.10 43.74
High 108.03 135.46 151.95
TOTAL 185.95 219.59 240.7
Of the fifty-five (55) identified educational institutions in the Bauang floodplain, eighteen (18) were
assessed to be exposed to High-level flooding for all three (3) rainfall return scenarios. Three (3) other
schools were found to be susceptible to flooding, projected to experience Medium-level flooding in the
5-year return period, and High-level flooding in the 25-year and 100-year rainfall scenarios. See Annex 12
for a detailed enumeration of schools in the Bauang floodplain.
Six (6) medical institutions were identified in the Bauang floodplain. Naguilian Health Center in Barangay
Cabaritan Sur and Naguilian District Hospital in Barangay Suguidan Norte were found to be highly prone
to flooding, computed to have High-level flooding in all three (3) rainfall scenarios. See Annex 13 for a
detailed enumeration of hospitals and clinics in the Bauang floodplain.
226
5.11 Flood Validation
In order to check and validate the extent of flooding in the different river systems, there is a need to
perform validation survey work. For this purpose, field personnel gathered secondary data regarding flood
occurrences in the respective areas within the major river systems in the Philippines.

From the flood depth maps produced by the Phil-LiDAR 1 Program, multiple points representing the
different flood depths for different scenarios were identified for validation.

The validation personnel then went to the specified points identified in a river basin to gather data regarding
the actual flood levels in each location. Data gathering was conducted through assistance from a local
DRRM office to obtain maps or situation reports about the past flooding events, or through interviews with
some residents with knowledge or experience of flooding in a particular area.

After which, the actual data from the field were compared with the simulated data to assess the accuracy
of the flood depth maps produced, and to improve on the results of the flood map. The points in the flood
map versus the corresponding validation depths are illustrated in Figure 192.
The flood validation survey was conducted in January 2017. The flood validation consists of one hundred
and ninety-seven (197) points randomly selected all over the Bauang floodplain. It attained an RMSE value
of 2.97. Table 151 presents a contingency matrix of the comparison. The validation points are found in
Annex 11.
The validation data were obtained on January 2017
227
Figure 191. Flood validation points of the Bauang River Basin
Figure 192. Flood map depth vs. actual flood depth for Bauang
228
Table 151. Actual flood depth vs. simulated flood depth in the Bauang River Basin
ARINGAY-BAU- Modeled Flood Depth (m)

ANG BASIN 0-0.20 0.21-0.50 0.51-1.00 1.01-2.00 2.01-5.00 > 5.00 Total
0-0.20 0 0 0 0 2 0 2
Actual Flood Depth (m)
0.21-0.50 39 21 20 25 50 20 175
0.51-1.00 7 3 1 6 2 0 19
1.01-2.00 1 0 0 0 0 0 1
2.01-5.00 0 0 0 0 0 0 0
> 5.00 0 0 0 0 0 0 0
Total 47 24 21 31 54 20 197
The overall accuracy generated by the flood model is estimated at 11.17%, with twenty-two (22) points
correctly matching the actual flood depths. In addition, there were sixty-eight (68) points estimated one
(1) level above and below the correct flood depths. Meanwhile, there were thirty-four (34) points and
seventy-three (73) points estimated two (2) levels above and below, and three (3) or more levels above
and below the correct flood depths, respectively. A total of one hundred and twenty-five (125) points
were overestimated, while a total of fifty (50) points were underestimated in the modeled flood depths of
Bauang.
Table 152. Summary of Accuracy Assessment in the Bauang River Basin

No. of Points %
Correct 22 11.17
Overestimated 125 63.45
Underestimated 50 25.38
Total 197 100
229
REFERENCES
Ang M.O., Paringit E.C., et al. 2014. DREAM Data Processing Component Manual. Quezon City, Philippines:
UP Training Center for Applied Geodesy and Photogrammetry.
Balicanta L.P., Paringit E.C., et al. 2014. DREAM Data Validation Component Manual. Quezon City,
Philippines: UP Training Center for Applied Geodesy and Photogrammetry.
Brunner, G. H. 2010a. HEC-RAS River Analysis System Hydraulic Reference Manual. Davis, CA: U.S.
ArmyCorps of Engineers, Institute for Water Resources, Hydrologic Engineering Center.
Lagmay A.F., Paringit E.C., et al. 2014. DREAM Flood Modeling Component Manual. Quezon City, Philippines:
UP Training Center for Applied Geodesy and Photogrammetry.
Paringit E.C, Balicanta L.P., Ang, M.O., Sarmiento, C. 2017. Flood Mapping of Rivers in the Philippines
Using Airborne Lidar: Methods. Quezon City, Philippines: UP Training Center for Applied Geodesy and
Photogrammetry.
Sarmiento C., Paringit E.C., et al. 2014. DREAM Data Acquisition Component Manual. Quezon City,
Philippines: UP Training Center for Applied Geodesy and Photogrammetry.
UP TCAGP 2016, Acceptance and Evaluation of Synthetic Aperture Radar Digital Surface Model (SAR DSM)
and Ground Control Points (GCP). Quezon City, Philippines: UP Training Center for Applied Geodesy and
Photogrammetry.
230
ANNEXES
Annex 1. Technical Specifications of the Pegasus LiDAR Sensor used in the Bauang
Floodplain Survey
Table A-1.1. Technical specifications of the Pegasus sensor

Parameter Specification
Operational envelope (1,2,3,4) 150-4000 m AGL, nominal
Laser wavelength 1064 nm
Horizontal accuracy (2) 1/5,500 x altitude, (m AGL)
Elevation accuracy (2) <5-35 cm, 1 σ
Effective laser repetition rate Programmable, 33-167 kHz
POS AV™ AP50 (OEM);
Position and orientation system 220-channel dual frequency GPS/GNSS/Galile-
o/L-Band receiver
Scan width (WOV) Programmable, 0-50˚
Scan frequency (5) Programmable, 0-70 Hz (effective)
Sensor scan product 1000 maximum
Dual divergence: 0.25 mrad (1/e) and 0.8 mrad
Beam divergence
(1/e), nominal
Roll compensation Programmable, ±5˚ (FOV dependent)
Up to 4 range measurements, including 1st, 2nd,
Range capture
3rd, and last returns
Up to 4 intensity returns for each pulse, including
Intensity capture
last (12 bit)
Video Camera Internal video camera (NTSC or PAL)
Image capture Compatible with full Optech camera line (optional)
12-bit Optech IWD-2 Intelligent Waveform Digitizer
Full waveform capture
(optional)
Data storage Removable solid state disk SSD (SATA II)
Power requirements 28 V; 900 W;35 A(peak)
Sensor: 260 mm (w) x 190 mm (l) x 570 mm (h); 23
kg
Dimensions and weight
Control rack: 650 mm (w) x 590 mm (l) x 530 mm
(h); 53 kg
Operating temperature -10˚C to +35˚C (with insulating jacket)
Relative humidity 0-95% no-condensing
231
Annex 2. NAMRIA Certification of Reference Points used in the LiDAR Survey
1. LUN-62
Figure A-2.1. LUN-62
232
2. LUN-176
233
3. LUN-3062
234
4. LUN-3047
235
5. LUN-3129
236
6. LU-94
Figure A-2.6. LU-94
237
Annex 3. Baseline Processing Reports of Control Points used in the LiDAR Survey
Table A-3.1. LU-94
238
Annex 4. The LiDAR Survey Team Composition
Table A-4.1. LiDAR Survey Team Composition

Data Acquisition
Component Designation Name Agency / Affiliation
Sub -Team
PHIL-LIDAR 1 Program Leader ENRICO C. PARINGIT, D.ENG UP-TCAGP
Data Acquisition Data Component Project

ENGR. CZAR JAKIRI SARMIENTO UP-TCAGP
Component Leader Leader – I
Chief Science Research
ENGR. CHRISTOPHER CRUZ UP-TCAGP
Survey Supervisor Specialist (CSRS)
FIELD TEAM
Senior Science Research
LOVELY GRACIA ACUNA UP-TCAGP
Specialist (SSRS)
LiDAR Operation Research Associate (RA) RENAN PUNTO UP-TCAGP
RA FAITH JOY SABLE UP-TCAGP
Ground Survey,
Data Download and RA KENNETH QUISADO UP-TCAGP
Transfer
PHILIPPINE AIR
Airborne Security SSG. OLIVER SACLOT
FORCE (PAF)
CAPT. MARK LAWRENCE TAN- ASIAN AEROSPACE
LiDAR Operation GONAN CORPORATION (AAC)
Pilot
CAPT. NEIL ACHILLES AGAWIN AAC
239
Annex 5. Data Transfer Sheets for the Bauang Floodplain Flights
Figure A-5.1. Data Transfer Sheet for the Bauang Floodplain
240
Annex 6. Flight Logs for the Flight Missions
1. Flight Log for 1153P Mission
Figure A-6.1. Flight Log for Mission 1153P
241
242
243
244
245
246
247
248
249
Annex 7. Flight Status Reports
Table A-7.1. Flight Status Report
NORTHERN MINDANAO
(February 25 – March 8, 2014)
FLIGHT DATE
AREA MISSION OPERATOR REMARKS
NO FLOWN
February
25, 2014 Survey voids of Block10A and
1153P BLOCK 10A 1BLK10AS056B F. SABLE 1 line Blk10B; renamed from
1151P
February Survey Block 10C with data

1155P BLOCK 10C 1BLK10C057A R. PUNTO 26, 2014 voids due to eye safety, laser
sets off; renamed from 1153P
February Survey Block 10G and 10D with
BLOCK 27, 2014 data voids due to eye safety,
1159P 1BLK10GD058A R. PUNTO
10GD laser sets off; renamed from
1157P
February
1161P BLK10D 1BLK10D058B F SABLE Survey Block 10D
27, 2014
February
1163P BLK10F 1BLK10F059A R PUNTO Finished Survey
28, 2014
BLOCK 10D March 2, Supplementary flight to cover
1171P & BLOCK 1BLK10CDS061A R. PUNTO 2014 voids for Block10D & Block10C;
10C renamed from 1169P
March 2, Supplementary flight to com-
1173P BLOCK 10D 1BLK10DGS061B F. SABLE 2014 plete Block 10D; renamed from
1171P
March 3, Supplementary flight to cover
1177P BLOCK 10C 1BLK10CS062B F. SABLE 2014 voids in Block 10C; renamed
from 1175P
BLOCK March 8,
1197P 1BLK10GCS067B R.PUNTO Mission Complete
10G,10C 2014
250
LAS BOUNDARIES PER FLIGHT
Flight No. : 1153P

Area: BLK 10A
Mission Name: 1BLK10BAS056B
Parameters: Altitude: 1200m; Scan Frequency: 30Hz;
Scan Angle: 25deg; Overlap: 30%
LAS
Figure A-7.1. Swath for Flight No. 1153P
Flight No. : 1155P

251
Area: BLK 10C

Mission Name: 1BLK10C057A
Parameters: Altitude: 1200 m; Scan Frequency: 30Hz;
Scan Angle: 25deg; Overlap: 30%
LAS
252
Flight No. : 1159P

Area: BLK 10G
Mission Name: 1BLK10GD058A
Scan Angle: 25 deg; Overlap: 30%
LAS
253
Flight No. : 1161P

Area: BLK10D
Mission name: 1BLK10D058B
Parameters: Altitude: 1200 m; Scan Frequency: 50 Hz;
Scan Angle: 15 deg; Overlap: %
Area covered: 188.09 km2
254
Flight No. : 1163P

Area: BLK10F
Mission name: 1BLK10F059A
Scan Angle: 15 deg; Overlap: 50 %
255
Flight No. : 1171P

Area: BLK 10C
Mission Name: 1BLK10CDS061A
LAS
256
Flight No. : 1173P

Area: BLK 10D
Mission Name: 1BLK10DS061B
LAS
257
Flight No. : 1177P

Area: BLK 10C
Mission Name: 1BLK10CS062B
LAS
258
Flight No. : 1197P

Area: BLK10GC
Mission name: 1BLK10GCS067B
Scan Angle: 15 deg; Overlap: 50 %
LAS
259
Annex 8. Mission Summary Reports
Table A-8.1. Mission Summary Report for Mission Blk10C

Flight Area La Union
Mission Name Blk10C
Inclusive Flights 1155P, 1171P, 1177P, 1197P
Range data size 62.86 GB
Base data size 22.2 MB
POS 695 MB
Image 98.6 GB
Transfer date March 08, 2014

Solution Status
Number of Satellites (>6) Yes
PDOP (<3) Yes
Baseline Length (<30km) No
Processing Mode (<=1) No

Smoothed Performance Metrics (in cm)
RMSE for North Position (<4.0 cm) 3.62
RMSE for East Position (<4.0 cm) 4.3
RMSE for Down Position (<8.0 cm) 6.55

Boresight correction stdev (<0.001deg) 0.000398
IMU attitude correction stdev (<0.001deg) 0.017218
GPS position stdev (<0.01m) 0.0267

Minimum % overlap (>25) 35.05%
Ave point cloud density per sq.m. (>2.0) 3.75
Elevation difference between strips (<0.20 m) Yes
Number of 1km x 1km blocks 505

Maximum Height 1038.87 m
Minimum Height 45.5 m
Classification (# of points)
Ground 270,659,389
Low vegetation 232,081,137
Medium vegetation 638,506,120
High vegetation 737,644,888
Building 23,046,935
Orthophoto Yes
Engr. Kenneth Solidum, Engr. Melanie Hingpit,
Processed by
Ailyn Biñas
260
Figure A-8.1. Solution Status
Figure A-8.2. Smoothed Performance Metric Parameters
261
Figure A-8.3. Best Estimated Trajectory
Figure A-8.4. Coverage of LIDAR data
262
Figure A-8.5. Image of Data Overlap
Figure A-8.6. Density map of merged LIDAR data
263
Figure A-8.7. Elevation difference between flight lines
264
Table A-8.2. Mission Summary Report for Mission Blk10D

Mission Name Blk10D
Inclusive Flights 1159P, 1161P, 1171P, 1173P
Range data size 75 GB
POS 695 MB
Image 68.3 GB

Solution Status
Number of Satellites (>6) Yes
PDOP (<3) Yes
Baseline Length (<30km) Yes




Maximum Height 401.8
Minimum Height 42.44
Ground 301,413,519
Building 30,496,828
Orthophoto
Engr. Kenneth Solidum, Engr. Angelo Carlo
Processed by Bongat, Engr. Benjamin Jonah Magallon, Engr.
Harmond Santos, Engr. Jeffrey Delica
265
266
267
268
Figure A-8.14 Elevation difference between flight lines
269
Table A-8.3. Mission Summary Report for Mission Blk10F

Mission Name Blk10F
Inclusive Flights 1163P
POS 216 MB
Image 26.6 GB
Transfer date February 28, 2014

Solution Status
Number of Satellites (>6) No
PDOP (<3) Yes
Baseline Length (<30km) No




Maximum Height 917.14
Minimum Height 54.72
Ground 313,178,442
Building 16,826,260
Orthophoto
Engr. Irish Cortez, Engr. Edgardo Gubatanga, Jr.,
Processed by
Ailyn Biñas
270
271
272
273
274
Table A-8.4. Mission Summary Report for Mission Blk10G

Mission Name Blk10G
Inclusive Flights 1159P, 1197P
POS 333 MB
Image 19.2 GB

Solution Status
Number of Satellites (>6) No
PDOP (<3) Yes
Baseline Length (<30km) Yes


Boresight correction stdev (<0.001deg) n/a
IMU attitude correction stdev (<0.001deg) n/a
GPS position stdev (<0.01m) n/a


Maximum Height 907.57 m
Minimum Height 63.17 m
Ground 87,019,657
Building 6,668,084
Orthophoto Yes
Engr. Angelo Carlo Bongat, Engr. Kenneth Solidum,
Processed by
Engr. Harmond Santos, Simonette Lat
275
276
277
278
279
Annex 9. Bauang Model Basin Parameters
Table A-9.1. Bauang Model Basin Parameters
Clark Unit
SCS Curve Number Loss Recession Baseflow
Hydrograph Transform
Basin
Number Initial Ab-
Curve Impervious Time of Con- Storage Coef- Initial Discharge Recession Threshold Ratio to
straction Initial Type
Number (%) centration (HR) ficient (HR) (cms) Constant Type Peak
(mm)
W1000 426.92 52.5 0 4.1347 2.9943 Discharge 0.7729 0.77141 Ratio to Peak 0.22199
W1070 348.315 35 0 2.0717 3.3811 Discharge 0.58396 1 Ratio to Peak 0.23
280

W1100 398.75 40.467 0 0.22153 2.8015 Discharge 0.10226 1 Ratio to Peak 0.15253
Clark Unit
Basin
Number Initial Ab-
(mm)
W610 174 40.336 0 9.1514 5.4524 Discharge 2.7478 1 Ratio to Peak 0.21316
281
Clark Unit
Basin
Number Initial Ab-
(mm)
Clark Unit
Basin
282
Number Initial Ab-

(mm)
Annex 10. Bauang Model Reach Parameters
Table A-10.1. Bauang Model Reach Parameters
Reach Muskingum Cunge Channel Routing

Num- Length Man- Side
ber Time Step Method Slope Shape Width
(m) ning’s n Slope
R1190 Automatic Fixed Interval 3354.8 0.004709 0.002321 Trapezoid 406 1
R170 Automatic Fixed Interval 7003.8 0.004856 0.004 Trapezoid 72.6 1

R180 Automatic Fixed Interval 2825 0.004294 0.004 Trapezoid 209 1


Reach Muskingum Cunge Channel Routing

Num- Length Manning’s Side
ber Time Step Method Slope Shape Width
(m) n Slope



283

284
Annex 11. Bauang Field Validation Points
Table A-11.1. Bauang Field Validation Points
Validation Coordinates Vali- Rain

Point Model
dation Return/
Num- Var Error (m) Event/Date
Points Scenar-
ber Lat Long (m)
(m) io
1 16.51108 120.32204 0.320 0.305 0.000 Mario/ September 18-22, 2014 5-Year
19 16.381372 120.35388 0.330 0.457 0.016 Mina/ November 24-27, 2007 5-Year
31 16.536744 120.32905 0.100 0.305 0.042 Pepeng/ October 2-5, 2009 5-Year
285
286
Typhoon Juan (Oct. 17-20,
106 16.514998 120.32929 0.050 1.219 1.367 5-Year
2010)
128 16.526032 120.3199 2.070 0.254 3.298 Lawin/ October 18-22, 2016 5-Year
287
Typhoon Juan (Oct. 17-20,
133 16.52746 120.32172 2.200 0.254 3.787 5-Year
2010)
288
RMSE 2.971563
289
Annex 12. Educational Institutions Affected by Flooding in Bauang Floodplain
Table A-12.1. Educational Institutions Affected by Flooding in the Bauang Floodplain
La Union
Bagulin
Rainfall Scenario
Building Name Barangay
5-year 25-year 100-year
BAGULIN CENTRAL ELEMENTARYSCHOOL Dagup Medium High High
SUYO NATIONAL HIGHSCHOOL Dagup High High High
Bauang
Rainfall Scenario
BALLAY COMMUNITYSCHOOL / BALLAY
Acao High High High
NATIONAL HIGH SCHOOL
BACCUIT ELEMEMTARYSCHOOL Baccuit Sur Medium High High
BACCUIT NATIONAL HIGHSCHOOL Baccuit Sur Medium Medium High
CABARITAN ELEMENTARYSCHOOL Ballay High High High
BOY-UTAN ELEMENTARYSCHOOL Boy-Utan High High High
CABALAYANGAN ELEMENTARYSCHOOL Bucayab
SUGUIDAN INTEGRATEDSCHOOL Cabalayangan Medium High
CALUMBAYA ELEM.SCHOOL Calumbaya
DON EULOGIO MEMORIAL NATIONAL
Calumbaya
HIGHSCHOOL
BAUANG LIBRARY Central East
BAUANG NORTH CENTRALSCHOOL Central East
SACRED HEARTSCHOOL Central East
STS. PETER AND PAUL LEARNING CENTER Central East
QUINAVITE ELEMENTARYSCHOOL Disso-Or Low
TUBAO CREDIT COOPERATIVE Disso-Or Low
YESHUAH THE MESIAH CHRISTIAN ACAD-
Disso-Or Low
EMY
YESHUAH THE MESSIAH CHRISTIAN
Disso-Or Low Low
ACADEMY
DAY CARE CENTER Parian Este Low Medium Medium
PARIAN ELEMENTARYSCHOOL Parian Este Low Medium Medium
PARIAN ESTE ELEMENTARYSCHOOL Parian Este Low Medium Medium
PUDOK PRIMARYSCHOOL Parian Este Low Medium Medium
Payocpoc Norte
GUERRERO ELEMENTARYSCHOOL
Este
Payocpoc Norte
PAYOCPOC ELEMENTARYSCHOOL Low Low
Este
LA UNION COLLEGES OFSCIENCE AND
Pugo High High High
TECHNOLOGY
KIDDIE CAMP PLAYSCHOOL Quinavite Low
290
ST. ANTHONY MONTESSORI EDUCATION-

Quinavite
AL NETWORK
SANTIAGO ELEMENTARYSCHOOL Santiago Low Low
BEEHIVE LEARNING CENTER Taberna High High High
PUGO ELEMENTARYSCHOOL Taberna High High High
Burgos
Rainfall Scenario
AGPAY ELEMENTARYSCHOOL Agpay
Naguilian
Rainfall Scenario
EASTERN NAGUILIAN NATIONAL HIGH- Ambaracao
SCHOOL Norte
Ambaracao
GUESSET ELEMENTARYSCHOOL
Norte
Ambaracao
GUESSET ELEMENTARYSCHOOL
Norte
DR. HERMOGGENES BELEN ELEMENTA-
Cabaritan Sur
RYSCHOOL
GIRLSCOUT OF THE PHILIPPINES REGION-
Cabaritan Sur High High High
AL OFFICE
NAGUILIAN DISTRICT OFFICE Cabaritan Sur High High High
NAGUILIAN ELEMENTARYSCHOOL Cabaritan Sur High High High
BARIQUIR ELEMENTARYSCHOOL Casilagan Low Medium High
UNIONSCHOOLS OF LA UNION Dal-Lipaoen High High High
BARAOAS SUR ELEMENTARYSCHOOL Daramuangan High High High
CASILAGAN ELEMENTARYSCHOOL Gusing Sur High High High
LOWER BIMMOTOBOT ELEMENTARY- Mamat-Ing
High High High
SCHOOL Norte
GUSING ELEMENTARYSCHOOL Mamat-Ing Sur High High High
GUSING NATIONAL HIGHSCHOOL Mamat-Ing Sur High High High
TUDDINGAN ELEMENTARYSCHOOL Mamat-Ing Sur High High
NAGUILIAN NATIONAL HIGHSCHOOL Natividad High High High
BURGOS CENTRALSCHOOL Sili
DON TOMAS R. MENDOZA ELEMENTARY-
Sili
SCHOOL
HOLY INFANT NINO MONTESSORI Suguidan Norte Medium High High
BATO ELEMENTARYSCHOOL Suguidan Sur
NAGSIDORISAN ELEMENTARYSCHOOL Suguidan Sur High High
San Fernando City

Rainfall Scenario
MASICONG ELEMENTARYSCHOOL Birunget Low Low
MAMAT-ING SUR ELEMENTARYSCHOOL Sacyud High High High
291
Annex 13. Medical Institutions Affected by Flooding in Bauang Floodplain
Table A-13.1. Medical Institutions Affected by Flooding in the Bauang Floodplain
La Union
Bauang
Rainfall Scenario
BAUANG RURAL HEALTH UNIT Central East
MUNICIPAL HEALTH CENTER Central East
GENERIKA DRUG STORE Quinavite Low Low Low
BARANGAY TABEMA HEALTH CENTER (TL) Taberna Medium High High
Naguilian
Rainfall Scenario
NAGUILIAN HEALTH CENTER Cabaritan Sur High High High
NAGUILIAN DISTRICT HOSPITAL Suguidan Norte High High High
292

LiDAR Surveys and Flood Mapping of Bauang River PDF

Uploaded by

Copyright:

Available Formats

LiDAR Surveys and Flood Mapping of Bauang River PDF

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

LiDAR Surveys and Flood Mapping of Bauang River PDF

Uploaded by

Copyright:

Available Formats

Hazard Mapping of the Philippines Using LIDAR (Phil-LIDAR 1)

LiDAR Surveys and Flood Mapping of Bauang River

For questions/queries regarding this report, contact:

Dr. Chelo Pascua

Enrico C. Paringit, Dr. Eng.

National Library of the Philippines

Figure 45. Bridge cross-section diagram of the Mamat-ing Bridge........................................................... 62

LIST OF ACRONYMS AND ABBREVIATIONS

CHAPTER 1: OVERVIEW OF THE PROGRAM AND

1.1 Background of the Phil-LiDAR 1 Program

1.2 Overview of the Bauang River Basin

Figure 1. Location map of the Bauang River Basin (in brown)

CHAPTER 2: LIDAR DATA ACQUISITION OF THE

2.1 Flight Plans

Table 1. Flight planning parameters for the Pegasus LiDAR system

2.2 Ground Base Stations

Figure 6. LUN-3047, as recovered in Barangay Nazareno, Agoo, La Union

Figure 8. NAMRIA reference point LU-94, as recovered by the field team

Date Surveyed Flight Number Mission Name Ground Control Points

February 25, 2014 1153P 1BLK10AS056B LUN-176 and LUN-3129

February 26, 2014 1155P 1BLK10C057A LUN-176 and LU-94

February 27, 2014 1159P 1BLK10GD058A LUN-62 and LUN-302

February 27, 2014 1161P 1BLK10D058B LUN-62 & LUN-3062

February 28, 2014 1163P 1BLK10F059A LUN-72 & LUN-3047

March 2, 2014 1171P 1BLK10CDS061A LUN-62 and LUN-3062

March 2, 2014 1173P 1BLK10DS061B LUN-62 and LUN-3062

March 3, 2014 1177P 1BLK10CS062B LUN-176 and LUN-3129

March 8, 2014 1197P 1BLK10GCS067B LUN-62 and LUN-3062

2.3 Flight Missions

TOTAL 2977.809 1906.437 308.936 1597.502 3263 26 7

Flying Scan Fre- Average Average

2.4 Survey Coverage

Figure 9. Actual LiDAR survey coverage for the Bauang floodplain

CHAPTER 3: LIDAR DATA PROCESSING OF THE BAUANG

3.1 Overview of the LIDAR Data Pre-Processing

These processes are summarized in the diagram in Figure 10.

Figure 10. Schematic diagram for the Data Pre-Processing Component

3.2 Transmittal of Acquired LiDAR Data

3.3 Trajectory Computation

Figure 11. Smoothed Performance Metric Parameters of Bauang Flight 1153P

Figure 12. Solution Status Parameters of Bauang Flight 1153P

3.4 LiDAR Point Cloud Computation

Table 12. Self-calibration results values for the Bauang flights

3.5 LiDAR Data Quality Checking

Table 13. List of LiDAR blocks for the Bauang floodplain

Figure 15. Image of data overlap for Bauang floodplain

3.6 LiDAR Point Cloud Classification and Rasterization

Table 14. Bauang classification results in TerraScan

3.7 LiDAR Image Processing and Orthophotograph Rectification

Figure 22. The Bauang floodplain, with available orthophotographs

Figure 23. Sample orthophotograph tiles for the Bauang floodplain

3.8 DEM Editing and Hydro-Correction

Table 15. LiDAR blocks with their corresponding areas

3.9 Mosaicking of Blocks

3.10 Calibration and Validation of Mosaicked LiDAR DEM

Table 17. Calibration statistical measures

Calibration Statistical Measures Value (meters)

Table 18. Validation statistical measures

Validation Statistical Measures Value (meters)