---
title: Sandwiched Volterra Volatility Model
url: https://www.emergentmind.com/topics/sandwiched-volterra-volatility-model
type: topic
---

# Sandwiched Volterra Volatility Model

Sandwiched Volterra Volatility (SVV) model is a stochastic volatility framework in which the asset price is driven by Brownian noise while the volatility factor is driven by a Gaussian Volterra process and constrained by a singular drift to remain between two prescribed deterministic barriers. In the one-factor formulation, the log-price and volatility satisfy
\[
S(t)=e^{X(t)},
\]
\[
X(t)=x_0+rt-\frac12\int_0^t Y^2(s)\,ds +\int_0^t Y(s)\bigl(\rho\,dB_1(s)+\sqrt{1-\rho^2}\,dB_2(s)\bigr),
\]
\[
Y(t)=y_0+\int_0^t b(s,Y(s))\,ds+\int_0^t \mathcal K(t,s)\,dB_1(s),
\]
with \(B_1,B_2\) independent Brownian motions, \(\rho\in(-1,1)\), and \(\mathcal K\) a Volterra kernel [2311.01228]. The class was introduced as a stochastic volatility model driven by an arbitrary Hölder continuous Gaussian Volterra process whose volatility is “sandwiched” between two arbitrary Hölder continuous functions chosen in advance, and it has since been developed for option pricing, implied-volatility asymptotics, and quadratic hedging [2209.10688][2209.13054][2311.01228].

## 1. Model specification and Volterra structure

The defining feature of SVV is that the instantaneous volatility is not postulated as an unconstrained diffusion or rough process, but as a process \(Y\) obtained by adding a singular drift to a Gaussian Volterra noise. In the multi-asset version, for each asset \(i=1,\dots,d\),
\[
S_i(t)=S_i(0)+\int_0^t \mu_i(s)S_i(s)\,ds+\int_0^t Y_i(s)S_i(s)\,dB_i^S(s),
\]
\[
Y_i(t)=Y_i(0)+\int_0^t b_i(s,Y_i(s))\,ds+\int_0^t \mathcal K_i(t,s)\,dB_i^Y(s),
\]
where \(B_i^S\) and \(B_i^Y\) are correlated Brownian motions and \(Z_i(t):=\int_0^t \mathcal K_i(t,s)\,dB_i^Y(s)\) is the Gaussian Volterra noise [2209.10688].

The kernel assumptions are standard Volterra regularity conditions. The kernel is square integrable,
\[
\int_0^T\int_0^T \mathcal K^2(t,s)\,ds\,dt<\infty,
\]
and satisfies a Hölder-type increment condition: there exists \(H\in(0,1)\) such that for every \(\lambda\in(0,H)\),
\[
\int_0^T(\mathcal K(t_2,s)-\mathcal K(t_1,s))^2\,ds \le C_\lambda |t_2-t_1|^\lambda,
\]
or, in the formulation used for hedging,
\[
\int_0^t(\mathcal K(t,u)-\mathcal K(s,u))^2\,du\le \ell_\lambda^2|t-s|^{2\lambda},\qquad 0\le s\le t\le T.
\]
These conditions imply that the Gaussian Volterra process
\[
Z(t)=\int_0^t \mathcal K(t,s)\,dB(s)
\]
has a modification with Hölder continuous paths of any order \(<H\) [2311.01228][2209.13054].

This architecture places SVV within the broader Volterra volatility literature, but with an additional confinement mechanism. General volatility-modulated Lévy-driven Volterra processes are of the form
\[
X(t)=\int_0^t g(t,s)\,\sigma(s)\,dL(s),
\]
or, in stationary form,
\[
X(t)=\int_{-\infty}^t g(t-s)\sigma(s^-)\,dL(s),
\]
with the kernel \(g\) encoding memory and the factor \(\sigma\) encoding stochastic volatility or intermittency [1205.3275][1602.02907]. SVV replaces the free volatility factor by a volatility state \(Y\) that is itself generated by a Volterra kernel and constrained by singular drift.

## 2. Sandwiching drift and boundary confinement

