HOME JOURNALS CONTACT

Journal of Environmental Science and Technology

Year: 2017 | Volume: 10 | Issue: 6 | Page No.: 320-324
DOI: 10.3923/jest.2017.320.324
Modeling of Extreme Rainfall Recurrence Time by Using Point Process Models
Nurtiti Sunusi , E.T. Herdiani and Nirwan

Abstract: Background and Objective: Period of extreme rainfall recurrence in certain location becomes a very interesting thing to be studied. In agriculture, uncertain emergence affects many things. The calendar of cropping pattern is closely related to the prediction of recurrence time of extreme rainfall. This study aimed to construct the mathematical model to find the time of extreme rainfall recurrence through point process modeling. Methodology: The assumption used in this study was the waiting time until the subsequent occurrence of subsequent rain follows the inverse Gaussian and the log normal distribution. To estimate the parameters model, the moment method was used. Results: The time of recurrence of subsequent extreme rainfall can be determined by calculating the length of waiting time since the last extreme rain event. Conclusion: The expected waiting time until the subsequent occurrence of the next rain, does not depends on the time difference since the appearance of the extreme last rainfall.

Fulltext PDF Fulltext HTML

How to cite this article
Nurtiti Sunusi, E.T. Herdiani and Nirwan , 2017. Modeling of Extreme Rainfall Recurrence Time by Using Point Process Models. Journal of Environmental Science and Technology, 10: 320-324.

Keywords: drought, inverse gaussian distribution, extreme rainfall recurrence, Point process modeling and log normal distribution

INTRODUCTION

Indonesia as an archipelago country became one of the most countries vulnerable to the negative impacts of climate change. As a country extremely vulnerable to climate change, Indonesia must adapt to these various environmental pressures. Indonesia has also become a front runner in the international limelight due to the country’s eagerness to reduce greenhouse emissions by roughly 26%1.

Due to global climate change, the territory of Indonesia is projected to experience changes in patterns and intensity rainfall that will increase the risk of flooding and drought in the season dry. This has an impact, such as prolonged drought, floods and the increasing frequency of extreme climatic events that affect health and community livelihoods. Risks of impact future losses will be suppressed if the projection of the weather and extreme climate can be done2. Extreme rainfall projection assessment in Indonesia expected can provide an overview of the expected extreme rainfall values. Forecasting was done to anticipate the adverse effects of extreme rainfall events. An important entity for the analysis of duration (include minimum inter event time)3, used to capture duration dependence is hazard rate which gives the instantaneous exit probability conditional on survival4,5.

Research on modeling of rain phenomenon has been studied by previous researchers4,6-9. Several special investigators studied the determination of time-lapse distributions between events, i.e., by Aldersley et al.2 and McShane et al.7. Although, various models of renewal have been studied by previous researchers, it has not been determined which model can best be used in the forecast of future extreme rainfall events. In this study, examined two types of time distribution between extreme rain events i.e., inverse Gaussian and log normal. These two distributions are usually used to describe non-negative positively skewed data. This corresponds to the phenomenon of extreme rain events where the data has extreme skewness properties. This study aims to construct the mathematical model to find the time of extreme rainfall recurrence through point process models.

MATERIALS AND METHODS

The study area was South Sulawesi during the period January-August, 2017. South Sulawesi province is located at 0°12’- 8° South latitude and 116° 48’-122°36’ East longitude. The total area is 62,482.54 km2. Administratively, the boundary of South Sulawesi province: North side with Central Sulawesi province, West side with Makassar Strait, East side with Bone Bay and South side with Flores Sea. The total area, 62,482.54 km2 (42% of the entire island of Sulawesi and 4.1% of the entire area of Indonesia). South Sulawesi has a strategic location in Eastern Indonesia enabling South Sulawesi to function as a service center for both Eastern Indonesia and international scale.

Variables and analytical procedures: The forecast analysis of the time of occurrence of extreme rainfall was done through the determination of waiting time until the occurrence of the next extreme rain. The assumption used in the process is the waiting time until the subsequent occurrence of subsequent rain follows the inverse Gaussian and the log normal model. To estimate the model parameters, the moment method was used.

RESULTS AND DISCUSSION

