Compound flooding, the concurrence of multiple flooding mechanisms such as storm surge, heavy rainfall, and riverine flooding, poses a significant threat to coastal communities. To mitigate the impacts of compound flooding, forecasts must represent the variability of flooding drivers over a wide range of spatial scales while remaining timely. One approach to develop these forecasts is through subgrid corrections, which utilize information at smaller scales to “correct” water levels and current velocities averaged over the model scale. Recent studies have shown that subgrid models can improve both accuracy and efficiency; however, existing models are not able to account for the dynamic interactions of hydrologic and hydrodynamic drivers and their contributions to flooding along the smallest flow pathways when using a coarse resolution. Here, we have developed a solver called CoaSToRM (Coastal Subgrid Topography Research Model) with subgrid corrections to compute compound flooding in coastal systems resulting from fluvial, pluvial, tidal, and wind-driven processes. A key contribution is the model’s ability to enforce all flood drivers and use the subgrid corrections to improve the accuracy of the coarse-resolution simulation. The model is validated for Hurricane Eta 2020 in Tampa Bay, showing improved prediction accuracy with subgrid corrections at 42 locations. Subgrid models with coarse resolutions (R2 = 0.70, 0.73, 0.77 for 3-, 1.5-, 0.75-km grids) outperform standard counterparts (R2 = 0.03, 0.14, 0.26). A 3-km subgrid simulation runs roughly 50 times faster than a 0.75-km subgrid simulation, with similar accuracy.
A Begmohammadi, D Wirasaet, N Lin, JC Dietrich, D Bolster, AB Kennedy (2023). “Subgrid Modeling for Compound Flooding in Coastal Systems.” Coastal Engineering Journal, 66(3), 434-451, DOI: 10.1080/21664250.2024.2373482.
The U.S. Atlantic and Gulf of Mexico coasts are vulnerable to storms, which can cause significant erosion of beaches and dunes that protect coastal communities. Real-time forecasts of storm-driven erosion are useful for decision support, but they are limited due to demands for computational resources and uncertainties in dynamic coastal systems and storm forcings. Current methods for coastal change forecasts are based on empirical calculations for wave run-up and conceptual models for erosion, which do not represent sediment transport and morphological change during the storm. However, with continued advancements in high-resolution geospatial data and computational efficiencies, there is an opportunity to apply morphodynamic models for forecasts of beach and dune erosion as a storm approaches the coast. In this study, we implement a forecast system based on a deterministic, dynamic model. The morphodynamic model is initialized with digital elevation models of the most up-to-date conditions and forced with hydrodynamics from wave and circulation model forecasts, and its predictions are categorized based on impact to the primary dune, defined in this study as the first ridge of sand landward of the beach. Results are compared spatially to the observed post-storm topography using changes to dune crest elevations and volumes, and temporally to the predicted total water level at the forecasted moment of dune impact.
With sea-level rise, flooding in coastal communities is now common during the highest high tides. Floods also occur at normal tidal levels when rainfall overcomes stormwater infrastructure that is partially submerged by tides. Data describing this type of compound flooding is scarce and, therefore, it is unclear how often these floods occur and the extent to which non-tidal factors contribute to flooding. We combine measurements of flooding on roads and within storm drains with a numerical model to examine processes that contribute to flooding in Carolina Beach, NC, USA — a community that chronically floods outside of extreme storms despite flood mitigation infrastructure to combat tidal flooding. Of the 43 non-storm floods we measured during a year-long study period, one-third were unexpected based on the tidal threshold used by the community for flood monitoring. We introduce a novel model coupling between an ocean-scale hydrodynamic model (ADCIRC) and a community-scale surface water and pipe flow model (3Di) to quantify contributions from multiple flood drivers. Accounting for the compounding effects of tides, wind, and rain increases flood water levels by up to 0.4 m compared to simulations that include only tides. Setup from sustained (non-storm) regional winds causes deeper, longer, more extensive flooding during the highest high tides and can cause floods on days when flooding would not have occurred due to tides alone. Rainfall also contributes to unexpected floods; because tides submerge stormwater outfalls on a daily basis, even minor rainstorms lead to flooding as runoff has nowhere to drain. As a particularly low-lying coastal community, Carolina Beach provides a glimpse into future challenges that coastal communities worldwide will face in predicting, preparing for, and adapting to increasingly frequent flooding from compounding tidal and non-tidal drivers atop sea-level rise.
Inundation models represent coastal regions with a grid of computational points, often with varying resolution of flow pathways and barriers. Models based on coarse grid solutions of shallow water equations have been improved recently via the use of subgrid corrections, which account for information (ground surface elevations, roughness characteristics) at smaller scales. In this work, numerical approaches of an established storm surge model are extended to include subgrid corrections. In an attempt to maintain continuity with existing users and results, model extensions were limited to those needed to provide basic subgrid capabilities, and included two major additions. First, a finite volume method is used to incorporate corrections to the mass and momentum equations using high-resolution ground surface elevations. Second, the no-slip condition imposed on the B-grid wet/dry interface in the model is modified to a slip condition to enable flows in channels with widths comparable to cell size. Numerical results demonstrate these numerical extensions can significantly enhance the accuracy of the model’s predictions of coastal flooding, with low additional computational cost.
The inland propagation of storm surge caused by tropical cyclones depends on large and small waterways to connect the open ocean to inland bays, estuaries, and floodplains. Numerical models for storm surge require these waterways and their surrounding topography to be resolved sufficiently, which can require millions of computational cells for flooding simulations on a large (ocean scale) computational domain, leading to higher demands for computational resources and longer wall-clock times for simulations. Alternatively, the governing shallow water equations can be modified to introduce subgrid corrections that allow coarser and cheaper simulations with comparable accuracy. In this study, subgrid corrections are extended for the first time to simulations at the ocean scale. Higher-level corrections are included for bottom friction and advection, and look-up tables are optimized for large model domains. Via simulations of tides, storm surge, and coastal flooding due to Hurricane Matthew in 2016, the improvements in water level prediction accuracy due to subgrid corrections are evaluated at 218 observation locations throughout 1500 km of coast along the South Atlantic Bight. The accuracy of the subgrid model with relatively coarse spatial resolution (ERMS = 0.41 m) is better than that of a conventional model with relatively fine spatial resolution (ERMS = 0.67 m). By running on the coarsened subgrid model, we improved the accuracy over efficiency curve for the model, and as a result, the computational expense of the simulation was decreased by a factor of 13.
Dunes and beaches are vulnerable to erosion during storm events. Numerical models can predict beach response to storms with fidelity, but their computational costs, the domain-specific knowledge necessary to use them, and the wide range of potential future storm and beach conditions can hinder their use in forecasting storm erosion for short- and long-term horizons. We develop an emulator, which is an efficient predictive model that behaves like a numerical model, to predict the morphologic response of the subaerial beach to storms. Specific emphasis is placed on providing antecedent beach states as an input to the emulator and predicting the post-storm profile shape. Training data include beach profiles at multiple stages in a nourishment life cycle to assess if such a framework can be applied in locations that nourish as a coastal defense policy. Development and application of the emulator is focused on Nags Head, North Carolina, which nourishes its beaches to mitigate hazards of storm waves, flooding, and erosion. A high-fidelity, process-based morphodynamic model is used to train the emulator with 1250 scenarios of sea-storms and beach profiles. The post-storm beach state is emulated with a parameterized power-law function fit to the eroded portion of the subaerial profile. When the emulator was tested for a sequence of real storms from 2019, the eroded beach profiles were predicted with a skill score of 0.66. This emulator is promising for future efforts to predict storm-induced beach erosion in hazard warnings or adaptation studies.
Predictions of storm surge and flooding require models with higher resolution of coastal regions, to describe fine-scale bathymetric and topographic variations, natural and artificial channels, flow features, and barriers. However, models for real-time forecasting often use a lower resolution to improve efficiency. There is a need to understand how resolution of inland regions can translate to predictive accuracy, but previous studies have not considered differences between models that both represent conveyance into floodplains and are intended to be used in real time. In this study, the effects of model resolution and coverage are explored using comparisons between forecast-ready and production-grade models that both represent floodplains along the US southeast coast, but with typical resolutions in coastal regions of 400 and 50 m, respectively. For two storms that impacted the US southeast coast, it is shown that, although the overall error statistics are similar between simulations on the two meshes, the production-grade model allowed a greater conveyance into inland regions, which improved the tide and surge signals in small channels and increased the inundation volumes between 40% and 60%. Its extended coverage also removed water level errors of 20–40 cm associated with boundary effects in smaller regional models.
Barrier islands are a primary coastal defense and often experience erosion during storms. When they fail due to storm-induced breaching, there can be significant changes to the small- and large-scale hydrodynamics and morphodynamics of the region. In this study, we explore the formation of a breach on Hatteras Island, North Carolina, during Isabel (2003) and the subsequent flooding into Pamlico Sound. Two-way coupling of high-fidelity, high-resolution numerical models for coastal erosion and flooding enables a better understanding of the formation of the breach, as well as scenarios of the breach’s effects on the circulation in the region. The breach connecting the ocean to the sound formed during the day of landfall. It is shown that, during the storm, overwash and inundation from the ocean led to deterioration of the beach and dunes, and then after the storm, the creation of channels through the island was sensitive to elevated water levels in the lagoon. Then flooding scenarios are considered in which the ground surface of the hydrodynamic model was (a) static, updated with the (b) pre- and post-storm observations, and updated dynamically with (c) erosion model predictions and (d) erosion model predictions with elevated lagoon-side water levels. The model results show that the breach has region-scale effects on flooding that extend 10 to 13 km into the lagoon, increasing the local water levels by as much as 1.5 m. These results have implications for similar island-lagoon systems threatened by storms.