This study focuses on soil thickness modeling of the Joshimath-Badrinath road section in the Joshimath region of Chamoli district, Uttarakhand, India, a crucial area for local transportation, pilgrimage, and military logistics. The study area, located between longitudes 79.46°E-79.62°E and latitudes 30.53°N-30.82°N, covers 213.43 km2. The terrain is rugged, with slopes ranging from 20 to over 60°, and features diverse vegetation and moderate climatic conditions. The methodology employed integrates soil thickness data from boreholes, a geological map, and a digital elevation model (DEM). Using the GIST model, which incorporates geomorphological and geological factors, soil thickness was modeled. Key factors include hillslope position, surface curvature, and slope gradient, with adjustments made to fit the study area's specific characteristics. Results reveal that the GIST model predicts soil thickness ranging from 0 to 30 m, with a tendency to underestimate mean and maximum values. The model's performance was evaluated against observed data, showing a moderate correlation (R2 = 0.7172), with a root mean square error (RMSE) of 3.4733 m and a mean absolute error (MAE) of 2.6260 m. The model tends to underestimate soil thickness, particularly in extreme values, as indicated by a positive residual skewness (0.6778). While the GIST model offers a reasonable estimate of soil thickness, it tends to fall short in its predictions and sometimes provides unrealistic values. Recognizing these limitations is crucial for refining the model and improving soil thickness assessments in related geotechnical and environmental scenarios.