# Valuing Catastrophe Bonds Involving Credit Risks

Authors: Jian Liu; Jihong Xiao; Lizhao Yan; Fenghua Wen
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

Catastrophe bonds are the most important products in catastrophe risk securitization market. For the operating mechanism, CAT bonds may have a credit risk, so in this paper we consider the influence of the credit risk on CAT bonds pricing that is different from the other literature. We employ the Jarrow and Turnbull method to model the credit risks and get access to the general pricing formula using the Extreme Value Theory. Furthermore, we present an empirical pricing study of the Property Claim Services data, where the parameters in the loss function distribution are estimated by the MLE method and the default probabilities are deduced by the US financial market data. Then we get the catastrophe bonds value by the Monte Carlo method.


---

## Body

## 1. Introduction

The securitization of catastrophe risk springing up in the early 1990s has created a direct link between the insurance industry and capital market. There are a variety of catastrophe risk securitization instruments, such as options, swaps, and bonds. Catastrophe bonds (CAT bonds), the largest issued and most successful instrument among them, have topped $6 billion in 2013 and are set to be the highest since 2007, according to Artemis, a deal tracker [1]. CAT bonds are usually insurance-linked and meant to raise money in case of a catastrophic event such as a hurricane or earthquake, so that insurance companies can hedge their exposure by transferring catastrophe risk to a wide pool of willing investors. Compared with other fixed income asset classes, CAT bonds offer high yield as the total returns were almost 9.5% in 2013, according to a Swiss Re index [1]. Furthermore, CAT bonds are not closely linked with the stock market or economic conditions, so they offer significant attractions to investors.The development of CAT bonds market depends on the reasonable prices, so the scientific pricing is the key problem to the field of CAT bonds research. As a kind of catastrophe risk securitization product,, the value of CAT bonds results from the probability of the catastrophe risk and the loss in the catastrophe, for CAT bonds have the dual properties of bonds and options. Moreover, for the highly skewed property of the catastrophe risk distribution, valuing CAT bonds has become very complicated.The main pricing models of CAT bonds, including Kreps model [2], LFC model [3, 4], Christofides model [5], and Wang two-factor model [6], are all using the quantitative methods to estimate essential factors of the price and then price the CAT bonds with them. Zimbidis et al. [7] use Extreme Value Theory to get the numerical results of CAT bonds prices under stochastic interest rates in an incomplete market framework. Z.-G. Ma and C.-Q. Ma [8] derive a bond pricing formula under stochastic interest rates with the losses following a compound nonhomogeneous Poisson process and find the numerical solution for the price of catastrophe risk bonds. Li et al. [9] study a representative agent-pricing model of the multievent CAT bonds by the data of catastrophic insured property losses of typhoon in China. Based on the idea of layered pricing, Xiao and Meng [10] use the extreme value model to discuss the pricing of the excess-of-loss reinsurance premium with different attachment points. Nowak and Romaniuk [11] use Monte Carlo simulation method to price the CAT bonds with different payoff functions. However, most of the literature in CAT bonds pricing research does not take the credit risk influence into account. Actually, the credit risk has the probability of existence for the operating mechanism of CAT bonds. So it is necessary to fully consider the credit risk in valuing research, which can improve the pricing validity.This paper presents a pricing model of CAT bonds with credit risks and conducts an empirical analysis of the US catastrophe market data using the Extreme Value Theory, where the parameters in the loss function distribution are estimated by the MLE method and the default probabilities are deduced by the US financial market data. Furthermore, based on the theoretical pricing formula, we get the CAT bonds value by the Monte Carlo method.


## 2. Operating Mechanism of CAT Bonds