The term “sandwiched” refers to the requirement that volatility remain between two deterministic Hölder-continuous bounds,
\[
\varphi_i(t) < Y_i(t) < \psi_i(t), \qquad t\in[0,T],
\]
with \(0<\varphi_i(t)<\psi_i(t)\) for all \(t\) [2209.10688]. The drift is designed to explode near the boundaries. In the option-pricing formulation, the standing assumptions are
\[
b_i(t,y)\ge \frac{c}{(y-\varphi_i(t))^{\gamma_i}} \quad \text{for } y \text{ near } \varphi_i(t),
\]
\[
b_i(t,y)\le -\frac{c}{(\psi_i(t)-y)^{\gamma_i}} \quad \text{for } y \text{ near } \psi_i(t),
\]
with
\[
\gamma_i > \frac{1}{H_i}-1.
\]
In the one-factor power-law analysis, the analogous exponents satisfy
\[
\gamma_1>\frac1H-1,\qquad \gamma_2>\frac1H-1,
\]
ensuring that the singular repulsion is strong enough relative to the kernel roughness parameter \(H\) [2209.10688][2311.01228].

Under these assumptions, the solution remains trapped in the corridor almost surely:
\[
\varphi(t)<Y(t)<\psi(t)\qquad \text{for all }t\in[0,T],\quad \text{a.s.}
\]
A sharper estimate states that for any \(\lambda\in(0,H)\),
\[
\varphi(t)+\frac{C_Y}{(1+\Lambda)^\beta}\le Y(t)\le \psi(t)-\frac{C_Y}{(1+\Lambda)^\beta},
\]
where \(\Lambda\) is the Hölder modulus of the Volterra process [2311.01228]. In the hedging paper, the corresponding theorem yields existence of a unique strong solution together with explicit pathwise lower and upper margins from the barriers, and finite moments of inverse distances to both boundaries [2209.13054].

This confinement is not a cosmetic addition. Because \(Y_i\) is bounded above and below away from zero, the asset price admits the explicit stochastic exponential representation
\[
S_i(t)=S_i(0)\exp\left\{ \int_0^t\left(\mu_i(s)-\frac{Y_i^2(s)}{2}\right)ds +\int_0^t Y_i(s)\,dB_i^S(s) \right\},
\]
and the model satisfies
\[
\mathbb E\left[\sup_{t\in[0,T]} S_i^r(t)\right]<\infty \qquad \text{for all } r\in\mathbb R
\]
[2209.10688]. A common misconception is therefore to treat SVV as merely another rough-volatility parametrization; the singular sandwiching drift is the structural device that prevents volatility from hitting zero or exploding.

## 3. Position within the Volterra stochastic calculus literature

SVV sits at the intersection of two earlier Volterra traditions: stochastic integration with respect to volatility-modulated Volterra processes, and infinite-dimensional representations of Volterra fields.

For Lévy-driven Volterra processes, the core object is
\[
X(t)=\int_0^t g(t,s)\sigma(s)\,dL(s),
\]
and a Malliavin-calculus-based integration theory defines anticipative integrals \(\int_0^t Y(s)\,dX(s)\) even when \(g\) is singular at the diagonal and the integrand is not adapted [1205.3275]. In the Brownian case, the operator
\[
\mathcal K_g(h)(t,s)=h(s)g(t,s)+\int_s^t \bigl(h(u)-h(s)\bigr)\,g(du,s)
\]
is central, and the integral is written as a Skorohod term plus a Malliavin-derivative correction [1205.3275]. This is directly relevant to nested Volterra structures, because the same paper emphasizes that the ability to integrate random functionals of \(X\) against \(X\) is especially important for models where volatility is itself driven by a Volterra process.

In the white-noise framework, Brownian-driven Volterra integration is extended from \(L^2\)-based calculus to the space \(G^*\) of Potthoff–Timpel distributions, and volatility can be inserted either by pointwise multiplication or by the Wick product,
\[
\int_0^t \Phi(s)\,dX_{\diamond\Sigma}(s) = \int_0^t K_g(\Phi)(t,s)\diamond \Sigma(s)\,\delta B(s) + \int_0^t D_s(K_g(\Phi)(t,s))\diamond \Sigma(s)\,ds.
\]
This provides a mathematically consistent formulation when the volatility is generalized and pointwise multiplication is not defined [1303.4625]. Although that paper does not use the SVV name, it supplies analytic machinery for Volterra kernels with a stochastic factor “sandwiched” inside the integrand.

A different branch of the literature represents Volterra processes as boundary traces of hyperbolic SPDEs. For volatility modulated Volterra processes and Lévy semistationary processes,
\[
X_{t_0}(t)=Y(t,0)
\]
where \(Y\) solves a hyperbolic SPDE driven by the same noise and volatility structure [1602.02907]. This suggests an infinite-dimensional state-space interpretation of layered Volterra volatility systems, especially when one seeks joint simulation of current values and forward curves.

