Accurate and continuous estimation of surface albedo is vital for assessing and understanding land–surface–atmosphere interactions. We developed a method for estimating instantaneous all-sky at-surface shortwave upwelling radiance and albedo over the Tibetan Plateau. The method accounts for the complex interplay of topography and atmospheric interactions and aims to mitigate the occurrence of data gaps. Employing an RTLSR-kernel-driven model, we retrieved surface shortwave albedo with a 1 km resolution, incorporating direct, isotropic diffuse; circumsolar diffuse; and surrounding terrain irradiance into the all-sky solar surface irradiance. The at-surface upwelling radiance and surface shortwave albedo estimates were in satisfactory agreement with ground observations at four stations in the Tibetan Plateau, with RMSE values of 56.5 W/m² and 0.0422, 67.6 W/m² and 0.0545, 98.6 W/m² and 0.0992, and 78.0 98.6 W/m² and 0.639. This comparison indicated an improved accuracy of at-surface upwelling radiance and surface albedo and significantly reduced data gaps. Valid observations increased substantially in comparison to the MCD43A2 data product, with the new method achieving an increase ranging from 40% to 200% at the four stations. Our study demonstrates that by integrating terrain, cloud properties, and radiative transfer modeling, the accuracy and completeness of retrieved surface albedo and radiance in complex terrains can be effectively improved.

Accurate knowledge of the at-surface solar irradiance (SSI) is essential for retrieving surface and atmospheric properties using satellite measurements of backscattered and reflected radiance. The latter is affected by surface-atmosphere interactions, including the effects of terrain. The SSI is affected by the same processes. This study proposes a method to estimate the components of instantaneous SSI: direct, isotropic and circumsolar diffuse, and terrain irradiance, which is expected to improve the simultaneous retrieval of aerosol optical depth (AOD) and surface reflectance. The method takes into account the coupled effects of topography and atmosphere by combining parameterization and the lookup table (LUT) approaches. The method was applied to rugged terrain over the Tibetan Plateau using Moderate Resolution Imaging Spectrometer (MODIS) atmosphere and surface data, the fifth generation European Centre for Medium-Range Weather Forecasts reanalysis (ERA5) data, Cloud-Aerosol Lidar With Orthogonal Polarization (CALIOP) aerosol data, and a digital elevation model (DEM). The results showed that the SSI estimates were in satisfactory agreement with ground observations at four stations over the Tibetan Plateau (TP) in 2018 with R2 values of 0.61, 0.44, 0.41, and 0.49, respectively, and root mean square error (RMSE) of 205.7, 176.9, 186.0, and 201.2 W/m2, respectively. Estimations of the diffuse irradiance were evaluated separately against the only available in situ observations at the Dali Station, and the results were better than our SSI estimates with R2, RMSE, and relative bias (BIAS) being 0.71, 94.98 W/m2, and 31%, respectively. The isotropic and circumsolar diffuse irradiances accounted for 37.57% and 7.68% of the total annual SSI, respectively, while diffuse irradiance accounted for 46.48% of the total annual SSI. Under clear skies, every 0.1 increase in AOD caused about a 35-W/m2 increase in diffuse irradiance and a decrease of about 25 W/m2 of SSI.