Available online at www.sciencedirect.com ScienceDirect

Procedia Environmental Sciences 24 ( 2015 ) 56 – 64

The 1st International Symposium on LAPAN-IPB Satellite for Food Security and Environmental Monitoring Monitoring and prediction of hydrological drought using a drought early warning system in Pemali-Comal river basin,

Waluyo Hatmokoa,c,*, Radhikaa, Bayu Raharjab, Daniel Tollenaarb, Ronald Vernimmenb

aResearch Center for Water Resources, Ministry of Public Works,Jl. Ir. H. Juanda 193, Bandung 40135, Indonesia bDeltares, Boussinesqweg 1, 2629 HV Delft, Netherlands cDoctoral Candidate, Parahyangan Catholic University, Jl. Ciumbuleuit No. 94, Bandung 40141, Indonesia

Abstract

Unlike meteorological drought, hydrological drought in Indonesia has not been routinely monitored. This paper discusses hydrological drought monitoring and prediction system in Pemali-Comal River Basin by utilizing a Drought Early Warning System (DEWS) based on Delft-FEWS software. The Standardized Runoff Index (SRI) is applied to river discharges data which is collected through a real-time telemetering system. In case of unavailability of river discharge data, TRMM satellite rainfall data is used to simulate river discharges. A preliminary prototype of DEWS shows that the characteristics of previous drought events can be evaluated, and the simulation shows the possibility of forecasting hydrological drought with time lag of 6 months. © 20152015 The The Authors. Authors. Published Published by Elsevier by Elsevier B.V ThisB.V. is an open access article under the CC BY-NC-ND license (Selectionhttp://creativecommons.org/licenses/by-nc-nd/4.0/ and peer-review under responsibility). of the LISAT-FSEM Symposium Committee. Selection and peer-review under responsibility of the LISAT-FSEM Symposium Committee Keywords: drought; hydrological drought; drought monitoring; drought prediction; Indonesia

1. Introduction

Drought is a natural disaster that threatens food security and having unique characteristics. It is a creeping phenomenon, difficult to identify at the beginning, the ending and the severity of drought events. Drought generally s classified into four categories which are meteorological drought, hydrological drought, agricultural drought, and socio-economic drought. Meteorological drought is defined as a lack of precipitation over a region for a period of time. Hydrological drought is related to a period with inadequate surface and subsurface water resources for

* Corresponding author. Tel.: +62-8122103205; fax: +62-22-2500163. E-mail address:[email protected], [email protected].

1878-0296 © 2015 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/). Selection and peer-review under responsibility of the LISAT-FSEM Symposium Committee doi: 10.1016/j.proenv.2015.03.009 Waluyo Hatmoko et al. / Procedia Environmental Sciences 24 ( 2015 ) 56 – 64 57

established water uses of a given water resources management system. Agricultural drought is associated with a period with declining soil moisture and consequent in crop failure. Finally, socio-economic drought is related to failure of water resources systems to meet the water demands and thus associating droughts with supply of and demand for an economic good [3]. To overcome the difficulties in evaluating the severity of drought events, identify its onset and ending, a drought index is applied. There are drought indices corresponding to meteorological drought, hydrological drought, and agricultural drought. Meteorological drought in Indonesia has been routinely monitored and predicted by BMKG using the Standardized Precipitation Index (SPI). For the purpose of water allocation for irrigation, public water supply, and other water users, hydrological drought indices usually perform better than meteorological drought indices, since these indices are directly related to discharges in the rivers and water levels in reservoirs or lakes. This paper discusses a hydrological drought monitoring and prediction system in Pemali-Comal River Basin, utilizing a Drought Early Warning System (DEWS) based on Delft-FEWS software, the same open shell system software used for meteorological drought index monitoring in Indonesia.

2. Methodology

2.1. Data and Location of the Study

The location of the study is at Pemali-Comal River Basin. Fig. 1 shows the location of Pemali Comal River Basin in Indonesia and in Central Province. The available data for drought indices computation is monthly river discharge at the irrigation weir from 1991 – 2013 for the six irrigation weirs in Pemali-Comal River Basin.

KOTA TEGAL KOTA PEKALONGAN   Sungapan     Kramat     Asemsiketek    PEMALANG  Sokawati  BATAN PEKALONGAN BREBES  TEGAL   Kaliwadas   Notog    Notog 





Fig. 1. Location of the study area, Pemali-Comal River Basin in Indonesia

The rivers and areal of the irrigation schemes of the six diversion structures is presented in Table 1. Diversion structure of the Notog Weir serves the largest irrigation area and is located in the district of Brebes. The smallest irrigation area is Asemsiketek in Pekalongan district.

