We investigate two-dimensional advective stirring by satellite-derived ocean surface currents in the north Bay of Bengal during December 2015–March 2016 using the framework of Lagrangian Coherent Structures (LCS). The forward and backward finite-time Lyapunov exponent (fFTLE and bFTLE, respectively) fields indicate that freshwater parcels from the river run-off are subjected to intense stirring by the ocean surface currents, which contain both geostrophic and wind-forced Ekman currents. Large-scale thermal fronts present in the satellite-based sea surface temperature (SST) fields are then shown to closely resemble attracting LCS (aLCS) that govern horizontal advective stirring on subseasonal time scales. Statistical estimates of the angle between the thermal fronts (of all spatial scales present in the SST data) and the aLCS from the entire north Bay of Bengal show a preference to align with each other. Significantly, in the presence of sufficiently strong aLCS in their neighbourhood, thermal fronts are shown to have a strong preference to align with them. Finite-time advection of fluid particles over around a week to compute the bFTLE fields results in the strongest alignment between the thermal fronts and neighbouring aLCS. The robustness of our results are established by presenting analyses of (i) geostrophic currents from altimeter data, and (ii) a reanalysis-based fine-scale salinity field. Attracting LCS, identified as sufficiently strong ridges in the bFTLE fields, are argued to represent reasonable estimates of the location and orientation of fronts (thermal or salinity) formed by advective stirring over subseasonal time scales. © 2019 Elsevier Ltd