JUser: :_load: Unable to load user with ID: 59
Thursday, 30 April 2020 11:10

Watershed Modelling of the Mindanao River Basin in the Philippines Using the SWAT for Water Resource Management

Written by
Rate this item
(1 Vote)


This study aims to simulate the watershed of the Mindanao River Basin (MRB) to enhance water resource management for potential hydropower applications to meet the power demand in Mindanao with an average growth of 3.8% annually. The soil and water assessment tool (SWAT) model was used with inputs for geospatial datasets and weather records at four meteorological stations from DOST-PAGASA. To overcome the lack of precipitation data in the MRB, the precipitation records were investigated by comparing the records with the global gridded precipitation datasets from the NCDC-CPC and the GPCC. Then, the SWAT simulated discharges with the three precipitation data were calibrated with river discharge records at three stations in the Nituan, Libungan and Pulangi rivers. Due to limited records for the river discharges, the model results were, then, validated using the proxy basin principle along the same rivers in the Nituan, Libungan, and Pulangi areas. The R2 values from the validation are 0.61, 0.50 and 0.33, respectively, with the DOSTPAGASA precipitation; 0.64, 0.46 and 0.40, respectively, with the NCDC-CPC precipitation; and 0.57, 0.48 and 0.21, respectively, with the GPCC precipitation. The relatively low model performances in Libungan and Pulangi rivers are mainly due to the lack of datasets on the dam and water withdrawal in the MRB. Therefore, this study also addresses the issue of data quality for precipitation and data scarcity for river discharge, dam, and water withdrawal for water resource management in the MRB and show how to overcome the data quality and scarcity.

Keywords: SWAT; Mindanao River Basin; Discharge; Watershed Modelling; Precipitation; Proxy Basin.


1.   Introduction

Among the developing countries, the Philippines faces a considerable challenge regarding development due to the continuous increase in electricity demands, with an annual average rate of increase of 4.3% [1]. The power demand of the Mindanao island group in the Philippines has increased by 3.8% annually over the past decades [2]. In April 2017, the maximum power peak demand in Mindanao reached approximately 1,696 MW [3]. However, the Mindanao water resources contributed 38%, or 1,947 GWh, of the gross power generation from hydropower in June 2017 [3]. Regardless of the current contribution of water resources to renewable energy, the power demand continues to outpace the supply. Thus, to address this emerging problem, assessment for a potential source of sustainable renewable energy is needed. The purpose of this study is to enhance water resource management for hydropower application in Mindanao to improve the electrification situation and support the implementation of the Renewable Energy Act of the Philippines.


* Corresponding author: This email address is being protected from spambots. You need JavaScript enabled to view it.


The renewable energy resources in the Philippines are geothermal, wind, solar, biomass, ocean and hydropower resources [1]. However, among these options, hydropower is more sustainable in this country due to the abundance of water resources in its 18 major basins [4]. Hydropower has the greatest potential, with an estimated contribution of 13.31% of the energy needs of the country [5]. Therefore, to maximize the utilization of water resources for hydropower, assessment of available water resources has to be carried out in the major basins of the country.

On the other hand, water resources are also utilized for irrigation for agricultural productivity; these irrigation systems cover 52.0% and 38.6% of the Philippines and Mindanao, respectively [6]. Hence, the Mindanao irrigation service covers a total of 20,212.71 ha, contributing to the primary income-generating agriculture industry [4]. Moreover, water resources play an important role in the community; for instance, only 82.6% and 86.8% of households have access to safe water supplies in 2011 in the Philippines and Mindanao, respectively. This low rate of access to potable water results in outbreaks of diseases carried by water [6].

The concerns of water resource management in Mindanao are severely affected by geological and hydrogeological hazards due to the physical environment. Mindanao is vulnerable to disasters induced by natural hazards such as storm surges, typhoons, earthquakes, tsunamis, droughts and floods [7, 8].

In 2011, Tropical Storm Washi (known as Tropical Storm Sendong in the Philippines) made landfall in the northern part of Mindanao and caused a heavy rain that led to overflow of the Cagayan Basin, resulting in calamitous flooding in Cagayan de Oro and Iligan City and in 14 provinces, with an estimated damage of 4.17 million USD to agriculture and 0.78 million USD to fisheries [9, 10]. In addition, Typhoon Bopha, caused damage in eastern Mindanao with an overall estimated cost to agriculture of 645 million USD [11]. Furthermore, flooding events occur due to extreme rainfall, tropical cyclones from Monsoon winds and the dynamic climate of tropical cyclones with low pressures [7].

This weather dynamic is very important to consider for the development of water resource management because of its direct influence on the watershed. For instance, high precipitation intensity may cause a flood because of the direct impact of precipitation on river runoff and slow ground absorption. Thus, more precipitation results in a higher possibility of flooding [12]. Therefore, the assessment of the sustainability of the water supply for hydropower application mainly depends on the characteristics of precipitation.


Figure 1. Study area of the Mindanao island group, Philippines: (a) the 17 regional administrative boundaries of the Philippines; (b) the population at the provincial level; (c) Mindanao, showing the population at the municipality/city level; and d) the major basins, weather stations and gridded precipitation points in Mindanao used in the SWAT simulations. The dashed lines in Figure 1(d) are the four points used for the estimation of rainfall patterns from the precipitation datasets.

Thus, this study carefully investigates the precipitation inputs in watershed modelling by comparing the observational data from the DOST-PAGASA with the global gridded precipitation datasets from the NCDC-CPC and GPCC. In addition, water management can include agricultural water footprint analysis by considering increasingly complex indicators, such as multiple climate variables, soil characteristics, and crop properties, to simulate the water cycle [13]. Therefore, weather-related events must be understood in terms of their possible effects on the watershed since hydropower generation is driven by river discharges. Moreover, SWAT considers other climate variables such as temperature, humidity and solar radiation, while soil characteristics and crop properties are some of the inputs to the hydrologic response unit (HRU) of SWAT [14].

Therefore, the objectives of this study are to evaluate the observed rainfall data with gridded global precipitation datasets to overcome the present data scarcity and to simulate the river discharges of the MRB for water resource management and future hydropower development. The literature review on watershed modelling with SWAT is presented in Section 2. The study area, material and method are presented in Section 3, 4 and 5, respectively. The results and discussion are described in Section 6 and 7, respectively, followed by conclusion in Section 8. 

2.   Watershed Modeling with SWAT in the Available Literature

The SWAT model was developed in the early 1990s by Jeff Arnold [14], and it has been recognized worldwide as an effective tool in water resource management for assessing the impact of the climate on water supplies and nonpoint sources of pollution in watersheds [15, 16]. Moreover, SWAT is a scientific tool used to evaluate streamflow, agricultural chemicals and sediment yield in a large basin [17]. For instance, SWAT was applied to a semi-arid climate in India for rainfall-runoff modelling of river basins [18]. SWAT was also applied to short-term climate data for the assessment of potential hydropower in Assam, India [19]. Correspondingly, the climate change impact on hydropower safety in Dak Nong, Vietnam, was carried out using SWAT [20].

