Multiple-Trigger Catastrophe Bond Pricing Model and Its Simulation Using Numerical Methods

: Investor interest in single-trigger catastrophe bonds (STCB) has the potential to decline in the future. It is triggered by the increasing trend of global catastrophe loss and intensity every year, which increases the probability that a claim of STCB will occur. To increase investor interest again, the issuance of multiple-trigger catastrophe bonds (MTCB) can be one solution. However, to issue MTCB, its pricing is more complex because it involves more factors than STCB. Therefore, this study aims to design a simple MTCB pricing model. The claim trigger indices used are actual loss and fatality. Then, a nonhomogeneous compound Poisson process is used to model actual losses and fatalities aggregate to consider catastrophe intensity. In addition, this study proposes numerical methods, namely the continuous distribution approximation method and the Nuel recursive method, to facilitate the application of the model. Finally, an analysis of the effect of catastrophe intensity and other factors on MTCB prices is also presented. This study is expected to help special-purpose vehicles as MTCB issuers in MTCB pricing.


Introduction
Catastrophe intensity in all countries worldwide has tended to increase in recent decades [1,2]. This trend of increasing catastrophe intensity is predicted to continue in the future [3]. To deal with it, the government of a country needs to improve its preparedness [4], one of which is the preparedness of catastrophe management funds. However, many countries have limited catastrophe management funds because they only rely on the budget and social aid [5]. Therefore, new fund sources to improve preparedness catastrophe management funds are needed.
To obtain new fund sources for preparedness catastrophe management, sponsors (governments, insurers, or reinsurers) can transfer their catastrophe risk to the capital market. A special mechanism is used here. It involves the role of investors to jointly bear the catastrophe risk of sponsors [6][7][8]. The involvement of investors is assessed by associating catastrophe risk with financial securities in the capital market. The associated catastrophe risk is generally in the form of a single index, namely the loss index or parameter index [9]. The loss index represents the actual catastrophe loss experienced by the insured or by the property claim service (PCS), while the parameter index represents the magnitude of the earthquake and the number of fatalities [10,11]. Among the financial securities in the capital market, bonds are essential securities and successfully relate the catastrophe risk of sponsors to investors [12][13][14][15][16][17][18][19]. For sponsors, the funds obtained from bonds are significant and quickly raised, and for investors, the yields generated from these bonds are high. These bonds are then referred to as catastrophe bonds.
Pricing is an essential step in issuing a catastrophe bond. Several studies have focused on modeling single-trigger catastrophe bond prices (STCB) with a loss trigger index. Chernobai et al. [20] designed an STCB price model in which the aggregate of the catastrophe losses is modeled by a nonhomogeneous compound Poisson process (NCPP) whose cumulative distribution function is left-truncated. Then, in their study, Ma and Ma [21] proposed a numerical method, a mixed approximation method, to design the STCB price model. Schmidt [22] carried out an STCB price modeling in which the aggregate of catastrophe losses is designed using a shot-noise process with time-nonhomogeneous drivers. Ma et al. [23] modeled a zero-coupon STCB price whose catastrophe intensity and aggregate catastrophe losses are designed using the Black Derman Toy model and NCPP. Then, Tang and Yuan [24] integrated a product pricing measure and a distorted probability measure to model the STCB price. Deng et al. [25] use the peaks over threshold (POT) method to design a drought STCB price model. Several studies also use parameter indices to model STCB prices. Zimbidis et al. [26] modeled an earthquake STCB pricing using the earthquake magnitude index. The model's maximum earthquake magnitude during the STCB period is designed with an extreme value theory (EVT) approach. Romaniuk [27] proposed using the Monte Carlo method and the iterative stochastic equation method to determine the earthquake STCB price with an earthquake magnitude index. Kurniawan et al. [28] applied a jump-diffusion process to model the flood STCB price by considering the number of flooded areas index.
Given that loss and the intensity of global catastrophes tend to increase every year, investor interest in purchasing STCB can decline in the future because the probability that a claim of STCB will occur is increasing. Therefore, efforts to attract investors to purchase catastrophe bonds are necessary. One of the efforts is to issue a multiple-trigger catastrophe bond (MTCB) [29,30]. MTCB can be more attractive to investors because it will lose its principal if two events of trigger claims occur. Another advantage is that MTCB can reduce the moral hazard of the sponsor [5]. In addition, MTCB can describe the actual catastrophe severity better than STCB. However, there are obstacles to issuing MTCB, where pricing is more complex than STCB [31]. It is because MTCB involves more factors in determining its price. Therefore, a study on MTCB price modeling needs to be conducted.
Based on the literature review that has been completed, two studies focus on MTCB price modeling. These studies were carried out by Reshetar [31] and Chao and Zou [5]. The MTCB claim trigger indices used are the actual loss and fatality indices. To describe the risk distribution, Reshetar [31] modeled it parametrically, while Chao and Zou [5] modeled it using a homogeneous compound Poisson process (HCPP). Then, the joint risk distribution of both models is designed with a copula approach. Finally, the use of both models is carried out using the Monte Carlo simulation.
The gap from the MTCB pricing model designed by Reshetar [31] is that the catastrophe intensity factor is not considered, even though this factor may affect MTCB prices. Then, there is no closed-form solution to the model, making it challenging to determine the model solution. Furthermore, Chao and Zou's [5] model has included the catastrophe intensity factor. It overcomes the gap from Reshetar's [31] model. However, another gap emerges: the catastrophe intensity is assumed to be constant. It is not appropriate because the catastrophe intensity tends to increase every year. In addition, the risk distribution of fatalities in this model is designed with a continuous distribution, even though this risk is a discrete random variable, so it is not appropriate. Finally, there is no closed-form solution in this model, so it is also challenging to determine the solution of its model.
Based on the gap related to the MTCB price modeling described, the MTCB pricing modeling that involves the increasing catastrophe intensity factor currently is very open to being performed. In addition, a tool to determine the model solution so that it can be determined easily also needs to be developed. Therefore, this study aims to design an MTCB pricing model that involves the increasing catastrophe intensity factor and develop a method to make its solution easier to determine. The MTCB price modeled in this study is zero-coupon MTCB and coupon-paying MTCB with an annual term. The actual loss and fatality indices are used as claim trigger indices on the MTCB. The NCPP is used to model the aggregate of actual losses and fatalities to account for the increasing catastrophe intensity. Then, this study proposes using the new numerical methods to find solutions to the model, namely the continuous distribution approximation method and the Nuel recursive method. The continuous distribution approximation method is used to calculate the cumulative distribution function (CDF) value of actual loss aggregate, while the Nuel recursive method is used to calculate the CDF value of fatality aggregate. Then, after we obtained the model, a simulation of its use is performed on storm catastrophe data in the United States. This simulation analyzes how the catastrophe intensity and other factors influence MTCB prices. This model has the advantage of involving the increasing catastrophe intensity, and the search for the model solution has been facilitated by the previously mentioned numerical method. Therefore, this model is expected to help special-purpose vehicles as MTCB issuers in pricing MTCB. The simulations provide a reference for investors in choosing the MTCB they want to purchase based on the influence of catastrophe intensity and other factors on its price.