## 4. Pricing, martingale measures, and Malliavin-based valuation

The pricing theory developed for SVV begins from the discounted asset
\[
\widetilde S_i(t):=e^{-\int_0^t \nu(s)\,ds}S_i(t),
\]
which satisfies
\[
d\widetilde S_i(t)=\widetilde \mu_i(t)\widetilde S_i(t)\,dt+Y_i(t)\widetilde S_i(t)\,dB_i^S(t), \qquad \widetilde \mu_i=\mu_i-\nu
\]
[2209.10688]. A central result is a full description of strict local martingale densities \(M\) such that \(M\widetilde S_i\) are local martingales. Choosing the თავისუფર directions \(\xi_l\equiv 0\), Novikov’s condition applies because \(Y_i\) is bounded away from zero, so an equivalent martingale measure exists; allowing nontrivial \(\xi_l\) yields many equivalent local martingale measures, hence market incompleteness [2209.10688].

The same paper proves Malliavin differentiability of both volatility and price. For volatility,
\[
D^j_u Y_i(t) = \ell_{d+i,j}\left( \mathcal K_i(t,u) + \int_u^t \mathcal K_i(s,u)\frac{\partial b_i}{\partial y}(s,Y_i(s)) e^{\int_s^t \frac{\partial b_i}{\partial y}(v,Y_i(v))\,dv}\,ds \right)\mathbf 1_{[0,t]}(u),
\]
and for the log-price \(X_i(t)=\log S_i(t)\), an explicit Malliavin derivative formula is also obtained [2209.10688].

These derivatives are used to price discontinuous claims by Malliavin integration by parts. The discontinuity is transferred from the payoff to a Malliavin weight, and explicit weighted formulas are derived in the one-dimensional case and for basket payoffs [2209.10688]. This is one of the model’s practical differentiators: the sandwiching bounds ensure that the weights involving \(1/Y_i\) remain integrable.

## 5. Short-maturity implied-volatility skew and power laws

A major theoretical development is the proof that SVV can reproduce the empirically observed short-maturity power law for the at-the-money implied-volatility skew. The target asymptotic is
\[
\left|\frac{\partial \widehat \sigma}{\partial \kappa}(\tau,0)\right| = O\!\left(\tau^{-1/2+H}\right),\qquad \tau\to 0,
\]
for some \(H\in(0,1/2)\) [2311.01228].

The analysis proceeds through Malliavin calculus. The first derivative of the volatility is
\[
D_s Y(t) = \mathcal K(t,s) + \int_s^t \mathcal K(u,s)\,b_y'(u,Y(u)) \exp\!\left\{\int_u^t b_y'(v,Y(v))\,dv\right\}du,
\]
and the paper proves the key regularity statement that for any \(t\in[0,T]\) and \(p\ge2\),
\[
Y(t)\in\mathbb D^{2,p}.
\]
An explicit formula is then obtained for \(D_rD_sY(t)\), together with a generalized Malliavin product rule for singular random terms [2311.01228].

The short-time skew constant is controlled by the leading small-time behavior of \(\mathbb E[D_sY(t)]\). Under the kernel conditions
\[
|\mathcal K(t,s)|\le C|t-s|^{-1/2+H},\qquad 0\le s<t\le T,
\]
with
\[
H\in\left(\frac16,\frac12\right),
\]
and the small-time integral asymptotic
\[
\tau^{-3/2-H}\int_0^\tau\int_s^\tau \mathcal K(t,s)\,dt\,ds \to K_Y,
\]
together with \(\rho\neq0\) and \(\rho K_Y<0\), the main theorem yields
\[
\lim_{\tau\to 0}\tau^{\frac12-H}\frac{\partial \widehat \sigma_{\mathrm{ATM}}}{\partial \kappa}
= -\frac{\rho}{\sigma(0)}\,K_Y,
\]
so that the SVV model reproduces the power law
\[
\left|\frac{\partial \widehat\sigma}{\partial \kappa}(\tau,0)\right|
\sim \tau^{-1/2+H}.
\]
The mechanism is transparent: \(D_sY(t)=\mathcal K(t,s)+\) higher-order terms, and a kernel singularity of order \(|t-s|^{-1/2+H}\) induces the same rough-volatility skew exponent [2311.01228].

