“…The time solutes and particles spend in the immobile zone is controlled by the following parameters: (a) the power-law residence time distribution of solute within the immobile zone, set by the power law slope, β S , (b) the rate of fine particle immobilization with the immobile zone, Λ IP , and (c) the powerlaw residence time distribution of particles in the immobile zone, set by the power-law slope, β IP .Following the fitting procedure outlined inDrummond, Schmadel, Kelleher, Packman, and Ward (2019), we performed several computational experiments with simulations and parameter sets constrained to match the conservative solute and fine particle breakthrough curves. We sampled the parameter space using a Latin Hypercube approach (N = 27,000; e.g.,Kelleher et al, 2019). TheBalanced mean square error (θ ; Bottacin-Busolin, Marion, Musner, Tregnaghi, & Zaramella, 2011) objective function was calculated for each simulation as:θ = 1 n…”