These analyses were performed by Carter Day, an undergraduate researcher in our team.
Like other spectral wave models, SWAN uses parameterizations to represent sources (e.g. wind), sinks (e.g. whitecapping, bottom friction, depth-limited breaking), and resonance (e.g. quadruplet and triad wave-wave interactions). Each parameterization is based on laboratory and experimental data to represent the transfer of energy to, from, and between waves. It is necessary for each term to represent its physical process, but it is also necessary for the terms to be calibrated collectively to represent their combined effects on wave evolution.
SWAN release version 41.31 was modified in two main ways: derivative computation was changed to use the Green-Gauss formula, and a new ‘package’ of wave physics (the so-called ST6 physics) was introduced. This package includes new parameterizations of wind input, whitecapping, swell dissipation, wind speed scaling, and other processes. The ST6 physics have been adopted by other wave models (e.g. NOAA’s WaveWatch III), and it will likely become the preferred physics package for SWAN. However, because the ST6 physics package has changes to so many parameterizations, it is necessary to quantify its effects on wave predictions during recent storms.
In this study, we simulate two recent hurricanes, Gustav (2008) and Florence (2018), and we compare wave predictions with the new ST6 physics package. Do the ST6 physics improve the SWAN wave predictions?
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.
One of the most unpredictable and deadly parts of a coastal storm is the storm surge, which can cause devastating flooding of coastal regions, and can result in loss of property and life. Storm surge is a result of winds pushing water from the nearshore ocean to rise above regular tide levels. Storm surge can have a short duration; elevated water levels are limited to when the storm winds are strongest at the coast, typically for a few hours as the storm makes landfall. This short duration is a challenge for predictions of when storm surge will start, how long it will persist, and which regions will be flooded.
This research will explore the arrival of storm effects for Florence, which made landfall along the North Carolina coast during September 2018. It is an ideal storm for this research as its track was shore-normal, and thus its coastal effects increased as it approached landfall. This will allow for investigating the most optimal switch by focusing on a single switch between a lower-resolution mesh to a higher-resolution mesh. The switches will be initiated by several triggers, including wind speeds and water levels at the coast and inland locations, and with several lead times, including near and several days before landfall. Model performance will be quantified via comparisons to observations of storm effects in the region, as well as to a single, high-resolution simulation for the full storm. It will be shown that switching from a coarse resolution mesh to a fine resolution mesh will lead to an increase in efficiency gains across all switching simulations with the most optimal switch time resulting in the most accurate predictions of water levels as compared to our full high-resolution simulation.
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.