Recent observational studies provided preliminary insights on the interannual variability of Bay of Bengal (BoB) Sea Surface Salinity (SSS), but are limited by the poor data coverage. Here, we describe the BoB interannual SSS variability and its driving processes from a regional eddy-permitting ocean general circulation model forced by interannually varying air-sea fluxes and altimeter-derived discharges of major rivers over the past two decades. Simulated interannual SSS variations compare favorably with both in situ and satellite data and are largest in boreal fall in three regions: the northern BoB, the coastal region off east India, and the Andaman Sea. In the northern BoB, these variations are independent from those in other regions an...