This paper develops a rigorous inferential framework for a class of Gaussian stochastic processes driven by white noise with constant drift, whose temporal evolution is governed by a Caputo fractional derivative of order α∈(1/2,1). The model belongs to the family of fractional Volterra processes, where memory is generated by the dynamics themselves rather than by correlated noise. We derive explicit analytical expressions for the mean, variance, and covariance structure of the solution, thereby characterizing in a precise manner how the fractional order α governs both variance growth and the strength of temporal dependence. In particular, the process exhibits correlated increments and a power-law variance scaling of order t2α−1, highlighting the dual role of α as a regularity and memory parameter. Building on this structural analysis, we address the statistical problem of estimating the parameter vector (μ,σ,α) from discrete-time observations. Two complementary procedures are proposed for the estimation of the fractional order: a variance-growth method based on log–log regression of empirical variances, and a wavelet-based estimator exploiting multi-scale scaling properties of the process. For the drift and diffusion parameters (μ,σ), we construct explicit Gaussian pseudo-maximum likelihood estimators derived from the Volterra covariance structure of the increment process. We establish unbiasedness, L2-convergence, strong consistency, and asymptotic normality for all estimators. Furthermore, we derive Berry–Esseen type bounds that quantify the rate of convergence toward the Gaussian law, providing sharp distributional approximations in a genuinely fractional and non-Markovian setting. A Monte Carlo study is carried out, using high-resolution Volterra discretizations, large-scale simulation budgets, covariance-structured linear algebra, and multi-scale diagnostic tools. The numerical experiments confirm the theoretical convergence rates, demonstrate the finite-sample reliability of the estimators, and illustrate the sensitivity of the process dynamics to the fractional order α: smaller values of α produce stronger memory effects and higher variability, while values closer to one lead to smoother and more stable trajectories. The proposed methodology unifies statistical inference for long-memory Gaussian processes with fractional differential stochastic dynamics, offering a coherent analytical and computational framework applicable in areas such as quantitative finance, anomalous diffusion in physics, hydrology, and engineering systems with hereditary effects.
Keddi et al. (Tue,) studied this question.