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.
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 countrys 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.,
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 τ. Its 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
View the Eq. 11 for interval [0, ∞). Due to f(τ)>0 all τ at [0,
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
Gamma model: Let probability density function of Gamma of random variable T is:
(13) |
Mean 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
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
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.