## 6. Markovian approximation, quadratic hedging, and computation

Because SVV is typically non-Markovian, direct computation of conditional expectations needed for hedging is difficult. The proposed remedy is a finite-dimensional Markovian approximation obtained by replacing \(\mathcal K\) with degenerate kernels
\[
\mathcal K_m(t,s)=\sum_{i=0}^m e_{m,i}(t)f_{m,i}(s)\mathbf 1_{\{s<t\}},
\]
which yields
\[
Z_m(t)=\sum_{i=0}^m e_{m,i}(t)\,U_{m,i}(t), \qquad U_{m,i}(t):=\int_0^t f_{m,i}(s)\,dB_1(s).
\]
The approximated state vector
\[
(X_m,Y_m,U_{m,0},\dots,U_{m,m})
\]
is finite-dimensional and Markovian [2209.13054].

Two approximation families are emphasized. For Hölder convolution kernels, Bernstein polynomials give
\[
\|\mathcal K-\mathcal K_m\|_{L^2([0,T])}\le C\,m^{-H/2}.
\]
For the rough fractional kernel
\[
\mathcal K(t,s)=\frac{(t-s)^{H-\frac12}}{\Gamma(H+\tfrac12)}\mathbf 1_{s<t},
\qquad H\in(0,\tfrac12),
\]
finite sums of exponentials yield
\[
\|\mathcal K-\mathcal K_m\|_{L^2([0,T])}\le C_H\,m^{-4H/5}.
\]
These kernel errors propagate to volatility and prices. In particular, for the discounted prices,
\[
\mathbb E\!\left[\sup_{t\in[0,T]}|X(t)-X_m(t)|^r\right] \le C\|\mathcal K-\mathcal K_m\|^r,
\]
with the right-hand side measured either in \(L^r([0,T]^2)\) or, for convolution kernels, in \(L^2([0,T])\) [2209.13054].

Quadratic hedging is formulated as
\[
\inf_u \mathbb E\left[\left(F-\int_0^T u(s)\,dX(s)\right)^2\right], \qquad F=f(X(T)).
\]
The optimal strategy is given by the non-anticipating derivative. For a partition \(\pi=\{t_k\}\),
\[
u_\pi = \sum_{k=0}^{n-1} \frac{\mathbb E\!\left[(\eta(t_{k+1})-\eta(t_k))\xi\,\big|\,\mathcal G_{t_k}\right]}
{\mathbb E\!\left[(\eta(t_{k+1})-\eta(t_k))^2\,\big|\,\mathcal G_{t_k}\right]}
\mathbf 1_{(t_k,t_{k+1}]},
\]
and
\[
\mathfrak D\xi := \lim_{|\pi|\to0}u_\pi \quad\text{in }L^2(\mathbb P\times[\eta]).
\]
For a European payoff \(F=f(X(T))\),
\[
u_\pi(t_k) = \frac{ \mathbb E\!\left[(X(t_{k+1})-X(t_k))F\,\big|\,\mathcal F_{t_k}\right] }
{ \mathbb E\!\left[(X(t_{k+1})-X(t_k))^2\,\big|\,\mathcal F_{t_k}\right] }.
\]
The Markovian approximation makes these conditional expectations state-dependent functions of \((X_m,Y_m,U_{m,0},\dots,U_{m,m})\), enabling nested Monte Carlo and least-squares Monte Carlo implementations [2209.13054].

The convergence theory extends to hedging. If \(\mathcal K_m\to\mathcal K\), then for general payoffs \(f=f_1+f_2\), with \(f_1\) globally Lipschitz and \(f_2\) of bounded variation,
\[
\mathbb E\!\left[\int_0^T|u_\pi(s)-u_\pi^m(s)|\,ds\right] \le C\,n\sqrt{|\pi|}\,\|\mathcal K-\mathcal K_m\|_{L^2([0,T])}^{1/5},
\]
while for globally Lipschitz payoffs the exponent \(1/5\) is replaced by \(1\) [2209.13054]. This establishes a direct chain from kernel approximation error to hedging error.

In this sense, SVV is simultaneously a model class and a methodological program: a confined Volterra volatility dynamics, a Malliavin-calculus-based pricing theory, a rough-kernel explanation of power-law skews, and a Markovian approximation scheme that restores computational tractability without discarding the underlying non-Markovian structure.

Source: https://www.emergentmind.com/topics/sandwiched-volterra-volatility-model