^{1}

^{1}

The main goal of this paper is to use the enlargement of filtration framework for pricing zero-coupon CAT bonds. For this purpose, we develop two models where the trigger event time is perfectly covered by an increasing sequence of stopping times with respect to a reference filtration. Hence, depending on the nature of these stopping times the trigger event time can be either accessible or totally inaccessible. When some of these stopping times are not predictable, the trigger event time is totally inaccessible, and very nice mathematical computations can be derived. When the stopping times are predictable, the trigger event time is accessible, and this case would be a meaningful choice for Model 1 from a practical point of view since features like seasonality are already captured by some quantities such as the stochastic intensity of the Poisson process. We compute the main tools for pricing the zero-coupon CAT bond and show that our constructions are more general than some existing models in the literature. We obtain some closed-form prices of zero-coupon CAT bonds in Model 2 so we give a numerical illustrative example for this latter.

Catastrophe (CAT) bonds are an alternative to traditional reinsurance intended to hedge against the risks incurred following a natural disaster such as earthquakes, pandemics, etc. The operation process is as follows. When an insurer or reinsurer wishes to cover a risk in a catastrophic geographical area, they are taken out by investors in the financial markets who perceive so many coupon exchanges. The interest of the issuer is to eliminate the risk of default and increase the available capacity. The operating principle is similar to that of conventional bonds. But in the case of the occurrence of triggering events, the investor loses all or part of the coupons or even the principal of the bond. In the absence of events after a given period, generally between 3 and 5 years, investors recover their initial stakes plus interests.

The first catastrophic bonds emerged in 1994, two years after Hurricane Andrew hit, which cost the US $20 billion. Insurers and reinsurers then realized that natural risks were going to be more and more expensive. CAT bonds were then created as an alternative and complementary to cover extreme risks. In a nutshell, the issuing of catastrophe (CAT) bonds has been and continuous to be essential for insurance companies. Recently, we are witnessing the occurrence of the Coronavirus disease COVID-19 pandemic that causes significant losses and the use of these products would be essential to cover against such losses.

The modeling of CAT bonds is somehow similar to that of credit risk. Indeed, most of the models interested in the construction of the investor’s trigger event time (by misuse of language, we may sometimes call it simply default time) τ are based on the structural approach used to price credit derivatives where τ is the first moment that the aggregate losses process L exceeds the threshold’s fixed value D (we simply call these models first hitting time models). Thus the investor loses all or part of his principal as soon as τ occurs before the maturity of the CAT bond. The advantage of this type of modeling is that the default time can coincide with the arrival times of catastrophic events. However, even if the aggregate losses approach the threshold (for example L t / D = 0.9999 , for all t ≥ 0 ) as long as the overshoot does not occur, the default mechanism cannot exist. This is not advantageous for the issuer of the CAT (for example a reinsurance company).

