The current study focuses on a probabilistic three-dimensional finite element analysis of geosynthetic-reinforced pile-supported embankments. The probabilistic analysis utilizes a finite element technique and was executed using Python scripting in Abaqus. A deterministic model is established for both unreinforced and geosynthetic-reinforced pile-supported embankments, employing a basic linear elastic perfectly plastic Mohr-Coulomb constitutive model. This numerical model adopted a simplified methodology that eliminates the need to simulate actual piles and subsoil. The results are validated across different embankment heights and mesh sizes. Based on results from Monte-Carlo simulations, we aimed to show how results can vary even by considering only one random variable. Probability density functions for various parameters, including pile efficacy, maximum differential settlement, and maximum strain in the reinforcement, are shown. The findings illustrate the ability of probabilistic approaches to predict the most likely occurrences for different governing parameters of geosynthetic-reinforced pile-supported embankments. They also highlight the potential of probabilistic approaches to expand further into a full-scale comprehensive study which can offer deeper insights into the stability of these embankments considering various sources of uncertainties.