Barrier islands are susceptible to erosion, overwash, and breaching during intense storms. However, these processes are not represented typically in large-domain models for storm surge and coastal inundation. In this study, we explore the requirements for bridging the gap between dune-scale morphodynamic and region-scale flooding models. A high-resolution XBeach model is developed to represent the morphodynamics during Hurricane Isabel (2003) in the North Carolina (NC) Outer Banks. The model domain is extended to more than 30 km of Hatteras Island and is thus larger than in previous studies. The predicted dune erosion is in good agreement with post-storm observed topography, and an ‘‘excellent’’ Skill Score of 0.59 is obtained on this large domain. Sensitivity studies show the morphodynamic model accuracy is decreased as the mesh spacing is coarsened in the cross-shore direction, but the results are less sensitive to the alongshore resolution. A new metric to assess model skill,
Water Overpassing Area (
WOA), is introduced to account for the available flow pathway over the dune crest. Together, these findings allow for upscaled parameterizations of erosion in larger-domain models. The updated topography, obtained from XBeach prediction, is applied in a region-scale flooding model, thus allowing for enhanced flooding predictions in communities along the Outer Banks. It is found that, even using a fixed topography in region-scale model, the flooding predictions are improved significantly when post-storm topography from XBeach is implemented. These findings can be generalized to similar barrier island systems, which are common along the U.S. Gulf and Atlantic coasts.
A Gharagozlou, JC Dietrich, A Karanci, RA Luettich, MF Overton (2020). “Storm-driven erosion and inundation of barrier islands from dune- to region-scales.” Coastal Engineering, 158, 103674, DOI: 10.1016/j.coastaleng.2020.103674
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.
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.
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.
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.
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.
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.
Predicting building damage as a function of hurricane hazards, building attributes, and the interaction between hazard and building attributes is a key to understanding how significant interaction reflects variation hazard intensity effect on damage based on building attribute levels. This paper develops multihazard hurricane fragility models for wood structure homes considering interaction between hazard and building attributes. Fragility models are developed for ordered categorical damage states (DS) and binary collapse / no collapse. Exterior physical damage and building attributes from rapid assessment in coastal Mississippi following Hurricane Katrina (2005), high-resolution numerical hindcast hazard intensities from the Simulating WAves Nearshore and ADvanced CIRCulation (SWAN+ADCIRC) models, and base flood elevation values are used as model input. Leave-one-out cross-validation (LOOCV) is used to evaluate model prediction accuracy. Eleven and forty-nine combinations of global damage response variables and main explanatory variables, respectively, were investigated and evaluated. Of these models, one DS and one collapse model met the rejection criteria. These models were refitted considering interaction terms. Maximum 3-s gust wind speed and maximum significant wave height were found to be factors that significantly affect damage. The interaction between maximum significant wave height and number of stories was the significant interaction term for the DS and collapse models. For every 0.3 m (0.98 ft) increase in maximum significant wave height, the estimated odds of being in a higher rather than in a lower damage state for DS model were found to be 1.95 times greater for one- rather than for two-story buildings. For every 0.3 m (0.98 ft) increase in maximum significant wave height, the estimated odds of collapse were found to be 2.23 times greater for one- rather than for two-story buildings. Model prediction accuracy was 84% and 91% for DS and collapse models, respectively. This paper does not consider the full hazard intensity experienced in Hurricane Katrina; rather, it focuses on single-family homes in a defined study area subjected to wind, wave, and storm surge hazards. Thus, the findings of this paper are not applicable for events with hazards that exceed those experienced in the study area, from which the models were derived.