A Univariate Extreme Value Analysis and Change Point Detection of Monthly Discharge in Kali Kupang, Central Java, Indonesia

Sandy Herho - University of Maryland, 8000 Regents Dr #237, College Park, 20742, United States

Citation Format:

DOI: http://dx.doi.org/10.30630/joiv.6.4.953


Kali Kupang plays an important role in the life of the people of Pekalongan and its surrounding areas. However, until recently, not many hydrological studies have been carried out in this area. This study presents how Extreme Value Analysis (EVA) can predict future extreme hydrological events and how a dynamic-programming-based change point detection algorithm can detect the abrupt transition in discharge events variability. Using the annual block maxima, we can predict the upper extreme discharge probability from the generalized extreme value distribution (GEVD) that best fits the data by using the Markov Chain Monte Carlo (MCMC) algorithm as a distribution fitting method. Metropolis-Hasting (MH) algorithm with 500 walkers and 2,500 samples for each walker is used to generate random samples from the prior distribution. As a result, this discharge data can be categorized as a Gumbel distribution (  = 6.818,  = 3.456, and  = 0). The recurrence intervals (RI) for this discharge data can be calculated through this distribution. The changepoint location of the annual standard deviation of this discharge data in the mid-1990s is detected by using the pruned exact linear time (PELT) algorithm. Despite some shortcomings, this study can pave the way for using data-driven algorithms, along with more traditional numerical and descriptive approaches, to analyze hydrological time-series data in Indonesia. This is crucial, considering an increasing number of hydro climatological disasters in the future as a consequence of global climate change.


Hydrological extreme; change point detection; block maxima; MCMC; PELT.

Full Text:



R. N. Poerbatjaraka, Riwayat Indonesia Djilid I. Jakarta: Yayasan Pembangunan, 1952.

M. Dirhamsyah, Pekalongan Yang (Tak) Terlupakan. Pekalongan: KPAD Kota Pekalongan, 2015.

F. N. Fritsch and R. E. Carlson, “Monotone piecewise cubic interpolation,†SIAM Journal on Numerical Analysis, vol. 17, no. 2, pp. 238–246, 1980.

C. B. Moler, Numerical computing with MATLAB. Society for Industrial and Applied Mathematics, 2011.

W. McKinney, “Data Structures for Statistical Computing in Python,†in Proceedings of the 9th Python in Science Conference (S. van der Walt and J. Millman, eds.), pp. 56 – 61, 2010.

P. Virtanen, R. Gommers, T. E. Oliphant, M. Haberland, T. Reddy, D. Cournapeau, E. Burovski, P. Peterson, W. Weckesser, J. Bright, S. J. van der Walt, M. Brett, J. Wilson, K. J. Millman, N. Mayorov, A. R. J. Nelson, E. Jones, R. Kern, E. Larson, C. J. Carey, I˙. Polat, Y. Feng, E. W. Moore, J. VanderPlas, D. Laxalde, J. Perktold, R. Cimrman, I. Henriksen, E. A. Quintero, C. R. Harris, A. M. Archibald, A. H. Ribeiro, F. Pedregosa, P. van Mulbregt, and SciPy 1.0 Contributors, “SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python,†Nature Methods, vol. 17, pp. 261–272, 2020.

S. S. Shapiro and M. B. Wilk, “An analysis of variance test for normality (complete samples),†Biometrika, vol. 52, no. 3/4, pp. 591–611, 1965.

A. Ghasemi and S. Zahediasl, “Normality tests for statistical analysis: a guide for non-statisticians,†International journal of endocrinology and metabolism, vol. 10, no. 2, p. 486, 2012.

S. Coles, An Introduction to Statistical Modeling of Extreme Values. London: Springer London, 2001.

I. Ologhadien, “Study of unbiased plotting position formulae for the generalized extreme value (gev) distribution,†European Journal of Engineering and Technology Research, vol. 6, no. 4, pp. 94–99, 2021.

J. Sampaio and V. Costa, “Bayesian regional flood frequency analysis with gev hierarchical models under spatial dependency structures,†Hydrological Sciences Journal, vol. 66, no. 3, pp. 422–433, 2021.

H. Tabari, “Extreme value analysis dilemma for climate change impact assessment on global flood and extreme precipitation,†Journal of Hydrology, vol. 593, p. 125932, 2021.

E. Towler, D. Llewellyn, A. Prein, and E. Gilleland, “Extreme-value analysis for the characterization of extremes in water resources: A generalized workflow and case study on new mexico monsoon precipitation,†Weather and Climate Extremes, vol. 29, p. 100260, 2020.

