Mixing of passive tracers in the Bay of Bengal, driven by altimetry derived daily geostrophic surface currents, is studied on subseasonal timescales. To begin with, Hovmöller plots, wavenumber-frequency diagrams and power spectra confirm the multiscale nature of the flow. Advection of latitudinal and longitudinal bands immediately brings out the chaotic nature of mixing in the Bay via repeated straining and filamentation of the tracer field. A principal finding is that mixing is local, i.e., of the scale of the eddies, and does not span the entire basin. Indeed, Finite Time Lyapunov Exponent (FTLE), Relative Dispersion (RD) and Finite Size Lyapunov Exponents (FSLE) maps in all seasons are patchy with minima scattered through the interior of the Bay. Further, FTLE, FSLE and RD maps show that the Bay experiences a seasonal cycle wherein rapid stirring progressively moves from the northern to southern Bay during pre and post monsoonal periods, respectively. The non-uniform stirring of the Bay is reflected in long tailed histograms of FTLEs, that become more stretched for longer time intervals. Quantitatively, advection for a week shows the mean FTLE lies near 0.15-0.16 day −1 , while extremes reach almost 0.5 day −1 . Averaged over the Bay, RD initially grows exponentially, this is followed by a powerlaw at scales between approximately 100 and 250 km, which finally transitions to an eddy-diffusive regime. These findings are confirmed by FSLEs; in addition, quantitatively, below 250 km, a scale dependent diffusion coefficient is extracted that behaves as a power-law with cluster size, while above 250 km, eddy-diffusivities range from 6 × 10 3 -10 4 m 2 /s. Finally, in concert with satellite salinity data, these Lagrangian tools are used to analyse a single post-monsoonal fresh water mixing event. Here, while stirring the salinity field at large scales, FTLEs and FSLEs allow the identification of transport barriers, and elucidate how individual eddies help preserve the identity of fresh water.