I like Richard's answer, but I think we can compute the mean and the variance of $\int_0^T W_t dt$ by ourselves using Ito's lemma. Let $f(W_t, t) = t W_t$.
$$
d( t W_t ) = W_t dt + t dW_t .
$$
Integrating both sides, and re-arranging the terms, we get
$$
\int_0^T W_t dt = T W_T - \int_0^T t dW_t \, .
$$
We'll be using Ito's isometry formula $\mathbb{E} \left[ \int_0^T f_t dW_t \int_0^T g_t dW_t \right] = \int_0^T \mathbb{E} \left[f_t g_t \right] dt$.
The integral $\int_0^T W_t dt$ is a Gaussian random variable with zero mean
$$
\mathbb{E} \left[ \int_0^T W_t dt \right] = T \mathbb{E} \left[ W_T \right] - \mathbb{E} \left[ \int_0^T t dW_t \right] = 0,
$$
and variance
$$
\mathbb{E} \left[ \left(\int_0^T W_t dt \right)^2 \right] = T^2 \mathbb{E} [W_T^2] - 2 T \mathbb{E} \left[W_T \int_0^T t dW_t \right] + \mathbb{E} \left[ \left(\int_0^T t dW_t \right)^2 \right]
$$
$$
= T^3-2 T \int_0^T t dt + \int_0^T t^2 dt = \frac{T^3}{3}.
$$
Hence, continuing with Richard's derivations, $Z(T)$ is a Gaussian random variable with mean $S(0)$ and variance $\sigma^2 T/3$.