Table 1. Description on irrigation diversion structures Diversion structures River Irrigation area (ha) District location of irrigation area Notog Weir Kali Pemali 26,952 Brebes Sungapan Weir Kali Waluh 7,277 Pemalang Sukowati Weir Kali Comal 9,005 Pemalang Kaliwadas Weir Kali Genteng 7,642 Pekalongan Asemsiketek Weir Kali Sengkarang 238 Pekalongan Kramat Weir Kali Sambong 1,176 Batang

58 Waluyo Hatmoko et al. / Procedia Environmental Sciences 24 ( 2015 ) 56 – 64

2.2. Hydrological Drought Indices

In order to enable identification of drought onset, duration and severity, a theory of run approach developed by Yevjevich [7] is applied. The drought indicator such as rainfall or discharge time-series Xt is truncated to a certain threshold level X0 that can be fixed or variable. Drought is defined at the time when the indicator value is below threshold level, or if the truncated value is negative, as presented in Fig. 2. The duration of a drought event is length of time between successive crosses of X0 when the truncated value is negative. Drought severity is the cumulative deviation from X0, represented by the area created below X0; and drought intensity is the average deviation from X0, or the ratio between drought severity and drought duration. The illustration in Fig. 2 shows three drought events. The highest drought severity is drought event 1; the longest drought duration is drought event 2; and event 3 represents drought with the highest intensity.

4.0

3.0

2.0

1.0

time 0.0 2 -1.0 1 Drought Index Drought 3

-2.0

-3.0 Severity = 5.0 Duration = 5 Intensity = 1.0 -4.0 Severity = 9.6 Severity = 7.5 Duration = 4 Duration = 2 Intensity = 2.4 Intensity= 3.75 -5.0

Fig. 2. The Theory of Run applied to the truncated series of drought indicator

Table 2. SPI and SRI values

Index Classification ≥ 2,00 Extremely wet 1,50 to 1,99 Very wet 1,00 to 1.49 Moderately wet -0,99 to 0,99 Near normal -1,00 to -1,49 Moderately dry -1,50 to -1,99 Very dry ≤ -2,00 Extremely dry

This Theory of Run approach is applied in many popular drought indices including Standardized Precipitation Index (SPI), as well as the Standardized Runoff Index (SRI), an application of SPI for river discharge. The SPI Index is calculated from several steps: 1) Compute of the average of 1, 3, 6, and 12 months of raw data, becoming Waluyo Hatmoko et al. / Procedia Environmental Sciences 24 ( 2015 ) 56 – 64 59 new time-series data of 1, 3, 6, and 12 months average; 2) For each time-series and each month, transform the data into Gamma Distribution; 3) Convert the time-series having Gamma distribution into standardized Normal Distribution having zero mean and variance of unity. This new series having standardized Normal Distribution is an SPI time-series corresponding to 1, 3, 6, or 12 months; 4) Analyse the SPI series according to the Theory of Run, identify the drought severities as negative accumulated departure from zero, drought durations, and drought intensities. The usual SPI drought level classifications according to WMO [6] and also applied to SRI are as follows.

2.3. Delft-FEWS

Delft-FEWS is an open environment for the application of various modelling tools build-up around a central database. A set of standard tools related to data handling is available. This includes modules for importing and exporting data, data validation, data interpolation (both filling gaps in time, as well as in space) and transformation of series (aggregation, disaggregation, and transformation). Delft-FEWS makes use of a so-called general adapter (Fig. 3), which allows a user to connect to its own models. Model input data such as rainfall or input discharges time series can be exported from Delft-FEWS and model results such as runoff, water levels and discharges time series can be imported in Delft-FEWS. For hydrological and hydraulic models, model wrappers are available that convert model result files to FEWS files. The TRMM data is part of the import time series data, and also part of pre- processing model input; the 3 hourly data are aggregated to daily time steps.

Fig. 3. (a) Concept of Delft-FEWS; (b) General Adapters in Delft-FEWS

Delft-FEWS is mostly applied for real-time forecasting purposes. The overall workflow structure for such forecasting application is composed of a set of workflows, each composed of one or more processes (see Fig. 4). The workflows are typically executed independently, with different time intervals. Fall back mechanisms can be applied to adjust the strategy in case one or more activities fail to provide data [1]. Delft-FEWS is also being used for several researches and management studies in water domain, as well as operational flood forecasting and warning system. J-FEWS is an example of the implementation of Delft-FEWS in City flood forecasting and warning system [2]. Application of the Delft-FEWS for drought early warning system in Indonesia had been developed in the framework of Joint Cooperation Programme (JCP) among research centres in Indonesia and in Netherlands, including Deltares, KNMI, Research Centre for Water Resources (Pusair) and BMKG.

