This dataset is an asset of the scientific manuscript "BIS-4D: Mapping soil properties and their uncertainties at 25m resolution in the Netherlands" (Helfenstein et al., 2024, under review). It contains maps of soil properties and their uncertainties at 25m resolution in the Netherlands obtained using the BIS-4D soil modelling and mapping platform. BIS-4D is based on well-established digital soil mapping practices. This dataset includes maps of predictions of the mean, 0.05, 0.50 (median) and 0.95 quantiles and the 90th prediction interval width (PI90) of clay content [%], silt content [%], sand content [%], bulk density (BD) [g/cm3], soil organic matter (SOM) [%], pH [KCl], total N (Ntot) [mg/kg], oxalate-extractable P (Pox) [mmol/kg] and cation exchange capacity (CEC) [mmol(c)/kg]. Prediction maps are available for the standard depth layers specified by the GlobalSoilMap initiative (0-5, 5-15, 15-30, 30-60, 60-100 and 100-200cm). For SOM, these prediction maps are available for the years 1953, 1960, 1970, ...