Stochastic Sewing in Analysis

• Stochastic sewing is a probabilistic extension of the classical sewing lemma that reconstructs one-parameter processes from random increments via conditional expectation and martingale techniques.
• It leverages mixed estimates on the defect using Lₘ-control and BDG-type inequalities to lower the regularity threshold below the deterministic requirement.
• The framework underpins innovative applications ranging from Itô integrals and SPDE mild solutions to regularization by noise and law-dependent jump SDE flows.

Searching arXiv for recent and foundational papers on stochastic sewing. Stochastic sewing is a probabilistic extension of the classical sewing lemma that reconstructs a one-parameter process from a two-parameter family of random increments $A_{s,t}$ by proving convergence of Riemann-type sums $\sum_{[s,t]\in\pi}A_{s,t}$ as the mesh $|\pi|\to0$. Its defining feature is that it replaces purely pathwise control of the defect

$\delta A_{s,u,t}:=A_{s,t}-A_{s,u}-A_{u,t}$

by a mixed estimate combining $L_m$-control of the defect with stronger control of its conditional expectation. In this way, stochastic sewing exploits adaptedness, martingale cancellations, and Burkholder–Davis–Gundy-type estimates to lower the effective regularity threshold from the deterministic $>1$ regime to a stochastic threshold in which the random part only needs to be better than order $1/2$ (Lê, 2018). Subsequent work has developed semigroup-adapted mild formulations for SPDEs, Banach-valued variants, extensions based on conditioning at earlier times to capture asymptotic decorrelation, multiparameter versions for random fields, and a Wasserstein-space analogue for flows of probability measures (Li et al., 2021, Lê, 2021, Matsuda et al., 2022, Bechtold et al., 2023, Alfonsi et al., 2024).

1. Foundational formulation and probabilistic mechanism

The basic stochastic sewing problem begins with an adapted two-parameter process $A_{s,t}$, typically $\mathcal F_t$-measurable, and asks when the random Riemann sums

$\sum_i A_{t_i,t_{i+1}}$

converge in $\sum_{[s,t]\in\pi}A_{s,t}$0 or in probability to an increment $\sum_{[s,t]\in\pi}A_{s,t}$1. In the one-parameter stochastic sewing lemma introduced by Lê, the central hypotheses are estimates on the defect $\sum_{[s,t]\in\pi}A_{s,t}$2 and on its predictable projection with respect to $\sum_{[s,t]\in\pi}A_{s,t}$3 (Lê, 2018).

In the form summarized in the lecture notes, one assumes exponents $\sum_{[s,t]\in\pi}A_{s,t}$4, $\sum_{[s,t]\in\pi}A_{s,t}$5 and constants $\sum_{[s,t]\in\pi}A_{s,t}$6 such that

$\sum_{[s,t]\in\pi}A_{s,t}$7

Under these bounds, the sewn increment exists and satisfies

$\sum_{[s,t]\in\pi}A_{s,t}$8

(Butkovsky, 14 Oct 2025).

The probabilistic improvement over deterministic sewing is tied to the decomposition of local errors into a predictable part and a martingale-difference part. Deterministic sewing requires direct summability of the full defect and thus a threshold $\sum_{[s,t]\in\pi}A_{s,t}$9. Stochastic sewing instead estimates the centered random part by BDG, so square-summability suffices; this is the source of the characteristic threshold $|\pi|\to0$0 for the random component (Lê, 2018, Butkovsky, 14 Oct 2025). The lecture notes emphasize that in stochastic problems one often does not need pathwise cancellation of $|\pi|\to0$1; conditional expectation is enough because randomness creates cancellations that are invisible in a deterministic framework (Butkovsky, 14 Oct 2025).

A useful simplification recorded in the notes is that, for the $|\pi|\to0$2-control, it is enough to show

$|\pi|\to0$3