Similarly, hydrological modelling of the Hoa Binh reservoir in Vietnam was conducted to optimize the utilization of flood control and hydropower generation. The results revealed a significant reduction in the peak flood downstream during the rainy season and a stable reservoir level during the dry season [21]. However, the hydrological model in the upper Mekong Basin identified a significant variation from the normal seasonal characteristic of river discharges since the hydropower began operation [22]. Moreover, water balance analysis in SWAT was used to quantify agricultural water demand for the sub-arid Mediterranean watershed [23]. SWAT modelling was carried out in a snowy area of Istanbul, in both Asia and Europe, to evaluate the water budgets of water resources in the context of uncertainties due to climate change and population growth in urbanized areas [24]. Land use change was evaluated using SWAT by simulating the streamflow of Murchison Bay in Uganda to further estimate sediment yield and nutrient loss for water resource management [25]. The groundwater analysis of the Taleghan Dam in Iran was also analysed by simulating the runoff river simulation using SWAT [26]. In addition, SWAT was used to demonstrate the importance of precipitation inputs as the main cause of uncertainties during the simulation of the Adige River basin in Italy, using multiple types of precipitation inputs [27]. Researchers found that monthly simulation produced better results than daily simulation in the ungauged Tonle Sap Basin in Cambodia [28]. Furthermore, SWAT was introduced to simulate runoff of the Mekong River to evaluate the hydrological application of tools in large basins [29].

In some parts of the Philippines, SWAT is utilized in different applications, such as for assessment of potential hydropower in the Visayas [30], Misamis Occidental [5], and the Agusan River basin [31], for predicting runoff in an ungauged watershed in Mabacan [32], and for simulating sediment yield in the Layawan watershed in Mindanao, to investigate land use change [33].

However, most of the applications in developing countries face a lack of precipitation and river discharge data, resulting in reliability issues in the validation process [34]. Therefore, we attempted to address how to overcome the data scarcity in the selected study area. Thus, this work used three types of precipitation datasets. Two are gridded datasets with a resolution of 0.5˚ latitude and 0.5˚ longitude, the NCDC-CPC dataset, and 1˚ latitude and 1˚ longitude, the GPCC dataset, as presented in Table 1. Hence, precipitation datasets were assigned to four stations to proportionally represent large areas of the MRB. Moreover, the calibration was conducted for 3 rivers: the Nituan, Libungan and Pulangi rivers. Then, validation of the calibrated models was conducted through the proxy basin to facilitate data scarcity. Finally, the SWAT interface was implemented in ArcGIS to carry out a watershed model of the MRB despite the limited hydrological datasets available for validation.

3.   Study Area

3.1. The Philippines and Mindanao

The Philippines is situated in Southeast Asia and includes three major zones: Luzon, Visayas and Mindanao, as shown in Figure 1(a). Mindanao is located in the southern Philippines and is the second-largest group of islands, after Luzon [35]. The Philippines had a population of 100,981,437 during the 2015 census [8] and includes a total of 7,107 islands with a land area of 300,000 km2 (Figure 1(b)) [36]. Moreover, the Philippines comprises 17 regions, 80 provinces, 143 cities, 1,491 municipalities, and 42,028 barangays (villages) [36]. Mindanao has a land area of 120,812 km2 and had a population of 24,135,775 during the 2015 census, as shown in Figure 1(c) [8], and it is subdivided into six administrative regions that further split into 27 provinces, 35 cities, and 422 municipalities [36].

The Philippines has 421 principal river basins and 18 major basins according to the National Water Regulatory Board (NWRB) [30]. Additionally, the Philippines has four types of climate, as defined by the spatial distribution of monthly rainfall, and experiences an average of 20 typhoons annually [7]. From 1990 to 2006, 520 disasters were induced by seven major natural hazards in this country, affecting 19,298,190 families (approximately 95 million people), who were repeatedly hit by natural hazards such as tropical cyclones, floods and landslides within the same period [37]. Considering these characteristics, water resource management is very challenging in the Philippines because of seasonal weather changes. Modelling is an alternative approach to account for the weather factors that influence the watershed.

3.2. Mindanao River Basin (MRB)

The MRB is the second-largest basin in the country [38], with a total area of 21,503 km2 [39]. It lies in four regions covering 72 municipalities and 1,732 villages in eight provinces, namely, Maguindanao, Lanao del Sur, Bukidnon, Sultan Kudarat, Davao del Sur, Davao del Norte, North Cotabato, and South Cotabato, as shown in Fig. 2 [39]. Due to the dependency on rain throughout the year, the MRB climate was classified as Types 3 and 4 under the modified Corona Climate Classification System of the PAGASA [4]. Moreover, the MRB includes major rivers, such as the Mindanao River and the Tamontaka River, which enters the sea of the lowest part of Cotabato City in Maguindanao [2]. The Pulangi River originates from Bukidnon Province. The Ambal-Simuay River has its waterhead in Lanao del Sur, and the Ala River navigates the Ala Valley in the south [2].

Moreover, the MRB is located at the coordinates of 124º47’35.71’ longitude and 7º12’17.06” latitude [4]. MRB has three vast marshes, namely, the Ligawasan, Ebpanan, and Libungan marshes, located within the central and lower parts of the basin. Thus, this large water resource in the MRB will be a potential asset to enhance hydropower development in support of the economic growth of the nearby regions. Therefore, the main reason to choose this study area is to maximize the application of the potential water resources for sustainable hydropower development in Mindanao.

Figure 2. (a) SAR-DEM for the MRB (10-m resolution), (b) land use and land cover map, (c) soil map and classification in the study area, and (d) slope classification in the HRU level

4.   Datasets

This study mainly used the available datasets from certain government agencies in the Philippines, as presented in Table 1. These datasets were requested from the corresponding listed agencies, but the available records are limited. The precipitation datasets from global gridded models were obtained online from the corresponding websites.

4.1. Digital Elevation Model

The synthetic aperture radar digital elevation model (SAR-DEM) with a 10-m resolution was obtained from the University of the Philippines Training Center for Applied Geodesy and Photogrammetry (UP-TCAGP) [30, 31]. These datasets were collected from point cloud data at a rate of 300 to 400 km2 per day at every sensor by the use of airborne light detection and ranging (LiDAR) technology and appended with SAR in some areas of concern [40]. This DEM was mostly used in the related previous studies because of its high resolution and accessibility [5, 38, 40]. The DEM was projected with the Universal Transverse Mercator (UTM) Zone 51 projection and World Geodetic System (WGS) 1984 as the horizontal datum, as displayed in Figure 2(a).

4.2. Administrative Boundaries

The Philippines administrative boundaries were obtained from a global administrative map and compared with other shapefiles from the Philippines GIS organization. The administrative shapefile was then projected with UTM 51N and WGS 1984, and then the MRB areas were overlaid onto the provincial boundaries, as shown in Figure 2(b), and municipal boundaries, as shown in Figure 2(c). The MRB lies in 72 municipalities of 6 provinces within 4 regions of Mindanao.

Table 1. Summary of the datasets used in the SWAT simulations of this study

4.3. Land Use and land Cover

