Soils with spatially varying shear strengths are modeled using random field theory and elasto-plastic finite element analysis to evaluate the extent to which spatial variability and cross-correlation in soil properties (c and ϕ) affect bearing capacity. The analysis is two dimensional, corresponding to a strip footing with infinite correlation length in the out-of-plane direction, and the soil is assumed to be weightless with footing placed on the soil surface. Theoretical predictions of the mean and standard deviation of bearing capacity, for the case where c and ϕ are independent, are derived using a geometric averaging model and then verified via Monte Carlo simulation. The standard deviation prediction is found to be quite accurate, while the mean prediction is found to require some additional semi-empirical adjustment to give accurate results for "worst case" correlation lengths. Combined, the theory can be used to estimate the probability of bearing-capacity failure, but also sheds light on the stochastic behaviour of foundation bearing failure.Key words: bearing capacity, probability, random fields, geometric averaging, cϕ soil, Monte Carlo simulation.