We study the compound flooding processes that occurred in Hurricane Florence (2018), which was accompanied by heavy precipitation, using a 3D creek-to-ocean hydrodynamic model. We examine the important role played by barrier islands in the observed compound surges in the coastal watershed. Locally very high resolution is used in some watershed areas in order to resolve small features that turn out to be critical for capturing the observed high water marks locally. The wave effects are found to be significant near barrier islands and have contributed to some observed over-toppings and breaches. Results from sensitivity tests applying each of the three major forcing factors (oceanic, fluvial, and pluvial) separately are succinctly summarized in a “dominance map” that highlights significant compound effects in most of the affected coastal watersheds, estuaries, and back bays behind the barrier islands. Operational forecasts based on the current model are being set up at NOAA to help coastal resource and emergency managers with disaster planning and mitigation efforts.

Holistic simulation approaches are often required to assess human impacts on a river-estuary-coastal system, due to the intrinsically linked processes of contrasting spatial scales. In this paper, a Semi-implicit Cross-scale Hydroscience Integrated System Model (SCHISM) is applied in quantifying the impact of a proposed hydraulic engineering project on the estuarine hydrodynamics. The project involves channel dredging and land expansion that traverse several spatial scales on an ocean-estuary-river-tributary axis. SCHISM is suitable for this undertaking due to its flexible horizontal and vertical grid design and, more importantly, its efficient high-order implicit schemes applied in both the momentum and transport calculations. These techniques and their advantages are briefly described along with the model setup. The model features a mixed horizontal grid with quadrangles following the shipping channels and triangles resolving complex geometries elsewhere. The grid resolution ranges from ~6.3 km in the coastal ocean to 15 m in the project area. Even with this kind of extreme scale contrast, the baroclinic model still runs stably and accurately at a time step of 2 min, courtesy of the implicit schemes. We highlight that the implicit transport solver alone reduces the total computational cost by 82%, as compared to its explicit counterpart. The base model is shown to be well calibrated, then it is applied in simulating the proposed project scenario. The project-induced modifications on salinity intrusion, gravitational circulation, and transient events are quantified and analyzed.