#### Abstract

Catastrophe events are attracting increased attention because of their devastating consequences. Aimed at the nonlinear dependency and tail characteristics of different triggered indexes of multiple-event catastrophe bonds, this paper applies Copula function and the extreme value theory to multiple-event catastrophe bond pricing. At the same time, floating coupon and principal payoff structures are adopted instead of fixed coupon and principal payoff structures, to reduce moral hazard and improve bond attractiveness. Furthermore, we develop a CIR-Copula-POT bond pricing model with CIR stochastic rate and estimate flood multiple-event triggered catastrophe bond price using Monte Carlo simulation method. Finally, we implement the sensitivity analysis to show how catastrophe intensity, maturity date, and the dependence affect the prices of catastrophe bonds.

#### 1. Introduction

Different kinds of natural disasters occur frequently in the world over the past several years. These low-frequency and high-losses catastrophic events have a serious influence on peoples life and the stability society. Traditionally, when the catastrophic events occur, the national finance and social aid would be used to compensate the catastrophe losses. When faced with the natural disaster losses, the insurance companies themselves cannot satisfy the demand of catastrophe risk due to the large financial pressure and the restrictions of business ability. In recent years, there appear some kinds of important insurance-linked securities (ILSs) in the international catastrophe insurance market. And catastrophe (CAT) bonds are one of the most prominent ILSs, which transfer the consequence of CAT financial risks from issuers to investors. CAT bonds not only improve the risk bearing capacity of insurance companies, but also bring more investment choices to capital market. Reasonable pricing is the critical point in the CAT bonds issuing and trading.

Recently, some research efforts have been devoted to the catastrophe bonds pricing. Lee and Yu [1] developed a contingent-claim model to price CAT bonds with consideration of default risk, basis risk, and moral hazard. Egami and Young [2] presented a method for pricing structured CAT bonds based on utility indifference pricing. Härdle and Cabrera [3] calibrated the government parameter index-triggered CAT bonds with Mexico earthquake data; they proved that mixing reinsurance with CAT bond can reduce exposure risk and default risk. Z. G. Ma and C. Q. Ma [4] proposed a mixed approximation method to find the numerical solution for the price of catastrophe risk bonds. Nowak and Romaniuk [5] implemented the CAT bond pricing model described by the two-factor Vasicek model. Moreover, they proposed an automated approach for decision-making in fuzzy environment with relevant examples presenting this method.

While the catastrophe risks have obvious thick tail features, it is more reasonable to use extreme value theory (EVT) to characterize the tail characteristic of catastrophic losses distribution. Zimbids et al. [6] studied the Greece earthquakes data using advanced techniques from the extreme value theory. Moreover, they evaluated the CAT bond price using Monte Carlo simulation techniques and stochastic iterative equations. Shao et al. [7] applied equilibrium pricing theory and EVT to construct a multiple-variable CAT bond for California earthquakes. Ma et al. [8] employed a doubly stochastic Poisson and Peak over Threshold (POT) to price zero-coupon catastrophe bonds.

Since one-triggering-event CAT bond is difficult to meet the diverse needs of investors, the multiple-event CAT bond starts to rise because of its advantages. Recently, several studies have mainly focused on multiple-event CAT bonds. For example, Woo [9] addressed multiple-event risk securitization as a good way of transferring terrorism risk to capital markets. Reshetar [10] developed a framework for pricing of a multiple-event coupon paying CAT bond. It should be noted that they do not take fully the fat tail features of catastrophic risks into account. Considering that, this paper tries to combine Copula function and EVT to model a multiple-event CAT bond. The details are evaluating the marginal distributions of flood catastrophe losses and deaths via EVT. Furthermore, this paper employs the Copula function to model a joint distribution of losses and deaths. To make the price have a more applied value, we describe the spot interest rate by the CIR stochastic interest model.

The remainder of the paper is organized as follows: Section 2 briefly describes the framework of CAT claim model and stochastic interest model. Section 3 presents an empirical analysis. Section 4 is devoted to Monte Carlo simulation and sensitivity analysis. Finally, Section 5 offers conclusions for this research.

