“…We adopted Monte Carlo sampling as the uncertainty propagation method. According to this method, a large number of sample points are randomly selected from the distribution of uncertain parameters and fed into a primary model to propagate uncertainty and, hence, quantify output variability, which is typically represented as a probability density function [ [64] , [65] , [66] ]. In the present study, we employed the generalized agglomeration model as the primary model, with being the output variable, and considered normal distributions in the uncertain parameters.…”