Credit Risk Modeling Approaches › Reduced-Form Models: Intensity-Based Approach

Reduced-form models take a fundamentally different approach from structural models. Rather than deriving default from balance sheet dynamics and capital structure, they model default as a random event governed by a hazard rate (or default intensity). Default occurs exogenously, without explicit modeling of economic mechanisms. Instead of asking why the firm defaults, reduced-form models focus on when: estimating the instantaneous default probability at any moment given that the firm has survived to that time.

This approach offers several practical advantages. First, it's easier to calibrate to market prices since the hazard rate connects directly to observable credit spreads. Second, the framework handles multiple credit events naturally (including rating downgrades and restructurings), not just binary default. Third, it doesn't require unobservable asset values, making implementation more straightforward.

The mathematical foundation of reduced-form models rests on the concept of a hazard rate, which originates from survival analysis in statistics and actuarial science. The default time $\tau$ is the first time a default-triggering event occurs. The survival probability to time $t$ is the likelihood that the obligor remains solvent through time $t$:

$$S(t) = P(\tau > t)$$

where:

• $S(t)$: probability of surviving (not defaulting) until time $t$
• $\tau$: random default time
• $P(\tau > t)$: probability that default occurs after time $t$

The survival probability function captures the declining likelihood of remaining solvent as time progresses. It starts at $S(0) = 1$ (the firm is solvent today) and decreases toward zero as the time horizon extends (eventually all firms default or cease to exist). The rate at which this probability declines at any moment is governed by the hazard rate.

The hazard rate has an intuitive interpretation as the instantaneous default probability per unit time: if the hazard rate is 5% per year, it means that at any moment, a firm that has survived faces approximately a 5% chance of defaulting over the next year. This probability is conditional on survival: it applies only to firms that haven't already defaulted.

The survival probability relates to the hazard rate through the exponential of the cumulative hazard. To see why, note that the survival probability must satisfy $\frac{dS(t)}{dt} = -\lambda(t) S(t)$ because the rate of decrease in survival probability equals the hazard rate times the probability of having survived so far. This differential equation simply states that the rate at which firms leave the surviving population equals the hazard rate times the population size. Solving this differential equation with initial condition $S(0) = 1$ gives:

$$S(t) = \exp\left(-\int_0^t \lambda(s) \, ds\right)$$

where:

• $S(t)$: survival probability to time $t$
• $\lambda(s)$: hazard rate (default intensity) at time $s$
• $s$: integration variable representing time from 0 to $t$
• $\int_0^t \lambda(s) \, ds$: cumulative hazard, the total default risk accumulated from time 0 to $t$
• $\exp(\cdot)$: exponential function
• The exponential function converts cumulative hazard to survival probability through the relationship $S(t) = e^{-\text{cumulative hazard}}$
• This relationship follows from solving the differential equation $\frac{dS(t)}{dt} = -\lambda(t) S(t)$ with initial condition $S(0) = 1$

The integral $\int_0^t \lambda(s) ds$ represents the cumulative hazard, capturing the total default risk accumulated from time 0 to time $t$. The exponential transformation converts this accumulated risk into a survival probability. Higher cumulative hazard means lower survival probability, with the exponential ensuring the probability remains between 0 and 1.

For a constant hazard rate $\lambda$, this simplifies to:

$$S(t) = e^{-\lambda t}$$

where the integral becomes simply $\lambda t$, giving exponential decay in survival probability. This exponential decay is the hallmark of the constant hazard model: the survival probability decreases at a rate proportional to itself.

The cumulative default probability is:

$$P(\tau \leq t) = 1 - S(t) = 1 - e^{-\lambda t}$$

where:

• $P(\tau \leq t)$: probability of defaulting by time $t$
• $\tau$: random default time
• $S(t) = e^{-\lambda t}$: survival probability with constant hazard
• $1 - S(t)$: complement gives cumulative default probability

The constant hazard rate model implies default arrives as a Poisson process with intensity $\lambda$. This connection provides deep insight into reduced-form model mechanics. The Poisson process describes events that occur randomly over time with a constant average rate, and the exponential distribution $P(\tau \leq t) = 1 - e^{-\lambda t}$ characterizes the waiting time until the first such event.

