The treatment of horizontal diffusion in the particle model deserves particular attention, because it illustrates a subtlety that is often overlooked in applied transport modelling. In a random-walk formulation, particles are displaced by both the resolved advective velocity and a stochastic term representing unresolved turbulent mixing. When the diffusion coefficient varies in space, as it does here, a naive implementation in which the stochastic displacement is drawn independently at each step systematically violates the well-mixed condition, the requirement that a uniformly distributed set of particles remains uniformly distributed over time. The consequence is a spurious drift of particles toward regions of low diffusivity, an artifact that could be mistaken for a physical accumulation zone. The correction introduced by Visser, which the authors adopt, adds an explicit gradient term to the particle displacement equation and evaluates the diffusivity at an intermediate position reached after applying half of the gradient drift. This two-stage procedure substantially reduces the bias, and its inclusion in the model design signals a careful attention to numerical fidelity that strengthens confidence in the resulting exposure maps and dispersion statistics.

The decision to represent the lagoon as a closed basin is a deliberate simplification with a clear evidential basis. Previous modelling work on the same system had shown that, even under idealized conditions with both inlets fully open, the exchange of water between the lagoon and the Adriatic Sea is remarkably weak, with the external sea level exerting a negligible influence on the internal circulation except in the immediate vicinity of the inlets themselves. By removing the open boundary altogether, the authors isolate the internal transport mechanisms, wind-driven circulation above all, and can attribute differences in plume behavior between scenarios directly to variations in the wind forcing rather than to changes in boundary exchange. This framing also carries a practical implication: if the lagoon behaves hydrodynamically like a nearly closed basin, then pollutants introduced into it are, to first order, retained within it, and the timescales over which they are redistributed internally become the relevant timescales for ecological exposure. The two inlets, whose operation is informal and unscheduled, therefore emerge not as reliable flushing conduits but as marginal contributors whose management status nonetheless shapes local conditions near the connections.

The six weekly wind scenarios were extracted from a forty-year historical record, which gives the scenario ensemble a degree of representativeness that single-case studies often lack. Rather than simulating one arbitrary week, the authors sampled the range of metocean conditions that the lagoon actually experiences, so that the contrast between scenarios reflects genuine climatological variability rather than an artifact of scenario selection. This approach sits between the two extremes commonly found in the literature: idealized forcing experiments, which are clean but may not correspond to real conditions, and long-term coupled simulations, which are realistic but computationally expensive and harder to interpret scenario by scenario. The weekly window itself is a meaningful choice for a shallow lagoon, since it is long enough for the wind-driven circulation to establish a coherent pattern and for particles to travel appreciable distances, yet short enough that the forcing can be treated as approximately stationary within each scenario.

The distinction between the two release configurations, instantaneous and continuous, addresses a question that is central to pollutant risk assessment but frequently glossed over in modelling studies. A single pulse release, such as an accidental discharge or an overflow event, produces a compact cloud of particles whose fate depends on the flow field encountered at the moment of release. A continuous release, such as the steady discharge from a sewage outfall, produces instead a stream of particles that sample the flow at many different phases of the circulation, and the resulting plume integrates over that variability. The finding that dispersion in some scenarios is sensitive to release timing, with particles either remaining close to the source or spreading more uniformly across the lagoon, indicates that the advective pathways within the lagoon are strongly organized by the wind regime. When the circulation is dominated by coherent, persistent flow structures, the moment of release largely determines which pathway a particle follows, and two releases separated by only hours or days can end up in entirely different regions of the basin. When the circulation is more dispersive, this sensitivity fades and the plume behaves more predictably. The authors attribute this sensitivity specifically to advection, which is a useful diagnostic: it means that the variability arises from the deterministic structure of the flow rather than from the stochastic diffusion term, and it suggests that even modest improvements in the timing of discharges, where such control is possible, could meaningfully reduce exposure in the most retention-prone conditions.

The exposure maps highlighting accumulation zones along the lagoon boundary carry direct management relevance. In a closed or nearly closed basin, material that is not flushed out must eventually be redistributed by the internal circulation, and the boundaries, where currents weaken and particles can strand in shallow areas or along the shoreline, act as natural sinks for floating matter. For floating pollutants, including oils, plastics, and buoyant sewage-associated debris, these boundary accumulation zones represent the locations where ecological contact and cleanup costs are likely to be concentrated. Identifying them in advance, and understanding how their location shifts between wind regimes, allows managers to prioritize monitoring stations, anticipate where visual or chemical impacts will appear after an event, and design interventions such as containment booms or targeted collection with some foreknowledge of where material will converge. The fact that the release point was chosen on the basis of previous field sampling that identified it as the most contaminated area of the lagoon further anchors the simulations in observed reality, ensuring that the modelled plumes originate from a location where pollutant inputs are known to occur rather than from a purely hypothetical source.

