Utilized SciPy’s stats module to fit six probability distributions to maximum yearly discharge data of four different sites. Checked applicability of Gaussian, Exponential, Gamma, Lognormal distributions and plotted histograms.