Effects of Model Resolution and Coverage on Storm-Driven Coastal Flooding Predictions

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.

A Thomas, JC Dietrich, CN Dawson, RA Luettich (2022). “Effects of Model Resolution and Coverage on Storm-Driven Coastal Flooding Predictions.” Journal of Waterway, Port, Coastal, and Ocean Engineering, 148(1), 04021046, DOI: 10.1061/(ASCE)WW.1943-5460.0000687.

Formation of a Barrier Island Breach and Its Contributions to Lagoonal Circulation

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.

A Gharagozlou, JC Dietrich, TC Massey, DL Anderson, JF Gorski, MF Overton (2021). “Formation of a barrier island breach and its contributions to lagoonal circulation.” Estuarine, Coastal and Shelf Science, 262, 107593, DOI: 10.1016/j.ecss.2021.107593.

Subgrid Corrections in Finite-Element Modeling of Storm-Driven Coastal Flooding

Coastal flooding models are used to predict the timing and magnitude of inundation during storms, both for real-time forecasting and long-term design. However, there is a need for faster flooding predictions that also represent flow pathways and barriers at the scales of critical infrastructure. This need can be addressed via subgrid corrections, which use information at smaller scales to ‘correct’ the flow variables (water levels, current velocities) averaged over the mesh scale. Recent studies have shown a decrease in run time by 1 to 2 orders of magnitude, with the ability to decrease further if the model time step is also increased. In this study, subgrid corrections are added to a widely used, finite-element-based, shallow water model to better understand how they can improve the accuracy and efficiency of inundation predictions. The performance of the model, with and without subgrid corrections, is evaluated on scenarios of tidal flooding in a synthetic domain and a small bay in Massachusetts, as well as a scenario with a real atmospheric forcing and storm surge in southwest Louisiana. In these tests we observed that the subgrid corrections can increase model speed by 10 to 50 times, while still representing flow through channels below the mesh scale to inland locations.

JL Woodruff, JC Dietrich, D Wirasaet, AB Kennedy, D Bolster, Z Silver, SD Medlin, RL Kolar (2021). “Subgrid corrections in finite-element modeling of storm-driven coastal flooding.” Ocean Modelling, 167, 101887, DOI: 10.1016/j.ocemod.2021.101887.

Methodology for Regional Multihazard Hurricane Damage and Risk Assessment

Hurricanes are devastating natural hazards that often cause damage to the built environment as a result of their loadings, which include storm surge, waves, and wind, often in combination. Modeling these hazards individually and their effects on buildings is a complex process because each loading component within the hazard behaves differently, affecting either the building envelope, the structural system, or the interior contents. Realistic modeling of hurricane effects requires a multihazard approach that considers the combined effects of wind, surge, and waves. Previous studies focused primarily on modeling these hazards individually, with less focus on the multihazard impact on the whole building system made up of the combination of the structure and its interior contents. The analysis resolution used in previous studies did not fully enable hurricane risk assessment through a detailed investigation of the vulnerability at the component-level or subassembly-level (a group of components such as interior contents, structural components, or nonstructural components). To address these research gaps, a robust multihazard hurricane risk analysis model that uses high-resolution hazard, exposure, and vulnerability models was developed. This model uses a novel approach to combine the storm surge and wave fragility functions with a suite of existing wind fragilities to account for structural damage and then combines them with another suite of flood-based fragilities to account for interior content damage. The proposed vulnerability model was applied to the state of North Carolina as an example of a regional-scale assessment to demonstrate the ability of the method to predict damage at the building level across this large spatial domain. This model enables better understanding of the damages caused by hurricanes in coastal regions, thereby setting initial post-impact conditions for community resilience assessment and investigation of recovery policy alternatives.

OM Nofal, JW van de Lindt, TQ Do, G Yan, S Hamideh, DT Cox, JC Dietrich (2021). “Methodology for Regional Multihazard Hurricane Damage and Risk Assessment.” Journal of Structural Engineering, 147(11), 04021185, DOI: 10.1061/(ASCE)ST.1943-541X.0003144.

Improving Coastal Flooding Predictions by Switching Meshes during a Simulation

Storm surge and coastal flooding predictions can require high resolution of critical flow pathways and barriers, typically with simulations using grids/meshes with millions of cells/elements to represent a coastal region. However, the cost of this resolution can slow forecasts during a storm. To add resolution when and where it is needed, previous studies have used adaptive mesh methods, which update resolution at single or multiple cells but which require hierarchies of and thresholds for refinement, and nesting methods, which update resolution at subdomains but which require additional simulations. This research proposes a middle way, in which predictions from a coarse mesh are mapped, mid-simulation, onto a fine mesh with increased resolution near the storm’s projected landfall location. The coarse and fine meshes are pre-developed, thus removing any refinement decisions during the simulation, the solution mapping uses a widely used framework, thus enabling an efficient interpolation, and the same simulation is continued, thus eliminating a separate full-domain simulation. For four historical storms, results show efficiency gains of up to 53 percent, with minimal accuracy losses relative to a static simulation.

