As an extension to the well-established stationary elliptic partial differential equation (PDE) with random continuous coefficients we study a time-dependent advection-diffusion problem, where may have spatial discontinuities. In subsurface flow model, randomness in parabolic account for insufficient measurements or uncertain material procurement, while discontinuities could represent transitio...