Accurate river ice mapping is crucial for predicting and managing floods caused by ice jams and for the safe operation of hydropower and water resource facilities. Although satellite multispectral images are widely used for river ice mapping, atmospheric contamination limits their effectiveness. This study developed river ice mapping models for the Han River in South Korea using atmospherically uncorrected Landsat-8 Operational Land Imager (OLI) multispectral reflectance data, addressing atmospheric influences with a Random Forest (RF) classification approach. The RF-based river ice mapping models were developed by implementing various combinations of input variables, incorporating the Landsat-8 multispectral top-of-atmosphere (TOA) reflectance, normalized difference indices for snow, water, and bare ice, and atmospheric factors such as aerosol optical depth, water vapor content, and ozone concentration from the Moderate Resolution Imaging Spectroradiometer observations, as well as surface elevation from the GLO-30 digital elevation model. The RF model developed using all variables achieved excellent performance in the classification of snow-covered ice, snow-free ice, and water, with an overall accuracy and kappa coefficient exceeding 98.4% and 0.98 for test samples, and higher than 83.7% and 0.75 when compared against reference river ice maps generated by manually interpreting the Landsat-8 images under various atmospheric conditions. The RF-based river ice mapping model for the atmospherically corrected Landsat-8 multispectral surface reflectance was also developed, but it showed very low performance under atmospheric conditions heavily contaminated by aerosol and water vapor. Aerosol optical depth and water vapor content were identified as the most important variables. This study demonstrates that multispectral reflectance data, despite atmospheric contamination, can be effectively used for river ice monitoring by applying machine learning with atmospheric auxiliary data to mitigate atmospheric effects.