---
title: Signature Volatility Model
url: https://www.emergentmind.com/topics/signature-volatility-model
type: topic
---

# Signature Volatility Model

Signature volatility models are stochastic-volatility models in which the instantaneous volatility is represented as a linear functional of the signature of a time-extended driving path, most commonly a Brownian motion or a more general primary continuous semimartingale. In the standard Brownian formulation, one sets $\widehat W_t=(t,W_t)$ and defines volatility through iterated integrals of $\widehat W$; in broader formulations, the same construction is applied to a primary process $\widehat X_t=(t,X_t)$ or to an augmented process that also carries the stock-price driver. The class is explicitly non-Markovian at the path level, yet it remains linearly tractable in signature coordinates, which is why the literature emphasizes universality, calibration, Fourier methods, and structural arbitrage conditions rather than a single parametric volatility factor [2503.17103], [2402.01820], [2207.13136].

## 1. Mathematical formulation

The canonical signature volatility model starts from two Brownian motions $(B,W)$ with
\[
d\langle B,W\rangle_t=\rho\,dt,\qquad \rho\in[-1,1],\quad \bar\rho:=\sqrt{1-\rho^2},
\]
and the time-extended path
\[
\widehat W_t=(t,W_t)\in\mathbb R^2.
\]
Its Stratonovich signature is
\[
\Sig(\widehat W)_{0,t}=\Bigl(1,\Sig^1(\widehat W)_{0,t},\Sig^2(\widehat W)_{0,t},\dots\Bigr),
\]
with
\[
\Sig^n(\widehat W)_{0,t}=\int_{0<u_1<\cdots<u_n<t}\circ\,d\widehat W_{u_1}\otimes\cdots\otimes\circ\,d\widehat W_{u_n}.
\]
Given a finite linear functional
\[
\ell\in\Bigl(\bigoplus_{n=0}^N(\mathbb R^2)^{\otimes n}\Bigr)^*,
\]
the instantaneous volatility is
\[
\sigma_t=\bigl\langle \ell,\Sig(\widehat W)_{0,t}\bigr\rangle,
\]
and under the risk-neutral measure the asset price satisfies
\[
\frac{dS_t}{S_t}=\sigma_t\,dB_t,
\qquad
dX_t=-\tfrac12\,\sigma_t^2dt+\sigma_t\,dB_t,\quad X_t=\log S_t,
\]
with $B_t=\rho W_t+\bar\rho W_t^\perp$ and $W^\perp$ independent of $W$ [2503.17103].

A broader formulation replaces the Brownian driver by a general primary process $X=(X^1,\dots,X^d)$ and its time-extension $\widehat X_t=(t,X_t)$. Then a truncated linear signature model takes the form
\[
S_n(\ell)_t=\ell_\emptyset+\sum_{1\le |I|\le n}\ell_I\,\langle e_I,\widehat{\mathbb X}_{0,t}\rangle,
\]
or, in the martingale representation relevant for one-asset pricing,
\[
S_n(\ell)_t=S_0+\int_0^t \sigma_s\,dX_s,\qquad
\sigma_s=\ell^1_\emptyset+\sum_{0<|I|\le n-1}\ell^1_I\,\langle e_I,\widehat{\mathbb X}_{0,s}\rangle.
\]
This makes the instantaneous volatility exactly a linear function of the signature [2207.13136].

The Sig-SDE literature places the same idea inside a controlled SDE:
\[
dX_t=b_t\bigl(\Sig_t^N(X);\theta_b\bigr)\,dt+\sigma_t\bigl(\Sig_t^N(X);\theta_\sigma\bigr)\,dW_t,
\]
with a truncated signature state $\Sig_t^N(X)$. In that setting, the volatility can be parameterized as a linear pre-activation of signature coordinates followed by a positive link function such as $\exp$ or $\softplus$ [2006.00218].

## 2. Tensor-algebra structure and universality

