Tidal response to sea‐level rise (SLR) varies in different coastal systems. To provide a generic pattern of tidal response to SLR, a systematic investigation was conducted using numerical techniques applied to idealized and realistic estuaries, with model results cross‐checked by analytical solutions. Our results reveal that the response of tidal range to SLR is nonlinear, spatially heterogeneous, and highly affected by the length and bathymetry of an estuary and weakly affected by the estuary convergence with an exception of strong convergence. Contrary to the common assumption that SLR leads to a weakened bottom friction, resulting in increased tidal amplitude, we demonstrate that tidal range is likely to decrease in short estuaries and in estuaries with a narrow channel and large low‐lying shallow areas.

Hurricane Sandy inflicted heavy damage in New York City and the New Jersey coast as the second costliest storm in history. A large-scale, unstructured grid storm tide model, Semi-implicit Eulerian Lagrangian Finite Element (SELFE), was used to hindcast water level variation during Hurricane Sandy in the mid-Atlantic portion of the U.S. East Coast. The model was forced by eight tidal constituents at the model’s open boundary, 1500 km away from the coast, and the wind and pressure fields from atmospheric model Regional Atmospheric Modeling System (RAMS) provided by Weatherflow Inc. The comparisons of the modeled storm tide with the NOAA gauge stations from Montauk, NY, Long Island Sound, encompassing New York Harbor, Atlantic City, NJ, to Duck, NC, were in good agreement, with an overall root mean square error and relative error in the order of 15–20 cm and 5%–7%, respectively. Furthermore, using large-scale model outputs as the boundary conditions, a separate sub-grid model that incorporates LIDAR data for the major portion of the New York City was also set up to investigate the detailed inundation process. The model results compared favorably with USGS’ Hurricane Sandy Mapper database in terms of its timing, local inundation area, and the depth of the flooding water. The street-level inundation with water bypassing the city building was created and the maximum extent of horizontal inundation was calculated, which was within 30 m of the data-derived estimate by USGS.

We study the compound flooding processes that occurred in Hurricane Florence (2018), which was accompanied by heavy precipitation, using a 3D creek-to-ocean hydrodynamic model. We examine the important role played by barrier islands in the observed compound surges in the coastal watershed. Locally very high resolution is used in some watershed areas in order to resolve small features that turn out to be critical for capturing the observed high water marks locally. The wave effects are found to be significant near barrier islands and have contributed to some observed over-toppings and breaches. Results from sensitivity tests applying each of the three major forcing factors (oceanic, fluvial, and pluvial) separately are succinctly summarized in a “dominance map” that highlights significant compound effects in most of the affected coastal watersheds, estuaries, and back bays behind the barrier islands. Operational forecasts based on the current model are being set up at NOAA to help coastal resource and emergency managers with disaster planning and mitigation efforts.