We assess erosion and flooding risk in the northern Gulf of Mexico by identifying interdependencies among oceanographic drivers and probabilistically modeling the resulting potential for coastal change. Wave and water level observations are used to determine relationships between six hydrodynamic parameters that influence total water level and therefore erosion and flooding, through consideration of a wide range of univariate distribution functions and multivariate elliptical copulas. Using these relationships, we explore how different our interpretation of the present‐day erosion/flooding risk could be if we had seen more or fewer extreme realizations of individual and combinations of parameters in the past by simulating 10,000 physically and statistically consistent sea‐storm time series. We find that seasonal total water levels associated with the 100 year return period could be up to 3 m higher in summer and 0.6 m higher in winter relative to our best estimate based on the observational records. Impact hours of collision and overwash—where total water levels exceed the dune toe or dune crest elevations—could be on average 70% (collision) and 100% (overwash) larger than inferred from the observations. Our model accounts for non‐stationarity in a straightforward, non‐parametric way that can be applied (with little adjustments) to many other coastlines. The probabilistic model presented here, which accounts for observational uncertainty, can be applied to other coastlines where short record lengths limit the ability to identify the full range of possible wave and water level conditions that coastal mangers and planners must consider to develop sustainable management strategies.

Predictions of total water levels, the elevation of combined tides, surge, and wave runup at the shoreline, are necessary to provide guidance on potential coastal erosion and flooding. Despite the importance of early warning systems for these hazards, existing real-time meteorological and oceanographic forecast systems at regional and national scales, until now, have lacked estimates of runup necessary to predict wave-driven overwash and erosion. To address this need, we present an approach that includes wave runup in an operational, national-scale modeling system. Using this system, we quantify the contribution of waves to potential dune erosion events along 4,700 km of U.S. Atlantic and Gulf of Mexico sandy coastlines for a one-year period. Dune erosion events were predicted to occur at over 80% of coastal locations, where waves dominated shoreline total water levels, representing 73% of the signal. This shows that models that neglect the wave component underestimate the hazard. This new, national-scale operational modeling system provides communities with timely, local-scale (0.5 km resolution) coastal hazard warnings for all wave conditions, allowing for rapid decision-making related to safety and emergency management. The modeling system also enables continued research into wave-driven processes at a broad range of coastal areas.