The model’s algebraic backbone is the extended tensor algebra and the shuffle product. In the Brownian formulation of Abi Jaber and Gérard, the signature lives in the extended tensor algebra $\widehat{\mathbb T}(\mathbb R^2)$, with coordinates indexed by words over the alphabet $\{1,2\}$. For $w\in\widehat{\mathbb T}$ one writes
\[
\langle w,\mathbb X_t\rangle=\sum_{n=0}^\infty\sum_{v\in V_n}w^v\,\mathbb X_t^v,
\]
whenever the series converges, and Itô’s formula can be written directly in signature coordinates:
\[
d\langle w_t,\mathbb X_t\rangle
=
\langle\dot w_t+w_t\shuffle e_1+\tfrac12 w_t\shuffle e_{22},\mathbb X_t\rangle\,dt
+
\langle w_t\shuffle e_2,\mathbb X_t\rangle\,dW_t.
\]
This algebraic representation is the source of the model’s tractability in pricing and hedging [2402.01820].

The universality claims in the literature are concrete. Signature volatility models are described as a “remarkably universal” class that includes, but is not limited to, Stein-Stein, Bergomi, and Heston models, together with path-dependent variants such as delayed equations and Volterra models [2402.01820]. In the 2022 calibration framework, the class is called universal in the sense that classical models can be approximated arbitrarily well and the parameters can be learned from all sources of available data by simple methods [2207.13136]. In the rough-path calibration study, the linear signature expansion is justified by the Universal Approximation Theorem for rough paths, and the paper explicitly states that any continuous path-functional $\sigma_t=f((X_s)_{s\le t})$ can be uniformly approximated on compact sets by a linear functional of the truncated signature [2507.23392].

Time augmentation is not merely cosmetic. One explicit statement is that augmenting $(t,X_t)$ guarantees that the signature uniquely encodes both spatial path and its parametrization, preventing aliasing of differently timed paths [2507.23392]. A related fact used in joint SPX–VIX calibration is that when the primary process is a polynomial diffusion, its truncated signature is again a finite-dimensional polynomial diffusion, so conditional expectations reduce to matrix exponentials [2301.13235].

## 3. Martingality, arbitrage, and moment explosions

A central issue is that the discounted price process is always a nonnegative local martingale, but it may fail to be a true martingale. For the finite-order Brownian signature volatility model with leading coefficient $\sigma^{2\cdots2}\neq 0$, the martingality criterion is
\[
S_t\text{ is a true martingale on }[0,T]
\Longleftrightarrow
\begin{cases}
N\text{ is odd},\ \rho\,\sigma^{2\cdots2}\le 0.
\end{cases}
\]
In particular, when $\sigma^{2\cdots2}>0$, one needs $\rho<0$ and odd $N$; equivalently,
\[
S_t\text{ is a martingale}\quad\iff\quad \deg(\ell)\text{ is odd and }\rho<0.
\]
The proof proceeds by a Girsanov change of measure that reduces martingality to non-explosion of the “signature SDE”
\[
X_t=\rho\int_0^t\bigl\langle \ell,\Sig(\widehat X)_{0,s}\bigr\rangle\,ds+W_t,
\qquad \widehat X=(t,X),
\]
and by the asymptotic estimate
\[
\langle\ell,\Sig(\widehat X)_{0,t}\rangle
\sim
\sigma^{2\cdots2}\,\frac{X_t^N}{N!}
\]
for large $|X|$ [2503.17103].

Once martingality is established, higher moments are governed by the correlation boundary
\[
\mathbb E[S_T^m]<\infty
\quad\Longleftrightarrow\quad
|\rho|>\sqrt{1-\tfrac1m},
\qquad m>1,
\]
equivalently
\[
\mathbb E[S_T^m]<\infty
\quad\Longleftrightarrow\quad
\rho^2>1-\tfrac1m.
\]
If $|\rho|<\sqrt{1-\tfrac1m}$ then $\mathbb E[S_T^m]=+\infty$ [2503.17103].

Low-order examples illustrate the theorem sharply. In the first-order case $N=1$,
\[
\sigma_t=\sigma^0+\sigma_1^1\,t+\sigma_2^1\,W_t,
\]
and $S_t$ is always a martingale, with no condition on $\rho$. In the cubic model $N=3$,
\[
\sigma_t=c_0+c_1W_t+c_2\frac{W_t^2}{2!}+c_3\frac{W_t^3}{3!},
\]
with $c_3\neq 0$, one has martingality if and only if $c_3\rho<0$, and moments of order $m$ are finite exactly when $|\rho|>\sqrt{1-1/m}$ [2503.17103].