#### 2. Valuation Framework

##### 2.1. Pricing Model for the CAT Bond

Throughout this paper, we use the following assumptions:(i) is a Poisson process with the intensity ;(ii) and denote the economy losses and the number of deaths, respectively;(iii) and are two sequences of independent and identically distribution random variables;(iv), and are mutually independent.

We consider a coupon paying CAT bond, namely, paying a certain percentage coupons to investors at the end of the year and returning a certain percentage of principal at maturity date. In this case of catastrophe event, both coupon and principal are at risk. We choose the catastrophe losses and death tolls as trigger indicators. When one of the indicators is triggered, the current and future coupons are paid in proportion to cumulative catastrophe losses. And the principal is also paid in proportion to the cumulative losses only when both indicators are triggered simultaneously. The structures of payoff are given by the following:

the coupon paying framework:

the principal paying framework:where and stand for the attachment point for property losses and deaths, respectively, is coupon value, is bond maturity date, and denotes the percentage of principal paying under different cases. Moreover, in order to avoid the denominator meaninglessness (the number of deaths equal to zero), we replace with in (1).

Compared with the previous researches about multiple-event CAT bond, we adopt floating coupon and principal payment structures to replace the fixed coupon and principal payment structures. This will prevent insurance companies from increasing catastrophe losses on purpose after catastrophe events, resulting in moral hazard.

Let denote the price of CAT bond; then can be calculated by the following formula:where is the price at the time of a risk-free interest rate.

##### 2.2. Interest Rate Model

Cox, Ingersoll, and Ross (CIR) [11] model can not only describe the mean-reverting characteristic of interest rate, but also guarantee nonnegative interest rate, which is just the characteristic in real interest rate markets. Therefore, in this paper, the spot interest rate is assumed to follow CIR model. The interest rate model can be expressed as where is the speed of mean-reverting, is an mean of interest rate in the long run, is the volatility of the interest rate, and is a standard Brownian process. In the CIR model, the risk-neutral pricing of zero-coupon bond is given by following equation: where

##### 2.3. Extreme Value Theory and Modeling

Extreme value theory provides two methods to portray the extreme value behavior of observations, namely, Block Maxima Method (BMM) model and Peak Over Threshold model. However, BMM model is only interested in the behavior of the sample maximum, which could cause vast valid data missing. To take advantage of data information, the POT model (see [12]), which considers all observations exceeding a certain threshold value, shall be introduced to discuss the data behavior.

Suppose a random variable with the distribution function and a high threshold ; can be viewed as the statistical extremes; then the excess distribution function of iswhich implies that

According to the theorem of PBdH (1975), as is large enough, the excess distribution in (10) can be approximated by the generalized Pareto distribution (GPD); that is, where and are, respectively, the shape parameter and the scale parameter.

Then substituting (11) into (10), we obtain

#### 3. Empirical Analysis

##### 3.1. Description of Data

Our data consists of flood events that are recorded in Global Archive of Large Flood Event, provided by Dartmouth College since 1985. And our study mainly considers the losses value exceeding 100 thousand dollars. Thus, the total of 827 pairs of observations for losses and deaths are picked out. Directly analyzing the data, we will find that the fitting performance of data is not very well. Aimed at improving fitting accuracy, the data is adjusted to logarithm method to eliminate the magnitude difference.

Before applying EVT, heavy-tailed characteristic of flood data should be discussed. In general, there are two methods to judge the heavy-tailed behavior: numerical method and the exponential quantile-quantile (Q-Q) plot. Here the two methods shall be used for analysis. Figure 1 obviously shows the tail of exponential QQ plots appears in a convex shape. Besides, the kurtosis values are bigger than three. Therefore, the graphs strongly suggest that the hypotheses that the data follows GPD is acceptable.

**(a)**

**(b)**

##### 3.2. POT Model for Marginal Distribution