3. Results and Discussions

3.1. Evaluation of the Previous Droughts

Time-series of monthly river discharges data from 6 irrigation weirs for the year of 1991 until 2013 were collected, screening for the outliers, and hydrological drought index of 1, 3, 6, and 12 months SRI-1, SRI-3, SRI-6, and SRI-12 were computed. The time series of SRI-1 at Notog Weir on Fig. 5a shows the previous extreme drought events of 60 Waluyo Hatmoko et al. / Procedia Environmental Sciences 24 ( 2015 ) 56 – 64

SRI-1 below -2, in the year of 1994 and 1997, and moderate drought events of SRI-1 below -1 in the year of 1993, 1996, 1999, 2000, 2002, 2003, 2006, 2007, 2008, and 2012. The longer average of hydrological drought indicator SRI-6 in Fig. 5b shows more stable fluctuation. It clearly identified drought in 1994 as the most extreme drought, followed by drought in the year of 1997, 2007 and 2008.

Fig. 4. Typical Delft FEWS forecasting workflow

3 3

2 2 -1 -6

1 1

0 0

-1 -1

-2 -2

-3 -3

Standardized Runoff Index SRI Standardized Runoff -4 Index SRI Standardized Runoff -4

-5 -5 1991 1993 1995 1997 1999 2001 2003 2005 2007 2009 2011 2013 1991 1993 1995 1997 1999 2001 2003 2005 2007 2009 2011 2013

Fig. 5. (a) SRI-1 at Notog Weir; (b) SRI-6 at Notog Weir

From this monthly drought index, annual characteristic of drought for each year: annual drought severity, duration and intensity of the drought can be evaluated, and compared with the field data on areal affected by drought. The case for Notog Weir, the annual drought severity based on both meteorological drought index SPI and hydrological drought index SRI is evaluated and together with the data on areal of rice field affected by drought is presented in Fig. 6. This figure shows that annual drought severity based on SRI correlates to the area affected by drought better than the meteorological drought SPI. The complete trend of correlation between drought severity and the data on areal affected by drought presented in Fig. 7 proves that hydrological drought index SRI performs much better than meteorological drought index SPI. The peak performance of the hydrological drought index SRI is at SRI-6 having Waluyo Hatmoko et al. / Procedia Environmental Sciences 24 ( 2015 ) 56 – 64 61

coefficient correlation of 94%, while the best performance of meteorological drought index SPI is at SPI-3 having correlation coefficient of only 39%.

25.00 6000

5000 20.00

4000

15.00

3000

10.00

Annual severity Annual severity of drought 2000 Areal affected by drought (ha) drought affected by Areal

5.00 1000

0.00 0 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 SRI-1 (r = 86%) SPI-1 (r = 33%) Affected by drought

Fig. 6. Annual drought severity and data on areal affected by drought

100%

94% 92% 93% 91% 80% 86%

60%

40% 39% 33% 20% 27%

0% Correlation coefficient with the areal by drought affectedcoefficient the areal with Correlation -5% -8% -20% 136912 Months

SRI SPI

Fig. 7. Correlation coefficient between annual drought severity and the areal affected by drought 62 Waluyo Hatmoko et al. / Procedia Environmental Sciences 24 ( 2015 ) 56 – 64

3.2. Monitoring of Hydrological Drought

Monitoring of current hydrological drought can be made in the same way of flood early warning system. Fig. 8 presents map of the weirs having drought warning level of below the specified threshold drought level. This kind of warning level map enables the decision makers as well as the stakeholders to take appropriate actions for drought mitigation. Evaluation of the hydrological droughts in several irrigation weirs is shown in Fig. 9, facilitating users to compare the onset, duration, severity, of the drought in several irrigation weirs simultaneously.

Fig. 8. Drought warning level in the 6 weirs

The process applied in the DEWS is the following: 1) The TRMM data is first imported from file transfer protocol “ftp://trmmopen.gsfc.nasa.gov//pub/merged/mergeIRMicro/“, and the binary data are converted into a format which can be read by Delft-FEWS; 2) The 3 hourly data are aggregated to daily or monthly data; 3) Bias correction for TMPA 3B42RT (result of validation between TRMM and ground stations as in Vernimmen et al. [5] is applied on monthly data; 4) Depending on model requirements data can subsequently be disaggregated again to daily data in order to get bias corrected daily data; and 5) For SPI computation the monthly TRMM data are used. The advantages of using TRMM are: meteorological drought index SPI can be calculated using TRMM rainfall data in case of no rainfall data from the ground stations; and the SRI is computed from distributed rainfall-runoff model Wflow utilizing TRMM rainfall data input. The limitation of TRMM is that it only available since the year of 2002.

