Well, your second equation, Z(t) = Z(t-1) + sqrt(a)ε(t-1) + sqrt(b) ε(t-2), looks like something similar to GARCH(1,2) or GARCH(2,1) model. Anyway, eventually, first, you can argue that the difference Z(t) - Z(t-1) is a normal random variable provided ε's are normal. And, second, applaying / using the Cramér's decomposition theorem