In this section it was expected to display the time lapse of the upcoming extreme rain by calculating the estimated waiting time until the next rainfall occurs using conditional probabilities. The concept of conditional probability has been introduced by previous researchers, including Gigerenzer and Hoffrage10, Sanfilippo et al.11 and Ferraes12. A review of the recurrence time in various phenomena has been done by some previous researchers using various time distribution models between events, such as Weibull, Pareto, Brownian Passage model and Normal log13. In the study it was shown that all types of time distribution between events used to model reappearance time gave the same result i.e., (the waiting time until the occurrence of the next extreme rain events) does not depend on the time difference t0 since the last extreme rainfall event12.

Consider Rt as the time of the last extreme rain event at one location, t0 is the time difference since the last extreme rain event. Furthermore, the probability of occurrence of extreme precipitation events as follows: If the time lapse t0 since the last extreme event will occur at intervals (t0, t0+Δt0). Suppose τ the waiting time until the occurrence of the next extreme rain events. This interval is a random variable T>0 the distribution F(t0) = P(T<t0).

Furthermore:

(1)

as extreme rainfall system reliability with F(t0) is probability that extreme rainfall system is broken off until time t0 and R(t0) represent probability of extreme rainfall system function at t0. Conditional probability:

(2)

is probability that extreme rainfall system will fail until time τ. It’s clear that:

Thus:

(3)

The differential of Eq. 3 to τ is defined in Eq. 4:

(4)

where, mf(τ|T>t0)dτ is probability that extreme rainfall system will occur in time interval (t0, t0+Δt0) with assumption that there is no extreme rainfall until time t0. Here described the distribution of waiting time model until the appearance of extreme rainfall next using inverse Gaussian model and Gamma model.

Inverse Gaussian model: Let the probability density function of inverse Gaussian of random variable T is defined as:

(5)

Then mean:

(6)

Cumulative distribution for random variable inverse Gaussian as follows:

(7)

So, based on Eq. 4, we have:

(8)

Since the cumulative distribution function for the random variable Gaussian inverse distributed in Eq. 7 is a constant C. Thus, the probability density function of extreme rainfall occurrence14 can be written as follows:

(9)

With:

Furthermore, the waiting time τ until the occurrence of the next extreme rainfall that maximizes the conditional probability density function for extreme rainfall in Eq. 9 is:

(10)

The critical point for the Eq. 10 is the solution of:

or:

(11)

Based on the Eq. 11, the critical point is obtained as follows:

(12)

Since τ is positive, the solution of Eq. 11 is . Further, it checked the maximum conditional probability criteria.

View the Eq. 11 for interval [0, ∞). Due to f’(τ)>0 all τ at [0, ), then based on the monotonous theorem, f increased on [0, ). Due to f’(τ)>0 all τ in (, ∞), then f decreased on [ , ∞). So it can be concluded that f() maximum.

Using Eq. 11, it is possible to predict the occurrence of an upcoming extreme rainfall for the Gaussian inverse model. In the theory of density function of conditional probability, it is known that the predicted time of occurrence of extreme rainfall that will come through the determination of waiting time depends on the elapse of time (elapsed time) t0. However, Eq. 12 shows that the expected waiting time until the occurrence of extreme rainfall in the future does not depend on the time difference t0 since the last extreme rainfall event.

Gamma model: Let probability density function of Gamma of random variable T is:

(13)

Mean of T is and variance of T is .

Cumulative distribution of random variable Gamma is:

(14)

So, based on Eq. 4 and 14, we have:

(15)

Since the cumulative distribution function for the Gamma distributed random variable in Eq. 14 is a constant C. Thus, the probability density function of extreme rainfall recurrence can be written as follows:

(16)

with:

Furthermore, the waiting time until the occurrence of the next extreme rainfall that maximizes conditional probability density function for extreme rainfall in Eq. 15 is:

(17)

The critical points for Eq. 17 is solution of:

or:

(18)

The above equation is true if:

(19)

Based on Eq. 19, the critical point is obtained as follows:

(20)

If α<1, then maximum occurs outside the definition range and for this case the conditional probability density for the occurrence of extreme rainfall is a monotonous function down to zero for t0→∞. However, Eq. 20 shows that the expected waiting time until the occurrence of extreme rainfall in the future does not depend on the time difference t0 since the last extreme rainfall event.