There are four main participants in CAT bonds: sponsor, special purpose vehicle (SPV), investor, and trustee. The sponsor usually is an insurance company or a reinsurance company which signs a reinsurance contract with SPV and pays the reinsurance premium to SPV. SPV connects the insurance market and the capital market and converts the reinsurance premium into CAT bonds which are issued to investors. Then the reinsurance premium and the funds coming from CAT bonds investors are deposited in the trust account of the trustee and they are invested in short-term securities with low risk. When a catastrophic event takes place and its loss is higher than the prespecified trigger, SPV will provide compensation, which is from the loss of CAT bonds investors, to the sponsor according to the contract. In this case, investors only receive a part of the principal and interest. If the catastrophic risk event does not occur or trigger during the term of the CAT bond, bonds investors will receive their principal plus a high-yield compensation for the catastrophic risk exposure [12, 13].It follows that SPV is the key to the CAT bonds issue and its operation is directly related to the expected effect of CAT bonds. SPV is the intermediary in the securitization of catastrophe risk and has a significant role in the risk diversification and financing funds. SPV converts the constant return of the high security assets into float return using swaps and forwards via the trustees, who manage the funds, in order to pay interests to CAT bonds investors. However, the final appreciation of the capital depends on the credit rating of the counterparty. If the credit rating of the counterparty is low, the funds may shrink. In this case, investors of CAT bonds will not receive the high-yield predetermined compensation. It means that the CAT bonds have the credit risk, so in this paper we consider the influence of the credit risk to CAT bonds pricing.


## 3. Valuation Framework

### 3.1. Pricing Expression of CAT Bonds