The broader arbitrage discussion is more structural. In the 2022 signature-based model, absence of arbitrage is tied to the existence of an equivalent local martingale measure under which the drift and quadratic-variation terms in the signature representation disappear, leaving a pure stochastic integral representation [2207.13136]. In the weighted-tensor-algebra formulation of Xodarev, the stochastic exponential
\[
S_t=s_0\exp\Bigl(\int_0^t\xi_s\,dB_s-\tfrac12\int_0^t\xi_s^2\,ds\Bigr),
\qquad
\xi_t=\langle \ell,\widehat W_t\rangle,
\]
is a true, uniformly integrable martingale under the exponential-integrability hypothesis H3, yielding NFLVR, while global solvability of the associated infinite-dimensional Riccati equation is separated from the question of transform non-explosion [2605.17142]. This suggests that “arbitrage-free” and “transform tractable” are related but distinct structural requirements.

## 4. Pricing, transforms, and hedging

Pricing theory for signature volatility models is built on transform methods. In the infinite-dimensional Brownian-signature model, if a time-dependent tensor-valued Riccati equation admits a solution and the relevant exponential process is a true martingale, then one obtains the joint characteristic functional of the log-price and integrated variance. This yields Fourier inversion formulas for European and path-dependent options and supports quadratic hedging formulas through martingale representation [2402.01820].

The pricing results extend beyond plain vanilla options. For European calls, Lewis’ formula is combined with the model characteristic function; for geometric Asian options, one inserts a time-dependent weight into the Riccati system; for quadratic hedging, one differentiates the Fourier integral under the integral sign and identifies the $Z$ and $Z^\perp$ processes in the martingale representation [2402.01820]. In the 2022 theory paper, a “sig-payoff” is any payoff that is linear in the signature of the price path, and its expectation reduces to expectations of finitely many signature coordinates, which are tractable when the primary process is polynomial [2207.13136].

A distinct line of work develops SPX–VIX pricing. Cuchiero, Gazzani, Möller, and Svaluto-Ferro show that when volatility is parametrized by a linear functional $\ell$ of the signature of a polynomial diffusion, one can write
\[
VIX_t^2=(1/\Delta)\,\ell^\top Q(t,\Delta)\,\ell,
\]
so $VIX^2$ is an explicit quadratic form in $\ell$, while $\log S_t$ and $VIX_t^2$ can both be represented as linear functionals of the signature of an augmented process. The calibration task then splits into an offline sampling stage and a standard optimization stage [2301.13235].

The hedging literature uses the same transform structure. For Fourier-invertible claims, semi-explicit quadratic hedges can be derived in a signature volatility model [2508.02759]. Malliavin calculus has also been specialized to finite linear combinations of time-extended Brownian motion signatures, yielding explicit formulas for the Malliavin derivative, Clark–Ocone representation, integration by parts, and Greeks for path-dependent options under signature volatility models [2604.22528].

## 5. Calibration methodologies and reported performance

Calibration exploits the linearity of volatility in signature coordinates, even though prices depend nonlinearly on the coefficients. Time-series calibration in the 2022 framework is a penalized least-squares problem on precomputed signature features, while implied-volatility-surface calibration is phrased as optimization over model prices computed from Monte Carlo samples of signature vectors [2207.13136]. In the Sig-SDE framework, risk-neutral calibration minimizes weighted squared pricing errors with $L^2$ regularization, and real-world calibration may proceed by Euler–Maruyama likelihood or by matching expected signature moments [2006.00218].

The rough-path calibration study uses a linear signature model
\[
\sigma_t=\langle \boldsymbol\beta,\mathrm{Sig}^{\le L}(X)_{0,t}\rangle
\]
and optimizes
\[
\mathcal L(\boldsymbol\beta)=\sum_{i=1}^N w_i\bigl(V_i^{\rm mkt}-V_i(\boldsymbol\beta)\bigr)^2,
\]
typically with $L=2$ or $3$, Monte Carlo precomputation of truncated signatures, and L-BFGS-B for the coefficient update [2507.23392].

The following reported results illustrate the range of calibration settings.

