The major drivers of soil variation in Saskatchewan at scales finer than the existing soil maps are parent material variance, slope position, and salinity. There is therefore a need to generate finer-scale parent material maps as part of updating soil maps in Saskatchewan. As spatially referenced soil point data are lacking in Saskatchewan, predictive soil mapping methods that disaggregate existing soil parent material maps are required. This study focused on investigating important environmental covariates to use in parent material disaggregation, particularly bare soil composite imagery (BSCI). Synthetic point observations were generated using an area-proportional approach based on existing soil survey polygons and a random forest model was trained with those synthetic observations to predict parent material classes. Including BSCI as environmental covariates increased model accuracy from 0.38 to 0.52 and the model Kappa score from 0.19 to 0.35 compared with models where it was not included. Models that included training points from all locations, regardless of whether BSCI was available, and included BSCI as environmental covariates had similar results to the BSCI model with an accuracy of 0.48 and a Kappa value of 0.30. Based on these results, BSCI is an important covariate for parent material disaggregation in the Saskatchewan Prairies. Future work to disaggregate soil classes based on slope position and salinity, and to combine those methods with parent material disaggregation is needed to generate detailed soil maps for the Canadian Prairies.