In POT model, if the losses after logarithm excess threshold , the exceeding part would be described by a generalized Pareto distribution. For the other part, the empirical function is suitable. In conclusion, the fitting distribution of the economic losses is given by

Since the logarithm of death number cannot be taken when its value is zero, we should take two steps to fit as follows. Firstly, for the deaths number exceeding one, we can also use the above-mentioned methods to give the following condition distribution of catastrophe death: Secondly, fit the data containing zero with (14). Then the unconditional distribution is defined as where can be evaluated with empirical distribution; namely, , where denotes the frequency of the number of deaths equaling zero.

##### 3.3. Threshold Selection, Parameter Estimates, and Testing

One of the main challenges in POT model is the selection of an optimal threshold for fitting the model. In practical application, mean excess function plot is usually used to set the threshold. Generally, the threshold is valid if the mean excess plot becomes roughly linear which starts from certain threshold level. For more details, see Embrechts et al. [13].

Figures 2 and 3 display the mean excess plots of flood losses and deaths, respectively. From Figure 2, we can see that the plot starts straight from the value 20.5 to the value 25 and seems to decline roughly after the value 25. Based on the mean excess plot, we might choose a threshold of 20.5. Furthermore, we choose a series of threshold values and use the maximum likelihood estimation getting a series of parameters to check the validity of the threshold. If the parameter around the threshold has stability property, in this sense the selected threshold is satisfactory. Then Figure 4 indicates that the shape parameter and scale parameter are more stable around the threshold 20.5. Similar conclusions can be drawn by considering the deaths data (Figure 5); the threshold for deaths is 5.8.

After setting the threshold values, the maximum likelihood method will be used to estimate the other parameters. Estimates are given in Table 1.

We use the diagnostic plots to check whether the GPD fits the data. Figures 6 and 7 are the diagnostic plots for losses and deaths, respectively. Almost all the dots of probability plots are on the same line. And the density plots fit well with the histograms. Consequently, it can be safely concluded that the chosen model gives a good fit for the loss and death data.

##### 3.4. Selection of Copula Function and Parameters Estimation

Copula function is called joint function, which connects joint distribution function with marginal distribution function. It usually studies the nonlinear relationship among variables. And Copula function mainly includes Elliptic Copula family and Archimedean Copula family. The later was used more often because the Elliptic Copula cannot describe the asymmetric relation of variables. Therefore, in this paper, three common Archimedean Copula functions, namely, Gumbel, Clayton, and Frank Copula, will be adopted to undertake related studies.

First of all, we plug the estimates of parameters shown in Table 1 into (13) and (15). Then, via probability integral transform, the distribution series whose values lie in are obtained. Furthermore, it follows from Box-test that the series are independent. Next, let the series be as the observations of Copula, and use maximum likelihood method to estimate the Copula parameters. As different Copula functions reflect different dependent patterns, it is very important to choose the suitable Copula function. To attain this goal, we adopt both the Kolmogorov-Smirnov (K-S) test and Q-Q plot to test for the appropriateness of the Copula model selection. As shown in Table 2, only Clayton Copula passes the test. And from Figure 8, we can observe that the Q-Q plot of Clayton Copula matches the straight line better than the Q-Q plots of Gumbel Copula and Frank Copula.

**(a)**

**(b)**

**(c)**

#### 4. Numerical Example and Stimulation

##### 4.1. Monte Carlo Simulation of CAT Bond

In this section, we will concentrate on valuing the price of CAT bond. Before the estimation, some related parameter values need to be set. According to the relative data of the real insurance market, the basic parameter values are given in Table 3.

Recalling that the explicit solutions of (1) and (2) are difficult to compute, we estimate the price of CAT bond by using Monte Carlo simulation techniques. The steps of simulation can be explicitly outlined as follows:

Simulate random numbers that follow Poisson distribution with intensity . And these values denote the times of flood in the first year, second year, , and the th year, respectively.

Generate pairs of random numbers from Clayton Copula function.

