The stochastic integral is defined as $$u_t = \int_{t-1}^t e^{-\kappa(t-s)}\int_0^s e^{-c(s-r)} \, dW(r) \, ds.$$ where $W(t)$ is a standard Brownian motion, $\kappa$ and $c$ are both positive.
I know this integral can be viewed as $$u_t = \int_{t-1}^t e^{-\kappa(t-s)}J_c(s) \, ds,$$ where $J_c(s) = \int_0^s e^{-c(s-r)} \, dW(r)$, is an Ornstein-Uhlenbeck process. But how do we handle this double integral and use Ito's isometry to get the variance of it?
Further, does this integral admit a Wold representation, that is, $$u_t = \sum_{j = 0}^\infty F_j \varepsilon_{t-j},$$ where $\varepsilon_t \sim \mathrm{iid}(0, \sigma^2)$.
Here is my calculations, which may be wrong, I am not so sure...
$$\begin{align} Var(u_t) &= E(u_t^2) = E\left(\int_{t-1}^{t} \int_{0}^{s}e^{-\kappa(t-s)-c(s-r)}dW(r) ds\right)^2 \\ &= E\left(\int_{0}^{t} \int_{r}^{t}e^{-\kappa(t-s)-c(s-r)}ds dW(r) \right)^2 \\ &= \int_{0}^{t} \left(\int_{r}^{t}e^{-\kappa(t-s)-c(s-r)}ds\right)^2 dr \\ &= \int_{0}^{t} e^{-2\kappa t + 2cr} \left(\int_{r}^{t}e^{(\kappa-c)s}ds\right)^2 dr \\ &= \dfrac{1}{(\kappa-c)^2}\int_{0}^{t} e^{-2\kappa t + 2cr} \left(e^{(\kappa-c)t}- e^{(\kappa-c)r} \right)^2 dr \\ &= \dfrac{1}{(\kappa-c)^2}\left(\int_{0}^{t} e^{-2c(t-r)} dr + \int_{0}^{t} e^{-2\kappa(t-r)} dr - 2\int_{0}^{t} e^{-(\kappa+c)(t-r)} dr\right)\\ &= \dfrac{1}{(\kappa-c)^2}\left(\dfrac{1-e^{-2ct}}{2c} + \dfrac{1-e^{-2\kappa t}}{2\kappa} - 2\dfrac{1-e^{-(\kappa+c)t}}{\kappa+c}\right) \end{align}$$ which indicates $u_t$ exhibits heteroskedasiticity. But still, I don't know if it is weakly stationary.