A Brief Explanation of Multiple-Trigger Catastrophe Bonds
Catastrophe bonds are bonds that are associated with catastrophe risk [32]. The issuer will obtain funds from the catastrophe bond if a predetermined claim trigger event occurs within the bond's life span [33,34]. In addition to a single-trigger index, catastrophe bonds can also use an independent multiple-trigger index [30]. A catastrophe bond that uses a multiple-trigger index is called a multiple-trigger catastrophe bond (MTCB). In the future, MTCB is potentially more attractive to investors than single-trigger catastrophe bonds (STCB). The increasing trend of global catastrophe loss and intensity every year causes the probability that a claim of STCB will occur higher. It is detrimental to investors because the possibility of losing the principal is also higher. On MTCB, the investor's principal is lost when two claim triggering events occur. It certainly benefits investors because moral hazard can be reduced so that the two events that trigger claims on MTCB can occur more slowly or even not occur at all. In addition, MTCB can also better describe the catastrophe severity that occurs than STCB. It is clear because the catastrophe severity in MTCB is measured by multiple-trigger indices, whereas in STCB, it is only measured by a singletrigger index.
A simple schematic of MTCB is described in this paragraph. The sponsor (government, insurer, or reinsurer) establishes a special-purpose vehicle (SPV) that acts as the issuer of the MTCB. Generally, SPV issues MTCB for one to five years [35,36]. Apart from acting as an MTCB issuer, SPV also acts as a claim fund provider for sponsors. After the SPV is established, the sponsor signs a protection contract and pays a premium to the SPV. The SPV then invests the proceeds from the sale of MTCB and sponsor premiums into short-term safe securities such as treasury bonds. The return from the investment is deposited into a trust account. These returns are then exchanged for floating returns based on the London Interbank Offered Rate (LIBOR) to immunize sponsors and investors from interest rate risk and default risk [10,35]. The funds raised usually exceed the required amount and have been arranged so that the payment of a claim for sponsors and coupons and MTCB principal for investors is guaranteed. If one of the trigger events occurs within the MTCB's life span, the investor will not receive the coupon back from that period until the maturity date. Then, if two trigger events occur within the MTCB's life span, the investor will receive the principal parts or not receive them at all [37]. Then, if the two trigger events do not occur within the MTCB's life span, the investor will receive the coupon and principal in total [5,31].