| Study | Setup | Reported result |
|---|---|---|
| [2507.23392] | Uncorrelated Heston, SIG $(L=3)$ vs ASV | RMSE in IV: $1.0\times10^{-4}$ for SIG, $1.2\times10^{-4}$ for ASV; calibration time: $90$ min vs $2$ s |
| [2207.13136] | S&P 500 market data, $n=3$, $M=10^6$ | worst-case error $\sim96$ bp initially; with weighted loss and/or $n=4$, reduced to $34$ bp |
| [2207.13136] | Time-dependent parameters, truncate $n=2$ | uniform absolute error $<12$ bp even at the shortest maturity, $<2$ bp at ATM |
| [2301.13235] | VIX only, $N=3$, 2-dimensional OU primary process | fits six maturities of VIX options to within bid–ask spreads, with relative futures-price errors $\lesssim10^{-3}$ |

These results are heterogeneous rather than directly comparable, because the objectives, data sets, and model classes differ. What is consistent across the cited papers is the computational strategy: expensive signature sampling is typically pushed offline or precomputed pathwise, while the online stage is a lower-dimensional optimization in linear signature coefficients [2301.13235], [2207.13136].

## 6. Extensions, deep architectures, and structural developments

The model class has expanded in two directions: richer representation spaces and sharper structural theory. The deep signature approach for non-Markovian stochastic volatility reformulates the asset dynamics as a rough stochastic differential equation and then represents the rough paths via linear or non-linear combinations of time-extended Brownian motion signatures. After truncation, the pair $(S_t,\widehat{\mathbb W}^{\le N}_t)$ becomes a finite-dimensional Markov diffusion, and the pricing error in the linear-signature case satisfies a factorial-decay bound of order
\[
O\!\Bigl(\frac{C^{N+1}}{(N+1)!}\Bigr)
\]
as $N\to\infty$ [2508.15237].

Machine-learning applications use signatures both as a model class and as a feature extractor. Sig-SDEs are presented as a framework that can be calibrated under both $\mathbb Q$ and $\mathbb P$, simulate future market scenarios, and provide theoretical guarantees for existence, uniqueness, and Euler–Maruyama convergence under global Lipschitz assumptions [2006.00218]. In hedging experiments under non-Markovian stochastic volatility, signature features in feedforward neural networks are reported to outperform LSTMs in most cases, while shallow regression based on a calibrated signature volatility model is described as more accurate and stable than direct linear signature regression across different payoffs and volatility dynamics [2508.02759].

The American-option literature has used signature methods in a different modular configuration. A 2025 hybrid framework combines a rolling R/S estimator for a time-varying Hurst parameter, gradient-boosted forecasting, regime switching between rough Bergomi and Heston simulation, Random Fourier Features for signature kernels, and a least-squares Monte Carlo primal–dual scheme. The reported implementation uses signature depth three and $D=128$, and the paper states that $N=32{,}000$ paths are feasible on a standard CPU [2508.07151].

Structural theory has also advanced. Xodarev’s 2026 treatment formulates the signature SDE on a weighted tensor algebra $T_w$, proves global existence and uniqueness of strong solutions under hypotheses H1–H3, separates NFLVR from Riccati non-explosion, characterizes market completeness through the density of truncated signature spans in $L^2(\mathcal F_T^S,\mathbb Q)$, and derives a Galtchouk–Kunita–Watanabe hedging-error decomposition with residual controlled by a weighted signature tail bound [2605.17142]. A common misconception is that truncation alone guarantees a viable asset-pricing model; the martingale criteria, exponential-integrability conditions, and transform solvability results show that admissibility depends on parity, sign conditions, and integrability, not merely on finite-dimensional approximation [2503.17103], [2605.17142]. Another misconception is that signature volatility models are intrinsically deep-learning constructions; in the core finance literature they are primarily linear functionals on signatures, with deep architectures appearing later as one implementation layer rather than the defining feature [2207.13136], [2006.00218].

Taken together, the literature presents signature volatility models as a unified pathwise framework: linear in signature coordinates, compatible with non-Markovian dynamics, able to recover classical stochastic-volatility families, and now supported by detailed results on martingality, moments, Fourier pricing, Malliavin Greeks, calibration, and completeness [2402.01820], [2503.17103], [2605.17142].

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