Assuming the model $f(\lambda) \sim G\lambda^{1-2H}$, as $\lambda \rightarrow 0 +$, for the spectral density of a covariance stationary process, we consider an estimate of $H \in (0, 1)$ which maximizes an approximate form of frequency domain Gaussian likelihood, where discrete averaging is carried out over a neighbourhood of zero frequency which degenerates slowly to zero as sample size tends to infinity. The estimate has several advantages. It is shown to be consistent under mild conditions. Under conditions which are not greatly stronger, it is shown to be asymptotically normal and more efficient than previous estimates. Gaussianity is nowhere assumed in the asymptotic theory, the limiting normal distribution is of very simple form, involving a variance which is not dependent on unknown parameters, and the theory covers simultaneously the cases $f(\lambda) \rightarrow \infty, f(\lambda) \rightarrow 0$ and $f(\lambda) \rightarrow C \in (0, \infty)$, as $\lambda \rightarrow 0$. Monte Carlo evidence on finite-sample performance is reported, along with an application to a historical series of minimum levels of the River Nile.