(Butkovsky, 14 Oct 2025). This reflects the fact that the random part can be propagated through the dyadic proof architecture using martingale inequalities rather than direct Hölder summability.

2. Lê’s stochastic sewing lemma and the Doob–Meyer-type structure

Lê’s 2018 paper gives the canonical one-parameter theorem. For an adapted $|\pi|\to0$4-integrable process $|\pi|\to0$5, with $|\pi|\to0$6, suppose there exist $|\pi|\to0$7 and $|\pi|\to0$8 such that for every $|\pi|\to0$9,

$\delta A_{s,u,t}:=A_{s,t}-A_{s,u}-A_{u,t}$0

and

$\delta A_{s,u,t}:=A_{s,t}-A_{s,u}-A_{u,t}$1

Then there exists a unique adapted $\delta A_{s,u,t}:=A_{s,t}-A_{s,u}-A_{u,t}$2-integrable process $\delta A_{s,u,t}:=A_{s,t}-A_{s,u}-A_{u,t}$3 with $\delta A_{s,u,t}:=A_{s,t}-A_{s,u}-A_{u,t}$4 such that $\delta A_{s,u,t}:=A_{s,t}-A_{s,u}-A_{u,t}$5 converges to $\delta A_{s,u,t}:=A_{s,t}-A_{s,u}-A_{u,t}$6 in $\delta A_{s,u,t}:=A_{s,t}-A_{s,u}-A_{u,t}$7 as $\delta A_{s,u,t}:=A_{s,t}-A_{s,u}-A_{u,t}$8, together with increment estimates of the form

$\delta A_{s,u,t}:=A_{s,t}-A_{s,u}-A_{u,t}$9

and

$L_m$0

(Lê, 2018).

The theorem also admits a simplification in which the centered estimate is replaced by the stronger but easier condition

$L_m$1

up to changing constants (Lê, 2018). No continuity of $L_m$2 is assumed; adaptedness is essential.

A distinctive aspect of Lê’s formulation is its Doob–Meyer-type decomposition. Under an additional estimate

$L_m$3

the sewn process decomposes uniquely as

$L_m$4

where $L_m$5 is a martingale and $L_m$6 is an adapted remainder satisfying quantitative increment estimates (Lê, 2018). This is not a finite-variation decomposition in the classical sense; rather, it is a sewing-based martingale/predictable splitting. The same paper shows that the corresponding partial Riemann sums

$L_m$7

converge in $L_m$8 to $L_m$9 and $>1$0 (Lê, 2018).

This framework subsumes standard objects from Itô calculus. For $>1$1, the stochastic sewing lemma yields the Itô integral as the limit of left-point Riemann sums. For $>1$2, it constructs quadratic variation as a sewn limit. It also recovers Itô’s formula by sewing the Taylor expansion increments of $>1$3 (Lê, 2018). The relevance of stochastic sewing is therefore not confined to abstract convergence theory; it supplies an increment-based architecture in which Itô integration, quadratic variation, and stochastic Taylor expansions appear as instances of a single principle.

3. Structural generalizations: Banach, mild, multiparameter, and Wasserstein formulations

Subsequent developments preserve the core idea while changing the ambient geometry.

Banach-space stochastic sewing

The Banach-valued extension replaces Hilbert or Euclidean target spaces by a separable Banach space $>1$4 with martingale type $>1$5. In this setting the stochastic threshold changes from $>1$6 to $>1$7. For a measurable $>1$8 with $>1$9 $1/2$0-measurable, one assumes

$1/2$1

and

$1/2$2

for a control $1/2$3. Then the Riemann sums converge to an adapted process $1/2$4, with quantitative analogues of the scalar estimates (Lê, 2021). In Hilbert spaces, martingale type $1/2$5 recovers the classical exponent $1/2$6. This extension makes stochastic sewing available in Besov, Sobolev, and Triebel–Lizorkin settings (Lê, 2021).

Mild stochastic sewing

For semigroup-driven evolution equations, the natural increment is not $1/2$7 but