A Thomas, JC Dietrich, M Loveland, A Samii, CN Dawson (2021). “Improving Coastal Flooding Predictions by Switching Meshes during a Simulation.” Ocean Modelling, 164, 101820, DOI: 10.1016/j.ocemod.2021.101820.

Subgrid Surface Connectivity for Storm Surge Modeling

Subgrid modeling to account for unresolved topography within the context of shallow water equations relies on the use of coarse grids for computational efficiency. However, excessively coarse grids can lead to artificial cross flows between hydrologically disconnected areas separated by physical barriers smaller than the grid size. An approach based on introducing cell and edge clones, consisting of connected groups of pixels in each cell, is able to systematically remove such artificial cross flows. Such an approach considers that the subgrid barriers permanently divide flow among clones and effectively restrict flow to a predetermined path. In this work, a simple algorithm, along with the use of an overtopping formula, is proposed to extend the clone approach to a scenario in which clones are allowed to be further split and merged as needed, depending on the surface elevation during a given runtime. The algorithm is intended for accommodating the possibility of the subgrid barriers being inundated and no-longer dividing the flow during an extreme event. The performance of the proposed algorithm is demonstrated through a series of idealized and more realistic test cases, showing considerable improvements over existing methodologies.

A Begmohammadi, D Wirasaet, Z Silver, D Bolster, AB Kennedy, JC Dietrich (2021). “Subgrid surface connectivity for storm surge modeling.” Advances in Water Resources, 153, 103939, DOI: 10.1016/j.advwatres.2021.103939.

Dynamic Load Balancing for Predictions of Storm Surge and Coastal Flooding

As coastal circulation models have evolved to predict storm-induced flooding, they must include progressively more overland regions that are normally dry, to where now it is possible for more than half of the domain to be needed in none or only some of the computations. While this evolution has improved real-time forecasting and long-term mitigation of coastal flooding, it poses a problem for parallelization in an HPC environment, especially for static paradigms in which the workload is balanced only at the start of the simulation. In this study, a dynamic rebalancing of computational work is developed for a finite-element-based, shallow-water, ocean circulation model of extensive overland flooding. The implementation has a low overhead cost, and we demonstrate a realistic hurricane-forced coastal flooding simulation can achieve peak speed-ups near 45% over the static case, thus operating now at 80−90% efficiency.

KJ Roberts, JC Dietrich, D Wirasaet, WJ Pringle, JJ Westerink (2021). “Dynamic load balancing for predictions of storm surge and coastal flooding.” Environmental Modelling & Software, 140, 105045, DOI: 10.1016/j.envsoft.2021.105045.

Downscaling of Real-Time Coastal Flooding Predictions for Decision Support

During coastal storms, forecasters and researchers use numerical models to predict the magnitude and extent of coastal flooding. These models must represent the large regions that may be affected by a storm, and thus, they can be computationally costly and may not use the highest geospatial resolution. However, predicted flood extents can be downscaled (by increasing resolution) as a post-processing step. Existing downscaling methods use either a static extrapolation of the flooding as a flat surface, or rely on subsequent simulations with nested, full-physics models at higher resolution. This research explores a middle way, in which the downscaling includes simplified physics to improve accuracy. Using results from a state-of-the-art model, we downscale its flood predictions with three methods: (1) static, in which the water surface elevations are extrapolated horizontally until they intersect the ground surface; (2) slopes, in which the gradient of the water surface is used; and (3) head loss, which accounts for energy losses due to land cover characteristics. The downscaling methods are then evaluated for forecasts and hindcasts of Hurricane Florence (2018), which caused widespread flooding in North Carolina. The static and slopes methods tend to over-estimate the flood extents. However, the head loss method generates a downscaled flooding extent that is a close match to the predictions from a higher-resolution, full-physics model. These results are encouraging for the use of these downscaling methods to support decision-making during coastal storms.

CA Rucker, N Tull, JC Dietrich, TE Langan, H Mitasova, BO Blanton, JG Fleming, RA Luettich Jr (2021). “Downscaling of Real-Time Coastal Flooding Predictions for Decision Support.” Natural Hazards, 107, 1341-1369, DOI: 10.1007/s11069-021-04634-8.