This project concerns the effect of freeze-thaw cycles on the pore-scale structure of nonaqueous phase liquid (NAPL) contaminants in water-saturated porous media. This problem is of critical importance to the entrapment of such contaminants in cold temperate, polar and high altitude regions, and has not been examined in the literature to date. This research work is conducted in three stages: (i) two-dimensional nondestructive visualisation of residual light non-aqueous phase liquid (LNAPL), and dense non-aqueous phase liquid (DNAPL), in porous media subject to successive freeze-thaw cycles; (ii) three-dimensional experiments on LNAPL in porous media subject to freeze-thaw, with quantification of phase volumes by X-ray micro-computed tomography (micro-CT); and (iii) the explanation of results by several pore scale mathematical and conceptual models. The two-dimensional cell experiments (using a monolayer of 0.5 mm diameter glass beads held between two glass sheets), and three-dimensional X-ray micro-CT experiments reveal a substantial mobilisation and rupture of ganglia during successive freeze-thaw cycles; this includes the detachment of smaller ganglia from larger ganglia and the mobilisation of NAPL in the direction of freezing front. The experiments also reveal significant shedding of numerous single/sub-singlet ganglia along narrow pore corridors, their entrapment in growing polycrystalline ice, and the coalescence of such small ganglia during thawing to form larger singlets. These changes were more predominant where the freezing commenced. The results of the experimental studies were interpreted by developing several mathematical and conceptual models, including freezing-induced pressure model, Darcy's law model, multipore ganglia model (rupture coefficient) and ice-snap off model.