The mountainous region of Asir is experiencing rapid and unsystematic urbanization leading to an increase in land surface temperatures (LST), which poses a challenge to human well-being and ecological balance. Therefore, it is necessary to study the interaction between land use and land cover (LULC)-induced urbanization and LST using advanced geostatistical techniques. In addition, understanding the urbanization process and urban density is essential for effective urban planning and management. The aim of this study was to investigate the interaction between the urbanization process, urban density and the associated LST. Using the Random Forest Algorithm, LULC mapping was conducted for the years 1990, 2000 and 2020. Metrics such as land cover change rate (LCCR), land cover index (LCI), landscape expansion index (LEI), mean landscape expansion index (MLEI) and area-weighted landscape expansion index (AWLEI) were used to understand urbanization processes and LULC changes. Convolutional kernels were used to model urban density, and the mono-window algorithm was applied to analyse LST in the selected years. In addition, the study assessed the Surface Urban Heat Island (SUHI) contribution index to LULC and used Generalized Additive Models (GAMs) in conjunction with Partial Dependence Plots (PDPs) to understand the relationship between urbanization processes, urban density and LST. In a detailed 30-year study, the application of the RF algorithm showed significant shifts in LULC with an overall validation accuracy of over 85%. Urban areas grew dramatically from 69.40 km2 in 1990 to 338.74 km2 in 2020, while water areas decreased from 1.51 to 0.54 km2. Dense vegetation increased from 43.36 to 52.22 km2, indicating positive ecological trends. The LST analysis showed a general warming, with the mean LST increasing from 40.51°C in 1990 to 46.73°C in 2020 and the highest temperature category (50-60°C) increasing from 0.78 to 33.35 km2. The built-up area of cities tripled between 1990 and 2020, with the Landscape Expansion Index reflecting significant growth in suburban areas. The modeling of urban density shows increasing urbanization in the centre, which will expand significantly to the east by 2020. The contribution of LULC to LST and the Urban Heat Island (SUHI) effect was evident, with built-up areas showing a constant temperature increase. GAMs confirmed a statistically significant relationship between urban density and LST, with different effects for different types of urban expansion. This comprehensive study quantitatively sheds light on the complicated dynamics of urbanization, land cover change and temperature variation and provides important insights for sustainable urban development.