This study explores the spatio-temporal distribution and trends on monthly, seasonal, and annual scales of rainfall in the central Punjab districts of Punjab province in Pakistan by using observation and satellite data products. The daily observed data was acquired from the Pakistan Metrological Department (PMD) between 1983 and 2020, along with one reanalysis, namely the Climate Hazard Infrared Group Precipitation Station (CHIRPS) and one satellite-based daily Precipitation Estimation from Remotely Sensed Information Using Artificial Neural Networks climate data record (PERSIANN-CDR) using the Google Earth Engine (GEE) web-based API platform to investigate the spatio-temporal fluctuations and inter-annual variability of rainfall in the study domain. Several statistical indices were employed to check the data similarity between observed and remotely sensed data products and applied to each district. Moreover, non-parametric techniques, i.e., Mann–Kendall (MK) and Sen’s slope estimator were applied to measure the long-term spatio-temporal trends. Remotely sensed data products reveal 422.50 mm (CHIRPS) and 571.08 mm (PERSIANN-CDR) mean annual rainfall in central Punjab. Maximum mean rainfall was witnessed during the monsoon season (70.5%), followed by pre-monsoon (15.2%) and winter (10.2%). Monthly exploration divulges that maximum mean rainfall was noticed in July (26.5%), and the minimum was in November (0.84%). The district-wise rainfall estimation shows maximum rainfall in Sialkot (931.4 mm) and minimum in Pakpattan (289.2 mm). Phase-wise analysis of annual, seasonal, and monthly trends demonstrated a sharp decreasing trend in Phase-1, averaging 3.4 mm/decade and an increasing tendency in Phase-2, averaging 9.1 mm/decade. Maximum seasonal rainfall decreased in phase-1 and increased Phase-2 during monsoon season, averaging 2.1 and 4.7 mm/decade, whereas monthly investigation showed similar phase-wise tendencies in July (1.1 mm/decade) and August (2.3 mm/decade). In addition, as district-wise analyses of annual, seasonal, and monthly trends in the last four decades reveal, the maximum declined trend was in Sialkot (18.5 mm/decade), whereas other districts witnessed an overall increasing trend throughout the years. Out of them, Gujrat district experienced the maximum increasing trend in annual terns (50.81 mm/decade), and Faisalabad (25.45 mm/decade) witnessed this during the monsoon season. The uneven variability and trends have had a crucial imprint on the local environment, mainly in the primary activities.