Soil data and soil maps are crucial for Earth system modeling, water management, agricultural production, and climate change studies, and reducing uncertainties in soil property and soil class maps improves their reliability. Here, we present a pruned Hierarchical Random Forest (pHRF) framework to map soil taxa and properties over the National Ecological Observatory Network (NEON) sites in the Contiguous United States (CONUS). The pHRF method reduces uncertainties in predictions compared to POLARIS v1, providing smaller prediction intervals for the distributions of soil properties. In addition, pHRF addresses two data imbalance issues in soil survey data— uneven spatial distribution of georeferenced soil observations, and secondly underrepresentation of certain soil taxa. Unlike traditional hierarchical soil classification, pHRF conditions the probabilities of finer taxonomic levels based on their parent levels and removes implausible predictions (identified as errors) using field-validated soil taxa, improving prediction intervals. To address the categorical imbalance, soil taxa belonging to minority parent soil taxa are predicted with their own models, without being overlooked compared to using a single model on all soil taxa. For spatial imbalance, each model dynamically adapts its spatial coverage, incorporating more neighboring soil data in areas where georeferenced soil observations are sparse. In data-scarce areas, fieldvalidated soil taxa are resampled to improve the representation of soil variation. The pHRF-derived soil classification showed out-of-bag scores above 0.7 at different taxonomic levels. The probabilistic map of soil series was then used to estimate soil properties, by linking them to a harmonized soil properties database. When evaluated against independent NEON measurements, pHRF performed better than POLARIS v1 for root zone properties (0–60 cm), particularly for sand, clay, and organic matter content. Specifically, pHRF reduced RMSE by 1.15 (sand%), 1.32 (clay%), and 0.21 (log-scaled organic matter%) while improving correlations. For pH, both models showed a reasonable fit (RMSE: ~0.70, correlation: 0.85). This approach presents a development in refining soil properties mapping, especially in its effectiveness in reducing uncertainties. Future work will focus on reducing uncertainties and correcting biases in soil property estimates. Keywords: Soil properties estimation; Pruned Hierarchical Random Forest (pHRF); POLARIS; Digital soil mapping; Soil covariates; Uncertainty reduction
