This paper investigates probabilistic failure envelopes of strip foundations on spatially variable soils with profiles of undrained shear strength su linearly increasing with depth using the lower bound random finite element limit analysis. The spatially variable su is characterised by a non-stationary random field with linearly increasing mean and constant coefficient of variation (COV) with depth. The deterministic uniaxial capacities and failure envelopes are firstly derived to validate numerical models and provide a reference for the subsequent probabilistic analysis. Results indicate that the random field parameters COVsu (COV of su) and Δ (dimensionless autocorrelation distance) have a considerable effect on the probabilistic normalised uniaxial capacities which alters the size of probabilistic failure envelopes. However, COVsu and Δ have an insignificant effect on the shape of probabilistic failure envelopes is observed in the V-H, V-M and H-M loading spaces, such that failure envelopes for different soil variabilities can be simply scaled by the uniaxial capacities. In contrast to COVsu and Δ, the soil strength heterogeneity index κ = μkB/μsu0 has the lowest effect on the probabilistic normalised uniaxial capacity factors but the highest effect on the shape of the probabilistic failure envelopes. A series of expressions are proposed to describe the shape of deterministic and probabilistic failure envelopes for strip foundations under combined vertical, horizontal and moment (V-H-M) loading.