Next Article in Journal
Research on Runoff Management of Sponge Cities under Urban Expansion
Previous Article in Journal
Source of Methanogens and Characteristics of Methane Emission in Two Wastewater Treatment Plants in Xi’an
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Alternative Source of Funding to Mitigate Flood Losses through Bonds: A Model for Pricing Flood Bonds in Indonesian Territory

1
Department of Mathematics, Faculty of Mathematics and Natural Sciences, Universitas Padjadjaran, Sumedang 45363, Indonesia
2
Doctoral Program of Mathematics, Faculty of Mathematics and Natural Sciences, Universitas Padjadjaran, Sumedang 45363, Indonesia
3
Faculty of Informatics and Computing, Universiti Sultan Zainal Abidin, Besut 22200, Malaysia
*
Author to whom correspondence should be addressed.
Water 2024, 16(15), 2102; https://doi.org/10.3390/w16152102
Submission received: 22 June 2024 / Revised: 23 July 2024 / Accepted: 23 July 2024 / Published: 25 July 2024

Abstract

:
Indonesia suffers significant economic losses from floods, and state budget allocations are often inadequate. Flood bonds provide an alternative funding source, but the pricing framework is complex due to simultaneous flood and financial risk considerations. Therefore, this study aims to model flood bond prices as an alternative flood funding in Indonesia. The model is formulated using the risk-neutral-pricing measure with the stochastic assumption of the force of interest. The claim trigger is represented as maximum rainfall, which is modeled as a continuous-stochastic process with a discrete-time index. Given the varying patterns of rainy and dry seasons, we assume both durations are dynamic. Then, we provide the approximate model solution for the government to estimate bond prices quickly. This estimation shows that the bond’s trigger point is proportional to the bond prices. Additionally, bond prices are proportional to the dry season duration and inversely proportional to the rainy season duration. We also show that using a stochastic force of interest yields significant differences from a constant one except for the constant as data average. This study can help the Indonesian government price flood bonds and provide more tools for related meteorological and climatological institutions to calculate the probability of future maximum rainfall.

1. Introduction

