Climate change may cause severe hydrological droughts, leading to water shortages which will require to be assessed using high-resolution data. Gravity Recovery and Climate Experiment (GRACE) satellite Terrestrial Water Storage (TWSA) estimates offer a promising solution to monitor hydrological drought, but its coarse resolution (1°) limits its applications to small regions of the Indus Basin Irrigation System (IBIS). Here we employed machine learning models such as Extreme Gradient Boosting (XGBoost) and Artificial Neural Network (ANN) to downscale GRACE TWSA from 1° to 0.25°. The findings revealed that the XGBoost model outperformed the ANN model with Nash Sutcliff Efficiency (NSE) (0.99), Pearson correlation (R) (0.99), Root Mean Square Error (RMSE) (5.22 mm), and Mean Absolute Error (MAE) (2.75 mm) between the predicted and GRACE-derived TWSA. Further, Water Storage Deficit Index (WSDI) and WSD (Water Storage Deficit) were used to determine the severity and episodes of droughts, respectively. The results of WSDI exhibited a strong agreement when compared with the Standardized Precipitation Evapotranspiration Index (SPEI) at different time scales (1-, 3-, and 6-months) and self-calibrated Palmer Drought Severity Index (sc-PDSI). Moreover, the IBIS had experienced increasing drought episodes, e.g., eight drought episodes were detected within the years 2010 and 2016 with WSDI of −1.20 and −1.28 and total WSD of −496.99 mm and −734.01 mm, respectively. The Partial Least Square Regression (PLSR) model between WSDI and climatic variables indicated that potential evaporation had the largest influence on drought after precipitation. The findings of this study will be helpful for drought-related decision-making in IBIS.

Pakistan’s agriculture and food production account for 27% of its overall gross domestic product (GDP). Despite ongoing advances in technology and crop varieties, an imbalance between water availability and demand, combined with robust shifts in drought propagation has negatively affected the agro-ecosystem and environmental conditions. In this study, we examined hydro-meteorological drought propagation and its associated impacts on crop yield across natural and human-disturbed agro-ecological zones (AEZs) in Pakistan. Multisource datasets (i.e., ground observations, reanalysis, and satellites) were used to characterize the most extensive, intense drought episodes from 1981 to 2018 based on the standardized precipitation evaporation index (SPEI), standardized streamflow index (SSFI), standardized surface water storage index (SSWSI), and standardized groundwater storage index (SGWI). The most common and intense drought episodes characterized by SPEI, SSFI, SSWSI, and SGWI were observed in years 1981–1983, 2000–2003, 2005, and 2018. SPEI yielded the maximum number of drought months (90) followed by SSFI (85), SSWSI (75), and SGWI (35). Droughts were frequently longer and had a slower termination rate in the human-disturbed AEZs (e.g., North Irrigated Plain and South Irrigated Plain) compared to natural zones (e.g., Wet Mountains and Northern Dry Mountains). The historical droughts are likely caused by the anomalous large-scale patterns of geopotential height, near-surface air temperature, total precipitation, and prevailing soil moisture conditions. The negative values (<−2) of standardized drought severity index (DSI) observed during the drought episodes (1988, 2000, and 2002) indicated a decline in vegetation growth and yield of major crops such as sugarcane, maize, wheat, cotton, and rice. A large number of low-yield years (SYRI ≤ −1.5) were recorded for sugarcane and maize (10 years), followed by rice (9 years), wheat (8 years), and cotton (6 years). Maximum crop yield reductions relative to the historic mean (1981–2017) were recorded in 1983 (38% for cotton), 1985 (51% for maize), 1999 (15% for wheat), 2000 (29% for cotton), 2001 (37% for rice), 2002 (21% for rice), and 2004 (32% for maize). The percentage yield losses associated with shifts in SSFI and SSWSI were greater than those in SPEI, likely due to longer drought termination duration and a slower termination rate in the human-disturbed AEZs. The study’s findings will assist policymakers to adopt sustainable agricultural and water management practices, and make climate change adaptation plans to mitigate drought impacts in the study region.