Land use and land cover data from Landsat 8 of 2010 with a 30-m resolution were obtained from the National Mapping and Resource Information Authority (NAMRIA) and validated on the ground in 2015 by the agency. The shapefile of this dataset also used the same projection as the DEM. This dataset is among the main components of the model structure. Thus, the HRU was determined using this dataset, and the reclassification results are presented in Figure 2(d). The HRU results reported the following figures: the total area of the watershed is 2,041,449.74 ha, including the comprising agriculture area (52.65%), bushland (23.78%), open forest (7.84%), closed forest (5.71%), marshland (4.01%), grassland (4.37%), water (1.12%), and built-up area (0.53%). The large agricultural land of the MRB with an area of 1,074,869.37 ha comprises perennial crops and annual crops.

4.4. Soil Type and Slope

The soil map and local soil type classification was obtained from the Bureau of Soil and Water Management

(BMWS). Additionally, this soil classification was used in the earlier studies conducted in Mindanao [30, 31]. This information is provided by a shapefile projected in UTM 51N. The reclassification results for the HRU, as reflected in Figure 2(e), is characterized by mountain soil (39.04%), clay (25.45%), sandy loam (15.87%), clay loam (16.21%), loam (2.39%) and silt loam (1.04%). Mountain soil is a local name, according to the BSWM. Moreover, the reclassified slopes in the study area are divided into five categories: 0-25 (74.07%), 25-50 (19.48%), 50-75 (5.25%), 75-100 (0.99%) and 100-9999 (0.21%), as shown in Figure 2(f).

4.5. Weather Records

The weather records for 22 years, as shown in Figure 3, were obtained from DOST-PAGASA. Three datasets, temperature, humidity and wind speed, are available at four DOST-PAGASA stations in Cotabato, General Santos, Davao and Malaybalay, for the period from 1995 to 2017, but solar radiation is available at only the General Santos station for 2016-2017. The weather dataset is a main input for the SWAT model [14, 17]. Thus, these weather records were applied to simulate the MRB watershed model. The four stations are shown in Figure 1(d) and were used


Figure 3. Monthly observed values for (a) temperature, (b) humidity, and (c) wind speed at the four stations within and outside the MRB, and (d) solar radiation available at only the General Santos station for 2016-2017

4.6. Precipitation

The precipitation is a sensitive input in SWAT modelling because of the direct effects on the streamflow output [27, 41]. However, the study area has only 2 weather stations, the Cotabato and Malaybalay stations, located within the domain of the MRB. Two other weather stations, General Santos and Davao, are located outside the MRB. These four stations are not enough to represent the large MRB. Therefore, to address this concern, the datasets from the global gridded precipitation model from NCDC and GPCC were compared with the observational datasets from DOSTPAGASA [12]. Thus, the precipitation records from the abovementioned four DOST-PAGASA stations were then used in a comparison with datasets from the closest points of two global gridded precipitation models to evaluate the quality of the rainfall dataset. The multiple precipitation types were used to address the concern on lack of access to quality inputs in a developing country [34].

The NCDC-CPC precipitation is described as a global daily spatial coverage with a resolution of 0.5° latitude and 0.5° longitude covering 1979 to 2017 [42], while the GPCC precipitation [43], a global daily land surface precipitation with a resolution of 1° latitude and 1° longitude, covers the period from 1982 to 2016. The comparison results of the precipitations are summarized in Table 2.

Moreover, the three precipitation datasets were applied for simulating the watershed of MRB to improve the results of simulated discharges. The MRB was simulated from 2000 to 2017 and assigned 3 years of warm-up. Then, datasets from 2005 to 2010 were used to evaluate the precipitation responses against the simulated discharges during calibration and validation of the models based on the available period of river discharge records, as shown in Figure 4.


NCDC-CPC and GPCC datasets at the a) General Santos, b) Cotabato, c) Davao, and d) Malaybalay stations. Note that the Cotabato and Malaybalay stations are located within the MRB, while the other two stations, General Santos and Davao, are outside the MRB (see Fig. 1(d))

4.7. River Discharge

According to the Bureau of Standards of the Department of Public Works and Highways (DPWH), these records of river discharges were acquired through the information communication centre (ICT) with technical assistance of the United States Agency for International Development (USAID). The data were collected in five regions, which include three stations in the study domain: the Nituan River from 2005 to 2010, the Libungan River from 2006 to 2008, and the Pulangi River from 2009 to 2010. These river discharge records are used in the calibration and validation processes of watershed modelling.

Table 2. Comparison of the precipitation data between the observations (DOST-PAGASA) and global gridded datasets (NCDC-CPC and GPCC) at four stations within and near the MRB in terms of statistical indices such as correlation coefficient (R), index of agreement (d), root mean square error (RMSE) Stations

5.   Method

As described, SWAT was designed for agricultural, non-point source pollution and runoff river flow research. However, it has many features; for example, it can model stream flow by validating the simulated discharge from measured discharges. The hydrological cycle in SWAT is controlled by the water balance equation presented in Eq. (1). Thus, this water balance equation drives the physics of SWAT, allowing it to model the watershed of a certain basin [24, 44].

Where SWt is the final soil water content (mm H2O), SW0 is the initial soil water content (mm H2O), Rday is the amount of precipitation on day i (mm H2O), Qsurf is the amount of surface runoff on day i (mm H2O), Ea is the amount of evapotranspiration on i (mm H2O), Wseep is the amount of water entering the vadose zone from the soil profile on day i (mm H2O), and Qgw is the amount of return flow on day i (mm H2O).

5.1. Analysis Procedure

This study applied the following procedure to model the MRB watershed with SWAT. Each step of the procedure was clearly stipulated to ensure the process during the model simulation. Moreover, the analysis procedure is summarized in Figure 5 to provide a clear overview of the methods being applied with the SWAT model.

First, input weather datasets were prepared to match required input formats. The geospatial datasets such as the DEM, land use and land cover were processed and unmasked within the basin boundaries. Then, the watershed is delineated by using the processed inputs of the geospatial datasets. The processed land use and land cover were used to generate the HRUs by ArcSWAT. The processed weather dataset was loaded and utilized for the entire modelling process. Then, initial values of model parameters are set up which can be adjusted after calibrating simulation results in SWAT-CUP interface.

Then, calibrations were carried out to improve model performance by adjusting the parameters based on the sensitivity analysis from the SWAT-CUP (SUFI2) outputs. The watershed of the MRB was delineated into 107 subbasins, as shown in Figure 6. Then, the calibration period was carried out from 2005 to 2010, as reflected in Figure 4. Due to a limited number of river gauges in the study area, the only rivers with a record were the Nituan River, from 2005 to 2010 the Libungan River, from 2006 to 2008, and the downstream Pulangi River, from 2009 to 2010. Therefore, these rivers were utilized for calibration: sub-basin 28 was assigned to the Nituan River (Figure 7), subbasin 40 was assigned to the Libungan River (Figure 8), and sub-basin 45 was assigned to the downstream of the Pulangi River (Figure 9). Then, validations were carried out using the proxy basin principle [45] due to the relatively short records of river discharges. For instance, the calibrated/fitted parameters of River A (the Nituan River) were applied to River B (the Libungan River and Pulangi River) to validate the simulated discharges against the observed discharges. Since the record data of stream flow is not enough to split into two equal parts therefore the same datasets of the Nituan, Libungan, and Pulangi were also used in the proxy basin validation.