$1/2$8

where $1/2$9 is a semigroup generated by an operator $A_{s,t}$0. The mild stochastic sewing lemma adapts the coboundary to

$A_{s,t}$1

and proves convergence of

$A_{s,t}$2

under $A_{s,t}$3-size and conditional-defect estimates in semigroup interpolation spaces (Li et al., 2021). A distinctive feature is the quantitative estimate

$A_{s,t}$4

which trades time regularity for spatial regularity through semigroup smoothing (Li et al., 2021). The authors explicitly note that this quantitative bound does not follow directly from earlier stochastic sewing results because the mild formulation depends essentially on interpolation and semigroup structure (Li et al., 2021).

Multiparameter stochastic sewing

In several parameters, one no longer controls a single defect. The multiparameter lemma for random fields on rectangles in $A_{s,t}$5 requires control of all “partial defects” obtained by cutting in subsets of coordinate directions. The theorem is formulated relative to a commuting multiparameter filtration and a function space $A_{s,t}$6. If $A_{s,t}$7, with $A_{s,t}$8, $A_{s,t}$9, $\mathcal F_t$0, $\mathcal F_t$1, and $\mathcal F_t$2, then the grid Riemann sums converge to sewing operators $\mathcal F_t$3 indexed by subsets $\mathcal F_t$4, uniquely characterized by additivity in the directions of $\mathcal F_t$5 and an alternating-remainder estimate (Bechtold et al., 2023).

The conceptual novelty is that multiple directions force a lattice of defects and conditional expectations. This distinguishes stochastic sewing in several parameters from deterministic hypercube sewing, such as the earlier multi-parameter deterministic theory on hypercubes (Harang, 2018). A plausible implication is that multiparameter stochastic sewing should be viewed not merely as a higher-dimensional reformulation but as a genuinely new algebra of conditional cancellations.

Wasserstein-space stochastic sewing

A further generalization lifts sewing from random variables to maps on probability measures

$\mathcal F_t$6

Here there is no intrinsic martingale structure at the level of measures, so the construction proceeds through a probabilistic representation by coupling operators $\mathcal F_t$7. Under a $\mathcal F_t$8-coupling family satisfying representation, growth, and stochastic sewing estimates, one obtains a flow $\mathcal F_t$9 on $\sum_i A_{t_i,t_{i+1}}$0 with

$\sum_i A_{t_i,t_{i+1}}$1

together with convergence rates of partition schemes in $\sum_i A_{t_i,t_{i+1}}$2 (Alfonsi et al., 2024). The paper emphasizes that the genuinely new regime is $\sum_i A_{t_i,t_{i+1}}$3, where classical sewing on Wasserstein space is not sufficient but the stochastic decomposition still yields convergence (Alfonsi et al., 2024).

4. Earlier-time conditioning and asymptotic decorrelation

A major extension is the 2022 result that replaces conditioning on the present $\sum_i A_{t_i,t_{i+1}}$4 by conditioning on a strictly earlier $\sum_i A_{t_i,t_{i+1}}$5, $\sum_i A_{t_i,t_{i+1}}$6, in order to exploit asymptotic decorrelation. The motivating observation is that some fractional Brownian increments become weakly correlated with the distant past in a way that is invisible if one conditions only at time $\sum_i A_{t_i,t_{i+1}}$7 (Matsuda et al., 2022).

The generalized theorem assumes $\sum_i A_{t_i,t_{i+1}}$8 is $\sum_i A_{t_i,t_{i+1}}$9-measurable, $\sum_{[s,t]\in\pi}A_{s,t}$00, and that for every $\sum_{[s,t]\in\pi}A_{s,t}$01,

$\sum_{[s,t]\in\pi}A_{s,t}$02

and

$\sum_{[s,t]\in\pi}A_{s,t}$03

Under

$\sum_{[s,t]\in\pi}A_{s,t}$04

