Landslides, prevalent in mountainous areas, are typically triggered by tectonic movements, climatic changes, and human activities. They pose catastrophic risks, especially when occurring near settlements and infrastructure. Therefore, detecting, monitoring, and predicting landslide deformations is essential for geo-risk mitigation. The mainstream of the previous studies have often focused on deterministic models for immediate landslide prediction. However, in most of them, the aspect of prediction uncertainties are not sufficiently addressed. This paper introduces an innovative probabilistic prediction method using a Variational Autoencoder (VAE) combined with Gated Recurrent Unit (GRU) to forecast landslide deformations from a generative standpoint. Our approach consists of two main elements: firstly, training the VAE-GRU model to maximize the variational lower bound on the likelihood of historical precipitation data; secondly, using the learned approximated posterior distribution to predict imminent deformations from a generative angle. To assess the prediction quality, we use four widely-used metrics: Prediction Interval Coverage Probability (PICP), Prediction Interval Normalized Average Width (PINAW), Coverage Width-Based Criterion (CWC), and Prediction Interval Normalized Root Mean Square Width (PINRW). The results demonstrate that our proposed VAE-GRU framework surpasses traditional state-of-the-art (SOTA) probabilistic deformation prediction algorithms in terms of accuracy and reliability.