Table 3. Summary of the fitted parameters during calibration of the model by SUFI2 of SWAT-CUP. The sub-basins 28, 40, and 45 indicate the Nituan, Libungan, and Pulangi basins, respectively. r_parameter means modifying the parameters by multiplying the existing value to 1+ the given value, and v_parameter means using the given value to replace the parameter

(see K.C. Abbaspour, 2015)

Table 4. General criteria for performance evaluation in a statistical test for watershed scale

5.2. SWAT Calibration

The SWAT-CUP program was established to support the SWAT tool to minimize concerns about uncertainties [46][47]. Moreover, a sensitivity analysis was included inside the SWAT-CUP to tune the parameters according to the recommended results from a number of iterations [48]. Hence, the range of parameters was also enumerated in the absolute values section of SWAT-CUP [47]. Furthermore, SWAT-CUP also incorporates the statistical formulas for the Nash-Sutcliffe coefficient (NSE), coefficient of determination (R2) and percent bias (PBIAS), as presented in Appendix A, to evaluate the model performance. Thus, this study used the sequential uncertainty fitting version 2 (SUFI2) in the SWAT-CUP program for the calibration of the model performances for the Nituan, Libungan and Pulangi rivers.

The calibration was executed with 500 simulations in every iteration, and we conducted 5 iterations per recommendation in the previous studies [46, 47, 49]. Model calibration does not guarantee the improvement in the model performance. However, it helps modellers evaluate uncertainties by elucidating the sensitivity of parameters that can be adjusted with the observational datasets to improve the statistical results and model fitness [46]. Here, the calibration process was executed cautiously, and the 11 sensitivity parameters shown in Table 3 were investigated, but it was observed that the model performances always depend on the quality of the weather inputs, especially the precipitation, and model structure.

In addition, there are two ways to adjust the parameters: by manual calibration in the ArcSWAT itself or by the use of the SWAT-CUP tool, as previously described. Regardless of these options, the purpose of the calibration process is to improve the model fitness and statistical indicators to ensure the quality of the model outputs. Thus, this study used the general evaluation criteria recommended for the watershed to evaluate the model performances [49, 50], as presented in Table 4.

6.   Results

6.1. Calibration

The calibration results were measured through the statistical indices shown in Table 5, depicting an R2 of 0.61 at the Nituan River, 0.50 at the Libungan River and 0.42 at the Pulangi River for the DOST-PAGASA model. The NCDC-CPC model has R2 values of 0.66 at the Nituan River, 0.49 at the Libungan River and 0.55 at the Pulangi River. Additionally, the R2 values of the GPCC model are 0.62 at the Nituan River, 0.51 at the Libungan River and 0.27 at the Pulangi River. The PBIAS at the Nituan River of the DOST-PAGASA model is better compared with that of the NCDC-CPC model, with 13.70, and the GPCC, with 16.30. Then, the percentage of uncertainty in the models were estimated by p-factor and r-factor at the Nituan River with 0.42 and 0.50, respectively, for the DOST-PAGASA model; 0.33 and 0.50, respectively, for the NCDC-CPC model; and 0.54 and 0.74, respectively, for the GPCC model. The p-factor and r-factor at the Libungan River are 0.08 and 0.71, respectively, for the DOST-PAGASA model; 0.17 and 1.04, respectively, for the NCDC-CPC model; and 0.06 and 0.85, respectively, for the GPCC model. The p-factor and r-factor at the Pulangi River are 0.25 and 0.20, respectively, for the DOST-PAGASA model; 0.04 and 0.10, respectively, for the NCDC-CPC model; and 0.08 and 0.34, respectively, for the GPCC model. These results are obtained in the 95% percentage of uncertainty (PPU) in the simulated discharge model.  

A negative NSE means that the model performance is unsatisfactory and is characterized by extreme values [50]. A negative statistical performance indicates that the observed average streamflow is better than the simulated streamflow. The simulated discharges at the Libungan and Pulangi rivers are underestimated due to a lack of information on the reservoir management of dams and the water withdrawal from agricultural irrigation in the study area. The system of water storage, release and distribution has a very significant effect on the behaviour of river discharge in a watershed [51]. Then, this has an excessive impact on model performance, as shown in Figures 8 and 9 In addition, sub-basins 40 and 45 have a wetland (marshland) component, and the elevation difference between upstream and downstream is high. Therefore, the terrain of the MRB has a large heterogeneous component that may not be able to be accounted for during the modelling process. The SWAT application is very challenging in a largescale model and in wetlands. The wetlands normally absorb the surface and subsurface water between the inlet and outlet at several points, even if the inlet is well-defined [52]. Although SWAT already employs the basic concept of wetlands (marshlands), its ability to emulate the riparian wetland-river interaction is still under-studied [53]. Furthermore, most of the statistical model results do not show a significant difference between them.


Figure 6. The delineated watershed of the MRB; the watershed was divided into 107 sub-basins. The monitoring points are the links between the rivers/stream networks or junctions. Sub-basins 28, 40, and 45 indicate the Nituan, Libungan, and Pulangi rivers, respectively


Figure 8. Calibration results of river discharge at the Libungan River with three different precipitation inputs: (a) the


Results of this study are similar to the previous studies conducted in the Philippines due to lack of accessibility of datasets [34]. For instance, the SWAT modelling results in Pagsanjan-Lumban Basin in the Philippines show an R2 and NSE of 0.42 and 0.22, respectively, because of the lack of information on dams and paddy areas (Marshland) [51]. The SWAT water balance simulates the seasonal average discharges of 6.53 m3/s for DOST-PAGASA, 5.26 m3/s for NCDC-CPC, and 4.73 m3/s for GPCC in Nituan River. The seasonal average simulated discharge in Libungan River are 5.08 m3/s for DOST-PAGASA, 3.97 m3/s for NCDC-CPC, and 4.37 m3/s for GPCC. The seasonal average simulated discharges in Pulangi River are 261.69 m3/s for DOST-PAGASA, 215.63 m3/s for NCDC-CPC, and 250.45 m3/s for GPCC. The observed seasonal average discharges are 6.20 m3/s in Nituan, 11.37 m3/s in Libungan, and 470.64 m3/s in Pulangi. Among the simulated discharge models, the DOST-PAGASA model has closer values with the observed river discharges. Therefore, the DOST-PAGASA simulated discharges in the river mouth of the Mindanao river basin which is located in sub-basin 42 has an average simulated discharge of 502.03 m3/s, with the peak and lowest discharges of 1,239 m3/s and 2.29 m3/s, respectively.

Libungan River from 2006 to 2008, and (c) Pulangi River from 2009 to 2010

6.2. Validation: Proxy Basin