From marginal distribution, calculate the observations about losses and deaths. They are denoted by ; that is, , .

Substitute the results of procedure (3) into (1) and (2), we can get the values, and into (4), one simulation price follows.

Repeat the procedures (1)–(4) for times, and calculate the mean value.

To determine the price, we implement Monte Carlo simulations and finally obtain the price of CAT bond as 57.4417 USD.

##### 4.2. Sensitivity Analysis for the Catastrophe Bond

###### 4.2.1. The Effects of Intensity and Maturity Time

The bond prices decrease as the intensity gets stronger; the discount values of coupons also show the inverse relationship with intensity from Table 4. The main reason is that the stronger intensity means the bigger probability that catastrophe event would cause losses and deaths above the attachment points. The result may lead the bond prices and coupon discount values to decrease.

As can be seen from Table 5, we can see that the bond prices tend to decrease with the maturity time, while the discount coupons increase with the maturity time. A possible explanation of this result may be related to the same reason with intensity effect. The effect of maturity time on bond price is totally in line with traditional bond pricing theory. Furthermore, investors can attain more coupons with the extension of maturity time; this will raise the prices of bond. Thus, the decrease of bond prices can be attributed to the results of combination of a default increase probability and an increase yielded by coupons.

As shown in Tables 4 and 5, due to the stronger intensity and the longer period of time, they both lead to a bigger probability of bonds default and, therefore, a lower price. While the two cases are different. The discount coupons increase more with the maturity time extension. However, the results of arrival intensity are totally opposite. The increase coupons can offset a part of the risk caused by the triggering bond. Thus, the impact of arrival intensity on the price of catastrophe bonds is greater than the extension of maturity time. In other words, the effect of intensity on the price of CAT bond is more pronounced than the effect of maturity time.

###### 4.2.2. The Effect of Dependence

As we know, the relationship between Kendall rank correlation coefficient and Clayton Copula parameter is . The value of increases with . In other words, the stronger positive dependence between losses and deaths implies the greater probability that the two indications are triggered at the same time. The relationship between and the price of bond is shown in Table 6. The bond prices decrease as the increases. Therefore, when designing multiple-event CAT bonds, the less relevant indicators should be given priority; then the design goals can be achieved.

#### 5. Conclusions

Catastrophe bonds triggered by multiple-event have a lower catastrophic risk than the one triggering event, while they have higher yields and more market potential than ordinary bonds. This paper designs a multiple-event triggering pricing model and combines the Copula function and the POT method in extreme value theory to study the pricing of catastrophe bonds. The pricing model not only retains the advantages of the previous multiple-event trigger model, but also has the characteristics of low risk and high return. Some improvements have been made, floating coupon and principal payment structures as replacements for the original fixed coupon and principal payment structures, so as to reduce moral hazard and improve bond attractiveness. In addition, we implement Monte Carlo simulation to price CAT bonds using Global Archive of Large Flood Event, provided by Dartmouth College since 1985. Finally, the sensitivity analysis about the parameters of pricing model is also conducted. The empirical studies reveal that the price has the inverse relationship with arrival intensity, maturity time, and Copula dependence coefficient. The effect of intensity on the price of CAT bond is more pronounced than the effect of maturity time.

In recent years, the frequency and loss affected by natural disasters in the world have been constantly expanding, which poses a more severe challenge to the traditional catastrophe insurance market. Faced with low-loss frequency and high-loss severity catastrophe risk, the capacity of insurance companies for catastrophe is very limited, making catastrophe risk not yet effectively dispersed in the insurance market. Catastrophe bonds and other derivatives emerged, connecting the insurance market with the capital market and well transferring the catastrophe risk to the capital market. The model put forward in this paper enriches existing research on catastrophe bond pricing, especially in multiple-event catastrophe bond pricing. Furthermore, our research not only provides theoretical guidance for insurers to price the multiple-event catastrophe bond, but also provides a low-risk investment products to investors, enriching their investment portfolios.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.