
Glioblastoma is one of the most dangerous tumors for patients in clinical practice at present, and since glioblastoma originates from the brain, it will have a serious impact on patients. Therefore, more effective clinical therapeutic targets are still needed at this stage. Kinesin family member 15 (KIF15) promotes proliferation in several cancers, but its effect on glioblastoma is unclear. In this study, differentially expressed gene analysis and network analysis were performed to identify critical genes affecting glioma progression. The samples were divided into a KIF15 high-expression group and KIF15 low-expression group, and the association between FIK15 expression level and clinical characteristics was summarized and analyzed by performing medical data analysis; the effect of KIF15 on glioblastoma cell proliferation was detected by employing colony formation and MTT assays. The effect of KIF15 on tumor growth in mice was determined. It was found that KIF15 was a potential gene affecting the progression of glioblastoma. In addition, KIF15 was highly expressed in glioblastoma tumor tissues, and KIF15 was correlated with tumor size, clinical stage and other clinical characteristics. After the KIF15 gene was knocked out, the proliferation ability of glioblastoma was significantly inhibited. KIF15 also contributed to the growth of glioblastoma tumors in mice. Therefore, we found KIF15 to be a promising clinical therapeutic target.
Citation: Leibo Wang, Xuebin Zhang, Jun Liu, Qingjun Liu. RETRACTED ARTICLE: Kinesin family member 15 can promote the proliferation of glioblastoma[J]. Mathematical Biosciences and Engineering, 2022, 19(8): 8259-8272. doi: 10.3934/mbe.2022384
[1] | Kaoutar MOUNIR, Isabelle LA JEUNESSE, Haykel SELLAMI, Abdessalam ELKHANCHOUFI . Spatiotemporal analysis of drought occurrence in the Ouergha catchment, Morocco. AIMS Environmental Science, 2023, 10(3): 398-423. doi: 10.3934/environsci.2023023 |
[2] | Zinabu A. Alemu, Emmanuel C. Dioha, Michael O. Dioha . Hydro-meteorological drought in Addis Ababa: A characterization study. AIMS Environmental Science, 2021, 8(2): 148-168. doi: 10.3934/environsci.2021011 |
[3] | Carolyn Payus, Siti Irbah Anuar, Fuei Pien Chee, Muhammad Izzuddin Rumaling, Agoes Soegianto . 2019 Southeast Asia Transboundary Haze and its Influence on Particulate Matter Variations: A Case Study in Kota Kinabalu, Sabah. AIMS Environmental Science, 2023, 10(4): 547-558. doi: 10.3934/environsci.2023031 |
[4] | Joseph A. Kazery, Dayakar P. Nittala, H. Anwar Ahmad . Meteorological effects on trace element solubility in Mississippi coastal wetlands. AIMS Environmental Science, 2019, 6(1): 1-13. doi: 10.3934/environsci.2019.1.1 |
[5] | Suhail Mahmud, Md Al Masum Bhuiyan, Nusrat Sarmin, Sanjida Elahee . Study of wind speed and relative humidity using stochastic technique in a semi-arid climate region. AIMS Environmental Science, 2020, 7(2): 156-173. doi: 10.3934/environsci.2020010 |
[6] | Tuan Syaripah Najihah, Mohd Hafiz Ibrahim, Nurul Amalina Mohd Zain, Rosimah Nulit, Puteri Edaroyati Megat Wahab . Activity of the oil palm seedlings exposed to a different rate of potassium fertilizer under water stress condition. AIMS Environmental Science, 2020, 7(1): 46-68. doi: 10.3934/environsci.2020004 |
[7] | Abebe Kebede Habtegebreal, Abebaw Bizuneh Alemu, U. Jaya Prakash Raju . Examining the Role of Quasi-biennial Oscillation on Rainfall patterns over Upper Blue Nile Basin of Ethiopia. AIMS Environmental Science, 2021, 8(3): 190-203. doi: 10.3934/environsci.2021013 |
[8] | John R. Christy . Examination of Extreme Rainfall Events in Two Regions of the United States since the 19th Century. AIMS Environmental Science, 2019, 6(2): 109-126. doi: 10.3934/environsci.2019.2.109 |
[9] | Khaled Hazaymeh, Quazi K. Hassan . Remote sensing of agricultural drought monitoring: A state of art review. AIMS Environmental Science, 2016, 3(4): 604-630. doi: 10.3934/environsci.2016.4.604 |
[10] | Volodymyr Ivanov, Viktor Stabnikov, Chen Hong Guo, Olena Stabnikova, Zubair Ahmed, In S. Kim, and Eng-Ban Shuy . Wastewater engineering applications of BioIronTech process based on the biogeochemical cycle of iron bioreduction and (bio)oxidation. AIMS Environmental Science, 2014, 1(2): 53-66. doi: 10.3934/environsci.2014.2.53 |
Glioblastoma is one of the most dangerous tumors for patients in clinical practice at present, and since glioblastoma originates from the brain, it will have a serious impact on patients. Therefore, more effective clinical therapeutic targets are still needed at this stage. Kinesin family member 15 (KIF15) promotes proliferation in several cancers, but its effect on glioblastoma is unclear. In this study, differentially expressed gene analysis and network analysis were performed to identify critical genes affecting glioma progression. The samples were divided into a KIF15 high-expression group and KIF15 low-expression group, and the association between FIK15 expression level and clinical characteristics was summarized and analyzed by performing medical data analysis; the effect of KIF15 on glioblastoma cell proliferation was detected by employing colony formation and MTT assays. The effect of KIF15 on tumor growth in mice was determined. It was found that KIF15 was a potential gene affecting the progression of glioblastoma. In addition, KIF15 was highly expressed in glioblastoma tumor tissues, and KIF15 was correlated with tumor size, clinical stage and other clinical characteristics. After the KIF15 gene was knocked out, the proliferation ability of glioblastoma was significantly inhibited. KIF15 also contributed to the growth of glioblastoma tumors in mice. Therefore, we found KIF15 to be a promising clinical therapeutic target.
The recent abnormal behavior of rainfall and air temperature has caused great impacts on the environment and human life. One of the most significant impacts is the unpredictable changes in seasonal patterns. In some areas, high rainfall will cause flooding; in others, low rainfall with high air temperatures will result in drought. Drought prediction is one of the biggest challenges for scientists and hydrologists, mainly due to its complex nature: these events are random and can fluctuate over time [1]. In predicting drought, the accuracy of rainfall data is crucial because rainfall is the main factor determining water availability in a region [2]. In addition, air temperature, which affects the rate of evaporation, also plays an important role in modeling and predicting drought phenomena [3,4,5].
Drought is one of the most devastating natural disasters, impacting water supply, agriculture, energy production, ecosystems, and society [6]. Drought has affected many parts of the world over the past few decades, such as in the Southeast United States [7], China [8,9], Brazil [10], and Pakistan [11]. Drought can be classified into four categories: Meteorological drought, which is the lack of rainfall to below-normal levels in a certain period of time, hydrological drought, the lack of water availability in and on the surface of the soil, agricultural drought, the reduced yield or agricultural production due to reduced water supply, and socioeconomic drought, related to demand and supply in a market for goods of economic value [12].
The measurement tool for drought is called the drought index, a single value that can describe the severity of drought. Meteorological drought indices that can be used to monitor drought conditions include the Palmer drought severity index (PDSI), which uses the water balance equation in the soil [13], and the standardized precipitation index (SPI) using the rainfall probability approach [14]. In recent years, new drought indices have been developed to improve the effectiveness of existing ones. One of these developments is the standardized precipitation evapotranspiration index (SPEI). The SPEI is a development of the SPI that only considers rainfall. SPEI adds potential evapotranspiration parameters in its calculation to describe drought better than relying on rainfall alone. This is a response to climate change and its effect on drought [15]. The time scales of SPEI calculation are the same as the SPI: The 1-month period is used for short-term drought recognition, the 3- and 6-month periods are used for seasonal drought recognition, the 12-month period is used for medium-term drought, and the 24- and 48-month periods are used for long-term drought assessment [14].
Drought monitoring using the 1-month SPEI drought index in Timor Island, East Nusa Tenggara, has shown that drought events in Kupang City spanned 94 months, with different intensity classifications: 63 months experienced "moderately dry" levels, 25 months had "severely dry" levels, and 6 months reached "extremely dry" levels. Meanwhile, in Kupang Regency, there were 93 months of drought, 62 months at the "moderately dry" level, 26 months at the "severely dry" level, and 5 months reaching the "extremely dry" level. South Central Timor Regency recorded 90 months of drought intensity, with 59 months at the "moderately dry" level, 25 months at the "severely dry" level, and 6 months reaching the "extremely dry" level. In North Central Timor Regency, 88 months of drought occurred, with 62 months at the "moderately dry" level, 20 months at the "severely dry" level, and 6 months at the "extremely dry" level. Malaka Regency recorded 95 months of drought, with 66 months at the "moderately dry" level, 25 months at the "severely dry" level, and 4 months at the "extremely dry" level. In Belu Regency, there were 87 months of drought, with 59 months at the "moderately dry" level, 23 months at the "severely dry" level, and 5 months at the "extremely dry" level.
The distribution of drought intensity on Timor Island, classified based on the 1-month SPEI drought level in each observation area, shows a very significant variation. In addition, there is a natural trend where more severe drought events, namely "extremely dry" and "very dry" events, tend to occur less frequently than less severe events, such as normal drought. In a stochastic process, this trend can be interpreted as a power law where the intensity of drought events decreases as their severity increases. The power law reflects that very severe drought events have a lower probability of occurrence than milder drought events. The large fluctuations in the tail of the power law distribution, which includes extreme but rare events, indicate that extreme events have a very low probability but can occur with very large intensity [16]. The power law process is one of the case models of the non-homogeneous Poisson process with the intensity function of the form $ \left(\frac{\beta }{\gamma }\right){\left(\frac{t}{\gamma }\right)}^{\left(\beta -1\right)} $ [17].
The non-homogeneous Poisson process is a commonly used model to model the number of events as a function of time [18]. Special case models of the non-homogeneous Poisson process have been widely used in various disciplines, including hydrometeorology, as shown in the study of Achcar et al. [19], where non-homogeneous Poisson process models, namely Weibull and Goel Okumoto with multiple variable points, were used to estimate the number of ozone levels exceeding the standard limit in Mexico City. Another study by Achar et al. [20] used the non-homogeneous Poisson process at the change point with the power law process model to analyze the drought period based on the SPI in Brazil. Ellahi et al. [21] used the non-homogeneous Poisson process model with a linear intensity function to assess the number of hydrological drought events using the SPI in Pakistan.
In addition, research on drought prediction using the SPEI index has also been carried out in many parts of the world. Ghasemi et al. [22] and Karbasi et al. [23] forecasted the SPEI 12 drought index in Iran; Dikshit et al. [24] predicted the size of drought using the SPEI on two different time scales (SPEI 1 and SPEI 3) in the New South Wales region, Australia; Affandy et al. [25] modeled and predicted meteorological drought measured by the SPEI with a time range of 1, 3, 6, and 12 months in Lamongan Regency, Indonesia.
The power law process can occur in various natural and artificial phenomena, covering several fields of science such as biology, economics, physics, chemistry, and computer science [26,27,28,29,30]. Statistical inference for the power law process is generally based on the maximum likelihood estimator (MLE) and its asymptotic properties. The MLE is used to find the parameter values in the power law process model that are most likely to yield the observed data [31]. A special characteristic of the power law process is the estimated value of the shape parameter $ \left(\beta \right) $, which can describe how an intensity can increase or decrease. If $ \beta > 1 $, the intensity of an event will increase; if $ \beta < 1 $, the intensity of an event will decrease; and if $ \beta = 1 $, the power law process reduces to a homogeneous Poisson process [17]. Two categories of data cases can be applied for parameter estimation in the power law process model: the time interval between events and the number of events observed in the specified interval [32]. In this study, the estimation method is applied to the second category, where the number of drought events based on 1-month SPEI in Timor Island is considered a random variable with a predetermined observed time interval. The fit test results using the acute Cramér-von Mises test showed that the intensity of drought events based on the 1-month SPEI on Timor Island fits the power law process model. The parameter estimation of the power law process using MLE in each observation area shows $ \left(\beta \right) > 1 $, with Kupang City having a value of 1.063; Kupang Regency 1.174; South Central Timor Regency 1.095; North Central Timor Regency 1.049; Malaka Regency 1.034, and Belu Regency 1.112. This indicates a possible increase in the intensity of drought events in these areas. Therefore, as a mitigation effort and early planning in the face of future drought events, this study aims to analyze short-term meteorological drought periods using the power law process to obtain an estimate of the duration of future drought events.
Figure 1 shows the area of study was conducted in the Timor Island region of East Nusa Tenggara (NTT), Indonesia, covering six regency/cities, namely Kupang City, Kupang Regency, South Central Timor Regency, North Central Timor Regency, Malaka Regency, and Belu Regency. Geographically, the six locations are located in the western part of the Timor Archipelago with coordinates of 9°14' N and 124°56' E. The data used are secondary in the form of monthly rainfall amounts and monthly average air temperatures obtained from NASA Power through the website https://power.larc.nasa.gov/data-access. Data were used as input data for the SPEI 1-month drought index calculation parameters These results were classified based on drought severity to obtain the intensity of drought events in each observation area. The intensity of drought events was measured based on the frequency of drought periods within a time span. The observation period was from January 1981 to December 2023, with a record length of 516 months.
The SPEI was designed to consider rainfall and potential evapotranspiration (PET) in determining drought. The SPEI drought index calculation is based on the deficit value between rainfall and PET [15]. PET can be calculated by the Thornthwaite method using average air temperature based on the following equation [33]:
$ PET = 16K{\left(\frac{10T}{I}\right)}^{\mathrm{m}} $ | (1) |
$ K $ is a correction factor based on the latitude position of the observation area, T is the monthly average air temperature $ \left(^ \circ {\rm{C}}\right) $, and I is the annual internal heat index obtained from the sum of $ i $ for 12 months in the following equation:
$ i = {\left(\frac{T}{5}\right)}^{\mathrm{1,514}} \; \text{and} \; I = \sum _{1}^{12}i $ | (2) |
The $ m $ is a coefficient that is a variable tied to $ I $ with the equation given by:
$ m = \mathrm{6, 75} \times {10}^{-7}\;\;{I}^{3}-\mathrm{7, 71}\times 1{0}^{-5}\;\;{I}^{2}+\mathrm{1,792}\times {10}^{-2}\;\;I+0.492 $ | (3) |
The deficit between rainfall and PET or climate water balance can be determined by the following equation:
$ {D}_{i} = {CH}_{i}-{\mathrm{P}\mathrm{E}\mathrm{T}}_{i} $ | (4) |
$ {D}_{i} $ is the value of climate water balance in month i, $ {CH}_{i} $ is the amount of rainfall in month $ i $, and $ PE{T}_{i} $ is PET in month $ i $ in mm. Next, the value of $ {D}_{i} $ is standardized based on the probability density function of the log-logistic distribution with three parameters to capture the deficit value, since it is likely that the moisture deficit in arid and semi-arid areas may be negative. For the two-parameter distribution as used in SPI, the variable D has a lower limit of zero $ (0 > D < \infty ) $, which means D can only take positive values, while for the three-parameter distribution used in SPEI, D can take values in the range $ (\gamma > D < \infty ) $, which means D can also take negative values [15]. The probability density function of the log-logistic distribution is given as:
$ f\left(D\right) = \frac{\beta }{\alpha }{\left(\frac{D-\gamma }{\alpha }\right)}^{\beta -1}{\left[1+{\left(\frac{D-\gamma }{\alpha }\right)}^{\beta }\right]}^{-2} $ | (5) |
The parameters $ \alpha $, $ \beta $, $ \gamma $ in the log-logistic distribution are calculated using the L-moment procedure. L-moment calculation of Pearson Ⅲ distribution parameters can be obtained through the following equation [34]:
$ \beta = \frac{2{W}_{1}-{W}_{0}}{6{W}_{1}-{W}_{0}-6{W}_{2}} $ | (6) |
$ \alpha = \frac{({W}_{0}-2{W}_{1})\beta }{Г\left(1+\frac{1}{\beta }\right)Г\left(1-\frac{1}{\beta }\right)} $ | (7) |
$ \gamma = {W}_{0}-\alpha Г\left(\frac{1+1}{\beta }\right)Г\left(\frac{1-1}{\beta }\right) $ | (8) |
$ Г\left(\beta \right) $ is the gamma distribution function of $ \beta $. W is the probability weighted moments (PWMs) obtained from the following equation:
$ {W}_{s} = \frac{1}{N}{\sum }_{i = 1}^{N}{\left(1-{F}_{i}\right)}^{s}{D}_{i} $ | (9) |
$ s $ is the number of PWMs, and $ {F}_{i} $ is a frequency estimator that can be calculated using the equation given by:
$ {F}_{i} = \frac{i-\mathrm{0, 35}}{N} $ | (10) |
$ i $ is the range of observations arranged in ascending order and N is the number of data used. The probability function of the distribution D over various time scales can be calculated using the following equation:
$ F\left(D\right) = {\left[1+{\left(\frac{\alpha }{D-\gamma }\right)}^{\beta }\right]}^{-1} $ | (11) |
Based on the probability function, the SPEI can be calculated using the following equation [35]:
$ SPEI = \left(t-\frac{{c}_{0}+{c}_{1}W+{c}_{2}{W}^{2}}{1+{d}_{1}W+{d}_{2}{W}^{2}+{d}_{3}{W}^{3}}\right) $ | (12) |
$ W = \sqrt{-2\mathrm{ln}\left(P\right)} \;\text{ for P } \le 0.5 $ | (13) |
$ W = \sqrt{-2\mathrm{ln}\left(1-P\right)} \;\text{ for P } \ge 0.5 $ | (14) |
P is the probability of exceeding the value of D, which is determined by the following equation:
$ P = 1-F\left(x\right) $ | (15) |
The coefficient value of McKee is as follows:
$ \begin{array}{l} {c}_{0} = 2.515517 \;\;\;\;\;\; {d}_{1} = 1.432788\\ {c}_{1} = 0.802853 \;\;\;\;\;\; {d}_{2} = 0.189269\\ {c}_{2} = 0.010328 \;\;\;\;\;\; {d}_{3} = 0.001308 \end{array} $ |
Drought occurs when SPEI reaches drought intensity with SPEI value ≤ –1. The classification of SPEI drought index values is based on Table 1 [36].
SPEI value | Classification |
$ \ge $2.00 | Extremely wet |
1.50–1.99 | Very wet |
1.00–1.49 | Moderately wet |
(–0.99)–0.99 | Normal |
(–1.00)–(–1.49) | Moderately dry |
(–1.50)–(–1.99) | Severely dry |
$ \le $ –2.00 | Extremely dry |
The non-homogeneous Poisson process is a stochastic process used to count the number of events in a given time interval, where the rate of events is not constant but depends on time. A process counting $ \left\{N\left(t\right), t\ge 0\right\} $ is said [37] to be a non-homogeneous Poisson process with intensity function $ \lambda \left(t\right) $, $ t\ge 0 $, if:
a. $ N\left(0\right) = 0 $,
b. $ \left\{N\left(t\right), t\ge 0\right\} $ has independent increment,
c. P$ \left\{N\left(t+h\right)-N\left(t\right) = 1\right\} = \lambda \left(t\right)+o\left(h\right) $, and
d. P$ \left\{N\left(t+h\right)-N\left(t\right)\ge 2\right\} = o\left(h\right) $, where h > 0 and $ o\left(h\right) $ is a small number satisfying the condition $ \underset{h\to 0}{\mathrm{lim}}\frac{o\left(h\right)}{h} = 0. $
The expected value, also known as the cumulative function of the non-homogeneous Poisson process $ \left\{N\left(t\right), t\ge 0\right\} $ with intensity function $ \lambda \left(t\right) $ is defined as:
$ m\left(t\right) = {\int }_{0}^{t}\lambda \left(t\right)dt $ | (16) |
Based on Eq 16, the average estimate of $ N\left(t\right) $ is given by the equation:
$ \widehat{m}\left(t\right) = E\left(N\left(t\right)\right) = {\int }_{0}^{t}\widehat{\lambda }\left(t\right)dt $ | (17) |
where $ E\left(.\right) $ is the expectation value. $ \left\{N\left(t\right), t\ge 0\right\} $, modeled as a non-homogeneous Poisson process, is expressed as:
$ P\left({N}_{\left(t\right)} = n\right) = \frac{\left[\int_0^t {} \lambda \left(t\right)dt\right]}{n!}{exp}^{-\int_0^t {} \lambda \left(t\right)dt}\;\;\;\;\;\;n = \mathrm{0, 1}, 2\dots n $ | (18) |
Based on Eq 16, for $ t, s > 0, N\left(t+s\right)-N\left(t\right) $ has the following expected value function:
$ m\left(t+s\right)-m\left(t\right) = {\int }_{0}^{t+s}\lambda \left(t\right)dt $ | (19) |
Thus, if $ \left\{N\left(t+s\right), t, s\ge 0\right\} $, based on Eq 18, it can be modeled as a non-homogeneous Poisson process as follows:
$ P\left({N}_{\left(t+s\right)}-{N}_{\left(t\right)} = n\right) = \frac{{(m\left(t+s\right)-m(t\left)\right)}^{n}}{n!}{exp}^{(-(m\left(t+s\right)-m\left(t\right))} $ | (20) |
A power law process is a special case of non-homogeneous Poisson process with intensity function given by [17]:
$ \lambda \left(t\right) = \left(\frac{\beta }{\gamma }\right)\left(\frac{t}{\gamma }\right)\genfrac{}{}{0pt}{}{\beta -1}{}, \gamma > 0, \beta > 1, t > 0 $ | (21) |
Meanwhile, the expectation value based on Eq 16 is given by:
$ m\left(t\right) = \left(\frac{t}{\gamma }\right)\genfrac{}{}{0pt}{}{\beta }{}, \gamma > 0, \beta > 1, t > 0 $ | (22) |
The intensity function of the power law process can be used to estimate the event rate at a given time. This is because the shape parameter $ \left(\beta \right) $ can describe how an intensity can increase or decrease. If β > 1, the intensity of an event will increase: if $ \beta < 1 $, the intensity of the event will decrease; and if $ \beta = 1 $, the power law process reduces to a homogeneous Poisson process [17].
There are several goodness-of-fit test procedures that can be used to test the suitability of the power law process model, including Kuiper's V Test, Watson's $ {U}^{2} $ Test, Anderson-Darling $ {A}^{2} $ Test, Shapiro-Wilk Test, and Cramér-von Mises Test. The Cramér-von Mises test uses the following hypothesis:
$ {H}_{0}: $ Event intensity fits the power law process model.
$ {H}_{1}: $ Event intensity does not fit the power law process model.
The Cramér-von Mises test statistic is expressed based on the following equation:
$ {C}_{R}^{2} = \frac{1}{12\left(n\right)}+\sum \limits_{i = 1}^{n}{\left(\overline {R}-\frac{2i-1}{2n}\right)}^{2} $ | (23) |
$ \overline {R} $ is the ratio power transformation given by Eq 29.
$ \overline {R} = {\left(\frac{{t}_{i}}{t}\right)}^{\overline {\beta }} $ | (24) |
$ \overline {\beta } $ is the unbiased estimator given by Eq 27.
$ \overline {\beta } = \frac{\left(n-2\right)}{\sum _{i = 1}^{n}\mathrm{ln}\left(\frac{t}{{t}_{i}}\right)} $ | (25) |
The $ {H}_{0} $ decision is accepted if the calculated value of the $ {C}_{R}^{2} $ test statistic is smaller than the critical value for the Cramér-von Mises test, which means that the power law process model is appropriate. If the value of the $ {C}_{R}^{2} $ test statistic is greater than the critical value for the Cramér-von Mises test, then $ {H}_{0} $ is rejected, which means that the model is not suitable, and a more suitable model needs to be used [17].
Suppose $ {t}_{1}, {t}_{2}, {t}_{3}, \dots , {t}_{n} $ are mutually independent random samples from a distribution with a joint probability density function $ f({t}_{1}, {t}_{2}, {t}_{3}, \dots , {t}_{n};\beta , \gamma ) $ with n representing the number of events occurring until time $ {t}_{i} $ for $ 0 < {t}_{1} < {t}_{2} < {t}_{3} < \dots < {t}_{n} $. If the joint likelihood function is expressed as a function of $ \beta $ and $ \gamma $, then the likelihood function is denoted as $ L({t}_{1}, {t}_{2}, {t}_{3}, \dots , {t}_{n};\beta , \gamma ). $ The likelihood function for the parameters $ \beta $ and $ \gamma $ is given as follows [17]:
$ L\left({t}_{1}, {t}_{2}, {t}_{3}, \dots , {t}_{n};\beta , \gamma \right) = \left(\prod \limits_{i = 1}^{n}\lambda \left({t}_{i};\beta , \gamma \right)\right)exp\left(-{\int }_{0}^{{t}_{n}}\lambda \left({t}_{i};\beta , \gamma \right)dt\right) $ | (26) |
Based on Eq 21, the likelihood function with intensity function $ \lambda \left({t}_{i};\beta , \gamma \right) $ in Eq 23 is:
$ L\left({t}_{i};\beta , \gamma \right) = \left(\prod \limits_{i = 1}^{n}\left(\frac{\beta }{\gamma }\right){\left(\frac{{t}_{i}}{\gamma }\right)}^{\genfrac{}{}{0pt}{}{\beta -1}{}}\right)exp\left(-{\int }_{0}^{{t}_{n}}\left(\frac{\beta }{\gamma }\right){\left(\frac{t}{\gamma }\right)}^{\genfrac{}{}{0pt}{}{\beta -1}{}}dt\right) $ | (27) |
Based on Eq 24, the logarithmic likelihood function $ \mathcal{L}\left(t;\beta , \gamma \right) = \mathrm{l}\mathrm{n}\left(L\left(t;\beta , \gamma \right)\right) $ is:
$ \mathcal{L}\left(t;\beta , \gamma \right) = n\mathrm{ln}\left(\beta \right)-n\beta \mathrm{ln}\left(\gamma \right)+(\beta -1)\sum \limits_{i = 1}^{n}\mathrm{ln}\left({t}_{i}\right)-{\left(\frac{t}{\gamma }\right)}^{\beta } $ | (28) |
Furthermore, Eq 25 is derived with respect to $ \beta $ and $ \gamma $ so that the maximum likelihood estimator is obtained as follows:
$ \widehat{\beta } = \frac{n}{\sum _{i = 1}^{n}\mathrm{ln}\left(\frac{t}{{t}_{i}}\right)} $ | (29) |
$ \widehat{\gamma } = \frac{t}{{n}^{\frac{1}{\widehat{\beta }}}} $ | (30) |
The SPEI is a drought index used to analyze meteorological drought conditions by considering the standardization of rainfall deficits with potential evapotranspiration (PET) or climate water balance. In this study, the SPEI calculation time scale used a 1-month period, which is adjusted to the needs of researchers to evaluate drought in the short term. SPEI defines a drought event as occurring when the SPEI value is below or equal to the -1 threshold, and the drought event ends when the SPEI value returns to positive. The index classifies drought levels into three main categories: dry, very dry, and extremely dry. The classification of drought levels by SPEI is based on the SPI classification table as follows [14,15,36]:
1) Moderately dry: Occurs when the SPEI value is between –1 and –1.49. This indicates mild drought that may affect water availability.
2) Severely dry: Occurs when the SPEI value is between –1.5 and –1.99. This indicates a more serious drought that can significantly impact agriculture, clean water, and ecosystems.
3) Extremely dry: Occurs when the SPEI value is below –2. This category represents the worst drought index and can cause major losses to agriculture, water availability, and the environment.
The SPEI calculation process uses Eqs 1–15. A time series plot of the calculated values of the 1-month SPEI for each observation area on Timor Island is shown in Figure 2.
Figure 2 shows that SPEI values close to 0 in each observation area indicate near-normal conditions, while positive or negative values indicate above or below-normal conditions. There are many negative SPEI values lower than or equal to –1, indicating that there are frequent droughts at dry, very dry, and extremely dry levels in each observation area of Timor Island. This result aligns with the research of Kuswanto et al. [38], which shows that very dry events are more common in the eastern region of NTT; in this case, Timor Island is the eastern region of NTT.
Furthermore, SPEI values lower than or equal to –1 were characterized to obtain drought intensity, duration, and severity during the observation period [14]. The characterization of 1-month SPEI values in each region of Timor Island can be seen in Table 2.
Observation area | Extremely dry index | Longest drought duration | Intensity of drought months by drought level | ||||
Value | Month of incident | Moderately dry | Severely dry | Extremely dry | Total | ||
Kupang City | –2.47 | August 1988 | 4 months | 63 months | 25 months | 6 months | 94 months |
Kupang | –2.74 | April 2016 | 4 months | 62 months | 26 months | 5 months | 93 months |
South Central Timor | –3.86 | August 1998 | 4 months | 59 months | 25 months | 6 months | 90 months |
North Central Timor | –2.61 | August 1998 | 4 months | 62 months | 20 months | 6 months | 88 months |
Malaka | –2.82 | August 2010 | 5 months | 66 months | 25 months | 4 months | 95 months |
Belu | –3.07 | August 2010 | 4 months | 59 months | 23 months | 5 months | 87 months |
The worst SPEI indices, with values lower than or equal to –2 (extremely dry conditions), have occurred throughout Timor Island. On average, the most severe SPEI drought indices occurred in April, except in Kupang Regency. The most extreme short-term meteorological drought ever recorded on Timor Island occurred in South Central Timor Regency, with a drought index of –3.86 in August 1998.
The longest droughts were as follows: In Kupang City for 4 consecutive months from August to November 1988; in Kupang Regency for 4 consecutive months, occurring from June to September 1998; in South Central Timor Regency for 4 consecutive months from May to September 1998 and again from February to March 2018; in North Central Timor Regency for 4 consecutive months, also in two different periods, namely from June to September 1996 and June to September 1998; in Malaka Regency for 5 consecutive months, occurring from August to November 2020; and in Belu Regency for 4 consecutive months, occurring from November 1997 to February 1998 and again from June to September 1998. In addition, the intensity of drought events in each region varies greatly. The distribution of the intensity of drought events can be seen in Figure 3.
Figure 3 illustrates the pattern of variance in the intensity of drought events, which is that extremely dry events tend to occur less frequently than severely dry and moderately dry ones. This phenomenon demonstrates the complexity of drought as a natural phenomenon that involves factors such as time distribution, scale, and varying intensity of occurrence. An appropriate and effective method is needed to understand and describe drought dynamics.
In the analysis of short-term meteorological drought using the 1-month SPEI on Timor Island, the natural trend where more severe events tend to be less frequent than weaker events can be interpreted with a stochastic process model, the power law process. This process can explain how the event's intensity is inversely proportional to its magnitude [16,18]. The power law process helps researchers understand the pattern of drought intensity and can be used to predict the likelihood of future drought events. Therefore, to better understand the spatial distribution of SPEI values on Timor Island, the SPEI index was mapped every month in 2023, as shown in Figure 4.
The 2023 drought distribution map shows that all areas of Timor Island were affected by drought, as seen from the orange and red colors. In 2023, there were 3 months of drought in Kupang City, Kupang, South Central Timor, and Malaka, and 2 months in North Central Timor and Belu. In January, all observed areas experienced drought at a moderately dry level; in October, almost all observed areas experienced drought at a severely dry level, except for Malaka Regency, which experienced a moderately dry level. In November, a moderately dry drought occurred in Malaka Regency, while other regions had returned to normal. In December, an extremely dry drought occurred in South Central Timor Regency, and a severely dry drought was observed in Kupang City and Kupang Regency.
The presentation of the drought distribution map only for 2023 is based on the need to provide up-to-date information on drought conditions on Timor Island. Although the range of observations covers the years 1981–2023, 2023 was chosen due to the relevance of the current information desired in this study. By focusing on that year, a more in-depth understanding of the current spatial distribution of SPEI values in the Timor Island region can be obtained.
We tested the suitability of the power law process model using the Cramér-von Mises test [17]. The hypothesis used is as follows:
$ {H}_{0}: $ The intensity of drought events based on the 1-month SPEI fits the power law process model.
$ {H}_{1}: $ The intensity of drought events based on the 1-month SPEI does not fit the power law process model.
The Cramér-von Mises statistic $ \left({C}_{R}^{2}\right) $ is obtained using Eqs 23–25. The results are presented in Table 3.
Observation area | $ {\boldsymbol{C}}_{\boldsymbol{R}}^{\bf{2}} $ | Critical value | Decision |
Kupang City | 0.082 | 0.22 | $ {H}_{0} $ accepted |
Kupang | 0.062 | 0.22 | $ {H}_{0} $ accepted |
South Central Timor | 0.089 | 0.22 | $ {H}_{0} $ accepted |
North Central Timor | 0.171 | 0.22 | $ {H}_{0} $ accepted |
Malaka | 0.043 | 0.22 | $ {H}_{0} $ accepted |
Belu | 0.022 | 0.22 | $ {H}_{0} $ accepted |
Based on Table 3, the $ {C}_{R}^{2} $ value of all observation areas is below the critical value determined based on the frequency of drought events, and $ {H}_{0} $ is accepted. These results indicate that the intensity of drought events based on 1-month SPEI in each region of Timor Island fits the power law process model.
The estimated values of the shape $ \left(\beta \right) $ and scale $ \left(\gamma \right) $ parameters in the power law process intensity function were obtained using the MLE method based on Eqs 29 and 30 using the time truncated estimate of the power law process [17]. Data used are drought frequency, time of occurrence, and observation time span. The parameter estimation results are presented in Table 4.
Observation area | Parameter | |
$ \widehat{\beta } $ | $ \widehat{\gamma } $ | |
Kupang City | 1.063 | 7.170 |
Kupang | 1.174 | 10.859 |
South Central Timor | 1.095 | 8.481 |
North Central Timor | 1.049 | 7.245 |
Malaka | 1.034 | 6.296 |
Belu | 1.112 | 9.282 |
Table 4 shows that the $ \widehat{\beta } $ parameter in each observation area was higher than 1. Based on the characteristics of the power law process, if the intensity function is greater than 1, the intensity of an event will increase [13]. Therefore, it can be concluded that the intensity of drought events in each region of Timor Island will increase, so it is necessary to estimate the frequency of future drought events.
Suppose the time value of drought occurrence in each observation area is $ {t}_{i} < {t}_{2} < \cdots < {t}_{n} < t $ so that $ N\left(t\right) = \{N, {t}_{i} < {t}_{2} < \cdots < {t}_{n};(0, t\left]\right\} $ with N expressing the frequency of months with drought occurrence in the time interval (0, t]. It is known that the total months observed from 1981–2023 are t = 516 months. The period of months to be estimated is the following 12 months, so if t + s, then 516 + 12 = 528, meaning that the last month of estimation is the 528th month. The initial time of estimation starts from the time after the last time of observation, so the estimation time interval becomes [517,528]. The number of months to be estimated is 12, so the possible value of n is n = 1, 2, 3, ..., 12. The interpretation of the observation time based on the range of 1981–2023 is as follows: The first month is January 2023, and the 516th month is December 2023. For the estimation months, the 517th month is January 2024, and the 528th month is December 2024.
The expected frequency m(t) of drought occurrence in the observation time [1,516] can be obtained by substituting the estimated values of parameters $ \widehat{\beta } $ and $ \widehat{\gamma } $ in each region of Timor Island into Eq 22. A comparison between expected values and real monthly frequency values is presented in Figure 5.
Based on Figure 4, the expected value shows a very good level of agreement with the actual value. This indicates that the power law process model accurately predicts the number of months with drought events. Furthermore, the estimation of the expected value $ \left(m\right(t+s\left)\right) $ of drought frequency in the future 12 months or in observations [517,528] is presented in Table 5.
Months | Kupang City | Kupang | South Central Timor | Nort Central Timor | Malaka | Belu |
517 (Jan 2024) | 94.207 | 93.244 | 90.195 | 87.957 | 95.393 | 87.374 |
518 (Feb 2024) | 94.400 | 93.456 | 90.386 | 88.136 | 95.584 | 87.562 |
519 (Mar 2024) | 94.594 | 93.668 | 90.577 | 88.314 | 95.775 | 87.750 |
520 (Apr 2024) | 94.788 | 93.879 | 90.768 | 88.493 | 95.965 | 87.938 |
521 (May 2024) | 94.981 | 94.091 | 90.959 | 88.671 | 96.156 | 88.126 |
522 (Jun 2024) | 95.175 | 94.304 | 91.150 | 88.850 | 96.347 | 88.315 |
523 (Jul 2024) | 95.369 | 94.516 | 91.342 | 89.028 | 96.538 | 88.503 |
524 (Aug 2024) | 95.563 | 94.728 | 91.533 | 89.207 | 96.729 | 88.691 |
525 (Sep 2024 | 95.756 | 94.940 | 91.724 | 89.385 | 96.920 | 88.879 |
526 (Oct 2024) | 95.950 | 95.152 | 91.916 | 89.564 | 97.111 | 89.067 |
527 (Nov 2024) | 96.144 | 95.365 | 92.107 | 89.743 | 97.301 | 89.256 |
528 (Dec 2024) | 96.338 | 95.577 | 92.299 | 89.921 | 97.492 | 89.444 |
The predicted results show that, in the next 12 months, the drought frequency in each region of Timor Island will increase by 2 months from the initial observation. For example, in Kupang City, the drought frequency in the 516-month observation is 94 months, and in the 528-month observation is 96 months. The same is true for every other observation area, where in the next 12 months, the frequency of drought events will increase by 2 months. Therefore, it is imperative to make early preparations and implement effective mitigation strategies to reduce the possible impacts of more frequent droughts in the future.
Using the homogeneous Poisson process, the probability value of the expected future frequency of drought can be obtained by substituting the expected value m(t+s) into the non-homogeneous Poisson process based on Eq 20, where m(t) is the expected value at the last time of observation. As mentioned, the power law process is a special case of the non-expected frequency of drought occurrence in the future, as presented in Tables 6–11. For example, if the expected frequency of drought in the next 12 months is 2 months, the probability values based on Tables 6–11 for each region of Timor Island are 0.264 for Kupang City, 0.254 for Kupang, 0.265 for South Central Timor, 0.269 for North Central Timor, 0.266 for Malaka, and 0.267 for Belu. In addition, the probability value of the expected frequency is 1, 3, and 12 months. For example, in Table 7, it is shown that the probability of drought frequency being 3 months in the future 12 months in Kupang City is 0.205. Meanwhile, in Kupang Regency, that probability value is 0.216. This interpretation also applies to each of the following estimates.
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.160 | |||||||||||
518 (Feb 2024) | 0.263 | 0.051 | ||||||||||
519 (Mar 2024) | 0.325 | 0.094 | 0.018 | |||||||||
520 (Apr 2024) | 0.357 | 0.138 | 0.036 | 0.007 | ||||||||
521 (May 2024) | 0.368 | 0.178 | 0.057 | 0.014 | 0.003 | |||||||
522 (Jun 2024) | 0.364 | 0.211 | 0.082 | 0.024 | 0.006 | 0.001 | ||||||
523 (Jul 2024) | 0.349 | 0.237 | 0.107 | 0.036 | 0.010 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.329 | 0.255 | 0.132 | 0.051 | 0.016 | 0.004 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.305 | 0.266 | 0.154 | 0.067 | 0.023 | 0.007 | 0.002 | 0.000 | 0.000 | |||
526 (Oct 2024) | 0.279 | 0.270 | 0.175 | 0.085 | 0.033 | 0.011 | 0.003 | 0.001 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.253 | 0.270 | 0.191 | 0.102 | 0.043 | 0.015 | 0.005 | 0.001 | 0.000 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.227 | 0.264 | 0.205 | 0.119 | 0.055 | 0.021 | 0.007 | 0.002 | 0.001 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.172 | |||||||||||
518 (Feb 2024) | 0.277 | 0.059 | ||||||||||
519 (Mar 2024) | 0.337 | 0.107 | 0.023 | |||||||||
520 (Apr 2024) | 0.363 | 0.154 | 0.043 | 0.009 | ||||||||
521 (May 2024) | 0.367 | 0.194 | 0.069 | 0.018 | 0.004 | |||||||
522 (Jun 2024) | 0.357 | 0.227 | 0.096 | 0.031 | 0.008 | 0.002 | ||||||
523 (Jul 2024) | 0.336 | 0.250 | 0.123 | 0.046 | 0.014 | 0.003 | 0.001 | |||||
524 (Aug 2024) | 0.311 | 0.264 | 0.149 | 0.063 | 0.021 | 0.006 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.283 | 0.270 | 0.172 | 0.082 | 0.031 | 0.010 | 0.003 | 0.001 | 0.000 | |||
526 (Oct 2024) | 0.254 | 0.270 | 0.191 | 0.101 | 0.043 | 0.015 | 0.005 | 0.001 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.226 | 0.264 | 0.205 | 0.120 | 0.056 | 0.022 | 0.007 | 0.002 | 0.001 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.200 | 0.254 | 0.216 | 0.137 | 0.070 | 0.030 | 0.011 | 0.003 | 0.001 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.158 | |||||||||||
518 (Feb 2024) | 0.261 | 0.050 | ||||||||||
519 (Mar 2024) | 0.323 | 0.093 | 0.018 | |||||||||
520 (Apr 2024) | 0.356 | 0.136 | 0.035 | 0.007 | ||||||||
521 (May 2024) | 0.367 | 0.175 | 0.056 | 0.013 | 0.003 | |||||||
522 (Jun 2024) | 0.364 | 0.209 | 0.080 | 0.023 | 0.005 | 0.001 | ||||||
523 (Jul 2024) | 0.351 | 0.235 | 0.105 | 0.035 | 0.009 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.331 | 0.253 | 0.129 | 0.049 | 0.015 | 0.004 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.308 | 0.265 | 0.152 | 0.065 | 0.022 | 0.006 | 0.002 | 0.000 | 0.000 | |||
526 (Oct 2024) | 0.283 | 0.270 | 0.172 | 0.082 | 0.031 | 0.010 | 0.003 | 0.001 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.257 | 0.270 | 0.189 | 0.100 | 0.042 | 0.015 | 0.004 | 0.001 | 0.000 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.231 | 0.265 | 0.203 | 0.116 | 0.053 | 0.020 | 0.007 | 0.002 | 0.000 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.149 | |||||||||||
518 (Feb 2024) | 0.250 | 0.045 | ||||||||||
519 (Mar 2024) | 0.313 | 0.084 | 0.015 | |||||||||
520 (Apr 2024) | 0.350 | 0.125 | 0.030 | 0.005 | ||||||||
521 (May 2024) | 0.366 | 0.163 | 0.049 | 0.011 | 0.002 | |||||||
522 (Jun 2024) | 0.367 | 0.197 | 0.070 | 0.019 | 0.004 | 0.001 | ||||||
523 (Jul 2024) | 0.358 | 0.224 | 0.093 | 0.029 | 0.007 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.342 | 0.244 | 0.116 | 0.042 | 0.012 | 0.003 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.322 | 0.259 | 0.139 | 0.056 | 0.018 | 0.005 | 0.001 | 0.000 | 0.000 | |||
526 (Oct 2024) | 0.299 | 0.267 | 0.159 | 0.071 | 0.025 | 0.008 | 0.002 | 0.000 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.276 | 0.271 | 0.177 | 0.087 | 0.034 | 0.011 | 0.003 | 0.001 | 0.000 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.251 | 0.269 | 0.192 | 0.103 | 0.044 | 0.016 | 0.005 | 0.001 | 0.000 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.158 | |||||||||||
518 (Feb 2024) | 0.261 | 0.050 | ||||||||||
519 (Mar 2024) | 0.323 | 0.092 | 0.018 | |||||||||
520 (Apr 2024) | 0.356 | 0.136 | 0.035 | 0.007 | ||||||||
521 (May 2024) | 0.367 | 0.175 | 0.056 | 0.013 | 0.003 | |||||||
522 (Jun 2024) | 0.364 | 0.209 | 0.080 | 0.023 | 0.005 | 0.001 | ||||||
523 (Jul 2024) | 0.351 | 0.235 | 0.104 | 0.035 | 0.009 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.332 | 0.253 | 0.129 | 0.049 | 0.015 | 0.004 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.308 | 0.265 | 0.152 | 0.065 | 0.022 | 0.006 | 0.002 | 0.000 | 0.000 | |||
526 (Oct 2024) | 0.283 | 0.270 | 0.172 | 0.082 | 0.031 | 0.010 | 0.003 | 0.001 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.257 | 0.270 | 0.189 | 0.099 | 0.042 | 0.015 | 0.004 | 0.001 | 0.000 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.232 | 0.266 | 0.203 | 0.116 | 0.053 | 0.020 | 0.007 | 0.002 | 0.000 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.156 | |||||||||||
518 (Feb 2024) | 0.258 | 0.049 | ||||||||||
519 (Mar 2024) | 0.321 | 0.090 | 0.017 | |||||||||
520 (Apr 2024) | 0.354 | 0.133 | 0.033 | 0.006 | ||||||||
521 (May 2024) | 0.367 | 0.173 | 0.054 | 0.013 | 0.002 | |||||||
522 (Jun 2024) | 0.365 | 0.206 | 0.077 | 0.022 | 0.005 | 0.001 | ||||||
523 (Jul 2024) | 0.353 | 0.232 | 0.102 | 0.034 | 0.009 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.334 | 0.251 | 0.126 | 0.047 | 0.014 | 0.004 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.311 | 0.264 | 0.149 | 0.063 | 0.021 | 0.006 | 0.001 | 0.000 | ||||
526 (Oct 2024) | 0.287 | 0.270 | 0.169 | 0.080 | 0.030 | 0.009 | 0.003 | 0.001 | 0.000 | |||
527 (Nov 2024) | 0.261 | 0.270 | 0.186 | 0.096 | 0.040 | 0.014 | 0.004 | 0.001 | 0.000 | 0.000 | ||
528 (Dec 2024) | 0.236 | 0.267 | 0.201 | 0.113 | 0.051 | 0.019 | 0.006 | 0.002 | 0.000 | 0.000 | 0.000 | 0.000 |
Note that a limitation of this study lies in the estimation of the power law process model, where the researchers stipulated that the value of the shape parameter (β) should be greater than 1 to indicate an increase in event intensity, as described in Rigdon and Basu’s study [17]. In addition, it should be noted that other modeling approaches or the use of different drought index datasets can be considered. For example, Ghasemi et al. used a Gaussian process regression model to forecast the SPEI drought index [22]; other drought indices, such as the Palmer drought severity index (PDSI) or the Z-Index, may also be used. However, it is important to highlight that future research could focus on modeling and comparing different results. Furthermore, the model used to analyze the SPEI drought index in Timor Island has a wide potential application in cases where extreme events are rare compared to common weaker events, such as in the study of earthquakes, extreme weather, temperature changes, and others.
The analysis of short-term meteorological drought events using SPEI for a 1-month period on Timor Island shows that extremely dry events are less frequent than very dry and dry events. The results of the power law process parameter estimation show a $ \beta > 1 $ value in all regions of Timor Island, specifically 1.063 for Kupang City, 1.174 for Kupang Regency, 1.095 for South Central Timor Regency, 1.049 for North Central Timor Regency, 1.034 for Malaka Regency, and 1.112 for Belu Regency. This indicates an increase in drought events in the future. In the next 12 months, the estimated duration of short-term meteorological droughts in all regions is 2 months, with the following probability values: 0.264 for Kupang City, 0.25 for Kupang, 0.265 for South Central Timor, 0.269 for North Central Timor, 0.265 for Malaka, and 0.266 for Belu.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors are grateful for financial support from the Directorate General of Higher Education, Ministry of Education and Culture, Research and Technology of Indonesia.
The authors stated that there is no conflict of interest for the study.
[1] |
X. X. Ke, Y. Pang, K. Chen, D. Zhang, F. Wang, S. Zhu, et al., Knockdown of arsenic resistance protein 2 inhibits human glioblastoma cell proliferation through the MAPK/ERK pathway, Oncol. Rep., 40 (2018), 3313-3322. https://doi.org/10.3892/or.2018.6777 doi: 10.3892/or.2018.6777
![]() |
[2] |
H. Y. Li, B. B. Lv, Y. H. Bi, FABP4 accelerates glioblastoma cell growth and metastasis through Wnt10b signalling, Eur. Rev. Med. Pharmacol. Sci., 22 (2018), 7807-7818. https://doi.org/10.26355/eurrev_201811_16405 doi: 10.26355/eurrev_201811_16405
![]() |
[3] |
T. Mashimo, K. Pichumani, V. Vemireddy, K. J. Hatanpaa, D. K. Singh, S. Sirasanagandla, et al., Acetate is a bioenergetic substrate for human glioblastoma and brain metastases, Cell, 159 (2014), 1603-1614. https://doi.org/10.1016/j.cell.2014.11.025 doi: 10.1016/j.cell.2014.11.025
![]() |
[4] |
A. Vartanian, S. Agnihotri, M. R. Wilson, K. E. Burrell, P. D. Tonge, A. Alamsahebpour, et al., Targeting hexokinase 2 enhances response to radio-chemotherapy in glioblastoma, Oncotarget, 7 (2016), 69518-69535. https://doi.org/10.18632/oncotarget.11680 doi: 10.18632/oncotarget.11680
![]() |
[5] |
B. Huang, T. A. Dolecek, Q. Chen, C. R. Garcia, T. Pittman, J. L. Villano, Characteristics and survival outcomes associated with the lack of radiation in the treatment of glioblastoma, Med. Oncol., 35 (2018), 74. https://doi.org/10.1007/s12032-018-1134-3 doi: 10.1007/s12032-018-1134-3
![]() |
[6] |
Z. Shboul, L. Vidyaratne, M. Alam, S. M. S. Reza, K. M. Iftekharuddin, Glioblastoma and survival prediction, Lect. Notes Comput. Sci., 10670 (2018), 358-368. https://doi.org/10.1007/978-3-319-75238-9_31 doi: 10.1007/978-3-319-75238-9_31
![]() |
[7] |
J. K. Sa, S. H. Kim, J. K. Lee, H. J. Cho, Y. J. Shin, H. Shin, et al., Identification of genomic and molecular traits that present therapeutic vulnerability to HGF-targeted therapy in glioblastoma, Neuro-oncology, 21 (2019), 222-233. https://doi.org/10.1093/neuonc/noy105 doi: 10.1093/neuonc/noy105
![]() |
[8] |
Z. C. Zhu, J. W. Liu, K. Li, J. Zheng, Z. Q. Xiong, KPNB1 inhibition disrupts proteostasis and triggers unfolded protein response-mediated apoptosis in glioblastoma cells, Oncogene, 37 (2018), 2936-2952. https://doi.org/10.1038/s41388-018-0180-9 doi: 10.1038/s41388-018-0180-9
![]() |
[9] |
M. Westphal, C. L. Maire, K. Lamszus, EGFR as a target for glioblastoma treatment: An unfulfilled promise, CNS Drugs, 31 (2017), 723-735. https://doi.org/10.1007/s40263-017-0456-6 doi: 10.1007/s40263-017-0456-6
![]() |
[10] |
N. Hirokawa, Y, Tanaka, Kinesin superfamily proteins (KIFs): Various functions and their relevance for important phenomena in life and diseases, Exp. Cell Res., 334 (2015), 16-25. https://doi.org/10.1016/j.yexcr.2015.02.016 doi: 10.1016/j.yexcr.2015.02.016
![]() |
[11] | N. Hirokawa, From electron microscopy to molecular cell biology, molecular genetics and structural biology: Intracellular transport and kinesin superfamily proteins, KIFs: Genes, structure, dynamics and functions, J. Electron Microsc., 60 (2011), 63-92. https://doi.org/10.1093/jmicro/dfr051 |
[12] |
S. S. Siddiqui, Metazoan motor models: Kinesin superfamily in C. elegans, Traffic, 3 (2002), 20-28. https://doi.org/10.1034/j.1600-0854.2002.30104.x doi: 10.1034/j.1600-0854.2002.30104.x
![]() |
[13] |
H. Miki, M. Setou, K. Kaneshiro, N. Hirokawa, All kinesin superfamily protein, KIF, genes in mouse and human, Proc. Natl. Acad. Sci. U.S.A., 98 (2001), 7004-7011. https://doi.org/10.1073/pnas.111145398 doi: 10.1073/pnas.111145398
![]() |
[14] |
Y. M. Lee, W. Kim, Kinesin superfamily protein member 4 (KIF4) is localized to midzone and midbody in dividing cells, Exp. Mol. Med., 36 (2004), 93-97. https://doi.org/10.1038/emm.2004.13 doi: 10.1038/emm.2004.13
![]() |
[15] |
Z. Shen, A. R. Collatos, J. P. Bibeau, F. Furt, L. Vidali, Phylogenetic analysis of the Kinesin superfamily from physcomitrella, Front. Plant Sci., 3 (2012), 230. https://doi.org/10.3389/fpls.2012.00230 doi: 10.3389/fpls.2012.00230
![]() |
[16] |
K. Tang, N. H. Toda, A microtubule polymerase cooperates with the kinesin-6 motor and a microtubule cross-linker to promote bipolar spindle assembly in the absence of kinesin-5 and kinesin-14 in fission yeast, Mol. Biol. Cell, 28 (2017), 3647-3659. https://doi.org/10.1091/mbc.e17-08-0497 doi: 10.1091/mbc.e17-08-0497
![]() |
[17] |
T. McHugh, H. Drechsler, A. D. McAinsh, N. J. Carter, R. A. Cross, Kif15 functions as an active mechanical ratchet, Mol. Biol. Cell, 29 (2018), 1743-1752. https://doi.org/10.1091/mbc.E18-03-0151 doi: 10.1091/mbc.E18-03-0151
![]() |
[18] |
M. Xu, D. Liu, Z. Dong, X. Wang, X. Wang, Y. Liu, et al., Kinesin-12 influences axonal growth during zebrafish neural development, Cytoskeleton, 71 (2014), 555-563. https://doi.org/10.1002/cm.21193 doi: 10.1002/cm.21193
![]() |
[19] |
J. Feng, Z. Hu, H. Chen, J. Hua, R. Wu, Z. Dong, et al., Depletion of kinesin-12, a myosin-ⅡB-interacting protein, promotes migration of cortical astrocytes, J. Cell. Sci., 129 (2016), 2438-2447. https://doi.org/10.1242/jcs.181867 doi: 10.1242/jcs.181867
![]() |
[20] |
H. Miki, M. Setou, K. Kaneshiro, N. Hirokawa, All kinesin superfamily protein, KIF, genes in mouse and human, Proc. Natl. Acad. Sci. U.S.A., 98 (2001), 7004-7011. https://doi.org/10.1073/pnas.111145398 doi: 10.1073/pnas.111145398
![]() |
[21] |
J. Wang, X. Guo, C. Xie, J. Jiang, KIF15 promotes pancreatic cancer proliferation via the MEK-ERK signalling pathway, Br. J. Cancer, 117 (2017), 245-255. https://doi.org/10.1038/bjc.2017.165 doi: 10.1038/bjc.2017.165
![]() |
[22] |
Y. Qiao, J. Chen, C. Ma, Y. Liu, P. Li, Y. Wang, et al., Increased KIF15 expression predicts a poor prognosis in patients with lung adenocarcinoma, Cell Physiol. Biochem., 51 (2018), 1-10. https://doi.org/10.1159/000495155 doi: 10.1159/000495155
![]() |
[23] |
G. Harris, D. Jayamanne, H. Wheeler, C. Gzell, M. Kastelan, G. Schembri, et al., Survival Outcomes of Elderly Patients With glioblastoma multiforme in their 75th year or older treated with adjuvant therapy, Int. J. Radiat. Oncol. Biol. Phys., 98 (2017), 802-810. https://doi.org/10.1016/j.ijrobp.2017.02.028 doi: 10.1016/j.ijrobp.2017.02.028
![]() |
[24] |
K. K. Jain, A critical overview of targeted therapies for glioblastoma, Front. Oncol., 8 (2018), 419. https://doi.org/10.3389/fonc.2018.00419 doi: 10.3389/fonc.2018.00419
![]() |
[25] |
Y. Zheng, N. Gao, Y. L. Fu, B. Y. Zhang, X. L. Li, P. Gupta, et al., Generation of regulable EGFRvⅢ targeted chimeric antigen receptor T cells for adoptive cell therapy of glioblastoma, Biochem. Biophys. Res. Commun., 507 (2018), 59-66. https://doi.org/10.1016/j.bbrc.2018.10.151 doi: 10.1016/j.bbrc.2018.10.151
![]() |
[26] |
M. Momeny, F. Moghaddaskho, N. K. Gortany, H. Yousefi, Z. Sabourinejad, G. Zarrinrad, et al., Blockade of vascular endothelial growth factor receptors by tivozanib has potential anti-tumour effects on human glioblastoma cells, Sci. Rep., 7 (2017), 44075. https://doi.org/10.1038/srep44075 doi: 10.1038/srep44075
![]() |
[27] |
E. G. Sturgill, S. R. Norris, Y. Guo, R. Ohi, Kinesin-5 inhibitor resistance is driven by kinesin-12, J. Cell. Biol., 213 (2016), 213-227. https://doi.org/10.1083/jcb.201507036 doi: 10.1083/jcb.201507036
![]() |
[28] |
C. Müller, D. Gross, V. Sarli, M. Gartner, A. Giannis, G. Bernhardt, et al., Inhibitors of kinesin Eg5: Antiproliferative activity of monastrol analogues against human glioblastoma cells, Cancer Chemother. Pharmacol., 59 (2007), 157-164. https://doi.org/10.1007/s00280-006-0254-1 doi: 10.1007/s00280-006-0254-1
![]() |
[29] |
A. I. Marcus, U. Peters, S. L. Thomas, S. Garrett, A. Zelnak, T. M. Kapoor, et al., Mitotic kinesin inhibitors induce mitotic arrest and cell death in Taxol-resistant and sensitive cancer cells, J. Biol. Chem., 280 (2005), 11569-11577. https://doi.org/10.1074/jbc.M413471200 doi: 10.1074/jbc.M413471200
![]() |
[30] |
D. W. Buster, D. H. Baird, W. Yu, J. M. Solowska, M. Chauviere, A. Mazurek, Expression of the mitotic kinesin Kif15 in postmitotic neurons: Implications for neuronal migration and development, J. Neurocytol., 32 (2003), 79-96. https://doi.org/10.1023/A:1027332432740 doi: 10.1023/A:1027332432740
![]() |
[31] |
B. Stangeland, A. A. Mughal, Z. Grieg, C. J. Sandberg, M. Joel, S. Nygard, et al., Combined expressional analysis, bioinformatics and targeted proteomics identify new potential therapeutic targets in glioblastoma stem cells, Oncotarget, 6 (2015), 26192-26215. https://doi.org/10.18632/oncotarget.4613 doi: 10.18632/oncotarget.4613
![]() |
[32] |
O. Rath, F. Kozielski, Kinesins and cancer, Nat. Rev. Cancer, 12 (2012), 527-539. https://doi.org/10.1038/nrc3310 doi: 10.1038/nrc3310
![]() |
[33] |
G. Bergnes, K. Brejc, L. Belmont, Mitotic kinesins: Prospects for antimitotic drug discovery, Curr. Top. Med. Chem., 5 (2005), 127-145. https://doi.org/10.2174/1568026053507697 doi: 10.2174/1568026053507697
![]() |
[34] |
Z. Z. Wang, J. Yang, B. H. Jiang, J. B. Di, P. Gao, L. Peng, et al., KIF14 promotes cell proliferation via activation of Akt and is directly targeted by miR-200c in colorectal cancer, Int. J. Oncol., 53 (2018), 1939-1952. https://doi.org/10.3892/ijo.2018.4546 doi: 10.3892/ijo.2018.4546
![]() |
[35] |
X. Zhao, L. L. Zhou, X. Li, J. Ni, P. Chen, R. Ma, Overexpression of KIF20A confers malignant phenotype of lung adenocarcinoma by promoting cell proliferation and inhibiting apoptosis, Cancer Med., 7 (2018), 4678-4689. https://doi.org/10.1002/cam4.1710 doi: 10.1002/cam4.1710
![]() |
[36] |
Y. Teng, B. Guo, X. Mu, S. Liu, KIF26B promotes cell proliferation and migration through the FGF2/ERK signaling pathway in breast cancer, Biomed. Pharmacother., 108 (2018), 766-773. https://doi.org/10.1016/j.biopha.2018.09.036 doi: 10.1016/j.biopha.2018.09.036
![]() |
[37] |
Q. Y. Wang, B. Han, W. Huang, C. J. Qi, F. Liu, Identification of KIF15 as a potential therapeutic target and prognostic factor for glioma, Oncol. Rep., 43 (2020), 1035-1044. https://doi.org/10.3892/or.2020.7510 doi: 10.3892/or.2020.7510
![]() |
1. | Nurtiti Sunusi, Nur Hikmah Auliana, Assessing SPI and SPEI for Drought Forecasting through the Power Law Process: A Case Study in South Sulawesi, Indonesia, 2025, 22150161, 103235, 10.1016/j.mex.2025.103235 |
SPEI value | Classification |
$ \ge $2.00 | Extremely wet |
1.50–1.99 | Very wet |
1.00–1.49 | Moderately wet |
(–0.99)–0.99 | Normal |
(–1.00)–(–1.49) | Moderately dry |
(–1.50)–(–1.99) | Severely dry |
$ \le $ –2.00 | Extremely dry |
Observation area | Extremely dry index | Longest drought duration | Intensity of drought months by drought level | ||||
Value | Month of incident | Moderately dry | Severely dry | Extremely dry | Total | ||
Kupang City | –2.47 | August 1988 | 4 months | 63 months | 25 months | 6 months | 94 months |
Kupang | –2.74 | April 2016 | 4 months | 62 months | 26 months | 5 months | 93 months |
South Central Timor | –3.86 | August 1998 | 4 months | 59 months | 25 months | 6 months | 90 months |
North Central Timor | –2.61 | August 1998 | 4 months | 62 months | 20 months | 6 months | 88 months |
Malaka | –2.82 | August 2010 | 5 months | 66 months | 25 months | 4 months | 95 months |
Belu | –3.07 | August 2010 | 4 months | 59 months | 23 months | 5 months | 87 months |
Observation area | $ {\boldsymbol{C}}_{\boldsymbol{R}}^{\bf{2}} $ | Critical value | Decision |
Kupang City | 0.082 | 0.22 | $ {H}_{0} $ accepted |
Kupang | 0.062 | 0.22 | $ {H}_{0} $ accepted |
South Central Timor | 0.089 | 0.22 | $ {H}_{0} $ accepted |
North Central Timor | 0.171 | 0.22 | $ {H}_{0} $ accepted |
Malaka | 0.043 | 0.22 | $ {H}_{0} $ accepted |
Belu | 0.022 | 0.22 | $ {H}_{0} $ accepted |
Observation area | Parameter | |
$ \widehat{\beta } $ | $ \widehat{\gamma } $ | |
Kupang City | 1.063 | 7.170 |
Kupang | 1.174 | 10.859 |
South Central Timor | 1.095 | 8.481 |
North Central Timor | 1.049 | 7.245 |
Malaka | 1.034 | 6.296 |
Belu | 1.112 | 9.282 |
Months | Kupang City | Kupang | South Central Timor | Nort Central Timor | Malaka | Belu |
517 (Jan 2024) | 94.207 | 93.244 | 90.195 | 87.957 | 95.393 | 87.374 |
518 (Feb 2024) | 94.400 | 93.456 | 90.386 | 88.136 | 95.584 | 87.562 |
519 (Mar 2024) | 94.594 | 93.668 | 90.577 | 88.314 | 95.775 | 87.750 |
520 (Apr 2024) | 94.788 | 93.879 | 90.768 | 88.493 | 95.965 | 87.938 |
521 (May 2024) | 94.981 | 94.091 | 90.959 | 88.671 | 96.156 | 88.126 |
522 (Jun 2024) | 95.175 | 94.304 | 91.150 | 88.850 | 96.347 | 88.315 |
523 (Jul 2024) | 95.369 | 94.516 | 91.342 | 89.028 | 96.538 | 88.503 |
524 (Aug 2024) | 95.563 | 94.728 | 91.533 | 89.207 | 96.729 | 88.691 |
525 (Sep 2024 | 95.756 | 94.940 | 91.724 | 89.385 | 96.920 | 88.879 |
526 (Oct 2024) | 95.950 | 95.152 | 91.916 | 89.564 | 97.111 | 89.067 |
527 (Nov 2024) | 96.144 | 95.365 | 92.107 | 89.743 | 97.301 | 89.256 |
528 (Dec 2024) | 96.338 | 95.577 | 92.299 | 89.921 | 97.492 | 89.444 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.160 | |||||||||||
518 (Feb 2024) | 0.263 | 0.051 | ||||||||||
519 (Mar 2024) | 0.325 | 0.094 | 0.018 | |||||||||
520 (Apr 2024) | 0.357 | 0.138 | 0.036 | 0.007 | ||||||||
521 (May 2024) | 0.368 | 0.178 | 0.057 | 0.014 | 0.003 | |||||||
522 (Jun 2024) | 0.364 | 0.211 | 0.082 | 0.024 | 0.006 | 0.001 | ||||||
523 (Jul 2024) | 0.349 | 0.237 | 0.107 | 0.036 | 0.010 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.329 | 0.255 | 0.132 | 0.051 | 0.016 | 0.004 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.305 | 0.266 | 0.154 | 0.067 | 0.023 | 0.007 | 0.002 | 0.000 | 0.000 | |||
526 (Oct 2024) | 0.279 | 0.270 | 0.175 | 0.085 | 0.033 | 0.011 | 0.003 | 0.001 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.253 | 0.270 | 0.191 | 0.102 | 0.043 | 0.015 | 0.005 | 0.001 | 0.000 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.227 | 0.264 | 0.205 | 0.119 | 0.055 | 0.021 | 0.007 | 0.002 | 0.001 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.172 | |||||||||||
518 (Feb 2024) | 0.277 | 0.059 | ||||||||||
519 (Mar 2024) | 0.337 | 0.107 | 0.023 | |||||||||
520 (Apr 2024) | 0.363 | 0.154 | 0.043 | 0.009 | ||||||||
521 (May 2024) | 0.367 | 0.194 | 0.069 | 0.018 | 0.004 | |||||||
522 (Jun 2024) | 0.357 | 0.227 | 0.096 | 0.031 | 0.008 | 0.002 | ||||||
523 (Jul 2024) | 0.336 | 0.250 | 0.123 | 0.046 | 0.014 | 0.003 | 0.001 | |||||
524 (Aug 2024) | 0.311 | 0.264 | 0.149 | 0.063 | 0.021 | 0.006 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.283 | 0.270 | 0.172 | 0.082 | 0.031 | 0.010 | 0.003 | 0.001 | 0.000 | |||
526 (Oct 2024) | 0.254 | 0.270 | 0.191 | 0.101 | 0.043 | 0.015 | 0.005 | 0.001 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.226 | 0.264 | 0.205 | 0.120 | 0.056 | 0.022 | 0.007 | 0.002 | 0.001 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.200 | 0.254 | 0.216 | 0.137 | 0.070 | 0.030 | 0.011 | 0.003 | 0.001 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.158 | |||||||||||
518 (Feb 2024) | 0.261 | 0.050 | ||||||||||
519 (Mar 2024) | 0.323 | 0.093 | 0.018 | |||||||||
520 (Apr 2024) | 0.356 | 0.136 | 0.035 | 0.007 | ||||||||
521 (May 2024) | 0.367 | 0.175 | 0.056 | 0.013 | 0.003 | |||||||
522 (Jun 2024) | 0.364 | 0.209 | 0.080 | 0.023 | 0.005 | 0.001 | ||||||
523 (Jul 2024) | 0.351 | 0.235 | 0.105 | 0.035 | 0.009 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.331 | 0.253 | 0.129 | 0.049 | 0.015 | 0.004 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.308 | 0.265 | 0.152 | 0.065 | 0.022 | 0.006 | 0.002 | 0.000 | 0.000 | |||
526 (Oct 2024) | 0.283 | 0.270 | 0.172 | 0.082 | 0.031 | 0.010 | 0.003 | 0.001 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.257 | 0.270 | 0.189 | 0.100 | 0.042 | 0.015 | 0.004 | 0.001 | 0.000 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.231 | 0.265 | 0.203 | 0.116 | 0.053 | 0.020 | 0.007 | 0.002 | 0.000 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.149 | |||||||||||
518 (Feb 2024) | 0.250 | 0.045 | ||||||||||
519 (Mar 2024) | 0.313 | 0.084 | 0.015 | |||||||||
520 (Apr 2024) | 0.350 | 0.125 | 0.030 | 0.005 | ||||||||
521 (May 2024) | 0.366 | 0.163 | 0.049 | 0.011 | 0.002 | |||||||
522 (Jun 2024) | 0.367 | 0.197 | 0.070 | 0.019 | 0.004 | 0.001 | ||||||
523 (Jul 2024) | 0.358 | 0.224 | 0.093 | 0.029 | 0.007 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.342 | 0.244 | 0.116 | 0.042 | 0.012 | 0.003 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.322 | 0.259 | 0.139 | 0.056 | 0.018 | 0.005 | 0.001 | 0.000 | 0.000 | |||
526 (Oct 2024) | 0.299 | 0.267 | 0.159 | 0.071 | 0.025 | 0.008 | 0.002 | 0.000 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.276 | 0.271 | 0.177 | 0.087 | 0.034 | 0.011 | 0.003 | 0.001 | 0.000 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.251 | 0.269 | 0.192 | 0.103 | 0.044 | 0.016 | 0.005 | 0.001 | 0.000 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.158 | |||||||||||
518 (Feb 2024) | 0.261 | 0.050 | ||||||||||
519 (Mar 2024) | 0.323 | 0.092 | 0.018 | |||||||||
520 (Apr 2024) | 0.356 | 0.136 | 0.035 | 0.007 | ||||||||
521 (May 2024) | 0.367 | 0.175 | 0.056 | 0.013 | 0.003 | |||||||
522 (Jun 2024) | 0.364 | 0.209 | 0.080 | 0.023 | 0.005 | 0.001 | ||||||
523 (Jul 2024) | 0.351 | 0.235 | 0.104 | 0.035 | 0.009 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.332 | 0.253 | 0.129 | 0.049 | 0.015 | 0.004 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.308 | 0.265 | 0.152 | 0.065 | 0.022 | 0.006 | 0.002 | 0.000 | 0.000 | |||
526 (Oct 2024) | 0.283 | 0.270 | 0.172 | 0.082 | 0.031 | 0.010 | 0.003 | 0.001 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.257 | 0.270 | 0.189 | 0.099 | 0.042 | 0.015 | 0.004 | 0.001 | 0.000 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.232 | 0.266 | 0.203 | 0.116 | 0.053 | 0.020 | 0.007 | 0.002 | 0.000 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.156 | |||||||||||
518 (Feb 2024) | 0.258 | 0.049 | ||||||||||
519 (Mar 2024) | 0.321 | 0.090 | 0.017 | |||||||||
520 (Apr 2024) | 0.354 | 0.133 | 0.033 | 0.006 | ||||||||
521 (May 2024) | 0.367 | 0.173 | 0.054 | 0.013 | 0.002 | |||||||
522 (Jun 2024) | 0.365 | 0.206 | 0.077 | 0.022 | 0.005 | 0.001 | ||||||
523 (Jul 2024) | 0.353 | 0.232 | 0.102 | 0.034 | 0.009 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.334 | 0.251 | 0.126 | 0.047 | 0.014 | 0.004 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.311 | 0.264 | 0.149 | 0.063 | 0.021 | 0.006 | 0.001 | 0.000 | ||||
526 (Oct 2024) | 0.287 | 0.270 | 0.169 | 0.080 | 0.030 | 0.009 | 0.003 | 0.001 | 0.000 | |||
527 (Nov 2024) | 0.261 | 0.270 | 0.186 | 0.096 | 0.040 | 0.014 | 0.004 | 0.001 | 0.000 | 0.000 | ||
528 (Dec 2024) | 0.236 | 0.267 | 0.201 | 0.113 | 0.051 | 0.019 | 0.006 | 0.002 | 0.000 | 0.000 | 0.000 | 0.000 |
SPEI value | Classification |
$ \ge $2.00 | Extremely wet |
1.50–1.99 | Very wet |
1.00–1.49 | Moderately wet |
(–0.99)–0.99 | Normal |
(–1.00)–(–1.49) | Moderately dry |
(–1.50)–(–1.99) | Severely dry |
$ \le $ –2.00 | Extremely dry |
Observation area | Extremely dry index | Longest drought duration | Intensity of drought months by drought level | ||||
Value | Month of incident | Moderately dry | Severely dry | Extremely dry | Total | ||
Kupang City | –2.47 | August 1988 | 4 months | 63 months | 25 months | 6 months | 94 months |
Kupang | –2.74 | April 2016 | 4 months | 62 months | 26 months | 5 months | 93 months |
South Central Timor | –3.86 | August 1998 | 4 months | 59 months | 25 months | 6 months | 90 months |
North Central Timor | –2.61 | August 1998 | 4 months | 62 months | 20 months | 6 months | 88 months |
Malaka | –2.82 | August 2010 | 5 months | 66 months | 25 months | 4 months | 95 months |
Belu | –3.07 | August 2010 | 4 months | 59 months | 23 months | 5 months | 87 months |
Observation area | $ {\boldsymbol{C}}_{\boldsymbol{R}}^{\bf{2}} $ | Critical value | Decision |
Kupang City | 0.082 | 0.22 | $ {H}_{0} $ accepted |
Kupang | 0.062 | 0.22 | $ {H}_{0} $ accepted |
South Central Timor | 0.089 | 0.22 | $ {H}_{0} $ accepted |
North Central Timor | 0.171 | 0.22 | $ {H}_{0} $ accepted |
Malaka | 0.043 | 0.22 | $ {H}_{0} $ accepted |
Belu | 0.022 | 0.22 | $ {H}_{0} $ accepted |
Observation area | Parameter | |
$ \widehat{\beta } $ | $ \widehat{\gamma } $ | |
Kupang City | 1.063 | 7.170 |
Kupang | 1.174 | 10.859 |
South Central Timor | 1.095 | 8.481 |
North Central Timor | 1.049 | 7.245 |
Malaka | 1.034 | 6.296 |
Belu | 1.112 | 9.282 |
Months | Kupang City | Kupang | South Central Timor | Nort Central Timor | Malaka | Belu |
517 (Jan 2024) | 94.207 | 93.244 | 90.195 | 87.957 | 95.393 | 87.374 |
518 (Feb 2024) | 94.400 | 93.456 | 90.386 | 88.136 | 95.584 | 87.562 |
519 (Mar 2024) | 94.594 | 93.668 | 90.577 | 88.314 | 95.775 | 87.750 |
520 (Apr 2024) | 94.788 | 93.879 | 90.768 | 88.493 | 95.965 | 87.938 |
521 (May 2024) | 94.981 | 94.091 | 90.959 | 88.671 | 96.156 | 88.126 |
522 (Jun 2024) | 95.175 | 94.304 | 91.150 | 88.850 | 96.347 | 88.315 |
523 (Jul 2024) | 95.369 | 94.516 | 91.342 | 89.028 | 96.538 | 88.503 |
524 (Aug 2024) | 95.563 | 94.728 | 91.533 | 89.207 | 96.729 | 88.691 |
525 (Sep 2024 | 95.756 | 94.940 | 91.724 | 89.385 | 96.920 | 88.879 |
526 (Oct 2024) | 95.950 | 95.152 | 91.916 | 89.564 | 97.111 | 89.067 |
527 (Nov 2024) | 96.144 | 95.365 | 92.107 | 89.743 | 97.301 | 89.256 |
528 (Dec 2024) | 96.338 | 95.577 | 92.299 | 89.921 | 97.492 | 89.444 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.160 | |||||||||||
518 (Feb 2024) | 0.263 | 0.051 | ||||||||||
519 (Mar 2024) | 0.325 | 0.094 | 0.018 | |||||||||
520 (Apr 2024) | 0.357 | 0.138 | 0.036 | 0.007 | ||||||||
521 (May 2024) | 0.368 | 0.178 | 0.057 | 0.014 | 0.003 | |||||||
522 (Jun 2024) | 0.364 | 0.211 | 0.082 | 0.024 | 0.006 | 0.001 | ||||||
523 (Jul 2024) | 0.349 | 0.237 | 0.107 | 0.036 | 0.010 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.329 | 0.255 | 0.132 | 0.051 | 0.016 | 0.004 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.305 | 0.266 | 0.154 | 0.067 | 0.023 | 0.007 | 0.002 | 0.000 | 0.000 | |||
526 (Oct 2024) | 0.279 | 0.270 | 0.175 | 0.085 | 0.033 | 0.011 | 0.003 | 0.001 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.253 | 0.270 | 0.191 | 0.102 | 0.043 | 0.015 | 0.005 | 0.001 | 0.000 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.227 | 0.264 | 0.205 | 0.119 | 0.055 | 0.021 | 0.007 | 0.002 | 0.001 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.172 | |||||||||||
518 (Feb 2024) | 0.277 | 0.059 | ||||||||||
519 (Mar 2024) | 0.337 | 0.107 | 0.023 | |||||||||
520 (Apr 2024) | 0.363 | 0.154 | 0.043 | 0.009 | ||||||||
521 (May 2024) | 0.367 | 0.194 | 0.069 | 0.018 | 0.004 | |||||||
522 (Jun 2024) | 0.357 | 0.227 | 0.096 | 0.031 | 0.008 | 0.002 | ||||||
523 (Jul 2024) | 0.336 | 0.250 | 0.123 | 0.046 | 0.014 | 0.003 | 0.001 | |||||
524 (Aug 2024) | 0.311 | 0.264 | 0.149 | 0.063 | 0.021 | 0.006 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.283 | 0.270 | 0.172 | 0.082 | 0.031 | 0.010 | 0.003 | 0.001 | 0.000 | |||
526 (Oct 2024) | 0.254 | 0.270 | 0.191 | 0.101 | 0.043 | 0.015 | 0.005 | 0.001 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.226 | 0.264 | 0.205 | 0.120 | 0.056 | 0.022 | 0.007 | 0.002 | 0.001 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.200 | 0.254 | 0.216 | 0.137 | 0.070 | 0.030 | 0.011 | 0.003 | 0.001 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.158 | |||||||||||
518 (Feb 2024) | 0.261 | 0.050 | ||||||||||
519 (Mar 2024) | 0.323 | 0.093 | 0.018 | |||||||||
520 (Apr 2024) | 0.356 | 0.136 | 0.035 | 0.007 | ||||||||
521 (May 2024) | 0.367 | 0.175 | 0.056 | 0.013 | 0.003 | |||||||
522 (Jun 2024) | 0.364 | 0.209 | 0.080 | 0.023 | 0.005 | 0.001 | ||||||
523 (Jul 2024) | 0.351 | 0.235 | 0.105 | 0.035 | 0.009 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.331 | 0.253 | 0.129 | 0.049 | 0.015 | 0.004 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.308 | 0.265 | 0.152 | 0.065 | 0.022 | 0.006 | 0.002 | 0.000 | 0.000 | |||
526 (Oct 2024) | 0.283 | 0.270 | 0.172 | 0.082 | 0.031 | 0.010 | 0.003 | 0.001 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.257 | 0.270 | 0.189 | 0.100 | 0.042 | 0.015 | 0.004 | 0.001 | 0.000 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.231 | 0.265 | 0.203 | 0.116 | 0.053 | 0.020 | 0.007 | 0.002 | 0.000 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.149 | |||||||||||
518 (Feb 2024) | 0.250 | 0.045 | ||||||||||
519 (Mar 2024) | 0.313 | 0.084 | 0.015 | |||||||||
520 (Apr 2024) | 0.350 | 0.125 | 0.030 | 0.005 | ||||||||
521 (May 2024) | 0.366 | 0.163 | 0.049 | 0.011 | 0.002 | |||||||
522 (Jun 2024) | 0.367 | 0.197 | 0.070 | 0.019 | 0.004 | 0.001 | ||||||
523 (Jul 2024) | 0.358 | 0.224 | 0.093 | 0.029 | 0.007 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.342 | 0.244 | 0.116 | 0.042 | 0.012 | 0.003 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.322 | 0.259 | 0.139 | 0.056 | 0.018 | 0.005 | 0.001 | 0.000 | 0.000 | |||
526 (Oct 2024) | 0.299 | 0.267 | 0.159 | 0.071 | 0.025 | 0.008 | 0.002 | 0.000 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.276 | 0.271 | 0.177 | 0.087 | 0.034 | 0.011 | 0.003 | 0.001 | 0.000 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.251 | 0.269 | 0.192 | 0.103 | 0.044 | 0.016 | 0.005 | 0.001 | 0.000 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.158 | |||||||||||
518 (Feb 2024) | 0.261 | 0.050 | ||||||||||
519 (Mar 2024) | 0.323 | 0.092 | 0.018 | |||||||||
520 (Apr 2024) | 0.356 | 0.136 | 0.035 | 0.007 | ||||||||
521 (May 2024) | 0.367 | 0.175 | 0.056 | 0.013 | 0.003 | |||||||
522 (Jun 2024) | 0.364 | 0.209 | 0.080 | 0.023 | 0.005 | 0.001 | ||||||
523 (Jul 2024) | 0.351 | 0.235 | 0.104 | 0.035 | 0.009 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.332 | 0.253 | 0.129 | 0.049 | 0.015 | 0.004 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.308 | 0.265 | 0.152 | 0.065 | 0.022 | 0.006 | 0.002 | 0.000 | 0.000 | |||
526 (Oct 2024) | 0.283 | 0.270 | 0.172 | 0.082 | 0.031 | 0.010 | 0.003 | 0.001 | 0.000 | 0.000 | ||
527 (Nov 2024) | 0.257 | 0.270 | 0.189 | 0.099 | 0.042 | 0.015 | 0.004 | 0.001 | 0.000 | 0.000 | 0.000 | |
528 (Dec 2024) | 0.232 | 0.266 | 0.203 | 0.116 | 0.053 | 0.020 | 0.007 | 0.002 | 0.000 | 0.000 | 0.000 | 0.000 |
Months | n | |||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
517 (Jan 2024) | 0.156 | |||||||||||
518 (Feb 2024) | 0.258 | 0.049 | ||||||||||
519 (Mar 2024) | 0.321 | 0.090 | 0.017 | |||||||||
520 (Apr 2024) | 0.354 | 0.133 | 0.033 | 0.006 | ||||||||
521 (May 2024) | 0.367 | 0.173 | 0.054 | 0.013 | 0.002 | |||||||
522 (Jun 2024) | 0.365 | 0.206 | 0.077 | 0.022 | 0.005 | 0.001 | ||||||
523 (Jul 2024) | 0.353 | 0.232 | 0.102 | 0.034 | 0.009 | 0.002 | 0.000 | |||||
524 (Aug 2024) | 0.334 | 0.251 | 0.126 | 0.047 | 0.014 | 0.004 | 0.001 | 0.000 | ||||
525 (Sep 2024 | 0.311 | 0.264 | 0.149 | 0.063 | 0.021 | 0.006 | 0.001 | 0.000 | ||||
526 (Oct 2024) | 0.287 | 0.270 | 0.169 | 0.080 | 0.030 | 0.009 | 0.003 | 0.001 | 0.000 | |||
527 (Nov 2024) | 0.261 | 0.270 | 0.186 | 0.096 | 0.040 | 0.014 | 0.004 | 0.001 | 0.000 | 0.000 | ||
528 (Dec 2024) | 0.236 | 0.267 | 0.201 | 0.113 | 0.051 | 0.019 | 0.006 | 0.002 | 0.000 | 0.000 | 0.000 | 0.000 |