Applications of Hybrid Algorithm (Successive Over Relaxation and Inverse Distance Weighting) for Interpolating Rainfall Data Obtained from a Dense Network of Meteorological Stations in Metro Manila, Philippines


A hybrid piecewise rainfall value interpolation algorithm was formulated using the commonly known Inverse Distance Weighting (IDW) and Gauss-Seidel variant Successive Over Relaxation (SOR) to interpolate rainfall values over Metro Manila, Philippines. Due to the fact that the SOR requires boundary values for its algorithm to work, the IDW method has been used to estimate rainfall values at the boundary. Iterations using SOR were then done on the defined boundaries to obtain the desired results corresponding to the lowest RMSE value. The hybrid method was applied to rainfall datasets obtained from a dense network of 30 stations in Metro Manila which has been collecting meteorological data every 5 minutes since 2012. Implementing the Davis Vantage Pro 2 Plus weather monitoring system, each station sends data to a central server which could be accessed through the website The stations are spread over approximately 625 sq km of area such that each station is approximately within 25 sq km from each other. The locations of the stations determined by the Metro Manila Development Authority (MMDA) are in critical sections of Metro Manila such as watersheds and flood-prone areas. Three cases have been investigated in this study, one for each type of rainfall present in Metro Manila: monsoon-induced (8/20/13), typhoon (6/29/13), and thunderstorm (7/3/15 & 7/4/15). The area where the rainfall stations are located is divided such that large measured rainfall values are used as part of the boundaries for the SOR. Measured station values found inside the area where SOR is implemented are compared with results from interpolated values. Root mean square error (RMSE) and correlation trends between measured and interpolated results are quantified. Results from typhoon, thunderstorm and monsoon cases show RMSE values ranged from 0.25 to 2.46 mm for typhoons, 1.55 to 10.69 mm for monsoon-induced rain and 0.01 to 6.27 mm for thunderstorms. R2 values, on the other hand, are 0.91, 0.89 and 0.76 for typhoons, monsoon-induced rain and thunderstorms, respectively. This study has shown that the method of approximating rainfall works and can be used in improved prediction, analysis and real time flood map generation.