Time-series analysis and statistical forecasting of daily rainfall in Kupang, East Nusa Tenggara, Indonesia

Sandy Herho (1) , Gisma Firdaus (2)
(1) Department of Geology, College of Computer, Mathematical, and Natural Sciences, University of Maryland, College Park, MD, 20742, United States
(2) Indonesian Sugar Research Institute (P3GI), Pasuruan, East Java, 67126, Indonesia
Fulltext View | Download


This pilot study presents a novel statistical time-series approach for analyzing daily rainfall data in Kupang, East Nusa Tenggara, Indonesia. By using the piecewise cubic hermite interpolation algorithm, we succeeded in filling in the null values in the daily rainfall time series. We then analyzed the monthly average and its pattern using the continuous wavelet transform (CWT) algorithm, which shows the strong annual pattern of rainfall in this region. In addition, we use the rainfall anomaly index (RAI) function to standardize daily rainfall as an indicator of dry/wet conditions in this region. Then we also use the daily RAI time-series objects from 1978 to 2020 for modeling and predicting daily RAI over the next year. The result is the root mean squared error (RMSE) of 0.8424041040593219. This Prophet model is also able to capture the linear trend of increasing drought throughout the study time period and the annual pattern of wet/dry conditions which is in accordance with previous study by Aldrian and Susanto (2003).

Article Details

How to Cite
S. Herho and G. Firdaus, “Time-series analysis and statistical forecasting of daily rainfall in Kupang, East Nusa Tenggara, Indonesia”, Int. J. Data. Science., vol. 3, no. 1, pp. 25-32, Jun. 2022.


[1] 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.
[2] L. Uieda, D. Tian, W. J. Leong, M. Jones, W. Schlitzer, L. Toney, M. Grund, J. Yao, Y. Magen, K. Materna, T. Newton, A. Anant, M. Ziebarth, and P. Quinn, J. aand Wessel, “PyGMT: A Python interface for the Generic Mapping Tools,” Oct. 2021.
[3] F. N. Fritsch and R. E. Carlson, “Monotone piecewise cubic interpolation,” SIAM Journal on Numerical Analysis, vol. 17, no. 2, pp. 238–246, 1980.
[4] C. B. Moler, Numerical computing with MATLAB. Society for Industrial and Applied Mathematics, 2011.
[5] 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.
[6] A. Grinsted, J. C. Moore, and S. Jevrejeva, “Application of the cross wavelet transform and wavelet coherence to geophysical time series,” Nonlinear Processes in Geophysics, vol. 11, no. 5/6, pp. 561–566, 2004.
[7] C. Torrence and G. P. Compo, “A practical guide to wavelet analysis,” Bulletin of the American Meteorological Society, vol. 79, no. 1, pp. 61–78, 1998.
[8] N. Salehnia, A. Alizadeh, H. Sanaeinejad, A. Z. M. Bannayan, and G. Hoogenboom, “Estimation of meteorological drought indices based on agmerra precipitation data and station-observed precipitation data,” Journal of Arid Land, vol. 9, p. 797–809, 2017.
[9] L. V. Povoa and J. T. Nery, precintcon: Precipitation Intensity, Concentration and Anomaly Analysis, 2016. R package version 2.3.0.
[10] S. J. Taylor and B. Letham, “Forecasting at scale,” The American Statistician, vol. 72, no. 1, pp. 37–45, 2018.
[11] T. Hastie and R. Tibshirani, “Generalized additive models; some applications,” in Generalized Linear Models (R. Gilchrist, B. Francis, and J. Whittaker, eds.), (New York, NY), pp. 66–81, Springer US, 1985.