Spatial-temporal rainfall assessments are integral to climate/hydrological modeling, agricultural studies, and water resource planning and management. Herein, we evaluate spatial-temporal rainfall trends and patterns in Pakistan for 1961-2020 using nationwide data from 82 rainfall stations. To assess optimal spatial distribution and rainfall characterization, twenty-seven interpolation techniques from geo-statistical and deterministic categories were systematically compared, revealing that the empirical Bayesian kriging regression prediction (EBKRP) technique was best. Hence, EBKRP was used to produce and utilize the first normal annual rainfall map of Pakistan for evaluating spatial rainfall patterns. While the largest under-prediction was estimated for Hunza (- 51%), the highest and lowest rainfalls were estimated for Malam Jaba in Khyber Pakhtunkhwa province (~ 1700 mm), and Nok-kundi in Balochistan province (~ 50 mm), respectively. A gradual south-to-north increase in rainfall was spatially evident with an areal average of 455 mm, while long-term temporal rainfall evaluation showed a statistically significant (p = 0.05) downward trend for Sindh province. Additionally, downward inter-decadal regime shifts were detected for the Punjab and Sindh provinces (90% confidence). These results are expected to help inform environmental planning in Pakistan; moreover, the rainfall data produced using the optimal method has further implications in several aforementioned fields.