GARCH(1,1) Stationarity, Unconditional Variance, and Shock Propagation
You model daily log returns $r_t$ using a GARCH(1,1) process:
$\sigma_t^2 = \omega + \alpha \, r_{t-1}^2 + \beta \, \sigma_{t-1}^2$
where $\omega$, $\alpha$, and $\beta$ are parameters and $r_t = \sigma_t z_t$ with $z_t \sim N(0,1)$ i.i.d.
- State the conditions on $\omega$, $\alpha$, $\beta$ that guarantee covariance stationarity.
- Derive the unconditional variance $E[\sigma_t^2]$ under stationarity.
- Explain how a volatility shock $r_{t-1}^2$ propagates forward through time and how the half-life of a shock relates to the persistence parameter $\alpha + \beta$.
Hints
- Think of the GARCH variance equation as a weighted average of three components. What constraint on the weights ensures the process does not explode?
- Under stationarity, $E[\sigma_t^2] = E[\sigma_{t-1}^2]$. Use this plus the law of iterated expectations (noting $E[r_t^2] = E[\sigma_t^2]$) to solve for the unconditional variance.
- Iterate the GARCH equation forward and track how the deviation $\sigma_{t+k}^2 - \bar{\sigma}^2$ evolves. Set $(\alpha + \beta)^h = 1/2$ and solve for $h$.
Worked Solution
How to Think About It: GARCH(1,1) is the workhorse volatility model on every trading desk. The intuition is simple: tomorrow's variance is a weighted average of a long-run baseline ($\omega$), yesterday's squared return (the "news" or shock, weighted by $\alpha$), and yesterday's variance forecast (the "memory," weighted by $\beta$). The key number is $\alpha + \beta$ -- that is the persistence. If it equals or exceeds 1, shocks never die out and there is no well-defined long-run variance. If it is less than 1, shocks decay geometrically and you get mean-reversion in volatility.
Quick Sanity Checks: For equity indices, typical fitted values are $\alpha \approx 0.05$-$0.10$, $\beta \approx 0.85$-$0.90$, giving $\alpha + \beta \approx 0.95$-$0.99$. This means volatility is highly persistent but stationary. The unconditional variance should be a finite positive number, and shocks should have a half-life on the order of weeks to months.
---
(a) Conditions for Covariance Stationarity
The conditional variance must be positive almost surely and finite in expectation. The conditions are:
- $\omega > 0$
- $\alpha \geq 0$
- $\beta \geq 0$
- $\alpha + \beta < 1$
The first three ensure $\sigma_t^2 > 0$ for all $t$. The critical condition is the fourth: $\alpha + \beta < 1$ guarantees that the variance process is mean-reverting and the unconditional moments exist.
To see why, note that GARCH(1,1) can be written as an ARMA(1,1) in $r_t^2$. Define $v_t = r_t^2 - \sigma_t^2$ (the variance "surprise"). Then:
$r_t^2 = \omega + (\alpha + \beta) \, r_{t-1}^2 + v_t - \beta \, v_{t-1}$
This is an ARMA(1,1) process for $r_t^2$ with AR coefficient $(\alpha + \beta)$. Covariance stationarity of an ARMA process requires the AR root to lie inside the unit circle, i.e., $\alpha + \beta < 1$.
---
(b) Unconditional Variance
Under stationarity, the unconditional expectation satisfies $E[\sigma_t^2] = E[\sigma_{t-1}^2]$. Take expectations of the GARCH equation:
$E[\sigma_t^2] = \omega + \alpha \, E[r_{t-1}^2] + \beta \, E[\sigma_{t-1}^2]$
Since $E[r_t^2] = E[E[r_t^2 \mid \mathcal{F}_{t-1}]] = E[\sigma_t^2]$, write $\bar{\sigma}^2 = E[\sigma_t^2]$:
$\bar{\sigma}^2 = \omega + \alpha \, \bar{\sigma}^2 + \beta \, \bar{\sigma}^2 = \omega + (\alpha + \beta) \bar{\sigma}^2$
Solving:
$\boxed{\bar{\sigma}^2 = \frac{\omega}{1 - \alpha - \beta}}$
This is well-defined and positive precisely when $\alpha + \beta < 1$ and $\omega > 0$.
Numerical check: With $\omega = 0.00001$, $\alpha = 0.08$, $\beta = 0.90$ (persistence $= 0.98$), the unconditional variance is $0.00001 / 0.02 = 0.0005$, corresponding to a daily vol of $\sqrt{0.0005} \approx 2.2\%$, or roughly $35\%$ annualized. That is in the right ballpark for a volatile equity index.
---
(c) Shock Propagation and Half-Life
Suppose at time $t-1$ there is a volatility shock -- a large squared return $r_{t-1}^2$ that pushes $\sigma_t^2$ above its long-run mean. How does this excess variance decay?
Iterate the GARCH equation forward. The deviation from the mean $\bar{\sigma}^2$ at horizon $k$ satisfies:
$E_t[\sigma_{t+k}^2] - \bar{\sigma}^2 = (\alpha + \beta)^k \, (\sigma_t^2 - \bar{\sigma}^2)$
So the impact of any shock decays geometrically at rate $(\alpha + \beta)$ per period. This is the persistence parameter.
Propagation mechanism: The shock enters through the $\alpha \, r_{t-1}^2$ term (direct impact). On subsequent days, it is carried forward through the $\beta \, \sigma_{t-1}^2$ term (the "memory" channel). Together, each day the remaining impact is multiplied by $(\alpha + \beta)$.
Half-life: The half-life $h$ is the number of periods until the shock's impact drops to half its initial size:
$(\alpha + \beta)^h = \frac{1}{2}$
$\boxed{h = \frac{\ln 2}{-\ln(\alpha + \beta)} = \frac{\ln 2}{\ln\!\left(\frac{1}{\alpha + \beta}\right)}}$
Example: If $\alpha + \beta = 0.98$, then $h = \ln 2 / \ln(1/0.98) = 0.693 / 0.0202 \approx 34$ trading days, or about 7 weeks. At persistence $0.95$, the half-life drops to about 14 days. At $0.99$ it climbs to 69 days. This is why equity vol feels "sticky" -- typical persistence is in the $0.95$-$0.99$ range.
Answer: (a) Require $\omega > 0$, $\alpha \geq 0$, $\beta \geq 0$, and $\alpha + \beta < 1$. (b) The unconditional variance is $\bar{\sigma}^2 = \omega / (1 - \alpha - \beta)$. (c) Shocks decay geometrically at rate $(\alpha + \beta)^k$, with half-life $h = \ln 2 \,/\, \ln(1/(\alpha + \beta))$.
Intuition
GARCH(1,1) captures two empirical facts about financial returns: volatility clusters (big moves follow big moves) and volatility mean-reverts (extreme vol eventually subsides). The $\alpha$ parameter controls how reactive the model is to new shocks, while $\beta$ controls how much memory the model has. Their sum $\alpha + \beta$ is the single most important number -- it determines how "sticky" volatility is. In practice, fitted persistence is almost always between 0.95 and 0.99, meaning a volatility spike takes weeks to months to fully dissipate.
The unconditional variance $\omega / (1 - \alpha - \beta)$ is the anchor the process reverts to, but when persistence is near 1, that reversion is so slow that vol can spend long stretches far above or below its long-run level. This is why even stationary GARCH models produce the kind of regime-like behavior you see in real markets. When an interviewer asks about GARCH, they want to hear you connect the math ($\alpha + \beta < 1$, geometric decay) to the trading reality (vol clustering, mean reversion speed, and the danger of near-unit-root persistence making your risk estimates stale).