Waluyo Hatmoko et al. / Procedia Environmental Sciences 24 ( 2015 ) 56 – 64 63

Fig. 9. Evaluation of drought events in some weirs

3.3. Prediction of Hydrological Drought

Drought prediction for several months ahead is an important aspect in drought mitigation planning, providing predicted inflow for better reservoir operations, gives warning sign for the farmers to adjust cropping calendar, and inform the people to start saving the water. A drought indicator is a variable to identify and assess drought conditions, such as precipitations and river discharges. A drought trigger is a threshold value of the drought indicator that distinguishes a drought category, and determines when drought response actions should begin or end [4]. Predicting the forthcoming drought events is identical with predicting the trigger of drought indicator. Therefore one should find the drought indices having highest correlation with the drought impact, and then predicting the value of the indicator for the next months, utilizing the accumulated information until this moment. A drought prediction has been simulated for the Western part of Pemali-Comal River Basin. The drought index SRI-LN-6 MJJASO, a six-month standardized runoff index using Log-Normal statistical distribution for the months of May-June-July-August-September-October is predicted based on SRI-LN-3 MAM, a three-months standardized runoff index using Log-Normal statistical distribution for the month of March-April-May, using the regression equation. The result in the Fig. 10b shows that the prediction is close to the observed data, especially for the extreme drought events in the year of 1991, 1994, and 2003, where both the predicted and the data are below -1.

64 Waluyo Hatmoko et al. / Procedia Environmental Sciences 24 ( 2015 ) 56 – 64

4 3

3 2 2 1 1

0 0 SRI SRI

-1 -1 -2 -2 -3

-4 -3 1991 1993 1995 1997 1999 2001 2003 2005 2007 2009 1991 1993 1995 1997 1999 2001 2003 2005 2007 2009 2011 SRI-1 SRI-6 MJJASO Data MJJASO Predicted MJJASO Fig. 10. (a) SRI-LN-6 MJJASO as a drought trigger; (b) Prediction of SRI-LN-6 MJJASO

4. Concluding Remarks

Hydrological drought index plays important part in water allocation management because it usually performs better than meteorological drought index. This fact is proved from their coefficient of correlation of annual drought severity with the data on area affected by drought. A preliminary prototype of Hydrological Drought Early Warning System utilizing Delft-FEWS is being developed. The system is designed to access hydrological drought indicators for example river discharges and reservoir water level by means of file transfer facility; offline manual input; telemetering system; and satellite access such as TRMM. The rainfall data from TRMM can be utilized in drought monitoring and prediction in two different ways: as input for meteorological drought index SPI, or becoming input for rainfall-runoff model producing river discharges that can be converted into hydrological drought index SRI. In offline mode, Delft-FEWS store and retrieve the hydrological drought indicators database enable the analysis and evaluation of the severity, duration and intensity from previous drought events. Current river discharges in the irrigation weirs can be monitored using on-line mode of Delft-FEWS, as well as several hydrological drought indicators giving information on the updated information of the ongoing drought events and the predicted future drought. Finally, simulation of drought prediction shows the possibility of forecasting hydrological drought with lead time of 6 months.

References

1. Gijsbers P. Opportunities and limitations of DelftFEWS as a scientific workflow tool for environmental modelling. International Environmental Modelling and Software Society (iEMSs). International Congress on Environmental Modelling and Software; 2010. 2. Gnting S, Adidarma W. Jakarta Flood Early Warning System (J-FEWS), Workshop on MCCOE Radar Meteorology/Climatology in Indonesia; 2013. 3. Mishra AK, Singh VP. A review of drought concepts. Journal of Hydrology 2010; 391(1-2): 202–16. doi:10.1016/j.jhydrol.2010.07.012. 4. Steinemann A. Drought indicators and triggers: a stochastic approach to evaluation. Journal of the American Water Resources Association 2003. 5. Vernimmen RRE, Hooijer A, Aldrian E, van Dijk AIJM. Evaluation and bias correction of satellite rainfall data for drought monitoring in Indonesia. Hydrology and Earth System Sciences 2012; 16 (1): 133–46. doi:10.5194/hess-16-133-2012. 6. WMO. Standardized Precipitation Index User Guide; 2012. 7. Yevjevich V. An Objective Approach to Definitions and Investigations of Continental Hydrologic Droughts. Hydrology Papers Colorado. State University Fort Collins, Colorado, (August); 1967.