This study elucidates the behavior of Markov-Chains Monte Carlo ensemble samplers for vadose zone inverse modeling by performing an in-depth comparison of four algorithms that use Affine-Invariant (AI) moves or Differential Evolution (DE) strategies to approximate the target density. Two Rosenbrock toy distributions, and one synthetic and one actual case study focusing on the inverse estimation of soil hydraulic parameters using HYDRUS-1D, are used to compare samplers in different dimensions d. The analysis reveals that an ensemble with N=d+1 chains evolved using DE-based strategies converges to the wrong stationary posterior, while AI does not suffer from this issue but exhibits delayed convergence. DE-based samplers regain their ergodic properties when using N≥2d chains. Increasing the number of chains above this threshold has only minor effects on the samplers’ performance, while initializing the ensemble in a high-likelihood region facilitates its convergence. AI strategies exhibit shorter autocorrelation times in the 7d synthetic vadose zone scenario, while DE-based samplers outperform them when the number of soil parameters increases to 16 in the actual scenario. All evaluation metrics degrade as d increases, thus suggesting that sampling strategies based only on interpolation between chains tend to become inefficient when the bulk of the posterior lays in increasingly small portions of the parameters’ space.