Let $B^H(t)$ be a fractional Brownian motion with Hurst parameter $H\in (0,1)$. We define fractional Gaussian noise as $X(t)=B^H(t+1)-B^H(t)$. We know the fBm has covariance $R(s,t)=E(B^H(t)B^H(s))=\frac{1}{2}(t^{2H}+s^{2H}-|t-s|^{2H})$. We can build up the covariance of increments of fBm from this and then even further we can ask about:
$$E((X(t)-X(s))(X(v)-X(u))$$
Is there any non horrible way of writing down what this covariance should be? I can write it in terms of a bunch of $R$s but nothing better.
Hint: $\mathbb{E}\left[\left(B_{t}^{H}-B_{s}^{H}\right)\left(B_{u}^{H}-B_{v}^{H}\right)\right]=\frac{1}{2} \mathbb{E}\left[\left(B_{t}^{H}-B_{v}^{H}\right)^{2}+\left(B_{s}^{H}-B_{u}^{H}\right)^{2}-\left(B_{t}^{H}-B_{u}^{H}\right)^{2}-\left(B_{s}^{H}-B_{v}^{H}\right)^{2}\right]=\frac{1}{2}\left(|t-v|^{2 H}+|s-u|^{2 H}-|t-u|^{2 H}-|s-v|^{2 H}\right)$