Regularity of Local Times Associated with Volterra–Lévy Processes and Path-Wise Regularization of Stochastic Differential Equations

We investigate the space-time regularity of the local time associated with Volterra–Lévy processes, including Volterra processes driven by α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}-stable processes for α∈(0,2]\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \in (0,2]$$\end{document}. We show that the spatial regularity of the local time for Volterra–Lévy process is P\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbb {P}}$$\end{document}-a.s. inverse proportional to the singularity of the associated Volterra kernel. We apply our results to the investigation of path-wise regularizing effects obtained by perturbation of ordinary differential equations by a Volterra–Lévy process which has sufficiently regular local time. Following along the lines of Harang and Perkowski (2020), we show existence, uniqueness and differentiability of the flow associated with such equations.


Introduction
Occupation measures and local times associated with d-dimensional paths ( p t ) t∈ [0,T ] have received much attention over the past decades from both in the analytical and the probabilistic community. The occupation measure essentially quantifies the amount of time the path p spends in a given set, i.e. for a Borel set A ∈ B(R d ) the occupation measure is given by where λ is the Lebesgue measure on R. The local time is given as the Radon-Nikodym derivative of the occupation measure with respect to the Lebesgue measure. The existence of the local time is generally not assured without some further knowledge of the path p, and the existence of the local time associated with the Weierstrass function, and other deterministic fractal like paths, is, to the best of our knowledge, still considered an open question. However, when ( p t ) t∈[0,T ] is a stochastic process, existence of the local time can often be proved using probabilistic techniques, and much research has been devoted to this aim, see, e.g. [14] and the references therein for a comprehensive overview. Knowledge of probabilistic and analytic properties of the local time becomes useful in a variety problems arising in analysis. For example, given a measurable path p with an existing local time, the following formula holds where * denotes convolution, and L : [0, T ] × R d → R + is the local time associated with p. Thus, analytical or probabilistic questions relating to the left hand side integral can often be answered with the knowledge of the probabilistic and analytic properties of the local time L.
In this article, we will study regularity properties of the local time associated with Volterra-Lévy processes given on the form where k(t, ·) ∈ L α ([0, t]) for all t ∈ [0, T ] with α ∈ (0, 2], and L is a Lévy process on a filtered probability space ( , F, P). In the case, when L = B is a Brownian motion, then joint regularity in time and space of the local time associated with Volterra processes has received some attention in recent years as this knowledge can be applied towards regularization of ordinary differential equations (ODEs) by noise [6,12,13,15], as discussed in detail below. Furthermore, in [13], the authors investigated the regularity of the local time associated with α-stable processes, i.e. when the kernel k ≡ 1, and L is an α-stable process. One goal of this article is therefore to extend these results to the general case of Volterra-Lévy processes, as well as apply this to the regularization by noise procedure. Towards this end, we formulate a simple local non-determinism condition for these processes, which will be used to determine the regularity of the local time. The regularity of the local time is then proved in Sobolev space, by application of the recently developed stochastic sewing lemma [19], similarly as done for Gaussian Volterra processes in [15]. By embedding, it follows that the local time is also contained in a wide range of Besov spaces.
As an application of our results on regularity of the local time, we show existence and pathwise uniqueness of stochastic differential equations (SDEs) of the form d dt even when b is a Besov-distribution (the exact regularity requirement of z and b will be given in Sect. (1.1) below). It is well known that certain stochastic processes provide a regularizing effect on SDEs on the form of (1.2). By this, we mean that if the process (z t ) t∈[0,T ] is given on some explicit form, (1.2) might be well posed, even when b does not satisfy the usual assumption of Lipschitz and linear growth. In fact, in [6], the authors show that if z is given as a sample path of a fractional Brownian motion with Hurst index H ∈ (0, 1), Eq. (1.2) is well posed and has a unique solution even when b is only a distribution in the generalized Besov-Hölder space C β with β < 1 2H − 2. More recently, Perkowski and one of the authors of the current article in [15] proved that there exists a certain class of continuous Gaussian processes with exceptional regularization properties. In particular, if z in (1.2) is given as a path of such a process, then a unique solution exists to (1.2) (where the equation is understood in the pathwise sense), for any b ∈ C β with β ∈ R. Moreover, the flow map ξ → x t (ξ ) is infinitely differentiable. We then say that the path z is infinitely regularizing. Not long after this result was published, Galeati and Gubinelli [12], showed that in fact almost all continuous paths are infinitely regularizing by using the concept of prevalence. Furthermore, the regularity assumption on b was proven to be inverse proportional to the irregularity of the continuous process z. In fact, this statement holds in a purely deterministic sense, see, e.g. [12,Thm. 1]. The main ingredient in this approach to regularization by noise is to formulate the ODE/SDE into a nonlinear Young equation, involving a nonlinear Young integral, as was first described in [6]. This reformulation allows one to construct integrals, even in the case when traditional integrals (Riemann, Lebesgue, etc.) do not make sense. A particular advantage of this theory is furthermore that the framework itself does not rely on any probabilistic properties of the processes, such as Markov or martingale properties. This makes this framework particularly suitable when considering SDEs where the additive stochastic process is of a more exotic type. As is demonstrated in the current paper, the framework is well suited to study SDEs driven by Volterra-Lévy processes, which is a class of processes difficult to analyse using traditional probabilistic techniques. We believe that this powerful framework can furthermore be applied towards analysing several interesting problems relating to ill-posed SDEs and ODEs in the future.
Historically, the investigation of similar regularising effects for SDEs with general Lévy noise seems to have received less attention compared to the case when the SDE (1.2) is driven by a continuous Gaussian process. Of course, the general structure of the Lévy noise excludes several techniques which has previously been applied in the Gaussian case. However, much progress has been made also on this front when the equation has jump type noise, and although several interesting results deserve to be mentioned, we will only discuss here some the most recent results and refer the reader to [4,9,10,18,30] for further results. In [22], Priola showed that (1.2) has a path-wise unique strong solution (in a probabilistic sense) when z = L is a symmetric α-stable process with α ∈ (0, 2) and b is a bounded β-Hölder continuous function of order β > 1 − α 2 . In [23], this result was put in the context of path-by-path uniqueness suggested by Davie [7]. More recently, in [8] the authors prove that the martingale problem associated with (1.2) is well posed, even when b is only assumed to be bounded and continuous, in the case when z = L is an α-stable process with α = 1 (being the critical case). Further in [2], the authors show strong existence and uniqueness of (1.2) when z = L is an one-dimensional α-stable process, and b ∈ C β with β > 1 2 − α 2 . Thus, allowing here for possibly distributional coefficients b when α is sufficiently large (i.e. greater than 1). Our results can be seen as an extension of the last result to a purely pathwise setting, and to the case of general Volterra-Lévy processes. Similarly as seen in the Gaussian case, the choice of Volterra kernel then dictates the regularity β ∈ R of the distribution b ∈ C β that can be considered to still obtain existence and uniqueness.

