Transport of reactive contaminant in a wetland flow with the effects of reversible and irreversible reactions on the bed surface

Wetlands are crucial for maintaining biodiversity in ecosystems. The reversible chemical reactions can significantly affect the pollutant reduction capacity in wetland flow. The spatial transport of contaminants in wetlands is further complicated by phase exchange kinetics due to reversible reactions in the wetland bed. This paper provides an analytical solution for studying the multi-dimensional concentration distribution of solutes in wetland flow, considering both reversible and irreversible boundary reactions. Unlike previous one-dimensional studies (Jiang et al. (2022), Zhan et al. (2024)), our approach offers a comprehensive view of wetland contaminant transport, accounting for complex interactions between different phases and environmental factors. The solute may undergo reversible sorptive phase exchange between the immobile (wetland bed phase) and the mobile phase (fluid phase). Multi-scale homogenization methods are used to determine the dispersion coefficient, mean concentration, longitudinal and vertical real concentration distributions, rate of variation of vertical concentration distribution under the influence of vegetation force, and reversible and irreversible reactions at a stationary wetland bed. The effects of retention parameter (theta), Damkohler number (Da), vegetation effect (alpha), dispersion time (tau), and first-order irreversible boundary reaction (beta') on solute mixing are thoroughly analyzed. The study's analytical results are validated by numerical results obtained using the finite difference method with a Shishkin mesh, and other comparisons are made with the limiting case of Taylor dispersivity. Results show that an increase in Damkohler number and vegetation force weakens contaminant transport in wetland flow, but an increase in D-a enhances contaminant concentration in the mobile phase. The presence of reversible phase exchange kinetics causes non-uniform vertical concentration variation across the wetland cross-section. Critical values where the multi-scale homogenization technique may fail are suggested as theta <= 0.5 and D-a >= 1. The study's findings have potential applications to pollution control in wetlands.