Several good examples using this approach include the model of Burnecki and Kukla [

In the first hitting time models, default occurs when the aggregate claims process exceeds a specified level (called threshold value). However, in that framework even if the trigger event time can coincide with the arrival time of catastrophe events this interesting property is not highlighted in most of the papers.

Another approach to model the investor’s trigger event mechanism is to evaluate all the losses until maturity, so the investor loses his principal as soon as the latter exceeds the value of the threshold fixed in advance, otherwise, he recovers its principal plus the coupons. But in this last setting, valuation and payment can be made by taking into account only what happens after maturity. While the default mechanism could even appear just a few months or years after the signing of the contract and for an investor it would be necessary to find out about his default as soon as possible than to wait until the end of the contract. Among others, it should be referred to Schmidt [

These two approaches are equivalent as soon as we are no longer interested in payment at hit and when the percentage of the principal lost by investor in case of the default mechanism is no more stochastic (this is generally the case in CAT bond modeling). Note also that in these two types of approaches, the calibration of the aggregate loss process is essential since the price of CAT bonds is affected by the frequency and also the severity of catastrophic losses.

In other terms, the CAT bond prices depend on the aggregate claims distribution. However in most of the cases, closed form of that distribution does not exist and one resorts to numerical approximations. For instance, [

To the best of our knowledge, only the paper of Jarrow [

In this paper, we introduced two models based on the enlargement of filtration theory for pricing zero-coupon CAT bonds. Our approaches will be focused on the information quantified in the CAT bonds. Indeed, we build two models that inherit some credit risk hybrid ones such as the ones developed in [

Our two models are general frameworks that may lead to a new perspective on modeling CAT bonds.

The paper is organized as follows. We first introduce notation and basic notions in Section 2. In Section 3, we present our two models in a general framework by computing the quantities of interest for the pricing of zero-coupon bonds, as the conditional expectations and the dual predictable projection. In the first model, we show that in the general case where some of the stopping times in the reference filtration are not predictable, our models cover some existing ones in the literature on CAT bond modeling. We also focus on a most meaningful case from the practical point of view, by which the stopping times in the reference filtration are predictable. In the last section, we establish a case study based on simulations of Model 2 where we deal with a particular aggregate loss process called the Shot-Noise process. We illustrate how the jumps in the (Shot-Noise) aggregate loss process induce the jumps in the CAT bonds pricing and compare our results with the ones obtained using a Compound Poisson process where the prices are always continuous.

In this section we recall, for the ease of the reader, some well-known results and definitions.

We consider a probability space ( Ω , G , ℙ ) and τ a positive random time defined on ( Ω , G ) . We introduce the right-continuous increasing default process A t = 1 l { τ ≤ t } associated with τ and we denote by A = ( A t ) t ≥ 0 the filtration (completed and right-continuous) generated by A. We recall that, for any process X, one has ∫ u t L s d A s = ∫ ] u , t ] L s d A s = L τ 1 l { u < τ ≤ t } .

Let ℍ a given filtration on Ω . There exists a unique ℍ -optional locally integrable variation process A o , ℍ , called the ℍ -dual optional projection of A, such that

E [ ∫ 0 ∞ L s d A s ] = E [ ∫ 0 ∞ L s d A s o , ℍ ]

for any bounded ℍ -optional process X such that E [ ∫ 0 ∞ | L s | d | A | s ] < ∞ .

There exists a unique ℍ -predictable locally integrable variation process A p , ℍ , called the dual predictable projection of A, such that

E [ ∫ 0 ∞ L s d A s ] = E [ ∫ 0 ∞ L s d A s p , ℍ ]

for any bounded ℍ -predictable process X such that E [ ∫ 0 ∞ | L s | d | A | s ] < ∞ . we shall sometimes call them ℍ -dual projections of τ. (see [

If τ is an ℍ -stopping time, the compensator of τ is by definition the unique ℍ -predictable increasing process J ℍ such that J 0 ℍ = 0 and A t − J t ℍ is an ℍ -martingale (see ( [

Lemma 2.1. For any ℍ -predictable bounded process H, the process

H τ 1 l { τ ≤ t } − ∫ 0 t ∧ τ H s d A s p , ℍ

isan ℍ -martingale.

Proof: This result follows from the fact that

H τ 1 l { τ ≤ t } = ∫ 0 t H s d A s = ( H ▪ A ) t ,

and, for H being ℍ -predictable, the ℍ -dual predictable projection of H ▪ A is H ▪ A p , ℍ (see [

Definition 2.2. A random time τ is said to avoid ℍ -stopping times, if for any finite ℍ -stopping time ξ , one has ℙ ( τ = ξ ) = 0 .

Definition 2.3. Let ϑ an ℍ -stopping.

● We say ϑ to be ℍ -predictable if there exists an increasing sequence of ℍ -stopping times ( ϑ i ) i ≥ 1 converging to ϑ such that ϑ i < ϑ on the set { ϑ i > 0 } ,for all i .If ϑ is ℍ -predictable, ( 1 l { ϑ ≤ t } , t ≥ 0 ) is a predictable process.

● We say that ϑ is accessible if 〚 ϑ 〛 ⊂ ∪ i 〚 ϑ i 〛 where ( ϑ i ) i ≥ 1 are ℍ -predi- ctablestopping times,with 〚 ϑ 〛 denotes the graph of ϑ (i.e. 〚 ϑ 〛 = { ( ω , t ) : ϑ ( ω ) = t } ).

● We say that ϑ is totally inaccessible if it avoids all ℍ -predictable stopping times (i.e., ℙ ( ϑ = ξ < ∞ ) = 0 for any ℍ -predictable stopping time ξ ).

We now work on a filtered probability space ( Ω , G , F , ℙ ) on which a random time τ is defined. We denote by Z the Azéma supermartingale (see ( [

Definition 2.4. Let G = ( G t ) t ≥ 0 be the progressive enlargement of F with τ , i.e., G = F ∨ A , which means that G t = ∩ ε > 0 G t + ε 0 , with G s 0 = F s ∨ A s for every s ≥ 0 (see, e.g., [

The filtration G is the smallest filtration satisfying the usual hypotheses containing F and turning out τ into a stopping time.

Definition 2.5. The F -predictable reduction of the compensator of τ

The process Λ given by

Λ t = ∫ 0 t 1 l { Z s − > 0 } d A s p , F Z s − (2.1)

is F -predictable and increasing,and,denoting by Λ τ the process Λ stopped at timeτ,

A t − Λ t ∧ τ = A t − Λ t τ = A t − ∫ 0 t ∧ τ d A s p , F Z s −

isa G -martingale (see ( [

The process Λ τ is the G -compensator of the default process A (we shall also say compensator of τ )and we call Λ the F -predictable reduction of the G -compensator ofτ.

If Λ is absolutely continuous with respect to the Lebesgue measure,i.e., Λ t = ∫ 0 t λ s d s ,then its derivative λ ,that is a non-negative F -predictable process,is called the F -intensity rate.

We recall that the random time τ avoids all F -stopping times (resp. all F -predictable stopping times) if and only if A o , F (resp. A p , F ) is continuous (see ( [

In a filtered probability space ( Ω , G , F , ℙ ) covering the market uncertainty, we consider an increasing sequence of F -stopping times ( θ i ) i , with θ 0 = 0 . We consider an increasing càdlàg process L (assumed to be independent of F ) to be the aggregate loss process related to a sequence of catastrophe events such that L 0 = 0 , L ∞ = ∞ .

Let τ be the catastrophe default time (known also as trigger event time) of a CAT bond contract.

Definition 3.1. A zero-coupon CAT bond with maturity T is a contract that pays a fixed amount (called the principal) P cat at time T if τ does not occur before T and a fraction δ (with 0 ≤ δ < 1 ) of the principal if τ occurs before T. Its payoff is then given by

ζ = P cat 1 l { τ > T } + δ P cat 1 l { τ ≤ T } .

We consider G to be the enlarged filtration of F with τ . Here G represents the extra information about the market. Hence, as one can see, the pricing of the CAT bond should be done in G . Fortunately, the enlargement of filtration theory offers a way to return the pricing in the reference filtration F . This requires the knowledge of some characteristics of the trigger event time τ such as the Azéma supermartingale and the dual predictable projection in F . Furthermore, these quantities depend on the model related to τ. In what follows, we introduce two models for which we give those quantities and discuss the pricing of zero-coupon CAT bonds.

We define the trigger event time τ as

τ = θ i on { L θ i − 1 ≤ D < L θ i } , for i ≥ 1 (3.1)

where D is a (finite) positive fixed amount representing a threshold value of the CAT bond. We consider the function Ψ related to the law of the aggregate loss process L as Ψ ( t , D ) : = ℙ ( L t ≤ D ) and note that Ψ ( ∞ , D ) = 1 hence, τ is almost surely finite.

An interesting feature of our approach is that the F -stopping times ( θ i ) i are not avoided by the trigger event time τ and these times can be fixed a priori so that they coincide with the catastrophe arrival times.

We start by computing the F -conditional survival law of the trigger event time τ which allows deriving the Azéma supermartingale Z associated with τ and then using the Doob-Meyer decomposition of Z for obtaining the dual predictable projection of τ . All these details related to the computation of the F -conditional survival law of τ can be seen in the proof of the following proposition.

Proposition 3.2. The Azéma supermartingale Z of the trigger event time τ is given by

Z t : = ℙ ( τ > t | F t ) = 1 − ∑ i = 1 ∞ 1 l { θ i ≤ t } [ Ψ ( θ i − 1 , D ) − Ψ ( θ i , D ) ] , for t ∈ ℝ + . (3.2)

Remark 3.3. Note the importance of the assumption (A) which implies that Z is non-increasing (this can be seen in (3.2)).

Proof: For all t , u ∈ ℝ + , one has

ℙ ( τ > u | F t ) = ∑ i = 1 ∞ ℙ ( θ i > u , L θ i − 1 ≤ D < L θ i | F t ) .

This implies that

ℙ ( τ > u | F t ) = ∑ i = 1 ∞ ℙ ( θ i > u , L θ i − 1 ≤ D | F t ) − ∑ i = 1 ∞ ℙ ( θ i > u , L θ i ≤ D | F t ) = ∑ i = 1 ∞ ℙ ( θ i > u , L θ i − 1 ≤ D | F t ) − ∑ i = 1 ∞ ℙ ( θ i − 1 > u , L θ i − 1 ≤ D | F t ) = ∑ i = 1 ∞ ℙ ( θ i > u ≥ θ i − 1 L θ i − 1 ≤ D | F t )

where the second equality is due to the fact that θ 0 = 0 .

Therefore, by using the tower property one obtains

ℙ ( τ > u | F t ) = ∑ i = 1 ∞ E [ 1 l { θ i > u ≥ θ i − 1 } ℙ ( L θ i − 1 ≤ D | F ∞ ) | F t ] = ∑ i = 1 ∞ E [ 1 l { θ i > u ≥ θ i − 1 } Ψ ( θ i − 1 , D ) | F t ] = E [ ∑ i = 0 ∞ 1 l { θ i ≤ u } Ψ ( θ i , D ) − ∑ i = 1 ∞ 1 l { θ i ≤ u } Ψ ( θ i − 1 , D ) | F t ] .

Since Ψ ( θ 0 , D ) = 1 , one has

∑ i = 0 ∞ 1 l { θ i ≤ u } Ψ ( θ i , D ) = 1 + ∑ i = 1 ∞ 1 l { θ i ≤ u } Ψ ( θ i , D )

Hence, it follows

ℙ ( τ > u | F t ) = 1 − E [ ∑ i = 1 ∞ 1 l { θ i ≤ u } [ Ψ ( θ i − 1 , D ) − Ψ ( θ i , D ) ] | F t ] .

If t ≥ u , one has

ℙ ( τ > u | F t ) = 1 − ∑ i = 1 ∞ 1 l { θ i ≤ u } [ Ψ ( θ i − 1 , D ) − Ψ ( θ i , D ) ]

which is due to the fact that the random variables 1 l { θ i ≤ u } and 1 l { θ i ≤ u } [ Ψ ( θ i − 1 , D ) − Ψ ( θ i , D ) ] are F t -measurable.

In particular, we have

Z t : = ℙ ( τ > t | F t ) = 1 − ∑ i = 1 ∞ 1 l { θ i ≤ t } [ Ψ ( θ i − 1 , D ) − Ψ ( θ i , D ) ] .

We can also easily check that

ℙ ( τ > u | F u ) = ℙ ( τ > u | F ∞ ) ,

which implies the immersion property of the model. □

In what follows, we compute the F -dual predictable projection of τ using the Doob-Meyer decomposition of Z. This decomposition depends on the nature of the F -stopping times ( θ i ) i .

Case where some of the ( θ i ) i are not predictable:

We denote, for any i , by Λ i the F -compensator of θ i , i.e., the F -predic- table increasing process Λ i , with Λ 0 i = 0 , such that ( 1 l { θ i ≤ t } − Λ t ∧ θ i i , t ≥ 0 ) is an F -martingale. If it exists, we denote the intensity rate of θ i by λ i (i.e., Λ t ∧ θ i i = ∫ 0 t ∧ θ i λ s i d s , for all t ≥ 0 ).

Proposition 3.4. The F -dual predictable projection A p , F of the trigger event time τ is given by

A t p , F = ∑ i = 1 ∞ ∫ 0 t ( Ψ ( θ i − 1 , D ) − Ψ ( s , D ) ) λ s i 1 l { θ i − 1 ≤ s < θ i } d s , for t ∈ ℝ + . (3.3)

Proof: This proposition can be proved as the same way as Proposition 3.15 of [

Case where the ( θ i ) i are predictable:

In this case, the Azéma supermartingale Z is predictable, and one has A t p , F = 1 − Z t , for all t ∈ ℝ + (i.e., m = 1 , almost surely). In addition, the trigger event time covered by the ( θ i ) i is accessible.

In the next subsections, we discuss the pricing of zero-coupon bonds by taking into account the two cases.

For simplicity, we assume a constant interest rate r and ℙ being the pricing measure. Hence, the price V t ( T ) of the zero-coupon CAT bond with principal P cat and maturity T at time t ≤ T is

V t ( T ) = e r t E [ P cat e − r T 1 l { T < τ } + e − r τ δ P cat 1 l { t < τ ≤ T } | G t ] .

According to ( [

V t ( T ) e − r t = 1 l { t < τ } P cat e − r T 1 Z t E [ Z T | F t ] + 1 l { t < τ } δ P cat 1 Z t E [ ∫ t T e − r u d A u p , F | F t ] , ∀ 0 ≤ t ≤ T . (3.4)

By consequence, one has

V 0 ( T ) = P cat e − r T E [ Z T ] + δ P cat E [ ∫ 0 T e − r u d A u p , F ] = P cat e − r T ( 1 − E [ A T p , F ] ) + δ P cat E [ ∫ 0 T e − r u d A u p , F ] (3.5)

where we have used the fact that Z = m − A p , F and since A 0 p , F = 0 then E [ m T ] = E [ m 0 ] = E [ Z 0 ] = 1 almost surely.

Note that the main quantity we need here for obtaining the price at time 0 of the zero-coupon CAT bond is the dual predictable projection of the trigger event time τ . According to equality (3.4), this quantity depends on the law of the aggregate losses, the catastrophe arrival times with their intensity rates, and the threshold level.

By combining (3.4) and (3.5), one obtains

V 0 ( T ) = P cat e − r T ( 1 − ∑ i = 1 ∞ ∫ 0 T E [ Q i ( s ) ] d s ) + δ P cat ∑ i = 1 ∞ ∫ 0 T e − r u E [ Q i ( u ) ] d u , (3.6)

where Q i ( t ) : = ( Ψ ( θ i − 1 , D ) − Ψ ( t , D ) ) λ t i 1 l { θ i − 1 ≤ t < θ i } , for i ≥ 1 .

A particular framework of our model

Here, we consider a particular case with the two F -stopping times θ 0 and θ 1 with θ 0 = 0 (i.e., i ∈ { 1,2 } ). Then from (3.2) the Aéma supermartingale Z has the following expression

Z t = { Ψ ( θ 1 , D ) if θ 1 ≤ t 1 if θ 1 > t .

From (3.4) the F -dual predictable projection A p , F of the trigger event time τ is given by

A t p , F = ∫ 0 t ( 1 − Ψ ( s , D ) ) λ s 1 1 l { 0 ≤ s < θ 1 } d s , for t ∈ ℝ + .

Therefore, the price V 0 ( T ) of the zero-coupon CAT bond that pays a principal P cat at T if τ does not occur before T and zero otherwise (i.e. δ = 0 ), is then given by

V 0 ( T ) = P cat e − r T ( 1 − ∫ 0 T E [ ( 1 − Ψ ( s , D ) ) λ s 1 1 l { 0 ≤ s < θ 1 } ] d s ) (3.7)

which is similar to the zero-coupon bond price at time 0 obtained in [

Comment 3.5. Note in this case that our model is a more general framework according to some existing CAT bond models in the literature (such as the ones in [

In this case, the main quantity for pricing zero-coupon bonds is the Azéma supermartingale Z, and since Z = 1 − A p , F , the equality (3.4) can be replaced by

V t ( T ) e − r t = 1 l { t < τ } P cat e − r T 1 Z t E [ Z T | F t ] − 1 l { t < τ } δ P cat 1 Z t E [ ∫ t T e − r u d Z u | F t ] , ∀ 0 ≤ t ≤ T .

In the case with zero interest rate, one has V 0 ( T ) = P cat ( 1 − δ ) E [ Z T ] , i.e.,

V 0 ( T ) = P cat ( 1 − δ ) E [ 1 − ∑ i = 1 ∞ 1 l { θ i ≤ T } ( Ψ ( θ i − 1 , D ) − Ψ ( θ i , D ) ) ] . (3.8)

Note in this case that the price depends on the law of the aggregate losses, the catastrophe arrival times, and the threshold level. In the particular case where the ( θ i ) i are deterministic, the computation is simply based on the computation of the aggregate loss distribution. Note that, in general, closed-form solutions are not easy to obtain for this distribution but numerical algorithms, such as Monte Carlo and Fourier transformation can be successfully used to estimate it. However, the problem of estimating that distribution arises when there is low available data, which is the case in general for CAT bonds.

We recall that given an F -survival process Z and a uniform random variable U in [ 0,1 ] independent from F , one can construct a random time τ associated with Z by extended the standard Cox construction as (see, e.g., in [

τ : = inf { t ≥ 0 : Z t ≤ U } . (3.9)

In this setting, immersion holds and we have

ℙ ( τ > t | F t ) = Z t .

We consider the survival process of the investor in the CAT bond introduced above to be defined by

Z t : = e − L t D , ∀ t ≥ 0 (3.10)

and we model the trigger event time of the investor as in (3.9). It is clear that this is non-increasing with respect to t and 0 ≤ Z t ≤ 1 , for all t ≥ 0 . Furthermore, Z is non-decreasing with respect to the threshold value D. This is a valid survival process in CAT bond modeling. Let us note that the definition (3.9) is equivalent to define to as

τ : = inf { t ≥ 0 : L t D > Θ } (3.11)

where Θ is a unit exponential random variable independent of the reference filtration F . We postulate that the F -stopping times ( θ i ) i are the jump times of the process L.

This model belongs to the so-called Generalized Cox model in credit risk (see

[

time τ can be found following [

We consider the F θ i -measurable non-negative random variable y i ^{1} to be the amount of losses at the catastrophe time θ i and we denote by μ the jump measure of the marked point process ( θ i , y i ) and ν its compensator (for simplicity, we suppose ν to be deterministic^{2}. We define the aggregate loss process as

L t : = ∑ i ≥ 1 1 l { θ i ≤ t } H ( t − θ i , y i ) = ∫ 0 t ∫ ℝ H ( t − s , x ) μ ( d s , d x ) , ∀ t ≥ 0 , (3.12)

where H is a function ℝ + × ℝ → ℝ + with

H ( t , x ) = H ( 0 , x ) + ∫ 0 t h ( s , x ) d s , ∀ t ≥ 0 , x ∈ ℝ , (3.13)

where h is a non-negative Borel function on ℝ + × ℝ . We assume that

∫ 0 T ∫ ℝ h 2 ( s , x ) ν ( d s , d x ) < ∞ , ∀ T < ∞ . (3.14)

The equality (3.13) guarantees the fact that L is increasing with respect to the time direction and the relation (3.14) insures the semi-martingale property of L (see, e.g., Lemma 2 of [

By using Proposition 3.15 of [

Z t ( u ) : = ℙ ( τ > u | F t ) = c ( u ) L t ( u ) , for u ≥ t (3.15)

where, c ( u ) = exp ( ∫ 0 u ∫ ℝ ( e − H ( u − s , x ) D − 1 ) ν ( d s , d x ) ) , ∀ u ∈ ℝ + and L t ( u ) = exp ( − ∫ 0 t ∫ ℝ H ( u − s , x ) D μ ( d s , d x ) − ∫ 0 t ∫ ℝ ( e − H ( u − s , x ) D − 1 ) ν ( d s , d x ) ) , for any u ∈ ℝ + which is an F -martingale. In particular, the survival function of τ is ℙ ( τ > u ) = c ( u ) and the Azéma supermatingale is Z t = Z t ( t ) .

Following the proposition 3.14 of [

A t p , F = ∫ 0 t ∫ ℝ Z s − ( e − H ( 0, x ) D − 1 ) ν ( d s , d x ) (3.16)

and the F -predictable reduction Λ of the compensator of τ satisfies

d Λ t = ∫ ℝ ( 1 − e − H ( 0 , x ) D ) ν ( d t , d x ) , Λ 0 = 0.

Furthermore, if v is continuous, the trigger event time τ admits an intensity rate λ . For instance, if ( θ i ) i are the jump times of a Poisson process N with deterministic intensity function λ N such that ν ( d t , d x ) : = F ( d x ) λ N ( t ) d t , where F is a distribution function, then

λ t = λ N ( t ) ∫ ℝ ( 1 − e − H ( 0, x ) D ) f ( x ) d x , for all t ≥ 0. (3.17)

Hence the intensity of the trigger event time depends on the threshold value D, the loss severity distribution function F, and the claim arrival intensity λ N . Note that the intensity rate of τ is time-varying if the intensity of the claim arrival is time-varying. Furthermore, this quantity is non-increasing with respect to D.

We consider a non-negative stochastic interest rate r adapted to a Brownian filtration F W independent of the Shot-Noise process L defined in (3.12). By denoting F L as the filtration generated by the L, we define the reference filtration F as F = F W ∨ F L . For simplicity, we consider that the principal P cat (that we assume to be equal to 1) is fully lost in case of occurrence of the trigger event, i.e., δ = 0 . In this setting, we have the following result.

Lemma 3.6. The price of the zero-coupon CAT bonds is given by

V t ( T ) = 1 l { t < τ } exp ( ∫ t T ∫ ℝ ( e − H ( T − s , x ) D − 1 ) ν ( d s , d x ) − 1 D ∫ 0 t ∫ ℝ [ H ( T − s , x ) − H ( t − s , x ) ] μ ( d s , d x ) ) Q t ( T ) ,

with Q t ( T ) : = E [ e − ∫ t T r s d s | F t W ] .

Proof: The price V t ( T ) of the zero-coupon CAT bond is given by

V t ( T ) = E [ P cat 1 l { τ > T } e − ∫ t T r s d s | G t ] .

Under the general framework of pricing defaultable zero-coupon bonds of ( [

V t ( T ) = 1 l { t < τ } 1 Z t E [ Z T e − ∫ t T r s d s | F t ] = 1 l { t < τ } 1 Z t E [ e − ∫ t T r s d s E [ Z T | F T W ∨ F t L ] | F t ]

where we have used the tower property and the fact that fact that e − ∫ t T r s d s is F T W ∨ F t L -measurable.

Since L is independent of F W (hence Z is independent of F W ) and F t L ⊂ F t , one obtains

V t ( T ) = 1 l { t < τ } 1 Z t E [ e − ∫ t T r s d s E [ Z T | F t L ] | F t ] = 1 l { t < τ } 1 Z t E [ Z T | F t L ] E [ e − ∫ t T r s d s | F t ] = 1 l { t < τ } 1 e − L t D E [ e − L T D | F t L ] E [ e − ∫ t T r s d s | F t W ] .

We have

E [ e − L T D | F t L ] = E [ e − L T D | F t ] = Z t ( T ) = c ( T ) L t ( T ) ,

where c ( T ) = exp ( ∫ 0 T ∫ ℝ ( e − H ( T − s , x ) D − 1 ) ν ( d s , d x ) ) and

L t ( T ) = exp ( − ∫ 0 t ∫ ℝ H ( T − s , x ) D μ ( d s , d x ) − ∫ 0 t ∫ ℝ ( e − H ( T − s , x ) D − 1 ) ν ( d s , d x ) ) . Hence, the result follows by replacing L t by its value given in (3.12). □

A very tractable Shot-Noise model called the stochastic discounting model (see [

L t = ∑ i ≥ 1 1 l { θ i ≤ t } y i e − α ( θ i − t ) , (3.18)

with ( y i ) some non-negative random variables and α a strictly positive parameter. In this case, the process L is a Markovian Shot-Noise process since H ( T − s , x ) : = x e α ( T − s ) = e α ( T − t ) H ( t − s , x ) (see [

V t ( T ) = 1 l { t < τ } exp ( ∫ t T ∫ ℝ + ( e − x e α ( T − s ) D − 1 ) ν ( d s , d x ) − 1 D ( e α ( T − t ) − 1 ) L t ) Q t ( T ) (3.19)

where Q t ( T ) : = E [ e − ∫ t T r s d s | F t W ] .

Here we say semi-closed form because of the quantity Q ( T ) but closed form can be obtained when the interest rate process is an affine process or a polynomial one.

For example, in the case where the interest rate follows a Cox-Ingersoll-Ross (CIR) model, i.e.,r verifies

d r t = γ r ( θ − r t ) d t + σ r t d W t , r 0 = x ,

where γ r , θ , and σ are positive parameters and W a Brownian motion, then Q t ( T ) = e A t ( T ) − B t ( T ) r t , where A and B verify

A t ( T ) = 2 γ r θ σ 2 ln ( 2 h e 1 2 ( γ r + h ) ( T − t ) h − γ r + e h ( T − t ) ( h + γ r ) ) and B t ( T ) = 2 ( e h ( T − t ) − 1 ) h − γ r + e h ( T − t ) ( h + γ r ) , where h = γ r 2 + 2 σ 2 .

Comment 3.7. One of the specificities of this model is that the arrival catastrophe events induce some jumps in the CAT bond prices at the time of occurrence of those events. Indeed let us consider the F -adapted process V ˜ ( T ) such that V t ( T ) 1 l { t < τ } = V ˜ t ( T ) 1 l { t < τ } . This process V ˜ ( T ) always exists and admits some negative jumps with sizes given by

Δ V ˜ θ i = V ˜ θ i − ( e − ( H ( T − θ i , y i ) D − H ( 0, y i ) D ) − 1 ) . (3.20)

This can be easily seen by using the results of example 4.2 of [

We now consider the F -stopping times ( θ i ) i to be the jumps times of a time-inhomogeneous Poisson process N with intensity function λ N and consider the compensator measure ν ( d t , d x ) : = F ( d x ) λ N ( t ) d t , where F is a distribution function. Therefore, the price of the CAT bond given in (3.19) becomes

V t ( T ) = 1 l { t < τ } exp ( ∫ t T ∫ ℝ + ( e − x e α ( T − s ) D − 1 ) F ( d x ) λ N ( s ) d s − 1 D ( e α ( T − t ) − 1 ) L t ) Q t ( T ) . (3.21)

Hence, the price of the zero-coupon CAT bond is affected by the threshold value D, the claim arrival intensity λ N , the interest rate uncertainty, the distribution of the shots F, and the speed of the growth α of the impulsion function H of the Shot-Noise. These two last parameters constitute the severities of the losses.

Furthermore, if λ N is constant, we have

V t ( T ) = 1 l { t < τ } exp ( λ N ∫ t T ( ϕ ( e α ( T − s ) D ) − 1 ) d s − 1 D ( e α ( T − t ) − 1 ) L t ) Q t ( T ) (3.22)

where ϕ ( u ) is the Laplace transform of the distribution F evaluated at u ∈ ℝ + . The Laplace transform is not always explicit but some approximations have been introduced for some distributions. For instance, if F is the log-normal distribution with parameters μ ∈ ℝ and σ > 0 then a closed-form of ϕ does not exist but numerical approximations can be used to compute ϕ (see, e.g., in [

[

ϕ ( J ( s ) ) = e − μ J ( s ) 1 + W ( σ 2 J ( s ) ) exp ( − 1 2 σ 2 W ( σ 2 J ( s ) ) 2 − 1 σ 2 W ( σ 2 J ( s ) ) )

where W is the Lambert function (see, e.g., [

When F is a Pareto distribution function, i.e., in the form

F ( x ) = 1 − b a b + x , x > 0 , a > 0 , b > 0.

Then ϕ can be setting using the result of [

ϕ ( J ( s ) ) = a ( b J ( s ) ) ( a − 1 ) / 2 exp ( b J ( s ) / 2 ) W − ( a + 1 ) / 2 , − a / 2 ( b J ( s ) )

where W κ , μ is the Whittaker function (see [

Comment 3.8.

● In the case where α = 0 , the ( θ i ) i being the jumps times of an homogeneous Poisson process N with intensity λ N and the shots ( γ i ) i being i.i.d and independent of N, the Shot-Noise process L defined in (3.18) is a Compound Poisson process.

● In this case the CAT bond pre-trigger price is continuous.

In this subsection, we compute the term structure of the pre-trigger price of a zero-coupon CAT bond with principal P cat = 1 and maturity of 3 years using the Shot-Noise model. As such, we use a threshold value D = 10000 and a CIR interest rate model with parameters ( r 0 = 0.0204 , θ = 0.0204 , γ r = 0.0884 , and σ = 0.0477 ) extracted from [

In Figure1, we represent the term structure of the zero-coupon CAT bond pre-trigger price (the blue curve) as a function of time in years along with the corresponding path of the aggregate losses (in red) generated from the Shot- Noise process L. We note two jumps in aggregate losses at times θ 1 = 1.104 and θ 2 = 1.971 with jump sizes respectively to equal 601.8668 and 582.0399. This automatically induces negative jumps at the same times, in the price of CAT bond with sizes corresponding respectively to 0.17989792 and 0.05843795 computed using the equality (3.20). One can be surprised at the difference in jump sizes during the two events. Indeed, despite the larger size of the second jump of the Shot-Noise, the size of the first price jump remains higher than that of the first. However, in terms of ratio, we can see that it is reasonable since the size of the first price jump represents 19.309152% of its last value before the jump while that of the second is 7.092837% of its last value just before the moment of the jumps. We generally observe an increasing trend for the CAT bond price in between the jump times. For more illustrations, FigureA1 of the Appendix shows different scenarios for the behavior of the time-varying CAT bond price according to the trajectory of the aggregate losses.

Now for a given trajectory of aggregated losses, we fix different threshold levels (5000, 9000, 15,000, 20,000) and compute for each Threshold the value of the zero-coupon CAT bond at time 0 ≤ t ≤ T . Figure2 presents the time-varying CAT bond price according to these levels of threshold. Unsurprisingly, the results show that the CAT bond price is non-decreasing with respect to the threshold values. Indeed, higher threshold value leads to low probability of exceeding it, which leads to the increase in the CAT price. This satisfies the results obtained in [

From a standpoint to show the importance of using the Shot-Noise process in our framework, we sample one path of the aggregate losses from a Compound Poisson process (CPP) using a log-normal distribution for the severity distribution and a homogeneous Poisson process with the same parameters as the ones used for the Shot-Noise (i.e., μ = 6.387 , σ = 0.153 and λ N = 0.5 ). We also compute the time-varying pre-trigger price of the CAT bond introduced above using the CPP and the same trajectory of the CIR interest rate. Figure3 represents the path of the CAT bond price (black curve) computed using the CPP with the corresponding trajectory of the CPP (green curve) and a trajectory of the CAT bond price obtained from the Shot-Noise (in blue) also along with the corresponding Shot-Noise path (in red). The results show two jumps of the compound Poisson process at times θ 1 = 0.549 and θ 2 = 1.674 with jump sizes respectively equal to 599.3870 and 608.2323. Despite the jumps of the CPP, we observe a continuity everywhere of its corresponding CAT bond price. The last finding shows that it is more convenient to use Shot-Noise processes in our framework than using a Compound Poisson process for the aggregate loss modeling. See FigureA3 of the Appendix for different scenarios.

In

In this paper, we have explored the abilities of two models related to credit risk modeling for CAT bond pricing. In Model 1, we have first investigated the case where the trigger event time is covered by totally inaccessible stopping times that make it to be totally inaccessible. We then study the second case where those stopping times are predictable hence the trigger event time is also predictable. We have derived all the quantities of interest for pricing zero-coupon CAT bonds. We also showed that, in some settings, the model generalizes some ones developed in the CAT bond modeling.

In Model 2 we started with the survival process of an investor in a CAT bond and then based on the results of enlargement of filtration, we construct the trigger event time. We have shown that this framework is related to the generalized Cox model [

While the two approaches offer a new perspective on modeling CAT bonds, some potential improvements have to be taken into account in the forthcoming studies.

In Model 1, we have shown that the zero-coupon CAT bond prices depend on the distribution of the aggregate losses. However, the estimation of this distribution could lead to some problems because of the low available data in the CAT bond framework. An adapted numerical approximation could be achieved for calibrating this model in real data.

In Model 2, we did not interested in the estimation of the Shot-Noise since several approaches to estimating that process exists in the literature (see, e.g., [

The authors declare no conflicts of interest regarding the publication of this paper.

Chaieb, Z. and Gueye, D. (2022) Pricing Zero-Coupon CAT Bonds Using the Enlargement of Filtration Theory: A General Framework. Journal of Mathematical Finance, 12, 582-605. https://doi.org/10.4236/jmf.2022.123031