Suppose the financial market is frictionless and arbitrage-free. The uncertainty in the market is characterized by the probability space(
Ω
,
ℱ
,
{
ℱ
t
}
t
≥
0
,
P
), where Ω is the state space, ℱ is the σ-algebra of measurable events, and P is the market probability measure. The maturity date of a CAT bond is T that is divided into n period with continuous trading interval Δ
t. An increasing filtration ℱ
t
⊂
ℱ, t
∈
[
0
,
T
]. The investor receives the predetermined interest at the end of each period and the interest of the current period and the whole principal in time T. The amount of the interest or the principal depends on whether the loss of the catastrophic event exceeds the trigger level. In this paper, when considering the credit risk of the CAT bond, the investment income also depends on the probability of default at each period. Suppose the probability of default at period i is λ
i and the recovery rate is a constant θ. Let H
i denote the payment amount of the CAT bond at period i which can be expressed as follows.When1
≤
i
≤
n
-
1,
(1)
H
i
=
{
C
i
,
I
i
≤
M
,
if
not
default
,
θ
C
i
,
I
i
≤
M
,
if
default
,
k
C
i
,
I
i
>
M
.Wheni
=
n,
(2)
H
n
=
{
C
n
+
F
,
I
n
≤
M
,
if
not
default
,
θ
(
C
n
+
F
)
,
I
n
≤
M
,
if
default
,
k
(
C
n
+
F
)
,
I
n
>
M
,
where F is the principal of the CAT bond, C
i is the interest of period i
,
k is the obtained ratio of interest and the principal when the CAT bond triggers, I
i is the loss amount of the catastrophic event at period i, and M is the amount of the loss trigger.Then, the price of the CAT bond at time 0 is(3)
V
=
E
[
∑
i
=
1
n
e
-
r
·
i
Δ
t
H
(
λ
i
,
I
i
)
]
,
where
(4)
H
(
λ
i
,
I
i
)
=
{
k
C
i
·
P
(
I
i
>
M
)
+
[
λ
i
θ
C
i
+
(
1
-
λ
i
)
C
i
]
·
P
(
I
i
≤
M
)
,
1
≤
i
≤
n
-
1
,
k
(
C
n
+
F
)
·
P
(
I
n
>
M
)
+
[
λ
n
θ
(
C
n
+
F
)
kk
+
(
1
-
λ
n
)
(
C
n
+
F
)
]
·
P
(
I
n
≤
M
)
,
i
=
n
.
Let
(5)
V
(
l
)
=
1
l
∑
j
=
1
l
∑
i
=
1
n
e
-
r
·
i
Δ
t
H
(
λ
i
(
j
)
,
I
i
(
j
)
)
.
Assuming expression (3) exists, we can approximate the price of the CAT bond as V
=
li
m
l
→
∞
V
(
l
) where the number of simulation paths is l. So the price of the CAT bond at time 0 can be calculated by Monte Carlo method. We will give more detail of the method in Section 4. Because P
(
I
i
>
M
)
=
1
-
P
(
I
i
≤
M
), so the key for the CAT bond pricing is to obtain the probability of default λ
i and the probability distribution of loss process I
i.


### 3.2. Credit Risk

This paper uses the methodology of Jarrow and Turnbull [14] to model the credit risk of the CAT bond. The credit risk can be described by the probability of default and the recovery value which are deduced by the arbitrage-free valuation techniques. This methodology is simple and fits the existing term structure of interest rate well. As the absolute priority rule is often violated and lots of other factors affect the payoff of bonds, such as the percentage of managerial ownership [14], modeling the actual payoff in default is a complicated problem. Therefore, we take the recovery rate in the event of default when the CAT bond does not trigger as an exogenously given constant. The recovery rate is denoted by θ and it is assumed to be the same for all bonds in a given credit ranking. Then, we can use the risk-free interest rate to discount the cash flow of the CAT bond but not use the discount rate involving risk premium [15]. In fact, the probability of default has reflected the credit risk of the CAT bond.Suppose the probability of default in time interval[
i
-
1
,
i
] is denoted by λ
i, and let λ
i
*
=
-
(
1
/
i
)
∑
t
=
1
i
ln
⁡
⁡
(
1
-
λ
t
); then we get ∏
t
=
1
i
(
1
-
λ
t
)
=
e
-
λ
i
*
·
i. If λ
i is small, then we can obtain
(6)
λ
i
*
≈
∑
t
=
1
i
λ
t
i
.
It means that λ
i
* is the average intensity of default.The interest rate of the risk bond in time interval[
i
-
1
,
i
] is denoted by r
i
* and the risk-free interest rate is r
i. For the no-arbitrage principle, parameters r
i, r
i
*, θ, λ
i
* meet the equation e
-
r
i
*
=
[
1
·
(
1
-
λ
i
*
)
+
θ
·
λ
i
*
]
e
-
r
i, and then we get
(7)
λ
i
*
=
1
-
e
r
i
-
r
i
*
1
-
θ
.
Consequently, we can obtain the average intensity of default λ
i
* by (6) and deduce the probability of default in each time interval [
i
-
1
,
i
] which is λ
i, i
=
1,2
,
…
,
n.


### 3.3. Distribution of Loss Function

The loss functionI
i is the maximum loss amount of the catastrophic event at period i and I
i
=
max
⁡
⁡
{
X
1
i
,
X
2
i
,
…
,
X
m
i
}, where m is the number of days at each period and X
1
i
,
X
2
i
,
…
,
X
m
i is the sequence of loss amount at each period. They are independent random variables with a common unknown distribution function. If there exist sequences of constants {
a
k
:
a
k
>
0
,
∀
k
∈
ℕ
}, {
b
k
}
k
∈
ℕ and a nondegenerate distribution function G
(
z
), such that
(8)
P
(
I
k
-
b
k
a
k
≤
z
)
⟶
G
(
z
)
,
as
k
⟶
∞
,
z
∈
ℝ
,
then by Fisher-Tipptt Gnedenko Theorem, G is a member of the generalized extreme value (GEV) family of distributions or von Mises Type Extreme Value distribution or the von Mises-Jenkinson type distribution [16] and its distribution function form can be denoted as
(9)
G
(
z
)
=
exp
⁡
⁡
{
-
[
1
+
ξ
(
z
-
μ
σ
)
]
-
1
/
ξ
}
.
It is defined on the set {
z
:
1
+
ξ
(
z
-
μ
)
/
σ
>
0
}, where the scale parameter satisfies σ
>
0, and the trail index satisfies -
∞
<
ξ
<
∞ and location parameter satisfies -
∞
<
μ
<
∞. If we got the estimated values of the three parameters, the distribution function of the maximum loss amount in a catastrophic event I
i can be gotten [17].There are different kinds of considerable estimation techniques for three parameters{
σ
,
ξ
,
μ
}, such as Generalized Method of Moments, the graphical techniques based on versions of probability plots, and maximum likelihood estimation (MLE) method. The MLE method is a classical estimation technique, and it is effective for the estimations of parameters {
σ
,
ξ
,
μ
} studied by many researchers [18–20]. We employ MLE method to estimate parameters. Suppose that the loss amount processes I
1
,
I
2
,
…
,
I
n are independent variables following GEV distribution; then the log-likelihood function for parameters μ
,
σ
,
ξ is
(10)
L
(
σ
,
ξ
,
μ
)
=
n
log
⁡
⁡
σ
-
(
1
+
1
ξ
)
∑
i
=
1
n
log
⁡
⁡
[
1
+
ξ
(
I
i
-
μ
σ
)
]
-
∑
i
=
1
n
[
1
+
ξ
(
I
i
-
μ
σ
)
]
-
1
/
ξ
,
where 1
+
ξ
(
(
I
i
-
μ
)
/
σ
)
>
0, as i
=
1
,
2
,
…
,
n.


## 3.1. Pricing Expression of CAT Bonds

Suppose the financial market is frictionless and arbitrage-free. The uncertainty in the market is characterized by the probability space(
Ω
,
ℱ
,
{
ℱ
t
}
t
≥
0
,
P
), where Ω is the state space, ℱ is the σ-algebra of measurable events, and P is the market probability measure. The maturity date of a CAT bond is T that is divided into n period with continuous trading interval Δ
t. An increasing filtration ℱ
t
⊂
ℱ, t
∈
[
0
,
T
]. The investor receives the predetermined interest at the end of each period and the interest of the current period and the whole principal in time T. The amount of the interest or the principal depends on whether the loss of the catastrophic event exceeds the trigger level. In this paper, when considering the credit risk of the CAT bond, the investment income also depends on the probability of default at each period. Suppose the probability of default at period i is λ
i and the recovery rate is a constant θ. Let H
i denote the payment amount of the CAT bond at period i which can be expressed as follows.When1
≤
i
≤
n
-
1,
(1)
H
i
=
{
C
i
,
I
i
≤
M
,
if
not
default
,
θ
C
i
,
I
i
≤
M
,
if
default
,
k
C
i
,
I
i
>
M
.Wheni
=
n,
(2)
H
n
=
{
C
n
+
F
,
I
n
≤
M
,
if
not
default
,
θ
(
C
n
+
F
)
,
I
n
≤
M
,
if
default
,
k
(
C
n
+
F
)
,
I
n
>
M
,
where F is the principal of the CAT bond, C
i is the interest of period i
,
k is the obtained ratio of interest and the principal when the CAT bond triggers, I
i is the loss amount of the catastrophic event at period i, and M is the amount of the loss trigger.Then, the price of the CAT bond at time 0 is(3)
V
=
E
[
∑
i
=
1
n
e
-
r
·
i
Δ
t
H
(
λ
i
,
I
i
)
]
,
where
(4)
H
(
λ
i
,
I
i
)
=
{
k
C
i
·
P
(
I
i
>
M
)
+
[
λ
i
θ
C
i
+
(
1
-
λ
i
)
C
i
]
·
P
(
I
i
≤
M
)
,
1
≤
i
≤
n
-
1
,
k
(
C
n
+
F
)
·
P
(
I
n
>
M
)
+
[
λ
n
θ
(
C
n
+
F
)
kk
+
(
1
-
λ
n
)
(
C
n
+
F
)
]
·
P
(
I
n
≤
M
)
,
i
=
n
.
Let
(5)
V
(
l
)
=
1
l
∑
j
=
1
l
∑
i
=
1
n
e
-
r
·
i
Δ
t
H
(
λ
i
(
j
)
,
I
i
(
j
)
)
.
Assuming expression (3) exists, we can approximate the price of the CAT bond as V
=
li
m
l
→
∞
V
(
l
) where the number of simulation paths is l. So the price of the CAT bond at time 0 can be calculated by Monte Carlo method. We will give more detail of the method in Section 4. Because P
(
I
i
>
M
)
=
1
-
P
(
I
i
≤
M
), so the key for the CAT bond pricing is to obtain the probability of default λ
i and the probability distribution of loss process I
i.


## 3.2. Credit Risk

This paper uses the methodology of Jarrow and Turnbull [14] to model the credit risk of the CAT bond. The credit risk can be described by the probability of default and the recovery value which are deduced by the arbitrage-free valuation techniques. This methodology is simple and fits the existing term structure of interest rate well. As the absolute priority rule is often violated and lots of other factors affect the payoff of bonds, such as the percentage of managerial ownership [14], modeling the actual payoff in default is a complicated problem. Therefore, we take the recovery rate in the event of default when the CAT bond does not trigger as an exogenously given constant. The recovery rate is denoted by θ and it is assumed to be the same for all bonds in a given credit ranking. Then, we can use the risk-free interest rate to discount the cash flow of the CAT bond but not use the discount rate involving risk premium [15]. In fact, the probability of default has reflected the credit risk of the CAT bond.Suppose the probability of default in time interval[
i
-
1
,
i
] is denoted by λ
i, and let λ
i
*
=
-
(
1
/
i
)
∑
t
=
1
i
ln
⁡
⁡
(
1
-
λ
t
); then we get ∏
t
=
1
i
(
1
-
λ
t
)
=
e
-
λ
i
*
·
i. If λ
i is small, then we can obtain
(6)
λ
i
*
≈
∑
t
=
1
i
λ
t
i
.
It means that λ
i
* is the average intensity of default.The interest rate of the risk bond in time interval[
i
-
1
,
i
] is denoted by r
i
* and the risk-free interest rate is r
i. For the no-arbitrage principle, parameters r
i, r
i
*, θ, λ
i
* meet the equation e
-
r
i
*
=
[
1
·
(
1
-
λ
i
*
)
+
θ
·
λ
i
*
]
e
-
r
i, and then we get
(7)
λ
i
*
=
1
-
e
r
i
-
r
i
*
1
-
θ
.
Consequently, we can obtain the average intensity of default λ
i
* by (6) and deduce the probability of default in each time interval [
i
-
1
,
i
] which is λ
i, i
=
1,2
,
…
,
n.


## 3.3. Distribution of Loss Function

The loss functionI
i is the maximum loss amount of the catastrophic event at period i and I
i
=
max
⁡
⁡
{
X
1
i
,
X
2
i
,
…
,
X
m
i
}, where m is the number of days at each period and X
1
i
,
X
2
i
,
…
,
X
m
i is the sequence of loss amount at each period. They are independent random variables with a common unknown distribution function. If there exist sequences of constants {
a
k
:
a
k
>
0
,
∀
k
∈
ℕ
}, {
b
k
}
k
∈
ℕ and a nondegenerate distribution function G
(
z
), such that
(8)
P
(
I
k
-
b
k
a
k
≤
z
)
⟶
G
(
z
)
,
as
k
⟶
∞
,
z
∈
ℝ
,
then by Fisher-Tipptt Gnedenko Theorem, G is a member of the generalized extreme value (GEV) family of distributions or von Mises Type Extreme Value distribution or the von Mises-Jenkinson type distribution [16] and its distribution function form can be denoted as
(9)
G
(
z
)
=
exp
⁡
⁡
{
-
[
1
+
ξ
(
z
-
μ
σ
)
]
-
1
/
ξ
}
.
It is defined on the set {
z
:
1
+
ξ
(
z
-
μ
)
/
σ
>
0
}, where the scale parameter satisfies σ
>
0, and the trail index satisfies -
∞
<
ξ
<
∞ and location parameter satisfies -
∞
<
μ
<
∞. If we got the estimated values of the three parameters, the distribution function of the maximum loss amount in a catastrophic event I
i can be gotten [17].There are different kinds of considerable estimation techniques for three parameters{
σ
,
ξ
,
μ
}, such as Generalized Method of Moments, the graphical techniques based on versions of probability plots, and maximum likelihood estimation (MLE) method. The MLE method is a classical estimation technique, and it is effective for the estimations of parameters {
σ
,
ξ
,
μ
} studied by many researchers [18–20]. We employ MLE method to estimate parameters. Suppose that the loss amount processes I
1
,
I
2
,
…
,
I
n are independent variables following GEV distribution; then the log-likelihood function for parameters μ
,
σ
,
ξ is
(10)
L
(
σ
,
ξ
,
μ
)
=
n
log
⁡
⁡
σ
-
(
1
+
1
ξ
)
∑
i
=
1
n
log
⁡
⁡
[
1
+
ξ
(
I
i
-
μ
σ
)
]
-
∑
i
=
1
n
[
1
+
ξ
(
I
i
-
μ
σ
)
]
-
1
/
ξ
,
where 1
+
ξ
(
(
I
i
-
μ
)
/
σ
)
>
0, as i
=
1
,
2
,
…
,
n.


## 4. Numerical Analysis

### 4.1. Data and Parameter Estimations

The analysis is based on the catastrophe data from Property Claim Services (PSC). The data covers all losses resulting from natural catastrophic events in the USA over the period 1985–2011, with 3628 original data. We get the series of annual maximum magnitude of catastrophe losses in the USA. In order to facilitate comparison, the data is adjusted for inflation given the time value of the capital, using the Consumer Price Index (CPI) provided by the US Department of Labor. Figure1 shows the adjusted annual maximum magnitude of losses. As described in Table 1, the mean of catastrophe losses is 7.046 billion, the standard deviation is 10.697 billion, the skewness is 2.54, and the kurtosis is 7.19. It means that the data is right-skewed and tail-dispersed.Table 1
Descriptive statistic of catastrophe losses data (billion dollars).
Sample size (n)
Minimum
Maximum
Mean
Standard deviation
Skewness
Kurtosis
Statistic
Std. error
Statistic
Std. error
27
0.168
4.7337
7.046
10.697
2.539
0.448
7.193
0.872Figure 1
Scatter plot of the annual maximum magnitude of catastrophe losses in the USA (1985–2011).Assume that the data is independent random variables with GEV distribution function. By the extRemes program package in R software, we get the Maximum Likelihood estimations of parameters in distribution function (10) as follows:
(11)
(
μ
,
σ
,
ξ
)
=
(
16.22
,
20.25
,
1.02
)
.
The covariance matrix of the three parameters estimations is
(12)
cov
⁡
=
[
23.40
25.99
-
0.53
25.99
38.36
0.14
-
0.53
0.14
0.11
]
.
Figure 2 shows the various plots for assessing the accuracy of the GEV model, fitted to the annual maximum magnitude of catastrophe losses of US data. Each set of plotted points is near-line in the probability plot and the quantile plot, so the fitted model is valid. The respective estimated curve in the return level plot is close to a straight line and the corresponding density plot seems consistent with the histogram of the data. Therefore, the estimated GEV distribution fits the real data well and the GEV model, deduced by the MLE method, is acceptable and valid. Then, we use the estimations of the model parameters to calculate the probability of the catastrophe loss by (9).Figure 2
Diagnostic plots for GEV fit to the annual maximum magnitude of catastrophe losses in the USA.Now, we consider the credit risk. The credit ratings of most CAT bonds are BB, while a CAT bond rated AA was issued in 2006 for the first time. So we assume that the credit rating of the CAT bond is BB. Furthermore, we take the treasury rates in the US market as the risk-free interest rate and take the corporation bonds rates as the risk rates. Suppose the recovery rateθ
=
40
% when the CAT bond defaults. Then, according to the US market data in July, 2013, the treasury rates of 6 months, 2 years, 3 years, and 5 years are 0.06%, 0.39%, 0.74%, and 1.61%,respectively, and the corresponding rates of corporation bonds rated BB are 3.05%, 5.24%, 6.83%,and 7.35%, respectively. We get the treasury rates and corporation bonds rates from the risk-free interest rates and the risk interest rates, respectively, over the time period of 1–5 years by the linear interpolation method as Table 2 shows. Moreover, we get the default rates of risk bonds in each period by expressions (6) and (7), just as Table 3 shows.Table 2
Risk-free interest rates and risk interest rates.
Time
1
2
3
4
5
Risk-free
0.17%
0.39%
0.74%
1.18%
1.61%
Risk
3.78%
5.24%
6.83%
7.09%
7.35%Table 3
Default rates of risk bonds.
Time period
0-1(
λ
1
)
1-2(
λ
2
)
2-3(
λ
3
)
3-4(
λ
4
)
4-5(
λ
5
)
Default rate
5.92%
9.86%
13.77%
8.72%
8.23%


### 4.2. Pricing Analysis

Consider a CAT bond with the durationn
=
5 years, the face value F
=
100 dollars, and the coupon rate 9%. The amount of trigger loss is assumed to be the mean of the loss amount; that is to say, M
=
7.046 billion dollars. The obtained ratio of the interest and the principal when the CAT bond triggers is k
=
0.8. Then, we get the probability of trigger P
(
I
>
M
)
=
0.2404 by expression (9) substituting the above-mentioned parameters. Then by (9), we deduce
(13)
z
=
μ
+
[
σ
(
1
-
(
-
log
⁡
⁡
G
(
z
)
)
-
ξ
)
]
.
As the function G
(
z
) is the distribution probability, it meets 0
≤
G
(
z
)
≤
1. So by generating random numbers on the interval [
0
,
1
], we simulate the variable z employing the Monte Carlo method with the number of simulation paths l
=
10000. For expression (4), we get the pricing result of the CAT bond which is V
=
118.86 dollars by the Matlab software. In the same way, the price of fixed income instrument without credit risk is calculated as 139.27 dollars with the same interest rates. The former is much less than the latter that attributes to the risk of catastrophe and the credit risk. For a higher risk, the CAT bond gives a higher yield.


## 4.1. Data and Parameter Estimations

The analysis is based on the catastrophe data from Property Claim Services (PSC). The data covers all losses resulting from natural catastrophic events in the USA over the period 1985–2011, with 3628 original data. We get the series of annual maximum magnitude of catastrophe losses in the USA. In order to facilitate comparison, the data is adjusted for inflation given the time value of the capital, using the Consumer Price Index (CPI) provided by the US Department of Labor. Figure1 shows the adjusted annual maximum magnitude of losses. As described in Table 1, the mean of catastrophe losses is 7.046 billion, the standard deviation is 10.697 billion, the skewness is 2.54, and the kurtosis is 7.19. It means that the data is right-skewed and tail-dispersed.Table 1
Descriptive statistic of catastrophe losses data (billion dollars).
Sample size (n)
Minimum
Maximum
Mean
Standard deviation
Skewness
Kurtosis
Statistic
Std. error
Statistic
Std. error
27
0.168
4.7337
7.046
10.697
2.539
0.448
7.193
0.872Figure 1
Scatter plot of the annual maximum magnitude of catastrophe losses in the USA (1985–2011).Assume that the data is independent random variables with GEV distribution function. By the extRemes program package in R software, we get the Maximum Likelihood estimations of parameters in distribution function (10) as follows:
(11)
(
μ
,
σ
,
ξ
)
=
(
16.22
,
20.25
,
1.02
)
.
The covariance matrix of the three parameters estimations is
(12)
cov
⁡
=
[
23.40
25.99
-
0.53
25.99
38.36
0.14
-
0.53
0.14
0.11
]
.
Figure 2 shows the various plots for assessing the accuracy of the GEV model, fitted to the annual maximum magnitude of catastrophe losses of US data. Each set of plotted points is near-line in the probability plot and the quantile plot, so the fitted model is valid. The respective estimated curve in the return level plot is close to a straight line and the corresponding density plot seems consistent with the histogram of the data. Therefore, the estimated GEV distribution fits the real data well and the GEV model, deduced by the MLE method, is acceptable and valid. Then, we use the estimations of the model parameters to calculate the probability of the catastrophe loss by (9).Figure 2
Diagnostic plots for GEV fit to the annual maximum magnitude of catastrophe losses in the USA.Now, we consider the credit risk. The credit ratings of most CAT bonds are BB, while a CAT bond rated AA was issued in 2006 for the first time. So we assume that the credit rating of the CAT bond is BB. Furthermore, we take the treasury rates in the US market as the risk-free interest rate and take the corporation bonds rates as the risk rates. Suppose the recovery rateθ
=
40
% when the CAT bond defaults. Then, according to the US market data in July, 2013, the treasury rates of 6 months, 2 years, 3 years, and 5 years are 0.06%, 0.39%, 0.74%, and 1.61%,respectively, and the corresponding rates of corporation bonds rated BB are 3.05%, 5.24%, 6.83%,and 7.35%, respectively. We get the treasury rates and corporation bonds rates from the risk-free interest rates and the risk interest rates, respectively, over the time period of 1–5 years by the linear interpolation method as Table 2 shows. Moreover, we get the default rates of risk bonds in each period by expressions (6) and (7), just as Table 3 shows.Table 2
Risk-free interest rates and risk interest rates.
Time
1
2
3
4
5
Risk-free
0.17%
0.39%
0.74%
1.18%
1.61%
Risk
3.78%
5.24%
6.83%
7.09%
7.35%Table 3
Default rates of risk bonds.
Time period
0-1(
λ
1
)
1-2(
λ
2
)
2-3(
λ
3
)
3-4(
λ
4
)
4-5(
λ
5
)
Default rate
5.92%
9.86%
13.77%
8.72%
8.23%


## 4.2. Pricing Analysis

Consider a CAT bond with the durationn
=
5 years, the face value F
=
100 dollars, and the coupon rate 9%. The amount of trigger loss is assumed to be the mean of the loss amount; that is to say, M
=
7.046 billion dollars. The obtained ratio of the interest and the principal when the CAT bond triggers is k
=
0.8. Then, we get the probability of trigger P
(
I
>
M
)
=
0.2404 by expression (9) substituting the above-mentioned parameters. Then by (9), we deduce
(13)
z
=
μ
+
[
σ
(
1
-
(
-
log
⁡
⁡
G
(
z
)
)
-
ξ
)
]
.
As the function G
(
z
) is the distribution probability, it meets 0
≤
G
(
z
)
≤
1. So by generating random numbers on the interval [
0
,
1
], we simulate the variable z employing the Monte Carlo method with the number of simulation paths l
=
10000. For expression (4), we get the pricing result of the CAT bond which is V
=
118.86 dollars by the Matlab software. In the same way, the price of fixed income instrument without credit risk is calculated as 139.27 dollars with the same interest rates. The former is much less than the latter that attributes to the risk of catastrophe and the credit risk. For a higher risk, the CAT bond gives a higher yield.


## 5. Conclusions

As catastrophes are small probability and high loss events and present high positive correlation among individuals, the traditional insurance company and reinsurance company may not spread risk completely. However, the securitization of catastrophe risk brings an effective solution to transfer and spread the catastrophe risk. In order to develop the CAT bond market, it is necessary to make effective pricing to CAT bond which is the most mature instrument in catastrophe securitization at the present. To make the pricing result accurate, the credit risk in the CAT bond should not be ignored. This paper builds the model of the CAT bond involving the credit risk. This model has the characteristics of an analytic method and numerical method and has a good practical feasibility. We employ the Jarrow and Turnbull method to model the credit risk and get access to the general pricing formula using the Extreme Value Theory. Furthermore, we present an empirical pricing study of the Property Claim Services data where the parameters in the loss function distribution are estimated by the MLE method and the default probabilities are deduced by the US financial market data. Consequently, we get the catastrophe bonds value by the Monte Carlo method which is lower than the price of fixed income instruments.Further research is directed to extend the model in more complex situations, for example, given suitable stochastic process to describe the trigger amount of loss, and it will perfect the model which we have discussed.



---

*Source: 563086-2014-04-17.xml*