© 2025 Lupiola, Bárcena, García-Alba and García. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Frontiers in Marine Science, 2025, 12, 1531684
Estuarine mixing and stratification
Potential energy anomaly budget
Small estuary
Estuarine hydrodynamics
Derivative of the potential energy anomaly
Driving mechanism
3D hydrodynamic model
The environmental variability of rivers and tides create complex mixing patterns, which modulate the estuarine ecosystem services. Therefore, a thorough understanding of these systems is vital, not only for their protection but also for their recovery. This study first applies a method to analyze the different mechanisms driving the mixing and stratification of the water column in Suances estuary, a small estuary with large intertidal zones, by using numerical modeling to calculate the potential energy anomaly (f) and its derivative (ft). Numerical results show that Suances estuary presents an ellipse of mixing and stratification variability driven by, firstly, the river flow (seasonal cycle monthly time scale), secondly, the tidal phase (ebb-flood cycle diurnal time scale) and, lastly, the tidal magnitude (spring-neap cycle fortnightly time scale). Furthermore, these results explicitly highlight how the driving mechanisms can vary for the same estuary geometry at different locations due to diurnal, fortnightly and seasonal changes in forcing. The predominant driving mechanisms in Suances estuary are straining (S) tending to stratify the estuary, river and tide advection (A) tending to mix it and non-linear straining (N), caused by contributions from intertidal areas that favor mixing or stratification according to the tidal cycles. Additionally, a threshold was found between the potential energy anomaly and the depth of the water column, confirming that there are limiting values of the potential energy anomaly depending on the depth that can develop in the estuary. This is especially significant in small and shallow estuaries, since maximum values of the potential energy anomaly will be obtained as a function of depth.