The　wide　production　and　application　of　engineered　nanomaterials　(ENMS)　inevitably　lead　to　their　release　in　the　groundwater　environment.　However,　the　release　and　transport　of　boron　nanoparticles　in　the　multi-dimensional　subsurface　remain　largely　unknown.　In　this　work,　a　multi-dimensional　numerical　simulator　for　the　transport　of　boron　nanoparticles　in　the　saturated　porous　media　was　first　developed　and　validated.　Hypothetical　scenarios　for　the　release　of　boron　nanoparticles　into　a　layered　two-dimensional　(2D)　and　heterogeneous　2D　saturated　porous　media　were　then　explored,　and　compared　with　the　fullerene　nanoparticles.　The　results　demonstrated　that　the　soil　heterogeneity　influenced　the　fate　of　nanoparticles,　with　high　permeable　layers　and　high　aqueous-phase　concentration.　Besides,　the　boron　nanoparticles　tend　to　accumulate　at　the　inlet　zones,　where　it　was　closer　to　a　nanoparticles　source.　Different　layers　of　interface　interaction　also　impact　the　concentration　of　nanoparticles.　In　general,　the　mobility　and　aqueous-phase　concentration　of　fullerene　nanoparticles　were　higher　than　those　of　the　boron　nanoparticles.　In　addition,　the　mobility　of　boron　nanoparticles　was　found　to　be　sensitive　to　release　concentration,　soil　porosity　and　nanoparticle　aggregate　size.