The highly urbanized estuary of San Francisco Bay is an excellent example of a location susceptible to flooding from both coastal and fluvial influences. As part of developing a forecast model that integrates fluvial and oceanic drivers, a case study of the Napa River and its interactions with the San Francisco Bay was performed. For this application we utilize Delft3D-FM, a hydrodynamic model that computes conservation of mass and momentum on a flexible mesh grid, to calculate water levels that account for tidal forcing, storm surge generated by wind and pressure fields, and river flows. We simulated storms with realistic atmospheric pressure, river discharge, and tidal forcing to represent a realistic joint fluvial and coastal storm event. Storm conditions were applied to both a realistic field-scale Napa river drainage as well as an idealized geometry. With these scenarios, we determine how the extent, level, and duration of flooding is dependent on these atmospheric and hydrologic parameters. Unsurprisingly, the model indicates that maximal water levels will occur in a tidal river when high tides, storm surge, and large fluvial discharge events are coincident. Model results also show that large tidal amplitudes diminish storm surge propagation upstream and that phasing between peak fluvial discharges and high tide is important for predicting when and where the highest water levels will occur. The interactions between tides, river discharge, and storm surge are not simple, indicating the need for more integrated flood forecasting models in the future.

Sea level rise (SLR) is placing both immediate and long‐term pressures on coastal communities to take protective actions. Projects in the United States, and in many locations throughout the world, generally involve local jurisdictions raising the elevation of shoreline protection elements, with limited or no analysis of the feedback between shoreline management decisions and the impacts to water levels regionally. Our study examines the impact of local shoreline development on regional flood risk and considers SLR scenarios up to 1.5 m using a large‐scale numerical model, as an example, for San Francisco Bay. Here we show that measures to prevent flooding along an embayment shoreline in one location or subregion may increase inundation elsewhere in the system. The network of interactions occurs not only within subbasins of the Bay but also across the greater geographic extent from one end of the Bay to the other, and local jurisdiction may have either reciprocal relationships with or asymmetric impacts on one other. Importantly, the nature of the interaction network is seen to evolve with SLR: interactions are purely subregional at current sea level but with higher sea level (e.g., 1 m of SLR), not only do the subregional interdependencies strengthen but also regional interdependences emerge.

Sea level rise increases the risk of storms and other short‐term water‐rise events, because it sets a higher water level such that coastal surges become more likely to overtop protections and cause floods. To protect coastal communities, it is necessary to understand the interaction among multiday and tidal sea level variabilities, coastal infrastructure, and sea level rise. We performed a series of numerical simulations for San Francisco Bay to examine two shoreline scenarios and a series of short‐term and long‐term sea level variations. The two shoreline configurations include the existing topography and a coherent full‐bay containment that follows the existing land boundary with an impermeable wall. The sea level variability consists of a half‐meter perturbation, with duration ranging from 2 days to permanent (i.e., sea level rise). The extent of coastal flooding was found to increase with the duration of the high‐water‐level event. The nonlinear interaction between these intermediate scale events and astronomical tidal forcing only contributes ∼1% of the tidal heights; at the same time, the tides are found to be a dominant factor in establishing the evolution and diffusion of multiday high water events. Establishing containment at existing shorelines can change the tidal height spectrum up to 5%, and the impact of this shoreline structure appears stronger in the low‐frequency range. To interpret the spatial and temporal variability at a wide range of frequencies, Optimal Dynamic Mode Decomposition is introduced to analyze the coastal processes and an inverse method is applied to determine the coefficients of a 1‐D diffusion wave model that quantify the impact of bottom roughness, tidal basin geometry, and shoreline configuration on the high water events.

Accurate and timely flood forecasts are critical for making emergency-response decisions regarding public safety, infrastructure operations, and resource allocation. One of the main challenges for coastal flood forecasting systems is a lack of reliable forecast data of large-scale oceanic and watershed processes and the combined effects of multiple hazards, such as compound flooding at river mouths. Offshore water level anomalies, known as remote Non-Tidal Residuals (NTRs), are caused by processes such as downwelling, offshore wind setup, and also driven by ocean-basin salinity and temperature changes, common along the west coast during El Niño events. Similarly, fluvial discharges can contribute to extreme water levels in the coastal area, while they are dominated by large-scale watershed hydraulics. However, with the recent emergence of reliable large-scale forecast systems, coastal models now import the essential input data to forecast extreme water levels in the nearshore. Accordingly, we have developed Hydro-CoSMoS, a new coastal forecast model based on the USGS Coastal Storm Modeling System (CoSMoS) powered by the Delft3D San Francisco Bay and Delta community model. In this work, we studied the role of fluvial discharges and remote NTRs on extreme water levels during a February 2019 storm by using Hydro-CoSMoS in hindcast mode. We simulated the storm with and without real-time fluvial discharge data to study their effect on coastal water levels and flooding extent, and highlight the importance of watershed forecast systems such as NOAA’s National Water Model (NWM). We also studied the effect of remote NTRs on coastal water levels in San Francisco Bay during the 2019 February storm by utilizing the data from a global ocean model (HYCOM). Our results showed that accurate forecasts of remote NTRs and fluvial discharges can play a significant role in predicting extreme water levels in San Francisco Bay. This pilot application in San Francisco Bay can serve as a basis for integrated coastal flood modeling systems in complex coastal settings worldwide.

A wave‐driven surface‐buoyant jet exiting a coral reef was studied in order to quantify the amount of water reentrained over the reef crest. Both moored observations and Lagrangian drifters were used to study the fate of the buoyant jet. To investigate in detail the effects of buoyancy and alongshore flow variations, we developed an idealized numerical model of the system. Consistent with previous work, the ratio of alongshore velocity to jet velocity and the jet internal Froude number were found to be important determinants of the fate of the jet. In the absence of buoyancy, the entrainment of fluid at the reef crest creates a significant amount of retention, keeping 60% of water in the reef system. However, when the jet is lighter than the ambient ocean water, the net effect of buoyancy is to enhance the separation of the jet from shore, leading to a greater export of reef water. Matching observations, our modeling predicts that buoyancy limits retention to 30% of the jet flow for conditions existing on the Moorea reef. Overall, the combination of observations and modeling we present here shows that reef‐ocean temperature gradients can play an important role in reef‐ocean exchanges.

Compound flooding, resulting from a combination of riverine and coastal processes, is a complex but important hazard to resolve along urbanized shorelines in the vicinity of river mouths. However, inland flooding models rarely consider oceanographic conditions, and vice versa for coastal flood models. Here, we describe the development of an operational, integrated coastal-watershed flooding model to address this issue of compound flooding in a highly urbanized estuarine environment, San Francisco Bay (CA, USA), where the surrounding communities are susceptible to flooding along the bay shoreline and inland rivers and creeks that drain to the bay. The integrated tributary-coastal forecast model (Hydro-Coastal Storm Modeling System, or Hydro-CoSMoS) was developed to provide water managers and other users with flood forecast information beyond what is currently available. Results presented here are focused on the interaction of the Napa River watershed and the San Pablo Bay at the northern end of San Francisco Bay. This paper describes the modeling setup, the scenario used in a tabletop exercise (TTE), and the assessment of the various flood forecast information products. Hydro-CoSMoS successfully demonstrated the capability to provide watershed and coastal flood information at scales and locations where no such information is currently available and was also successful in showing how tributary flows could be used to inform the coastal storm model during a flooding scenario. The TTE provided valuable feedback on how to guide continued model development and to inform what model outputs and formats are most useful to end-users.