Geotechnical engineering problems are characterized by many sources of uncertainty. Some of these sources are connected to the uncertainties of soil properties involved in the analysis. In this paper, a numerical procedure for a probabilistic analysis that considers the spatial variability of soil properties is presented to study the response of spatially random soil. The approach integrates a commercial finite difference method and random field theory into the framework of a probabilistic analysis. Two-dimensional non-Gaussian random fields are generated based on a Karhunen-$Lo{grave{e}}ve$ expansion in a fashion consistent with a specified marginal distribution function and an autocorrelation function. A Monte Carlo simulation is then used to determine the statistical response based on the random fields. A series of analyses were performed to study the effects of uncertainty due to the spatial heterogeneity on the settlement and bearing capacity of a rough strip footing. The simulations provide insight into the application of uncertainty treatment to the geotechnical problem and show the importance of the spatial variability of soil properties with regard to the outcome of a probabilistic assessment.