F. De Paola, M. Giugni, F. Pugliese, A. Annis, and F. Nardi, “Gev parameter estimation and stationary vs. non-stationary analysis of extreme rainfall in african test cities,†Hydrology, vol. 5, no. 2, 2018.

Y. Liu, Y. Hao, Y. Fan, T. Wang, X. Huo, Y. Liu, and T.-C. J. Yeh, “A nonstationary extreme value distribution for analysing the cessation of karst spring discharge,†Hydrological Processes, vol. 28, no. 20, pp. 5251–5258, 2014.

G. Bocharov. (2021) pyextremes. [Online]. Available: https://georgebv.github.io/pyextremes/

I. Jeliazkov and E. Lee, “Mcmc perspectives on simulated likelihood estimation,†Advances in Econometrics, vol. 26, pp. 3–39, 12 2010.

E. Gilleland and R. W. Katz, “extremes 2.0: An extreme value analysis package in r,†Journal of Statistical Software, vol. 72, no. 8, pp. 1–39, 2016.

G. Casella and R. L. Berger, Statistical inference. Belmont, CA: Cengage Learning, 2021.

D. Gamerman and H. F. Lopes, Markov chain Monte Carlo: stochastic simulation for Bayesian inference. Boca Raton, FL: CRC press, 2006.

W. K. Hastings, “Monte Carlo sampling methods using Markov chains and their applications,†Biometrika, vol. 57, pp. 97–109, 04 1970.

D. Foreman-Mackey, D. W. Hogg, D. Lang, and J. Goodman, “emcee: The MCMC hammer,†Publications of the Astronomical Society of the Pacific, vol. 125, pp. 306–312, mar 2013.

J. Goodman and J. Weare, “Ensemble samplers with affine invariance,†Communications in applied mathematics and computational science, vol. 5, no. 1, pp. 65–80, 2010.

R. Killick, P. Fearnhead, and I. A. Eckley, “Optimal detection of changepoints with a linear computational cost,†Journal of the American Statistical Association, vol. 107, no. 500, pp. 1590–1598, 2012.

D. T. Myers, D. L. Ficklin, S. M. Robeson, R. P. Neupane, A. Botero-Acosta, and P. M. Avellaneda, “Choosing an arbitrary calibration period for hydrologic models: How much does it influence water balance simulations?,†Hydrological Processes, vol. 35, no. 2, p. e14045, 2021.

N. Singh and P. Chinnasamy, “Non-stationary flood frequency analysis and attribution of streamflow series: a case study of periyar river, india,†Hydrological Sciences Journal, vol. 66, no. 13, pp. 1866–1881, 2021.

R. V. Rocha and F. A. S. Filho, “Mapping abrupt streamflow shift in an abrupt climate shift through multiple change point methodologies: Brazil case study,†Hydrological Sciences Journal, vol. 65, no. 16, pp. 2783–2796, 2020.

K. R. Ryberg, G. A. Hodgkins, and R. W. Dudley “Change points in annual peak streamflows: Method comparisons and historical change points in the united states,†Journal of Hydrology, vol. 583, p. 124307, 2020.

M. Sakizadeh and L. Chua, “Environmental impact of karkheh dam in the southern part of iran on groundwater quality by intervention and trend analysis,†Environmental Monitoring and Assessment, vol. 192, no. 683, 2020.

P. Gogumalla, S. Dubey, S. Tripathi, and D. Chandniha, “Trend and change point detection of precipitation in urbanizing districts of uttarakhand in india,†Indian Journal of Science and Technology, vol. 7, pp. 1573–1582, 10 2014.

C. Truong, L. Oudre, and N. Vayatis, “Selective review of offline change point detection methods,†Signal Processing, vol. 167, p. 107299, 2020.

E. Aldrian and R. D. Susanto, “Identification of three dominant rainfall regions within indonesia and their relationship to sea surface temperature,†International Journal of Climatology, vol. 23, no. 12, pp. 1435–1452, 2003.

B. Masitho, “Dinamika politik pembangunan pada masa orde baru (studi tentang industrialisasi ketergantungan dan peran modal jepang),†PERSPEKTIF, vol. 3, no. 2, 2014.

R. Thombs and A. Jorgenson, “Manufacturing the urban rift: A cross-national study of urbanization, manufacturing, and co2 emissions, 2000-2013,†Human Ecology Review, vol. 25, 12 2019.