Access to enough data is a common problem in a developing country, and the Philippines is into exception [34]. Due to the challenges of limited access to river discharge records and a limited number of river gauges located in the study area, model validation was used to create a proxy basin. This principle is not new; it was introduced under the hierarchical scheme for systematic testing of hydrological simulation [45]. It was explained that streamflow at Basin C is to be selected as ungagged, and two basins within the region will be selected as gauged rivers, for example, Basins A and B. Then, the models will be calibrated in one of these basins and validated in another basin within the region. For instance, the model is calibrated in Basin A and validated in Basin B, or vice versa. In addition, the proxy basin is useful if the available streamflow record in a basin is not insufficient for equal split-sample and only if two validation results are acceptable and identical [45]. Moreover, the proxy basin was also used for model development of ungagged basins and at regional scales of watershed models [54]. The validation of the calibrated model was also carried out with SUFI2 of SWAT-CUP at the same rivers and in the same period as mentioned in the calibration section. However, the conventional way to split the available dataset is not applicable due to the insufficient length of the data record for river discharge and inconsistency of the duration, as shown in Figure 10. Therefore, from the proxy basin, the resulting parameters from the Nituan River were applied to the Libungan and Pulangi rivers, and vice versa, in this study to conduct the validation of the calibrated model and address the issues on data scarcity.

The validation results for the DOST-PAGASA model are similar to the R2 results of 0.61 and 0.61 at the Nituan River and 0.50 and 0.50 at the Libungan River during calibration and validation. In contrast, the model performance at the Pulangi River decreased from 0.42 to 0.33 for the DOST-PAGASA model, 0.55 to 0.40 for the NCDC-CPC model, and 0.27 and 0.21 for the GPCC model. Moreover, the PBIAS at the Nituan River did not change much, remaining at the good and satisfactory level according to the general ratings for the watershed provided in Table 4. In contrast, the NSE values for all the sub-basins remain unsatisfactory with negative values, as depicted in Table 6.  The simulated discharges at Libungan and Pulangi remain underestimated. In addition, the p-factor and r-factor of all the models did not change significantly from calibration to validation, even though proxy basins were applied. With these results, the calibration and validation satisfied only the Nituan River model. The Libungan and Pulangi rivers are the subject of more in-depth studies. Thus, the model performance will improve only if the dam management and irrigation water withdrawal will be accounted for in the model. Although datasets for dams and irrigation are not available; however, these results may be useful for understanding the role of dams and irrigation systems, especially downstream of rivers. In addition, the purpose of validation is to evaluate the fitness of the model after tuning the parameters during the calibration process. The proxy basin validation was carried out to evaluate the application of the SWAT model in a large basin with limited hydrological datasets. Therefore, regional modelling of large-scale basins with limited datasets does not accurately account for all the heterogeneous characteristics of the watershed.

Table 6. Summary of the statistical results for the validation of the river discharges of the models by proxy basin using the fitted parameters from the calibration results at sub-basins 28, 40, and 45, which indicate the Nituan, Libungan, and Pulangi basins, respectively


































































































7.   Discussions 

7.1. Effects of Precipitation on Simulated Discharge

Since river discharge is basically characterized by precipitation patterns, the observed precipitation inputs were carefully examined by evaluating the global gridded datasets and ground-truth dataset discussed in the previous precipitation section. The results of the correlation between the global gridded datasets and the observational dataset are 0.46 for the NCDC-CPC data and 0.78 for the GPCC data at General Santos, 0.90 for the NCDC-CPC data and 0.83 for the GPCC at Cotabato, 0.95 for the NCDC-CPC data and 0.63 for the GPCC data at Davao, and 0.92 for the NCDC-CPC for the GPCC and 0.52 for the GPCC for the GPCC at the Malaybalay station.

These strong correlations of the precipitation at the 3 stations indicate that the precipitation datasets have similar intensity and seasonal precipitation patterns in Mindanao. Thus, both global gridded precipitation models, either NCDC-CPC or GPCC, and DOST-PAGASA produce a trend of simulated discharges. However, among them, the precipitation model from DOST-PAGASA produced higher simulated discharges, as observed in the peak and base flow in Figure 10. In addition, the characteristics of simulated discharges physically influence the precipitation behaviours, as presented in Figures 7 to 9, respectively.

The three precipitation models proved that a high intensity of rainfall usually occurs in the months of June, October, and November, as reflected in the peak flow of the discharges in Figure 10. Moreover, the simulated peak discharge is slightly higher than the observed peak discharge in June for the Nituan and Libungan rivers, unlike that at the Pulangi River, which has a lower simulated peak discharge in June and November. Therefore, this situation might be affected by the dam system upstream of the Pulangi River. Thus, precipitation patterns for the Nituan River are valuable information for the nearest basin and for the upstream sub-basin to create a scenario for monthly seasonal discharge. However, for the Libungan and Pulangi rivers, the precipitation pattern is not beneficial for hydropower analysis because the simulated model is underestimated, and the area of interest is in wetlands with low elevations. Then, this terrain is not suitable for hydropower development; the ground is soft, and the elevation difference is not enough to generate power. Therefore, for the purpose of this study, utilization of the Nituan River will be suggested for the planning and development of hydropower in the study area.

7.2. Quality of Observation

As stated in the previous sections, the number of river gauges and river discharge records are very limited in MRB, and there are only 3 gauges with inconsistent records. Hence, the methods and types of instruments used for data collection were not mentioned in the source. This also affects the model performances; for instance, s shorter period of recorded data is most likely to produce low statistical index results, as observed for the Libungan and Pulangi rivers. Moreover, the ideal calibration and validation require enough river gauges data to split a dataset into two equal periods.

The quality of the river discharge has substantial effects on model calibration and validation. Aside from the length of records and the number of gauges available in the study domain, the method of collecting and processing the data are very serious factors to be considered. As explained previously, these datasets are part of a pilot project to support the technical capability of responsible agencies in certain regions. In short, a transitional process for transferring technology know-how from the service provider to the agency might lead to transitional development. Thus, the datasets duration is inconsistent among the stations, and the number of gauges is very limited despite the large size of this basin, with hundreds of streams and rivers. Therefore, as an alternative, using the regional watershed model and proxy basin process is an efficient way to implement model validation in the MRB.

8.   Conclusion

The main purpose of this study was to contribute to the improvement in water resource management for hydropower applications. Watershed modelling in the MRB in Philippines was carried out using SWAT. Since precipitation is a critical input that has a direct influence on river discharges, measured precipitation dataset within the the MRB is desirable as many as possible. However, only 2 DOST-PAGASA weather stations are available in MRB while another 2 stations are located close to MRB, despite the large area. Therefore, global precipitation gridded datasets from NCDC and GPCC were utilized to investigate the quality of precipitation in the MRB. Each precipitation dataset was individually used in SWAT simulations. Then, the simulated discharge was calibrated at 3 river gauges in the Nituan, Pulangi, and Libungan rivers using SWAT-CUP (SUF1). Moreover, due to limited short records from the 3 river gauges, the proxy basin process was applied for the validation of calibrated models of the same rivers. The models of the Nituan River provide better results compared with those of the Libungan and Pulangi rivers, even though a calibration was executed, and a proxy basin was also applied.