there exists a unique adapted process $\sum_{[s,t]\in\pi}A_{s,t}$05 such that the Riemann sums converge in $\sum_{[s,t]\in\pi}A_{s,t}$06, with conditional and moment sewing estimates

$\sum_{[s,t]\in\pi}A_{s,t}$07

and

$\sum_{[s,t]\in\pi}A_{s,t}$08

(Matsuda et al., 2022).

The singular factor $\sum_{[s,t]\in\pi}A_{s,t}$09 is the key new ingredient. It quantifies dependence on the past and turns asymptotic decorrelation into a sewing estimate. The theorem also yields a convergence rate

$\sum_{[s,t]\in\pi}A_{s,t}$10

for some $\sum_{[s,t]\in\pi}A_{s,t}$11 (Matsuda et al., 2022). A companion corollary allows singularities at the origin, with factors such as $\sum_{[s,t]\in\pi}A_{s,t}$12 and $\sum_{[s,t]\in\pi}A_{s,t}$13, which is needed in local-time applications (Matsuda et al., 2022).

The proof retains the standard dyadic/BDG architecture but estimates the predictable part by iteratively reapplying the conditional estimate across dyadic levels, rather than by a direct triangle inequality. This iterative use of earlier-time conditioning is precisely what removes an extra technical restriction that would otherwise appear (Matsuda et al., 2022).

For fractional Brownian motion $\sum_{[s,t]\in\pi}A_{s,t}$14, represented through Mandelbrot–Van Ness,

$\sum_{[s,t]\in\pi}A_{s,t}$15

the paper derives estimates such as

$\sum_{[s,t]\in\pi}A_{s,t}$16

for $\sum_{[s,t]\in\pi}A_{s,t}$17 (Matsuda et al., 2022). This is the template instance of a defect that fails the original $\sum_{[s,t]\in\pi}A_{s,t}$18-condition but becomes tractable once conditioned on a remote past $\sum_{[s,t]\in\pi}A_{s,t}$19. A plausible interpretation is that this result enlarges stochastic sewing from a local-adapted principle to a framework capable of encoding long-memory Gaussian structure.

5. Principal applications

The applications of stochastic sewing are broad but structurally unified: one identifies a local increment whose defect has a small predictable component and a square-summable random component, then uses the sewn limit as an integral, additive functional, or effective flow.

Stochastic integrals and rough Gaussian drivers

For Brownian motion, stochastic sewing reconstructs the Itô integral from left-point Riemann sums and yields quadratic variation and Itô’s formula (Lê, 2018). The same technique extends to fractional Brownian motion in regimes where deterministic sewing is insufficient. In the 2022 extension, if $\sum_{[s,t]\in\pi}A_{s,t}$20 and $\sum_{[s,t]\in\pi}A_{s,t}$21 is bounded measurable, then for

$\sum_{[s,t]\in\pi}A_{s,t}$22

the theorem applies with

$\sum_{[s,t]\in\pi}A_{s,t}$23

and the sums

$\sum_{[s,t]\in\pi}A_{s,t}$24

converge in $\sum_{[s,t]\in\pi}A_{s,t}$25 for every $\sum_{[s,t]\in\pi}A_{s,t}$26, with

$\sum_{[s,t]\in\pi}A_{s,t}$27

(Matsuda et al., 2022). The significance is that no Hölder smoothness of $\sum_{[s,t]\in\pi}A_{s,t}$28 is required.

For Stratonovich-type sums along fBM with $\sum_{[s,t]\in\pi}A_{s,t}$29, the same paper considers

$\sum_{[s,t]\in\pi}A_{s,t}$30

proves convergence of the midpoint/trapezoidal sums in $\sum_{[s,t]\in\pi}A_{s,t}$31, and derives

$\sum_{[s,t]\in\pi}A_{s,t}$32

