Concentration fields of solutes in porous media often exhibit large fluctuations, driven by physical and chemical heterogeneity from the pore to Darcy scale. For many applications, ranging reactive transport modeling toxicology, knowledge mean concentrations is not sufficient, quantifying concentration variability necessary. The probability density function (PDF) quantifies frequency occurrence...