Understanding porous fluid flow is key for many geology applications. Traditional methods cannot resolve cases with sharp discontinuities in hydraulic/mechanical properties across those layers. Here we present a new space-time method that can handle such discontinuities. This approach is coupled with trace element transport. Our study reveals that the layering of rocks significantly influences the formation of fluid-rich channels and the material distribution adjacent to discontinuities.