In the present study, we estimated the aerosol radiative forcing and heating rates near Yala Glacier, Nepal (28.21°N, 85.61°E; 4900 masl), using in situ black carbon (BC) mass concentration measurements, satellite data sets, and model simulations. The real-time ambient BC mass concentration was continuously measured using an Aethalometer (AE-33) from October 2016 to May 2017. The Optical Properties of Aerosols and Clouds (OPAC) model was used to simulate the aerosol optical properties in conjunction with the in situ measurements and satellite data sets. Outputs from OPAC and the satellite data sets were used as inputs for the Santa Barbara Discrete Ordinate Radiative Transfer Atmospheric Radiative Transfer (SBDART) model to estimate the radiative forcing. The in situ measurements showed that the BC mass concentration peaked during the pre-monsoon season (707.9 ± 541.8 ng m–3), which was corroborated by the higher aerosol optical depth (AOD) values during this season (0.058 ± 0.002). The diurnal cycle of the BC mass concentration exhibited a night-time low and afternoon high, which were influenced by the boundary layer dynamics and valley wind flow pattern. The Concentration Weighted Trajectory (CWT) analysis indicated diverse source regions, including northern Asia, the Indo-Gangetic Plain (IGP), and parts of Nepal and Bangladesh. The Moderate Resolution Imaging Spectroradiometer (MODIS)-derived AOD and Ångström exponent (AE), and the OPAC-simulated single-scattering albedo (SSA) and asymmetry parameter (AP) over the study site were estimated to be 0.048 ± 0.009 and 1.32 ± 0.01, and 0.938 ± 0.019 and 0.710 ± 0.042, respectively, during the study period. The mean radiative forcing during the study period for the top of the atmosphere, surface and atmosphere were 3.4, –0.5 and 3.9 W m–2, respectively. Higher atmospheric forcing was observed in the pre-monsoon season, leading to changes in the heating rates.