Transport　processes　in　the　subsurface　are　strongly　influenced　by　the　heterogeneity　of　soil.　Understanding　the　non-Fickian　behavior　of　solute　transport　in　such　environments　is　critically　important　for　the　accurate　prediction　of　the　spatial　and　temporal　distribution　of　mass.　Traditional　local　continuum　models　(e.g.,　classical　differential　equation-based　approaches)　encounter　inherent　limitations　when　dealing　with　strong　heterogeneities,　fractures,　and　other　discontinuous　structures　due　to　their　reliance　on　spatial　derivative　continuity　assumptions.　This　study　develops　a　peridynamic　non-Fickian　transport　model　to　consider　the　non-Fickian　and　non-local　effects　simultaneously.　The　model　is　formulated　in　the　framework　of　bond-based　peridynamic　(PD)　by　introducing　the　concept　of　dual　phase　lag　(DPL)　effects,　including　flux　lag　τJ　quantifying　inertial　collision　delays　from　pore　structure　irregularities　and　storage　lag　τC　reflecting　interfacial　adsorption–desorption　effects,　to　provide　an　alternative　numerical　framework　for　addressing　these　complexities.　The　developed　dual-phase　lag　peridynamic　(DPL-PD)　framework　is　validated　against　the　analytical　solution.　The　effects　of　phase　lag　and　heterogeneity　were　discussed.　A　series　of　case　studies,　including　solute　transport　in　layered　soil,　heterogeneous　soil,　and　porous　media　with　discontinuity,　demonstrate　the　capacity　of　DPL-PD　model　to　capture　the　non-Fickian　and　non-local　transport　of　solute　in　heterogeneous　and　discontinuous　porous　media.　Furthermore,　these　cases　quantitatively　assess　the　influence　of　permeable　and　impermeable　inclusions,　as　well　as　fractures,　on　solute　transport.　The　results　show　that　the　developed　dual-phase　lag　peridynamic　framework　can　offer　an　alternative　to　classical　local　formulations　for　investigating　non-Fickian　and　non-local　transport　in　media　with　physically　realistic　microstructures.　©　2025