Scientists Build Powerful New Model to Tame the Storm-Struck Bay of Bengal

The northern Bay of Bengal is one of the most dangerous stretches of ocean on Earth, a funnel-shaped sea where some of the deadliest tropical cyclones in recorded history have pushed lethal walls of water into the low-lying delta of Bangladesh. Now, researchers have delivered a new weapon in the fight against that menace: a two-dimensional coupled wave-hydrodynamic model that, for the first time for this region, rigorously calibrates and validates water levels, significant wave heights, and current velocities together within a single unified modeling framework. The work, published in Discover Oceans, offers a statistically robust baseline for understanding how tides, waves, and currents interact along one of the world’s most densely populated and hazard-exposed coastlines.

Led by Md. Khairul Amin and G. M. Jahid Hasan of the Department of Civil Engineering at the Military Institute of Science and Technology in Dhaka, the study tackles a notorious modeling gap. Previous efforts along the Bangladesh coast typically validated either water level or wave height in isolation, leaving current velocity largely untested. The new model, built on the widely used Delft3D suite, closes that gap by coupling a hydrodynamic component with the SWAN spectral wave model and confronting both against real observations at multiple coastal stations. The payoff is impressive: Pearson correlation coefficients ranging from 0.79 to 0.99 and Nash-Sutcliffe efficiency values of 0.72 to 0.99 across all tested parameters.

The challenge the model confronts is extraordinary. The Bay of Bengal’s distinctive geometry, a wide continental shelf in the west that narrows sharply toward the southeast, combined with the enormous freshwater discharge of the Ganges-Brahmaputra-Meghna river system, creates a uniquely complicated hydrodynamic environment. Bangladesh’s coastal zone covers roughly 47,201 square kilometers and is home to approximately 46 million people, much of it sitting only a few meters above mean sea level. Tides here are predominantly semi-diurnal and amplify toward the head of the bay, while monsoon winds, river discharge, and tropical cyclones all modulate water levels, wave activity, and current patterns in ways that defy simple prediction.

Technically, the model domain stretches from 18 degrees north to 23 degrees north and from 83 degrees east near Vishakhapatnam, India, to 95 degrees east near Gwa beach in Myanmar. The hydrodynamic component employs a curvilinear-orthogonal grid of roughly 27,000 cells with variable resolution, refined to 100-200 meters along the intricately sculpted coastline near the Meghna Estuary and coarsened to 10-15 kilometers near the open ocean boundaries. Bathymetry was drawn from the GEBCO 2020 dataset at 15 arc-second resolution, while tidal forcing at the open boundaries came from the TPXO 8.0 global inverse tide model using 13 principal harmonic constituents. Upstream river boundaries were driven by hourly discharge measurements from the Padma and Upper Meghna rivers.

Atmospheric forcing derives from the European Centre for Medium-Range Weather Forecasts’ ERA5 reanalysis, with hourly 10-meter wind fields and mean sea level pressure at 0.25-degree resolution. The SWAN wave component runs on a separate rectangular grid and incorporates third-generation physics, including depth-induced breaking via the Battjes-Janssen formulation, JONSWAP bottom friction, nonlinear triad interactions, and a phase-decoupled refraction-diffraction approximation. Crucially, the researchers employed online two-way coupling between the wave and flow models, allowing radiation stress gradients, enhanced bed shear stress, wave boundary layer streaming, turbulence mixing, and Stokes drift to be exchanged dynamically between components at every coupling interval, capturing the feedback between waves and currents that is essential during energetic monsoon conditions.

Calibration against observed tide levels and ERA5 wave heights during June 2020, a period spanning a full spring-neap tidal cycle under active monsoon conditions, produced tide-level correlations of 0.88 to 0.99 and Nash-Sutcliffe values of 0.91 to 0.99. Independent validation for October 2020 under contrasting post-monsoon conditions confirmed the model’s robustness, with tide-level correlations holding between 0.86 and 0.99. Perhaps most notably, simulated current velocities at Payra Port achieved a correlation coefficient and Nash-Sutcliffe efficiency of 0.80 each, reasonably capturing the temporal rhythm and dominant directions of the semi-diurnal tidal currents, though the authors note a strongly negative percentage bias that reflects the metric’s known sensitivity to reversing tidal flows rather than a genuine failure of the model.

A year-long simulation for 2020 then revealed the striking spatial texture of the region’s coastal hydrodynamics. At Sandwip, an island adjacent to the tide-dominated Meghna Estuary, water levels ranged from about minus 1.8 meters to nearly plus 2.0 meters relative to mean sea level, and depth-averaged currents commonly reached 1.0 to 1.3 meters per second. At Kuakata on the more open south-central coast, the tidal range was somewhat smaller and currents generally stayed below 0.5 meters per second during neap tides, peaking near 1.0 to 1.1 meters per second at spring tide. Wave heights at Sandwip generally remained below 1.0 meter with a late-August peak of roughly 1.7 meters, while Kuakata’s waves stayed between 0.3 and 0.7 meters most of the time, with mean wave periods of 6 to 11 seconds that are longer and more variable than Sandwip’s estuarine 4.5 to 6.5 seconds.

