## Some sample path properties of multifractional Brownian motion

## Paul Balança

e-mail: paul.balanca@gmail.com

École Centrale Paris Laboratoire MAS, ECP Grande Voie des Vignes 92295 Châtenay-Malabry, France

Abstract: The geometry of the multifractional Brownian motion (mBm) is known to present a complex and surprising form when the Hurst function is greatly irregular. Nevertheless, most of the literature devoted to the subject considers sufficiently smooth cases which lead to sample paths locally similar to a fractional Brownian motion (fBm). The main goal of this paper is therefore to extend these results to a more general frame and consider any type of continuous Hurst function. More specifically, we mainly focus on obtaining a complete characterization of the pointwise Hölder regularity of the sample paths, and the Box and Hausdorff dimensions of the graph. These results, which are somehow unusual for a Gaussian process, are illustrated by several examples, presenting in this way different aspects of the geometry of the mBm with irregular Hurst functions.

AMS 2000 subject classifications: 60G07, 60G17, 60G22, 60G44.

Keywords and phrases: 2-microlocal analysis, Box dimension, Hausdorff dimension, Hölder regularity, multifractional Brownian motion.

## 1. Introduction

↦

The multifractional Brownian motion (mBm) has been independently introduced by Peltier and Lévy Véhel [29] and Benassi et al. [6] as a natural extension of the well-known fractional Brownian motion (fBm). The main idea behind these two works was to drop the stationary assumption on the process, and allow the Hurst exponent to change as time passes. In this way, the mBm is parametrized by a function t → H t ( ), usually continuous, and happens to be an interesting stochastic model for non-stationary phenomena (e.g. signal processing, traffic on internet, terrain modelling, . . . ).

Since its introduction, several authors have investigated sample path properties of the mBm. For instance, Benassi et al. [6] and Peltier and Lévy Véhel [29] have respectively studied the law of iterated logarithm and the Hölder regularity of its trajectories. In the latter, the Box and Hausdorff dimensions of the graph have also been determined. Moreover, the fine covariance structure have been analysed and refined by Ayache et al. [3], Herbin [15] and Stoev and Taqqu [35]. Finally, the local time and the Hausdorff dimension of the level sets have recently been considered by Boufoussi et al. [10]. Several multiparameter extensions, which generalise the classic Lévy fractional Brownian motion and fractional Brownian sheet, have also been introduced and investigated by Herbin [15] Meerschaert et al. [24] and Ayache et al. [4].

Fine geometric properties obtained in the aforementioned works usually rely on a key as-

sumption H 0 on the Hurst function:

$$H \text{ is a } \beta \text{-Holder continuous function such that } \sup _ { t \in R } H ( t ) < \beta. \quad \quad ( \mathcal { H } _ { 0 } )$$

Owing to this hypothesis, the study of the H 0 -mBm is usually easier and leads to results closely related to their counterparts on fractional Brownian motion. This structure of the mBm is due to the form of its covariance which, under H 0 , is locally equivalent to fBm's one with Hurst exponent H t ( ):

$$\forall u, v \in B ( t, \rho ) ; \ \mathbb { E } \left [ ( X _ { u } - X _ { v } ) ^ { 2 } \right ] \asymp | u - v | ^ { 2 H ( t ) },$$

for any t ∈ R \ { 0 } .

The study of the more general multifractional Brownian motion, i.e. when the assumption H 0 does not hold, has only been recently considered by Herbin [15] and Ayache [1]. Both have investigated the Hölder regularity of the sample paths, and more precisely the pointwise and local Hölder exponents, which are defined by

$$\forall t \in { \mathbf R } ; \quad \alpha _ { X, t } = \sup \left \{ \alpha \, \colon \lim \sup _ { \rho \to 0 } \sup _ { u, v \in B ( t, \rho ) } \frac { | X _ { u } - X _ { v } | } { \rho ^ { \alpha } } < \infty \right \} \quad \quad ( 1. 1 )$$

and

$$\forall t \in { \mathbf R } ; \quad \widetilde { \alpha } _ { X, t } = \sup \left \{ \alpha \, \colon \lim \sup _ { \rho \to 0 } \sup _ { u, v \in B ( t, \rho ) } \frac { | X _ { u } - X _ { v } | } { | u - v | ^ { \alpha } } < \infty \right \}.$$

These two exponents aim to characterise the local Hölder behaviour of the trajectories at t , by respectively comparing its oscillations to a power of ρ or | u -v | . For instance, it is wellknown that a fractional Brownian motion B H satisfies with probability one and for all t ∈ R , α B H ,t = ˜ α B H ,t = H .

The general form of the Hölder regularity of the multifractional Brownian motion has first been obtained by Herbin [15], proving that for all t &gt; 0

$$\alpha _ { X, t } = H ( t ) \wedge \alpha _ { H, t } \quad \text{and} \quad \widetilde { \alpha } _ { X, t } = H ( t ) \wedge \widetilde { \alpha } _ { H, t } \quad \text{a.s.}$$

The previous result only holds for a fixed t &gt; 0, and thus not uniformly on the sample paths. We also note that the geometric properties of a generic mBm are not a simple extension of those on fBm, but display a more complicated structure where the fine geometry of the Hurst function intervene. Ayache [1] has recently extended this study, obtaining a uniform characterisation of the local Hölder exponent and proving that its pointwise exponent can behave randomly as time passes. The resulting process, with unusual sample path properties, is sometimes called the irregular multifractional Brownian motion in the literature.

The main goal of this work is therefore to push further the study of this irregular mBm. In Section 2, we first discuss the existence of an alternative deterministic representation of the fractional Brownian field, which will then be used to study more precisely the geometric properties of the mBm. Section 3 is devoted to the 2-microlocal and Hölder regularity of the trajectories (Theorem 3.1 and Proposition 3.2), and therefore extends the results obtained by Herbin [15] and Ayache [1]. Hausdorff and Box dimensions of the graph, and images of fractal sets, are investigated in Section 4 (Theorems 4.1, 4.2 and 4.3), where further connections with

the geometry of the Hurst function are presented. Appendix 5 gathers a few deterministic results related to the 2-microlocal frontier and the fractal dimensions which are used along this work. Finally, several examples are also given in these different sections to illustrate the main results and some particular aspects of the geometry of the irregular multifractional Brownian motion.

## 2. Deterministic representation

The question of the equivalence of the different mBm representations is known to be non-trivial. On one side, Peltier and Lévy Véhel [29] have introduced the mBm as the following stochastic integral

$$\forall t \in \mathbf R ; \ \ X _ { t } = \frac { 1 } { \Gamma ( H ( t ) + \frac { 1 } { 2 } ) } \int _ { \mathbf R } \left [ ( t - u ) _ { + } ^ { H ( t ) - 1 / 2 } - ( - u ) _ { + } ^ { H ( t ) - 1 / 2 } \right ] \mathrm d W _ { u }.$$

Whereas on the other hand, Benassi et al. [6] have considered a real-harmonizable definition:

$$\forall t \in \mathbf R ; \ \widetilde { X } _ { t } = \int _ { \mathbf R } \frac { e ^ { i u t } - 1 } { | u | ^ { H ( t ) + 1 / 2 } } \, \mathrm d \widehat { W } _ { u }.$$

Stoev and Taqqu [35] have proved that the covariance structure is slightly different between these two definitions. Considering this fact, they have suggested to name multifractional Brownian motion any process X which has the following form

$$\forall t \in \mathbf R ; \ \ X _ { t } = a ^ { + } B ^ { + } ( t, H ( t ) ) + a ^ { - } B ^ { - } ( t, H ( t ) ),$$

↦

where ( a + , a -) ∈ R 2 \ { 0 0 , } and the fractional Brownian fields ( t, H ) → B ± ( t, H ) are defined by

$$B ^ { \pm } ( t, H ) = \frac { 1 } { \Gamma ( H + \frac { 1 } { 2 } ) } \int _ { \mathbf R } \left [ ( t - u ) _ { \pm } ^ { H - 1 / 2 } - ( - u ) _ { \pm } ^ { H - 1 / 2 } \right ] \mathrm d W _ { u },$$

for all t ∈ R and H ∈ (0 , 1).

In the next proposition, we present an alternative deterministic representation for these fractional Brownian fields which will be very useful for the study of sample path properties.

Proposition 2.1. Suppose ( B t ) t ∈ R is a Brownian motion such that B t := W ([0 , t ]) , where W is a Wiener measure. Then, for all t ∈ R and any H ∈ (0 1) , , the fractional Brownian field B ± ( t, H ) is almost surely equal to:

- · if ( t, H ) ∈ R × ( 1 2 , 1 ) ,

$$B ^ { \pm } ( t, H ) = \frac { \pm 1 } { \Gamma ( H - \frac { 1 } { 2 } ) } \int _ { \mathbf R } B _ { u } \left [ ( t - u ) _ { \pm } ^ { H - 3 / 2 } - ( - u ) _ { \pm } ^ { H - 3 / 2 } \right ] \mathrm d u. \quad \quad ( 2. 2 )$$

- · if ( t, H ) ∈ R × ( 0 , 1 2 ] and T is a fixed number such that T ∈ ±∞ ( , t ) ,

$$B ^ { \pm } ( t, H ) & = \frac { \pm ( H - \frac { 1 } { 2 } ) } { \Gamma ( H + \frac { 1 } { 2 } ) } \int _ { \mathbf R } ( B _ { u } - B _ { t } \mathbf 1 _ { \{ t \geq \pm T \} } ) ( t - u ) _ { \pm } ^ { H - 3 / 2 } - B _ { u } ( - u ) _ { \pm } ^ { H - 3 / 2 } \mathbf d u \\ & \pm B _ { t } ( t - T ) _ { \pm } ^ { H - 1 / 2 },$$

or equivalently,

$$\cdots \, \Psi ^ { \pm } ( t, H ) & = \frac { \pm 1 } { \Gamma ( H + \frac { 1 } { 2 } ) } \frac { \mathrm d } { \mathrm d s } \left ( \int _ { [ \pm T, s ] } B _ { u } ( s - u ) _ { \pm } ^ { H - \frac { 1 } { 2 } } \mathrm d u \right ) \{ ( t ) - ( 0 ) \} \\ & \pm \frac { H - \frac { 1 } { 2 } } { \Gamma ( H + \frac { 1 } { 2 } ) } \int _ { ( \pm \infty, \pm T ] } B _ { u } \left [ ( t - u ) _ { \pm } ^ { H - 3 / 2 } - ( - u ) _ { \pm } ^ { H - 3 / 2 } \right ] \mathrm d u. \quad \ ( 2. 4 ) \\ \cdots \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm. \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm + \, u } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \ \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ \ \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \cdots \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ { \pm } \, \Psi ^ \ \Ps }$$

Proof. Without any loss of generality, we may only prove the equality for the term B + ( t, H ). Let first consider the case t ∈ R and H ∈ ( 1 2 , 1). Using Ito's lemma on the product B u ( t -u ) H -1 / 2 + , we obtain for any x ≤ t

$$B _ { t } ( t - t ) ^ { H - 1 / 2 } = B _ { x } ( t - x ) ^ { H - 1 / 2 } + \int _ { x } ^ { t } ( t - u ) ^ { H - 1 / 2 } \mathrm d B _ { u } - ( H - \frac { 1 } { 2 } ) \int _ { x } ^ { t } B _ { u } ( t - u ) ^ { H - 3 / 2 } \mathrm d u.$$

Hence, subtracting the case t = 0,

$$\lambda _ { x } \left [ ( t - x ) _ { + } ^ { H - 1 / 2 } - ( - x ) _ { + } ^ { H - 1 / 2 } \right ] & = - \int _ { x } ^ { t } \left [ ( t - u ) _ { + } ^ { H - 1 / 2 } - ( - u ) _ { + } ^ { H - 1 / 2 } \right ] \mathrm d B _ { u } \\ & \quad + ( H - \frac { 1 } { 2 } ) \int _ { x } ^ { t } B _ { u } \left [ ( t - u ) _ { + } ^ { H - 3 / 2 } - ( - u ) _ { + } ^ { H - 3 / 2 } \right ] \mathrm d u.$$

$$- \, \frac { 1 } { 2 } ) \int _ { x } \, B _ { u } \left [ \left ( t - u \right ) _ { + } ^ { H - 3 / 2 } - \left ( - u \right ) _ { + } ^ { H - 3 / 2 } \right ] \mathrm d u.$$

The Brownian motion is known to satisfy lim u →±∞ | B u | / u | | 1 / 2+ ε = 0. Moreover, we observe that [ ( t -x ) H -1 / 2 + - -( x ) H -1 / 2 + ] ∼ -∞ ( H -1 ) t H -1 / 2 ( -x ) H -3 / 2 . Therefore, as 3 -H &gt; 1 ,

$$x \rightarrow - \infty$$

2 2 2 B x [ ( t -x ) H -1 / 2 + - -( x ) H -1 / 2 + ] a.s. - - - - -→ 0 .

Similarly, using the dominated convergence theorem,

$$\int _ { x } ^ { t } B _ { u } \left [ ( t - u ) _ { + } ^ { H - 3 / 2 } - ( - u ) _ { + } ^ { H - 3 / 2 } \right ] \mathrm d u \xrightarrow { a. s. } \int _ { - \infty } ^ { t } B _ { u } \left [ ( t - u ) _ { + } ^ { H - 3 / 2 } - ( - u ) _ { + } ^ { H - 3 / 2 } \right ] \mathrm d u,$$

Finally, owing to the L 2 -continuity of the stochastic integral,

$$\int _ { x } ^ { t } \left [ ( t - u ) _ { + } ^ { H - 1 / 2 } - ( - u ) _ { + } ^ { H - 1 / 2 } \right ] \mathrm d B _ { u } \xrightarrow { L ^ { 2 } ( \Omega ) } \int _ { - \infty } ^ { t } \left [ ( t - u ) _ { + } ^ { H - 1 / 2 } - ( - u ) _ { + } ^ { H - 1 / 2 } \right ] \mathrm d B _ { u },$$

which proves the expected equality.

Let now consider the second case H ∈ (0 , 1 2 ). The calculus is quite similar, with only a slight modification in the application of Ito's lemma between T and t -ε .

$$& ( B _ { t - \varepsilon } - B _ { t } ) \varepsilon ^ { H - 1 / 2 } + B _ { t } ( t - T ) ^ { H - 1 / 2 } \\ & = B _ { T } ( t - T ) ^ { H - 1 / 2 } + \int _ { T } ^ { t - \varepsilon } ( t - u ) ^ { H - 1 / 2 } \mathrm d B _ { u } - ( H - \frac { 1 } { 2 } ) \int _ { T } ^ { t - \varepsilon } ( B _ { u } - B _ { t } ) ( t - u ) ^ { H - 3 / 2 } \mathrm d u, \\ \hat { \nu } & \hat { \nu } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim } \hat { \sim }$$

for any ε &gt; 0. Owing to the Hölder continuity of the sample paths of B , the previous expression converges when ε → 0. In addition, the calculus presented in the case H &gt; 1 2 holds as well on the interval [ x, T ], therefore proving Equation (2.3). Finally, the second form (2.4) of the fractional Brownian field is obtained using an alternative representation of the RiemannLiouville fractional derivative of a function on the interval [ T, t ] (see the book of Samko et al. [33] for more information on the subject).

This type of representation has already been exhibited and studied by Picard [32] for the fractional Brownian motion, even though his proof is based on a different approach (approximation of the Brownian motion with Lipschitz functions). The case H &gt; 1 2 has also been considered by Takashima [36].

↦

As presented in the next simple lemma and owing to the Hölder continuity of the sample paths of Brownian motion, we obtain a deterministic representation of the fractional Brownian field ( t, H ) → B t, H ( ) which is continuous in both variables. It will be very useful in the two next sections to study in the sample path properties of the multifractional Brownian motion.

