Precise assessment of drought and its impact on the natural ecosystem is an arduous task in regions with limited climatic observations due to sparsely distributed in situ stations, especially in the hyper-arid region of Kingdom of Saudi Arabia (KSA). Therefore, this study investigates the application of remote sensing techniques to monitor drought and compare the remote sensing-retrieved drought indices (RSDIs) with the standardized meteorological drought index (Standardized Precipitation Evapotranspiration Index, SPEI) during 2001–2020. The computed RSDIs include Vegetation Condition Index (VCI), Temperature Condition Index (TCI), and Vegetation Health Index (VHI), which are derived using multi-temporal Landsat 7 ETM+, Landsat 8 OLI/TIRS satellites, and the Google Earth Engine (GEE) platform. Pearson correlation coefficient (CC) is used to find the extent of agreement between the SPEI and RSDIs. The comparison showed CC values of 0.74, 0.67, 0.57, and 0.47 observed for VHI/SPEI-12, VHI/SPEI-6, VHI/SPEI-3, and VHI/SPEI-1, respectively. Comparatively low agreement was observed between TCI and SPEI with CC values of 0.60, 0.61, 0.42, and 0.37 observed for TCI/SPEI-12, TCI/SPEI-6, TCI/SPEI-3, and TCI/SPEI-1. A lower correlation with CC values of 0.53, 0.45, 0.33 and 0.24 was observed for VCI/SPEI-12, VCI/SPEI-6, VCI/SPEI-3, and VCI/SPEI-1, respectively. Overall, the results suggest that VHI and SPEI are better correlated drought indices and are suitable for drought monitoring in the data-scarce hyper-arid regions. This research will help to improve our understanding of the relationships between meteorological and remote sensing drought indices.

Satellite precipitation products (SPPs) are undeniably subject to uncertainty due to retrieval algorithms and sampling issues. Many research efforts have concentrated on merging SPPs to create high-quality merged precipitation datasets (MPDs) in order to reduce these uncertainties. This study investigates the efficacy of dynamically weighted MPDs in contrast to those using static weights. The analysis focuses on comparing MPDs generated using the “dynamic clustered Bayesian averaging (DCBA)” approach with those utilizing the “regional principal component analysis (RPCA)” under fixed-weight conditions. These MPDs were merged from SPPs and reanalysis precipitation data, including TRMM (Tropical Rainfall Measurement Mission) Multi-satellite Precipitation Analysis (TMPA) 3B42V7, PERSIANN-CDR, CMORPH, and the ERA-Interim reanalysis precipitation data. The performance of these datasets was evaluated in Pakistan’s diverse climatic zones—glacial, humid, arid, and hyper-arid—employing data from 102 rain gauge stations. The effectiveness of the DCBA model was quantified using Theil’s U statistic, demonstrating its superiority over the RPCA model and other individual merging methods in the study area The comparative performances of DCBA and RPCA in these regions, as measured by Theil’s U, are 0.49 to 0.53, 0.38 to 0.45, 0.37 to 0.42, and 0.36 to 0.43 in glacial, humid, arid, and hyper-arid zones, respectively. The evaluation of DCBA and RPCA compared with SPPs at different elevations showed poorer performance at high altitudes (>4000 m). The comparison of MPDs with the best performance of SPP (i.e., TMPA) showed significant improvement of DCBA even at altitudes above 4000 m. The improvements are reported as 49.83% for mean absolute error (MAE), 42.31% for root-mean-square error (RMSE), 27.94% for correlation coefficient (CC), 40.15% for standard deviation (SD), and 13.21% for Theil’s U. Relatively smaller improvements are observed for RPCA at 13.04%, 1.56%, 10.91%, 1.67%, and 5.66% in the above indices, respectively. Overall, this study demonstrated the superiority of DCBA over RPCA with static weight. Therefore, it is strongly recommended to use dynamic variation of weights in the development of MPDs.