“…The Monte Carlo method was adopted in this section to carry out the statistical analysis of different parameters (e.g., retardation factor, degradation of contaminants, the porosity of CCL, and crack width) (Jacoboni and Reggiani 1983). The ranges of the half-width of the crack (b), the crack retardation factor (R f ), the half-life of organic contaminants in cracks (t f, 1/2 ), and the porosity of CCL (θ) was assumed to be 0-25 mm (Basnett and Brungard 1992;Xie et al 2019), 1-19 (Tang et al 1981;Yan et al 2021a), 10-90 years (Yan et al 2021a), 0.2-0.44 (Musso et al 2020), respectively. The average values of these parameters b, R f , t f,1/2 , θ are assumed to be 12.5 mm, 10, 50 years, 0.32, respectively, for the reference scenario.…”