Lemma 2.1. Let ( B t ) t ∈ R be a continuous Brownian motion and B t, H ( ) be the fractional Brownian field defined in Proposition 2.1.

↦

Then, with probability one for all ( t 0 , H 0 ) ∈ R × (0 , 1) , the field ( t, H ) → B t, H ( )( ω ) is continuous at ( t 0 , H 0 ) .

Proof. As previously, we only need to consider the component B + ( t, H ). The continuity on the domains R × (0 , 1 2 ) and R × ( 1 2 , 1) is a simple consequence of the dominated convergence theorem: the deterministic integral is split in two parts on which the theorem can be applied with different bounds.

To threat the case H 0 = 1 2 , we observe that the representation (2.3) in fact holds for all H ∈ (0 1). Then, similarly to the previous case, the dominated convergence theorem implies , that the first term converges to zero when H → 1 2 whereas the second one simply converges to B t 0 .

## 2.1. Well-balanced case

As observed by Stoev and Taqqu [35], the introduction of the mBm is more delicate in the well-balanced case a + = a -. Indeed, when H = 1 2 , the aforementioned definition (2.1) of the fractional Brownian field B t, H ( ) implies that B , ( · 1 2 ) = 0 almost surely. Hence, as presented in [35], to obtain a well-defined mBm, one needs to change the normalisation term, and more specifically, add a diverging factor: 1 / H ( -1 2 ).

Considering this slightly modified definition of the well-balanced mBm, we observe that it corresponds to a particular deterministic expression.

Proposition 2.2. Suppose ( B t ) t ∈ R + is a continuous Brownian motion and a + = a -= 1 . Then, for all ( t, H ) ∈ R × (0 , 1) ,

$$B ( t, H ) = \frac { 1 } { \Gamma ( H + \frac { 1 } { 2 } ) } \int _ { \mathbf R } B _ { u } \left [ ( t - u ) ^ { < H - 3 / 2 > } - ( - u ) ^ { < H - 3 / 2 > } \right ] \mathrm d u \ \ a. s. \quad \ \ ( 2. 5 )$$

where for any x ∈ R , x &lt;α&gt; := sign( x ) | x | α .

Proof. Let ( t, H ) ∈ R × (0 , 1), δ &gt; 0 and T ± = t ∓ δ . Owing to the representation obtained in Proposition 2.1 and using the well-balanced re-normalisation,