The model also reproduced recognizable large-scale patterns, including the persistent anticlockwise residual circulation around Sandwip Island documented in earlier studies and the pronounced tidal asymmetry near major estuarine systems. Snapshot comparisons showed flood-tide currents reaching 0.69 meters per second and ebb currents 0.38 meters per second during the monsoon neap-tide snapshot, compared with up to 1.11 and 0.90 meters per second respectively during the post-monsoon spring-tide snapshot, with the strongest currents consistently concentrated in shallow nearshore zones and estuarine channels where bathymetry and channel geometry amplify flow. Mean wave directions clustered predominantly between 160 and 210 degrees, a southerly to southwesterly approach consistent with the monsoon wind field and swell arriving from the southern Indian Ocean. The authors are careful to note that the two snapshots differ simultaneously in tidal phase and season, so they illustrate spatial structure rather than isolating a pure seasonal effect.

The study is candid about its limitations. GEBCO bathymetry cannot fully resolve the fine-scale estuarine channels and shoals of the Meghna Estuary, and ERA5’s relatively coarse resolution may smooth localized nearshore wind gradients and coastline-induced variability, meaning wave validation is effectively a comparison against reanalysis fields rather than fully independent observations. Current-velocity validation rests on a single station, and the model is accordingly best suited for regional-scale assessment of tide-wave-current variability rather than absolute sediment-flux estimates or site-specific engineering predictions. The researchers recommend future evaluation against satellite altimetry wave products and in situ buoy records, along with higher-resolution hydrographic surveys and nested local models for areas of specific engineering interest.

Even with those caveats, the significance for Bangladesh is hard to overstate. A physically consistent, statistically validated coupled model of the northern Bay of Bengal provides exactly the platform needed for storm surge forecasting, disaster preparedness, early warning systems, coastal infrastructure design, and long-term climate change impact assessment in a region where sea-level rise and intensifying cyclones threaten millions. By demonstrating that a single coupled framework can simultaneously simulate tides, waves, and currents with good-to-excellent skill, the Dhaka team has established a regional baseline onto which future studies of compound flooding, sediment transport, morphodynamics, and coastal protection interventions can build. For the 46 million people living along this vulnerable deltaic edge, better models mean better warnings, and better warnings mean lives saved.

The statistical metrics used to judge the model deserve a brief explanation for readers outside the field. Pearson’s correlation coefficient measures how closely simulated and observed values move together in time, while the Nash-Sutcliffe efficiency compares the model’s errors against the simple baseline of using the observed mean; a value of one indicates a perfect match, and values above roughly 0.7 are generally considered good for coastal hydrodynamic simulations. The fact that the coupled framework achieved efficiencies as high as 0.99 for tide levels, across both monsoon and post-monsoon conditions, indicates that the model captures not only the timing of the semi-diurnal tidal signal but also its amplitude through spring-neap cycles.

The choice of a full annual simulation spanning 2020 is itself methodologically significant. Many earlier regional studies focused on single cyclone events or short calibration windows, which can overstate skill because a model tuned to one tidal or wind regime may fail when conditions shift. By testing against the June monsoon period, when southwesterly winds and heavy Ganges-Brahmaputra-Meghna discharge dominate, and again in October as conditions relaxed, the researchers exposed the model to the two most contrasting hydroclimatic states the coast experiences.

The two-way coupling approach also reflects a broader trend in coastal oceanography. Earlier generations of surge models treated waves and water levels as separate problems, but studies of cyclone events in the Bay of Bengal have repeatedly shown that wave setup and radiation stress gradients can materially alter coastal water levels during storms. By exchanging radiation stresses, bed shear stress, and turbulence terms between the wave and flow components continuously, the new framework provides a physically consistent foundation for the compound-flooding and storm-surge applications that regional disaster planners increasingly demand.

Subject of Research: Coupled wave-hydrodynamic modeling of tides, waves, and currents in the northern Bay of Bengal

Article Title: Two-dimensional coupled wave-hydrodynamic modeling of the northern Bay of Bengal

Article References: Amin, M. K., & Hasan, G. M. J. (2026). Two-dimensional coupled wave-hydrodynamic modeling of the northern Bay of Bengal. Discover Oceans, 3(1), Article 54. https://doi.org/10.1007/s44289-026-00167-9

Image Credits: AI Generated

DOI: 10.1007/s44289-026-00167-9

Keywords: Bay of Bengal, Bangladesh coast, Delft3D, SWAN wave model, tidal hydrodynamics, storm surge, coastal modeling, Meghna Estuary, monsoon, ERA5 reanalysis, current velocity, climate resilience