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.

Rainwater harvesting (RWH), which is the collection and storage of rainwater for multiple purposes, is gaining recognition in water supply issues. Selection of harvesting sites is the most critical factor in RWH projects. The objective of this study is to develop a suitability map of RWH sites for a basin in Saudi Arabia. The method used, constitute the identification and assigning weights to criteria, and generation of suitability map using Analytical Hierarchy Process (AHP). Eight appropriate criteria were considered. Results showed that excellent and good sites covered about 40.6% of the total available sites. Sensitivity analysis showed that the curve number (CN), slope, rainfall, and soil were the most influential criteria. The maximum increase in the percentage area of excellent sites was 92% while good and moderate classes decreased by 43 and 53%, respectively. The developed suitability maps provide useful information to the decision maker for use in water management.