Main Results
We present here the main results to be proven in this article. The first result provides a simple condition to show regularity of the local time associated with Volterra-Lévy processes.

Corollary 2
There exists a class of Volterra-Lévy processes z t = t 0 k(t, s) dL s such that for each t ∈ [0, T ], its associated local time L t is a test function. More precisely, See Example 32, (iv) for proof of this corollary. Inspired by [6,13,15], we apply the result on regularity of the local time to prove regularization of SDEs by Volterra-Lévy noise. Since we will allow the coefficient b in (1.2) to be distributional-valued, it is not a priori clear what we mean by a solution. Indeed, since the integral t 0 b(x s ) ds is not well defined in a Riemann or Lebesgue sense if b is truly distributional, it is not a priori clear how to make sense of (1.2). We therefore begin with the following definition of a solution, which is in line with the definition of pathwise solutions to SDEs used in [6,13,15].

Definition 3
Consider a Volterra-Lévy process z given as in (1.1) with measurable paths, and associated local time L. Let b ∈ S (R d ) be a distribution such that b * L ∈ C γ ([0, T ]; C 2 (R d )) for some γ > 1 2 . Then, for any ξ ∈ R d we say that x is a solution to and θ solves the nonlinear Young equation Here,L t (z) = L t (−z) where L is the local time associated with (z t ) t∈[0,T ] , and the integral is interpreted in the nonlinear Young sense, described in Lemma 38.
where the solution is interpreted in sense of Definition 3. Moreover, if β > n +1−κ for some n ∈ N, then the flow mapping ξ → x t (ξ ) is n-times continuously differentiable.

Structure of the Paper
In Sect. 2, we recall some basic aspects from the theory of occupation measures, local times, and Sobolev/Besov distribution spaces. Section 3 introduces a class of Volterra processes where the driving noise is given as a Lévy process. We show a construction of such processes, even in the case of singular Volterra kernels, and discuss conditions under which the process is continuous in probability. Several examples of Volterra-Lévy processes are given, including a rough fractional α-stable process, with α ∈ [1, 2). In Sect. 4, we provide some sufficient conditions for the characteristics of Volterra-Lévy processes such that their associated local time exists and is P-a.s. contained in a Hölder-Sobolev space of positive regularity. At last, we apply the concept of local times in order to prove regularization by noise for SDEs with additive Volterra-Lévy processes. Here, we apply the framework of nonlinear Young equations and integration, and thus, our results can truly be seen as pathwise, in the "rough path" sense. An appendix is included in the end, where statements and proofs of some auxiliary results are given.