under $\sum_{[s,t]\in\pi}A_{s,t}$33, plus the symmetry condition $\sum_{[s,t]\in\pi}A_{s,t}$34 when $\sum_{[s,t]\in\pi}A_{s,t}$35 and $\sum_{[s,t]\in\pi}A_{s,t}$36 (Matsuda et al., 2022).

Additive functionals and regularization by noise

Lê’s original paper applies stochastic sewing to additive functionals of Markov processes and SDEs with irregular drift, including Brownian and fractional Brownian settings (Lê, 2018). The lecture notes revisit this theme as one of the central applications. For Brownian motion, they state the averaging estimate

$\sum_{[s,t]\in\pi}A_{s,t}$37

obtained by choosing the germ

$\sum_{[s,t]\in\pi}A_{s,t}$38

and applying stochastic sewing together with heat-semigroup smoothing (Butkovsky, 14 Oct 2025). For fractional Brownian motion, the notes record

$\sum_{[s,t]\in\pi}A_{s,t}$39

for $\sum_{[s,t]\in\pi}A_{s,t}$40, linking stochastic sewing to Catellier–Gubinelli-type regularization by noise (Butkovsky, 14 Oct 2025).

The Banach-valued theorem is designed precisely for singular additive functionals. For

$\sum_{[s,t]\in\pi}A_{s,t}$41

with $\sum_{[s,t]\in\pi}A_{s,t}$42 distributional in space, the paper proves that for suitable $\sum_{[s,t]\in\pi}A_{s,t}$43,

$\sum_{[s,t]\in\pi}A_{s,t}$44

and gives the estimate

$\sum_{[s,t]\in\pi}A_{s,t}$45

(Lê, 2021). This yields a regularization-by-noise principle in Besov scales.

Local times

The 2022 fractional stochastic calculus paper gives a discretization formula for the local time $\sum_{[s,t]\in\pi}A_{s,t}$46 of one-dimensional fBM. For

$\sum_{[s,t]\in\pi}A_{s,t}$47

the Riemann sums converge in $\sum_{[s,t]\in\pi}A_{s,t}$48 to

$\sum_{[s,t]\in\pi}A_{s,t}$49

(Matsuda et al., 2022). It also provides an integrated convergence statement in $\sum_{[s,t]\in\pi}A_{s,t}$50 and a crossing-count representation of local time (Matsuda et al., 2022). In several parameters, the multiparameter stochastic sewing lemma is used to control Fourier transforms of occupation measures of Gaussian sheets, leading to Sobolev/Bessel regularity of local times and then to regularization-by-noise results for planar SDEs (Bechtold et al., 2023).

SPDEs and averaging

The mild stochastic sewing lemma is applied to SPDEs

$\sum_{[s,t]\in\pi}A_{s,t}$51

with $\sum_{[s,t]\in\pi}A_{s,t}$52 generating a bounded analytic semigroup, $\sum_{[s,t]\in\pi}A_{s,t}$53 a trace-class fBM with $\sum_{[s,t]\in\pi}A_{s,t}$54, and $\sum_{[s,t]\in\pi}A_{s,t}$55 a fast random environment (Li et al., 2021). For the stochastic convolution, the paper studies

$\sum_{[s,t]\in\pi}A_{s,t}$56

and proves an $\sum_{[s,t]\in\pi}A_{s,t}$57-estimate for the mild integral

$\sum_{[s,t]\in\pi}A_{s,t}$58

(Li et al., 2021). Together with the corresponding drift estimate, this yields uniform $\sum_{[s,t]\in\pi}A_{s,t}$59-bounds, stability under perturbation, and a fractional averaging principle for SPDEs in random environments (Li et al., 2021).

Law-dependent jump SDEs

The Wasserstein-space theorem is applied to a McKean–Vlasov jump SDE on a Hilbert space, yielding a flow $\sum_{[s,t]\in\pi}A_{s,t}$60 on $\sum_{[s,t]\in\pi}A_{s,t}$61, weak existence of solutions, and uniqueness of marginal laws among solutions with $\sum_{[s,t]\in\pi}A_{s,t}$62-continuous marginals (Alfonsi et al., 2024). The one-step Euler map satisfies the abstract $\sum_{[s,t]\in\pi}A_{s,t}$63-sewing property with

