The paper examines three-dimensional (3D) analyses of the load bearing capacity of square and strip footings on a spatially variable cohesive-frictional (c-u) soil. Embedded, rigid and smooth footings are considered. Cohesion and friction angle are modelled using random fields. The vertical scale of fluctuation (SOF) is based on earlier studies of natural soils. For the square footing problem, the horizontal SOF is the subject of a parametric study. For the strip footing, only one value of the horizontal SOF is assumed, and the parametric study concerns the modelled length of the foundation. To generate the fields of strength parameters, the Fourier series method is employed. For both problems analyzed, the influence of the parameters on the probabilistic characteristics of the bearing capacity of the soil is investigated. For the square footing problem, both the mean value of bearing capacity and its coefficient of variation increase with an increasing horizontal fluctuation scale, but they stabilize at higher values of SOF. For the strip footing problem, it appears that 2D modelling is a very conservative approach, and for the precise assessment of the probability of foundation failure, the footing should be modelled in 3D. Keywords FLAC3D Á Random field Á Reliability index Á Scale of fluctuation List of symbols Roman a Lower bound of bounded interval b Upper bound of bounded interval CoV Coefficient of variation c Cohesion