Notation
For a fixed T > 0, we will denote by x t the evaluation of a function at time t ∈ [0, T ], and write x s,t = x t − x s . For some n ∈ N, we define n T := {(s 1 , . . . , s n ) ∈ [0, T ] n | s 1 ≤ · · · ≤ s n }.
To avoid confusion, the letter L will be used to denote a Lévy process, while L will be used to denote the local time of a process. For γ ∈ (0, 1) and a Banach space E, the space C

Occupation Measures and Local Times, and Distributions
This section is devoted to give some background on the theory of occupation measures and local times, as well as definitions of Sobolev and Besov spaces, which will play a central role throughout this article.

Occupation Measure and Local Times
The occupation measure associated with a process (x t ) t∈ [0,T ] gives information about the amount of time the process spends in a given set. Formally, we define the occupation measure μ associated with where λ denotes the Lebesgue measure. The Local time L associated with x is then the Radon-Nikodym derivative with of μ with respect to the Lebesgue measure(as long as this exists). We therefore give the following definition.
Definition 5 Consider a process x : [0, T ] → R d be a process, and let μ denote the occupation measure of x. If there exists a function L : then we say that L is the local time associated with the process (x t ) t∈[0,T ] .

Remark 6
The interpretation of the local time L t (z) is the time spent by the process x : [0, T ] → R d at a given point z ∈ R d . Thus, the study of this object has received much attention from people investigating both probabilistic and path-wise properties of stochastic processes. For purely deterministic processes (x t ) t∈[0,T ] , the local time might still exist; however, as discussed in [15], if x is a Lipschitz path, there exists at least two discontinuities of the mapping z → L t (z). On the other hand, it is well known (see [14]) that the local time associated with the trajectory of a one-dimensional Brownian motion is 1 2 -Hölder regular in its spatial variable (a.s.). More generally, for the trajectory of a fractional Brownian motion with Hurst index H ∈ (0, 1), we know that its local time L is contained in H κ (a.s.) for κ < 1 2H − d 2 , while still preserving Hölder regularity in time. This clearly shows that the more irregular the trajectory of the fractional Brownian motion is, the more regularity we obtain in the local time associated with this trajectory. In this case, the regularity of the local time can therefore be seen as an irregularity condition. This heuristic has recently been formalized in [13]. There, the authors show that if the local time associated with a continuous path (x t ) t∈[0,T ] is regular (i.e. Hölder continuous or better) in space, then x is truly rough, in the sense of [11]. More recently, the authors of [15] showed that the local time associated with trajectories of certain particularly irregular Gaussian processes (for example the log-Brownian motion) is infinitely differentiable in space, and almost Lipschitz in time. In the current article, we will extend this analysis to Lévy processes.
The next proposition will be particularly interesting towards applications in differential equations, and which we will use in subsequent sections.

Proposition 7 (Local time formula) Let b be a measurable function, and suppose (x t ) t∈[0,T ] is a process with associated local time L. Then, the following formula holds for any
A proof of this statement follows directly from the definition of the local time, see [14,Thm. 6.4] for further details.