The key property of the Poisson process is the memoryless property. Given that a firm has survived to time $t$, the probability distribution of remaining time to default is the same as if we started fresh at time zero. This may seem counterintuitive (shouldn't older firms be closer to default?), but it captures the idea that conditional on survival, the firm faces the same ongoing default intensity. The waiting time until default follows an exponential distribution with mean $1/\lambda$, so a firm with a 2% annual hazard rate has an expected time to default of 50 years, though actual default could occur much sooner or later due to the high variance of exponential distributions.

The survival probability declines from 100% to 70.8% over 10 years, while cumulative default probability rises from 0% to 29.2%. The upward-sloping hazard-rate term structure, starting at 2% annually and rising to 4%, reflects the common empirical pattern that credit uncertainty increases with maturity. The continuous curves change slope at years 1, 2, 3, and 5 as the piecewise-constant hazard rate steps upward.

The primary advantage of reduced-form models is their ability to be calibrated directly to market data. As we discussed in Part II Chapter 13 on credit default swaps, CDS spreads reflect the market's assessment of default risk. The CDS market provides a direct window into market-implied default probabilities, making reduced-form models highly practical for trading and risk management applications.

For a CDS with spread $s$ and recovery rate $R$, no-arbitrage pricing requires that the present value of premium payments (premium leg) equals the present value of protection payments (protection leg). This equilibrium condition reflects the fact that a fairly priced CDS has zero initial value: neither party pays upfront because the expected benefits equal the expected costs. This equilibrium allows us to back out the implied hazard rate from observable CDS spreads.

The premium leg represents the periodic payments made by the protection buyer.

$$\text{Premium\_Leg} = s \sum_{i=1}^n \Delta t_i \cdot S(t_i) \cdot D(t_i)$$

where:

• $s$: CDS spread (annualized premium rate paid by protection buyer)
• $n$: number of premium payment periods
• $i$: index for payment periods, running from 1 to $n$
• $\Delta t_i$: length of payment period $i$ in years
• $t_i$: time of the $i$-th payment
• $S(t_i)$: survival probability to time $t_i$ (premium only paid if no default has occurred)
• $D(t_i)$: discount factor to time $t_i$, equal to $e^{-r t_i}$ for continuous discounting
• The sum captures expected present value of all premium payments

The summation reflects that premiums are only paid if the reference entity has not defaulted. Each term multiplies the premium rate by the payment period length, the probability of survival to that date, and the appropriate discount factor. The survival probability weighting removes premiums scheduled after default: if default occurs before payment date $t_i$, that premium is never paid.

$$\text{Protection\_Leg} = (1-R) \int_0^T \lambda(t) S(t) D(t) \, dt$$

where:

• $R$: recovery rate (fraction of notional recovered upon default)
• $(1-R)$: loss given default (fraction not recovered)
• $T$: CDS maturity
• $\lambda(t)$: hazard rate (default intensity) at time $t$
• $S(t)$: survival probability to time $t$
• $D(t)$: discount factor for payment at time $t$
• $\lambda(t) S(t) dt$: probability of default occurring in the infinitesimal interval $[t, t+dt]$
• The integral sums expected present value of protection payments across all possible default times from 0 to $T$

where $D(t)$ is the discount factor to time $t$. The protection leg integral captures the expected payout to the protection buyer. Default can occur at any time between now and maturity, and when it does, the protection seller pays $(1-R)$ times the notional. The integrand $\lambda(t) S(t)$ represents the probability density of defaulting at exactly time $t$: the hazard rate times the probability of having survived to that point.

For a flat hazard rate $\lambda$, we can derive an approximate relationship between CDS spread and hazard rate. Under the assumptions that discount rates are small and survival probabilities do not vary too much over the CDS tenor.

$$\begin{aligned} \text{Premium\_Leg} &\approx s \sum_{i=1}^n \Delta t_i \cdot S(t_i) && \text{(sum of premium payments weighted by survival)} \\ &\approx s \cdot T \cdot \bar{S} && \text{(approximate as constant average survival } \bar{S} \text{)} \end{aligned}$$

Similarly, the protection leg becomes:

$$\begin{aligned} \text{Protection\_Leg} &\approx (1-R) \int_0^T \lambda S(t) \, dt && \text{(integral of loss payment hazard)} \\ &\approx (1-R) \lambda T \bar{S} && \text{(approximate with constant } \bar{S} \text{)} \end{aligned}$$

Setting these equal (no-arbitrage condition):

$$s \cdot T \cdot \bar{S} = (1-R) \lambda T \bar{S}$$

Canceling $T \bar{S}$ from both sides gives:

$$s \approx \lambda (1 - R)$$

where:

• $s$: CDS spread (annualized premium rate)
• $\lambda$: constant hazard rate (default intensity per year)
• $R$: recovery rate
• $(1-R)$: loss given default (LGD)
• The approximation holds when discount rates and survival probabilities are relatively stable over the CDS tenor

This relationship shows that CDS spreads compensate for two components of credit risk: the probability of default ($\lambda$) and the severity of loss when default occurs ($1-R$). The approximation becomes increasingly accurate for shorter maturities and lower hazard rates where discounting effects are minimal.

Solving for the hazard rate provides a direct link between observable CDS spreads and default intensity:

$$\lambda \approx \frac{s}{1 - R}$$

where:

• $\lambda$: implied hazard rate (annual default intensity)
• $s$: observed CDS spread
• $R$: assumed recovery rate
• The formula shows that CDS spreads reflect both default risk ($\lambda$) and loss severity ($1-R$)

This inversion provides a direct estimate of market-implied default intensity from a CDS spread under the flat-hazard approximation, making the reduced-form model straightforward to calibrate. The hazard rate interpretation also helps explain why CDS spreads differ across credits: investment-grade companies have low spreads because their hazard rates are low, while distressed companies have high spreads because their default intensities are high.

The visualization reveals how recovery rate assumptions affect hazard rate inference. For a given CDS spread, a lower recovery rate implies a lower hazard rate because the expected loss per default is higher. At 500 bps, assuming 30% recovery yields an implied hazard rate of approximately 7.1%, while assuming 50% recovery yields approximately 10%. This sensitivity highlights the importance of accurate recovery rate assumptions when calibrating reduced-form models, as the same CDS spread can imply substantially different default probabilities depending on the assumed recovery.

For more realistic dynamics, the hazard rate itself can be modeled as a stochastic process. Rather than assuming a constant default intensity, we allow it to fluctuate over time in response to changing economic conditions. A firm's credit quality varies with business cycles, competitive pressures, management decisions, and countless other factors. Capturing this time variation requires modeling the hazard rate as a random process.

A common specification is the Cox-Ingersoll-Ross (CIR) process, which ensures the hazard rate remains positive and mean-reverts to a long-term level:

$$d\lambda_t = \kappa(\theta - \lambda_t)dt + \sigma\sqrt{\lambda_t}dW_t$$

where:

• $\lambda_t$: hazard rate (default intensity) at time $t$
• $d\lambda_t$: infinitesimal change in the hazard rate over time $dt$
• $\kappa$: speed of mean reversion (how quickly $\lambda_t$ returns to its long-term level $\theta$)
• $\theta$: long-term mean level of the hazard rate
• $\sigma$: volatility parameter controlling the magnitude of random fluctuations
• $W_t$: standard Brownian motion driving the random shocks
• $dt$: infinitesimal time increment
• $dW_t$: infinitesimal Brownian increment
• $\kappa(\theta - \lambda_t)$: mean reversion term that pulls $\lambda_t$ toward $\theta$ when $\lambda_t < \theta$ and pushes it down when $\lambda_t > \theta$
• $\sigma\sqrt{\lambda_t}$: square root diffusion term that ensures the hazard rate remains non-negative (volatility scales with the level)

The CIR specification has several desirable properties that make it well-suited for modeling credit risk. The square root diffusion $\sigma\sqrt{\lambda_t}$ ensures non-negativity: as $\lambda_t$ approaches zero, the volatility shrinks proportionally, making it mathematically impossible for the process to go negative. This is essential since a negative default probability has no economic meaning. The mean reversion term $\kappa(\theta - \lambda_t)$ pulls the hazard rate toward its long-term mean $\theta$, capturing the empirical observation that credit quality tends to revert over economic cycles. When $\lambda_t < \theta$, the drift is positive (credit quality deteriorates toward normal), and when $\lambda_t > \theta$, the drift is negative (credit quality improves from distressed levels). This mean-reverting behavior reflects the idea that extremely good or extremely bad credit conditions are temporary states that tend to normalize over time.