Lack of enough and qualified data is a common problem in a developing country. Due to the challenges of limited access to river discharge records and a limited number of river gauges located in the study area, model validation based on a proxy basin principle was applied in this study; for instance, the model was calibrated in Basin A and validated in nearby Basin B, or vice versa, to overcome the data scarcity in the study area. The proxy basin was also used for model development of ungagged basins at regional scales of watershed models. Therefore, this approach of calibration and validation of watershed modelling based on the proxy basin principal with various precipitation inputs demonstrates a method of watershed modelling for regions with insufficient precipitation and discharge data in developing countries.

The underestimated results of the Libungan and Pulangi rivers are mainly due to a lack of information on dams and irrigation water withdrawal. The study area has a vast wetland (marshland) and is characterized by a high elevation difference between the upstream and downstream, contributing to the uncertainty in the models and weak performance of the models in wetlands with limited datasets. Thus, these findings will be applicable to only the upstream side of the MRB for identifying potential sites for hydropower development.

In addition, comparison results for the precipitation inputs may be useful references to improve the meteorological measurements by adding additional weather stations in the regions. The trend of simulated discharges against precipitation inputs at 3 stations demonstrates the monthly seasonal characteristics of a watershed. Thus, this information can be used to determine which month might have a high potential for hydropower generation. Finally, this study specifically discussed the importance of meteorological agencies improving data collection and the application of the collected data in large areas; additionally, this study discussed the very significant role of river discharge, dam management and agriculture water withdrawal data for watershed analysis.

9.   Funding 

This research was partly supported by the Grant-in-Aid for Scientific Research (17K06577) from the Japan Society for the Promotion of Science (JSPS), Japan. The first author is supported by The Project for Human Resources Development Scholarship (JDS), Japan. 

10.  Conflicts of Interest

 The authors declare no conflict of interest. 

11.  References 

  • Ayson, L.G.; Tamang, J.T.; Barisco, C.A.; Sinocruz, M.O. Undersecretary, J.T. Tamang, E. Carmencita, A. Bariso, A. Director, M.O. Sinocruz, Philippine Energy Plan 2012-2030, Department of Energy. 223 (2010). Available online:
  • Bangsamoro Development Agency, Final Report Chapter 5. Existing Conditions of Flood and Disaster Management in Bangsamoro, 2014. Available online:
  • Department of Energy, Power Supply and Demand Highlights Total Non-coincidental Peak Demand, Philippines in MW, 2017. Available online:
  • University of the Philippines, Development of Climate-Responsive Integrated River Basin Master Plan, College of Forestry and Natural Resources, 2015. Available online: CC%20MRB_Inception%20Report.pdf.
  • Tarife, Rovick P., Anacita P. Tahud, Ellen Jane G. Gulben, Haroun Al Raschid Christopher P. Macalisang, and Ma. Teresa T. Ignacio. “Application of Geographic Information System (GIS) in Hydropower Resource Assessment: A Case Study in

Misamis Occidental, Philippines.” International Journal of Environmental Science and Development 8, no. 7 (2017): 507–511.


  • Bangsamoro Development Agency,   Bangsamoro            Development           Plan,        1st           ,          2015.       Available       online:
  • Otieno, J.A. Scenario study for flood hazard assessment in the lower Bicol floodplain Philippine using a 2D flood model : based on the 1988 flood event caused by typhoon Yonning : a case study for the flood hazard assessment, WP 4500 SLARIM and ITC research project, (2004) 109. Available online:
  • National Economic and Development Authority, Philippine Development Plan 2017-2022, 2017.
  • Ramos, B.T. Final Report on the Effects and Emergency Management re Tropical Storm “Sendong” (Washi), (2012) 1–35.
  • United Nation ESCAP, Assessment Report of the damages caused by Tropical Storm Washi, 2012.
  • Wendt, U.J. Typhoon Bopha and People Displacements in the Philippines, (2013) 33–46.
  • Cabrera, Jonathan, and Han Lee. “Impacts of Climate Change on Flood-Prone Areas in Davao Oriental, Philippines.” Water 10, no. 7 (July 4, 2018): 893. doi:10.3390/w10070893.
  • Novoa, Vanessa, Ramón Ahumada-Rudolph, Octavio Rojas, Katia Sáez, Francisco de la Barrera, and José Luis Arumí.

“Understanding Agricultural Water Footprint Variability to Improve Water Management in Chile.” Science of The Total Environment 670 (June 2019): 188–199. doi:10.1016/j.scitotenv.2019.03.127.

  • Neitsch, S. L., J. G. Arnold, J. R. Kiniry, and J. R. Williams. "Soil and Water Assessment (SWAT) Tool theoretical documentation version 2009. Texas Water Resources Institute." Texas AgriLife Research and USDA Agriculural Research Service, Temple, Texas, USA (2011): 1–647.
  • Arnold, J. G., and N. Fohrer. “SWAT2000: Current Capabilities and Research Opportunities in Applied Watershed Modelling.” Hydrological Processes 19, no. 3 (2005): 563–572. doi:10.1002/hyp.5611.
  • Adu, Joy, and Muthukrishna Vellaisamy Kumarasamy. “Assessing Non-Point Source Pollution Models: A Review.” Polish Journal of Environmental Studies 27, no. 5 (May 30, 2018): 1913–1922. doi:10.15244/pjoes/76497.
  • W. Gassman, M. R. Reyes, C. H. Green, and J. G. Arnold. “The Soil and Water Assessment Tool: Historical Development, Applications, and Future Research Directions.” Transactions of the ASABE 50, no. 4 (2007): 1211–1250. doi:10.13031/2013.23637.
  • Shanbor Kurbah, and Dr. Manoj Kumar Jain. “Rainfall-Runoff Modeling of a River Basin Using SWAT Model.” International Journal of Engineering Research And V6, no. 12 (December 28, 2017). doi:10.17577/ijertv6is120111.
  • Kusre, B.C., D.C. Baruah, P.K. Bordoloi, and S.C. Patra. “Assessment of Hydropower Potential Using GIS and Hydrological Modeling Technique in Kopili River Basin in Assam (India).” Applied Energy 87, no. 1 (January 2010): 298–309. doi:10.1016/j.apenergy.2009.07.019.
  • Bang, Ho Quoc, Nguyen Hong Quan, and Vo Le Phu. "Impacts of climate change on catchment flows and assessing its impacts on hydropower in Vietnam’s central highland region." Glob. Perspect. Geogr 1 (2013): 1-8.
  • Ngo, Long Le, Henrik Madsen, and Dan Rosbjerg. “Simulation and Optimisation Modelling Approach for Operation of the Hoa Binh Reservoir, Vietnam.” Journal of Hydrology 336, no. 3–4 (April 2007): 269–281. doi:10.1016/j.jhydrol.2007.01.003.
  • Räsänen, Timo A., Paradis Someth, Hannu Lauri, Jorma Koponen, Juha Sarkkula, and Matti Kummu. “Observed River Discharge Changes Due to Hydropower Operations in the Upper Mekong Basin.” Journal of Hydrology 545 (February 2017): 28–41. doi:10.1016/j.jhydrol.2016.12.023.
  • Rivas-Tabares, David, Ana M. Tarquis, Bárbara Willaarts, and Ángel De Miguel. “An Accurate Evaluation of Water Availability in Sub-Arid Mediterranean Watersheds through SWAT: Cega-Eresma-Adaja.” Agricultural Water Management 212 (February 2019): 211–225. doi:10.1016/j.agwat.2018.09.012.
  • Cuceloglu, Gokhan, Karim Abbaspour, and Izzet Ozturk. “Assessing the Water-Resources Potential of Istanbul by Using a Soil and Water Assessment Tool (SWAT) Hydrological Model.” Water 9, no. 10 (October 24, 2017): 814. doi:10.3390/w9100814.
  • Anaba, Listowel Abugri, Noble Banadda, Nicholas Kiggundu, Joshua Wanyama, Bernie Engel, and Daniel Moriasi.

