Understanding　of　heat　and　mass　transfer　in　fractured　porous　media　is　vital　for　the　effective　analysis　and　management　of　subsurface　systems.　We　have　developed　a　transient　analytical　model　that　incorporates　advection,　diffusion,　degradation,　and　adsorption　processes　of　organic　contaminants　in　fractured　porous　media　under　non-isothermal　conditions.　The　solution　is　obtained　through　Laplace　inversion.　The　analytical　model　is　validated　using　experimental　data,　field　test,　and　numerical　solution.　We　investigate　the　influence　of　various　factors,　including　thermal　properties,　flow　conditions,　and　fracture　characteristics　on　organic　contaminants　transport　in　fractured　porous　media.　Our　findings　show　that　when　the　temperature　at　the　inlet　exceeds　the　initial　temperature　of　the　matrix,　contaminant　transport　into　the　matrix　is　enhanced　due　to　the　thermal　gradient.　In　contrast,　contaminant　transport　in　the　fractures　is　comparatively　impeded.　For　example,　the　contaminant　flux　for　the　case　with　a　temperature　difference　(ΔT)　of　−15　K　at　the　fracture　outlet　is　1.8　times　greater　than　that　observed　at　ΔT　=　20　K.　Accurate　prediction　of　contaminant　transport　requires　consideration　of　non-isothermal　diffusion,　especially　in　the　case　involving　elevated　ST　and　|ΔT|.　The　contaminant　flux　at　the　fracture　outlet　for　the　case　with　ΔT　=　20　K　is　68.24　%　of　the　value　without　considering　thermal　diffusion　effects.　Increasing　the　adsorption　capacity　of　the　matrix　(Rm)　reduces　the　diffusive　mass　into　the　matrix　from　the　fracture,　while　an　increase　in　half-life　of　organic　contaminant　(t1/2,m)　significantly　increases　the　mass　flux　at　the　outlet.　Increasing　t1/2,m　from　1　year　to　200　years　results　in　a　rise　of　the　outlet　flux　by　factors　of　82.3　and　70.4　for　the　case　with　Rm　=　10　and　Rm　=　80,　respectively.　When　the　fracture　aperture　reaches　a　threshold　value　(e.g.,　7　mm),　the　mass　transfer　at　the　interface　between　the　fracture　and　matrix　becomes　negligible　compared　to　the　dominant　flow　within　the　fracture.　©　2025