The hydro-meteorological variables of extreme rainfall are not easy to explain due to unexpected changes in climate and varied usage of water with growing population. Regional rainfall frequency analysis is the one such method that is useful for the requirement of more accurate estimates of rainfall yearly or desineally for the regions having lack of fresh water resources. The series of Annual Maximum Monthly Rainfall Totals (AMMRT) has been used for the seven sites of northern Punjab, Pakistan using L-moments. The results of different test, the run test, lag-1 correlation and Mann-Whitney U test illustrate that the data series of the seven sites of northern Punjab were found random and independently and identically distributed and have no serial correlation. Heterogeneity measure exposed that the region is homogeneous and discordancy measure gives the evidence that no site is discordant among the seven. The result of goodness of fit test including L-moment Ratio diagrams, ZDIST statistic and Mean Absolute Deviation Index exposed the Pearson Type III (PE3), Generalized Normal (GNO) and Generalized Extreme Value(GEV) are best suitable of the regional distribution for the quantiles estimation. The quantiles estimates obtained for different return periods. A linear regression model was developed with good fit between the at site characteristics and the mean of the AMMRT of the sites. The estimates of the study may be used for the estimation of the rainfall quantiles of the seven sites for different return periods. The estimates will be useful to design future preventive measures for the harmful impact of hydro meteorological events at these sites in Punjab Pakistan.