Identifying the changes in precipitation and temperature at a regional scale is of great importance for the quantification of climate change. This research investigates the changes in precipitation and surface air temperature indices in the seven irrigation zones of Punjab Province during the last 50 years; this province is a very important region in Pakistan in terms of agriculture and irrigated farming. The reliability of the data was examined using double mass curve and autocorrelation analysis. The magnitude and significance of the precipitation and temperature were visualized by various statistical methods. The stations’ trends were spatially distributed to better understand climatic variability across the elevation gradient of the study region. The results showed a significant warming trend in annual T_min (minimum temperature) and T_mean (mean temperature) in different irrigation zones. However, T_max (maximum temperature) had insignificant variations except in the high elevation Thal zone. Moreover, the rate of T_min increased faster than that of T_max, resulting in a reduction in the diurnal temperature range (DTR). On a seasonal scale, warming was more pronounced during spring, followed by that in winter and autumn. However, the summer season exhibited insignificant negative trends in most of the zones and gauges, except in the higher-altitude Thal zone. Overall, Bahawalpur and Faisalabad are the zones most vulnerable to warming annually and in the spring, respectively. Furthermore, the elevation-dependent trend (EDT) indicated larger increments in T_max for higher-elevation (above 500 m a.s.l.) stations, compared to the lower-elevation ones, on both annual and seasonal scales. In contrast, the T_min showed opposite trends at higher- and lower-elevation stations, while a moderate increase was witnessed in T_mean trends from lower to higher altitude over the study region. An increasing trend in DTR was observed at higher elevation, while a decreasing trend was noticed at the lower-elevation stations. The analysis of precipitation data indicated wide variability over the entire region during the study period. Most previous studies reported no change or a decreasing trend in precipitation in this region. Conversely, our findings indicated the cumulative increase in annual and autumn precipitation amounts at zonal and regional level. However, EDT analysis identified the decrease in precipitation amounts at higher elevation (above 1000 m a.s.l.) and increase at the lower-elevation stations. Overall, our findings revealed unprecedented evidence of regional climate change from the perspectives of seasonal warming and variations in precipitation and temperature extremes (T_max and T_min) particularly at higher-elevation sites, resulting in a variability of the DTR, which could have a significant influence on water resources and on the phenology of vegetation and crops at zonal and station level in Punjab.

The performance evaluation of satellite-based precipitation products at local and regional scales is crucial for modification in satellite-based precipitation retrieval algorithms, as well as for the provision of guidance during the selection of substitute precipitation data. This study evaluated the performances of three Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis (TMPA) products (3B42V6, 3B42RT and 3B42V7) with a reference to rain gauge observations in the Hunza River basin, northern Pakistan. Multi-spatial (pixel and basin) and temporal (daily, monthly, seasonal and annual) resolutions were considered for performance evaluation of TMPA products. Results revealed that the spatial pattern of observed precipitation over the basin was adequately captured by 3B42V7 but misplaced by 3B42V6 and 3B42RT. All TMPA products were unable to capture the intense precipitation events. On the daily time scale, the performance of TMPA products was very poor over both spatial scales. 3B42V6 underestimated the precipitation (31.25% and 44.27% on pixel and basin scales, respectively). By contrast, 3B42RT significantly overestimated the precipitation (47.91% and 38.62% on pixel and basin scales, respectively), while 3B42V7 showed overestimation (17.30%) on pixel scale and slight underestimation (6.24%) on the basin scale. On the seasonal scale, TMPA products showed significant biases with observed precipitation data. We found that the TMPA products performed relatively better on monthly and annual time scales and overall performance of 3B42V7 product was better than that of 3B42V6 and 3B42RT. The bias in 3B42V7 was improved by 85.90% compared with 3B42V6 and by 116.16% compared with 3B42RT. Thus, it is concluded that the TMPA products were unreliable to capture the intense precipitation events and retain high errors on daily and seasonal scales. Therefore, caution should be considered while using these precipitation estimates as a substitute data in hydrology, meteorology and climatology studies in Hunza River basin. However, due to the reasonable performance of monthly and annual 3B42V7 estimates, these can be used as an acceptable substitute for applications in the region.

In recent years, heatwaves have been reported frequently by literature and the media on the Tibetan Plateau. However, it is unclear how alpine vegetation responds to the heatwaves on the Tibetan Plateau. This study aimed to identify the heatwaves using long-term meteorological data and examine the impact of heatwaves on vegetation growth rate with remote sensing data. The results indicated that heatwaves frequently occur in June, July, and August on the Tibetan Plateau. The average frequency of heatwaves had no statistically significant trends from 2000 to 2020 for the entire Tibetan Plateau. On a monthly scale, the average frequency of heatwaves increased significantly (p < 0.1) in August, while no significant trends were in June and July. The intensity of heatwaves indicated a negative correlation with the vegetation growth rate anomaly (ΔVGR) calculated from the normalized difference vegetation index (NDVI) (r = −0.74, p < 0.05) and the enhanced vegetation index (EVI) (r = −0.61, p < 0.1) on the Tibetan Plateau, respectively. Both NDVI and EVI consistently demonstrate that the heatwaves slow the vegetation growth rate. This study outlines the importance of heatwaves to vegetation growth to enrich our understanding of alpine vegetation response to increasing extreme weather events under the background of climate change.

Spatial and temporal precipitation data acquisition is highly important for hydro-meteorological applications. Gridded precipitation products (GPPs) offer an opportunity to estimate precipitation at different time and resolution. Though, the products have numerous discrepancies that need to be evaluated against in-situ records. The present study is the first of its kind to highlight the performance evaluation of gauge based (GB) and satellite based (SB) GPPs at annual, winter, and summer monsoon scale by using multiple statistical approach during the period of 1979–2017 and 2003–2017, respectively. The result revealed that the temporal magnitude of all the GPPs was different and deviate up to 100–200 mm with overall spatial pattern of underestimation (GB product) and overestimation (SB product) from north to south gradient. The degree of accuracy of GB products with observed precipitation decreases with the increase in the magnitude of precipitation and vice versa for SB precipitation products. Furthermore, the observed precipitation revealed the positive trend with multiple turning points during the period 1979–2005. However, the gentle increase with no obvious break point has been detected during the period of 2005–2017. The large inter-annual variability and trends slope of the reference data series were well captured by Global Precipitation Climatology Centre (GPCC) and Tropical Rainfall Measuring Mission (TRMM) products and outperformed the relative GPPs in terms of higher R² values of ≥ 0.90 and lower values of estimated RME ≤ 25% at annual and summer monsoon season. However, Climate Research Unit (CRU) performed better during winter estimates as compared with in-situ records. In view of significant error and discrepancies, regional correction factors for each GPPs were introduced that can be useful for future concerned projects over the study region. The study highlights the importance of evaluation by the careful selection of potential GPPs for the future hydro-climate studies over the similar regions like Punjab Province.