$$B ( t, H ) = \frac { 1 } { \Gamma ( H + \frac {$$

Without any restriction, we may assume that 0 / B t, δ ∈ ( ). Then, let us observe the component of the integral corresponding to the interval [ t -δ, t + ]. It is equal to δ

$$\int _ { B ( t, \delta ) } ( B _ { u } - B _$$

For any ε &gt; 0, we note that ∫ B t,δ ( ) \ B t,ε ( ) B t t ( -u ) &lt;H -3 / &gt; 2 d u = 0. Furthermore, the other part converges as well to

$$\begin{array} {$$

since B is Hölder continuous at t . Hence, the integral converges when ε → 0 and we obtain the representation presented in Equation (2.5).

On the contrary to the formulas presented in Proposition 2.1, we observe that the representation (2.5) is an improper deterministic integral when H ≤ 1 2 .

Also note that we may check that when H = 1 2 , the corresponding Gaussian process is, up to a constant, a Brownian motion. More precisely, Stoev and Taqqu [35] have shown that it has the following stochastic integral representation

$$B ( t, \frac { 1 } { 2 } ) = \frac { 1 } { \Gam$$

As proved in the next proposition, the well-balanced case is in fact a classic mBm defined with respect to the latter Brownian motion.

Proposition 2.3. Suppose B is a continuous Brownian motion and B t, H ( ) denotes the wellbalanced fractional Brownian field. Define the Brownian motion ˜ B by

$$\forall t \in \mathbf R ; \ \ \widetilde { B } _ {$$

Then, for all t ∈ R and H ∈ (0 , 1) ,

$$B ( t, H ) = \frac { 1 } { ( H - \frac { 1 } { 2 } ) \tan \left ( ( H + \frac { 1 } { 2 } ) \frac { \pi } { 2 } \right ) } \left ( \widetilde { B } ^ { + } ( t, H ) - \widetilde { B } ^ { - } ( t, H ) \right ) \ \ a. s.$$

where ˜ B ± ( t, H ) designates the fractional Brownian field with respect to ˜ B . Note that the normalising function is continuous at H = 1 2 .

Proof. There exist several ways to prove this result. One might directly obtain the expression of the fractional Brownian field ˜ B ± ( t, H ) using Propositions 2.1 and 2.2, but the calculus is unfortunately a bit tricky.

The proof happens to be simpler if we consider the Fourier representation of B t, H ( ) obtained by Stoev and Taqqu [35]. For all t ∈ R and H ∈ (0 , 1),

$$B ^ { \pm } ( t, H ) = \frac { 1 } { 2 \sqrt { 2 \pi } } \int _ { \mathbf R } \frac { e ^ { i t \xi } - 1 } { | \xi | ^ { H + 1 / 2 } } e ^ { \mp i \text{sign} ( \xi ) ( H + 1 / 2 ) \pi / 2 } \text{d} \widehat { W } ( \xi ),$$

where ̂ W is a complex Gaussian measure such that ̂ W is a complex Wiener measure on R + and ̂ W A ( ) = ̂ W ( -A ) for any A ∈ B ( R + ). The well-balanced fractional Brownian field then has the following representation

$$B ( t, H ) = \frac { 2 \cos \left ( ( H + \frac { 1 } { 2 } ) \frac { \pi } { 2 } \right ) } { \sqrt { 2 \pi } ( H - \frac { 1 } { 2 } ) } \int _ { \mathbf R } \frac { e ^ { i t \xi } - 1 } { | \xi | ^ { H + 1 / 2 } } \text{d} W ( \xi ).$$

We introduce a slightly modified complex Gaussian measure ˜ W , given by

$$\forall A \in \mathcal { B } ( \mathbf R _ { + } ) ; \ \ W ( A ) \coloneqq i W ( A ) \ \text{ and } \ W ( - A ) \coloneqq W ( A ) = - i W ( - A ).$$

Then, let ˜ B ± ( t, H ) denote the fractional Brownian field defined with respect to ˜ W . Since ˜ W (d ξ ) = i sign( ξ W ) ̂ (d ξ ),

$$\widetilde { B } ^ { + } ( t, H ) - \widetilde { B } ^ { - } ( t, H ) & = \frac { 2 \sin \left ( ( H + \frac { 1 } { 2 } ) \frac { \pi } { 2 } \right ) } { \sqrt { 2 \pi } } \int _ { \mathbf R } \frac { e ^ { i t \xi } - 1 } { | \xi | ^ { H + 1 / 2 } } ( - i ) \, \text{sign} ( \xi ) \, \text{d} \, \widetilde { W } ( \xi ) \\ & = ( H - \frac { 1 } { 2 } ) \tan \left ( ( H + \frac { 1 } { 2 } ) \frac { \pi } { 2 } \right ) \, B ( t, H ).$$

To conclude the proof, we need to check that the field ˜ B ± ( t, H ) corresponds to the same process defined in the proposition, i.e. the integral with respect to the Brownian motion ˜ B . For this purpose, we simply need to consider the previous equality when H = 1 2 . We first note that lim H → 1 / 2 ( H -1 2 ) tan ( ( H + 1 2 ) π 2 ) = 1. Hence,

$$\widetilde { B } _ { t } \coloneqq \frac { 1 } { 2 } B ( t, \frac { 1 } { 2 } ) = \frac { 1 } { 2 } \left ( \widetilde { B } _ { + } ( t, \frac { 1 } { 2 } ) - \widetilde { B } _ { - } ( t, \frac { 1 } { 2 } ) \right ).$$

According to Equation (2.1), this last equality proves that ˜ B is the Brownian motion with respect to which are defined the fractional Brownian fields, and thus we clearly obtain the expected representation of B t, H ( ).

The equality in law between the two processes can be obtained more directly using Theorem 4.1 from [35]. Nevertheless, Proposition 2.3 indicates precisely how these two fractional Brownian fields are related, and proves that the well-balanced one is simply a classic field with a change of Brownian motion. Besides, this statement allows us to simply the study of the sample path properties of the mBm since we do not need to consider separately the specific well-balanced case a + = a -.

## 3. Hölder and 2-microlocal regularity

Following the convention initiated by [35], we will call multifractional Brownian motion a process X which has the following form:

$$\forall t \in \mathbf R ; \ \ X _ { t } = a ^ { + } B ^ { + } \left ( t, H ( t ) \right ) + a ^ { - } B ^ { - } \left ( t, H ( t ) \right ),$$

/negationslash

↦

where ( a -, a + ) ∈ R 2 \ { (0 0) , } . As previously said, we may also suppose without any loss of generality that a + = a -. The Hurst function t → H t ( ) ∈ (0 1) will be assume to be at least a , deterministic càdlàg function.

/negationslash

We may nevertheless note that several of the following results remain valid if H ( ) is replaced · by a stochastic process with càdlàg sample paths. This last model has already been considered by Papanicolaou and Sølna [28] and Ayache and Taqqu [2], where the latter has obtained some general sample path properties. For any càdlàg function H ( ), · the set { t ∈ R : ∆ H t ( ) = 0 } is countable, and therefore we can easily see from the definition of the mBm that with probability one,

/negationslash

$$\forall t \in \mathbf R ; \ \Delta H ( t ) \neq 0 \ \iff \ \Delta X _ { t } \neq 0.$$

/negationslash

Furthermore, ∆ X t is a centered Gaussian variable whose variance is proportional to ∆ H t ( ) 2 .

In this section, we will use specific functional spaces, called 2-microlocal spaces, to characterize the fine sample path properties of the mbm. The use of the unusual framework is particularly motivated by the deterministic fractional integrals and derivatives of Brownian motion exhibited in Proposition 2.1. Indeed, the latter spaces behave particularly well under the action of fractional operators, which will allow us to deduce many properties on the mBm from the geometry of the Brownian motion.

Historically, the 2-microlocal formalism has been introduced by Bony [8] to study the singularities of the solutions of PDEs, and has only been recently introduced in a stochastic frame by Herbin and Lévy Véhel [16] to characterise the fine regularity of processes. Several characterizations of the 2-microlocal spaces have been investigated in the literature, including the following time-domain one due to Kolwankar and Lévy Véhel [22] and Seuret and Lévy Véhel [34].

Definition 3.1. Let t ∈ R , s ′ ≤ 0 and σ ∈ (0 , 1) such that σ -s ′ / ∈ N . A function f : R → R belongs to the 2-microlocal space C σ,s ′ t if there exist C &gt; 0 , ρ &gt; 0 and a polynomial P t such that

$$| (.$$

$$\iota ^ { \prime } \quad \left | \left ( f ( u ) - P _ { t } ( u ) \right ) - \left ( ( f ( v ) - P _ { t } ( v ) ) \right | \leq C | u - v | ^ { \sigma } \left ( | u - t | + | v - t | \right ) ^ { - s ^ { \prime } },$$

$$( 3. 1 )$$

for all u, v ∈ B t, ρ ( ) .

In this work, it will only be necessary to consider the case σ ∈ (0 , 1). Nevertheless, Definition 3.1 can be extended σ ∈ R Z \ . We refer to [34, 5] for a more complete presentation on the subject.

As previously outlined, a very useful properties of the 2-microlocal spaces is their stability under the application of fractional integrations and derivations:

$$\forall \alpha > 0 ; \ \ f \in C _ { t } ^ { \sigma, s ^ { \prime } } \ \iff \ I _ { \pm } ^ { \alpha } f \in C _ { t } ^ { \sigma + \alpha, s ^ { \prime } },$$

where I α ± f denotes the fractional integral of f of order α : ( I α ± f ) ( ) := t ∫ R ( t -u ) α -1 ± f ( u ) d u . Following the presentation of these functional spaces, we can introduce a related regularity tool, named the 2-microlocal frontier and defined by

$$\forall s ^ { \prime } \in \mathbf R ; \quad \sigma _ { f, t } ( s ^ { \prime } ) = \sup \{ \sigma \in \mathbf R \, \colon f \in C _ { t } ^ { \sigma, s ^ { \prime } } \}.$$

↦

We observe that when using the 2-microlocal formalism, the local regularity at t is characterised by a function, not a single coefficient. Owing to inclusion properties of the 2-microlocal spaces, the map s ′ → σ f,t ( s ′ ) is well-defined and satisfies several interesting properties:

- · σ f,t ( ) is a concave non-decreasing function; ·

## · σ f,t ( ) has left and right derivatives between 0 and 1. ·

Furthermore, if the modulus of continuity of f is such that ω f ( h ) = O (1 / | log( h ) ), the pointwise | and local Hölder exponents can be retrieved from the frontier:

$$\alpha _ { f, t } = - \inf \{ s ^ { \prime } \, \colon \sigma _ { f, t } ( s ^ { \prime } ) \geq 0 \} \ \text{ and } \ \widetilde { \alpha } _ { f, t } = \sigma _ { f, t } ( 0 ).$$

Finally, owing to the stability of the 2-microlocal spaces, the application of a fractional integration simply corresponds to a translation of the 2-microlocal frontier along the vertical axis:

$$\forall s ^ { \prime } \in \mathbf R ; \ \ \sigma _ { I _ { \pm } ^ { \alpha } f, t } ( s ^ { \prime } ) = \sigma _ { f, t } ( s ^ { \prime } ) + \alpha.$$

The 2-microlocal frontier of some Gaussian processes has been determined by Herbin and Lévy Véhel [16]. For instance, a fractional Brownian motion B H , H ∈ (0 1) , satisfies with probability one and for all t ∈ R :

$$\forall s ^ { \prime } \in \mathbf R ; \ \ \sigma _ { B ^ { H }, t } ( s ^ { \prime } ) = \left ( H + s ^ { \prime } \right ) \wedge H.$$

To investigate the 2-microlocal and Hölder regularity of the multifractional Brownian motion, we first need to properly study the behaviour of the fractional Brownian field and its successive partial derivatives. Hence, in the next proposition, we consider uniformly on the variables t and H the 2-microlocal frontier of ∂ k H B t, H ( ) along the time axis.

↦

Proposition 3.1. Suppose B t, H ( ) is a fractional Brownian field. For any H ∈ (0 1) , and any k ∈ N , let B H,k denotes the process t → ∂ k H B t, H ( ) . Then, with probability one, for all H ∈ (0 , 1) , k ∈ N and t ∈ R

$$\forall s ^ { \prime } \in \mathbf R ; \quad \sigma _ { B ^ { H, k }, t } ( s ^ { \prime } ) = \left ( H + s ^ { \prime } \right ) \wedge H.$$

Proof. Since we use the deterministic representation obtained in Proposition 2.1, we can set a fix ω ∈ Ω in the proof. For the sake of readability, let us first consider H &gt; 1 2 and remove the normalisation constant 1 / Γ ( H -1 2 ) . Then, for any k ∈ N , the partial derivative with respect to H is equal to

$$\partial _ { H } ^ { k } \, B ^ { \pm } ( t, H ) = \int _ { \mathbf R } B _ { u } \left [ ( t - u ) _ { \pm } ^ { H - 3 / 2 } \log ^ { k } ( t - u ) _ { \pm } - ( - u ) _ { \pm } ^ { H - 3 / 2 } \log ^ { k } ( - u ) _ { \pm } \right ] \mathrm d u.$$

Suppose ( T + , T -) is a non-empty open interval. Since the 2-microlocal frontier is a purely local property of the sample paths, we may restrict our study to t ∈ ( T + , T -) with any loss of generality. Then, we note that the component

↦

$$t \mapsto \int _ { ( \mp \infty, T _ { \pm } ] } B _ { u } \left [ ( t - u ) _ { \pm } ^ { H - 3 / 2 } \log ^ { k } ( t - u ) _ { \pm } - ( - u ) _ { \pm } ^ { H - 3 / 2 } \log ^ { k } ( - u ) _ { \pm } \right ] \mathrm d u$$

is a smooth function, and therefore does not have any influence on the 2-microlocal regularity. Hence, we focus on the study of the remaining integral

↦

$$t \in ( T _ { + }, T _ { - } ) \mapsto \int _ { [ T _ { + }, T _ { - } ] } B _ { u } ( t - u ) _ { \pm } ^ { H - 3 / 2 } \log ^ { k } ( t - u ) _ { \pm } \mathrm d u,$$

when H &gt; 1 2 . Similarly, the case H &lt; 1 2 is reduced to the analysis of

↦

$$t \in ( T _ { + }, T _ { - } ) \mapsto \frac { \mathrm d } { \mathrm d s } \left ( \int _ { [ T _ { \pm }, s ] } B _ { u } ( s - u ) _ { \pm } ^ { H - \frac { 1 } { 2 } } \log ^ { k } ( t - u ) _ { \pm } \mathrm d u \right ) ( t ).$$

If we define the process Y t = B t 1 t ∈ [ T + ,T -] , we observe that the two previous formulas correspond respectively to fractional-like integral and derivative of Y on the real axis. Compared to the usual definition, there is only a slight modification in the kernel with a logarithmic term.

Furthermore, when H &lt; 1 2 , the fractional-like derivative is defined as the classic derivative of a integral of order H + 1 2 . Since we know the behaviour of the 2-microlocal regularity under derivation, we may restrict our study to the fine analysis of the following process

↦

$$t \in ( T _ { + }, T _ { - } ) \longmapsto ( I _ { \pm } ^ { \alpha, k } Y ) ( t ) \coloneqq \int _ { \mathbf R } Y _ { u } \, ( t - u ) _ { \pm } ^ { \alpha - 1 } \log ^ { k } ( t - u ) _ { \pm } \mathrm d u,$$

for any α ∈ (0 1) and , k ∈ N .

We expect the 2-microlocal frontier of this integral to behave similarly to the one of a classic fractional integral of Y , since the only difference lies in a logarithmic term. This easiest way to prove this statement is to compute the Fourier transform of the corresponding pseudodifferential operator and show it has an effect equivalent to the fractional integration in the Fourier space.

The computation of Fourier transform of the tempered distribution u α -1 ± log k u ± is similar to the classic calculus on u α -1 ± . We first consider the Laplace transform, with /Rfractur ( z ) &gt; 0. Since this transform is analytic on the previous domain, we can restrict ourselves to the case that z &gt; 0. Then,

$$\psi _ { z } \, \dots \, \dots \, \dots, \\ \forall z > 0 ; \quad \int _ { \mathbf R _ { + } } u ^ { \alpha - 1 } \log ^ { k } ( u ) \, e ^ { - z u } \, d u & = z ^ { - \alpha } \int _ { \mathbf R _ { + } } v ^ { \alpha - 1 } \left ( \log v - \log z \right ) ^ { k } e ^ { - v } \, d v \\ & = z ^ { - \alpha } \sum _ { j = 0 } ^ { k } ( - 1 ) ^ { j } \log ^ { j } ( z ) \, C _ { k } ^ { j } \underbrace { \int _ { \mathbf R _ { + } } v ^ { \alpha - 1 } \log ^ { k - j } ( v ) \, e ^ { - v } } _ { F ( \alpha, k - j ) }.$$

As the last term is also analytic, it implies the equality for any z ∈ C such that /Rfractur ( z ) &gt; 0.

To obtain the Fourier transform of u α -1 ± log k u ± as a tempered distribution, we simply have to consider z = iξ + ε and let ε → 0. The scalar product with a test function converges, and proves that the generalized Fourier transform is equal to

$$& \mathcal { F } _ { g } ( u _ { + } ^ { \alpha - 1 } \log ^ { k } u _ { + } ) ( \xi ) \\ & = e ^ { - i \, \text{sign} ( \xi ) \frac { \pi } { 2 } \alpha } | \xi | ^ { - \alpha } \frac { 1 } { \sqrt { 2 \pi } } \sum _ { j = 0 } ^ { k } ( - 1 ) ^ { j } ( \log | \xi | + i \, \text{sign} ( \xi ) \frac { \pi } { 2 } ) ^ { j } \cdot C _ { k } ^ { j } F ( \alpha, k - j ).$$

Therefore, the Fourier transform corresponds to a weighted sum of elements | ξ | -α log n | ξ | , with a phase correction term.

To express the action of these pseudo-differential operators on the regularity of the Brownian motion, we rely on the original definition of the 2-microlocal spaces in the Fourier space and

make use of the arguments presented by Jaffard [18] to prove the stability. More precisely, the aforementioned definition is based on the Littlewood-Paley decomposition of tempered distributions. Shortly, consider φ ∈ S ( R ) such that its Fourier transform satisfies

$$\widehat { \phi } ( \xi ) = \begin{cases} 1 & \text{si} & | \xi | \leq 1 / 2 ; \\ 0 & \text{si} & | \xi | \geq 1, \end{cases}$$

and define for any tempered distribution f ∈ S ′ ( R ) and for all j ∈ N , φ j ( x ) = 2 j φ (2 j x ), ψ j = φ j +1 -φ j , S f 0 = φ ∗ f and ∆ j f = ψ j ∗ f . Then, the Littlewood-Paley decomposition tells us that f can be expressed as following

$$f = S _ { 0 } f + \sum _ { j = 0 } ^ { \infty } \Delta _ { j } f,$$

where we observe that every term in the sum belongs to L 2 ( R ). Roughly speaking, this decomposition corresponds to the application of band-pass filters in the Fourier space. As presented by Bony [8], Meyer [25], f ∈ C σ,s ′ t if and only there exists C &gt; 0 such that for all j &gt; 0,

$$| \Delta _ { j } f ( u ) | \leq C \, 2 ^ { - j \sigma } \left ( 2 ^ { - j } + | u - t | \right ) ^ { - s ^ { \prime } }.$$

Observe that, up to a simple translation, we may assume that t = 0. Then, define the rescaled function U j ( x ) = (∆ j f )(2 -j u ). Its Fourier transform ̂ U j is carried by { ξ : 1 2 &lt; ξ | | &lt; 2 } and U j is such that | U j ( u ) | ≤ C 2 -j ( σ -s ′ ) ( 1+ | u | ) -s ′ . Applying the operator | ξ | -α log k | ξ | in the Fourier space on f , we obtain a corresponding collection of rescaled functions ( V j ) j&gt; 0 which satisfy for every j &gt; 0,

$$\widehat { V } _ { j } ( \xi ) = \widehat { U } _ { j } ( \xi ) \cdot \left [ \sum _ { \ell = 0 } ^ { k } c ( \ell, k ) \, | \xi | ^ { - \alpha } \log ^ { \ell } | \xi | \cdot j ^ { k - \ell } \, 2 ^ { - j \alpha } \right ],$$

where c /lscript, ( k ) is a constant only depending on /lscript and k . Note that similarly to the classic fractional integral, the Fourier multipliers ̂ K α,/lscript ( ξ ) := | ξ | -α log /lscript | ξ | coincide with Schwartz functions on 1 2 &lt; ξ | | &lt; 2. Coming back into the time space, V j is given by

$$V _ { j } = \sum _ { \ell = 0 } ^ { k } j ^ { k - \ell } 2 ^ { - j \alpha } \, c ( \ell, k ) K _ { \alpha, \ell } \ast U _ { j },$$

where K α,/lscript ∈ S ( R ) for every /lscript . The convolution with a Schwartz function K α,/lscript preserves the polynomial decay, hence for every j &gt; 0, | V j ( u ) | ≤ C j k 2 -j ( σ + α -s ′ ) ( 1 + | u | ) -s ′ . Rescaling by a factor 2 , the latter upper bound leads to the characterization of 2-microlocal spaces, proving j that

$$\forall s ^ { \prime }, \sigma \in \mathbf R ; \ \ Y \in C _ { t } ^ { \sigma, s ^ { \prime } } \Longrightarrow I _ { \pm } ^ { \alpha, k } Y \in C _ { t } ^ { \sigma - \varepsilon, s ^ { \prime } }$$

for any ε &gt; 0. Note that the property is slightly weaker than the classic stability under the action of fractional integration due to the logarithmic terms.

The converse case is obtained similarly,

$$\forall s ^ { \prime }, \sigma \in \mathbf R ; \ \ I _ { \pm } ^ { \alpha, k } Y \in C _ { t } ^ { \sigma, s ^ { \prime } } \Longrightarrow Y \in C _ { t } ^ { \sigma - \varepsilon, s ^ { \prime } }.$$

The same property holds for any weighted sum of I α,k -Y and I α,k + Y , and therefore, the 2microlocal frontier of I α,k ± Y satisfies for any H ∈ (0 , 1) and k ∈ N

$$\forall s ^ { \prime } \in \mathbf R ; \ \ \sigma _ { I _ { \pm } ^ { \alpha, k } Y, t } ( s ^ { \prime } ) = \sigma _ { Y, t } ( s ^ { \prime } ) + H - \frac { 1 } { 2 } = ( s ^ { \prime } + H ) \wedge H.$$

It concludes the proof, since we have shown previously that the 2-microlocal frontier of B H,k is equal to the regularity I α,k ± Y .

↦

Ayache [1] has obtained a similar result on the field ( t, H ) → B + ( t, H ) based on a multiresolution analysis of the latter.

As outlined in the introduction, the Hölder regularity of a general multifractional Brownian motion has already been investigated by Herbin [15] and Ayache [1]. In particular, the latter has proved that with probability one, the local Hölder exponent is given by

$$\forall t \in \mathbf R \, \sum \{ 0 \} ; \ \widetilde { \alpha } _ { X, t } = H ( t ) \wedge \widetilde { \alpha } _ { H, t }.$$

Note that the behaviour at t = 0 is slightly different as the form of the variance of increments differs at this point.

↦

Let us first obtain a estimate for the general 2-microlocal frontier. For this purpose, we need to introduce the multiplicity of the component h → ( B t, h ( ) -B t, H ( ) ) of the fractional field at ( t, H ) ∈ R × (0 1), i.e. ,

/negationslash

$$m _ { t, H } = \inf \{ k \in \mathbf N \, \sum \{ 0 \} \colon \partial _ { H } ^ { k } B ( t, H ) \neq 0 \}.$$

Theorem 3.1. Suppose X is a multifractional Brownian motion with Hurst function H . Then, with probability one, the 2-microlocal frontier of X at any t ∈ R satisfies

$$\sigma _ { X, t } ( s ^ { \prime } ) \geq \left ( s ^ { \prime } + H ( t ) \right ) \wedge H ( t ) \wedge \sigma _ { H, t } ( s ^ { \prime } + H ( t ) ) \wedge m _ { t, H ( t ) } \sigma _ { H, t } ( s ^ { \prime } / m _ { t, H ( t ) } ).$$

for all s ′ ≥ -α X,t . Furthermore, when m t,H t ( ) = 1 , the frontier is equal to

$$\forall s ^ { \prime } \geq - \alpha _ { X, t } ; \ \ \sigma _ { X, t } ( s ^ { \prime } ) = \left ( s ^ { \prime } + H ( t ) \right ) \wedge H ( t ) \wedge \sigma _ { H, t } ( s ^ { \prime } ).$$

Proof. Without any loss of generality, we can assume that ˜ α H,t &gt; 0, since otherwise the lower bound is equal to 0. Then, let ε &gt; 0, ρ &gt; 0, t ∈ R and u, v be in the neighbourhood B t, ρ ( ). We first estimate the size of the increments X u -X v :

$$B ( u, H _ { u } ) - B ( v, H _ { v } ) = \left ( B ( v, H _ { u } ) - B ( u, H _ { u } ) \right ) + \left ( B ( v, H _ { v } ) - B ( v, H _ { u } ) \right ).$$

Using a simple integration by parts, we observe that the first term is equal to

$$B ( v, H _ { u } ) - B ( u, H _ { u } ) = B ( v, H _ { t } ) - B ( u, H _ { t } ) + \int _ { H _ { t } } ^ { H _ { u } } \left ( \partial _ { H } B ( v, h ) - \partial _ { H } B ( u, h ) \right ) \text{d} h.$$

Owing to Lemma 3.1, for any ε &gt; 0, there exists c 1 &gt; 0 such that

$$\forall u, v \in B ( t, \rho ), \ \forall h \in [ H _ { t }, H _ { u } ] ; \ \left | \partial _ { H } B ( v, h ) - \partial _ { H } B ( u, h ) \right | \leq c _ { 1 } | u - v | ^ { h - \varepsilon }.$$

Therefore, the integral term satisfies

$$\left | \int _ { H _ { t } } ^ { H _ { u } } ( \partial _ { H } B ( v, h ) - \partial _ { H } B ( u, h ) ) d h \right | \leq c _ { 1 } | H _ { t } - H _ { u } | \cdot | v - u | ^ { H _ { t } \wedge H _ { u } - \varepsilon },$$

and is thus negligible in front of the first one. Still using Lemma 3.1, we know that for all u, v ∈ B t, ρ ( )

$$^ { \cdot \cdot \cdot } _ { \cdot } \,, & \quad | B ( v, H _ { t } ) - B ( u, H _ { t } ) | \leq c _ { 2 } | u - v | ^ { \sigma } \left ( | u - t | + | v - t | \right ) ^ { - s ^ { \prime } },$$

where s ′ ≥ -H t and σ &lt; ( H t + s ′ ) ∧ H t . This inequality corresponds to the first term in the lower bound.

Using Taylor's formula, we can obtain an estimate of the second component ( B v, H ( v ) -B v, H ( u )).

$$B ( v, H _ { v } ) - B ( v, H _ { u } ) = \sum _ { k = 1 } ^ { m - 1 } \frac { \partial _ { H } ^ { k } B ( v, H _ { u } ) } { k! } ( H _ { v } - H _ { u } ) ^ { k } + \int _ { H _ { u } } ^ { H _ { v } } \frac { \partial _ { H } ^ { m } B ( v, h ) } { ( m - 1 )! } ( H _ { v } - h ) ^ { m - 1 } d h,$$

where m denotes the coefficient m t,H t previously introduced. There exists a neighbourhood B t, ρ ( ) of t such that for every k ∈ { 1 , . . . , m -1 } and for all v ∈ B t, ρ ( ), | ∂ k H B v, H ( u ) | ≤ c 3 | v - | t H t -ε . Hence the term k = 1 is dominant in the sum and for all u, v ∈ B t, ρ ( )

$$| \partial _ { H } B ( v, H _ { u } ) | \cdot | H _ { v } - H _ { u } | \leq c _ { 4 } | u - v | ^ { \sigma } \left ( | u - t | + | v - t | \right ) ^ { - s + H _ { t } - \varepsilon },$$

for any s ′ ≥ -α H,t and σ &lt; σ H,t ( s ′ ). This inequality leads to the second term in the lower bound.

Finally, let consider the last integral. ∂ m H B v, h ( ) can be simply bound by a constant, implying that the integral is upper bounded by | H u -H v | m . The third lower bound is a direct consequence of this estimate.

To end the proof, we study the case m t,H t ( ) = 1. The lower bound is clearly equal to ( H t + s ′ ) ∧ H t ∧ σ H,t ( s ′ ). Let set s ′ ≥ -α X,t . Suppose first that ( H t + s ′ ) ∧ H &lt; σ t H,t ( s ′ ). Then, the term B v, H ( u ) -B u,H ( u ) will dominate in the increment X u -X v , therefore implying the inequality σ X,t ( s ′ ) ≤ ( H t + s ′ ) ∧ H t . In the other case, we observe that

$$B ( v, H _ { v } ) - B ( v, H _ { u } ) = \partial _ { H } B ( v, H _ { v } ) ( H _ { u } - H _ { v } ) + \int _ { H _ { v } } ^ { H _ { u } } \partial _ { H } ^ { 2 } B ( v, h ) ( H _ { v } - h ) \text{d} h,$$

There exists ρ &gt; 0 and c 5 &gt; 0 such that for all v ∈ B t, ρ ( ), | ∂ H B v, H ( v ) | &gt; c 5 , implying that the regularity of this component is σ H,t ( s ′ ) Furthermore, as previously noted, the second term is upper-bounded by ( H u -H v ) 2 , and thus is negligible in front of the first one. Therefore, the increment X u -X v behave similarly to H u -H v , proving that σ X,t ( s ′ ) ≤ σ H,t ( s ′ ).

As observed in the previous proof, the study of the general form of the 2-microlocal frontier is quite technical, in particular because of the interaction which may appear between the oscillations of the Hurst function H and the derivatives of the fractional Brownian field B t, H ( ). Nevertheless, a complete estimate can be obtained for the pointwise Hölder regularity.

Proposition 3.2. Suppose X is a multifractional Brownian motion with Hurst function H . Then, with probability one, the pointwise exponent is equal to

$$\forall t \in \mathbf R ; \quad \alpha _ { X, t } = H ( t ) \wedge m _ { t, H ( t ) } \alpha _ { H, t }.$$

Proof. Let us proceed similarly to the proof of Theorem 3.1 and study the increment X u -X t = ( B u,H ( u ) -B u,H ( t ) ) + ( B u,H ( t ) -B t, H ( t ) ) . Owing to Lemma 3.1, the pointwise regularity of the second term is exactly H t ( ).

Then, we consider a slightly different Taylor expansion on the first increment,

$$B ( u, H _ { u } ) - B ( u, H _ { t } ) = \sum _ { k = 1 } ^ { m } \frac { \partial _ { H } ^ { k } B ( u, H _ { t } ) } { k! } ( H _ { u } - H _ { t } ) ^ { k } + \int _ { H _ { t } } ^ { H _ { u } } \frac { \partial _ { H } ^ { m + 1 } B ( u, h ) } { m! } ( H _ { u } - h ) ^ { m } d h.$$

For every k ∈ { 1 , . . . , m - } 1 , there exists a simple upper bound : | ∂ k H B u,H ( t ) || H u -H t | k ≤ c 1 | u - | t H t + kα H,t -ε . Hence, the pointwise regularity of this type of term is strictly greater than H t ( ), they are therefore negligible. In the case k = m , we note that there exists ρ &gt; 0 and c 2 &gt; 0 such that for all u ∈ B t, ρ ( ), | ∂ k H B u,H ( t ) | ≥ c 2 . Hence, the pointwise exponent of this component is exactly m α · X,t . Finally, similarly to the previous proof, the last term can be simply upper bounded by c 3 | H u -H v | m +1 ≤ c 3 | u - | t ( m +1) α H,t . This last estimate is sufficient to conclude the proof.

The previous result constitutes an extension of the estimate obtained by Ayache [1]. As pointed out by the latter, it shows that the multifractional Brownian motion can have a nondeterministic pointwise regularity. More precisely, Ayache [1] has proved that for any ε &gt; 0, with positive probability,

$$\dim _ { \mathrm H } \{ t \in [ a, b ] \, \colon m _ { t, H ( t ) } > 1 \} \geq 1 - \inf _ { u \in [ a, b ] } H ( u ) - \varepsilon.$$

In the view of the result obtained in Proposition 3.2, it would be interesting to study the Hausdorff dimension of the following sets

$$\forall k \in \mathbf N ; \ \ E _ { k } = \left \{ t \in [ a, b ] \, \colon m _ { t, H ( t ) } = k \right \}.$$

This question seems to be much more difficult than the previous estimate, since it involves the successive derivatives ∂ k H B t, H ( ) of the fractional field and thus the correlation existing between them. Finally, we also note that in terms of regularity, the behaviour at t = 0 is specific, since the multiplicity m 0 ,H (0) is equal to infinity.

To end this section, we present a couple of examples of Hurst functions which illustrate this particular geometry of the multifractional Brownian motion.

Example 1. Suppose H : R → ( 1 2 , 3 4 ) is an Hurst function such that for all t ∈ R , the 2-microlocal local frontier of H at t is equal to

$$\forall s ^ { \prime } \in \mathbf R ; \ \ \sigma _ { H, t } ( s ^ { \prime } ) = \left ( s ^ { \prime } + \frac { 3 } { 8 } \right ) \wedge \frac { 3 } { 8 }.$$

This kind of function can be easily constructed from the sample paths of a fractional Brownian motion with parameter H = 3 8 .

Then, the 2-microlocal frontier of the corresponding multifractional Brownian motion satisfies with probability one:

- · if m t,H t = 1,

$$\forall s ^ { \prime } \in \mathbf R ; \ \ \sigma _ { X, t } ( s ^ { \prime } ) = \sigma _ { H, t } ( s ^ { \prime } ) = \left ( s ^ { \prime } + \frac { 3 } { 8 } \right ) \wedge \frac { 3 } { 8 } ;$$

- · if m t,H t &gt; 1,

$$\forall s ^ { \prime } \in \mathbf R ; \ \ \sigma _ { X, t } ( s ^ { \prime } ) = \left ( s ^ { \prime } + H ( t ) \right ) \wedge \frac { 3 } { 8 }.$$

These two cases are illustrated by Figure 1a.

To obtain the complete determination of 2-microlocal frontier, we use Theorem 3.1, observing that it implies the inequality

$$\sigma _ { X, t } ( s ^ { \prime } ) & \geq \left ( s ^ { \prime } + H ( t ) \right ) \wedge H ( t ) \left ( s ^ { \prime } + \frac { 3 } { 8 } + H ( t ) \right ) \wedge \frac { 3 } { 8 } \wedge \left ( s ^ { \prime } + \frac { 3 m _ { t, H ( t ) } } { 8 } \right ) \wedge \frac { 3 m _ { t, H ( t ) } } { 8 } \\ & = \left ( s ^ { \prime } + H ( t ) \right ) \wedge \frac { 3 } { 8 },$$

since H t ( ) ∈ ( 1 2 , 3 4 ). Then, the upper bound is a consequence of the pointwise exponent obtained in Proposition 3.2 and the following property of the 2-microlocal frontier:

$$\forall s ^ { \prime } \in \mathbf R ; \ \sigma _ { X, t } ( s ^ { \prime } ) \leq \liminf _ { s \to t } \widetilde { \alpha } _ { X, t },$$

proved in the Appendix (Lemma 5.1).

′

Figure 1: 2-microlocal frontiers of multifractional Brownian motion

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0317v1-figure-f09512d7-e4d4-4e0b-b93b-fb74adace17c.jpg)

In the following second example, we illustrate the particular Hölder regularity which may appear at t = 0.

↦

Example 2. Consider the Hurst function H : t → 3 4 + sin t ( 1 t 2 ) , which locally corresponds to the "chirp" function. Then, the 2-microlocal frontier of the multifractional Brownian motion at 0 is equal to

$$\forall s ^ { \prime } \in \mathbf R ; \ \ \sigma _ { X, 0 } ( s ^ { \prime } ) = \left ( s ^ { \prime } + \frac { 3 } { 4 } \right ) \wedge \left ( \frac { s ^ { \prime } } { 3 } + \frac { 7 } { 1 2 } \right ) \wedge \frac { 3 } { 4 }.$$

Figure 1b gives an illustration of this particular and hybrid 2-microlocal frontier.

To obtain this estimate, we need to refine the proof of Theorem 3.1. Let first recall that the 2-microlocal frontier of H at 0 is equal to

$$\forall s ^ { \prime } \in \mathbf R ; \ \ \sigma _ { H, 0 } ( s ^ { \prime } ) = \frac { 1 } { 3 } ( s ^ { \prime } + 1 ).$$

If we consider the regularity of the term ∂ H B v, H ( 0 )( H v -H v ), we know there exists for any ε &gt; 0 a sequence u n → 0 such that | ∂ H B u ,H ( n 0 ) | ≥ u H 0 -ε n . Then, owing to form of the Hurst function, there exists v n → 0 such that | H u ( n ) -H v ( n ) | ≥ | u n -v n | σ ( u n + v n ) -s ′ where σ &lt; σ H, 0 ( s ′ ). Combining the two, we obtain that the 2-microlocal frontier of the component ∂ H B v, H ( 0 )( H v -H v ) is

$$\forall s ^ { \prime } \geq - \frac { 7 } { 1 2 } ; \ \ \sigma ( s ^ { \prime } ) = \sigma _ { H, 0 } ( s ^ { \prime } + H ( 0 ) ) = \left ( \frac { s ^ { \prime } } { 3 } + \frac { 7 } { 1 2 } \right ).$$

Furthermore, we know that m 0 ,H (0) = ∞ , inducing that the other terms in the Taylor expansion are negligible.

Wenote that the 2-microlocal frontier obtained in this example is consistent with the estimate of the Hölder exponents proved in Proposition 3.2.

## 4. Fractal dimensions of the graph and images

In this section, we aim to extend our work on the sample path geometry of mBm by studying the fractal dimension of its graph and the images of sets. As previously, we will be particularly interested in the study of the irregular case in order to generalize the classic results of Peltier and Lévy Véhel [29].

We first need to recall a few notations on fractal dimensions. For any set A ⊂ R d , we respectively denote by dim H A and dim B A (or dim B A ) the Hausdorff and the lower (or upper) Box dimensions of A . In the case dim B A = dim B A , we denote by dim B A the previous quantity. The reader might refer to the book of Falconer [13] for the precise definition and the basic properties of these fractal dimensions.

As we investigate the fractal geometry of the multifractional Brownian motion, which is a purely non-stationary process, it is more relevant to consider the localised fractal dimension. For all A ⊂ R d and x ∈ R d , it is defined by

$$\dim _ { *, x } A = \lim _ { \rho \to 0 } \dim _ { * } \left ( A \cap B ( x, \rho ) \right ),$$

↦

where dim ∗ denotes any fractal dimension. The previous limit always exists as the right term is decreasing when ρ → 0. We note that the map x → dim ∗ ,x A is upper semi-continuous. Finally, we also observe that the fractal dimension can always be retrieve from the localised values: dim ∗ A = sup x ∈ R d dim ∗ ,x A .

In the first part of this section, we investigate the local fractal dimension of the graph of the mBm. For any function f , its graph on the set V will be denoted Gr( f, V ), i.e.

$$\text{Gr} ( f, V ) = \{ ( t, f ( t ) ) \colon t \in V \}.$$

In the case V = R , we will simply use the notation Gr( f ). To simplify, we will write dim ∗ ,t Gr( f ) for the local fractal dimension at ( t, f ( )). The latter is therefore equivalently defined by t

$$\dim _ { *, t } \text{Gr} ( f, V ) \coloneqq \lim _ { \wedge \to \theta } \dim _ { * } \text{Gr} ( f, V \cap B ( t, \rho ) ).$$

$$\rho \rightarrow 0$$

Several deterministic properties satisfied by the local graph dimension are presented in the Appendix 5, including some connections with the 2-microlocal frontier.

The graph dimension of an H 0 -multifractional Brownian motion has already been determined by Peltier and Lévy Véhel [29], proving that with probability one

$$\forall t \in \mathbf R ; \quad \dim _ { H, t } \mathbf G r ( f ) = \dim _ { B, t } \mathbf G r ( f ) = 2 - H ( t ).$$

In the next result, we investigate the extension of this formula to the Box dimension of the graph of a general multifractional Brownian motion.

Theorem 4.1. Suppose X is a multifractional Brownian motion. Then, with probability one, it satisfies

$$\forall t \in \mathbf R \, \sum \{ 0 \} ; \ \underline { \dim } _ { B, t } \, \text{Gr} ( X ) = \left ( 2 - H ( t ) \right ) \vee \underline { \dim } _ { B, t } \, \text{Gr} ( H )$$

and

$$\forall t \in \mathbf R \, \sum \{ 0 \} ; \ \overline { \dim } _ { B, t } \, \text{Gr} ( X ) = \left ( 2 - H ( t ) \right ) \vee \overline { \dim } _ { B, t } \, \text{Gr} ( H ).$$

Hence, the graph has a local Box dimension at t when the previous two quantities are equal.

Proof. We first obtain the upper bound for both the lower and upper Box dimensions. Let t ∈ R and ρ &gt; 0. Owing to Lemma 3.1, with probability one, there exist c 1 ( ω ) and ε &gt; 0 such that for all u, v ∈ B t, ρ ( ),

$$| X _ { u } - X _ { v } | \leq c _ { 1 } | u - v | ^ { H ( t ) - \varepsilon } + c _ { 1 } | H ( u ) - H ( v ) |.$$

Let N δ ( H ) denotes the smallest number of balls of diameter δ necessary to cover the graph of H on the interval [ t -ρ, t + ]. Up to a constant, we can assume that this cover is organised ρ on columns of size δ , as usually presented in the case of a graph covering. Owing to the upper bound of the increment X u -X v , we know that by adding /ceilingleft δ H t ( ) -ε δ -1 /ceilingright balls to every column covering, we obtain a covering of the graph of the process X . Hence, there exists a constant c 2 such that

$$\forall \delta > 0 ; \ \ N _ { \delta } ( X ) \leq c _ { 2 } \left ( \rho \delta ^ { H ( t ) - \varepsilon - 2 } + N _ { \delta } ( H ) \right )$$

and therefore

$$\frac { \log ( N _ { \delta } ( X ) ) } { - \log ( \delta ) } & \leq \frac { \log ( c _ { 2 } \rho \delta ^ { H ( t ) - \varepsilon - 2 } + c _ { 2 } N _ { \delta } ( H ) ) } { - \log ( \delta ) } \\ & \leq \max \left ( \frac { \log ( c _ { 3 } \delta ^ { H ( t ) - \varepsilon - 2 } ) } { - \log ( \delta ) }, \frac { \log ( c _ { 3 } N _ { \delta } ( H ) ) } { - \log ( \delta ) } \right ).$$

Considering the inferior and superior limits of this inequality, we respectively obtain the upper bound of the lower and upper Box dimensions.

To prove the lower bound 2 -H t ( ), we use classic potential methods that have already been presented by Peltier and Lévy Véhel [29] and Herbin et al. [17] for the multifractional Brownian

motion. Briefly, for any interval [ a, b ] ⊂ R + , there exists C &gt; 0 such that E [ X u -X v ] 2 ≥ C u | - | v 2 H , where H = min u ∈ [ a,b ] H u ( ). This estimate is sufficient, using a standard potential argument to show that almost surely,

$$2 - \min _ { u \in [ a, b ] } H ( u ) \leq \dim _ { H } G r ( X, [ a, b ] ) \leq \underline { \dim } _ { B } G r ( X, [ a, b ] ) \leq \overline { \dim } _ { B } G r ( f, [ a, b ] ).$$

Considering a, b ∈ Q , we obtain the inequality with probability one for any rational interval. Then, we note that Equation (4.1) holds as well if we use in the limit rational intervals containing t . Hence, owing to the continuity of H , with probability one and for all t ∈ R , we get 2 -H t ( ) ≤ dim B ,t Gr( X ) and 2 -H t ( ) ≤ dim B ,t Gr( X ).

To prove the remaining lower bound, let us first consider the case m t,H t ( ) = 1. The constant ρ can be chosen small enough such that for all ( v, h ) in neighbourhood of ( t, H ( )), t | ∂ H B v, h ( ) | &gt; c 4 , where c 4 is a positive constant. Furthermore, for all u, v ∈ B t, ρ ( ),

$$X _ { u } - X _ { v } = B ( u, H ( u ) ) - B ( v, H ( v ) ) + \int _ { H ( v ) } ^ { H ( u ) } \partial _ { H } B ( v, h ) \, \mathrm d h.$$

Hence, in the neighbourhood of t ,

$$c _ { 4 } | H ( u ) - H ( v ) | \leq \left | \int _ { H ( v ) } ^ { H ( u ) } \partial _ { H } B ( v, h ) \, d h \right | \leq | X _ { u } - X _ { v } | + c _ { 5 } | u - v | ^ { H ( t ) - \varepsilon }.$$

Using the covering argument described at the beginning of the proof, we obtain that the inequality dim B ,t Gr( H ) ≤ ( 2 -H t ( ) ) ∨ dim B ,t Gr( X ) and dim B ,t Gr( H ) ≤ ( 2 -H t ( ) ) ∨ dim B ,t Gr( X ). Since we already know that the lower and upper Box dimensions are greater than 2 -H t ( ), it induces the expected lower bound.

↦

Finally, let us extend the inequality to any t ∈ R \ { 0 } and any value of m t,H t ( ) . We have already noted that the function t → dim ∗ ,t Gr( H ) is upper semi-continuous for every dimension considered. As a consequence, there exists a deterministic countable set E ⊂ R \ { 0 } such that for all t , there is a sequence ( t n n ) ∈ N ∈ E satisfying

$$\dim _ { *, t } \text{Gr} ( H ) = \lim _ { n \rightarrow \infty } \dim _ { *, t _ { n } } \text{Gr} ( H ),$$

/negationslash

Since for every s ∈ E ∂ , H B s, H s ( ( )) is a centered Gaussian variable with a positive variance, with probability one and for all s ∈ E , ∂ H B s, H s ( ( )) = 0. Hence, as the local dimension is upper semi-continuous,

$$\underline { \dim } _ { B, t } \, \text{Gr} ( X ) & \geq \lim _ { n \in N } \underline { \dim } _ { B, t _ { n } } \, \text{Gr} ( X ) \\ & = \lim _ { n \in N } \left ( 2 - H ( t _ { n } ) \right ) \vee \underline { \dim } _ { B, t _ { n } } \, \text{Gr} ( H ) = \left ( 2 - H ( t ) \right ) \vee \underline { \dim } _ { B, t } \, \text{Gr} ( H ).$$

The argument holds as well for the upper Box dimension, therefore uniformly proving the result.

Interestingly, we have obtained an expression of the Box dimension of the graph which is similar to the behaviour observed on a fractional Brownian motion with variable drift (see the recent works [11, 30]). This common property happens to extend as well to the Hausdorff

dimension of the graph, recently obtained by Peres and Sousi [31] on the fBm with variable drift.

The Hausdorff dimension of the graph is a bit more complicate to analyse. For this purpose, we have to introduce a more general dimension called the parabolic Hausdorff dimension . We first need to define a parabolic metric /rho1 H on R 2 , with H &gt; 0:

$$\varrho _ { H } \left ( ( u, x ) \, ; ( v, y ) \right ) \coloneqq \max \left ( | u - v | ^ { H }, | x - y | \right ).$$

For any set A ⊂ R 2 , we denote by dim H ( A /rho1 ; H ) the parabolic Hausdorff dimension of A . It is defined similarly to the classic Hausdorff dimension using covering balls relatively to the metric /rho1 H , i.e. it corresponds to the infimum of s ≥ 0 for which

$$\lim _ { \delta \to 0 } \inf \left \{ \sum _ { i = 0 } ^ { \infty } \text{diam} ( O _ { i } \, ; \, \varrho _ { H } ) ^ { s } \, \colon ( O _ { i } ) _ { i \in N } \text{ is a $\delta$-cover of $A$} \right \} < \infty$$

Note that given the form of the parabolic metric, a /rho1 H -ball B (( t, x ) , δ ) has the following form: B (( t, x ) , δ ) = [ t -δ 1 /H , t + δ 1 /H ] × [ x -δ, x + ]. δ

The parabolic Hausdorff dimension has first been considered by Taylor and Watson [37] for the study of polar sets of the heat equation. Recently, it has been proved to be useful by Khoshnevisan and Xiao [20] to analyse the geometry of the images of Brownian motion. Finally, as previously outlined, Peres and Sousi [31] have shown it is also a natural idea for the study of fractional Brownian motion with variable drift. Note that the convention used in the latter work is slightly different from ours.

Theorem 4.2. Suppose X is a multifractional Brownian motion. Then, with probability one, it satisfies

$$\forall t \in \mathbf R \, \sum \{ 0 \} ; \quad \dim _ { H, t } \mathbf G r ( X ) = 1 + H ( t ) \left ( \dim _ { H, t } \left ( \mathbf G r ( H ) \, ; \, \varrho _ { H ( t ) } \right ) - 1 \right ).$$

Proof. To begin with, we prove that with probability one,

$$\forall t \in \mathbf R \, \sum \{ 0 \} ; \ \dim _ { H, t } ( \mathbf G r ( X ) \, ; \, \varrho _ { H ( t ) } ) = \dim _ { H, t } ( \mathbf G r ( H ) \, ; \, \varrho _ { H ( t ) } )$$

The sketch of the proof of this equality is similar to Theorem 4.1, as it uses as well an ad-hoc covering argument. Let us set δ &gt; 0, ρ &gt; 0 and γ &gt; dim H ,t ( Gr( H ) ; /rho1 H t ( ) ) . Suppose ( O i ) i ∈ N is a δ -cover of Gr( H, t [ -ρ, t + ]) with of ρ /rho1 H t ( ) -balls and such that

$$\sum _ { i = 0 } ^ { \infty } \delta _ { i } ^ { \gamma } < \infty \ \text{ where } \delta _ { i } \coloneqq \text{diam} ( O _ { i } \, ; \varrho _ { H ( t ) } ) \ \forall i \in \mathbf N.$$

We observe that | X u -X v | ≤ c 1 | u - | v H t ( ) -ε + c 1 | H u ( ) -H v ( ) | for all u, v ∈ B t, ρ ( ). For any i ∈ N , we need at most 2 δ 1 -ε/H t ( ) i · δ -1 i /rho1 H t ( ) -balls of diameter δ i to cover the part of the graph related to the term | u - | v H t ( ) -ε . Hence, from ( O i ) i ∈ N , we construct a δ -cover ( V k ) k ∈ N of Gr( X, t [ -ρ, t + ]) which satisfies ρ

$$\sum _ { k = 0 } ^ { \infty } \text{diam} ( V _ { k } \, ; \, \varrho _ { H ( t ) } ) ^ { s } \leq c _ { 2 } \sum _ { i = 0 } ^ { \infty } \delta _ { i } ^ { - \varepsilon / H ( t ) } \cdot \delta _ { i } ^ { s }.$$

The last series converges when s ≥ γ + ε/H t ( ). Hence, considering the limits ε → 0, ρ → 0 and γ → dim H ,t ( Gr( H ) ; /rho1 H t ( ) ) , we obtain dim H ,t ( Gr( X ) ; /rho1 H t ( ) ) ≤ dim H ,t ( Gr( H ) ; /rho1 H t ( ) ) .

↦

To prove the other side inequality, we proceed similarly to the proof of Theorem 4.1. We first assume that m t,H t ( ) = 1. In this case, the estimate of the increments presented in Theorem 4.1 and the previous reasoning similarly imply that dim H ,t ( Gr( H ) ; /rho1 H t ( ) ) ≤ dim H ,t ( Gr( X ) ; /rho1 H t ( ) ) . Then, we use the upper semi-continuity of the function t → dim H ,t ( Gr( H ) ; /rho1 H t ( ) ) , obtained in Lemma 5.5, to extend the inequality to any t ∈ R \ { 0 } .

The upper bound of dim H ,t Gr( X ) is then a consequence of the deterministic Lemma 5.4, with H 1 = 1 and H 2 = H t ( ). With probability one, for all t ∈ R , we have dim H Gr( X, t [ -ρ, t + ]) ρ ≤ 1+ H t ( ) ( dim H ( Gr( X, t [ -ρ, t + ]); ρ /rho1 H t ( ) ) -1 . Therefore, when ) ρ → 0, we obtain the first upper bound with probability one: for all t ∈ R \{ } 0 , dim H ,t Gr( X ) ≤ 1+ H t ( ) ( dim H ,t ( Gr( H ) ; /rho1 t ) -1 . )

To prove the lower bound, we make use of potential arguments on the parabolic and euclidean metrics. We first set t ∈ R \ { 0 } . Then, let ρ &gt; 0 and γ &lt; dim H ( Gr( X,B t,ρ ( )) ; /rho1 H t ( ) ) . Using Frostman's lemma applied on the parabolic metric /rho1 t , there exists a probability measure µ on B t, ρ ( ) such that

$$\mathcal { E } _ { \gamma } ( \mu _ { H } ) = \iint _ { B ( t, \rho ) ^ { 2 } } \frac { \mu ( \mathrm d u ) \, \mu ( \mathrm d v ) } { \left ( | u - v | ^ { H ( t ) } + | H ( u ) - H ( v ) | \right ) ^ { \gamma } } < \infty.$$

↦

Suppose µ X denotes the image of the measure µ by the map u → ( u, X u ). To prove the lower on the dimension of the graph of X , we need to study s -energy of µ X , i.e.

$$\mathcal { E } _ { s } ( \mu _ { X } ) = \iint _ { \mathbf R ^ { 4 } } \frac { \mu _ { X } ( \mathrm d x ) \, \mu _ { X } ( \mathrm d y ) } { \| x - y \| ^ { s } } = \iint _ { B ( t, \rho ) ^ { 2 } } \frac { \mu ( \mathrm d u ) \, \mu ( \mathrm d v ) } { \left ( | u - v | ^ { 2 } + | X _ { u } - X _ { v } | ^ { 2 } \right ) ^ { - s / 2 } },$$

where s &gt; 1. Let us set ε &gt; 0. To study the previous energy, we need to split the integral into two components, respectively denoted E 1 s and E 2 s :

$$\mathcal { E } _ { s } ^ { 1 } + \mathcal { E } _ { s } ^ { 2 } = \iint _ { B ( t, \rho ) ^ { 2 } } \frac { \mu ( \mathrm d u ) \, \mu ( \mathrm d v ) } { \left ( | u - v | ^ { 2 } + | X _ { u } - X _ { v } | ^ { 2 } \right ) ^ { - s / 2 } } \left ( \mathbf 1 _ { \left \{ | H ( u ) - H ( v ) | \leq | u - v | ^ { H ( t ) ( 1 - \epsilon ) } \right \} } \\ + \mathbf 1 _ { \left \{ | H ( u ) - H ( v ) | > | u - v | ^ { H ( t ) ( 1 - \epsilon ) } \right \} } \right ).$$

Let us prove E 1 s is finite. For this purpose, we estimate its expectation

$$\mathbb { E } [ \mathcal { E } _ { s } ^ { 1 } ] = \iint _ { B ( t, \rho ) ^ { 2 } } \mathbb { E } \left [ \left ( | u - v | ^ { 2 } + | X _ { u } - X _ { v } | ^ { 2 } \right ) ^ { - s / 2 } \right ] \mathbf 1 _ { \left \{ | H ( u ) - H ( v ) | \leq | u - v | ^ { H ( t ) ( 1 - t ) } \right \} } \mu ( \mathrm d u ) \, \mu ( \mathrm d v ).$$

Using the notation σ 2 u,v := E [ X u -X v ] 2 , the inner term is equal to

$$\mathbb { E } \left [ \left ( | u - v | ^ { 2 } + | X _ { u } - X _ { v } | ^ { 2 } \right ) ^ { - s / 2 } \right ] & = \frac { 2 } { \sigma _ { u, v } \sqrt { 2 \pi } } \int _ { \mathbf R _ { + } } \left ( x ^ { 2 } + | u - v | ^ { 2 } \right ) ^ { - s / 2 } \exp \left ( - \frac { x ^ { 2 } } { 2 \sigma _ { u, v } ^ { 2 } } \right ) \mathrm d x \\ & = c _ { 1 } \int _ { \mathbf R _ { + } } \left ( r \sigma _ { u, v } ^ { 2 } + | u - v | ^ { 2 } \right ) ^ { - s / 2 } r ^ { - 1 / 2 } e ^ { - r / 2 } \, \mathrm d r.$$

Splitting the integral into two terms, and using classic approximations (see e.g. [13]), the previous term is upper bounded by

$$\int _ { 0 } ^ { | u - v | ^ { 2 } / \sigma _ { u, v } ^ { 2 } } | u - v | ^ { - s } r ^ { - 1 / 2 } \, \mathrm d r + \int _ { | u - v | ^ { 2 } / \sigma _ { u, v } ^ { 2 } } ^ { \infty } r ^ { - s / 2 } \sigma _ { u, v } ^ { - s } r ^ { - 1 / 2 } \, \mathrm d r \leq c _ { 2 } | u - v | ^ { 1 - s } \sigma _ { u, v } ^ { - 1 }.$$

We recall that Herbin [15] has proved that

$$\sigma _ { u, v } ^ { 2 } = \mathbb { E } [ X _ { u } - X _ { v } ] ^ { 2 } \asymp | u - v | ^ { 2 H ( t ) } + | H ( u ) - H ( v ) | ^ { 2 }$$

in the neighbourhood of t . Hence, on the domain { | H u ( ) -H v ( ) | ≤ | u - | v H t ( )(1 -ε ) } , we observe that

$$| u - v | & = \left ( | u - v | ^ { 2 H ( t ) ( 1 - \varepsilon ) } \right ) ^ { 1 / 2 H ( t ) ( 1 - \varepsilon ) } \\ & \geq c _ { 3 } \left ( | u - v | ^ { 2 H ( t ) ( 1 - \varepsilon ) } + | H ( u ) - H ( v ) | ^ { 2 } \right ) ^ { 1 / 2 H ( t ) ( 1 - \varepsilon ) } \\ & \geq c _ { 4 } \, \sigma _ { u, v } ^ { 1 / H ( t ) ( 1 - \varepsilon ) }. \\.$$

The expectation of E 1 s is therefore upper bounded by

$$\mathbb { E } [ \mathcal { E } _ { s } ^ { 1 } ] & \leq c _ { 4 } \iint _ { B ( t, \rho ) ^ { 2 } } \sigma _ { u, v } ^ { - 1 + ( 1 - s ) / H ( t ) ( 1 - \varepsilon ) } \, \mu ( \mathrm d u ) \, \mu ( \mathrm d v ) \\ & \leq c _ { 5 } \iint _ { B ( t, \rho ) ^ { 2 } } \left ( | u - v | ^ { H ( t ) } + | H ( u ) - H ( v ) | \right ) ^ { - 1 + ( 1 - s ) / H ( t ) ( 1 - \varepsilon ) } \mu ( \mathrm d u ) \, \mu ( \mathrm d v ).$$

The term E [ E 1 s ] is thus finite when -1 + (1 -s /H t ) ( )(1 -ε ) &gt; -γ , i.e. when

$$s < 1 + ( \gamma - 1 ) \left ( H ( t ) ( 1 - \varepsilon ) \right ) \text{ \ where } \gamma < \dim _ { H } \left ( \text{Gr} ( X, B ( t, \rho ) ) \, ; \, \varrho _ { H ( t ) } \right ).$$

Hence, considering the limits on γ and ε , the upper bound obtained on s corresponds to the expected estimate.

Let now consider the second term E 2 s . We know that m t,H t ( ) = 1 almost surely. Recall that the increment X u -X v has the following form

$$X _ { u } - X _ { v } = B ( u, H ( u ) ) - B ( v, H ( u ) ) + \int _ { H ( v ) } ^ { H ( u ) } \partial _ { H } B ( v, h ) \, \mathrm d h.$$

If ρ is sufficiently small, there exists a positive constant c 1 &gt; 0 such that for all s ∈ B t, ρ ( ) and any h ∈ [min B t,ρ ( ) H u , ( ) max B t,ρ ( ) H u ( )], we have | ∂ H B v, h ( ) | ≥ c 1 . Furthermore, there also exists c 2 &gt; 0 such that

$$\forall u, v \in B ( t, \rho ) ; \ \left | B ( u, H ( u ) ) - B ( v, H ( u ) ) \right | \leq c _ { 2 } | u - v | ^ { H ( t ) ( 1 - \varepsilon / 2 ) }.$$

Using the previous estimates,

$$| X _ { u } - X _ { v } | & \geq \left | \int _ { H ( v ) } ^ { H ( u ) } \partial _ { H } B ( v, h ) \, \mathrm d h \right | - \left | B ( u, H ( u ) ) - B ( v, H ( u ) ) \right | \\ & \geq c _ { 1 } | H ( u ) - H ( v ) | \left ( 1 - c _ { 2 } | H ( u ) - H ( v ) | ^ { - 1 } \cdot | u - v | ^ { H ( t ) ( 1 - \varepsilon / 2 ) } \right )$$

Hence, on the domain { | H u ( ) -H v ( ) | &gt; u | - | v H t ( )(1 -ε ) } , we obtain

$$| X _ { u } - X _ { v } | \geq c _ { 1 } | H ( u ) - H ( v ) | \left ( 1 - c _ { 2 } | u - v | ^ { H ( t ) \varepsilon / 2 } \right ) \geq c _ { 3 } | H ( u ) - H ( v ) |,$$

when ρ is sufficiently small. Therefore, the energy term E 2 s satisfies

$$\mathcal { E } _ { s } ^ { 2 } & \leq c _ { 4 } \iint _ { B ( t, \rho ) } | H ( u ) - H ( v ) | ^ { - s } \mathbf 1 _ { \{ | H ( u ) - H ( v ) | > | u - v | ^ { H ( v ) ( 1 - s ) } \} } \, \mu ( \mathrm d u ) \, \mu ( \mathrm d v ) \\ & \leq c _ { 5 } \iint _ { B ( t, \rho ) ^ { 2 } } \left ( | u - v | ^ { H ( t ) } + | H ( u ) - H ( v ) | \right ) ^ { - s } \mu ( \mathrm d u ) \, \mu ( \mathrm d v ).$$

The last integral is finite when s &lt; γ . But, since γ &gt; 1, we remark that γ (1 -H t ( )) &gt; 1 -H t ( ) and thus γ &gt; 1+ H t ( )( γ -1), showing that E 2 s is always finite when s &lt; 1+( γ -1) H t ( ). Therefore, we have proved that the s -energy E s ( µ X ) is almost surely finite for all s &lt; 1+( γ -1) ( H t ( )(1 -ε ) ) . Considering the limits ε → 0 and γ → dim H ,t ( Gr( X ) ; /rho1 H t ( ) ) , we obtain

$$\dim _ { H, t } \| \text{Gr} ( X ) \geq 1 + H ( t ) \| ( \dim _ { H, t } \| ( \text{Gr} ( H ) \, ; \, \varrho _ { H ( t ) } ) - 1 ) \ \text{ a.s.}$$

↦

The lower bound has been obtained for a fixed t ∈ R \ { 0 } . To extend the latter uniformly on R \ { 0 } , we proceed similarly to Theorem 4.1. Shortly, as proved in Lemma 5.5, the map t → dim H ,t ( Gr( H ) ; /rho1 H t ( ) ) is upper semi-continuous. Hence, there exists a countable set E such that for any t , there is sequence ( t n n ) ∈ N ∈ E satisfying

$$\dim _ { H, t } \left ( \text{Gr} ( H ) \, ; \varrho _ { H ( t ) } \right ) = \lim _ { n \rightarrow \infty } \dim _ { H, t _ { n } } \left ( \text{Gr} ( H ) \, ; \varrho _ { H ( t _ { n } ) } \right ).$$

The equality on the dimension of the graph holds almost surely for all s ∈ E . Therefore, with probability one, for all t ∈ \{ 0 } ,

$$\dim _ { H, t } \text{Gr} ( X ) & \geq \limsup _ { n \in \mathbb { N } } \dim _ { H, t _ { n } } \text{Gr} ( X ) \\ & = 1 + \limsup _ { n \in \mathbb { N } } \, H ( t _ { n } ) \left ( \dim _ { H, t _ { n } } \left ( \text{Gr} ( H ) \, ; \varrho _ { H ( t _ { n } ) } \right ) - 1 \right ) \\ & = 1 + H ( t ) \left ( \dim _ { H, t } \left ( \text{Gr} ( H ) \, ; \varrho _ { H ( t ) } \right ) - 1 \right ).$$

This last inequality concludes the proof, since the upper bound is already uniform on R .

Note that Theorem 4.2 does not contradict the result obtained by Peltier and Lévy Véhel [29] on H 0 -mBms. Indeed, under this assumption H 0 and due to the form of the parabolic metric /rho1 H t ( ) , covering the graph of H in the neighbourhood of t is equivalent to covering a straight line, implying that the local dimension satisfies: dim H ,t ( Gr( H ) ; /rho1 H t ( ) ) = H t ( ) -1 . Hence, the application of Theorem 4.2 proves that with probability one

$$\forall t \in \mathbf R \, \sum \{ 0 \} ; \ \dim _ { H, t } \text{Gr} \, X = 2 - H ( t ),$$

where X is an H 0 -multifractional Brownian motion.

It is interesting to observe that the Hausdorff and Box dimensions of the graph of mBm and fBm with variable drift have a similar form despite the rather different Gaussian structures

of these two processes. Indeed, the first one is a centered process whose incremental variance satisfies

$$\mathbb { E } [ X _ { u } - X _ { v } ] ^ { 2 } \asymp | u - v | ^ { 2 H ( t ) } + | H ( u ) - H ( v ) | ^ { 2 },$$

whereas the latter still has the covariance of fBm, but a mean equal to the drift f ( ). t

As a consequence of these similarities, the examples constructed by Peres and Sousi [31] can also be applied to the mBm, proving that the graph of the latter may have different Hausdorff and Box dimensions.

↦

Example 3. Peres and Sousi [31] have computed the parabolic Hausdorff dimension of the celebrated deterministic construction of McMullen [23]. Hence, there exists a continuous function H : R → [1 / 2 1) such that for all , t ∈ R , ˜ α X,t = log(2) / log(3) := θ ,

$$\dim _ { H, t } ( \text{Gr} ( H ) \, ; \, \varrho _ { 1 / 2 } ) = \log _ { 2 } ( 5 ^ { 2 \theta } + 1 ) \ \text{ and } \ \dim _ { B, t } \text{Gr} ( H ) = 2 - \theta.$$

Then, considering the mBm X with Hurst function H , we first observe that with probability one, for all t &gt; 0, ˜ α X,t = θ . Furthermore, as a consequence of the properties of H , at every t &gt; 0 such that H t ( ) = 1 / 2,

$$\dim _ { H, t } G r ( X ) = \frac { 1 + \log _ { 2 } \left ( 5 ^ { 2 \theta } + 1 \right ) } { 2 } < \dim _ { B, t } G r ( X ) = 2 - \theta.$$

Hence, at these particular times, the local Hausdorff dimension is strictly smaller that the Box one. In addition, as shown by Peres and Sousi [31], the local Hausdorff dimension is such that dim H ,t Gr( X &gt; ) max dim { H ,t Gr( H , ) 3 / 2 } , meaning the simple and straightforward lower bound on the Hausdorff dimension is not optimal.

This example illustrates the kind of unusual fractal geometry the mBm can present. It clearly shows that the parabolic Hausdorff dimension is the relevant concept, as in this case, dim H ,t Gr( H ) happens to be insufficient to characterise entirely the Hausdorff dimension of the graph. We also note that in contrary to fBm, the fractal dimension of the graph can be completely unrelated to the local Hölder exponent.

As presented in the following second example, more classic sample path properties also occur.

Example 4. Suppose H is a sample path of a fractional Brownian motion with Hurst exponent α . Up to a rescaling, it is assumed to be valued in the interval (0 , 1). Then, using some classic estimates (see e.g. [13]), we easily obtain that

$$\forall t \in \mathbf R ; \ \dim _ { H, t } \left ( \mathbf G r ( H ) \, ; \, \varrho _ { H ( t ) } \right ) = \max \left \{ \frac { 1 } { H ( t ) }, 1 + \frac { 1 - \alpha } { H ( t ) } \right \}.$$

Then, owing to Theorem 4.2, a mBm X with Hurst function H satisfies with probability one:

$$\forall t \in { \mathbf R } \, \left \{ 0 \right \} ; \quad \dim _ { H, t } { \mathbf G } ( X ) = \dim _ { B, t } { \mathbf G } ( X ) = \max \left \{ 2 - \alpha, 2 - H ( t ) \right \} = 2 - \widetilde { \alpha } _ { X, t }.$$

This example displays a more common fractal geometry, where the Hausdorff and Box dimensions are both equal to 2 -˜ α X,t . Nevertheless, the local Hölder exponent may still have an irregular behaviour, depending on the value of H t ( ).

In the second part of this section, we are interested in studying the fractal dimension of an image X F ( ), where F can be any fractal set. A large literature has already investigated the

question of the geometry of images of Gaussian processes. It is a well-known result (see e.g. [19]) that for any fixed Borel set F ∈ (0 , ∞ ),

$$\dim _ { H } B ^ { H } ( F ) = \min \left \{ 1, \frac { \dim _ { H } F } { H } \right \} \text{ a.s.}$$

Mountford [27], Khoshnevisan et al. [21], Wu and Xiao [38] have generalized the previous result to ( N,d )-(fractional) Brownian sheet, and in particular, obtained uniform results on F under some assumption on the dimensions N and d .

Adapting classic arguments from Kahane [19], we can easily extend Equation (4.3) to an H 0 -mBm, i.e. for any fixed Borel set F ∈ (0 , ∞ ) and all t ∈ (0 , ∞ )

$$\dim _ { H, X ( t ) } X ( F ) = \min \left \{ 1, \frac { \dim _ { H, t } F } { H ( t ) } \right \} \text{ a.s.}$$

In the irregular case, the statement is not as simple and straightforward. Furthermore, the parabolic Hausdorff dimension again happens to be the relevant tool to describe the geometry of X F ( ).

Theorem 4.3. Suppose X is a multifractional Brownian motion. Then, with probability one,

$$\forall t \in \mathbf R \, \sum \{ 0 \} ; \ \dim _ { H, X ( t ) } X ( F ) = \min \left \{ 1, \dim _ { H, t } \left ( \mathbf G r ( H, F ) \, ; \varrho _ { t } \right ) \right \}.$$

Proof. The proof of this result is similar to the work of Peres and Sousi [31] on the drifted fBm. Let us first remark that

$$\dim _ { H, X ( t ) } X ( F ) \leq \dim _ { H, t } \text{Gr} ( X, F ).$$

↦

This inequality is a simple consequence of the property of the Hausdorff under the application of Lipschitz function (here, the projection p : ( t, X ( )) t → X t ( )). Then, using Lemma 2.1 from [31], we obtain dim H ,t Gr( X,F ) ≤ dim H ,t ( Gr( H,F ) ; /rho1 t ) , thus proving the upper bound.

To obtain the lower bound, we use a classic potential argument. Owing to Frostman's lemma, for any ρ &gt; 0 and any γ &lt; dim H ( Gr( H,F ∩ B t, ρ ( )) ; /rho1 H t ( ) ) , there exists a probability measure µ on B t, ρ ( ) such that

$$\iint _ { B ( t, \rho ) ^ { 2 } } \frac { \mu ( \mathrm d u ) \, \mu ( \mathrm d v ) } { \left ( | u - v | ^ { H ( t ) } + | H ( u ) - H ( v ) | \right ) ^ { \gamma } } < \infty.$$

Recall, as proved by Herbin [15] the covariance of the mBm has the following form

$$\sigma _ { u, v } ^ { 2 } \coloneqq \mathbb { E } [ X _ { u } - X _ { v } ] ^ { 2 } \asymp | u - v | ^ { 2 H ( t ) } + | H ( u ) - H ( v ) | ^ { 2 }$$

↦

Suppose µ X denotes the image of the measure µ under the map u → X u . Then, for any s ∈ (0 , 1), the expectation of its s -energy is equal to

$$\mathbb { * } _ { 1 } & \sim \mathbb { * } _ { 2 } \\ \mathbb { E } \left [ \iint _ { \mathbf R ^ { 2 } } \frac { \mu _ { X } ( \mathrm d x ) \, \mu _ { X } ( \mathrm d y ) } { \| x - y \| ^ { s } } \right ] & = \iint _ { \mathbf R ^ { 2 } } \mathbb { E } [ | X _ { u } - X _ { v } | ^ { - s } ] \mu ( \mathrm d u ) \, \mu ( \mathrm d v ) \\ & = c _ { 2 } \iint _ { \mathbf R ^ { 2 } } \sigma _ { u, v } ^ { - 2 s } \, \mu ( \mathrm d u ) \, \mu ( \mathrm d v )$$

Therefore, owing to the form of the covariance, the integral is finite when s ≤ γ . Hence, considering the limit ρ → 0, we obtain dim H ,X t ( ) X F ( ) ≥ dim H ,t ( Gr( H,F ) ; /rho1 t ) , which concludes the proof.

The local parabolic Hausdorff dimension dim H ,t ( Gr( H,F ) ; /rho1 t ) might not seem to be a very intuitive notion at first. To have a better intuition of the result described in Theorem 4.3, we may first note, similarly to Theorem 4.2, it does not contradict Equation (4.4) in the case of an H 0 -multifractional Brownian motion. Indeed, when H t ( ) &lt; α ˜ H,t , due to the property of the parabolic metric /rho1 H t ( ) , it is completely equivalent to cover F and Gr( H,F ) in the neighbourhood of t . Hence, in this case

$$\dim _ { H, t } \left ( \text{Gr} ( H, F ) \, ; \, \varrho _ { t } \right ) & = \lim _ { \rho \rightarrow 0 } \dim _ { H } \left ( \text{Gr} ( H, F \cap B ( t, \rho ) ) \, ; \, \varrho _ { t } \right ) \\ & = \lim _ { \rho \rightarrow 0 } H ( t ) ^ { - 1 } \dim _ { H } \left ( F \cap B ( t, \rho ) \right ) = H ( t ) ^ { - 1 } \dim _ { H, t } F.$$

The next example shows that different behaviours may occur when the mBm is irregular.

Example 5. Let H be a trajectory of fractional Brownian motion defined on the interval [0 , 1], with Hurst index α ∈ (0 1). , Up to a rescaling and a translation, the sample path satisfies H ([0 , 1]) ⊂ (0 , 1).

We first assume that α &lt; H t ( ). Owing to Equation (4.3), there exists a Borel set F such that dim H ,H t ( ) H F ( ) = α -1 dim H ,t F . Furthermore, since the proof of this statement is based on Frostman's lemma, there exists a probability measure on F ∩ B t, ρ ( ) such that ∫∫ R 2 | H u ( ) -H v ( ) | -γ µ (d u µ ) (d v ) &lt; ∞ where γ &lt; dim H ,H t ( ) H F ( ). Therefore, still using Frostman's lemma, we obtain that α -1 dim H ,t F ≤ dim H ,t ( Gr( H,F ) ; /rho1 t ) . The other side inequality is simply obtained by observing that ˜ α X,t = α , inducing the upper bound using classic covering arguments. Hence, we obtain

$$\dim _ { H, X ( t ) } X ( F ) = \dim _ { H, H ( t ) } H ( F ) = \frac { 1 } { \alpha } \dim _ { H, t } F \text{ \ a.s. }$$

For the second example, we consider the set F = H -1 { 1 2 } , which might be assumed to be non-empty. Monrad and Pitt [26] have proved that dim H F = 1 -α . In this case, as the Hurst function H is constant on the set F , we easily prove that the dimension of Gr( H,F ) is characterised by the dimension on F . Therefore,

$$\dim _ { H, X ( t ) } X ( F ) = \dim _ { H, t } ( \text{Gr} ( H, F ) \, ; \, \varrho _ { t } ) = \frac { 1 } { H ( t ) } \dim _ { H, t } F = \frac { 1 - \alpha } { H ( t ) } \text{ \ a.s. } \text{ \quad (4.6)}$$

The two previous examples show that the estimate obtained in Theorem 4.3 can lead to very different results, depending on the properties of the set F chosen. In particular, we observe that the Hausdorff geometry of images displays a richer structure than the local Hölder regularity. Indeed, in the irregular case, we know that the later only depends on H ( ) · sample path properties, whereas Equations (4.5) and (4.6) show a more complex geometry of the image X F ( ). In the first case, we observe that the formula depends only on the local exponent of H and the Hausdorff dimension of F , whereas in the second one, Equation (4.6) presents a case where the Hausdorff dimension follows the classic formula (4.3), even though the mBm is irregular at t .

To conclude this section, we discuss the question of obtaining the Hausdorff dimension of the level sets X -1 ( { x } ). This problem has already been investigated by Boufoussi et al. [9] for the H 0 -multifractional Brownian motion, proving in particular that for all t ∈ R \ { 0 } ,

$$\dim _ { H, t } X ^ { - 1 } ( \{ X _ { t } \} ) = 1 - H ( t ) \text{ a.s.}$$

It is quite natural to wonder if the previous result can be extended to the irregular mBm. The upper bound can easily extended, since it is a direct consequence of the Hölder regularity of the process at t . In addition, we note that the proof from [9] of the lower bound stands as well in the irregular case. As a consequence, we easily obtain bounds for all t ∈ R \ { 0 } ,

$$1 - H ( t ) \leq \dim _ { H, t } X ^ { - 1 } ( \{ X _ { t } \} ) \leq 1 - \min ( H ( t ), \widetilde { \alpha } _ { H, t } ) \ \text{ a.s.}$$

Lower and upper bounds are not equal in the irregular case, and the remaining question is whether the lower one can be improved. The sketch of the proof of this inequality is rather classic and relies on the study of the Hölder continuity of the local time (see e.g. the seminal works of Berman [7] and Geman and Horowitz [14] on the subject). The key point to prove the Hölder regularity is a property called the local nondeterminism (LND) . In the case of the mBm, it states that for every n ∈ N , there exists c n &gt; 0 such that for all ( t 1 , . . . , t n , t ) ∈ ( ε, ∞ ) n +1

$$\text{Var} ( X _ { t } | X _ { t _ { 1 } }, \dots, X _ { t _ { n } } ) \geq c _ { n } \min _ { 1 \leq k \leq n } | t - t _ { k } | ^ { 2 H ( t ) }.$$

To obtain a better lower bound for the local dimension of the level set, we would need to improve the exponent in the right term. Unfortunately, the next lemma proves that it is optimal under a mild assumption on the Hurst function.

Lemma 4.1. Let t ∈ R and suppose there exists a sequence s i → t such that H s ( i ) = H t ( ) . Then, for any n ∈ N , any ε &gt; 0 and any c n &gt; 0 , there exists ( t 1 , . . . , t n ) ∈ R n such that

$$\text{Var} ( X _ { t } \left | X _ { t _ { 1 } }, \dots, X _ { t _ { n } } ) \leq c _ { n } \min _ { 1 \leq k \leq n } \left | t - t _ { k } \right | ^ { 2 H ( t ) - \varepsilon }.$$

Proof. For all ( t 1 , . . . , t n ) ∈ R n and any k ∈ { 1 , . . . , n } , we note that Var( X t | X ,...,X t 1 t n ) ≤ Var( X t | X t k ) ≤ E [ ( X t -X t k ) 2 ] . Using elements t k ∈ { s i } i ∈ N and the estimate of the variance of the increments, we obtain that

$$\text{Var} ( X _ { t } | X _ { t _ { 1 } }, \dots, X _ { t _ { n } } ) \leq C \min _ { 1 \leq k \leq n } | t - t _ { k } | ^ { 2 H ( t ) }.$$

Let set ρ &gt; 0. We may assume that t k ∈ B t, ρ ( ) for every k ∈ { 1 , . . . , n } , implying that

$$\text{Var} ( X _ { t } | X _ { t _ { 1 } }, \dots, X _ { t _ { n } } ) \leq C \rho ^ { \varepsilon } \min _ { 1 \leq k \leq n } | t - t _ { k } | ^ { 2 H ( t ) - \varepsilon }.$$

Using ρ sufficiently small, it proves that the property of local nondeterminism does not stand with the exponent 2 H t ( ) -ε , for any ε &gt; 0.

Hence, the classic exponent 2 H t ( ) in the LND property of the mBm is most of the time optimal. Intuitively, it can be understood as following: in the neighbourhood of t , the Hurst function can influence the Hölder regularity, but is still a deterministic component and thus does

not bring any randomness which could increase the conditional variance. As the consequence, the latter is only a result of the local fBm form of the multifractional Brownian motion.

Therefore, the estimate of the Hausdorff dimension of the level sets can not be improved using the classic methods based on the LND property, and it remains an open question to determine the precise form of the fractal dimension in terms of the geometric properties of the Hurst function H .

## 5. Appendix: a few deterministic properties

We gather in the appendix a few deterministic properties related to 2-microlocal analysis and local graph dimensions and which have been used through the paper. We begin with a simple upper bound on the 2-microlocal which extends the classic property of lower semi-continuity on the local Hölder exponent.

Lemma 5.1. Suppose f is a continuous function and t ∈ R such that lim inf u → t ˜ α f,t &lt; 1 . Then, the 2-microlocal frontier of f at t satisfies

$$\forall s ^ { \prime } \in \mathbf R ; \ \ \sigma _ { f, t } ( s ^ { \prime } ) \leq \liminf _ { u \to t } \widetilde { \alpha } _ { f, u }.$$

Proof. Let us set ρ &gt; 0 and ε &gt; 0. There exists s ∈ B t, ρ ( ) such that ˜ α f,s ≤ lim inf u → t ˜ α f,u + . ε Then, we can find u , v n n → s such that

$$| f ( u _ { n } ) - f ( v _ { n } ) | \geq n \, | u _ { n } - v _ { n } | ^ { \tilde { \alpha } _ { f, s } + \varepsilon }.$$

Therefore, for any s ′ ∈ R , we have

$$\frac { | f ( u _ { n } ) - f ( v _ { n } ) | } { | u _ { n } - v _ { n } | ^ { \widetilde { \alpha } _ { f, s } + \varepsilon } \left ( | u _ { n } - t | + | v _ { n } - t | \right ) ^ { - s ^ { \prime } } } \geq \frac { n } { \left ( | u _ { n } - t | + | v _ { n } - t | \right ) ^ { - s ^ { \prime } } } \longrightarrow _ { n \to + \infty } + \infty,$$

$$\iota \, \sigma _ { f, t } ( s ^ { \prime } ) \leq \liminf _ { u \to t } \widetilde { \alpha } _ { f, u } + \varepsilon, \, \text{for any } \varepsilon > 0.$$

proving that σ f,t ( s ) ≤ lim inf u → t ˜ α f,u + , for any ε ε &gt;

In the next two lemmas, we present an extension of a well-known property connecting the local dimension of the graph to the local Hölder exponent: dim B ,t Gr( f ) ≤ 2 -˜ α f,t .

Lemma 5.2. Suppose f is a continuous function. Then, for all t ∈ R ,

$$\underline { \dim } _ { B, t } \, \text{Gr} ( f ) \leq \overline { \dim } _ { B, t } \, \text{Gr} ( f ) \leq 2 - \left ( \sigma _ { f, t } ( 1 ) \wedge 1 \right ).$$

Proof. We need to construct a local cover of the graph which is sufficiently efficient. Let first set t ∈ R and σ &lt; σ f,t (1) ∧ 1. Then, there exists ρ &gt; 0 and C &gt; 0 such that for all u, v ∈ B t, ρ ( ),

$$| f ( u ) - f ( v ) | \leq C | u - v | ^ { \sigma } \left ( | u - t | + | v - t | \right ) ^ { - 1 }.$$

The cover is constructed similarly to the original proof. Hence, if N δ denotes the number of

δ -boxes necessary to cover Gr( f, [ t -ρ, t + ]), we have ρ

$$N _ { \delta } \leq C \delta ^ { - 1 } + \sum _ { k = - \lceil \rceil \rho / \delta \rceil } ^ { \lceil \rho / \delta \rceil } \sup _ { u, v \in t + [ k \delta, ( k + 1 ) \delta \rceil } | f ( u ) - f ( v ) | \cdot \delta ^ { - 1 }$$

$$N _ { \delta } & \leq C \delta ^ { - 1 } + \sum _ { k = - \lceil \rho / \delta \rceil } ^ { \lceil \rho / \delta \rceil } \sup _ { u, v \in t + \lfloor k \delta, ( k + 1 ) \delta \rceil } | f ( u ) - f ( v ) | \cdot \delta ^ { - 1 } \\ & \leq C \delta ^ { - 1 } + \sum _ { k = - \lceil \rho / \delta \rceil } ^ { \lceil \rho / \delta \rceil } C \delta ^ { \sigma } ( k \delta ) ^ { - 1 } \cdot \delta ^ { - 1 } \\ & = C \delta ^ { - 1 } + 2 C \delta ^ { \sigma - 2 } \sum _ { k = 1 } ^ { \lceil \rho / \delta \rceil } k ^ { - 1 } \leq 2 C \delta ^ { \sigma - 2 } ( 1 + \log \lceil \rho / \delta \rceil ). \\ \cdot 1. \text{ Then}.$$

since 2 -σ &gt; 1. Then,

$$\overline { \dim } _ { B } \, \text{Gr} ( f, [ t - \rho, t + \rho ] ) = \lim \sup _ { \delta \to 0 } \frac { N _ { \delta } } { - \log \delta } \leq 2 - \sigma,$$

for all σ &lt; σ f,t 0 (1) ∧ 1, which proves to the expected inequality.

Lemma 5.3. Suppose f is a continuous function. Then, for all t ∈ R

$$\dim _ { H, t } \text{Gr} ( f ) \leq 2 - \left ( \sigma _ { f, t } ( + \infty ) \wedge 1 \right ).$$

Proof. As the estimate we want to obtain differs from the previous lemma, we need to adopt a constructing procedure which is slightly different. Let t ∈ R and σ &lt; σ f,t (+ ∞ ∧ ) 1. Then, there exist s ′ ≥ 0, ρ &gt; 0 and C &gt; 0 such that for all u, v ∈ B t, ρ ( )

$$| f ( u ) - f ( v ) | \leq C | u - v | ^ { \sigma } \left ( | u - t | + | v - t | \right ) ^ { - s ^ { \prime } }.$$

Let us set δ &gt; 0 and γ &gt; 1. The δ -cover of Gr( f, [ t -ρ, t + ]) is constructed as following. There ρ exists n ∈ N such that for all k ≥ n , k -γ &lt; δ . To cover the set [ t -ρ, t + ], we first choose ρ successive intervals ( I k ) k ≥ n of size k -γ in the neighbourhood of t . The rest is simply covered by intervals of size δ . The first type of intervals cover a set whose length is equal to:

$$\sum _ { k = n } ^ { + \infty } k ^ { - \gamma } \geq \int _ { n } ^ { + \infty } x ^ { - \gamma } \mathrm d x = C _ { \gamma } n ^ { 1 - \gamma }.$$

Hence, the number M δ of intervals of size δ which are necessary is upper bounded by

$$M _ { \delta } \leq \delta ^ { - 1 } ( 2 \rho - C _ { \gamma } n ^ { 1 - \gamma } ).$$

Let now consider the cover ( O i ) i ∈ N of Gr( f, [ t -ρ, t + ]) which can be defined from the previous ρ intervals (still using the same idea of covering by columns). It satisfies for any s &gt; 0

$$\sum _ { i \in \mathbb { N } } | O _ { i } | ^ { s } \leq \sum _ { k = n } ^ { + \infty } ( k ^ { - \gamma } ) ^ { s } \cdot k ^ { \gamma } \sup _ { u, v \in I _ { k } } | f ( u ) - f ( v ) | + \sum _ { k = 1 } ^ { M _ { \delta } } \delta ^ { s } \cdot \delta ^ { - 1 } \sup _ { u, v \in \widetilde { I } _ { k } } | f ( u ) - f ( v ) |.$$

Let us estimate the first term,

$$\sum _ { k = n } ^ { + \infty } ( k ^ { - \gamma } ) ^ { s } \cdot k ^ { \gamma } \sup _ { u, v \in I _ { k } } | f ( u ) - f ( v ) | \leq \sum _ { k = n } ^ { + \infty } ( k ^ { - \gamma } ) ^ { \sigma + s - 1 } | x _ { k } - t | ^ { - s ^ { \prime } }$$

$$\sum _ { k = n } ^ { \sum _ { k } ^ { \prime } } ( \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \mu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^{\ 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2} \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nu ^ { 2 } \nolimits _ { k = n } ^ { \sum _ { k = n } ^ { + \infty } ( k ^ { - 1 } ) ^ { \gamma ( \sigma + s - 1 ) - s ^ { \prime } ( \gamma - 1 ) }.$$

where we have used that d t, I ( k ) ∼ k 1 -γ . The previous series converges if s satisfies γ σ ( + s -1) -s ′ ( γ -1) &gt; 1, i.e. if s &gt; 1 -σ + s ′ (1 -γ -1 ) + γ -1 .

Let now consider the second term,

$$\sum _ { k = 1 } ^ { M _ { \delta } } \delta ^ { s } \cdot \delta ^ { - 1 } \sup _ { u, v \in \widetilde { I } _ { k } } | f ( u ) - f ( v ) | \leq \sum _ { k = 1 } ^ { M _ { \delta } } \delta ^ { \sigma + s - 1 } \, | \widetilde { x } _ { k } - t | ^ { - s ^ { \prime } }$$

$$\leq \delta ^ { \sigma + s - 1 } \sum _ { \L } ^ { M _ { \delta } } \left ( C _ { \gamma } n ^ { 1 - \gamma } + k \delta \right ) ^ { - s ^ { \prime } }.$$

$$\sum _ { k = 1 } ^ { \mu } \mu ^ { \mu } \sum _ { \substack { \mu \mu \nu \\ u, v \in \widetilde { I } _ { k } } } \nu \sum _ { \substack { \nu \nu \\ k = 1 } } ^ { \mu } \mu ^ { \mu \nu } \mu ^ { \mu } \\ \leq \delta ^ { \sigma + s - 1 } \sum _ { k = 1 } ^ { M _ { \delta } } ( C _ { \gamma } n ^ { 1 - \mu }$$

The sum can upper bounded by

$$\int _ { 0 } ^ { M _ { \delta } } \left ( C _ { \gamma } n ^ { 1 - \gamma } + x \delta \right ) ^ { - s ^ { \prime } } \mathrm d x & \leq \delta ^ { - 1 } \int _ { 0 } ^ { \rho } \left ( C _ { \gamma } n ^ { 1 - \gamma } + v \right ) ^ { - s ^ { \prime } } \mathrm d v \\ & \leq c _ { 2 } \, \delta ^ { - 1 } \cdot n ^ { ( \gamma - 1 ) ( s ^ { \prime } - 1 ) } \\ & \leq c _ { 3 } \, \delta ^ { - 1 } \cdot \delta ^ { - ( \gamma - 1 ) ( s ^ { \prime } - 1 ) / \gamma }.$$

Combining the previous estimate, we obtain

$$\sum _ { k = 1 } ^ { M _ { \delta } } \delta ^ { s } \cdot \delta ^ { - 1 } \sup _ { u, v \in \widetilde { I } _ { k } } | f ( u ) - f ( v ) | \leq c _ { 4 } \, \delta ^ { \sigma + s - 2 - ( \gamma - 1 ) ( s ^ { \prime } - 1 ) / \gamma }.$$

Hence, the latter converges as δ → 0 when s &gt; 2 -σ +( γ -1)( s ′ -1) /γ . Therefore, combining the two conditions we have obtained, the Hausdorff dimension satisfies

$$\dim _ { H } \text{Gr} ( f, [ t - \rho, t + \rho ] ) \leq \min \{ 1 - \sigma + s ^ { \prime } ( 1 - \gamma ^ { - 1 } ) + \gamma ^ { - 1 }, 2 - \sigma + ( \gamma - 1 ) ( s ^ { \prime } - 1 ) / \gamma \},$$

for all γ &gt; 1 and s ′ ∈ R . Considering the limit γ → 1, we obtain dim H Gr( f, [ t -ρ, t + ]) ρ ≤ 2 -σ , which implies the expected inequality.

↦

$$\forall s ^ { \prime } \in \mathbf R ; \ \ \sigma _ { f, 0 } ( s ^ { \prime } ) = \frac { s ^ { \prime } + \alpha } { 1 + \beta }.$$

We note that the two new upper bounds presented in the previous Lemmas are consistent with the original one, since ˜ α f,t = σ f,t (0) ≤ σ f,t (1) ≤ σ f,t (+ ∞ ). Furthermore, it proves to be the optimal ones for a classic deterministic function called the chirp function : f : x → | x | α sin( | x | -β ), with α, β &gt; 0. As proved by Echelard [12], its 2-microlocal frontier at 0 is equal to

It is an easy calculus to check that the local Hausdorff and Box dimensions at zero are given by

$$\dim _ { H, 0 } \, \text{Gr} ( f ) = 1 \quad \text{and} \quad \dim _ { B, 0 } \, \text{Gr} ( f ) = 2 - \left ( \frac { 1 + \alpha } { 1 + \beta } \right ) \wedge 1,$$

Hence, they are equal to the upper bounds presented in Lemmas 5.2 and 5.3, on the contrary to the usual Hölder estimate 2 -α/ (1 + β ) which is not optimal.

To end this section, we prove two useful lemmas related to the parabolic Hausdorff dimension. The first one constitutes an extension of an inequality proved by Peres and Sousi [31] (Lemma 2 1). .

Lemma 5.4. For every A ⊂ R 2 and all H &gt; H 1 2 &gt; 0 ,

$$\dim _ { H } ( A \, ; \, \varrho _ { H _ { 2 } } ) + \frac { 1 } { H _ { 2 } } - \frac { 1 } { H _ { 1 } } \ \leq \ \dim _ { H } ( A \, ; \, \varrho _ { H _ { 1 } } ) \ \leq \ 1 + \frac { H _ { 2 } } { H _ { 1 } } ( \dim _ { H } ( A \, ; \, \varrho _ { H _ { 2 } } ) - 1 ).$$

↦

As a consequence, if H ( ) · is a positive continuous function, the map t → dim H ( A /rho1 ; H t ( ) ) is also continuous.

Proof. Suppose A ⊂ R 2 and H &gt; H 1 2 . Let us first prove the upper bound on dim H ( A /rho1 ; H 1 ) . For any γ &gt; dim H ( A /rho1 ; H 2 ) , there exists a /rho1 H 2 -cover ( O i ) i ∈ N of A such that

$$\sum _ { i = 0 } ^ { \infty } \delta _ { i } ^ { \gamma } < \infty \ \text{ where } \ \delta _ { i } \coloneqq \text{diam} ( O _ { i } \, ; \, \varrho _ { H _ { 2 } } ) \ \forall i \in \mathbf N.$$

We denote by ρ i the quantity δ H /H 1 2 i . Without any loss of generality, we may assume that for every i ∈ N , O i is a rectangle of size δ 1 /H 2 i × δ i (as a simple consequence of the definition of /rho1 H 2 ). We want to construct a /rho1 H 1 -cover of A from ( O i ) i ∈ N . Since ρ 1 /H 1 i = δ 1 /H 2 i , we need to split O i along the space axis to obtain /rho1 H 1 -balls. More precisely, there will be at most 2 δ /ρ i i = 2 δ 1 -H /H 1 2 i resulting /rho1 H 1 -balls of diameter ρ i . Hence, we obtain a /rho1 H 1 -cover ( V k ) k ∈ N of A such that

$$\sum _ { k = 0 } ^ { \infty } \text{diam} \left ( V _ { k } \, ; \varrho _ { H _ { 1 } } \right ) ^ { s } \leq 2 \sum _ { i = 0 } ^ { \infty } \delta _ { i } ^ { 1 - H _ { 1 } / H _ { 2 } } \cdot \delta _ { i } ^ { s H _ { 1 } / H _ { 2 } }.$$

The last sum is finite when 1 + H /H 1 2 ( s -1) &gt; γ , i.e. when s &gt; 1+( γ -1) H /H 2 1 . Considering γ → dim H ( A /rho1 ; H 2 ) , we obtain the expected upper bound.

The lower bound is obtained using a similar reasoning. For any s &gt; dim H ( A /rho1 ; H 1 ) , there exists a /rho1 H 1 -cover ( V k ) k ∈ N of A such that

$$\sum _ { k = 0 } ^ { \infty } \rho _ { k } ^ { s } < \infty \ \text{ where } \ \rho _ { k } \coloneqq \text{diam} \left ( V _ { k } \, ; \, \varrho _ { H _ { 1 } } \right ) \ \forall k \in \mathbf N.$$

Similarly, we may assume that V k is a rectangle of size ρ 1 /H 1 k × ρ k for every k ∈ N . As H &gt; H 1 2 , we note that ρ 1 /H 2 &lt; ρ 1 /H 1 k . Hence, to obtain a /rho1 H 2 -cover of A , we split the previous balls along the time axis, in at most 2 ρ 1 /H 1 -1 /H 2 k sub-balls. Therefore, we obtain a /rho1 H 2 -cover ( O i ) i ∈ N of A such that

$$\sum _ { i = 0 } ^ { \infty } \text{diam} \left ( O _ { i } \, ; \, \varrho _ { H _ { 2 } } \right ) ^ { \gamma } \leq 2 \sum _ { k = 0 } ^ { \infty } \rho _ { k } ^ { 1 / H _ { 1 } - 1 / H _ { 2 } } \cdot \rho _ { k } ^ { \gamma }.$$

↦

This sum is finite when γ +1 /H 1 -1 /H 2 &gt; s , inducing the lower bound on dim H ( A /rho1 ; H 1 ) . The continuity of the map t → dim H ( A /rho1 ; H t ( ) ) is straightforward as we observe that lower and upper bounds converge to dim H ( A /rho1 ; H 1 ) when H 2 → H 1 .

The previous lemma helps us to extend the property of upper semi-continuity to the local parabolic Hausdorff dimension.

↦

Lemma 5.5. Suppose A ⊂ R 2 and H ( ) · is a positive continuous function. Then, the map t → dim H ,t ( A /rho1 ; H t ( ) ) is upper semi-continuous, i.e. for all t ∈ R

$$\dim _ { H, t } ( A \, ; \, \varrho _ { H ( t ) } ) \geq \lim _ { s \to t } \sup _ { s \to t } \dim _ { H, s } ( A \, ; \, \varrho _ { H ( s ) } ).$$

Proof. For any t ∈ R , the local parabolic Hausdorff dimension satisfies

$$\dim _ { H, t } ( A \, ; \varrho _ { H ( t ) } ) \geq \lim _ { s \to t } \sup _ { s \to t } \dim _ { H, s } ( A \, ; \varrho _ { H ( t ) } ).$$

This inequality is straightforward extension of the usual one on the Hausdorff dimension. Let δ &gt; 0. Owing to the uniform bounds obtained in Lemma 5.4, there exists c δ &gt; 0 such that for all s ∈ B t, δ ( ) and all ρ &gt; 0,

$$\left | \dim _ { H } ( A \cap B ( s, \rho ) \, ; \, \varrho _ { H ( t ) } ) - \dim _ { H } ( A \cap B ( s, \rho ) \, ; \, \varrho _ { H ( s ) } ) \right | \leq c _ { \delta } | H ( s ) - H ( t ) |.$$

↦

Hence, considering the limit ρ → 0, we obtain | dim H ,s ( A /rho1 ; H t ( ) ) -dim H ,s ( A /rho1 ; H s ( ) ) | ≤ c δ | H s ( ) -H t ( ) , | proving, jointly with the first inequality, the upper semi-continuity of the map t → dim H ,t ( A /rho1 ; H t ( ) ) .

## References

- [1] A. Ayache. Continuous Gaussian multifractional processes with random pointwise Hölder regularity. J. Theoret. Probab. , 26(1):72-93, 2013.
- [2] A. Ayache and M. S. Taqqu. Multifractional processes with random exponent. Publ. Mat. , 49(2):459-486, 2005.
- [3] A. Ayache, S. Cohen, and J. Lévy Véhel. The covariance structure of multifractional brownian motion, with application to long range dependence. In Proceedings of the Acoustics, Speech, and Signal Processing, 2000 , volume 6, pages 3810-3813, Washington, DC, USA, 2000. IEEE Computer Society.
- [4] A. Ayache, N.-R. Shieh, and Y. Xiao. Multiparameter multifractional brownian motion: local nondeterminism and joint continuity of the local times. Ann. Inst. H. Poincaré Probab. Statist , 2011.
- [5] P. Balança. Fine regularity of Lévy processes and linear (multi)fractional stable motion. Preprint , 2013. arXiv:1302.3140 .
- [6] A. Benassi, S. Jaffard, and D. Roux. Elliptic Gaussian random processes. Rev. Mat. Iberoamericana , 13(1):19-90, 1997.
- [7] S. M. Berman. Local nondeterminism and local times of Gaussian processes. Indiana Univ. Math. J. , 23:69-94, 1973.
- [8] J.-M. Bony. Second microlocalization and propagation of singularities for semilinear hyperbolic equations. In Hyperbolic equations and related topics (Katata/Kyoto, 1984) , pages 11-49. Academic Press, Boston, MA, 1986.

- [9] B. Boufoussi, M. Dozzi, and R. Guerbaz. On the local time of multifractional Brownian motion. Stochastics , 78(1):33-49, 2006.
- [10] B. Boufoussi, M. Dozzi, and R. Guerbaz. Sample path properties of the local time of multifractional Brownian motion. Bernoulli , 13(3):849-867, 2007.
- [11] P. H. A. Charmoy, Y. Peres, and P. Sousi. Minkowski dimension of brownian motion with drift. 2012.
- [12] A. Echelard. Analyse 2-microlocale et application au débruitage . PhD thesis, Université de Nantes, 2007. http://tel.archives-ouvertes.fr/tel-00283008/fr/ .
- [13] K. Falconer. Fractal geometry . John Wiley &amp; Sons Inc., Hoboken, NJ, second edition, 2003. Mathematical foundations and applications.
- [14] D. Geman and J. Horowitz. Occupation densities. Ann. Probab. , 8(1):1-67, 1980.
- [15] E. Herbin. From N parameter fractional Brownian motions to N parameter multifractional Brownian motions. Rocky Mountain J. Math. , 36(4):1249-1284, 2006.
- [16] E. Herbin and J. Lévy Véhel. Stochastic 2-microlocal analysis. Stochastic Process. Appl. , 119(7):2277-2311, 2009.
- [17] E. Herbin, B. Arras, and G. Barruel. From almost sure local regularity to almost sure Hausdorff dimension for Gaussian fields. Preprint , 2012.
- [18] S. Jaffard. Pointwise smoothness, two-microlocalization and wavelet coefficients. Publ. Mat. , 35(1):155-168, 1991. Conference on Mathematical Analysis (El Escorial, 1989).
- [19] J.-P. Kahane. Some random series of functions , volume 5 of Cambridge Studies in Advanced Mathematics . Cambridge University Press, Cambridge, second edition, 1985.
- [20] D. Khoshnevisan and Y. Xiao. Brownian motion and thermal capacity. Preprint , 2012.
- [21] D. Khoshnevisan, D. Wu, and Y. Xiao. Sectorial local non-determinism and the geometry of the Brownian sheet. Electron. J. Probab. , 11:no. 32, 817-843, 2006.
- [22] K. M. Kolwankar and J. Lévy Véhel. A time domain characterization of the fine local regularity of functions. J. Fourier Anal. Appl. , 8(4):319-334, 2002.
- [23] C. McMullen. The Hausdorff dimension of general Sierpiński carpets. Nagoya Math. J. , 96:1-9, 1984.
- [24] M. Meerschaert, D. Wu, and Y. Xiao. Local times of multifractional Brownian sheets. Bernoulli , 14(3):865-898, 2008.
- [25] Y. Meyer. Wavelets, vibrations and scalings , volume 9 of CRM Monograph Series . American Mathematical Society, Providence, RI, 1998.
- [26] D. Monrad and L. D. Pitt. Local nondeterminism and Hausdorff dimension. In Seminar on stochastic processes, 1986 (Charlottesville, Va., 1986) , volume 13 of Progr. Probab. Statist. , pages 163-189. Birkhäuser Boston, Boston, MA, 1987.
- [27] T. S. Mountford. Uniform dimension results for the Brownian sheet. Ann. Probab. , 17(4): 1454-1462, 1989.
- [28] G. C. Papanicolaou and K. Sølna. Wavelet based estimation of local Kolmogorov turbulence. In Theory and applications of long-range dependence , pages 473-505. Birkhäuser Boston, Boston, MA, 2003.
- [29] R. F. Peltier and J. Lévy Véhel. Multifractional brownian motion: Definition and preliminary results. Rapport de recherche INRIA , (2645):1-39, 1995.
- [30] Y. Peres and P. Sousi. Brownian motion with variable drift: 0-1 laws, hitting probabilities and Hausdorff dimension. Math. Proc. Cambridge Philos. Soc. , 153(2):215-234, 2012.
- [31] Y. Peres and P. Sousi. Dimension of fractional brownian motion with variable drift. 2013.
- [32] J. Picard. Representation formulae for the fractional Brownian motion. Séminaire de

Probabilités , XLIII:3-70, 2011.

- [33] S. G. Samko, A. A. Kilbas, and O. I. Marichev. Fractional integrals and derivatives . Gordon and Breach Science Publishers, Yverdon, 1993.
- [34] S. Seuret and J. Lévy Véhel. A time domain characterization of 2-microlocal spaces. J. Fourier Anal. Appl. , 9(5):473-495, 2003.
- [35] S. A. Stoev and M. S. Taqqu. How rich is the class of multifractional Brownian motions? Stochastic Process. Appl. , 116(2):200-221, 2006.
- [36] K. Takashima. Sample path properties of ergodic self-similar processes. Osaka J. Math. , 26(1):159-189, 1989.
- [37] S. J. Taylor and N. A. Watson. A Hausdorff measure classification of polar sets for the heat equation. Math. Proc. Cambridge Philos. Soc. , 97(2):325-344, 1985.
- [38] D. S. Wu and Y. M. Xiao. Dimensional properties of fractional Brownian motion. Acta Math. Sin. (Engl. Ser.) , 23(4):613-622, 2007.
Paul Balança
Some sample path properties of multifractional Brownian motion
The geometry of the multifractional Brownian motion (mBm) is known to present
a complex and surprising form when the Hurst function is greatly irregular.
Nevertheless, most of the literature devoted to the subject considers
sufficiently smooth cases which lead to sample paths locally similar to a
fractional Brownian motion (fBm). The main goal of this paper is therefore to
extend these results to a more general frame and consider any type of
continuous Hurst function. More specifically, we mainly focus on obtaining a
complete characterization of the pointwise H\"older regularity of the sample
paths, and the Box and Hausdorff dimensions of the graph. These results, which
are somehow unusual for a Gaussian process, are illustrated by several
examples, presenting in this way different aspects of the geometry of the mBm
with irregular Hurst functions