Experimental and numerical investigation of microparticle transport and deposition in the context of permafrost thaw

In natural porous environments, soil particle migration during flow plays an important role in soil stability and pollutant transport by affecting soil mechanical properties and water quality. In northern areas, permafrost degradation alters the subsurface connection pathways leading to mass movements and rearrangement of the soil. To date, few models have included the influence of temporal and spatial variations of flow velocity and porous media heterogeneity on the transport and deposition of suspended particles.

In this study, laboratory column experiments and a numerical model were used to investigate these issues. The laboratory column experiments were carried out under different flow rates and the effect of porous media heterogeneity was investigated using different grain size distributions. The soil columns were reconstituted from several samples taken in the studied site, the Tasiapik Valley, located in the discontinuous permafrost zone near Umiujaq, Nunavik, Québec. During the experiments, the spatio-temporal distribution of the porosity and the hydraulic conductivity was monitored using X-ray computed tomography imaging (CT-SCAN). Using the pore water velocity computed from the groundwater flow solution, the advection–dispersion transport equation with a first-order kinetic term for particle deposition was solved using the finite element model Heatflow/Smoker. The dependency of the attachment kinetics on the pore water velocity and on the porous media heterogeneity was included. The model was tested and validated with an analytical solution and calibrated with the experimental data. Our simulations highlight the roles of hydrodynamic conditions and soil characteristics on particle transport and deposition mechanisms and the susceptibility of the porous medium to thermo-suffosion in permafrost environments.