“…where f pre,i (x) and f post,i (x) are the probability density functions of ith IHA in pre-and post-impact periods, respectively. Since the IHAs cover a variety of indicators such as flow, duration, number, and rate, Vahedberdi et al (2022) suggested applying Gaussian kernel function to simulate the density function of IHA based on the theory of Gaussian Kernel Density Estimation (KDE), which has well fitted the distribution of the indicators in practice (Wglarczyk, 2018). The Gaussian kernel function is given as follows (Sheather, 2004;Silverman, 1986):…”