Evolutionary integral equations as appearing in the theory of linear parabolic viscoelasticity are studied in the presence of white noise. It is shown that the stochastic convolution leads to regular solutions, and that under suitable assumptions the samples are HH older-continuous. These results are put in a wider perspective by consideration of equations with fractional derivatives which are ...