MCTB Price Modeling Assumptions
The following are the MTCB price modeling assumptions used in this study:

Annual Force of Interest
This study uses the annual force of interest to determine the present value of the annual coupon and principal of MTCB. The present value of one unit in year is expressed as follows [38][39][40][41]: where represents the annual force of interest in year , = 1,2, … , . In this study, the annual force of interest in year is modeled using the autoregressive integrated moving average (ARIMA) model. The ARIMA model of is stated as follows [42]: where , , and respectively represent autoregressive, differentiation, and movingaverage order, ∇ represents the differentiation operator, where ∇ = ∇(∇ −1 ) given ∇ = − −1 , and represents a random error. In Equation (2), the ARIMA model can also write as ARIMA ( , , ). The assumptions in modeling the annual force of interest with ARIMA are as follows: (a) is independent and identically normally distributed with zero mean and constant variance (~. . . (0, 2 )).

Aggregate of Actual Losses and Fatalities
The process of the number of catastrophes until time , { : ∈ [0, ]}, in this study is defined as a nonhomogeneous Poisson process with catastrophe intensity . Meanwhile, the aggregate of actual losses and fatalities are modeled by the nonhomogeneous compound Poisson process (NCPP), and both are assumed to be independent. The NCPP for the aggregate of actual losses is stated as follows: where { , = 1,2, … , } is a sequence of independent and identically distributed random variables of actual catastrophe losses with the distribution function cumulative (CDF) and probability density functions (PDF), which are expressed as (•) and (•), respectively. Note that { , = 1,2, … , } and { : ∈ [0, ]} are assumed to be independent. The CDF of is expressed as follows: where * (•) represents the convolution function of the -th fold of , which is expressed as follows: where * 0 ( ) = { 1 ; ≥ 0 0 ; < 0 , and * 1 ( ) = ( ) [43]. Then, the NCPP for the aggregate of fatalities is stated as follows: where { , = 1,2, … , } is a sequence of independent and identically geometric distributed random variables of catastrophe fatalities with the CDF and probability mass functions (PMF), which are expressed as (•) and (•), respectively. Note that { , = 1,2, … , } and { : ∈ [0, ]} are assumed to be independent. The CDF of is expressed as follows: where * (•) represents the convolution function of the -th fold of , which is expressed as follows: where * 0 ( ) = { 1 ; ≥ 0 0 ; < 0 , and * 1 ( ) = ( ) [43].
The first time the threshold value of the actual loss aggregate is exceeded is expressed as follows: where is the threshold value of the actual loss aggregate . Note that ≤ is equivalent to > , and > is equivalent to ≤ [44]. The first time the threshold value of the fatality aggregate is exceeded is expressed as follows: where is the threshold value of the fatality aggregate . Note that ≤ is equivalent to > , and > is equivalent to ≤ . The first time when one of the threshold values of the aggregate of actual losses and fatalities is exceeded is expressed as follows: The first time when both threshold values of the aggregate of actual losses and fatalities are exceeded is expressed as follows:

MTCB Price Modeling
The MTCB price modeled in this study is zero-coupon MTCB and coupon paying MTCB. First, the zero-coupon MTCB price is modeled in advance. If both threshold values of the aggregate of actual losses and fatalities occur within the MTCB's life span, then the principal paid to investors is the same as the proportion, whereas otherwise, the principal paid to investors is intact. Mathematically, the structure of principal payment ( ) to investors on the maturity date is expressed as follows: where represents the principal, is the proportion of principal payments, represents the year of maturity, and represents the first time both threshold values of the aggregate of actual losses and fatalities occurred. The zero-coupon MTCB pricing model is designed as the present value of the expectations. Mathematically, the zerocoupon MTCB price is expressed as follows: Proof. See Appendix A. □ Next is coupon-paying MTCB price modeling. If one of the threshold values of the aggregate of actual losses and fatalities occurs in a year within the MTCB's life span, then the coupon for that year until the maturity date will not be paid. If both threshold values of the aggregate of actual losses and fatalities occur within the MTCB's life span, the principal will be paid in proportion. If both threshold values of the aggregate of actual losses and fatalities do not occur within the MTCB's life span, the coupon and principal will be paid in full. The structure of coupon payments in year ( ) is mathematically expressed as follows: where represents the coupon, and represents the first time one of the threshold values of the aggregate of actual losses and fatalities occurred. The principal structure on the maturity date ( ) is the same as shown in Equation (13). The coupon-paying MTCB pricing model is designed as the sum of the present value of the expectations in each year in MTCB's life span and the present value of the expectation on the maturity date. Mathematically, the coupon-paying MTCB price ( ′) is expressed as follows: Proof. See Appendix B. □

The Numerical Methods in the MTCB Pricing Model
Calculation of the price of MTCB using the model in Equations (14) and (16) is difficult to do analytically. In general, the calculation of the CDF value of the aggregate of actual losses and fatalities is very complex, except for some cases where there are already closed solutions. This complexity is found in the calculation of the convolution function values * (•) and * (•) for high -folds. Therefore, numerical methods can be a tool to solve this problem. In this study, the CDF value of and is calculated using the continuous distribution approximation method and the Nuel recursive method, respectively.

The Continuous Distribution Approximation Method
The continuous distribution approximation method is the approach used to calculate the CDF value of the continuous NCPP [43]. Two of the continuous distributions used to approximate the distribution of the NCPP are the inverse-Gaussian (IG) distribution and the gamma-inverse-Gaussian (GIG) distribution. These two approaches were proposed by Chaubey et al. [45], and Reijnen et al. [46] later developed a rule of thumb for both.
Suppose that the skewness of on the aggregate of actual losses is expressed as ( ). Then, suppose that the mean, variance, skewness, and kurtosis of are expressed as ( ), ( ), ( ), and ( ), respectively. The PDF of can be approximated by the PDF of the IG distribution if ( ) ∈ (5, 15) and ( ) ∈ (1.5, 50). The PDF of approximated by the PDF of the IG distribution ( (•)) is expressed as follows:

The Nuel Recursive Method
The Nuel recursive method proposed by Nuel [47] is a numerical method used to calculate the CDF value of NCPP whose individual random variables are geometrically distributed with parameter ∈ (0,1) . For = 0 , (0) = − , and for = 1 , (b) Determine the value of with the following equation: (c) Determine the value of with the following equation: (d) Determine the value of and . If ≥ 0, then = −1 , and = .
(e) Determine the value of ( ). If = , then the process is complete, and If ≠ , then is restated as = + 1, and the process returns to step (b).

Data Description
The data used for the simulation of the MTCB price model are as follows:

ARIMA Parameter Estimation of the Annual Force of Interest Data and Its Forecasting
The stationarity assumption of is tested first. The test is carried out using the Augmented Dickey-Fuller (ADF) test. Briefly, with a significance level of 0.05, the test statistic value of , −3.4135, is greater than the critical value, −3.5560. Therefore, is not stationary, so is differentiated. differentiation is performed so that ∇ is obtained. Briefly, the new test statistic value of ∇ , −3.9161, is smaller than the critical value, −3.4135. Therefore, ∇ is stationary. The next step is to determine the autoregressive order and the moving-average order . The order and are determined using the partial autocorrelation function (PACF) and autocorrelation function (ACF) diagrams. Lags cut off from the PACF and ACF diagrams are selected as autoregressive order and moving-average order, respectively. In this study, the maximum lag considered is ten. The PACF and ACF diagrams from lag one to lag ten are presented in Figure 1.  Figure 1 shows that PACF and ACF are cut off at lags 8 and 3, respectively. It indicates that the and orders are 8 and 3, respectively. Therefore, the ARIMA (8, 1, 3) model is chosen to candidate the best ARIMA models. Next, the parameters of the model are estimated. This parameter estimation is completed by using the maximumlikelihood (ML) method. The parameter estimation results are presented in Table 1.
After the estimated parameters of the ARIMA model (8, 1, 3) presented in Table 1 were obtained, a test of the assumption that ~. . . (0, 2 ) is carried out. The test that is independent is carried out by the Ljung-Box (LB) test, while the test that ~(0, 2 ) is carried out by the Kolmogorov-Smirnov (KS) test. Briefly, with a significance level of 0.01, the LB and KS test statistic values are 1.0626 and 0.1685, respectively. Both test statistics are smaller than critical values, namely 8.1132 and 0.2653, respectively. Therefore, ~. . . (0, 2 ) . After that, checking the feasibility of the ARIMA (8, 1, 3) model based on the error is carried out. The error checked is the mean absolute percentage error (MAPE). The MAPE of the ARIMA model (8, 1, 3) is 17.6181%. This MAPE value is in the 10% to 20% interval, indicating that the ARIMA (8, 1, 3) model is good for practical forecasting of annual forces of interest in the next few years [48]. The results of force of interest forecasts for the next five years, obtained practically using the ARIMA (8, 1, 3) model, are presented in Table 2.

Fitting Data Distributions of Catastrophe Loss and Fatality to Theoretical Distributions
The fitting test is carried out first between the actual catastrophe loss data distribution and the theoretical distribution. The actual catastrophe loss data distribution generally has a heavy tail [21]. Therefore, the actual catastrophe loss data distribution is fitted to the heavy-tailed theoretical distribution. Some of the heavy-tailed theoretical distributions considered are presented in Table 3.
Parameter estimation of theoretical distributions is completed first. The method used is the ML method. The parameter estimation results of all theoretical distributions are presented in Table 4. After the parameter estimators of the theoretical distributions are obtained, the fitting test between the actual catastrophe loss data distribution and the theoretical distributions is carried out using the KS test. The level of significance used is 0.01. As shown in Table 5, the KS test statistic value of each theoretical distribution appears to be smaller than the critical value, 0.1407. It indicates that all theoretical distributions are suitable for approaching the actual catastrophe loss data distribution. However, the most suitable among the theoretical distributions must be selected. The most suitable theoretical distribution is the smallest KS test statistic value. The Weibull distribution appears to have the smallest KS test statistic value, 0.0800. Therefore, the Weibull distribution is chosen as the most suitable theoretical distribution to describe the actual catastrophe loss data distribution.
Next is the fitting test between the catastrophe fatality data and the geometric distribution. Geometric distribution parameter estimation is completed first. The estimation is carried out using the ML method. The parameter estimator of the geometric distribution obtained is = 0.0618. The fitting test used is the KS test. In short, with a significance level of 0.01, the KS test statistic value obtained, 0.1325, is less than the critical value, 0.1407. It indicates that the geometric distribution is suitable for describing the catastrophe fatality data distribution.

Parameter Estimation for the Annual Number of Catastrophe Data and Its Forecasting
The annual number of catastrophes modeling is carried out similarly to the annual force of interest modeling with ARIMA. In short, with a significance level of 0.01 for each statistical test performed, the ARIMA model (3,1,9) satisfies all assumptions in the ARIMA modeling. The model also has an MAPE in the 10% to 20% interval, namely 14.7897%. Therefore, the ARIMA model (3,1,9), whose parameter estimates are presented in Table 6, is practically good for forecasting the annual catastrophe intensity in the next few years [48]. The catastrophe intensity forecast results in year , = 1, 2, 3, 4, 5 are presented in Table 7.  15.1417 Table 7 shows that the catastrophe intensity at the year , = 1, 2, 3, 4, 5 tends to increase. We must be aware of this.

Determination of the Continuous Distribution to Approximate (•)
To determine which continuous distribution is used to approximate (•) , the skewness of and the kurtosis of are determined first. Previously, is Weibull distributed with parameters = 0.7253 and = 1.8058 . With this information, the skewness of can be determined, which is 3.2973. It is constant for all . Then, to determine the kurtosis of , information on the catastrophe intensity in year , = 1, 2, 3, 4, 5 is required. Then, using the information obtained in Table 7, the kurtosis of for = 1, 2, 3, 4, 5 is presented in Table 8.

Pricing of MTCB
Suppose that  Table 2. (g) The annual catastrophe intensity can be seen in Table 7.
Applying the continuous distribution approximation method and the Nuel recursive method in Equations (14) and (16), the price of zero-coupon MTCB and coupon-paying MTCB is 0.8759 USD and 0.9308 USD, respectively.

Analysis of the Effect of Catastrophe Intensity on MTCB Prices
Suppose that for = 1, 2, 3, = ∈ [14,20]. With the same assumptions as points (a) to (f) in Section 4.6, the visualization of the MTCB price for each is presented in Figure 2.   Figure 3. Figure 3 shows that the MTCB term and the price MTCB, both zero-coupon MTCB and coupon-paying MTCB, have a negative relationship. It means that the longer the MTCB term, the cheaper the MTCB price.  Figure 4. Figure 4 shows that both threshold values of MTCB and its price, both zero-coupon MTCB and coupon-paying MTCB, have a positive relationship. It means that the higher both threshold values of MTCB, the higher its price. It can also be seen that the larger both threshold values of MTCB, the brighter the green and purple colors.

The Increasing Catastrophe Intensity
The findings of this study are in line with Hoeppe [3], who stated that the catastrophe intensity worldwide, including in the United States, due to changes in weather and climate will tend to continue to increase in the future. From the model simulation results, the forecast for the annual catastrophe storm intensity in the United States carried out in Section 4.4 will tend to increase in the next five years (see Table 7). According to Zhang et al. [49], this is partly due to rising sea temperatures, particularly along the east coast of the United States, due to climate change. Warmer sea temperatures on the east coast of the United States can sustain hurricanes and increase the frequency of storm events.

Considerations for Investors in Purchasing MTCB
The increasing catastrophe intensity should be considered by investors when purchasing MTCB so that the probability that the coupon and principal are lost is low. The following is further information regarding the relationship between catastrophe intensity and MTCB prices that investors can consider when purchasing MTCB. Based on the analysis of the effect of catastrophe intensity on MTCB prices carried out in Section 4.7.1, the higher the catastrophe intensity, the lower the MTCB price. It is logically appropriate that if the catastrophe intensity is high, then the probability that MTCB claims occur so that investors lose their coupons and principal is also high. It causes investor demand for MTCB in the capital market to decrease so that the price will also decrease. Therefore, in purchasing MTCB, investors should pay attention to the catastrophe intensity so that it is not too high to reduce the probability of losing the coupon and principal.
In addition to considering the catastrophe intensity, investors must also consider the MTCB term when purchasing it. It is also so that the probability of losing the coupon and principal is low. The following is further information regarding the relationship between the MTCB term and its price that investors can consider when purchasing MTCB. Based on the analysis of the effect of the MTCB term on its prices carried out in Section 4.7.2, the longer the MTCB term, the lower its price. It makes sense that if the MTCB term is long, then the probability that MTCB claims occur is high. It is in line with the high probability that the investor will lose the coupon and principal, which is also high. It causes investor demand for MTCB with a long term in the capital market to decrease so that the price will also decrease. Therefore, in purchasing the MTCB, investors should pay attention to its term so that it is not too long to minimize the probability of the coupon and principal being lost.
Finally, when purchasing MTCB, investors should consider both threshold values. The probability of receiving the total coupon and principal is undoubtedly the goal. The following is further information regarding the relationship between threshold values and MTCB prices that investors can consider when purchasing MTCB. Based on the analysis of the influence of both MTCB threshold values on its price carried out in Section 4.7.3, the larger both MTCB threshold values are, the higher its price. It follows intuition, where if both MTCB threshold values are rising, then the probability that MTCB claims will occur is low. It is proportional to the loss of the investor's coupon and principal probability, which is also low. It certainly causes investor demand for MTCB with high threshold values in the capital market to increase its price. Therefore, in choosing the MTCB to purchase, investors should pay attention to both threshold values so that they are not too small to avoid the high probability of losing coupons and principal.

Conclusions
This study develops a zero-coupon and coupon-paying MTCB pricing model involving the increasing catastrophe intensity factor. The involvement of these factors is essential considering that the catastrophe intensity in countries worldwide tends to increase every year. The increasing catastrophe intensity is included in the modeling aggregate of losses and fatalities via NCPP. There is no closed-form solution to the MTCB pricing model because the CDF values of and are difficult to compute analytically in most cases. Therefore, we propose using numerical methods, the continuous distribution approximation method, and the Nuel recursive method to facilitate its calculation. With these numerical methods, the model solution can be determined simpler and more efficient.
The model simulation results show that the catastrophe intensity and other factors such as the term and both threshold values of MTCB affect the MTCB price. The catastrophe intensity and the MTCB term have a negative relationship with the MTCB price, where the greater the catastrophe intensity and the longer the MTCB term, the lower the MTCB price. The reason is that the greater the catastrophe intensity and the longer the MTCB term, the greater the probability that the claim will occur. It causes investor interest to purchase MTCB to decrease so that the price also decreases. Meanwhile, both threshold values of MTCB have a positive relationship with the price, where the larger both threshold values of MTCB, the higher the price. The reason is that the larger both threshold values of the MTCB, the smaller the probability that the claim will occur. It causes investor interest to purchase it to increase so that the price of MTCB becomes expensive.
This study can help practitioners/professionals, including sponsors, SPV, and investors. For sponsors and SPV, this study can help manage their catastrophe risk through MTCB (particularly in the pricing section), which can interest investors in the future. For investors, the results of model simulations related to the effect of catastrophe intensity, term, and threshold value on MTCB price can contribute to purchasing MTCB with the risk of claim events they want. For academics, as a further suggestion study, the factor of the correlation rate between aggregate losses and fatalities may affect the MTCB price. The involvement of this factor in MTCB price modeling can be used as an opportunity to develop the MTCB model in this study in the future.