Remark 8
It is readily seen that, formally, the local time can be expressed in the following way for ξ ∈ R d and (s, where δ is the Dirac distribution.

Remark 9
For future reference, we also recall here that the Dirac distribution δ is contained in the in-homogeneous Sobolev space H − d 2 − for any > 0 (See, e.g. [3, Remark 1.54]).

Besov Spaces and Distributions
Before introducing the notion of Besov spaces, we give a definition of the Paley-Littlewood blocks, which plays a central role in the construction of these spaces.

Remark 12
We will work with regularity of the local time in the Sobolev space H κ . However, towards applications to regularization by noise in SDEs, we will also encounter Besov spaces, through Young's convolution inequality. We therefore give a definition of these spaces here. Of course, through Besov embedding, for any p, q ∈ [2, ∞] and κ ∈ R, (e.g. [3,Prop. 2.20]), and thus, our results imply that the local time is also included in these Besov spaces. We will, however, not specifically work in this setting to avoid extra confusion, but refer the reader to [12,13] for a good overview of regularity of the local time associated with Gaussian processes in such spaces.

Volterra-Lévy Process
In this section, we give a brief introduction on Lévy processes and stochastic integral for a Volterra kernel with respect to a Lévy process. General references for this part are [28,Chp. 4] and [1,Chp. 2,Chp. 4]. In Sect. 3.1, we give the definition of Volterra-Lévy processes (with possibly singular kernels) and obtain the associated characteristic function. Particularly, our framework includes Volterra processes driven by symmetric α-stable noise. In the end, we provide several examples of Volterra-Lévy processes, including the fractional α-stable process.
We begin to provide a definition of Lévy processes, as well as a short discussion on a few important properties.

Definition 13 (Lévy process)
Let T > 0 be fixed. We say that a càdlàg and (F t )adapted stochastic process (L t ) t∈[0,T ] defined on a complete probability space ( , F, (F t ) t∈[0,T ] , P), and which satisfies the usual assumptions is a Lévy process if the following properties hold: (ii) L has independent and stationary increments. (iii) L is continuous in probability, i.e. for all > 0, and all s > 0,

Remark 14
A known description of Lévy process is Lévy-Khintchine formula: for a d-dimensional Lévy process L , the characteristic function ψ of L verifies that for t ≥ 0, there exists a vector a ∈ R d , a positive definite symmetric d × d matrix σ and a Lévy measure ν such that the characteristic function is given by Here, the triple (a, σ, ν) is called the characteristic of the random variable L 1 .
The typical examples for Lévy processes are the case when the Lévy triplet is given by (0, σ, 0), resulting in a Brownian motion. Another typical example is when the characteristic triplet is given by (0, 0, ν) and the Lévy measure ν defines an α-stable process. We provide the following definition for this class of processes.
has the following characteristic function with α ∈ (0, 2] and some positive constant c α , then we say (L t ) t≥0 is a standard α-stable process.
We now move on to the construction of Volterra-Lévy processes, given of the form Of course in the case when (L t ) t∈[0,T ] is a Gaussian process, or even a square integrable martingale, the construction of such a stochastic integral is by now standard, and z is constructed as an element in L 2 ( ) given that k(t, ·) ∈ L 2 ([0, t]) for all t ∈ [0, T ], see, e.g. [24]. However, in the case when L is not square integrable, then the construction of z as a stochastic integral is not as straight forward. However, several articles discuss also this construction in the case of α-stable processes, which would be sufficient for our purpose. The next remark gives only a brief overview on this construction, and we therefore ask the interested reader to consult the given references for further details on the construction.

Remark 16
Consider a symmetric α-stable process L with α ∈ (0, 2 . In particular, in [27,Chp. 3.6] it is shown that the Volterrastable process below is well-defined and exists in L p ( ) for any p < α, given that the kernel In fact, in the case when L is a symmetric α-stable process, it is known that for any 0 where p,α,d means that they differ up to a constant depending on p, α and d (recall that d is the dimension of L ). See, e.g. [26] and the references therein for more details on this relation and the construction of such integrals.
The above discussion yields the following definition of the Volterra-Lévy process.
Definition 17 (Volterra-Lévy process) Fix T > 0, and let (L t ) t∈[0,T ] be a Lévy process as given in Definition 13. For a given kernel k : 2 T → R with the property that where the integral is constructed in L p ( ) sense for p ≤ β, as discussed above. Then, we call the stochastic process (z t ) t∈[0,T ] a Volterra-Lévy process, where L is the associated Lévy process to z and k is called the Volterra kernel.
Based on this definition, it is readily to have the following known property of Volterra-Lévy process, see, e.g. [25].

Proposition 18 Let (L t ) t∈[0,T ] be a Lévy process on a probability space
for any t ∈ [0, T ], then the Volterra-Lévy process (z t ) t∈[0,T ] given by is well defined as an element of L p ( ) for any 0 < p < β. For 0 ≤ s ≤ t ≤ T , the characteristic function of z is given by

4)
and the conditional characteristic function is given by Everything we have introduced so far only relates to the probabilistic properties of Volterra-Lévy process without any details regarding its sample path behaviour. Towards the goal of proving regularity of the local time associated with (z t ) t∈[0,T ] , as done in Sect. 4, we require that the process z is continuous in probability.

Remark 19
Recall that given a sufficiently regular Volterra kernel k, continuity in probability can for example be obtained from the fact L p -convergence entails Pconvergence using, e.g. Markov's inequality. Indeed, suppose (z t ) t∈[0,T ] is a Volterra-Lévy process, as given in Definition 17. It is readily checked that z is continuous in probability if there exists a p > 0 such that This particular fact will be used in Examples 22 and 23 below.
Below we provide three examples of different types of Volterra processes driven by Lévy noise.

Example 20 (Brownian motion)
Suppose L is a Brownian motion with values in R d . Then, it is well known that z t = t 0 k(t, s) dL s is well-defined in L 2 ( ) as a Wiener integral. The sample paths of such processes are clearly measurable, and depending on the regularity of the kernel k, the process may also be (Hölder) continuous.

Example 21 (Square-integrable martingale case)
Then, we know z t = t 0 k(t, s) dL s , t ≥ 0 is well-defined according to Proposition 18 (this is also clear from classical martingale theory, e.g. [1]).
The next example considers the case of α-stable processes and gives explicit conditions on the Volterra kernel k so that the resulting Volterra process is continuous in probability. (Standard α-stable case) Let α ∈ (0, 2), and suppose L is a standard αstable process, as defined in Definition 13. Assume k : 2 T → R is a Volterra kernel such that k(t, ·) ∈ L α ([0, t]) for all t ∈ [0, T ], and that

Example 22
where by assumption H (t, s) → 0 when s → t for all t ∈ [0, T ], and the integrability of k yields out that t s |k(t, r )| α dr → 0 as s → t. Therefore, we conclude by Remark 19 that (z t ) t∈[0,T ] is continuous in probability.
With the above preparation at hand, we can then construct fractional α-stable processes and give a representation of its characteristic function. We summarize this in the following example.
Example 23 (Fractional α-stable process) Let L be an α-stable process with α ∈ (0, 2], and consider the Volterra kernel k(t, s) = (t − s) H − 1 α , H ∈ (0, 1). Then, the process z t = t 0 k(t, s) dL s is called a fractional α-stable process (of Riemann-Liouville type) and specifically if α = 2, then L is a Brownian motion and z is a fractional Brownian motion. Note that in this case k(t, ·) ∈ L α ([0, t], ds) for any H ∈ (0, 1). There is a more detailed study of fractional processes of this type in [27,Chapter 7]. An application of Proposition (3.4) yields that the characteristic function associated with the fractional α-stable process z is given by Furthermore, it is simply checked that the function H defined as in (3.6) with the kernel k(t, s) = |t − s| H − 1 α satisfies H (t, s) → 0 when s → t, and thus by the same argument as used in Example 22, it is readily checked that the fractional α-stable process (z t ) t∈[0,T ] is continuous in probability.

Regularity of the Local Time Associated with Volterra-Lévy Processes
This section is devoted to prove space-time regularity of the local time associated with Volterra-Lévy processes, as defined in Sect. 3. We begin to give a notion of local non-determinism for these processes and provide a few examples of specific processes which satisfy this property.

Local Non-determinism Condition for Volterra-Lévy Process
The following definition of a local non-determinism condition can be seen as an extension of the concept of strong local non-determinism used in the context of Gaussian processes, see, e.g. [13,15,29].

Remark 25
The elementary example of a Volterra kernel is k(t, s) := 1 [0,t] (s) for any 0 ≤ s ≤ t ≤ T . In this case, the Volterra-Lévy process is just given as the Lévy process itself, i.e. z t = L t . If we let L be a standard d-dimensional α-stable process, condition (4.1) fulfils for ζ = 1. Hence, a standard d-dimensional α-stable process is (α, 1)-LND, which coincides with the conclusion in [20, Proposition 4.5, Example (a)].

Remark 26
There already exists several concepts of local non-determinism, but, as far as we know, most of them are given in terms of a condition on the variance of certain stochastic processes. The only exception we are aware of is the definition of Nolan in [20] for α-stable processes, where a similar condition is stated in L p spaces, with p = α (see [20,Definition 3.3]). Of course working with general α-stable processes, we do in general not have finite variance, and thus, the standard definitions of such a concept are not applicable. On the other hand, in the case when α = 2, we have finite variance, and then, the above criterion would be very similar to the condition for strong local non-determinism for Gaussian Volterra processes, as discussed for example in [29]. Working with the conditional characteristic function of Volterra αstable processes, we see, however, that this condition in some sense is what needs to be replaced in order to prove existence and regularity of local times associated with these processes.
It is readily seen that the Volterra α-stable process satisfies (4.1), with ζ depending on the choice of kernel k. The condition is, however, somewhat more general, as we only require the processes to behave similarly to Volterra α-stable processes. Let us provide an example to discuss some interesting process that satisfies the LND condition.

Example 27 (Volterra kernel)
As two examples of Volterra kernel that we are interested most, we give a specific discussion here. The first one usually relates to fractional type processes, for instance, fractional Brownian motion and fractional stable processes. As we will see later, the second one makes the corresponding Volterra-Lévy process an infinitely regularising process, similarly to the Gaussian counterpart discussed in [15]. As we mentioned, the Lévy process L does not have to be Gaussian type processes. For non-Gaussian type L , we mostly consider α-stable processes or the processes which has similar behaviour to stable processes. Since the condition (4.1) only focuses on the characteristic function ψ of L , there is a large class of jump processes which can be studied here.

Example 29 (Stable type processes)
Fix an α ∈ (0, 2). Given a kernel k : 2 T → R with k(t, ·) ∈ L α ([0, t]) and satisfying for some ζ > 0, the following inequality Then, the following list of processes satisfies the LND condition in Definition 24: α-stable processes. In this case, the corresponding characteristic function ψ is given by By Jensen's inequality, it follows that We additionally assume α ∈ (0, 1) (see [17] Example 1.5). This process is not really a stable process, but the small size jumps of this process have similar behaviour to stable processes. Since |ξ | α log(2 + |ξ |) ≥ |ξ | α for ξ ∈ R d , then Therefore, z t = t 0 k(t, r ) dL r , t > 0, is (α, ζ )-LND. The following theorem shows the regularity of the local time associated with Volterra-Lévy processes which is (α, ζ )-LND according to Definition 24.

Regularity of the Local Time
With the concept of local non-determinism at hand, we are now ready to prove the regularity of the local time associated with Volterra-Lévy processes, and thus also proving Theorem 1. The following theorem provides a proof of Theorem 1, as well as giving P-a.s. bounds for the Fourier transform of the occupation measure and the local time.  Proof We will follow along the lines of the proof of [15,Theorem 17], but adapt to the case of Lévy processes. To this end, we will apply the stochastic sewing lemma from [19], which is provided in Lemma 37 for self-containedness.
A Fourier transform of the occupation measure μ s,t (dx) yields t s e i ξ,z r dr . Note that this coincides with the Fourier transform of the local time L s,t (x) whenever L exists. Our first goal is therefore to show that for any p ≥ 2, the following inequality holds for some λ ≥ 0 and γ ∈ ( 1 2 , 1) To this end, the stochastic sewing lemma (see Lemma 37) will provide us with this information. We begin to define

and for a partition P[s, t] of [s, t] define
At last, we will need to control the term δ u A ξ s,t L p ( ) . To this end, using Proposition 18, we know that where E is defined as in (3.5). Therefore, it is readily checked that Of course, moments of the complex exponential E 0,s,r (ξ ) is bounded by 1, i.e. for any r ∈ [s, t], E 0,s,r (ξ ) L p ( ) ≤ 1, and therefore, it follows that Using the fact that z is (α, ζ )-LND for some ζ ∈ (0, 1), and using that (r − s) ζ ≥ (r − u) ζ we obtain the estimate Note in particular that this holds for any p ≥ 2. By the property of the exponential function, we have that for any Since 1 + |ξ | α (1 + |ξ | 2 ) α 2 for all α ∈ (0, 2], applying this relation in (4.3), and assuming that 0 ≤ ηζ < 1 2 , it follows that for all p ≥ 2 there exists a γ > 1 2 and λ < α Thus, both conditions of (A.1) in Lemma 37 are satisfied. By a simple addition and subtraction of the integrand A We will now show that μ s, By Minkowski's inequality, we have and by Jensen's inequality it follows that This implies that Furthermore, for fixed > 0 then e i ξ,z r − e i ξ,z u L p ( ) ≤ (e i ξ,z r − e i ξ,z u )1 |e i ξ,zr −e i ξ,zu |> L p ( ) + ≤ P(|e i ξ,z r − e i ξ,z u | > ) + since e i ξ,z r − e i ξ,z u L q ( ) ≤ 1 for any q. We conclude that for any > 0 (4.5) Since this holds for any partition P, letting the mesh tend to 0, and using the assumption that z is continuous in probability, it follows that e i ξ,z r is continuous in probability (see, e.g. the continuous mapping theorem), and thus Since (4.5) holds for any partition (also for partitions with infinitesimal mesh), we conclude that for any > 0 and since can be chosen arbitrarily small, we conclude that μ s,t (ξ ) = A ξ s,t in L p ( ). We move on to estimate the Sobolev norm L s,t H κ for some appropriate κ ∈ R. We begin to observe that By Minkowski's inequality, it follows that and then use the bound from (4.4) to observe that Choosing κ = λ − d 2 − for some arbitrarily small > 0, it follows that Recalling that λ < α 2ζ , since > 0 could be chosen arbitrarily small, we obtain that for any κ < α 2ζ − d 2 there exists a γ > 1 2 such that Since p ≥ 2 can be chosen arbitrarily large, we conclude by Kolmogorov's theorem of continuity that there exists a set ⊂ of full measure such that for all ω ∈ there exists a C(ω) > 0 such that In particular, this implies that for almost all ω ∈ , μ(ω) ∈ L 2 (R d ) and thus the local time L(ω) (given as the density of μ) exists, and our claim follows.

Remark 31
In the case, when α = 2, then X is a Gaussian process, and Theorem 30 provides the same regularity of the associated local time for a Gaussian Volterra process as proven in for example [15] (or without considering the joint time regularity, as shown in [5,14,21]). This theorem can therefore be seen as an extension of this work to the class of Volterra-Lévy processes.
local time L has compact support, P-a.s. and thus L ∈ C γ ([0, T ]; D(R d )), P-a.s. where D(R d ) denotes the space of test functions on R d . This proves in particular Corollary 2.

Regularization of ODEs Perturbed by Volterra-Lévy Processes
With the knowledge of the spatio-temporal regularity of the local time associated with a Volterra-Lévy process, we can solve additive SDEs with possibly distributionalvalued drift's coefficients. The goal of this section is to prove Theorem 4. To this end, we will recall some of the tools from the theory of nonlinear Young integrals and corresponding equations. This theory for construction of integrals and equations is by now well known (see, e.g. [6,16] and more recently [12,15] for an overview), but for the sake of self-containedness, we have included some short versions of proofs in the appendix. We also mention that conditions for existence and uniqueness of nonlinear Young equations can be stated in more general terms than what is used here. We choose to work with a simple set of conditions to provide a clear picture of the regularising effect in SDEs driven by Volterra-Lévy noise, in contrast to the full generality which could be accessible. More general conditions for existence and uniqueness of nonlinear Young equations can for example be found in [12].
γ , and satisfies the following inequalities for (s, t) ∈ T 2 and ξ,ξ ∈ R d Then, for any ξ ∈ R d there exists a unique solution to the equation Here, the integral is interpreted as the nonlinear Young integral described in for any partition P of [0, t].
Proof See proof in Appendix B.
From here on, all analysis is done pathwise. That is, we now consider a subset ⊂ of full measure such that for all ω ∈ the local time L(ω) associated with a Volterra-Lévy process is contained in C γ T H κ , for γ and κ as given through Theorem 30. With a slight abuse of notation, we will write L = L(ω).
Before moving on to prove existence and uniqueness of ODEs perturbed by Volterra-Lévy processes, we will need a technical proposition on the convolution of the local time with certain (possibly distributional) vector fields.
A combination of Lemma 33 and Proposition 34 provides the existence and uniqueness of ODEs perturbed by (α, ζ )-LND Volterra-Lévy processes. The following corollary and proposition can be seen as proof of Theorem 4.
Corollary 35 (SDEs driven by stable Volterra processes) Let (z t ) t∈[0,T ] be a Volterra-Lévy process which is continuous in probability, and (α, ζ )-LND according to definition 24 for some ζ ∈ (0, α d ) and α ∈ (0, 2]. Suppose b ∈ H β for some β ∈ R such that the following inequality holds β + α 2ζ − d 2 ≥ 2. Then, for any ξ ∈ R d there exists a unique solution y ∈ C γ T (R d ) to the equation Here, the integral and solution is interpreted pathwise in sense of Lemma 33 by setting , and it follows directly that conditions (i)-(iii) of Lemma 33 are satisfied, and thus, a unique solution to (5.4) exists.
Additionally to existence and uniqueness, the authors of [15] provided a general program to prove higher order differentiability of the flow mapping ξ → y t (ξ ). We will here apply this program in order to show differentiability of flows associated with ODEs perturbed by sample paths of a Volterra-Lévy process. It is well known that if b ∈ C k for some k ≥ 1, the flow ξ → y t (ξ ) where y is the solution to the ODE is k-times differentiable. Translating this to the abstract framework of nonlinear Young equations; let y be the solution to where ∈ C γ T C κ for some κ > 1 γ . Then, the flow ξ → y t (ξ ) is κ times differentiable. Recall that in our setting represents the convolution between the (possibly distributional) vector field b and the local time associated with the irregular path of a Volterra-Lévy process. We therefore provide a proposition to highlight the relationship between the regularity of the vector field b, the regularity of the local time associated with a Volterra-Lévy process, and the differentiability of the flow. probability, and (α, ζ )-LND for some α ∈ (0, 2] and ζ ∈ (0, α d ). Suppose b ∈ H β for some β ∈ R such that β + α 2ζ − d 2 ≥ 1 + n for some integer n ≥ 1. Let y(ξ ) ∈ C γ T (R d ) denote the solution to (5.4) starting in ξ ∈ R d . Then, the flow map ξ → y · (ξ ) is n-times Fréchet differentiable.

Proposition 36 Let (z t ) t∈[0,T ] be a Volterra-Lévy process taking values in R d which is continuous in
Proof This result for abstract Young equations was proven in [15, Thm. 2], but we give a short outline of the proof here. Denote by θ = β + α 2ζ − d 2 , and since θ ≥ 2 it follows that there exists a unique solution to (5.4). We will prove the differentiability of the flow by induction and begin to show the existence of the first derivative. It is readily checked that the first derivative of the flow ξ → y · (ξ ) needs to satisfy the equation where the integral is understood in sense of the nonlinear Young integral in Lemma 38,by setting 1 s,t (u s ) = ∇ s,t (y s (ξ ))u s where u s = ∇ y s (ξ ) ∈ R d×d . Since (5.5) is a linear equation, existence and uniqueness can be simply verified following along the lines of the proof of Lemma 33.
Funding Open Access funding enabled and organized by Projekt DEAL.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Appendix B. Nonlinear Young integration and equations
This section is devoted to give the necessary background regarding the nonlinear young integral, and Young equations. We begin to prove the existence of the nonlinear Young integral.
for some γ ∈ (0, 1). Furthermore, suppose y : Then, there exists a unique function I ∈ C γ T satisfying I s,t = I t − I s given by I t := lim |P|→0 I P . We then define t s dr (y r ) = I s,t .
Moreover, we have that δ (y) C γ +η C γ T C 1 y C η , and it follows from [11,Lemma 4 Proof To prove this, we make use of the classical sewing lemma from the theory of rough paths (see [11,Lemma 4.2] and thus invoking the assumption in (B.1), we obtain that |δ u Ξ s,t | |t − u| γ |y s − y u | |t − s| γ +η .
Since ∈ C γ T , and due to the assumption that η + γ > 1, it follows that (B.4) is satisfied, and thus our claim follows from the sewing lemma ( [11,Lemma 4.2]).

Remark 39
Of course, the nonlinear Young integral coincides with the classical Young integral if the abstract integrand s,t (y s ) is for example given by s,t (y s ) = y s X s,t , for some γ -Hölder continuous path x. Furthermore, if b is a measurable function, and z is path of finite p-variation, then set s,t (y s ) = t s b(y s + z r ) dr . In this case, it is readily checked that the integral coincides with the classical Riemann integral See [12] for a comprehensive introduction and discussion of the nonlinear integral.
For self-containedness, we include a proof of Lemma 33. The existence and uniqueness of these equations has been proven in [6,15], and we refer to these references for a full account on these results.

Proof of Lemma 33
This proof follows along the lines of [15,Lemma 30], and thus, we only give here a shorter recollection of the most important details. Let β ∈ ( 1 κ , γ ) where we recall that κ > 1 γ by assumption, and let S T : C Let B T (ξ ) ⊂ C β T (R d ) be a unit ball centered at ξ ∈ R d . In order to prove existence and uniqueness of (5.2), we will begin to show that there exists a τ > 0 such that S τ leaves the unit ball B τ (ξ ) invariant. In the second step, we will show that there exists a τ > 0 such that the solution map S τ is a contraction on the unit ball B τ (ξ ). It then follows by Picard's fixed point theorem that a unique solution exists in the unit ball Bτ (ξ ) forτ = τ ∧ τ . In the end, since ξ → (ξ ) is globally bounded, we can iterate the solution to the intervals [kτ , (k + 1)τ ∧ T ] for k ∈ N. We begin to show the invariance. By application of (B.3) and (i) in (5.1), it follows that for y ∈ B τ (ξ ) Using that y ∈ B τ (ξ ), and thus in particular y C β τ ≤ 1, it follows that S τ (y) C β τ C γ τ C κ τ γ −β .
By choosing τ > 0 sufficiently small, we obtain S τ (y) C β τ ≤ 1, and thus, S τ leaves the unit ball B τ (ξ ) invariant. We continue to prove the contraction property. Applying (B.3), it follows from Lemma 38 that for y, z ∈ B τ (ξ ) we have for some β > 1. We may assume that z 0 = y 0 = ξ . For the first term on the right hand side, it follows by assumption (B.1) (ii) that where we have used that y 0 −z 0 = 0. For the second term on the right hand side of (B.6), we appeal to the proof of the nonlinear Young integral in Lemma 38, and we will need to show that the action of the δ-operator on the integrand Ξ s,t := s,t (y s ) − s,t (z s ) is sufficiently regular and has a contractive property. That is, we will prove that for (s, u, t) ∈ 3 τ , the following inequality holds |δ u Ξ s,t | |t − s| μ y − z C β . By the fundamental theorem of calculus, it follows that Ξ s,t = 1 0 ∇ s,t (ρ y s + (1 − ρ)z s )dρ(y s − z s ).
By addition and subtraction of ∇ u,t (ρ y s + (1 − ρ)z s )(y u − z u ) inside the above integral, invoking (i) of (5.1) and using that κ ≥ 1, we begin to observe that Due to the assumption that β ∈ ( 1 κ , γ ) it follows that β(κ − 1) + γ > 1. Combining (B.8) and (B.9), it follows that for y, z ∈ B τ (ξ ) with z 0 = y 0 = ξ we set β = β(κ − 1) + γ and we have (B.10) Thus, inserting (B.7) and (B.10) into the right hand side of (B.6), we obtain the inequality By choosing τ > 0 small enough, it is clear that the solution map S τ is a contraction on the ball B τ (ξ ). Note in particular that the contraction bound is independent on the initial data, due to the assumption of boundedness of the derivatives of (recall that C γ τ C γ b ([0, τ ]) when γ ∈ (0, 1)). It follows that S τ ∧τ is a contraction and leaves the ball B τ (ξ ) invariant, and it follows by Picard's fixed point theorem that there exists a unique solution to (5.2) on in B τ (ξ ). By standard procedures, one can now iterate the solution to the whole interval [0, T ], and we ask the patient reader to consult [11,Section 8.3] for further details on this part. At last we note that the solution is indeed contained in the space C γ T . Indeed, assume y ∈ C β T satisfies (5.2). Using the inequality in (B.4), the following inequality holds |y s,t | = | t s dr (y r )| | s,t (y s )| + C γ T C κ y C β T |t − s| γ +β y, ,T |t − s| γ , and it follows that y ∈ C γ T . This concludes our proof.