Understanding the pixel-scale hydrology and the spatiotemporal distribution of regional precipitation requires high precision and high-resolution precipitation data. Satellite-based precipitation products have coarse spatial resolutions (~10 km–75 km), rendering them incapable of translating high-resolution precipitation variability induced by dynamic interactions between climatic forcing, ground cover, and altitude variations. This study investigates the performance of a downscaled-calibration procedure to generate fine-scale (1 km × 1 km) gridded precipitation estimates from the coarser resolution of TRMM data (~25 km) in the Indus Basin. The mixed geographically weighted regression (MGWR) and random forest (RF) models were utilized to spatially downscale the TRMM precipitation data using high-resolution (1 km × 1 km) explanatory variables. Downscaled precipitation estimates were combined with APHRODITE rain gauge-based data using the calibration procedure (geographical ratio analysis (GRA)). Results indicated that the MGWR model performed better on fit and accuracy than the RF model to predict the precipitation. Annual TRMM estimates after downscaling and calibration not only translate the spatial heterogeneity of precipitation but also improved the agreement with rain gauge observations with a reduction in RMSE and bias of ~88 mm/year and 27%, respectively. Significant improvement was also observed in monthly (and daily) precipitation estimates with a higher reduction in RMSE and bias of ~30 mm mm/month (0.92 mm/day) and 10.57% (3.93%), respectively, after downscaling and calibration procedures. In general, the higher reduction in bias values after downscaling and calibration procedures was noted across the downstream low elevation zones (e.g., zone 1 correspond to elevation changes from 0 to 500 m). The low performance of precipitation products across the elevation zone 3 (>1000 m) might be associated with the fact that satellite observations at high-altitude regions with glacier coverage are most likely subjected to higher uncertainties. The high-resolution grided precipitation data generated by the MGWR-based proposed framework can facilitate the characterization of distributed hydrology in the Indus Basin. The method may have strong adoptability in the other catchments of the world, with varying climates and topography conditions.

Here, we provided a comprehensive analysis of long-term drought and climate extreme patterns in the agro ecological zones (AEZs) of Pakistan during 1980–2019. Drought trends were investigated using the standardized precipitation evapotranspiration index (SPEI) at various timescales (SPEI-1, SPEI-3, SPEI-6, and SPEI-12). The results showed that droughts (seasonal and annual) were more persistent and severe in the southern, southwestern, southeastern, and central parts of the region. Drought exacerbated with slopes of −0.02, −0.07, −0.08, −0.01, and −0.02 per year. Drought prevailed in all AEZs in the spring season. The majority of AEZs in Pakistan’s southern, middle, and southwestern regions had experienced substantial warming. The mean annual temperature minimum (Tmin) increased faster than the mean annual temperature maximum (Tmax) in all zones. Precipitation decreased in the southern, northern, central, and southwestern parts of the region. Principal component analysis (PCA) revealed a robust increase in temperature extremes with a variance of 76% and a decrease in precipitation extremes with a variance of 91% in the region. Temperature and precipitation extremes indices had a strong Pearson correlation with drought events. Higher temperatures resulted in extreme drought (dry conditions), while higher precipitation levels resulted in wetting conditions (no drought) in different AEZs. In most AEZs, drought occurrences were more responsive to precipitation. The current findings are helpful for climate mitigation strategies and specific zonal efforts are needed to alleviate the environmental and societal impacts of drought.

Accurate global monitoring of carbon dioxide (CO₂) is essential for understanding climate change and informing policy decisions. This study compares column-averaged dry-air mole fractions of CO₂ (XCO₂) between ACOS_L2_Lite_FP V9r for Japan’s Greenhouse Gases Observing Satellite (GOSAT), OCO-2_L2_Lite_FP V10r for the USA’s Orbiting Carbon Observatory-2 (OCO-2), and IAPCAS V2.0 for China’s Carbon Dioxide Observation Satellite (TANSAT) collectively referred to as GOT, with data from the Total Carbon Column Observing Network (TCCON). Our findings are as follows: (1) Significant data quantity differences exist between OCO-2 and the other satellites, with OCO-2 boasting a data volume 100 times greater. GOT shows the highest data volume between 30–45°N and 20–30°S, but data availability is notably lower near the equator. (2) XCO₂ from GOT exhibits similar seasonal variations, with lower concentrations during June, July, and August (JJA) (402.72–403.74 ppm) and higher concentrations during December, January, and February (DJF) (405.74–407.14 ppm). XCO₂ levels are higher in the Northern Hemisphere during March, April, and May (MAM) and DJF, while slightly lower during JJA and September, October, and November (SON). (3) The differences in XCO₂ (ΔXCO₂) reveal that ΔXCO₂ between OCO-2 and TANSAT are minor (−0.47 ± 0.28 ppm), whereas the most significant difference is observed between GOSAT and TANSAT (−1.13 ± 0.15 ppm). Minimal differences are seen in SON (with the biggest difference between GOSAT and TANSAT: −0.84 ± 0.12 ppm), while notable differences occur in DJF (with the biggest difference between GOSAT and TANSAT: −1.43 ± 0.17 ppm). Regarding latitudinal variations, distinctions between OCO-2 and TANSAT are most pronounced in JJA and SON. (4) Compared to TCCON, XCO₂ from GOT exhibits relatively high determination coefficients (R² > 0.8), with GOSAT having the highest root mean square error (RMSE = 1.226 ppm, <1.5 ppm), indicating a strong relationship between ground-based observed and retrieved values. This research contributes significantly to our understanding of the spatial characteristics of global XCO₂. Furthermore, it offers insights that can inform the analysis of differences in the inversion of carbon sources and sinks within assimilation systems when incorporating XCO₂ data from satellite observations.