D. Pumo, E. Arnone, A. Francipane, D. Caracciolo, and L. Noto, “Potential implications of climate change and urbanization on watershed hydrology,†Journal of Hydrology, vol. 554, pp. 80–99, 2017.

E.-S. Chung, K. Park, and K. S. Lee, “The relative impacts of climate change and urbanization on the hydrological response of a korean urban watershed,†Hydrological Processes, vol. 25, no. 4, pp. 544–560, 2011.

J. Sheng and J. P. Wilson, “Watershed urbanization and changing flood behavior across the los angeles metropolitan region,†Natural Hazards, vol. 48, no. 1, pp. 41–57, 2009.

P. E. Kindermann, W. S. Brouwer, A. van Hamel, M. van Haren, R. P. Verboeket, G. F. Nane, H. Lakhe, R. Prajapati, and J. C. Davids, “Return level analysis of the hanumante river using structured expert judgment: A reconstruction of historical water levels,†Water, vol. 12, no. 11, p. 3229, 2020.

M. Bashirgonbad, N. A. Moghaddam, S. S. Khalighi, M. Mahdavi, E. Paque, and M. Lang, “Comparative study of frequency analysis and hydro-climatic methods for estimating maximum flood discharge (case study: Bakhtiary watershed),†2018.

F. E. Eregno, V. Nilsen, R. Seidu, and A. Heistad, “Evaluating the trend and extreme values of faecal indicator organisms in a raw water source: a potential approach for watershed management and optimizing water treatment practice,†Environmental Processes, vol. 1, no. 3, pp. 287–309, 2014.

C. Kabeja, R. Li, J. Guo, D. E. R. Rwatangabo, M. Manyifika, Z. Gao, Y. Wang, and Y. Zhang, “The impact of reforestation induced land cover change (1990–2017) on flood peak discharge using hec-hms hydrological model and satellite observations: a study in two mountain basins, china,†Water, vol. 12, no. 5, p. 1347, 2020.

A. Sarminingsih, A. Rezagama, et al., “Simulation of rainfall-runoff process using hec-hms model for garang watershed, semarang, indonesia,†in Journal of Physics: Conference Series, vol. 1217, p. 012134, IOP Publishing, 2019.

A. Saleh, R. Ghobad, and R. Noredin, “Evaluation of hec-hms methods in surface runoff simulation (case study: Kan watershed, iran),†Advances in Environmental Biology, pp. 1316–1322, 2011.

X. Chu and A. Steinman, “Event and continuous hydrologic modeling with hec-hms,†Journal of Irrigation and Drainage Engineering, vol. 135, no. 1, pp. 119–124, 2009.

J. Searcy and C. Hardison, “Double-mass curves. manual of hydrology—part 1—general surfacewater techniques, methods and practices of the geological survey,†US Geological Survey water-supply paper. Waterlow and Sons, 1962.

M. R. Zaghiyan, S. Eslamian, A. Gohari, and M. S. Ebrahimi, “Temporal correction of irregular observed intervals of groundwater level series using interpolation techniques,†Theoretical and Applied Climatology, vol. 145, no. 3, pp. 1027–1037, 2021.

E. L. Dan, M. Dîns¸oreanu, and R. C. Mures¸an, “Accuracy of six interpolation methods applied on pupil diameter data,†in 2020 IEEE International Conference on Automation, Quality and Testing, Robotics (AQTR), pp. 1–5, IEEE, 2020.

I. Azizan, S. A. B. A. Karim, and S. S. K. Raju, “Fitting rainfall data by using cubic spline interpolation,†in MATEC Web of Conferences, vol. 225, p. 05001, EDP Sciences, 2018.

R. T. Clarke, “Hydrological prediction in a non-stationary world,†Hydrology and Earth System Sciences, vol. 11, no. 1, pp. 408–414, 2007.

L. Cheng, A. AghaKouchak, E. Gilleland, and R. W. Katz, “Non-stationary extreme value analysis in a changing climate,†Climatic change, vol. 127, no. 2, pp. 353–369, 2014.

D. P. Solomatine, “Data-driven modelling: paradigm, methods, experiences,†in Proc. 5th international conference on hydroinformatics, pp. 1–5, 2002.

P. R. Shukla, J. Skeg, E. C. Buendia, V. Masson-Delmotte, H.-O. Pörtner, D. Roberts, P. Zhai, R. Slade, S. Connors, S. van Diemen, et al., “Climate change and land: an ipcc special report on climate change, desertification, land degradation, sustainable land management, food security, and greenhouse gas fluxes in terrestrial ecosystems,†2019.