Stratified and stirred: Monsoon freshwater in the Bay of Bengal

Abstract

Submesoscale ocean dynamics and instabilities, with characteristic scales 0.1–10 km, can play a critical role in setting the ocean’s surface boundary layer thickness and associated density stratification. Submesoscale instabilities contribute to lateral stirring and tracer dispersal. These dynamics are investigated in the Bay of Bengal, motivated by the upper ocean’s potentially coupled interactions with Monsoon winds and convection. The region’s excess precipitation and runoff generates strong salinity gradients that typically set density fronts and stratification in the upper 50 m. Since we cannot synoptically measure currents containing fast-evolving and oscillating components across the submesoscale range, we instead analyze passive tracer distributions (spice ⌘ density-compensated temperature (T) and salinity (S) anomalies), identifying signatures of flows and testing dynamical theories. The analysis is based on over 9000 vertical profiles of T and S measured along ⇠4800 km of ship tracks in the Bay of Bengal during ASIRI and MISO-BOB expeditions in 2013, 2015, and 2018. Observations in the surface mixed layer reveal ⇠1 km scale-selective correlation of surface T and S, with compensation reducing cross-front density gradients by ⇠50%. Using a process study ocean model, we show this is caused by submesoscale instabilities slumping fronts, plus surface cooling over the resultant enhanced salinity stratification, potentially thwarting the forward cascade of energy. In the stratified interior, we present a spectral analysis of horizontal spice variance statistics from wavenumber k ⇠0.01 cpkm to ⇠1 cpkm. At scales <10 km, stratified layers that are closer to the surface exhibit redder passive tracer spectra (power spectra k−3, gradient spectra k−1) than predicted by quasi-geostrophic or frontogenetic theories. Complimentary observations reveal spice patterns with multiple, parallel, ⇠10 m thin layers, crossing isopycnals with O(10−4) slopes, coherent over at least 30–80 km, with coincident layers of stratification anomalies. Comparison with shear measurements, and a numerical process study, suggest that both submesoscale sheared eddies, and thin near-inertial waves, form such layers. Fast formation timescales and large aspect ratios suggest they enhance horizontal mixing by shear dispersion, reducing variance at ⇠1–10 km scales.