The ecological and economic context of the Narta Lagoon amplifies the significance of these results. The lagoon supports migratory birds, European eels, and local aquaculture activities that employ a substantial share of the nearby workforce, so water quality is not merely an environmental concern but a livelihood issue. The documented pressures on the system, including waste disposal serving Vlorë and surrounding towns, salt production, leather tanning, and incomplete wastewater treatment, mean that the pollutant loading simulated in this study is not an abstract scenario but a plausible representation of ongoing inputs. Shallow lagoons are especially sensitive in this respect because their limited depth means that floating material remains in close contact with the biologically productive surface layer, and their restricted exchange means that dilution rates are low. A pollutant that would be rapidly dispersed in an open coastal setting can persist in a lagoon of this type for weeks, repeatedly exposing the same benthic and pelagic communities.

The finding that advective flow plays the dominant role in determining plume behavior also has implications for how monitoring and modelling of such systems should be designed. If mixing were controlled primarily by diffusion, then coarse descriptions of the flow field would suffice, and dispersion could be parameterized with simple bulk coefficients. The dominance of advection implies instead that the spatial structure of the wind-driven circulation, including its gyres, shear zones, and stagnation points, is what controls where pollutants go. This raises the bar for the hydrodynamic model: the velocity fields must resolve the relevant circulation features accurately, because errors in the flow structure translate directly into errors in predicted exposure. It also suggests that wind forecasting, rather than tidal prediction, is the more valuable operational input for anticipating pollutant behavior in this lagoon, a somewhat counterintuitive conclusion for a coastal water body connected to the sea.

The observation that the northern inlet is prone to seasonal closure through sedimentation driven by littoral transport, and is manually reopened by fishermen roughly twice per year, adds a human dimension to the physical problem. The intermittency of this connection means that the lagoon’s effective degree of isolation varies through the year in an unmanaged way, and the southern inlet’s manually operated sluice gates introduce a further element of discretionary control. The authors’ emphasis on the importance of properly managing the two inlets follows logically from their results: if internal retention is already high when the lagoon is treated as closed, then informal or infrequent inlet operation can only prolong stagnation, and a more deliberate management regime, with scheduled openings timed to coincide with conditions favorable to exchange, could reduce the residence time of both water and pollutants. The study’s framing of the lagoon as a closed basin thus serves a dual purpose: it is a defensible modelling simplification supported by prior evidence, and it is also a conservative bounding case whose results quantify the worst-case retention that inlet management would seek to alleviate.

More broadly, the study contributes to a growing body of work demonstrating that semi-enclosed coastal systems cannot be managed by analogy with open coastlines. The combination of high-resolution hydrodynamic modelling, Lagrangian particle tracking, and statistically grounded dispersion descriptors offers a transferable template that can be applied to other shallow lagoons facing similar pressures, many of which are distributed along the Mediterranean and Adriatic coasts where urbanization, tourism, and aquaculture converge on ecologically fragile water bodies. The scenario-based design, in which a small number of representative forcing conditions is analyzed in depth, strikes a practical balance between computational cost and physical realism, and the explicit comparison of release configurations provides a dimension of insight that single-release studies cannot offer. As climate variability alters wind regimes and sea level modifies inlet dynamics in the coming decades, the kind of scenario library constructed here could be extended and reweighted to reflect future conditions, extending the usefulness of the modelling framework well beyond the present analysis of a single Albanian lagoon.

Subject of Research: Lagrangian particle-tracking simulation of pollutant plume dispersion in the Narta Lagoon under variable metocean conditions

Article Title: Assessing lagoon plume dispersion under variable metocean conditions

Article References: Scovenna, M., & Besio, G. (2026). Assessing lagoon plume dispersion under variable metocean conditions. Frontiers in Environmental Engineering, 5, Article 1882143. https://doi.org/10.3389/fenve.2026.1882143

Image Credits: AI Generated

DOI: 10.3389/fenve.2026.1882143

Keywords: coastal lagoon, Lagrangian particle tracking, plume dispersion, metocean conditions, Narta Lagoon, water quality, hydrodynamic modelling, advection, pollutant transport, Delft3D, exposure maps, Albania

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