$\sum_{[s,t]\in\pi}A_{s,t}$64

which is sufficient for the stochastic-sewing construction of the law flow (Alfonsi et al., 2024).

6. Proof architecture, related sewing paradigms, and conceptual scope

Across formulations, the proof strategy is remarkably stable. The standard ingredients are dyadic refinement, decomposition of partition errors into localized defect terms, a martingale inequality of BDG or Burkholder type for the centered component, and direct control of the predictable part by conditional expectations (Lê, 2018, Li et al., 2021, Matsuda et al., 2022, Butkovsky, 14 Oct 2025). In the Banach-valued theory, BDG is replaced by martingale type $\sum_{[s,t]\in\pi}A_{s,t}$65, changing the random threshold from $\sum_{[s,t]\in\pi}A_{s,t}$66 to $\sum_{[s,t]\in\pi}A_{s,t}$67 (Lê, 2021). In several parameters, repeated one-dimensional BDG arguments require a commuting multiparameter filtration (Bechtold et al., 2023). In the mild theory, semigroup smoothing provides the additional “seesaw” gain from time regularity to spatial regularity (Li et al., 2021). In the earlier-time extension, the predictable estimate is iterated across dyadic levels to exploit the singular factor $\sum_{[s,t]\in\pi}A_{s,t}$68 (Matsuda et al., 2022).

The broader sewing landscape includes deterministic hypercube sewing and multiplicative sewing. The hypercube theory extends the deterministic lemma to fields on $\sum_{[s,t]\in\pi}A_{s,t}$69, leading to multi-parameter Young integration and hyperbolic equations driven by Hölder fields (Harang, 2018). The multiplicative sewing lemma, by contrast, constructs an evolution system from an almost multiplicative family of operators $\sum_{[s,t]\in\pi}A_{s,t}$70, with quantitative weak-approximation applications to SDEs (Hocquet et al., 2022). These theories are adjacent rather than identical. Stochastic sewing is additive at the level of increments, but its proof and applications often interface with semigroup or operator structures that resemble multiplicative sewing. This suggests a common algebraic motif: local consistency plus a summability mechanism produce a global object. What differs is the summability mechanism—pathwise Hölder regularity in deterministic sewing, stochastic cancellation in stochastic sewing, and near-multiplicativity in multiplicative sewing.

A recurrent misconception is that stochastic sewing is merely a probabilistic proof of the deterministic sewing lemma. The literature indicates a stronger claim. The method is specifically designed for situations in which deterministic sewing fails because the defect is not regular enough pathwise, yet the conditional expectation or martingale structure is sufficiently regular to force convergence (Lê, 2018, Butkovsky, 14 Oct 2025). Another misconception is that the framework is limited to semimartingales. The applications to fractional Brownian motion, Gaussian sheets, SPDEs with semigroup structure, and Wasserstein flows show that stochastic sewing is effective precisely in non-semimartingale and non-Markovian settings where conditional cancellations remain exploitable (Matsuda et al., 2022, Li et al., 2021, Bechtold et al., 2023, Alfonsi et al., 2024).

In that sense, stochastic sewing functions as a unifying technique for converting local conditional smoothing into global integral, flow, or averaging estimates. The lecture notes place under this umbrella regularization by noise for Brownian and fractional Brownian SDEs, well-posedness of SPDEs with irregular drift, averaging operators and local times, and numerical algorithms such as Euler schemes (Butkovsky, 14 Oct 2025). A plausible synthesis is that stochastic sewing has evolved from a lemma about convergence of random Riemann sums into a general methodology for stochastic analysis whenever local increments are irregular pathwise but regular after conditioning.