This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract
This study was conducted to determine the effect of the climate change to the level of groundwater recharging. This research was conducted on the watershed of Bangga by using the Soil Water Balance of MockWyn-UB model. Input data compose of evapotranspiration, monthly rainfall, watershed area, canopy interception, heavy rain factor and the influence of climate change factors (rainfall and temperature). The conclusion of this study indicates that there is a decreasing trend in annual groundwater recharge observed from 1995 to 2011. The amount of groundwater recharge varied linearly with monthly rainfall and between 3% to 25% of the rainfall. This result implies that rain contributed more than groundwater recharge to runoff and evaporation and the groundwater recharge and Bangga River discharge depends largely on the rainfall. In order to increase the groundwater recharge in the study area, reforestation programmes should be intensified.
Groundwater as set forth in the Indonesian Government Regulation number 43 of 2008 [1] was defined as the water contained in the soil or rock layers below the soil surface. Groundwater is one of the sources of water supply to the lives of human beings and animals on earth. Accumulation and spread of groundwater is determined by various factors, such as rainfall, morphology and geology [2]. The average height of rainfall in Indonesia is between 1,000 mm to 6,000 m per year. There are certain areas that have a rainfall of less than 1,000 mm but spread is very limited and only 0.9% of the total area of Indonesia. Areas with rainfall are 1,000–1,500 mm per year only covers an area of less than 4%.
Global climate change that will be encountered in addition to increase or reduce rainfall in some regions, the increasing in air temperature, can also be associated with changes in weather patterns, wind patterns, humidity, and solar radiation. The decline in rainfall as input variables watershed due to irregularities global climate will affect groundwater recharging and seasonal dynamics. In general, the impact is very simple: The higher rainfall will produce greater groundwater recharging and declining rainfall will reduce the recharging of groundwater.
The evidence on climate change has been reported systematically by official sources, including: The Intergovernmental Panel on Climate Change (IPPC), the United Nations Framework Convention on Climate Change (UNFCCC) and the World Wide Fund (WWF) Indonesia. IPCC in the 3rd report states that the global average temperature is projected to rise to 1.4–5.8°C between 1990 and 2100. A scenario of climate change [3] predicted that the temperature would rise between 1.3°C and 4.6°C until 2100 with the trend of 0.1°C - 0.4°C per year.
According to Rao and Al-Wadany (1995): 1) Global climate change may occur because of the increase in CO2 and other gases in the atmosphere of active radiation. Temperature is expected to rise with the increase of CO2 and other gases. 2) Based on climate models, where rainfall and temperature shows significant changes in the future. Increased CO2 in the atmosphere and changes in forest cover is the main reason suggested for climate change. Changes in precipitation and temperature obviously affect the groundwater [4].
The result research from Tung and Haith (1995) states that the effect of global warming don’t happening only in the discharge flow in most watersheds but will vary from time and space. Because of seasonal changes in rainfall and temperature, river discharge may decline for several months and increases in others. Watershed in different locations will see different patterns of climate change, and watersheds with different physical characteristics can respond in different ways. The flow of river flow is strongly influenced by the recharging of groundwater [5].
The study of Yates and Strzepek (1998) states that changes in rainfall and temperature could have serious consequences on regional water resources throughout the basin. Hydrological model showed a strong response among discharge as the dependent variable with rain and evapotranspiration as an independent factor. The flow of river flow is strongly influenced by the recharging of groundwater [6].
The research from Fu et al. (2007) states that 30% increase in precipitation causes a 50% increase in flow when normal temperatures compared to only 20–30% increase in the flow if the annual average temperature of 1.5°C higher than normal. In contrast, a decrease in rainfall of 20% produces about 25–30% less flow when the temperature is normal but 45% reduce in the flow if the temperature is 1.5°C higher than normal. Thus, the increase and decrease in rainfall affect groundwater recharge [7].
According to research Lorena et al. (2010) states that in simulated river flow and evapotranspiration reflect different variations in rainfall. The positive trend of rainfall resulted in the increase of surface water and groundwater, whereas evapotranspiration only modestly affected. Comparison of the effects of rainfall trend towards surface water and groundwater shows that the increase in surface water is 3 times larger, implying ground water system has little effect on climate change [8].
Studies conducted in Indonesia on groundwater ranges in estimating the potential for groundwater, groundwater discharge capacity and utilization of groundwater. However, none of the above studies are trying to see the effects of climate change on groundwater recharge estimates. In this context, this study was conducted in order to determine the effect that may occur from climate change on the level of recharging of groundwater as an important parameter known to develop or to use groundwater without damaging the environment or negative effect. This study aims to answer the problems by investigating the effects of climate change on the groundwater recharge.
2. Material and Methods
2.1. Description of Study
This research was conducted in the catchment area of Bangga that is tributary Palu River. According to administrative located in village Bangga, Sigi District of Central Sulawesi of Indonesia. Catchment area of Bangga as geographically located between 01° 15′07″ - 01°21′30″ south latitude and 119°49′20″ - 119°56′05″ east longitude. Area of catchment Bangga is 65.90 km2 and the long of main river around 15.50 km. For more details, the location of the research is presented in the Fig. 1.
2.2. Model Description
2.2.1. Soil water balance model
The Soil Water Balance (SWB) model is part of the model MockWyn-UB [9] used to simulate groundwater recharging. This model considers the monthly hydro climatology components such as rainfall, evapotranspiration, infiltration coefficient and factors groundwater recession. Model of MockWyn-UB recommended using if the location of the study of climate changes. Detection of presence or absence of climate change can use the model Mann Kendall. Model inputs, parameters and outputs are presented in Table 1.
In this model, the rain that falls in the watershed is divided into three parts, namely in the forest, mixed farms and open land, using the equation:
(1)
(2)
(3)
Where: PDAS, PHT, PKC, PLT = total rain, rain forest, mixed farm and open land respectively; LDAS, LHT, LKC, LLT = total area, forest area, mix farm area and open land respectively.
Net rainfall based on land cover and vegetation canopy interception by using the results of research [10].
(4)
(5)
(6)
(7)
Potential evapotranspiration for each month is calculated from Penman Monteith method [9, 11, 12, 14, 15, 17]:
(8)
Actual evapotranspiration (ETa) is divided into two parts:
(9)
(10)
Difference between TPN with the monthly evapotranspiration ETo,
(11)
Accumulated potential water loss (APWL) is divided into two parts:
(12)
(13)
Soil Moisture (SM) is divided into two parts:
In the wet months or TPN > ETo, SM value for each month is equal to field capacity In the dry months or TPN < ETo, SM value is calculated by the equation
The volume of groundwater is calculated by the Eq. (18).
2.2.2. Climate change
Detection of climate change using a model Mann - Kendall [9, 11–13, 17–20]:
(19)
(20)
(21)
Where Xj and Xk is the data value of the data “j” and “k”, j > k.
After the detection of whether there is a trend of improvement or deterioration in the rain with the Mann-Kendall, so to determine the amount the trend to use methods of non-parametric Sen’s [21, 22] assuming the trend linear, the procedure is started Eq. (19)–(21). Both methods are combined so-called Makesens method
(22)
Where: Q is the slope and B is a constant.
To obtain a slope estimates Q in Eq. (22), it first needs to be calculated slope for all data with the equation:
(23)
where j > k.
If there is “n” value “ Xj” in a time series, it is obtained as N = n(n– 1)/2 slope estimation Qi . Sens slope estimate is the median of N values Qi . N value of Qi is ranked from small to large, with an estimated Sens is:
(24)
To obtain estimates of “B” in Eq. (22), the value of “n” data from the difference (Xi–Q.ti) is calculated. The median value is the estimated of “B”.
2.3. Data Collection and Analysis
2.3.1. Data collection and analysis for climate change
The data required to analyze climate change in the form of secondary data: 1) Data daily rain of station Up Bangga and Down Bangga available from 1980 to 2011; 2) Data climatology of Bora station available from 1980 to 2011; 3) Potential Evapotranspiration data available from 1980 to 2011 [9].
Fig. 3 presents the monthly rainfall in Bangga watershed for the observation period 1995 to 2011. During this period the average rain is 110 mm, the largest rainfall occurred in September 1995 of 345 mm and the smallest occurred in April 2002 amounted to 0.00 mm. It can be seen that a decline in the monthly rainfall trends over the period.
Fig. 4 presents the potential evapotranspiration for the observation period from 1995 to 2011. During this period the average rain is 120 mm/mon, the largest rainfall occurred in March 1996 of 160 mm/mon and the smallest occurred in February 2009 amounted to 92 mm/mon. It can be seen that a decline in the monthly rainfall trends over the period.
2.3.2. Data collection and analysis for model
The data required to analyze the model in the form of secondary data and primary data. For secondary data are required: 1) Data monthly rain of station Up Bangga and Down Bangga available from 1995 to 2011); 2) Data temperature of Bora station available from 1995 to 2011); 3) Data Potential Evapotranspiration available from 1995 to 2011; 4) Maps satellite imagery; 5) Maps the earth in such a scale of 1:50,000 and 6) Soil type maps, scale 1:50,000. 7) The primary data obtained by directly sampling the soil in the research site to obtain the percentage fraction of land, SM content, soil bulk density and soil type conducted in accordance vegetation land cover (mixed garden and forest) to make a hole in the soil profile depth of the surface soil, topsoil and subsoil below (SM storage) [9].
3. Results and Discussion
The SWB model is part of the model MockWyn-UB which considers climate change in the form of rain correction factor (αP ) and temperature correction factor (αT) as input. Thus, before using this model needs to be analyzed in advance whether there is a climate change study site by using Mann-Kendall models. Furthermore, the projected changes are analyzed by the method of non-parametric Sen’s by Sutapa (2015) [9] there has been a climate change in Bangga watershed as mention in Table 2. Thus this model can be used SWB model calculation results are presented in Table 3.
Climate change is occurred if the values of Z test more or less than zero. Otherwise, if the value of Z test is equal to zero is not occurred climate change. A confidence value (α) used in the calculation of Mann-Kendall is 0.001; 0.01; 0.05 and 0.1. According to the table of normal standard “Z” the values are: Z0,001 = 3.292; Z0,01 = 2.576; Z0,05 = 1.96; Z0,1 = 1.645.
Signs of significance in the calculation of the Mann-Kendall categorized into five categories, namely: three star (***), two (**), one (*), the (+) and empty (..) which indicates the level of confidence (α) = 0.001; 0.01; 0.05 and 0.1 respectively. If there is no sign (blank) means a significant level (α) of more than 0.1 or can be said to be insignificant [9].
Fig. 5 presents the relationship between monthly rainfall (R) and potential evapotranspiration (ETo) with groundwater (GW) as a result of modeling taking into account climate change in year 1995. It can be seen that the potential evapotranspiration almost the same throughout the month, then the amount of ground water is only affected by the amount of monthly rainfall. This model may give good performances because of the variation of the ground water or linear follow the up and down of monthly rainfall. The mean effect of evapotranspiration (ETo) and rainfall (R) to be groundwater is 58% and 42%, where the influence of evapotraspiration ranged from 37% to 68%, while rainfall ranges between 32% to 54%.
Fig. 6 presents the results of the calculation model of groundwater in the period 1995 to 2011. At the beginning of the simulation in 1995 resulted in ground water is high enough then dropped dramatically in all months. In January (year 1995 to 2011) the value of the largest groundwater occurred in 2000, while in others nearly flat. In February and March the value of groundwater fluctuation but not too big. Period May–August (year 1995 to 2011) the largest fluctuation occurred in 2010. Finally, for the period September–December (year 1995 to 2011) the value of groundwater occurred in 1995, while for the others there is a fluctuation but not too big.
Fig. 7 presents variation of yearly groundwater for period 1995 to 2011. During this period the average groundwater is 147 mm/y, the largest groundwater is 676 mm/y in year 1995 and the smallest occurred in year 2006 amounted to 33 mm/y. It can be seen that there is a downward trend in groundwater from the beginning to the end of the simulation.
Fig. 8 and Fig. 9 present variation of yearly rainfall (R), evapotranspiration (ETo), groundwater (GW) and percentage of rainfall to be groundwater for period 1995 to 2011. During this period the average rainfall to be groundwater is 9.9%, the largest rainfall to be groundwater is 25.60% in year 1995 and the smallest occurred in year 2005 amounted to 3.24%. It can be seen that the fluctuation of the rain changes into groundwater.
From the results that have been presented can be seen that groundwater recharging period 1995 to 2011, in January decline of simulation start in 1995 until 1998, raised significantly in 2000 before finally almost the same until 2011. In February decreased very sharply from the beginning of the simulation until 1998, then slightly increased in 2000, 2007, 2009 and other years is almost the same. In March a very sharp decline from the beginning of the simulation until 1998 and then fluctuated until the end of the simulation (2011). In April there is a decrease until 1996, there is an increase in 1999, 2001, 2007, 2008 and 2009, while in others fluctuate. In May decline began early simulation and an increase in 1999, 2001 and 2007. In June decline began early simulation and an increase in 2000, 2007 and 2010. In July a decline from early simulation until 1997 and increase in 1998, 2008 and 2010. In August a decline from early simulation to 1997 and increased in 1998, 2007 and 2010. In September, October, November and December, there was a very sharp decline from the beginning of the simulation until 1996 and fluctuate for other years.
The annual simulation results indicate a trend decline in groundwater recharge from early 1995 until the end of the simulation in 2011. From 1998 to 2000, 2007, 2008, 2010 and 2011 were above the average annual groundwater for the period 1995 to 2011. The amount of groundwater is very linear with the amount of monthly rainfall is happening. It means that the higher the rainfall that occurred then recharging groundwater will be even greater, and vice versa.
Percentage rain to fill ground water is very small, ranging from 3% till 25%. This indicates that the rain that falls on the watershed Bangga to be a lot more runoff and evaporation. If observed Fig. 8, it is clear that in general the amount of annual rainfall is less than evaporation and no recharging groundwater. This is unlikely to happen. That is the amount of annual rainfall, evaporation annual and annual groundwater cannot be compared. What can be done is to compare a monthly basis in accordance with this simulation. It can be seen from Fig. 5 where the value is greater than the monthly rainfall evaporation and groundwater recharge occurs. This occurs and goes for the other years during the period 1995 to 2011.
4. Conclusions
The conclusion from this study is going on a downward trend in the annual groundwater recharge from early 1995 until the end of the simulation in 2011. For 1998 till 2000, 2007, 2008, 2010 and 2011 recharging groundwater is above the average annual groundwater for the simulation period (1995 to 2011). The amount of groundwater recharging linearly with the mount of monthly rainfall is happening. It means that the higher the rainfall that occurred then recharging groundwater will be even greater, and vice versa. Percentage of rain to fill ground water is very small, ranging from 3% to 25%. This indicates that the rain that falls on the Bangga watershed to be a lot more runoff and evaporation.
Acknowledgments
The author thanks very much to the River Regional Department Office of Sulawesi III in Palu and Watershed Management Department Office of Palu Poso in Center Sulawesi Province of Indonesia which prepared the hydro-climatology data and land use map for supporting this research. The author would also like to thank the Faculty of Engineering, University of Tadulako, Palu, Central Sulawesi of Indonesia for providing technical support.
Nomenclature
B
Constant
ETo
Potential evapotranpiration, mm/d
ETa
Actual evapotranspiration, mm/d
es
Saturated water vapor pressure, kPa
ea
Actual water vapor pressure, kPa
I
Infiltration
In
Coefficient of infiltration
k
Groundwater flow recession factor
LDAS,
Total area in watershed, km2
LHT
Forest area, km2
LKC
Mix farm area, km2
LLT
Open land area, km2
n
The amount of observational data
PDAS
Total rain in watershed, mm
PHT
Rain in forest, mm
PKC
Rain in mixed farm, mm
PLT
Rain in open land, mm
PNTHT
Net rain forest, mm
PNTKC
Net rain mixed farm, mm
PNTLT
Net rain open land, mm
Q
Slope
Rn
Net radiation, cal/cm2/d
S
Variants
TPN
Total net rain, mm
U2
Wind velocity at a height of 2 m above the ground, m/s
Vn
Volume of groundwater, mm/mon
Z
Statistical value
αP
Precipitation correction factor
αT
Temperature correction factor
σs
Root of Variants
Δ
The slope of the vapor pressure of water against temperature, kPa/°C
γ
psychrometric constant, kPa/°C
References
1. The Government of the Republic of Indonesia. The Indonesian Government Regulation number 43 of 2008 on Groundwater. Indonesia: 2008.
2. Bisri M. Groundwater. 1st edMalang, Indonesia: Univ. of Brawijaya Press (UB Press); 2012. p. 11–21.
3. WWF Indonesia and IPCC. Climate change in Indonesia, implications for humans and nature. Indonesia: 1999.
4. Rao A, Al-Wagdany A. Effects of climate change in Wabah River basin. J Hydrol Eng. 1995;121:207–215.
5. Tung CP, Haith DA. Global-warming effetcs on New York streamflows. J Water Resour Plann Manage. 1995;121:216–225.
6. Yates D, Strzepek K. Modeling the Nile basin under climate change. J Hydrol Eng. 1998;3:98–108.
7. Fu G, Barber M, Chen S. Impacts of climate change on regional hydrological regims in the Spokane River watershed. J Hydrol Eng. 2007;12:452–461.
8. Lorena L, Leonardo VN, Enrique RV, Goffreda LL. Basin-scale water resources assessment in Oklahoma under synthetic climate change scenario using a fully distributed hydrologic model. J Hydrol Eng. 2010;15:107–122.
9. Sutapa IW. Modeling discharge of Bangga watershed under climate change. Appl Mech Mater. 2015;776:133–138.
10. Dunne , Leopold . Interception. Lecture 3A of Geography. IHE Delft; The Netherlands: 1978.
11. Sutapa IW. Application model Mann-Kendall and Sen’s (Makesens) for detecting climate change. Infrast J Civil Eng Univ Tadulako. 2014;4:31–40.
12. Sutapa IW. Long-term trend climatology in Sigi, Central Sulawesi province. In : National Seminar on Civil Engineering; 28 February 2015; Narotama Univ., Surabaya. p. 267–277.
13. Sutapa IW. Application of non parametric test to detect trend rainfall in watershed Palu, Central Sulawesi, Indonesia. Int J Hydrol Sci Technol. 2016;6:238–253.
14. Hadisusanto Nugroho. Hydrology Application. 1st edMalang, Indonesia: Jogja Mediutama; 2011. p. 95–105.
15. Richard AG. Crop evapotranspiration-guidelines for computing crop water requirement-FAO irrigation and drainage paper No. 56. Rome: Food Agriculture Organization of the United Nation; 1998. chapter 2:
16. Nasution , Syaifullah Djazim. Spatial analysis of the North Coast regional drought index (North Coast) in West Java. J Air Indonesia. 2005;1:235–243.
17. Sutapa IW. Effect of climate change modeling to discharge [dissertation]. Malang, Indonesia: Univ. of Brawijaya; 2013. 113–114.
18. Onoz B, Bayazit M. The power of statistical test for trend detection. Turkish J Eng Env Sci. 2002;27:247–251.
19. Wang S, Zhang Z, Sun G, et al. Long-term streamflow response to climatic variability in the Loess Plateau, China. J Am Water Resour Assoc. 2008;44:1098–1107.
20. Hakan A, Savaş K, Osman S. Trend analysis of hydrometeorological parameters in climate regions of Turkey. In : BALWOIS 2010; Ohrid, Republic of Macedonia. 2010.
21. Salmi T, Määttä A, Anttila P, Ruoho-Airola T, Amnell T. Detecting trends of annual of atmospheric pollutants by the Mann-Kendall test and Sen’s slope estimates. Helsingki: Finnish Meteorological Institute; 2002.
22. Deo RC, McAlpines CA, Syktus J, McGowan HA, Phinn Q. On Australian heat waves: Time series analysis of extreme temperature event in Australia, 1950–2005. In : MODSIM 2007 International Congress on Modelling and Simulation; 10–13 December 2007; Canberra, Australia. The Modelling & Simulation Society of Australia and NZ Inc.; 2007. p. 626–636.
Fig. 1
Map of location.
Fig. 2
Flowchart of soil water balance model.
Fig. 3
Variation of monthly rainfall data in year 1995–2011.
Fig. 4
Variation of monthly evapotranspiration data in year 1995–2011.
Fig. 5
Relationship between GW, R versus ETo in year 1995.
Fig. 6
Variation of monthly groundwater January-December.
Fig. 7
Variation of yearly groundwater in period 1995–2011.
Fig. 8
Variation of yearly rainfall (R), evapotranspiration (ETo), groundwater (GW).
Fig. 9
Variation of percentage rainfall (R) to be groundwater (GW).
Fig. 10
Relationship between groundwater (GW) versus rainfall (R) and potential evapotranpiration (ETo) base line 1995–2011 and projection to 2031.
Table 1
Summary of SWB Model, Inputs, Parameters, Outputs
(1)
(2)
(a) Inputs Data
LDAS, LHT, LKC, LLT
Total area, forest area, area mix farm and open land area respectively (km2)