2.1 Study area
The Magat River Basin has been chosen as the optimal location for this study. Designated as a forest reservation area through Proclamation 573 on June 26, 1969 (Elazegui & Combalicer, 2004), it holds critical watershed status in region 2 due to its role in maintaining ecological balance and its economic significance to the local population and the Philippines at large. Situated in the northern part of the Philippines, it spans major portions of Nueva Vizcaya, and parts of Quirino and Isabela provinces in the Cagayan Valley region, Philippines covering a total area of 5,156 square kilometers.
The climate within the Magat watershed is categorized as Type I and Type III according to Corona's classification. The western section falls under Type I climate, characterized by a dry season from December to May and a wet season from June to November. This section is exposed to the Southwest Monsoon, receiving a substantial amount of rainfall brought by tropical cyclones occurring from June to September. The eastern section, on the other hand, experiences a Type III climate. The type III climate season is characterized by a relatively dry period from November to April, but there is no pronounced maximum rain period (Tattao, 2010).
2.2 Description of the WEAP model
The Water Evaluation and Planning (WEAP) tool is a user-friendly application that employs a comprehensive approach to water resource planning and policy analysis. Utilizing the water balance principle to simulate hydrological processes (Asghar et al., 2019), this model demonstrates its versatility across individual river basins and complex basin systems, as seen in the study of Metobwa et al. in 2018. Offering a thorough evaluation of factors including hydrology, land use, hydrogeology, climate, water quality, and water allocation, the WEAP model functions as a conceptual model, enabling the representation of the physical system, as outlined by Li et al. in 2015.
Previously, the WEAP model has been proven successful in applications to agricultural and urban catchments worldwide in particular, simulating climate (Joyce et al., 2005; Mehta et al., 2013; Ougougdal et al., 2020), water supply and demands (Yao et al., 2021; Agarwal et al., 2018), and population growth (Arsiso et al., 2017).
2.3 Model development
This water supply and demand study generally followed the methodological framework shown in Fig. 1. The data on water resources which include hydrologic variables and water demand were gathered and used as inputs in the WEAP model. Considering this local data, the WEAP model was calibrated and validated until it satisfactorily mimics the hydrological processes in the Magat River watersheds. The simulation of scenarios considering socio-economic, technological, climate changes, and forest loss was done to assess its impacts on water supply and demand in the area.
Furthermore, the hydrological model is semi-theoretical, continuous, semi-distributed, and deterministic. Given its semi-theoretical nature, the model requires calibration and verification. Notably, the WEAP system does not offer automatic calibration for the hydrological model, necessitating a manual implementation of the calibration process. The standard method calculates water demand by multiplying the activity level with the water use rate across various sectors, including industrial, municipal, and agricultural. This method applies to all sectors except for agriculture. To allocate resources among demands, WEAP utilizes a one-period linear programming routine, as outlined by Letcher et al. (2007) and Van Cauwenbergh et al. (2008).
2.4 Model calibration and validation
The parameters governing runoff generation from climate inputs underwent calibration and validation using historical streamflow observations from gauging stations in the Magat River Basin. Additionally, the calibration parameters of the WEAP model were manually adjusted, drawing on data collected, existing literature, and expert knowledge. Model performance evaluation utilized statistical indices, including the coefficient of determination (R2), Nash-Sutcliffe Efficiency (NSE), and percent bias (PBIAS).
2.5 Data input and data collection
To simulate monthly hydrological processes in the WEAP hydrological model, meteorological/climate data are essential. This includes monthly information on precipitation, temperature, wind speed, and relative humidity, alongside Digital Elevation Model (DEM), land use, and streamflow data. The DEM and land use data were procured from the National Mapping and Resources Information Authority (NAMRIA). Meanwhile, climate and streamflow data were obtained from the Philippine Atmospheric, Geophysical, and Astronomical Services Administration (PAGASA) and the National Irrigation Administration (NIA), respectively. Demand-related data, encompassing population water use rates, water consumption, agricultural demand, and population growth rates, were sourced from the National Water Resources Board (NWRB) and the Philippine Statistics Authority (PSA).
2.6 Schematic description of Magat sub-basin WEAP model
The Magat River Basin WEAP model structure is divided into seven (7) hydrological catchments, namely Alimit, Subwatershed 2, Ibulao, Lamut, Subwatershed 1, Matuno, and Sta. Fe with corresponding percent areas of 14.82%, 10.83%, 15.13%, 7.64%, 14.61%, 14.70%, and 22.28% in the Magat Watershed, respectively. Specifically, there is one (1) municipality that may have been considered domestic, while one (1) is accounted for agricultural water use along the downstream of the Magat River Basin. Figure 2 shows the schematic illustration of the model at Magat River Basin.
2.7 Future scenarios development
Watershed-based water management represents a collaborative strategy that brings together all water stakeholders to formulate planned, concerted, and consensus-driven actions for the common good. This participatory approach involves individuals from diverse backgrounds, even when their interests may diverge. The WEAP system facilitates the creation of different current scenarios and their simulation in the future, addressing various "what if" questions, as discussed by Sieber & Purkey in 2015. These simulations can be compared with the existing situation. In addition to the reference scenario, four scenarios were developed for this study: the inclusion of additional Local Government Units (LGUs) as water users, changes in population growth rate, adoption of new irrigation techniques, and consideration of climate change. These scenarios aim to predict the future gap between water supply and demand. The identified management intervention options are designed to meet both current and future water demands at a minimal cost. The following section outlines the four scenarios developed:
Baseline scenario (BS). The initial baseline scenario developed represents the present water supply and demand conditions within the watershed. This scenario forms the foundation for a more in-depth analysis of the current state and for making comparisons with other simulated scenarios. The data from the current accounts, covering the years 2000 to 2020, serves as the foundational information for the model. From this, various scenarios are created to investigate possible changes in the system in the years beyond the current accounts year, including 2025, 2030, 2050, and 2080, effectively extending the assessment of the current situation into the future. Furthermore, the domestic water demand for the baseline year was based on the actual population of the LGU who tapped Magat River as a source of their domestic water. The water users in the reference scenario are the current water users of the water resources in the basin which include the NIA-Magat River Integrated Irrigation System service areas and the LGU of Alfonso Lista. Also, the industry, mining, and other LGUs were added in the reference scenario as water users of the groundwater resources.
Additional LGUs as water users. In this set of scenarios, the purpose of the Magat Dam for domestic water use was extended to LGU Santiago City outside the Magat River watershed. The baseline scenario for domestic water use is LGU-Alfonso Lista in Ifugao.
Population growth scenario. Under this scenario, the population growth was projected using the arithmetic method for 2025, 2030, 2050, and 2080. The population growth rate in Region 2 from 2000 to 2021 was 3.23 percent (PSA, 2020). There are four scenarios under this category, namely 2025, 2030, 2050, and 2080 population. The baseline scenario is the current domestic water usage of Alfonso Lista.
Irrigation system improvement scenario. The primary livelihood for the people in the Magat watershed is agriculture, but this fragile ecosystem is experiencing growing degradation due to a combination of several factors, including adverse climate changes that are already impacting local farmers. The impact of increased conveyance efficiency because of irrigation system improvement was evaluated at 95 percent. According to the Food and Agriculture Organization (FAO 1989), paving canals with concrete would achieve a conveyance efficiency of 95 percent. This technique allows for the improvement of conveyance efficiency due to the improvement measures. According to Alejo and Balderama (2021), the actual average conveyance efficiency of Magat is 76 percent. The water users from a certain LGU, specifically for irrigation and domestic purposes, are the same across scenarios. This is to allow assessment of the impact of improved conveyance efficiencies on water demand and unmet water demand. There were three scenarios used under conveyance efficiency, namely, 95 percent conveyance efficiency which is categorized as high, conveyance efficiency of 87 percent as moderate, and 80 percent conveyance efficiency as low. The alternate wetting and drying, with assumed water savings of 15 percent (low), 25 percent (moderate), and 30 percent (high) based on literature and the NIA Master Plan for 2020–2030, was also simulated under the irrigation system improvement scenario.
Climate change scenario. By using the Climate Information Risk Analysis Matrix (CLIRAM) tool, PAGASA was able to provide the downscaled projected changes in climate variables, particularly for rainfall and temperature in both mid (2036–2065) and late 21st century (2070–2099). Table 1 as shown below quantifies the projections on the rainfall and mean temperature. These were inputted into the WEAP model after it was considered acceptable.
Table 1
Projected changes in seasonal Rainfall and mean temperature mid (2036–2065) and late (2070–2099) 21st century.
Projected Changes in Seasonal Rainfall in the Mid-21st Century (2036–2065). |
Month | Moderate Emission (RCP4.5) | High Emission (RCP8.5) |
Lower Bound | Median Bound | Upper Bound | Lower Bound | Median Bound | Upper Bound |
DJF | 3.5 | 11.5 | 48.4 | -1.5 | 12.4 | 34.8 |
MAM | 1.1 | 10.3 | 23.5 | 1.1 | 7 | 17.3 |
JJA | -27.7 | -17.2 | 0.3 | -24.2 | -2.6 | 25.7 |
SON | -4.4 | 3 | 11.9 | -2.3 | 11 | 16.2 |
Projected Changes in Seasonal Mean Temperature in the Mid-21st Century (2036–2065). |
Month | Moderate Emission (RCP4.5) | High Emission (RCP8.5) |
Lower Bound | Median Bound | Upper Bound | Lower Bound | Median Bound | Upper Bound |
DJF | 1 | 1.2 | 1.5 | 1.1 | 1.6 | 1.8 |
MAM | 0.9 | 1.2 | 1.7 | 1.2 | 1.7 | 2.3 |
JJA | 1 | 1.3 | 2 | 1.3 | 1.6 | 2.5 |
SON | 1 | 1.1 | 1.9 | 1.3 | 1.6 | 2.3 |
Projected Changes in Seasonal Rainfall in the Late-21st Century (2070–2099). |
Month | Moderate Emission (RCP4.5) | High Emission (RCP8.5) |
Lower Bound | Median Bound | Upper Bound | Lower Bound | Median Bound | Upper Bound |
DJF | 0.2 | 16.5 | 41.7 | -10.2 | 32.3 | 51.1 |
MAM | -11.2 | 2.4 | 27.3 | -6.7 | 9.7 | 28.4 |
JJA | -26.8 | -12.9 | 6.2 | -31.2 | -21.6 | 15.2 |
SON | -5.2 | 4.4 | 12.2 | -5.1 | 5.7 | 24.7 |
Projected Changes in Seasonal Mean Temperature in the Late-21st Century (2070–2099). |
Month | Moderate Emission (RCP4.5) | High Emission (RCP8.5) |
Lower Bound | Median Bound | Upper Bound | Lower Bound | Median Bound | Upper Bound |
DJF | 1.1 | 1.5 | 2.2 | 1.9 | 2.9 | 3.5 |
MAM | 1.2 | 1.7 | 2.6 | 2.4 | 2.9 | 4 |
JJA | 1.4 | 1.6 | 2.7 | 2.8 | 3.3 | 4.5 |
SON | 1.3 | 1.5 | 2.7 | 2.6 | 3.1 | 4.3 |
Forest loss scenario. The land use change scenario suggests a notable increase in agricultural land in the basin. This change can be attributed to the decrease in grassland which was converted into agricultural land. The forested land increased from 90,516.48 ha to 170,490.93 ha from 2006 to 2010. This may be attributed to the reforestation initiatives, people-participated reforestation, and afforestation led by NIA-MARIIS DRD. The urban areas also increased five-fold, possibly due to the increase of urbanization in the provinces where the watershed is located. Likewise, the water bodies in the basin are also seen expanding from 5,629.39 ha to 9,415.75 ha.
Based on the land-use change from 2010 to 2015, there is an increase in agricultural land by 21 percent due to the conversion of grasslands and forested lands into agricultural areas. On the other hand, the forested lands and water bodies of the basin are observed to be declining. Both areas have a 13 percent decrease from 2010 to 2015. There are three scenarios under this category, namely, 2030, 2050, and 2080 forest loss scenario.