Understanding rainfall’s statistical distribution is crucial for effective water resource management, disaster mitigation, and climate adaptation in tropical regions. This study identifies the best-fit probability distributions for monthly rainfall in the Lake Toba region, Indonesia, based on long-term data from 34 rain gauge stations. Ten commonly used probability distributions were evaluated, with parameters estimated via Maximum Likelihood Estimation (MLE). The Kolmogorov–Smirnov (KS) test was applied to assess model goodness-of-fit at each station and month. Results indicate that the Generalized Extreme Value (GEV), Gamma, and Weibull distributions consistently provide the best fit for most stations and regencies, while Exponential and Inverse Gaussian distributions perform poorly. Spatial analysis reveals notable variation in best-fit models among regencies, emphasizing the influence of local topography and microclimate. These results highlight the need to select flexible probability models for hydrological planning and climate risk assessment in complex tropical regions. The findings provide valuable references for rainfall modeling and bias correction elsewhere.