Numerical weather- and climate prediction models rely on soil data to accurately model land surface processes. However, as soil data are produced using soil profiles and maps with multiple sources of uncertainty, wide discrepancies prevail in global soil datasets. Comparison of four commonly used soil datasets in Earth system climate models, i.e., Food and Agriculture Organization soil data, Harmonized World Soil Database, Global Soil Dataset for Earth System Model, and global gridded soil information system SoilGrids, yields widespread differences in southern Africa. This study investigates the simulated land-atmosphere interactions in southern Africa in the context of the uncertainties from applying different global soil datasets. We conducted ensemble simulations using the fully coupled Weather Research and Forecasting Hydrological Modeling system (WRF-Hydro) incorporated with each of the global soil datasets mentioned above. Model simulations were performed at 4-km convection-permitting scale from January 2015 to June 2016. By quantifying model\u27s internal variability and comparing the modeling results, results show that the simulated temperature, soil moisture, and surface energy fluxes are largely impacted by soil texture differences. For instance, changes in soil texture and associated hydrophysical parameters result in large differences in air temperature up to 1.7°C and surface heat flux up to 25 W/m2, and disparities in averaged surface soil moisture differ up to 0.1 m3/m3 in austral summer months. Differences in soil texture characteristics also regulate local climatic conditions differently in the wet and dry seasons as well as in different climatic regions. Furthermore, the thermodynamic differences in surface energy fluxes caused by soil texture demonstrate physical feedback perspective on atmospheric processes, resulting in distinct changes in planetary boundary layer height. This study demonstrates the non-negligible impact of soil data on land surface-atmosphere coupled modeling and highlights the need for consistent consideration of modeling uncertainties from soil data in modeling applications