“Application of SWAT to Assess the Effects of Land Use Change in the Murchison Bay Catchment in Uganda.” Computational Water, Energy, and Environmental Engineering 06, no. 01 (2017): 24–40. doi:10.4236/cweee.2017.61003.

  • Hosseini, M., M. S.M. Amin, A. M. Ghafouri, and M. R. Tabatabaei. “Application of Soil and Water Assessment Tools Model for Runoff Estimation.” American Journal of Applied Sciences 8, no. 5 (May 1, 2011): 486–494. doi:10.3844/ajassp.2011.486.494.
  • Tuo, Ye, Zheng Duan, Markus Disse, and Gabriele Chiogna. “Evaluation of Precipitation Input for SWAT Modeling in Alpine Catchment: A Case Study in the Adige River Basin (Italy).” Science of The Total Environment 573 (December 2016): 66–82. doi:10.1016/j.scitotenv.2016.08.034.
  • Ang, Raksmey, and Chantha Oeurng. “Simulating Streamflow in an Ungauged Catchment of Tonlesap Lake Basin in Cambodia Using Soil and Water Assessment Tool (SWAT) Model.” Water Science 32, no. 1 (April 2018): 89–101. doi:10.1016/j.wsj.2017.12.002.
  • Rossi, C.G.; Srinivasan, R.; Jirayoot, K.; Le Duc, T.; Souvannabouth, P.; Binh, N.; Gassman, P.W. Hydrologic evaluation of the lower mekong river basin with the soil and water assessment tool model, Int. Agric. Eng. J. 18 (2009) 1–13.
  • Jason, J.; Garcia, S.; Marie, A.; De La Serna, L.; Fesalbon, M.A; Silapan, J.R. Estimation of Hydropower Potential Energy Using Gis and Swat Hydrologic Model in Western Visayas, (2015) 1–8.
  • Cuasay, J.L.; Agno, G.C; D.A. Malonzo, K.M.; May Fesalbon, R.A; Inocencio,L.V; Rosario, M.O.; Ang, C. Evaluation of Climate Forecast System Reanalysis and local weather station data as input for run-of-river hydropower assessment in Agusan River Basin, Philippines, (2014).
  • Tolentino, Arlene B., and Victor B. Ella. “Assessment of SWAT Model Applicability and Performance for Predicting Surface Runoff in an Ungauged Watershed in the Philippines.” IAMURE International Journal of Ecology and Conservation 17, no. 1 (January 29, 2016). doi:10.7718/ijec.v17i1.1067.
  • Palao, Leo Kris M., Moises M. Dorado, Kharmina Paola A. Anit, and Rodel D. Lasco. “Using the Soil and Water Assessment Tool (SWAT) to Assess Material Transfer in the Layawan Watershed, Mindanao, Philippines and Its Implications on Payment for Ecosystem Services.” Journal of Sustainable Development 6, no. 6 (May 27, 2013). doi:10.5539/jsd.v6n6p73.
  • Tan, Mou Leong, Philip W. Gassman, Raghavan Srinivasan, Jeffrey G. Arnold, and XiaoYing Yang. “A Review of SWAT Studies in Southeast Asia: Applications, Challenges and Future Directions.” Water 11, no. 5 (May 1, 2019): 914. doi:10.3390/w11050914.
  • Ministry of Economy Trade and Industry, Study on Infrastructure Development in Mindanao, Philippines Final Report Ministry of Economy , Trade and Industry Commissioned to : Ernst & Young ShinNihon LLC Table of Contents, 2017. [36] Philippine Statistic Authority, The Philippines in figures: 2015, 2015. doi:ISSN-1655-2539.
  • Lindfield, M.; Singru, R.N. Republic of the Philippines national urban assessment, 2014.
  • University of the Philippines, Mindanao River: Dream Ground Survey Report, 2017.
  • Bangsamoro Development Agency, Comprehensive Capacity Development Project for the Bangsamoro Development Plan for the Bangsamoro Final Report Sec, 2016.
  • University of the Philippines, Buayan-Malungon River Basin: DREAM Flood Forecasting and Flood hazard Mapping, (2015) 76. Available online:
  • Eini, Mohammad Reza, Saman Javadi, Majid Delavar, José A.F. Monteiro, and Mohammad Darand. “High Accuracy of Precipitation Reanalyses Resulted in Good River Discharge Simulations in a Semi-Arid Basin.” Ecological Engineering 131 (June 2019): 107–119. doi:10.1016/j.ecoleng.2019.03.005.
  • Chen, Mingyue, Wei Shi, Pingping Xie, Viviane B. S. Silva, Vernon E. Kousky, R. Wayne Higgins, and John E. Janowiak.

“Assessing Objective Techniques for Gauge-Based Analyses of Global Daily Precipitation.” Journal of Geophysical Research 113, no. D4 (February 29, 2008): 1-13. doi:10.1029/2007jd009132.

  • Schneider, Udo, Tobias Fuchs, Anja Meyer-Christoffer, and Bruno Rudolf. "Global precipitation analysis products of the GPCC." Global Precipitation Climatology Centre (GPCC), DWD, Internet Publikation 112 (2008).
  • Winchell, M.; Srinivasan, R.; Di Luzio, M.; Arnold, J. ArcSWAT Interface For SWAT2012, User’s Guide, Texas Agrilife Res. United States Dep. Agric. Agric. Reseach Serv. (2013). Available online:
  • Klemeš, V. “Operational Testing of Hydrological Simulation Models.” Hydrological Sciences Journal 31, no. 1 (March 1986): 13–24. doi:10.1080/02626668609491024.
  • Abbaspour, Karim, Saeid Vaghefi, and Raghvan Srinivasan. “A Guideline for Successful Calibration and Uncertainty Analysis for Soil and Water Assessment: A Review of Papers from the 2016 International SWAT Conference.” Water 10, no. 1 (December 22, 2017): 6. doi:10.3390/w10010006.
  • Abbaspour, K.C. SWAT-CUP: SWAT Calibration and Uncertainty Programs- A User Manual, Department of Systems Analysis, Integrated Assessment and Modeling (SIAM), EAWAG. Swiss Federal Institute of Aqualtic Science and Technology, Duebendorf, Switzerland. User Manual (2015) 100p.
  • Abbaspour, K.C., E. Rouholahnejad, S. Vaghefi, R. Srinivasan, H. Yang, and B. Kløve. “A Continental-Scale Hydrology and Water Quality Model for Europe: Calibration and Uncertainty of a High-Resolution Large-Scale SWAT Model.” Journal of Hydrology 524 (May 2015): 733–752. doi:10.1016/j.jhydrol.2015.03.027.
  • Aqnouy, Mourad, Jamal Eddine Stitou El Messari, Hilal Ismail, Abdelmounim Bouadila, Jesús Gabriel Moreno Navarro,

