“…This model is created on basis of previous microscopic models which described pore-scale flow patterns and/or conservative solutes transport in biofilms and open pore spaces system (e.g., Deng et al, 2013;Karimifard et al, 2021), and the previous microscale models that coupled the three interconnected processes but typically neglected or simplified the flow within biofilms (e.g., Chowdhury et al, 2020;Knutson et al, 2005;Tang et al, 2013;Thullner & Baveye, 2008). The modeling framework uses a realistic 2-D pore network structures (Figure 1b) that are obtained from the experimental work of Weng (2020) and Zhou et al (2020). The 2-D model domain has a size of 8 cm [length] × 6 cm [width], and such size is similar with the characterizations of the HZ water-sediment interfaces reported in recent modeling and field investigations (e.g., Briggs et al, 2018;Chowdhury et al, 2020).…”