CONCLUSION

Analysis of the occurrence of extreme rainfall is very important to determine the period of recurrence of extreme rain. Although many researchers have conducted research in this field, it has not found the most appropriate distribution to state the time between extreme rain events. In this study we have discussed the time between extreme rain events that distributed inverse Gaussian and Gamma. The results obtained show that the expected waiting time until the occurrence of extreme rainfall in the future does not depend on the time difference t0 since the last extreme rainfall event.

SIGNIFICANCE STATEMENTS

This study provides another way to determine the approximate time of recurrence of extreme rainfall using the concept of conditional probability. Based on the characteristics of extreme rainfall, this study provides significant input to the possibility of appropriate distribution for the time between two extreme events. Based on the formula obtained, both models (inverse Gaussian and Gamma) show that the chances of occurrence of the subsequent rainfall extremes do not depend on the time difference since the appearance of the last extreme.

ACKNOWLEDGMENT

The authors would like to thank the Ministry of Technology Research and Higher Education Indonesia for financial funding PUPT research grant 005/SP2H/LT/DRPM/IV/ 2017 through decentralization research.

REFERENCES

  • World Bank, 2010. Indonesia: Climate change. The World Bank, Washington, DC. http://go.worldbank.org/HQQFW5MV70.


  • Aldersley, A., S.J. Murray and S.E. Cornell, 2011. Global and regional analysis of climate and human drivers of wildfire. Sci. Total Environ., 409: 3472-3481.
    CrossRef    PubMed    Direct Link    


  • Nojumuddin, N.S., F. Yusof and Z. Yusop, 2016. Determination of minimum inter-event time for storm characterisation in Johor, Malaysia. J. Flood Risk Manage.
    CrossRef    


  • Jose, K.K. and B. Abraka, 2013. A counting process with gumbel inter-arrival times for modeling climate data. J. Environ. Statist., 4: 1-13.
    Direct Link    


  • Sunusi, N., S. Darwis, W. Triyoso and I.W. Mangku, 2010. Study of earthquake forecast through hazard rate analysis. Int. J. Applied Math. Stat., 17: 96-103.
    Direct Link    


  • Smith, J.A. and A.F. Karr, 1985. Statistical inference for point process models of rainfall. Water Resour. Res., 21: 73-79.
    CrossRef    Direct Link    


  • McShane, B., M. Adrian, E.T. Bradlow and P.S. Fader, 2008. Count models based on weibull interarrival times. J. Bus. Econ. Statist., 26: 369-378.
    Direct Link    


  • Caballero, W. and A. Rahman, 2013. Variability in rainfall temporal patterns: A case study for New South Wales, Australia. J. Hydrol. Environ. Res., 1: 41-48.
    Direct Link    


  • Ghassabi, Z., G.A. Kamali, A. Mahkati, S. Hajam and N. Javaheri, 2017. Analysis of the temporal distribution of rainfall storms in South West of Iran. Iran Water Res. J., 9: 99-107.
    Direct Link    


  • Gigerenzer, G. and U. Hoffrage, 1995. How to improve Bayesian reasoning without instruction: Frequency formats. Psychol. Rev., 102: 684-704.
    CrossRef    Direct Link    


  • Sanfilippo, P.G., A.W. Hewitt and D.A. Mackey, 2017. The importance of conditional probability in diagnostic reasoning and clinical decision making: A primer for the eye care practitioner. Ophthalmic Epidemiol., 24: 81-89.
    PubMed    Direct Link    


  • Ferraes, S.G., 2003. The conditional probability of earthquake occurrence and the next large earthquake in Tokyo, Japan. J. Seismol., 7: 145-153.
    CrossRef    Direct Link    


  • Sunusi, N., S. Darwis, W. Triyoso and I. Mangku, 2008. The Brownian Passage Time (BPT) model for earthquake recurrence models. Far East J. Math. Sci., 29: 711-718.


  • Sunusi, N., E.T. Herdiani, F. Nawawi and Giarno, 2015. Exploratory analysis of rainfall occurrence in South Sulawesi region using spatial point process. J. Math. Statist., 11: 113-118.
    CrossRef    Direct Link    

  • © Science Alert. All Rights Reserved