Variability of seasonal thaw depth in permafrost regions: a stochastic modeling approach