The physical setting explains why this region is so difficult to simulate. The Bay of Bengal is the largest semi-enclosed bay in the global ocean, and its northern boundary is closed, so the entire basin responds intensely to the twice-yearly reversal of the monsoon winds. From May to September, strong southwesterly winds dominate; from November to February, weaker northeasterlies take over. During the transitional months, westerly winds over the equatorial Indian Ocean generate disturbances that propagate toward the bay as equatorial and coastal Kelvin waves and westward-moving Rossby waves, reshaping the thermocline and upper-ocean stratification on seasonal to interannual timescales. On top of this remote forcing sits an extraordinary local input: the Ganges–Brahmaputra–Meghna river system discharges more than 4.4 × 10¹¹ cubic meters of freshwater each year, spreading a buoyant, sediment-laden lens across the shelf that dominates the hydrography up to roughly 21.5°N. This freshwater cap produces intense haline stratification, a barrier layer between the halocline and the isothermal layer, and, in winter, temperature inversions in which warm subsurface water becomes trapped between cooler surface and deeper layers.

Those inversions are a signature feature of wintertime hydrography in the northern bay, and the CTD profiles captured them in remarkable detail. Nine of the twelve coastal stations showed temperature inversions, and several stations, including stations 4, 5, 6, 9, and 12, displayed multiple inversion layers stacked at different depths. At stations 5 and 6, a sharp thermal inversion appeared near 15 meters in the in situ data, yet the model completely missed the feature. At stations 9 and 10, small but distinct mid-depth temperature drops were present in the observations but absent from the model output. Station 3 revealed a mixed layer only 4 to 5 meters deep, while the model showed no clear mixed-layer signature at all. The largest discrepancy occurred at station 12, where the observed mixed layer reached nearly 19 meters but the model profile ended around 11 meters. Only a handful of stations, such as stations 2 and 8, showed broadly similar profile shapes between observations and model, and even there the agreement was only approximate.

The coastal observations also revealed how shallow the winter mixed layer really is on this shelf. Mixed-layer depths from the CTD casts ranged from just 4 to 12 meters, with a mean of 5.18 meters at the northeastern stations and 10.31 meters at the northwestern stations. The authors attribute the shallower northeastern mixed layers to stronger near-surface stratification, while the relatively deeper northwestern values suggest somewhat greater vertical mixing or weaker surface stabilization near the river mouths. A time series of mixed-layer depth averaged over the northern bay showed the layer deepening slightly to a peak of about 11.4 meters around 8 February before shoaling again through late February and March. The model, by contrast, barely varied at all, holding steady between 9.8 and 9.9 meters throughout the entire three-month period and capturing none of the observed fluctuation. In station-by-station comparisons, the model generally overestimated mixed-layer depth, exceeding observations by roughly 3 to 5 meters at stations 1, 5, and 6, while underestimating at stations 11 and 12.

Offshore, the picture was somewhat better but still flawed. The team compared three Argo float profiles, located at 17.50°N 88.50°E, 18.50°N 90.50°E, and 19.50°N 92.50°E, against the model for January through March. Below about 80 to 100 meters, the model performed reasonably well, reproducing the depth of the well-developed thermocline and the general stratification, and both datasets showed a gradual deepening of the thermocline from January to March. In the upper 50 meters, however, the model diverged sharply from the floats. Argo revealed strong salinity-driven stratification and wintertime temperature inversions that the model did not fully reproduce, and collocated comparisons at 5 meters depth showed a consistent warm near-surface temperature bias alongside a tendency to underestimate salinity. Density differences followed the same pattern, with negative density biases in freshwater-dominated regions. The authors caution that because only three Argo sites were available over a short seasonal window, these bias patterns should be read as preliminary and location-specific rather than a complete regional assessment.

