Abstract: Due to the geological tectonic movements such as stratum uplift, fractures and folds, natural soils often exhibit rotated anisotropy. In addition, a large number of field trials show that soil parameters will show obvious non-stationarity with the increase of buried depth. It is necessary to study the non-stationary random field simulation method for rotated anisotropy of soil parameters. On the basis of the non-stationary random field simulation method proposed by Griffiths, through the coordinate transformation of the autocorrelation function of the soil parameters, the non-stationary random field and the rotated anisotropy correlation structure are merged, and the rotated anisotropy non-stationary random field simulation method for soil parameters is proposed. A saturated clay slope is taken as an example to study the influences of the spatial variability of undrained shear strength on the slope stability. The results show that the random field which ignores the non-stationary characteristics of soil parameters will underestimate the reliability of the slope. The rotated anisotropy correlation structure has a significant impact on the reliability of the slope, and this effect gradually weakens with the increase of the parameter variability.