Bounab Loubna, and Mohammed Reda Aoulad Mansour. "Assessment of the SWAT Model and the Parameters Affecting the Flow Simulation in the Watershed of Oued Laou (Northern Morocco)." Journal of Ecological Engineering 20, no. 4 (2019): 104–113.

  • Moriasi, D.SA; Gitau, M.W.; Pai, N.; Daggupati, P. “Hydrologic and Water Quality Models: Performance Measures and Evaluation Criteria.” Transactions of the ASABE 58, no. 6 (December 30, 2015): 1763–1785. doi:10.13031/trans.58.10715.
  • Ligaray, Mayzonee, Minjeong Kim, Sangsoo Baek, Jin-Sung Ra, Jong Chun, Yongeun Park, Laurie Boithias, Olivier Ribolzi,

Kangmin Chon, and Kyung Cho. “Modeling the Fate and Transport of Malathion in the Pagsanjan-Lumban Basin, Philippines.” Water 9, no. 7 (June 22, 2017): 451. doi:10.3390/w9070451.

  • Rezaeianzadeh, Mehdi, Latif Kalin, and Mohamed Hantush. “An Integrated Approach for Modeling Wetland Water Level: Application to a Headwater Wetland in Coastal Alabama, USA.” Water 10, no. 7 (July 2, 2018): 879. doi:10.3390/w10070879. [53] Rahman, Mohammed M., Julian R. Thompson, and Roger J. Flower. “An Enhanced SWAT Wetland Module to Quantify

Hydraulic Interactions Between Riparian Depressional Wetlands, Rivers and Aquifers.” Environmental Modelling & Software 84 (October 2016): 263–289. doi:10.1016/j.envsoft.2016.07.003.

  • Rientjes, T. H. M., B. U. J. Perera, A. T. Haile, P. Reggiani, and L. P. Muthuwatta. “Regionalisation for Lake Level Simulation – the Case of Lake Tana in the Upper Blue Nile, Ethiopia.” Hydrology and Earth System Sciences 15, no. 4 (April 8, 2011): 1167–1183. doi:10.5194/hess-15-1167-2011.
  • Legates, David R., and Gregory J. McCabe. “A Refined Index of Model Performance: a Rejoinder.” International Journal of Climatology 33, no. 4 (April 18, 2012): 1053–1056. doi:10.1002/joc.3487.
  • Lee, Han Soo. “Evaluation of WAVEWATCH III Performance with Wind Input and Dissipation Source Terms Using Wave

Buoy Measurements for October 2006 Along the East Korean Coast in the East Sea.” Ocean Engineering 100 (May 2015): 67– 82. doi:10.1016/j.oceaneng.2015.03.009.

  • Kumar, Nirmal, Sudhir Kumar Singh, Prashant K. Srivastava, and Boini Narsimlu. “SWAT Model Calibration and Uncertainty Analysis for Streamflow Prediction of the Tons River Basin, India, Using Sequential Uncertainty Fitting (SUFI-2) Algorithm.” Modeling Earth Systems and Environment 3, no. 1 (March 3, 2017). doi:10.1007/s40808-017-0306-z.
  • N. Moriasi, J. G. Arnold, M. W. Van Liew, R. L. Bingner, R. D. Harmel, and T. L. Veith. “Model Evaluation Guidelines for

Systematic Quantification of Accuracy in Watershed Simulations.” Transactions of the ASABE 50, no. 3 (2007): 885–900. doi:10.13031/2013.23153.

  • Gupta, Hoshin V., Harald Kling, Koray K. Yilmaz, and Guillermo F. Martinez. “Decomposition of the Mean Squared Error and NSE Performance Criteria: Implications for Improving Hydrological Modelling.” Journal of Hydrology 377, no. 1–2 (October 2009): 80–91. doi:10.1016/j.jhydrol.2009.08.003.
  • Narsimlu, Boini, Ashvin K. Gosain, Baghu R. Chahar, Sudhir Kumar Singh, and Prashant K. Srivastava. “SWAT Model Calibration and Uncertainty Analysis for Streamflow Prediction in the Kunwari River Basin, India, Using Sequential Uncertainty Fitting.” Environmental Processes 2, no. 1 (February 6, 2015): 79–95. doi:10.1007/s40710-015-0064-8.






























Appendix A: Equations of the Statistical Indices Used to Evaluate the Model Performance

The index of agreement between two variables is computed by the following formula [55, 56]:



d is the index of agreement of two variables predicted and observed

Pi is the predicted value in a sample

?? is the observed values in a sample  is the mean value of the observed samples

n is the number of observations

The correlation coefficient shows the strength of the relations between two variables by:



Rxy is the correlation coefficient of the linear relationship between two variables, such as x and y:

Xi is the value of the x-variable of a sample  is the mean value of the x-variables

?? is the value of the y-variable of a sample  is the mean value of the y-variables

The coefficient of determination between two model simulations and measured values is [47, 57]:

2     ∑??=1[(? ?,? − Ǭ?)(??,? − Ǭ?)]2                            (3)


R2 is the coefficient of determinants

?? is the measured discharge

Ǭm is the mean of the measured discharge

?? is the simulated discharge

Ǭ? is the average of the simulated discharge

The root mean square error (RMSE) is used to measure the absolute fitness between the observed and the modelled results:



RMSE is the root mean square error of the samples

Pi is the predicted values in a sample ?? is the observed values in a sample n is the number of observations

The Nash-Sutcliffe coefficient is calculated as follows [57, 58]:



NSE is the Nash-Sutcliffe coefficient

?? is the measured discharge

Ǭ? is the mean of the measured discharge

?? is the simulated discharge

PBIAS measures the model fitness in terms of average tendency, and a small value of PBIAS indicates a better model fitness. PBIAS is calculated by [47, 57]:



?? is the measured discharge ?? is the simulated discharge n is the number of observations

RSR is a standardization of the RMSE to measure how well the model results fit with the observed values, and a lower value of RSR indicates a better model fitness. RSR is calculated by the following equation [47, 57]:



?? is the measured discharge

Ǭ? is the mean of measured discharge

?? is the simulated discharge n is the number of observations

The Kling-Gupta efficiency (KGE) is used to examine the decomposition of the Nash-Sutcliffe efficiency, with a value close to 1 indicating a better performance. KGE is calculated by the following equation [47, 59]:

??? = 1 − √(? − 1)2 + (? − 1)2 + (? − 1)2                 (8)


r is a linear regression coefficient between the simulated and the observed data α is a ratio of standard deviation between the simulated and measured data (α β is a ratio of the means between simulated and measured data .

The r-factor is the thickness of the 95% predicted uncertainties (95PPU) [60]:



????,97.5% is the upper boundary of the 95PPU (95% predicted uncertainties)

????,2.5% is the lower boundary of the 95PPU (95% predicted uncertainties) ???? is the observed standard deviation.

Read 770 times Last modified on Friday, 01 May 2020 05:20

Login to post comments