The broader oceanographic context of the 2020 winter helps explain the observed structures. Sea-level anomaly remained negative across the entire northern bay, ranging from about −0.32 to −0.36 meters in January and February and deepening to roughly −0.44 meters by March, with the strongest negative anomalies concentrated near the mouths of the Ganges–Brahmaputra system. Surface currents showed only weak variability, but subsurface currents told a different story: at 55 meters depth, currents shifted direction considerably and strengthened through the season, while at 21 meters they remained weak and variable. These depth-dependent current changes align closely with the thermal features seen in vertical sections, including warm subsurface layers of 25.4 to 26.6 °C trapped between cooler 23 to 24 °C waters in January and February, a cold intrusion between 90.4°E and 91.6°E that strengthened into March, and a downward penetration of warm water between 89.4°E and 90.2°E in February. Hovmöller diagrams confirmed a sudden cooling of roughly 2 to 2.5 °C around 21.3 to 21.5°N in February, followed by broad warming to about 26 °C in March as the winter monsoon weakened and wind stress declined from peaks of 0.04 to 0.054 pascals in the northeastern bay.

Why does a global model with a roughly 9-kilometer horizontal resolution and 5-meter vertical spacing near the surface struggle so much here? The authors point to a combination of factors rather than a single culprit. The model’s discrete vertical levels inherently smooth thin mixed layers, sharp thermohaline gradients, and narrow inversion layers that the CTD, sampling at roughly 0.2-meter intervals, resolves easily. Its global configuration and tripolar ORCA12 grid may not adequately represent regional subsurface processes, and the satellite datasets used to set initial conditions are often poorly suited to coastal and estuarine zones, forcing the use of pseudo-observation methods that can miss delicate features. Complex deltaic bathymetry, tide–river interactions, and the sheer magnitude of freshwater discharge add further difficulty. Importantly, the authors note that the region’s observational scarcity itself weakens the data assimilation that models depend on, although their analysis does not establish a direct causal link between sparse data and model bias. The result is a model that treats the northern shelf as an averaged, interpolated approximation of a far more intricate reality.

The observational gap itself is sobering. Despite the Argo program’s more than two million temperature–salinity profiles since 1999, fewer than 50 collocated Argo observations were reported for the extreme northern bay between 2015 and 2019, because the floats simply cannot sample the shallow Bangladesh shelf where freshwater input, sediment loading, and boundary currents are strongest. The in situ record consists almost entirely of short campaign-based surveys: a 2018 Nansen cruise with 38 CTD stations, a winter 2016 campaign with 15 casts, the 2020 dataset of 12 profiles used here, and roughly 60 nearshore stations near Cox’s Bazar collected in 2021–2022. To the authors’ knowledge, none of these CTD data are openly accessible, making their twelve profiles, to their knowledge, the only open-source coastal CTD data available for the region. Satellites offer no easy substitute, since they measure mainly the surface and often lack reliable algorithms for optically complex coastal waters, and validating satellite products requires exactly the kind of in situ data that is missing.

The stakes extend well beyond academic model evaluation. Two of the four major fishing grounds of the Bay of Bengal lie along the northeastern Bangladesh coast, suggesting a possible connection between episodic upwelling, nutrient supply, and fisheries distribution, yet the relationships among monsoon forcing, stratification, and nutrient delivery remain poorly understood precisely because strong stratification may limit vertical nutrient transfer. The bay’s freshwater-dominated upper ocean also shapes air–sea heat fluxes, cyclone response, and the monsoon system itself, so errors in simulating its mixed layer and barrier layer can propagate into regional climate predictions. The authors argue that continuous, large-scale, integrated monitoring of the northern bay is essential, not only for characterizing its subsurface structure but for building a regional model capable of accurate forecasting. Until such sustained observations exist, they conclude, the northern Bay of Bengal will remain what their title calls it: a blind spot, where even the best global models cannot see the delicate, fast-changing structures that define this extraordinary shelf sea.

Subject of Research: Winter upper-ocean hydrography and model performance in the data-scarce northern Bay of Bengal

Article Title: Inadequate in situ observations make the northern Bay of Bengal a blind spot for models

Article References: Inadequate in situ observations make the northern Bay of Bengal a blind spot for models. (n.d.). https://doi.org/10.1007/s44289-026-00148-y

Image Credits: AI Generated

DOI: 10.1007/s44289-026-00148-y

Keywords: Bay of Bengal, oceanography, CTD profiles, Argo floats, mixed layer depth, thermal inversion, barrier layer, monsoon, ocean modeling, CMEMS reanalysis, freshwater stratification, in situ observations