Indonesia has a high rainfall danger with an annual average cumulative rainfall of 1500 to 4000 mm/year [1]. High rainfall generally occurs during the rainy season, which occurs from October to March on average. The geographical location at the equator, precisely between two continents (Asia and Australia) and two oceans (Indian and Pacific), causes sunlight throughout the year, high evaporation of seawater, and the rapid formation of rain clouds. Monsoon winds and diverse mountain topography also increase rainfall intensity through orographic processes.
The high rainfall danger in Indonesia makes this country vulnerable to weather catastrophes, especially floods. Based on data from The International Disaster Database (accessed on 2 February 2024 at https://public.emdat.be/) as visually presented in Figure 1, there are, on average, eight flood catastrophes yearly. This average is the second highest worldwide after China. The highest frequency of floods occurred in 2020, with 25 annual catastrophes, while the lowest was in 2004, with one flood catastrophe. According to Figure 1, the frequency of flood catastrophes tends to increase each year, as indicated by the positive slope of the linear trend line in Figure 1. Meteorologically, this increasing trend is due to the intensification of the La Niña phenomenon, which increases rainfall, especially during the rainy season. Increased tropical storm activity in the western Pacific Ocean also brings extreme rainfall [2]. Climatologically, global climate change has caused an increase in sea surface temperature by 0.8 °C since 1901. It increases evaporation and atmospheric moisture, resulting in higher and more intense rainfall in Indonesia [3].
Unfortunately, the increasing tendency of flood frequency in Indonesia is still slightly mitigated. Hence, the increasing frequency tendency positively correlates with the losses experienced [4], where one flood disaster causes average economic losses of 50.7 million dollars. In 2020 alone, economic losses due to flooding in Indonesia reached more than 2.1 billion USD or around 31.5 trillion IDR (exchange rate 1 USD = 15,000 IDR) (accessed on 2 February 2024 on National Disaster Management Agency of the Republic of Indonesia, https://data.bnpb.go.id/). These economic losses include damage to infrastructure, damage to homes, decreased agricultural productivity, and significant disruption to business activities. According to a report of The World Bank Group & Asian Development Bank [5], the intensity and frequency of extreme rainfall in Indonesia will increase further in the future because of climate change. It could cause annual economic losses to increase to 4.3 billion USD by 2030.
To fund flood mitigation and recovery, the Indonesian government has generally relied on two sources of funds: the state budget (central or regional) and social assistance (from the private sector, international institutions, communities, and other donors). These funds are usually used for infrastructure projects such as building embankments, normalizing rivers, drainage systems, and providing catastrophe management tools. However, this funding is often insufficient to overcome all significant flooding problems because the amount is limited [6,7]. Therefore, new sources of funding to deal with floods are needed.
Efforts to increase sources of funds to mitigate flood catastrophes in Indonesia have been started through insurance. Based on data from the International Disaster Database (Accessed on 18 July 2024 on https://public.emdat.be/), since 2000, one flood catastrophe has received insurance benefits of 7.2 million dollars on average. However, this benefit’s average amount is only 14.27 percent of the actual losses on average, so it is not significant enough to cover it. For example, in the 2013 Jakarta Flood, insurance benefits only covered 3.33 percent of actual losses, while in the 2014 Tangerang Flood, insurance benefits only covered 9.99 percent of actual losses. These data can also be seen in Figure 2. Therefore, the insurance mechanism still needs to be effective in funding flood mitigation in Indonesia, and other sources of funds are needed.
The World Bank and several other countries have tried to identify alternative sources of funding for catastrophe management (including floods) [8,9]. As a result of this identification, insurance-linked securities (ILSs) became a new mechanism [10,11]. Through this mechanism, the community is invited to play a role in flood funding as ILS investors. Communities as investors not only pursue financial profits but also support social goals to increase preparedness and resilience to catastrophes. Given the high risk of flooding, they will be given significant compensation from this investment in return. Among existing securities, bonds are the most successful securities in ILS. Bonds have a clear interest (coupon) distribution schedule, fast liquidity, and moderate financial risk. Investors and insureds prefer these advantages in this scheme. These are, after this, referred to as flood bonds. Many countries have issued flood bonds. Under historical data on flood bond issuance taken on 1 March 2024, from ARTEMIS (https://www.artemis.bm/deal-directory), the UK used catastrophe bonds worth 150 million USD to transfer flood risks in April 2007. Then, in May 2018, Japan used catastrophe bonds worth 320 million USD in funding measures to mitigate the consequences of typhoons, floods, and earthquakes. The United States used 300 million USD from flood bonds in April 2019.
Based on the previous paragraph, flood bonds can be used as an alternative source of funds for flood management in Indonesia. However, to issue it, one thing that must be conducted is to study the pricing framework in Indonesia. The flood bond pricing framework is complex to formulate because prices are simultaneously influenced by financial and flood risk variables, e.g., interest rates, rainfall, duration of the rainy season, and duration of the dry season. Therefore, this research aims to model flood bond prices in Indonesia. Considering that the rainy season and dry season patterns in Indonesia are not always constant, we consider it in the designed model. To model the trigger event, we use the maximum rainfall trigger index to quantify it. This index has the advantage that the measurement is quickly carried out, objective, and transparently based on the flood catastrophe that occurred. These advantages minimize administrative costs in assessing claims and the potential for moral hazard in claims. We model this index as a continuous stochastic process with a discrete-time index. Then, we model flood bond prices using a risk-neutral-pricing measure assuming a stochastic force of interest. After the modeling was conducted, we applied the model to rainfall and interest rate data in Indonesia. Finally, we also analyze the influence of rainfall and financial risks on flood bond prices, e.g., trigger point of rainfall, bond term, force of interest, and time to buy the bond. This study can increase the tools available for calculating the probability of maximum rainfall in the future for meteorological or climatological institutions in this country. Apart from that, this study can certainly be used by flood bond issuers in Indonesia in determining flood bond prices that better reflect the real situation based on the inconstant duration of rainy and dry seasons and the stochastic force of interest.

2. Related Research

Before reviewing articles on flood bond pricing modeling, we describe some research on catastrophe bond pricing modeling of general catastrophes. This research became popular in the 2000s. Cox and Pedersen [12] used equilibrium pricing theory to develop a method for pricing catastrophe bonds based on a model of interest rate structure and catastrophe risk probability. This model has been used as a basis for modeling in subsequent research. Lee and Yu [13] developed a catastrophe bond pricing model considering moral hazards, basis risks, and stochastic interest rates. Albrecher et al. [14] examined quasi-Monte Carlo methods and variance reduction algorithms for pricing catastrophe bonds involving integrals with discontinuous and infinite-dimensional integrands. Zimbidis et al. [15] modeled earthquake catastrophe bond prices using an extreme value theory (EVT) approach and presented numerical results using Monte Carlo simulations and stochastic iterative equations. Jarrow [16] provided a simple solution for a catastrophe bond (CAT) pricing model using a model for the evolution of the Libor term structure of interest rates.
Research on modeling catastrophe bond prices has become increasingly popular since 2011. Nowak and Romaniuk [17] developed the model of Cox and Pedersen [12] by describing interest rate dynamics using the Cox–Ingersoll–Ross and Hull–White models. Ma and Ma [18] applied the approximation method introduced by Chaubey et al. [19] to determine the solution of the model by Lee and Yu [13]. Liu et al. [20] used the Jarrow and Turnbull model and EVT methods to model catastrophe bond prices with credit risk involved. Ma et al. [21] modeled catastrophe bond prices using a doubly stochastic Poisson process with intensity function proposed by Black et al. [22]. Georgiopoulos [23] proposed a stochastic optimization model for pricing catastrophe bonds. Giuricich and Burnecki [24] modeled catastrophe bond prices using a left-truncated heavy-tailed distribution of a nonhomogeneous compound Poisson process. Deng et al. [25] modeled catastrophe bond prices for drought catastrophe financing using a homogeneous compound Poisson process. Then, Ibrahim et al. [26] modeled catastrophe bond prices with an aggregate loss trigger index and a number of fatalities using a nonhomogeneous compound Poisson process. Sukono et al. [27] developed the model of Ibrahim et al. [26] involving the inflation rate and the correlation rate between trigger indices. Both are modeled using vector autoregression and copula models. Tang et al. [28] developed the model of Ibrahim et al. [26] by modeling the trigger indices using the peaks-over-threshold (POT) method. Then, Manathunga and Deng [29] modeled pandemic bond prices using the stochastic logistic growth model.
After providing a brief previous research overview on modeling catastrophe bond prices in general, we review research with the main topic of modeling flood bond prices until 4 June 2024. After searching in the Scopus, Springer, and Taylor databases, four previous studies discussed this topic. Chen et al. [30] modeled flood bond prices with a single-period maturity. They used the trigger indemnity (actual loss) index to determine the trigger event. They designed the index using a homogeneous compound Poisson process. The model they studied did not have a closed-form solution, so the solution was approximated by the Wang-double-factor transformation method. Once the model was designed, they applied it using data on extreme flood losses in China from 1961 to 2009. In the application, they found that the trigger point of flood bonds was proportional to their price.
Chao and Zou [31] modeled multi-period flood bond prices using parametric trigger and indemnity indices. The parametric and indemnity trigger indices were represented by the number of fatalities and aggregate losses from extreme flooding, respectively. They designed the marginal distributions of both indices using a combination of a homogeneous compound Poisson process and peaks-over threshold method. Then, the joint distribution of the two indices was modeled with an Archimedean copula. To model flood bond prices, they used a risk-neutral-pricing measure and a stochastic interest rate based on the Cox–Ingersoll–Ross (CIR) model. The model they designed did not have a closed-form solution, so the solution was approached using the Monte Carlo method. After designing the model, they applied it to flood loss and fatality data provided by Dartmouth College since 1981. In the application, they found that flood intensity, period, and correlation between indices had an inverse relationship with flood bond prices.
Li et al. [32] designed a flood bond pricing model with a trigger index for the number of dangerous floods that occur. Dangerous flood was defined as flooding accompanied by rainfall that exceeds a specified threshold value. The number of these dangerous floods was modeled with compound distributions. This compound distribution consisted of the dangerous flood frequency distribution represented by the Poisson distribution and the rainfall distribution represented by the generalized Pareto distribution (GPD). The flood bond pricing model was modeled using a risk-neutral-pricing measure and stochastic interest rates based on the Vasicek model. The bond price model did not have a closed-form solution, so the solution was approximated by the Wang transformation method. The model was then applied to flood catastrophe data in China from 2004 to 2019, which was accessed in the Ministry of Water Resources of the People’s Republic of China.
Ibrahim et al. [33] designed a flood bond pricing model with a closed-form solution. A closed-form solution was obtained because they assumed the interest rates were constant. The model was designed with an indemnity trigger index. The index was modeled using a homogeneous compound Poisson process. To estimate the solution of the model, they applied the Normal, Normal-Power 2, Gamma, and Gamma Inverse-Gaussian approximation methods described by Chaubey et al. [19] and Reijnen et al. [34]. After the model was designed, they applied it to flood data in Indonesia from 2010 to 2021, which was accessed at the National Catastrophe Management Agency of the Republic of Indonesia. On that application, they found that the loss an investor might experience on a flood bond was half the possible gain.
We identify the gaps based on a review of the research carried out. Then, we made these gaps novelties for our research. These novelties are as follows:
a.
This study uses a parametric trigger index represented by the maximum rainfall in each dekadal (10 days) each year.
b.
This study considers the duration of the rainy and dry seasons each year, which are assumed to be not constant. This inconstancy is adjusted to the actual conditions regarding the varying patterns of the country’s rainy and dry seasons.
c.
This study examines the effect of buying time on flood bond prices so that investors can use it as a guide about when it is the right time to buy bonds based on the claim probability and price they expect.

3. Flood Bond: Simply and Brief Explanation

Flood bonds, one type of catastrophe bond, are financial securities linked to flood risk belonging to a sponsor (government or other institution, insurer, or reinsurer). This linking is conducted by sponsors by sharing some or all the flood risk they bear with investors in the capital market. These linking aims to obtain a broader coverage capacity to face the risk of enormous flood losses. The enormous risk of major flood losses provides significant compensation to investors of these securities in return. If the claim trigger event for a flood bond does not occur within its life span, the investor receives in full of all periodic coupons and the redemption value. However, they must be prepared to risk losing some or all of it if the triggering event occurs [35].
In simple terms, flood bonds have three main entities: sponsors, special purpose vehicles (SPVs), and investors. The sponsor is the insured (government or another institution, insurer, or reinsurer) of the flood bond. They share flood risks with investors in the capital market through bonds, form special purpose vehicles (SPVs) that will issue bonds, and pay premiums to the SPV. Once an SPV is established and issued bonds, they collect bond funds from investors and premiums from sponsors. These funds are then managed by investing them in short-term safe securities, e.g., treasury bonds. The investment results will later be paid as coupons to investors.
In flood bonds, the trigger event is measured by various indices. Several indexes are commonly used, namely the indemnity index, parametric index, model index, and industry index. The indemnity index is measured by actual losses from flooding experienced by sponsors. Then, the parametric index is measured through specific parameters related to flood catastrophes, e.g., maximum rainfall and the number of fatalities. Meanwhile, the model index is measured based on risk analysis results from risk modeling agencies, e.g., Risk Management Solution (RMS) and Applied Insurance Research (AIR) Worldwide. Finally, the industry index is measured by the losses incurred by the industry, e.g., the issuer. Combinations between trigger indices are also possible. This study uses a parametric index, maximum rainfall, as the trigger index of flood bonds. This index has the advantage that it measures quickly, objectively, and transparently based on the flood that occurred [36,37,38]. These advantages minimize administrative costs in assessing claims and the potential for moral hazard in claims.

4. The Model

4.1. Assumptions

The following are the assumptions used in modeling the flood bond prices in this study:
a.
The term of a flood bond is T years, where T is a positive integer.
b.
All variables are modeled in ( Ω , F , F k , P ) , where Ω is the set of states of the world, F is the sigma-algebra on Ω , F k F is increasing filtration with k 0 , T , and P is a probability measure (also called the probability set function) on F . Suppose that A is a random variable, and A F . Then, P A A = P ω Ω A ω A .
c.
The face value of the flood bond is K , where this value is paid at the end of year T .
d.
Coupons from flood bonds are paid annually in the amount of C at the end of each year t with t = 1 ,   2 , ,   T .
e.
One month is divided into three dekadals: the 1st to the 10th is called the first dekadal, the 11th to the 20th is called the second dekadal, and the 21st to the end of the month is called the third dekadal. In other words, there are 36 dekadals in one year.
f.
Every year begins with the first rainy season. Then, the season continues with the dry season. Finally, the year ends with the second rainy season. It is adjusted to the duration of rainy and dry seasons in Indonesia in reality [39]. For other countries, it can be reformulated.
g.
The duration of rainy and dry seasons is discrete with dekadal units. The start of the first rainy season in year t is the 1st dekadal. Then, the start of the dry season in year t is the d t -th dekadal, and the start of the second rainy season in year t is the r t -th dekadal. In other words, the duration of the first rainy season is 1,2 , , d t 1 , the duration of the dry season is d t , d t + 1 , , r t 1 , and the duration of the second rainy season is r t , r t + 1 , , 36 .
h.
The start of the dry season d t is determined as the 1st dekadal with cumulative rainfall of fewer than fifty millimeters and is followed by the next two dekadals, each of which is less than fifty millimeters of cumulative precipitation. Then, the start of the second rainy season r t is determined as the 1st dekadal with cumulative rainfall of more than or equal to fifty millimeters and is followed by the next two dekadals, each of which is more than or equal to fifty millimeters of cumulative precipitation.
i.
Maximum rainfall at each dekadal does not affect each other. In reality, this assumption may not hold. However, assumptions are made to facilitate the modeling process, likewise with the use of assumptions on the following random variable.
j.
Maximum rainfall values over the m -th dekadal in the rainy season in a year, denoted by R m with m = 1 ,   2 , , d t 1 , r t , r t + 1 , ,   36 , are assumed to have identical probability distribution characteristics. In other words, P R m u = P R u = F R u , for all m = 1 ,   2 , , d t 1 , r t , r t + 1 , ,   36 . Then, maximum rainfall values over the n -th dekadal in the dry season of a year, denoted by D n with n = d t ,   d t + 1 , , r t 1 , are assumed to have identical probability distribution characteristics. In other words, P D n u = P D u = F D u , for all n = d t ,   d t + 1 , , r t 1 .
k.
The annual force of interest I t for all t = 1 ,   2 ,   ,   T in one year is assumed to not affect the annual force of interest in other years and have the same distribution characteristics. In other words, P I t i = P I i = F I i , for all t = 1 ,   2 ,   ,   T .

4.2. Maximum Rainfall Trigger Index Model

Maximum rainfall from the c 1 -th dekadal to c 2 -th dekadal in year t with 1 c 1 c 2 36 and t = 1 , 2 , , T is designed as follows:
a.
When c 1 and c 2 are in the first rainy season interval, the value is equal to the maximum values of R c 1 , R c 1 + 1 , , and R c 2 .
b.
When c 1 is in the first rainy season interval, and c 2 is in the dry season interval, the value is equal to the maximum value of R c 1 , R c 1 + 1 , , R d t 1 , D d t , D d t + 1 , , and D c 2 .
c.
When c 1 is in the first rainy season interval, and c 2 is in the second rainy season interval, the value is equal to the maximum value of R c 1 , R c 1 + 1 , , R d t 1 , D d t , D d t + 1 , , D r t 1 , R r t , R r t + 1 , , and R c 2 .
d.
When c 1 and c 2 are in the dry season interval, the value is equal to the maximum values of D c 1 , D c 1 + 1 , , and D c 2 .
e.
When c 1 is in the dry season interval, and c 2 is in the second rainy season interval, the value is equal to the maximum value of D c 1 , D c 1 + 1 , , D r t 1 , R r t , R r t + 1 , , and R c 2 .
f.
When c 1 and c 2 are in the second rainy season interval, the value is equal to the maximum values of R c 1 , R c 1 + 1 , , and R c 2 .
Mathematically, we denote it by M c 1 , c 2 , t formulated as follows:
M c 1 , c 2 , t = max R c 1 , , R c 2 : c 1 = 1 , 2 , d t 1 , c 2 = 1 , 2 , , d t 1 max R c 1 , , R d t 1 , D d t , , D c 2 : c 1 = 1 , 2 , d t 1 , c 2 = d t , d t + 1 , , r t 1 max R c 1 , , R d t 1 , D d t , , D r t 1 , R r t , , R c 2 : c 1 = 1 , 2 , d t 1 , c 2 = r t , r t + 1 , , 36 max D c 1 , , D c 2 : c 1 = d t , d t + 1 , , r t 1 , c 2 = d t , d t + 1 , , r t 1 max D c 1 , , D r t 1 , R r t , , R c 2 : c 1 = d t , d t + 1 , , r t 1 , c 2 = r t , r t + 1 , , 36 max R c 1 , , R c 2 : c 1 = r t , r t + 1 , , 36 , c 2 = r t , r t + 1 , , 36
After M c 1 , c 2 , t is defined, we then formulate its cumulative distribution function as follows:
P M c 1 , c 2 , t u = P max R c 1 , , R c 2 u : c 1 = 1 , 2 , d t 1 , c 2 = 1 , 2 , , d t 1 P max R c 1 , , R d t 1 , D d t , , D c 2 u : c 1 = 1 , 2 , d t 1 , c 2 = d t , d t + 1 , , r t 1 P max R c 1 , , R d t 1 , D d t , , D r t 1 , R r t , , R c 2 u : c 1 = 1 , 2 , d t 1 , c 2 = r t , r t + 1 , , 36 P max D c 1 , , D c 2 u : c 1 = d t , d t + 1 , , r t 1 , c 2 = d t , d t + 1 , , r t 1 P max D c 1 , , D r t 1 , R r t , , R c 2 u : c 1 = d t , d t + 1 , , r t 1 , c 2 = r t , r t + 1 , , 36 P max R c 1 , , R c 2 u : c 1 = r t , r t + 1 , , 36 , c 2 = r t , r t + 1 , , 36 .
Note that for each m = 1 ,   2 , , d t 1 , r t , r t + 1 , , 36 and n = d t ,   d t + 1 , , r t 1 , it holds as follows [40]:
max R m , D n u R m u , D n u .
Then, because the maximum rainfall values in each dekadal do not affect each other, the following equation holds:
P R m u , D n u = P R m u P D n u ,
for each m = 1 ,   2 , , d t 1 , r t , r t + 1 , , 36 and n = d t ,   d t + 1 , , r t 1 . Based on Equations (2) and (3), the CDF of M c 1 , c 2 , t in Equation (1) can be rewritten as follows:
P M c 1 , c 2 , t u = P R c 1 u P R c 2 u : c 1 = 1 ,   2 , d t 1 ,   c 2 = 1 ,   2 , , d t 1                                             P R c 1 u P R d t 1 u P D d t u P D c 2 u : c 1 = 1 ,   2 , d t 1 ,   c 2 = d t ,   d t + 1 , , r t 1                         P R c 1 u P R d t 1 u P D d t u P D r t 1 u P R r t u P R c 2 u : c 1 = 1 ,   2 , d t 1 ,   c 2 = r t ,   r t + 1 , , 36                                           P D c 1 u P D c 2 u : c 1 = d t ,   d t + 1 , , r t 1 ,   c 2 = d t ,   d t + 1 , , r t 1 P D c 1 u P D r t 1 u P R r t u P R c 2 u : c 1 = d t ,   d t + 1 , , r t 1 ,   c 2 = r t ,   r t + 1 , , 36                   P R c 1 u P R c 2 u : c 1 = r t ,   r t + 1 , , 36 ,   c 2 = r t ,   r t + 1 , , 36                                     .
Then, recall that for every m , R m has an identical distribution, and for every n , D n also has an identical distribution. These are as mentioned in point j in Section 4.1. Therefore, Equation (4) can be rewritten as follows:
P M c 1 , c 2 , t u = F R c 2 c 1 + 1 u : c 1 = 1 , 2 , d t 1 , c 2 = 1 , 2 , , d t 1 F R d t c 1 u F D c 2 d t + 1 u : c 1 = 1 , 2 , d t 1 , c 2 = d t , d t + 1 , , r t 1 F R d t c 1 + c 2 r t + 1 u F D r t d t u : c 1 = 1 , 2 , d t 1 , c 2 = r t , r t + 1 , , 36 F D c 2 c 1 + 1 u : c 1 = d t , d t + 1 , , r t 1 , c 2 = d t , d t + 1 , , r t 1 F D r t c 1 u F R c 2 r t + 1 u : c 1 = d t , d t + 1 , , r t 1 , c 2 = r t , r t + 1 , , 36 F R c 2 c 1 + 1 u : c 1 = r t , r t + 1 , , 36 , c 2 = r t , r t + 1 , , 36
To determine the survival distribution function (SDF) of M c 1 , c 2 , t , we can use Equation (5) to obtain the following equation:
P M c 1 , c 2 , t > u = 1 P M c 1 , c 2 , t u .
Then, we formulate the maximum rainfall from the c 1 -th dekadal in year t 1 to the c 2 -th dekadal in year t 2 with 1 c 1 , c 2 36 and 1 t 1 < t 2 T as follows:
M c 1 , t 1 , c 2 , t 2 = max M c 1 , 36 , t 1 , M 1,36 , t 1 + 1 , , M 1 , c 2 , t 2 .
Conveniently, the CDF of M c 1 , t 1 , c 2 , t 2 can be formulated as follows:
P M c 1 , t 1 , c 2 , t 2 u = P max M c 1 , 36 , t 1 , M 1,36 , t 1 + 1 , , M 1 , c 2 , t 2 u , = P M c 1 , 36 , t 1 u , M 1,36 , t 1 + 1 u , , M 1 , c 2 , t 2 u , = P M c 1 , 36 , t 1 u P M 1,36 , t 1 + 1 u P M 1 , c 2 , t 2 u .
In other words, the CDF of M c 1 , t 1 , c 2 , t 2 is determined based on the CDF of M c 1 , c 2 , t in Equation (5). To determine the SDF of M c 1 , t 1 , c 2 , t 2 , we can use Equation (6) to obtain the following equation:
P M c 1 , t 1 , c 2 , t 2 > u = 1 P M c 1 , t 1 , c 2 , t 2 u .

4.3. Flood Bond Price Model

Suppose that a flood bond with a term of T years is issued at the beginning of the c 1 -th dekadal in the year t 1 . In this study, annual coupons are paid to investors at the end of each year. In other words, the t -th annual coupon is paid at every c c 1 -th dekadal in year t c 1 , t , where c c 1 and t c 1 , t are functions, each expressed as follows:
c c 1 = 36 : c 1 = 1 c 1 1 : c 1 = 2 , 3 , , 36 ,
and
t c 1 , t = t 1 + t 1 : c 1 = 1 , t = 1 , 2 , , T t 1 + t : c 1 = 2 , 3 , , 36 , t = 1 , 2 , , T .
Then, the redemption value of the flood bonds is paid to investors at the end of the year T. In other words, the redemption value is paid on the c c 1 -th dekadal in year t c 1 , T . The annual coupon payment scheme and redemption value are intact if the maximum rainfall up to the payment date does not exceed the trigger point and in proportion if the maximum rainfall up to the payment date exceeds the trigger point. Mathematically, the annual coupon payment scheme and redemption value are denoted as P t and are formulated as follows:
P t = C 1 M c 1 , t 1 , c c 1 , t c 1 , t ω ( u ) + λ 1 C 1 M c 1 , t 1 , c c 1 , t c 1 , t > ω ( u ) : t = 1 , 2 , , T 1 ( C + K ) 1 M c 1 , t 1 , c c 1 , t c 1 , t ω ( u ) + ( λ 1 C + λ 2 K ) 1 M c 1 , t 1 , c c 1 , t c 1 , t > ω ( u ) : t = T , = C 1 M c 1 , t 1 , 36 , t 1 + t 1 ω ( u ) + λ 1 C 1 M c 1 , t 1 , 36 , t 1 + t 1 > ω ( u ) : c 1 = 1 , t = 1 , 2 , , T 1 C 1 M c 1 , t 1 , c 1 1 , t 1 + t ω ( u ) + λ 1 C 1 M c 1 , t 1 , c 1 1 , t 1 + t > ω ( u ) : c 1 = 2 , 3 , , 36 , t = 1 , 2 , , T 1 ( C + K ) 1 M c 1 , t 1 , 36 , t 1 + t 1 ω ( u ) + ( λ 1 C + λ 2 K ) 1 M c 1 , t 1 , 36 , t 1 + t 1 > ω ( u ) : c 1 = 1 , t = T ( C + K ) 1 M c 1 , t 1 , c 1 1 , t 1 + t ω ( u ) + ( λ 1 C + λ 2 K ) 1 M c 1 , t 1 , c 1 1 , t 1 + t > ω ( u ) : c 1 = 2 , 3 , , 36 , t = T .
where ω is the trigger point, and λ 1 and λ 2 are, respectively, the proportion of annual coupon payment and redemption value when the maximum rainfall trigger point until the payment date is exceeded,
1 M c 1 , t 1 , 36 , t 1 + t 1 ω u = 1 : u M c 1 , t 1 , 36 , t 1 + t 1 ω 0 : u M c 1 , t 1 , 36 , t 1 + t 1 > ω ,
1 M c 1 , t 1 , c 1 1 , t 1 + t ω u = 1 : u M c 1 , t 1 , c 1 1 , t 1 + t ω 0 : u M c 1 , t 1 , c 1 1 , t 1 + t > ω ,
1 M c 1 , t 1 , 36 , t 1 + t 1 > ω u = 1 : u M c 1 , t 1 , 36 , t 1 + t 1 > ω 0 : u M c 1 , t 1 , 36 , t 1 + t 1 ω ,
and
1 M c 1 , t 1 , c 1 1 , t 1 + t > ω u = 1 : u M c 1 , t 1 , c 1 1 , t 1 + t > ω 0 : u M c 1 , t 1 , c 1 1 , t 1 + t ω .
To determine the price of flood bonds, we use the summation of the present value of the expected annual coupon and the expected redemption value based on the risk-neutral pricing measure Q . For example, now is the beginning of the c 1 -th dekadal. The present value of a payment of one currency unit in the next t years ahead is formulated as follows:
B c 1 , t = E Q e s = 1 t I s 1 c 1 36 + I s + 1 c 1 36 F c 1 36 ,
where E Q is the statistical expectation under Q (risk-neutral-pricing measure). In Equation (9), the form 1 c 1 36 indicates the time of compounding one currency unit with the force of interest I s in year s . Then, the form c 1 36 indicates the time of compounding one currency unit with the force of interest I s + 1 in year s + 1 . Based on the present value concept in Equation (9), the flood bond price is designed as the sum of the present values of expected P t with t = 1,2 , , T . Mathematically, it is formulated as follows:
V c 1 , t 1 , T = t = 1 T E Q P t e s = 1 t I s 1 c 1 36 + I s + 1 c 1 36 F c 1 36 .
Based on Cox and Pedersen [12], after converting the physical probability measure P to the risk-neutralized measure Q , it is assumed that the original distribution characteristics of the maximum rainfall process remain intact. Events dependent on rainfall risk variables and events dependent on financial variables have no effect on one another under the risk-neutralized measure Q (for more detail, see Cox and Pedersen [12]). Consequently, Equation (10) can be rewritten as follows:
V c 1 , t 1 , T = t = 1 T E P P t F t E Q e s = 1 t I s 1 c 1 36 + I s + 1 c 1 36 F c 1 36 = t = 1 T E P P t F c 1 36 B c 1 , t ,
where B c 1 , t is determined using Equation (9), and E P P t F t is determined based on Equations (6) and (7) as follows:
E P P t F t = C P M c 1 , t 1 , 36 , t 1 + t 1 ω + λ 1 C P M c 1 , t 1 , 36 , t 1 + t 1 > ω : c 1 = 1 , t = 1 , 2 , , T 1 C P ( M c 1 , t 1 , c 1 1 , t 1 + t ω ) + λ 1 C P M c 1 , t 1 , c 1 1 , t 1 + t > ω : c 1 = 2 , 3 , , 36 , t = 1 , 2 , , T 1 ( C + K ) P M c 1 , t 1 , 36 , t 1 + t 1 ω + ( λ 1 C + λ 2 K ) P M c 1 , t 1 , 36 , t 1 + t 1 > ω : c 1 = 1 , t = T ( C + K ) P ( M c 1 , t 1 , c 1 1 , t 1 + t ω ) + ( λ 1 C + λ 2 K ) P M c 1 , t 1 , c 1 1 , t 1 + t > ω : c 1 = 2 , 3 , , 36 , t = T .

5. Numerical Simulation

5.1. Data Description

The primary data used in this study are daily rainfall data in Bandung Regency, West Java Province, Indonesia, from 1 January 1982 to 31 December 2023. These data were obtained on 1 March 2024 openly at the Prediction of Worldwide Energy Resource (POWER) provided by the National Aeronautics and Space Administration (NASA) of the United States of America at the following link: https://power.larc.nasa.gov/data-access-viewer/. From these data, we processed it to produce the following data:
a.
Data on cumulative rainfall per dekadal in Bandung Regency, West Java Province, Indonesia, from the 1st dekadal in 1982 to the last dekadal in 2023. The data size is 1.512 and has units of mm/dekadal.
b.
Data for the start of the dry and the second rainy seasons from 1982 to 2023. Data have dekadal units. Both data have the same size, namely 42. The dry season starts in the 16th dekadal every year on average. Then, the earliest dry season occurs in the 8th dekadal and the latest occurs in the 22th dekadal. The second rainy season, on average, starts in the 30th dekadal every year. Then, the earliest the second rainy season occurs is in the 23th dekadal and the latest is in the 35th dekadal. A visualization of both data is given in Figure 3.
c.
Maximum rainfall data per dekadal in the rainy and dry seasons from the 1st dekadal in 1982 to the last dekadal in 2023. The data sizes for the rainy and dry seasons are 915 and 597, respectively. Both data are in units of mm/day. The average maximum rainfall in the dry season is 8.26 mm/day. Then, the maximum rainfall in the dry season is at least 0 mm/day and at most 64.68 mm/day. Meanwhile, the average maximum rainfall in the rainy season is 26.94 mm/day. Then, the maximum rainfall in the rainy season is at least 1.20 mm/day and at most 226.27 mm/day. To determine the shape of the distribution, we analyzed the skewness and kurtosis of the two datasets. The data skewness in the dry and rainy seasons is 1.81 and 2.48, respectively. Both skewness values are positive, which indicates that both data distributions are skewed to the left. Then, the data kurtosis in the dry and rainy seasons is 6.17 and 35.92, respectively. Both kurtosis values are more than three, which indicates that both data distributions are sharper than the normal distribution (leptokurtic). From the results of the skewness and kurtosis analysis of the two data, the conclusion is that the distribution form of the two data is a heavy-right-tail distribution. Visually, this shape can be seen in the histogram of the two data in Figure 3.
Apart from that, the second primary data are data on the annual effective interest rate in Indonesia from 1991 to 2023. These data are then transformed into a force of interest with the equation I t = ln ( 1 + I t ) , where I t is the force of annual interest and I t is the annual effective interest rate. The data size is 33 and was openly obtained on 1 March 2024 from Bank Indonesia at the following link: https://www.bi.go.id/id/statistik/indikator/BI-Rate.aspx. The average force of interest in each year is 0.092329. Then, the force of interest has an average deviation of 0.053137 from the average. The skewness and kurtosis of the data are 2.247412 and 7.234458, respectively. In other words, the data distribution is skewed to the left and is sharper than a normal (leptokurtic) curve. We can also say that the data distribution is a heavy-right-tail distribution. A histogram of the data is given in Figure 3.

5.2. Fitting Rainfall and Force of Interest Data Distribution

The distribution analysis results in Section 5.1 show that the distributions of maximum rainfall data per dekadal in the rainy and dry seasons are heavy-right-tailed. The distribution of force of interest data also has the same form. Therefore, this section will match the distribution of these data with a theoretical distribution with heavy-right-tailed characteristics. Many theoretical distributions have it, e.g., gamma, beta, Burr, Pareto, Dagum, and inverse Gaussian distributions. Before matching, the parameters of the theoretical distributions are estimated first. Several methods can be used in this assessment, e.g., the maximum likelihood (ML) method. Then, to check the match between the distribution of the three data and the theoretical one, several tests can be used, e.g., Kolmogorov–Smirnov (KS) and Anderson–Darling (AD) tests. The KS test is nonparametric and does not depend on certain distribution assumptions. The AD test is a modification of the KS test, which emphasizes the tail of the distribution. The theoretical distribution that most closely fits the empirical distribution has the lowest KS and AD statistical values, denoted by D α and A α , where α represents the significance level. The KS and AD statistical values are determined as follows [41]:
D = sup F n x F ( x ) : x X
and
A = 1 n n = 1 n 2 n 1 log F x n + log 1 F x n n + 1 n ,
where n represents the data size, F ( x ) is the theoretical distribution and F n x is the empirical distribution.
In this study, we estimated the parameters of theoretical distributions and fitted them to three data distributions using the EasyFit software version 5.5. In more detail, parameter estimation was carried out using the ML method, while fitness tests were carried out using KS and AD tests. In summary, using the software, we found the following:
a.
The Dagum distribution ( α = 4.2444 , β = 18.175 , κ = 0.15306 ) is the most fit compared to other theoretical distributions to describe the distribution of maximum rainfall data in the dry season. In more detail, α and κ represent shape parameters, while β is a scale parameter. The statistical values of the KS and AD tests are D = 0.02213 and A = 0.2513 , respectively. The critical values with significance level α = 0.05 from the KS and AD tests are D 0.05 = 0.05558 and A 0.05 = 2.5018 , respectively. Therefore, this distribution was chosen to represent the distribution of maximum rainfall data in the dry season. Mathematically, this is written as follows:
F D u = 1 + u β α κ = 1 + u 18.175 4.2444 0.15306 .
b.
The Burr distribution ( α = 2.8735 , β = 32.627 , κ = 1.9002 ) is the most fit compared to other theoretical distributions to describe the distribution of maximum rainfall data in the rainy season. In more detail, α and κ represent shape parameters, while β is a scale parameter. The statistical values of the KS and AD tests are D = 0.02114 and A = 0.38158 , respectively. The critical values with significance level α = 0.05 from the KS and AD tests are D 0.05 = 0.04489 and A 0.05 = 2.5018 , respectively. Therefore, this distribution was chosen to represent the distribution of maximum rainfall data in the rainy season. Mathematically, this is written as follows:
F R u = 1 1 + u β α κ = 1 1 + u 32.627 2.8735 1.9002 .
c.
The Inverse-Gaussian distribution ( λ = 0.27875 , μ = 0.09233 ) is the more fit compared to other theoretical distributions to describe the distribution of the annual force of interest data, where λ and μ are the shape and the mean parameters, respectively. The statistical values of the KS and AD tests are D = 0.0798 and A = 0.36667 , respectively. The critical values with significance level α = 0.05 from the KS and AD tests are D 0.05 = 0.2342 and A 0.05 = 2.5018 , respectively. Therefore, this distribution was chosen to represent the annual force of interest data distribution. Mathematically, this is written as follows:
F I i = Φ λ i i μ 1 + e 2 λ μ Φ λ i i μ + 1 , = Φ 0.27875 i i 0.09233 1 + e 2 ( 0.27875 ) 0.09233 Φ 0.27875 i i 0.09233 + 1 ,
where Φ represents the cumulative distribution function of the standard normal random variable.

5.3. Forecasting Duration Ranges of Rainy and Dry Seasons

This section contains the forecasting stages for the start of the dry and the second rainy seasons in the following years. Many forecasting models can be used. In this research, we use a simple moving-average time series model. The selection of this model is based on data that do not have extreme fluctuations (stable) [42]. This stability can be seen from the visualization of both data in Figure 3. This stability can also be shown visually through the Box and Whisker diagram of the two data in Figure 4, which shows no outliers.
Mathematically, the simple moving-average models of order p for forecasting the start of the dry and the second rainy seasons are formulated as follows [42]:
d t + 1 = 1 p j = 1 p d t j + 1 + e t
and
r t + 1 = 1 p j = 1 p r t j + 1 + e t ,
where t = p + 1 , p + 2 , , T , T represents data size, e t and e t represents error. To determine the order of p , we conduct experiments for p = 1 ,   2 ,   3 ,   4 , and 5 . The one selected is the one with the best accuracy. In this study, accuracy is measured by mean-absolute-percentage error (MAPE), which is calculated using the following equation:
M A P E p = 1 T p t = p + 1 T r t r ^ t r t × 100 %
or
M A P E p = 1 T p t = p + 1 T d t d ^ t d t × 100 % ,
where r ^ t and d ^ t are the forecasts of r t and d t , respectively. In summary, we find that the MAPE of the single moving-average model of orders 1, 2, 3, 4, and 5 of r_t is 25%, 21.7%, 20%, 18%, and 18.3%, respectively. Then, the MAPE of the single moving average model of order 1, 2, 3, 4, and 5 of d t is 12.13%, 10.4%, 9.31%, 9.33%, and 9.4%, respectively. Based on these MAPE values, r t and d t are modeled with a single moving average of order 4 and 3, respectively.

5.4. Estimating Flood Bond Price

This section contains estimates of flood bond prices. In this study, the flood bond pricing model in Equation (11) does not have a closed-form solution. In more detail, this form results from the absence of a solution to B c 1 , t in the equation. Therefore, we determine the solution numerically using the Monte Carlo method. Because only the form B c 1 , t does not have a closed-form solution, computing the numerical solution of Equation (11) can be conducted quickly. Numerically, the solution of the flood bond price model in Equation (11) can be approximated by the following equation:
V ^ c 1 , t 1 , T = 1 W w = 1 W t = 1 T E P P t F c 1 36 e s = 1 t i s , w 1 c 1 36 + i s + 1 , w c 1 36 ,
where W represents the number of simulations, and i s , w is the random value of I s generated in the w -th simulation with s = 1,2 , , t and w = 1 ,   2 , , W . To start estimation, the values of the variables in Equation (13) are determined first. We assume that flood bonds are issued at the beginning of the c 1 = 1 st dekadal in year t 1 = 2023 and have a term of T = 3 years. Then, the annual coupon and redemption value are assumed to be C = 0.1 USD and K = 1 USD, respectively. Then, the trigger point for the maximum rainfall index is assumed to be ω = 100 mm/day. When the trigger point is exceeded, the annual coupon and redemption value are paid in the proportions λ 1 = 0 and λ 2 = 0.5 , respectively. By conducting W = 500,000 simulations, we visually provide all the estimation results in the histogram in Figure 5. The description of all estimation results is as follows:
a.
The average of the estimation results is 0.9339297. In other words, the estimated price of flood bonds in Bandung Regency is V ^ 1,2023,3 = 0.9339297 USD.
b.
The standard deviation of the estimation results is 0.07329027. In other words, the average deviation of flood bond price estimates from the mean is 0.07329027 USD.
c.
The range of the estimated results is 0.7066545. This value is so large. Hence, many simulations are needed to obtain the no bias estimation. It is why we conducted 500,000 simulations.

6. Discussion

This section contains an analysis of the effect of rainfall and financial risks on flood bond prices. The rainfall and financial risks analyzed include trigger points, term, stochastic and constant forces of interest, and time to buy bonds. The first is analyzing the effect of trigger points on flood bond prices. To simplify the analysis, we assume that all variables except the trigger point are the same as those used in Section 5.4. Meanwhile, we assume the trigger point is a set U = { 50 ,   75 ,   100 ,   125 ,   150 ,   175 ,   200 } with units of mm/day. For each member of U , we determine the price of the flood bond corresponding to it. Briefly, we visually provide flood bond prices for each U member in Figure 6.
Figure 6 illustrates that trigger points positively correlate with flood bond prices. The higher the trigger point, the higher the price of flood bonds, and the lower the trigger point, the lower the price of flood bonds. This result is validated logically, where if the trigger point is high, the chance of claiming from flood bonds will be low. This results in high investor demand for these bonds. Hence, the price becomes expensive. Then, if the trigger point is low, the chance of a claim from flood bonds will be high. This results in low investor demand for these bonds. Hence, the price becomes cheap. The positive correlation results between trigger points and flood bond prices obtained in this study align with the results conducted by Chao and Zou [31] and Ibrahim et al. [33].
Next is an analysis of the effect of the stochastic and constant force of interest on flood bond prices. We assume all variable values used are the same as in Section 5.4 except for the term. We assume the term is a set T = { 1 ,   2 ,   3 ,   4 ,   5 ,   6 } with units of year(s). For each member of T , we determine the price of a flood bond assuming a stochastic and constant force of interest. After that, we determine the differences. In more detail, the constant force of interest is assumed to be 0.09233, which is the average of the data from 1990 to 2023. Briefly, we visually provide the results of the analysis in Figure 7.
Figure 7 illustrates that the estimated price of flood bonds, assuming a stochastic and constant force of interest, significantly differs if the estimated constant value is far from its annual average. However, the choice of a constant force of interest in Equation (13) of its annual average is very reliable to use because the difference is not significant, although this is not analytical proof. Next is an analysis of the effect of the term of a flood bond on its price. This analysis can be carried out through Figure 7, where it appears that the term of the flood bond has a negative correlation with its price. The longer the term of a flood bond, the lower the price, and the shorter the term of a flood bond, the higher the price. It is validated logically in the same way as traditional bonds [43], where if the term of the flood bond is long, the chance of a claim occurring will be high. This results in low investor demand for these bonds. Hence, the price becomes cheap. Then, if the term of the flood bond is short, the chance of a claim occurring will be low. This results in high investor demand for these bonds. Hence, the price becomes expensive.
Next is an analysis of the effect of the buying time of flood bonds on its price. To analyze it, we assume that the variable values in the model are the same as in Section 5.4 except for c 1 . We let c 1 be the set C 1 = 1 ,   2 , ,   36 . For each c 1 in C 1 , we calculate the price of the flood bond corresponding to it. The flood bond price calculation used two force of interest assumptions: stochastic and constant. In summary, the analysis results are given visually in Figure 8.
Figure 8 illustrates that flood bonds tend to have the same price at the 1st to 16th dekadal, increase at the 17th dekadal, tend to be the same again until the 28th dekadal, increase at the 29th dekadal, and tend to be the same again until the final dekadal. After this was analyzed further, we found the cause. There is a difference in the duration of the dry and rainy seasons at each buying time for the flood bonds. Flood bonds with the more extended dry season duration have the highest prices and vice versa. This cause can be visually seen in Figure 8. The more extended duration of the dry season is related to the buying time of flood bonds in the 29th to 36th dekadals. As a result, bond prices in this dekadal range are the highest compared to other dekadals. Then, the shortest duration of the dry season is related to the buying time of flood bonds in the 1st to 16th dekadals. As a result, flood bond prices in this dekadal range have the lowest prices compared to other dekadals. Rationally, it is validated because the longer the dry season duration, the lower the chance of high rainfall. Then, the lower the chance, the higher the demand for flood bonds. Hence, the price becomes high. On the contrary, when the duration of the dry season is short, the same logic applies. Then, flood bonds with a long rainy season duration have low bond prices, and vice versa. This can be visually seen in Figure 8. The longest duration of the rainy season is owned by the buying time of flood bonds in the 1st to 16th dekadals. As a result, flood bond prices in this dekadal range have the lowest prices compared to other dekadals. Then, the shortest duration of the rainy season is owned by the buying time of flood bonds in the 29th to 36th dekadals. As a result, flood bond prices in this dekadal range are the highest compared to other dekadals. Rationally, it is validated because the longer the duration of the rainy season, the higher the chance of heavy rainfall. Then, the higher the chance, the lower the demand for flood bonds. As a result, the price is low. On the contrary, when the duration of the rainy season is short, the same logic applies.

7. Conclusions

This study aims to model flood bond prices in Indonesia, which assumes that the durations of the rainy and dry seasons are not constant. This inconstancy is adjusted to the actual conditions regarding the varying patterns of the country’s rainy and dry seasons. To model the trigger event, we use the maximum rainfall trigger index to quantify it. This index has the advantage that it measures quickly, objectively, and transparently based on the flood catastrophe that occurred. These advantages minimize administrative costs in assessing claims and the potential for moral hazard in claims. We model this index as a continuous stochastic process with a discrete-time index. Then, we model flood bond prices using a risk-neutral-pricing measure assuming a stochastic force of interest.
After the model was designed, the numerical simulation was conducted using rainfall and interest rate data in Indonesia. The simulation shows that the maximum rainfall trigger point is proportional to the price of flood bonds. In other words, the higher the trigger point, the more expensive the flood bond will be, and vice versa. It is validated logically, where if the trigger point is high, the chance of claiming from flood bonds will be low. This results in high investor demand for these bonds. Hence, the price becomes expensive. This trigger point sensitivity to flood bond prices aligns with the results conducted by Chao and Zou [31] and Ibrahim et al. [33]. Then, using a stochastic force of interest provides a significant difference compared to using a constant force of interest if the estimated constant value is far from the annual average of the data. Suppose a constant force of interest is still to be used. In that case, the annual average is better used for results closer to estimates with a stochastic force of interest, even though this is not analytically evident. Then, the duration of the rainy season is inversely proportional to the price of flood bonds, while the duration of the dry season is proportional to the price of flood bonds. In other words, the longer the dry season, the lower the chance of heavy rainfall. Then, the lower the chance, the higher the demand for flood bonds. Hence, the price becomes high. On the contrary, the longer the rainy season, the higher the chance of heavy rainfall. Then, the higher the chance, the lower the demand for flood bonds. Hence, the price becomes low.
This study can increase the tools available for calculating the probability of maximum rainfall in the future for related meteorological or climatological institutions in this country. Apart from that, flood bond issuers in Indonesia can undoubtedly use this study to determine flood bond prices that better reflect the actual situation in the form of inconstant duration of rainy and dry seasons and stochastic force of interest. The designed model has shortcomings that can be used in future research. The designed model does not consider the intensity of flood catastrophes. This consideration can better describe the flood conditions in Indonesia. Then, our model still needs to consider the inflation rate. By involving it, investors’ return estimates will be more accurate because they have considered reducing factors.

Author Contributions

Conceptualization, S. and R.A.I.; methodology, M.H. and R.A.I.; software, M.H. and R.A.I.; validation, S., M.H. and J.N.; formal analysis, R.A.I. and S.; investigation, S.; resources, S.; data curation, J.N.; writing—original draft preparation, M.D.J. and J.N.; writing—review and editing, N.Z.; visualization, M.D.J.; supervision, S.; project administration, M.H.; funding acquisition, S. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC for this paper were funded by Universitas Padjadjaran, grant number 2123/UN6.3.1/PT.00/2024.

Data Availability Statement

Data are contained within the article.

Acknowledgments

Thanks to Universitas Padjadjaran for providing the RDPD grant for the 2024 fiscal year.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Nirwansyah, A.W.; Braun, B. Tidal Flood Risk on Salt Farming: Evaluation of Post Events in the Northern Part of Java Using a Parametric Approach. Geosciences 2021, 11, 420. [Google Scholar] [CrossRef]
  2. Wang, B.; Xiang, B.; Lee, J.-Y. Subtropical High Predictability Establishes a Promising Way for Monsoon and Tropical Storm Predictions. Proc. Natl. Acad. Sci. USA 2013, 110, 2718–2722. [Google Scholar] [CrossRef] [PubMed]
  3. Cai, W.; Borlace, S.; Lengaigne, M.; van Rensch, P.; Collins, M.; Vecchi, G.; Timmermann, A.; Santoso, A.; McPhaden, M.J.; Wu, L.; et al. Increasing Frequency of Extreme El Niño Events Due to Greenhouse Warming. Nat. Clim. Chang. 2014, 4, 111–116. [Google Scholar] [CrossRef]
  4. Laksonoa, B.C.; Wulansari, I.Y.; Permatasari, N. Small Area Estimation of Household Expenditure on Insurance Programs for Minimizing the Impact of Natural Disasters in West Java, Indonesia. Stat. J. IAOS 2023, 39, 729–743. [Google Scholar] [CrossRef]
  5. The World Bank Group; Asian Development Bank. Climate Risk Profile: Indonesia; Asian Development Bank: Mandaluyong, Philippines, 2021. [Google Scholar]
  6. Setiawan, E.P.; Wutsqa, D.U.; Abadi, A.M.; Kusuma, E. Pricing Indonesian Earthquake Catastrophe Bond Based on Depth and Magnitude. In Proceedings of the the 3rd International Conference on Science, Mathematics, Environment, and Education, AIP Conference Proceedings, Surakarta, Indonesia, 27–28 July 2021; pp. 1–9. [Google Scholar]
  7. Gunardi; Setiawan, E.P. Valuation of Indonesian Catastrophic Earthquake Bonds with Generalized Extreme Value (GEV) Distribution and Cox-Ingersoll-Ross (CIR) Interest Rate Model. In Proceedings of the 1st International Conference on Actuarial Science and Statistics, AIP Conference Proceedings, Bandung, Indonesia, 21–23 October 2014; pp. 1–14. [Google Scholar]
  8. Ibrahim, R.A.; Sukono; Napitupulu, H.; Ibrahim, R.I. Earthquake Bond Pricing Model Involving the Inconstant Event Intensity and Maximum Strength. Mathematics 2024, 12, 786. [Google Scholar] [CrossRef]
  9. Sukono; Napitupulu, H.; Riaman; Ibrahim, R.A.; Johansyah, M.D.; Hidayana, R.A. A Regional Catastrophe Bond Pricing Model and Its Application in Indonesia’s Provinces. Mathematics 2023, 11, 3825. [Google Scholar] [CrossRef]
  10. Burnecki, K.; Giuricich, M.N. Stable Weak Approximation at Work in Index-Linked Catastrophe Bond Pricing. Risks 2017, 5, 64. [Google Scholar] [CrossRef]
  11. Lee, J.-P.; Yu, M.-T. Valuation of Catastrophe Reinsurance with Catastrophe Bonds. Insur. Math. Econ. 2007, 41, 264–278. [Google Scholar] [CrossRef]
  12. Cox, S.H.; Pedersen, H.W. Catastrophe Risk Bonds. North Am. Actuar. J. 2000, 4, 56–82. [Google Scholar] [CrossRef]
  13. Lee, J.; Yu, M. Pricing Default-Risky CAT Bonds with Moral Hazard and Basis Risk. J. Risk Insur. 2002, 69, 25–44. [Google Scholar] [CrossRef]
  14. Albrecher, H.; Hartinger, J.; Tichy, R.F. QMC Techniques for CAT Bond Pricing. Monte Carlo Methods Appl. 2004, 10, 197–211. [Google Scholar] [CrossRef]
  15. Zimbidis, A.A.; Frangos, N.E.; Pantelous, A.A. Modeling Earthquake Risk via Extreme Value Theory and Pricing the Respective Catastrophe Bonds. ASTIN Bull. 2007, 37, 163–183. [Google Scholar] [CrossRef]
  16. Jarrow, R.A. A Simple Robust Model for Cat Bond Valuation. Financ. Res. Lett. 2010, 7, 72–79. [Google Scholar] [CrossRef]
  17. Nowak, P.; Romaniuk, M. Pricing and Simulations of Catastrophe Bonds. Insur. Math. Econ. 2013, 52, 18–28. [Google Scholar] [CrossRef]
  18. Ma, Z.-G.; Ma, C.-Q. Pricing Catastrophe Risk Bonds: A Mixed Approximation Method. Insur. Math. Econ. 2013, 52, 243–254. [Google Scholar] [CrossRef]
  19. Chaubey, Y.P.; Garrido, J.; Trudeau, S. On the Computation of Aggregate Claims Distributions: Some New Approximations. Insur. Math. Econ. 1998, 23, 215–230. [Google Scholar] [CrossRef]
  20. Liu, J.; Xiao, J.; Yan, L.; Wen, F. Valuing Catastrophe Bonds Involving Credit Risks. Math. Probl. Eng. 2014, 2014, 1–6. [Google Scholar] [CrossRef]
  21. Ma, Z.; Ma, C.; Xiao, S. Pricing Zero-Coupon Catastrophe Bonds Using EVT with Doubly Stochastic Poisson Arrivals. Discret. Dyn. Nat. Soc. 2017, 2017, 3279647. [Google Scholar] [CrossRef]
  22. Black, F.; Derman, E.; Toy, W. A One-Factor Model of Interest Rates and Its Application to Treasury Bond Options. Financ. Anal. J. 1990, 46, 33–39. [Google Scholar] [CrossRef]
  23. Georgiopoulos, N. Pricing Catastrophe Bonds with Multistage Stochastic Programming. Comput. Manag. Sci. 2017, 14, 297–312. [Google Scholar] [CrossRef]
  24. Giuricich, M.N.; Burnecki, K. Modelling of Left-Truncated Heavy-Tailed Data with Application to Catastrophe Bond Pricing. Phys. A Stat. Mech. its Appl. 2019, 525, 498–513. [Google Scholar] [CrossRef]
  25. Deng, G.; Liu, S.; Li, L.; Deng, C. Research on the Pricing of Global Drought Catastrophe Bonds. Math. Probl. Eng. 2020, 2020, 3898191. [Google Scholar] [CrossRef]
  26. Ibrahim, R.A.; Sukono; Napitupulu, H. Multiple-Trigger Catastrophe Bond Pricing Model and Its Simulation Using Numerical Methods. Mathematics 2022, 10, 1363. [Google Scholar] [CrossRef]
  27. Sukono; Ibrahim, R.A.; Saputra, M.P.A.; Hidayat, Y.; Juahir, H.; Prihanto, I.G.; Halim, N.B.A. Modeling Multiple-Event Catastrophe Bond Prices Involving the Trigger Event Correlation, Interest, and Inflation Rates. Mathematics 2022, 10, 4685. [Google Scholar] [CrossRef]
  28. Tang, Y.; Wen, C.; Ling, C.; Zhang, Y. Pricing Multi-Event-Triggered Catastrophe Bonds Based on a Copula–POT Model. Risks 2023, 11, 151. [Google Scholar] [CrossRef]
  29. Manathunga, V.; Deng, L. Pricing Pandemic Bonds under Hull–White and Stochastic Logistic Growth Model. Risks 2023, 11, 155. [Google Scholar] [CrossRef]
  30. Chen, J.; Liu, G.; Yang, L.; Shao, Q.; Wang, H. Pricing and Simulation for Extreme Flood Catastrophe Bonds. Water Resour. Manag. 2013, 27, 3713–3725. [Google Scholar] [CrossRef]
  31. Chao, W.; Zou, H. Multiple-Event Catastrophe Bond Pricing Based on CIR-Copula-POT Model. Discret. Dyn. Nat. Soc. 2018, 2018, 5068480. [Google Scholar] [CrossRef]
  32. Li, J.; Cai, Z.; Liu, Y.; Ling, C. Extremal Analysis of Flooding Risk and Its Catastrophe Bond Pricing. Mathematics 2022, 11, 114. [Google Scholar] [CrossRef]
  33. Ibrahim, R.A.; Sukono, S.; Napitupulu, H.; Ibrahim, R.I.; Johansyah, M.D.; Saputra, J. Estimating Flood Catastrophe Bond Prices Using Approximation Method of the Loss Aggregate Distribution: Evidence from Indonesia. Decis. Sci. Lett. 2023, 12, 179–190. [Google Scholar] [CrossRef]
  34. Reijnen, R.; Albers, W.; Kallenberg, W.C.M. Approximations for Stop-Loss Reinsurance Premiums. Insur. Math. Econ. 2005, 36, 237–250. [Google Scholar] [CrossRef]
  35. Loubergé, H.; Kellezi, E.; Gilli, M. Using Catastrophe-Linked Securities to Diversify Insurance Risk: A Financial Analysis of Cat Bonds. J. Insur. Issues 1999, 22, 125–146. [Google Scholar]
  36. Anggraeni, W.; Supian, S.; Sukono; Halim, N.B.A. Earthquake Catastrophe Bond Pricing Using Extreme Value Theory: A Mini-Review Approach. Mathematics 2022, 10, 4196. [Google Scholar] [CrossRef]
  37. Zhang, N.; Huang, H. Assessment of World Disaster Severity Processed by Gaussian Blur Based on Large Historical Data: Casualties as an Evaluating Indicator. Nat. Hazards 2018, 92, 173–187. [Google Scholar] [CrossRef]
  38. Finken, S.; Laux, C. Catastrophe Bonds and Reinsurance: The Competitive Effect of Information-Insensitive Triggers. J. Risk Insur. 2009, 76, 579–605. [Google Scholar] [CrossRef]
  39. Monika, P.; Ruchjana, B.N.; Abdullah, A.S. GSTARI-X-ARCH Model with Data Mining Approach for Forecasting Climate in West Java. Computation 2022, 10, 204. [Google Scholar] [CrossRef]
  40. Salvadori, G.; De Michele, C.; Kottegoda, N.T.; Rosso, R. Univariate Extreme Value Theory. In Extremes in Nature; Springer: Dordrecht, The Netherlands, 2007; Volume 56, pp. 1–112. [Google Scholar]
  41. Purwani, S.; Ibrahim, R.A. Using Simple Fixed-Point Iterations to Estimate Generalized Pareto Distribution Parameters. IAENG Int. J. Appl. Math. 2024, 54, 194–204. [Google Scholar]
  42. Hyndman, R.J. Moving Averages. In International Encyclopedia of Statistical Science; Springer: Berlin/Heidelberg, Germany, 2011; pp. 866–869. [Google Scholar]
  43. Kwok, Y.-K. Interest Rate Derivatives: Bond Options, LIBOR and Swap Products. In Mathematical Models of Financial Derivatives; Springer: Berlin/Heidelberg, Germany, 2008; pp. 441–505. [Google Scholar]
Figure 1. Frequency of annual flood catastrophes in Indonesia from 2000 to 2023.
Figure 1. Frequency of annual flood catastrophes in Indonesia from 2000 to 2023.
Water 16 02102 g001
Figure 2. Comparison of insurance benefits and actual losses due to floods in Indonesia since 2000 and several case examples.
Figure 2. Comparison of insurance benefits and actual losses due to floods in Indonesia since 2000 and several case examples.
Water 16 02102 g002
Figure 3. Data visualization of the start of dry and second rainy seasons from 1982 to 2023 (a), histogram of maximum daily rainfall data per dekadal in the dry (b) and rainy (c) seasons, and histogram of annual force of interest data (d).
Figure 3. Data visualization of the start of dry and second rainy seasons from 1982 to 2023 (a), histogram of maximum daily rainfall data per dekadal in the dry (b) and rainy (c) seasons, and histogram of annual force of interest data (d).
Water 16 02102 g003
Figure 4. Box and Whisker diagram from data for the start of the dry (a) and the second rainy (b) seasons.
Figure 4. Box and Whisker diagram from data for the start of the dry (a) and the second rainy (b) seasons.
Water 16 02102 g004
Figure 5. Histogram of the results of 500,000 simulations of flood bond price estimations.
Figure 5. Histogram of the results of 500,000 simulations of flood bond price estimations.
Water 16 02102 g005
Figure 6. Visualization of the effect of trigger points on flood bond prices.
Figure 6. Visualization of the effect of trigger points on flood bond prices.
Water 16 02102 g006
Figure 7. Visualization of flood bond prices with the assumption of stochastic and constant force of interest (a) and the absolute difference among them (b).
Figure 7. Visualization of flood bond prices with the assumption of stochastic and constant force of interest (a) and the absolute difference among them (b).
Water 16 02102 g007
Figure 8. Visualization of the effect of buying time (a), duration of the dry season (b), and duration of the rainy season (c) on flood bond prices assuming a stochastic and constant force of interest.
Figure 8. Visualization of the effect of buying time (a), duration of the dry season (b), and duration of the rainy season (c) on flood bond prices assuming a stochastic and constant force of interest.
Water 16 02102 g008aWater 16 02102 g008b
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sukono; Hidayanti, M.; Nahar, J.; Ibrahim, R.A.; Johansyah, M.D.; Zamri, N. An Alternative Source of Funding to Mitigate Flood Losses through Bonds: A Model for Pricing Flood Bonds in Indonesian Territory. Water 2024, 16, 2102. https://doi.org/10.3390/w16152102

AMA Style

Sukono, Hidayanti M, Nahar J, Ibrahim RA, Johansyah MD, Zamri N. An Alternative Source of Funding to Mitigate Flood Losses through Bonds: A Model for Pricing Flood Bonds in Indonesian Territory. Water. 2024; 16(15):2102. https://doi.org/10.3390/w16152102

Chicago/Turabian Style

Sukono, Monika Hidayanti, Julita Nahar, Riza Andrian Ibrahim, Muhamad Deni Johansyah, and Nurnadiah Zamri. 2024. "An Alternative Source of Funding to Mitigate Flood Losses through Bonds: A Model for Pricing Flood Bonds in Indonesian Territory" Water 16, no. 15: 2102. https://doi.org/10.3390/w16152102

APA Style

Sukono, Hidayanti, M., Nahar, J., Ibrahim, R. A., Johansyah, M. D., & Zamri, N. (2024). An Alternative Source of Funding to Mitigate Flood Losses through Bonds: A Model for Pricing Flood Bonds in Indonesian Territory. Water, 16(15), 2102. https://doi.org/10.3390/w16152102

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop