Predicting site-specific storm wave run-up

Storm wave run-up causes beach erosion, wave overtopping, and street flooding. Extreme runup estimates may be improved, relative to predictions from general empirical formulae with default parameter values, by using historical storm waves and eroded profiles in numerical runup simulations. A climatology of storm wave run-up at Imperial Beach, California is developed using the numerical model SWASH, and over a decade of hindcast spectral waves and observed depth profiles. For use in a local flood warning system, the relationship between incident wave energy spectra E(f) and SWASH-modeled shoreline water levels is approximated with the numerically simple integrated power law approximation (IPA). Broad and multi-peaked E(f) are accommodated by characterizing wave forcing with frequency-weighted integrals of E(f). This integral approach improves runup estimates compared to the more commonly used bulk parameterization using deep water wave height H0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_0$$\end{document} and deep water wavelength L0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_0$$\end{document} Hunt (Trans Am Soc Civ Eng 126(4):542–570, 1961) and Stockdon et al. (Coast Eng 53(7):573–588, 2006. https://doi.org/10.1016/j.coastaleng.2005.12.005). Scaling of energy and frequency contributions in IPA, determined by searching parameter space for the best fit to SWASH, show an H0L0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_0L_0$$\end{document} scaling is near optimal. IPA performance is tested with LiDAR observations of storm run-up, which reached 2.5 m above the offshore water level, overtopped backshore riprap, and eroded the foreshore beach slope. Driven with estimates from a regional wave model and observed βf\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _f$$\end{document}, the IPA reproduced observed run-up with <30%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$<30\%$$\end{document} error. However, errors in model physics, depth profile, and incoming wave predictions partially cancelled. IPA (or alternative empirical forms) can be calibrated (using SWASH or similar) for sites where historical waves and eroded bathymetry are available.