Total water levels (TWLs), including the contribution of wind waves, associated with tropical cyclones (TCs) are among the most damaging hazards faced by coastal communities. TC-induced economic losses are expected to increase because of stronger TC intensity, sea level rise, and increased populations along the coasts. TC intensity, translation speed, and distance to the coast affect the magnitude and duration of increased TWLs and wind waves. Under climate change, the proportion of high-intensity TCs is projected to increase globally, whereas the variation pattern of TC translation speed also depends on the ocean basin and latitude. There is an urgent need to improve our understanding of the linkages between TC characteristics and TWL components. In the past few years, hurricanes Matthew (2016), Dorian (2019), and Isaias (2020) propagated over the South Atlantic Bight (SAB) with similar paths but resulted in different coastal impacts. We combined in situ observations and numerical simulations with the Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) modeling system to analyze the extreme TWLs under the three TCs. Model verification showed that the TWL components were well reproduced by the present model setup. Our results showed that the peak storm surge and the peak wave runup depended mainly on the TC intensity, the distance to the TC eye, and the TC heading direction. A decrease in TC translation speed primarily led to longer exceedance durations of TWLs, which may result in more severe economic losses. Wave-dependent water level components (i.e., wave setup and wave swash) were found to dominate the peak TWL within the near-TC field. Our results also showed that in specific conditions, the prestorm wave runup associated with the TC-induced swell may lead to TWLs higher than at the peak of the storm. This was the case along the SAB during Hurricane Isaias. Isaias's fast TC translation speed and the fact that its swell was not blocked by any islands were the main factors contributing to these peak TWLs ahead of the storm peak.

An immersed boundary‐finite element with soft‐body dynamics has been implemented to study steady flow over a finite patch of submerged flexible aquatic vegetation. The flow structure interaction model can resolve the flow interactions with flexible vegetation, and hence the reconfiguration of vegetation blades to ambient flow. Flow dynamics strongly depend on two dimensionless parameters, namely vegetation density and Cauchy number (defined as the ratio of the fluid drag force to the elastic force). Five different flow patterns have been identified based on vegetation density and Cauchy number, including the limited reach, swaying, “monami” A, “monami” B with slow moving interfacial wave, and prone. The “monami” B pattern occurred at high vegetation density and is different from “monami” A, in which the passage of Kelvin‐Helmholtz billows strongly affects the vegetation interface. With soft‐body dynamics, blade‐to‐blade interactions can also be resolved. At high vegetation density, the hydrodynamic interactions play an important role in blade‐to‐blade interactions, where adjacent vegetation blades interact via the interstitial fluid pressure. At low vegetation density, direct contacts among vegetation blades play important roles in preventing unphysical penetration of vegetation blades.

Cohesive sediment forms flocs of various sizes and structures in the natural turbulent environment. Understanding flocculation is critical in accurately predicting sediment transport and biogeochemical cycles. In addition to aggregation and breakup, turbulence also reshapes flocs toward more stable structures. An Eulerian–Lagrangian framework has been implemented to investigate the effect of turbulence on flocculation by capturing the time-evolution of individual flocs. We have identified two floc reshaping mechanisms, namely breakage-regrowth and restructuring by hydrodynamic drag. Surface erosion is found to be the primary breakup mechanism for strong flocs, while fragile flocs tend to split into fragments of similar sizes. Aggregation of flocs of sizes comparable to or greater than the Kolmogorov scale is modulated by turbulence with lower aggregation efficiency. Our findings highlight the limiting effects of turbulence on both floc size and structure.

Considerable uncertainty exists about how to represent wave bottom friction in coastal systems like reefs where orbital excursions are similar to roughness element size. Here, interactions between waves and large bottom roughness were investigated using Large Eddy Simulations of oscillatory flow over infinite hemisphere arrays. Wave amplitude, period, and hemisphere spacing were varied to investigate the dependence of kinematics and dynamics on dimensionless parameters. The net effect of topography on the oscillatory flow was assessed using a spatially and phase‐averaged Navier‐Stokes framework. Dynamics depended strongly on Keulegan‐Carpenter number (KC), the ratio of wave orbital excursion to roughness element size. For, flow separation was weak, form drag was small, stress gradients were negligible, and the main effect of topography on the flow was the inertial force associated with acceleration around roughness elements. For, strong flow separation occurred, and both drag and inertial forces were important. Phase‐dependent dispersive stresses were the main mechanism for vertical momentum transfer between the canopy layer and the overlying water column. Friction factors based on the drag force increased withKCfor, different from previously proposed empirical curves, but approached these curves for highKC. Friction factors based on the total force decreased with increasingKCand were consistent with previously proposed curves. These results highlight the importance of distinguishing the total force on the bottom, the drag force that removes energy from the flow, and the shear stress above the canopy layer, which were very different for the parameter range in this study.

In the coastal ocean, interactions of waves and currents with large roughness elements, similar in size to wave orbital excursions, generate drag and dissipate energy. These boundary layer dynamics differ significantly from well-studied small-scale roughness. To address this problem, we derived spatially and phase-averaged momentum equations for combined wave–current flows over rough bottoms, including the canopy layer containing obstacles. These equations were decomposed into steady and oscillatory parts to investigate the effects of waves on currents, and currents on waves. We applied this framework to analyse large-eddy simulations of combined oscillatory and steady flows over hemisphere arrays (diameter $D$ ), in which current ( $U_c$ ), wave velocity ( $U_w$ ) and period ( $T$ ) were varied. In the steady momentum budget, waves increase drag on the current, and this is balanced by the total stress at the canopy top. Dispersive stresses from oscillatory flow around obstacles are increasingly important as $U_w/U_c$ increases. In the oscillatory momentum budget, acceleration in the canopy is balanced by pressure gradient, added-mass and form drag forces; stress gradients are small compared to other terms. Form drag is increasingly important as the Keulegan–Carpenter number $KC=U_wT/D$ and $U_c/U_w$ increase. Decomposing the drag term illustrates that a quadratic relationship predicts the observed dependences of steady and oscillatory drag on $U_c/U_w$ and $KC$ . For large roughness elements, bottom friction is well represented by a friction factor ( $f_w$ ) defined using combined wave and current velocities in the canopy layer, which is proportional to drag coefficient and frontal area per unit plan area, and increases with $KC$ and $U_c/U_w$ .