An Alternative Source of Funding to Mitigate Flood Losses through Bonds: A Model for Pricing Flood Bonds in Indonesian Territory
Abstract
:1. Introduction
2. Related Research
- 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
4. The Model
4.1. Assumptions
- a.
- The term of a flood bond is years, where is a positive integer.
- b.
- All variables are modeled in , where is the set of states of the world, is the sigma-algebra on , is increasing filtration with , and is a probability measure (also called the probability set function) on . Suppose that is a random variable, and . Then, .
- c.
- The face value of the flood bond is , where this value is paid at the end of year .
- d.
- Coupons from flood bonds are paid annually in the amount of at the end of each year with
- 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 is the 1st dekadal. Then, the start of the dry season in year is the -th dekadal, and the start of the second rainy season in year is the -th dekadal. In other words, the duration of the first rainy season is , the duration of the dry season is , and the duration of the second rainy season is .
- h.
- The start of the dry season 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 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 -th dekadal in the rainy season in a year, denoted by with , are assumed to have identical probability distribution characteristics. In other words, for all . Then, maximum rainfall values over the -th dekadal in the dry season of a year, denoted by with , are assumed to have identical probability distribution characteristics. In other words, for all .
- k.
- The annual force of interest for all 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, , for all .
4.2. Maximum Rainfall Trigger Index Model
- a.
- When and are in the first rainy season interval, the value is equal to the maximum values of , , , and .
- b.
- When is in the first rainy season interval, and is in the dry season interval, the value is equal to the maximum value of , , , , , , , and .
- c.
- When is in the first rainy season interval, and is in the second rainy season interval, the value is equal to the maximum value of , , , , , , , , , , and .
- d.
- When and are in the dry season interval, the value is equal to the maximum values of , , , and .
- e.
- When is in the dry season interval, and is in the second rainy season interval, the value is equal to the maximum value of , , , , , , , and .
- f.
- When and are in the second rainy season interval, the value is equal to the maximum values of , , , and .
4.3. Flood Bond Price Model
5. Numerical Simulation
5.1. Data Description
- 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.
5.2. Fitting Rainfall and Force of Interest Data Distribution
- a.
- The Dagum distribution () 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 and , respectively. The critical values with significance level from the KS and AD tests are and , respectively. Therefore, this distribution was chosen to represent the distribution of maximum rainfall data in the dry season. Mathematically, this is written as follows:
- b.
- The Burr distribution 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 and , respectively. The critical values with significance level from the KS and AD tests are and , respectively. Therefore, this distribution was chosen to represent the distribution of maximum rainfall data in the rainy season. Mathematically, this is written as follows:
- c.
- The Inverse-Gaussian distribution 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 and , respectively. The critical values with significance level from the KS and AD tests are and , respectively. Therefore, this distribution was chosen to represent the annual force of interest data distribution. Mathematically, this is written as follows:
5.3. Forecasting Duration Ranges of Rainy and Dry Seasons
5.4. Estimating Flood Bond Price
- a.
- The average of the estimation results is 0.9339297. In other words, the estimated price of flood bonds in Bandung Regency is 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
7. Conclusions
Author Contributions
Funding
Data Availability Statement
Acknowledgments
Conflicts of Interest
References
- 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]
- 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]
- 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]
- 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]
- The World Bank Group; Asian Development Bank. Climate Risk Profile: Indonesia; Asian Development Bank: Mandaluyong, Philippines, 2021. [Google Scholar]
- 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]
- 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]
- 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]
- 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]
- Burnecki, K.; Giuricich, M.N. Stable Weak Approximation at Work in Index-Linked Catastrophe Bond Pricing. Risks 2017, 5, 64. [Google Scholar] [CrossRef]
- Lee, J.-P.; Yu, M.-T. Valuation of Catastrophe Reinsurance with Catastrophe Bonds. Insur. Math. Econ. 2007, 41, 264–278. [Google Scholar] [CrossRef]
- Cox, S.H.; Pedersen, H.W. Catastrophe Risk Bonds. North Am. Actuar. J. 2000, 4, 56–82. [Google Scholar] [CrossRef]
- 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]
- Albrecher, H.; Hartinger, J.; Tichy, R.F. QMC Techniques for CAT Bond Pricing. Monte Carlo Methods Appl. 2004, 10, 197–211. [Google Scholar] [CrossRef]
- 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]
- Jarrow, R.A. A Simple Robust Model for Cat Bond Valuation. Financ. Res. Lett. 2010, 7, 72–79. [Google Scholar] [CrossRef]
- Nowak, P.; Romaniuk, M. Pricing and Simulations of Catastrophe Bonds. Insur. Math. Econ. 2013, 52, 18–28. [Google Scholar] [CrossRef]
- Ma, Z.-G.; Ma, C.-Q. Pricing Catastrophe Risk Bonds: A Mixed Approximation Method. Insur. Math. Econ. 2013, 52, 243–254. [Google Scholar] [CrossRef]
- 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]
- Liu, J.; Xiao, J.; Yan, L.; Wen, F. Valuing Catastrophe Bonds Involving Credit Risks. Math. Probl. Eng. 2014, 2014, 1–6. [Google Scholar] [CrossRef]
- 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]
- 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]
- Georgiopoulos, N. Pricing Catastrophe Bonds with Multistage Stochastic Programming. Comput. Manag. Sci. 2017, 14, 297–312. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- 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]
- 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]
- Manathunga, V.; Deng, L. Pricing Pandemic Bonds under Hull–White and Stochastic Logistic Growth Model. Risks 2023, 11, 155. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- 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]
- Reijnen, R.; Albers, W.; Kallenberg, W.C.M. Approximations for Stop-Loss Reinsurance Premiums. Insur. Math. Econ. 2005, 36, 237–250. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- Hyndman, R.J. Moving Averages. In International Encyclopedia of Statistical Science; Springer: Berlin/Heidelberg, Germany, 2011; pp. 866–869. [Google Scholar]
- 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]
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. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
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
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 StyleSukono, 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 StyleSukono, 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