Journal cover
Journal topic
**Natural Hazards and Earth System Sciences**
An interactive open-access journal of the European Geosciences Union

Journal topic

- About
- Editorial board
- Articles
- Special issues
- Highlight articles
- Manuscript tracking
- Book reviews
- Subscribe to alerts
- Peer review
- For authors
- For reviewers
- EGU publications
- Imprint
- Data protection

- About
- Editorial board
- Articles
- Special issues
- Highlight articles
- Manuscript tracking
- Book reviews
- Subscribe to alerts
- Peer review
- For authors
- For reviewers
- EGU publications
- Imprint
- Data protection

**Research article**
20 Aug 2019

**Research article** | 20 Aug 2019

Spatial distribution of water level impacting back-barrier bays

- U.S. Geological Survey, Woods Hole, Massachusetts 02543, USA

- U.S. Geological Survey, Woods Hole, Massachusetts 02543, USA

**Correspondence**: Alfredo L. Aretxabaleta (aaretxabaleta@usgs.gov)

**Correspondence**: Alfredo L. Aretxabaleta (aaretxabaleta@usgs.gov)

Abstract

Back to toptop
Water level in semi-enclosed bays, landward of barrier islands, is mainly driven by offshore sea level fluctuations that are modulated by bay geometry and bathymetry, causing spatial variability in the ensuing response (transfer). Local wind setup can have a complementary role that depends on wind speed, fetch, and relative orientation of the wind direction and the bay. Bay area and inlet geometry and bathymetry primarily regulate the magnitude of the transfer between open ocean and bay. Tides and short-period offshore oscillations are more damped in the bays than longer-lasting offshore fluctuations, such as a storm surge and sea level rise. We compare observed and modeled water levels at stations in a mid-Atlantic bay (Barnegat Bay) with offshore water level proxies. Observed water levels in Barnegat Bay are compared and combined with model results from the Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) modeling system to evaluate the spatial structure of the water level transfer. Analytical models based on the dimensional characteristics of the bay are used to combine the observed data and the numerical model results in a physically consistent approach. Model water level transfers match observed values at locations inside the bay in the storm frequency band (transfers ranging from 50 %–100 %) and tidal frequencies (10 %–55 %). The contribution of frequency-dependent local setup caused by wind acting along the bay is also considered. The wind setup effect can be comparable in magnitude to the offshore transfer forcing during intense storms. The approach provides transfer estimates for locations inside the bay where observations were not available, resulting in a complete spatial characterization. An extension of the methodology that takes advantage of the ADCIRC tidal database for the east coast of the United States allows for the expansion of the approach to other bay systems. Detailed spatial estimates of water level transfer can inform decisions on inlet management and contribute to the assessment of current and future flooding hazard in back-barrier bays and along mainland shorelines.

Download & links

How to cite

Back to top
top
How to cite.

Aretxabaleta, A. L., Ganju, N. K., Defne, Z., and Signell, R. P.: Spatial distribution of water level impacting back-barrier bays, Nat. Hazards Earth Syst. Sci., 19, 1823–1838, https://doi.org/10.5194/nhess-19-1823-2019, 2019.

1 Introduction

Back to toptop
Back-barrier bays or coastal lagoons are common features along the coast of the United States. Their depths are usually on the order of a few meters and their horizontal extents are on the order of several tens of kilometers. They are often surrounded by highly populated areas and susceptible to intense human and environmental stressors. During storms, a surge and larger-than-normal waves combine to inundate low-elevation areas, resulting in hazards to coastal communities. Both hurricanes and winter storms affect coastal populations, infrastructure, and natural resources along the coastal bays of the United States (Nicholls et al., 2007, 2014; Rahmstorf, 2017; Wahl et al., 2017).

Hazard assessments consist of a characterization of the spatial and temporal extent of damaging physical events and the determination of the specific characteristics of those events (Ludwig et al., 2018). While flooding on the mainland side of back-barrier bays has severe socioeconomic implications, most coastal hazard evaluations (Gornitz et al., 1994; Thieler and Hammar-Klose, 1999; Klein and Nicholls, 1999; Kunreuther et al., 2000; Neumann et al., 2015; Vitousek et al., 2017) have focused on open-coast areas. Vulnerability evaluation of coastal areas around back-barrier bays requires extensive knowledge of the main hazard sources and their physical controls.

Water level in the bays is mainly driven by offshore sea level fluctuations with additional effects from local wind and wave setups. The bay exchange with the ocean usually occurs through narrow inlets. The size of the inlet determines the frictional effects and the amount of dampening offshore fluctuations encounter (Keulegan, 1967). Tides and short-period offshore oscillations tend to be more dampened in the bays than longer-lasting offshore fluctuations, such as a storm surge and sea level rise.

Bay water level fluctuations are linked to offshore forcing, especially at low frequencies, while wind acting directly over the bay is more connected to current fluctuations in the bay (Garvine, 1985). Chuang and Swenson (1981) determined that water level changes at subtidal frequencies in Lake Pontchartrain were controlled by coupled coastal ocean–bay fluctuations. Wong and Wilson (1984) studied subtidal sea level fluctuations in Great South Bay and again found them primarily driven by bay–shelf coupling. In Delaware Bay, a bay–inlet system with a relatively large opening, Wong and DiLorenzo (1988) showed that remote effects dominate over local effects and that fluctuations at both tidal and subtidal frequencies in connected bays of the Delaware Bay system were forced by shelf sea level.

More recently, Aretxabaleta et al. (2014) analyzed water level data in Barnegat Bay and Great South Bay before and after Hurricane Sandy and demonstrated that the offshore–bay transfer was not significantly altered by the geomorphologic changes caused by the storm. Aretxabaleta et al. (2017) described observed changes in both tidal amplitude and bay water level transfer from offshore in Great South Bay and connected bays and related the changes to the dredging of nearby inlets and the changing size of a breach across Fire Island caused by Hurricane Sandy. They also introduced an analytical model, based on the Chuang and Swenson (1981) approach but extended to interconnected bays, that incorporated bay and inlet dimensions and matched the observed transfer of offshore water level fluctuations into the bay system.

In this study, we combine an analysis of observed water levels in Barnegat Bay with the results of numerical models and an analytical description of the system to characterize the spatial characteristics of the bay response to offshore fluctuations. The observations provide detailed information at five locations in the bay, while the numerical simulations can expand the analysis to the entire bay system. The analytical model allows for the evaluation of the importance of the dominant factors affecting water level in bays. The combined approach can be used to provide consistent spatial maps of offshore water level impact into back-barrier bays. The method will be useful for coastal hazard assessment, assisting in the management of nuisance flooding (Moftakhari et al., 2018), providing spatial differences in vulnerability to perigean spring tides (king tides), and planning for flooding in response to storms of different durations. The final hazard estimates will be included as part of the U.S. Geological Survey Coastal Change Hazard portal (USGS, 2018) in an effort to expand the total water level predictions (Aretxabaleta et al., 2019).

2 Regional description

Back to toptop
The Barnegat Bay–Little Egg Harbor (BBLEH) estuary is a back-barrier bay along the coast of New Jersey (Fig. 1). It is a shallow (average depth around 1.5 m) bay connected to the ocean through three openings: Little Egg Inlet in the south, Barnegat Inlet in the center, and the Point Pleasant Canal, which is a much smaller connection in the north of the bay. Offshore tidal amplitudes decrease slightly from 0.7 m in northern New York Bay to 0.6 m in central New Jersey. The southern sub-embayment (Little Egg Harbor) is more connected to the open ocean, with tidal amplitudes ranging between 0.2 and 0.5 m, while its northern part (Barnegat Bay) has less exchange and tidal amplitudes are smaller than 0.2 m (Chant, 2001; Defne and Ganju, 2015).

3 Observational and model data

Back to toptop
Water level observations from five stations in the BBLEH system (Table 1) and from two external coastal stations are used to determine transfer from ocean to bay. The bay stations started recording in October 2007, while Sandy Hook and Atlantic City are long-term NOAA water level stations, operational since 1910 and 1911, respectively (Table 1). Wind observations were obtained from the National Data Buoy Center (NDBC) buoy 44065 (New York Harbor Entrance) for the period 2008–2018.

We used numerical simulations of Barnegat Bay for the period March–September 2012 (Defne and Ganju, 2015; data available from Defne and Ganju, 2018) and October–December 2012 (USGS, 2019) to obtain the spatial character of the water level response. The simulations used the Coupled Ocean–Atmosphere–Wave–Sediment Transport modeling system (COAWST; Warner et al., 2010). The model resolution ranged from 40 to 200 m, with the higher resolution located near complex geometry and around the inlets. The model is forced at the boundaries with tides from the ADCIRC tidal database for the western North Atlantic Ocean (Szpilka et al., 2016) and open-ocean forcing from subtidal water level and velocity from the ESPreSSO model (Wilkin and Hunter, 2013; http://www.myroms.org/espresso/, last access: 15 August 2019) and COAWST US east coast forecast (Warner et al., 2010; https://catalog.data.gov/dataset/coawst-forecast-system-usgs-us-east-coast-and-gulf-of, last access: 15 August 2019). Defne and Ganju (2015) showed that the numerical model solution had sufficient flow and elevation skill to characterize bay dynamics under normal and storm conditions.

The impact in the bay of offshore forcing can be evaluated spectrally by estimating transfer functions in frequency space between observed water levels offshore (input) and in the bays (output). The transfer functions are calculated using a Hanning window with overlapping (50 %) data segments with lengths of 29 d to provide estimates near the main tidal frequencies (Aretxabaleta et al., 2017). The uncertainty envelopes for the transfer function are estimated using the Bendat and Piersol (1986) formulation.

4 Analytical water level models

Back to toptop
The impact of ocean sea level fluctuations in the bay can be explored with an analytical model of a generic bay system (Fig. 2), consisting of multiple interconnected sub-embayments connected to the offshore by three separate inlets: Little Egg Inlet, Barnegat Inlet, and the Point Pleasant Canal. The model assumes that the bay water level responds as a level surface in each sub-embayment to ocean fluctuations, as local forcing in the bay is not included. The formulation is an extension of the approach proposed by Chuang and Swenson (1981) for a single inlet connecting to a bay and expanded by Wong and DiLorenzo (1988) to two connected bays and to multiple bays and inlets by Aretxabaleta et al. (2017). An analytical solution can be found for the entire system, with expressions for all the connections in the system. The model solves the along-channel depth-averaged momentum equation based on the balance between frictional effects and the elevation gradient between the offshore and bay areas and the continuity equation for the bay–channel system based on the changing volume of the bays as water flows through the inlets. The model also allows the estimation of the effect of the breach in Mantoloking during Hurricane Sandy. An analytical solution can be found by dividing the entire system into five sub-embayments (based on constrictions inside the bay system), resulting in a system of equations that includes 13 equations and unknowns (Appendix A).

$$\begin{array}{}\text{(1)}& \begin{array}{rl}{\displaystyle \frac{\partial}{\partial t}}\left(\begin{array}{l}{u}_{\mathrm{1}}\\ {u}_{\mathrm{2}}\\ {u}_{\mathrm{3}}\\ {u}_{\mathrm{4}}\\ {u}_{\mathrm{5}}\\ {u}_{\mathrm{6}}\\ {u}_{\mathrm{7}}\\ {u}_{\mathrm{8}}\end{array}\right)& =g\left(\begin{array}{l}\mathrm{1}/{L}_{\mathrm{1}}\\ \mathrm{1}/{L}_{\mathrm{2}}\\ \mathrm{1}/{L}_{\mathrm{3}}\\ \mathrm{1}/{L}_{\mathrm{4}}\\ \mathrm{1}/{L}_{\mathrm{5}}\\ \mathrm{1}/{L}_{\mathrm{6}}\\ \mathrm{1}/{L}_{\mathrm{7}}\\ \mathrm{1}/{L}_{\mathrm{8}}\end{array}\right)\left(\begin{array}{l}{\mathit{\eta}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{LEI}}-{\mathit{\eta}}_{\mathrm{1}}\\ {\mathit{\eta}}_{\mathrm{1}}-{\mathit{\eta}}_{\mathrm{2}}\\ {\mathit{\eta}}_{\mathrm{2}}-{\mathit{\eta}}_{\mathrm{3}}\\ {\mathit{\eta}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{BI}}-{\mathit{\eta}}_{\mathrm{3}}\\ {\mathit{\eta}}_{\mathrm{3}}-{\mathit{\eta}}_{\mathrm{4}}\\ {\mathit{\eta}}_{\mathrm{4}}-{\mathit{\eta}}_{\mathrm{5}}\\ {\mathit{\eta}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{breach}}-{\mathit{\eta}}_{\mathrm{5}}\\ {\mathit{\eta}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{PPC}}-{\mathit{\eta}}_{\mathrm{5}}\end{array}\right)\\ & -\left(\begin{array}{l}{u}_{\mathrm{1}}{r}_{\mathrm{1}}/{h}_{\mathrm{1}}\\ {u}_{\mathrm{2}}{r}_{\mathrm{2}}/{h}_{\mathrm{2}}\\ {u}_{\mathrm{3}}{r}_{\mathrm{3}}/{h}_{\mathrm{3}}\\ {u}_{\mathrm{4}}{r}_{\mathrm{4}}/{h}_{\mathrm{4}}\\ {u}_{\mathrm{5}}{r}_{\mathrm{5}}/{h}_{\mathrm{5}}\\ {u}_{\mathrm{6}}{r}_{\mathrm{6}}/{h}_{\mathrm{6}}\\ {u}_{\mathrm{7}}{r}_{\mathrm{7}}/{h}_{\mathrm{7}}\\ {u}_{\mathrm{8}}{r}_{\mathrm{8}}/{h}_{\mathrm{8}}\end{array}\right)\end{array}\end{array}$$

$$\begin{array}{}\text{(2)}& \begin{array}{rl}& \left(\begin{array}{l}{A}_{\mathrm{1}}\\ {A}_{\mathrm{2}}\\ {A}_{\mathrm{3}}\\ {A}_{\mathrm{4}}\\ {A}_{\mathrm{5}}\end{array}\right){\displaystyle \frac{\partial}{\partial t}}\left(\begin{array}{l}{\mathit{\eta}}_{\mathrm{1}}\\ {\mathit{\eta}}_{\mathrm{2}}\\ {\mathit{\eta}}_{\mathrm{3}}\\ {\mathit{\eta}}_{\mathrm{4}}\\ {\mathit{\eta}}_{\mathrm{5}}\end{array}\right)=\\ & \left(\begin{array}{l}{h}_{\mathrm{1}}{W}_{\mathrm{1}}{u}_{\mathrm{1}}-{h}_{\mathrm{2}}{W}_{\mathrm{2}}{u}_{\mathrm{2}}\\ {h}_{\mathrm{2}}{W}_{\mathrm{2}}{u}_{\mathrm{2}}-{h}_{\mathrm{3}}{W}_{\mathrm{3}}{u}_{\mathrm{3}}\\ {h}_{\mathrm{3}}{W}_{\mathrm{3}}{u}_{\mathrm{3}}+{h}_{\mathrm{4}}{W}_{\mathrm{4}}{u}_{\mathrm{4}}-{h}_{\mathrm{5}}{W}_{\mathrm{5}}{u}_{\mathrm{5}}\\ {h}_{\mathrm{5}}{W}_{\mathrm{5}}{u}_{\mathrm{5}}-{h}_{\mathrm{6}}{W}_{\mathrm{6}}{u}_{\mathrm{6}}\\ {h}_{\mathrm{6}}{W}_{\mathrm{6}}{u}_{\mathrm{6}}+{h}_{\mathrm{7}}{W}_{\mathrm{7}}{u}_{\mathrm{7}}+{h}_{\mathrm{8}}{W}_{\mathrm{8}}{u}_{\mathrm{8}}\end{array}\right)\end{array}\end{array}$$

Here *g* is the gravitational acceleration, *A*_{m} is the surface area of
sub-embayment *m*; *η*_{m} the sea level in the *m* sub-embayment; *η*_{o}
the sea level in the ocean; *h*_{n} the water depth; and *W*_{n} the width, *L*_{n} the length, and *r*_{n} the
linear drag coefficient of channel *n*. *ϕ*_{LEI}, *ϕ*_{BI}, *ϕ*_{breach}, and *ϕ*_{PPC} are the linear frequency-dependent relationships
between the water levels at offshore proxy stations (Sandy Hook or Atlantic
City) and the water level just offshore of Little Egg Inlet, Barnegat Inlet,
the breach at Mantoloking caused by Hurricane Sandy, and the Point Pleasant Canal.

Assuming $\mathit{\eta}=\stackrel{\mathrm{\u0303}}{\mathit{\eta}}{e}^{i\mathit{\omega}t}$ and $u=\stackrel{\mathrm{\u0303}}{u}{e}^{i\mathit{\omega}t}$,
where $\stackrel{\mathrm{\u0303}}{\mathit{\eta}}$ and $\stackrel{\mathrm{\u0303}}{u}$ represent the magnitude of the water
level and velocity oscillations, respectively, to reduce the size of
the equations, we can define ${K}_{n}=\frac{{h}_{n}{W}_{n}g}{{L}_{n}\left(\frac{{r}_{n}}{{h}_{n}}+i\mathit{\omega}\right)}$ for *n*=1, …, 8, as the
relative contribution of each sub-embayment based on its geometric and
frictional characteristics. Then, with the proper rearrangement, it yields the following equation.

$$\begin{array}{}\text{(3)}& {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}={\displaystyle \frac{{\displaystyle \begin{array}{c}{K}_{\mathrm{4}}\left({\mathit{\eta}}_{\mathrm{o}}\right){\mathit{\varphi}}_{\mathrm{BI}}+{\displaystyle \frac{\frac{{K}_{\mathrm{1}}{K}_{\mathrm{2}}{K}_{\mathrm{3}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{LEI}}}{i\mathit{\omega}{A}_{\mathrm{1}}+{K}_{\mathrm{1}}+{K}_{\mathrm{2}}}}{i\mathit{\omega}{A}_{\mathrm{2}}+{K}_{\mathrm{2}}+{K}_{\mathrm{3}}-\frac{{K}_{\mathrm{2}}{K}_{\mathrm{2}}}{i\mathit{\omega}{A}_{\mathrm{1}}+{K}_{\mathrm{1}}+{K}_{\mathrm{2}}}}}\\ +\frac{\frac{{K}_{\mathrm{5}}{K}_{\mathrm{6}}\left({K}_{\mathrm{7}}{\mathit{\varphi}}_{\mathrm{breach}}+{K}_{\mathrm{8}}{\mathit{\varphi}}_{\mathrm{PPC}}\right){\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}}{i\mathit{\omega}{A}_{\mathrm{5}}+{K}_{\mathrm{6}}+{K}_{\mathrm{7}}+{K}_{\mathrm{8}}}}{i\mathit{\omega}{A}_{\mathrm{4}}+{K}_{\mathrm{5}}+{K}_{\mathrm{6}}-\frac{{K}_{\mathrm{6}}{K}_{\mathrm{6}}}{i\mathit{\omega}{A}_{\mathrm{5}}+{K}_{\mathrm{6}}+{K}_{\mathrm{7}}+{K}_{\mathrm{8}}}}\end{array}}}{{\scriptscriptstyle \begin{array}{c}i\mathit{\omega}{A}_{\mathrm{3}}+{K}_{\mathrm{3}}+{K}_{\mathrm{4}}+{K}_{\mathrm{5}}-\frac{{K}_{\mathrm{3}}{K}_{\mathrm{3}}}{i\mathit{\omega}{A}_{\mathrm{2}}+{K}_{\mathrm{2}}+{K}_{\mathrm{3}}-\frac{{K}_{\mathrm{2}}{K}_{\mathrm{2}}}{i\mathit{\omega}{A}_{\mathrm{1}}+{K}_{\mathrm{1}}+{K}_{\mathrm{2}}}}\\ -{\scriptscriptstyle \frac{{K}_{\mathrm{5}}{K}_{\mathrm{5}}}{i\mathit{\omega}{A}_{\mathrm{4}}+{K}_{\mathrm{5}}+{K}_{\mathrm{6}}-\frac{{K}_{\mathrm{6}}{K}_{\mathrm{6}}}{i\mathit{\omega}{A}_{\mathrm{5}}+{K}_{\mathrm{6}}+{K}_{\mathrm{7}}+{K}_{\mathrm{8}}}}}\end{array}}}}\end{array}$$

The solution for the water level of the central sub-embayment can be used to recursively calculate the solutions for the rest of the sub-embayments as follows.

$$\begin{array}{}\text{(4)}& {\displaystyle}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}={\displaystyle \frac{{K}_{\mathrm{3}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}+\frac{{K}_{\mathrm{2}}{K}_{\mathrm{1}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{LEI}}}{i\mathit{\omega}{A}_{\mathrm{1}}+{K}_{\mathrm{1}}+{K}_{\mathrm{2}}}}{i\mathit{\omega}{A}_{\mathrm{2}}+{K}_{\mathrm{2}}+{K}_{\mathrm{3}}-\frac{{K}_{\mathrm{2}}{K}_{\mathrm{2}}}{i\mathit{\omega}{A}_{\mathrm{1}}+{K}_{\mathrm{1}}+{K}_{\mathrm{2}}}}}\text{(5)}& {\displaystyle}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}={\displaystyle \frac{{K}_{\mathrm{5}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}+\frac{{K}_{\mathrm{6}}\left({K}_{\mathrm{7}}{\mathit{\varphi}}_{\mathrm{breach}}+{K}_{\mathrm{8}}{\mathit{\varphi}}_{\mathrm{PPC}}\right){\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}}{i\mathit{\omega}{A}_{\mathrm{5}}+{K}_{\mathrm{6}}+{K}_{\mathrm{7}}+{K}_{\mathrm{8}}}}{i\mathit{\omega}{A}_{\mathrm{4}}+{K}_{\mathrm{5}}+{K}_{\mathrm{6}}-\frac{{K}_{\mathrm{6}}{K}_{\mathrm{6}}}{i\mathit{\omega}{A}_{\mathrm{5}}+{K}_{\mathrm{6}}+{K}_{\mathrm{7}}+{K}_{\mathrm{8}}}}}\text{(6)}& {\displaystyle}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}={\displaystyle \frac{{K}_{\mathrm{1}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{LEI}}+{K}_{\mathrm{2}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}}{i\mathit{\omega}{A}_{\mathrm{1}}+{K}_{\mathrm{1}}+{K}_{\mathrm{2}}}}\text{(7)}& {\displaystyle}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}={\displaystyle \frac{\left({K}_{\mathrm{7}}{\mathit{\varphi}}_{\mathrm{breach}}+{K}_{\mathrm{8}}{\mathit{\varphi}}_{\mathrm{PPC}}\right){\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}+{K}_{\mathrm{6}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}}{i\mathit{\omega}{A}_{\mathrm{5}}+{K}_{\mathrm{6}}+{K}_{\mathrm{7}}+{K}_{\mathrm{8}}}}\end{array}$$

The resulting expressions include all the sub-embayment and offshore exchanges under the same assumptions of the Chuang and Swenson (1981) model (e.g., no local influences, no overtopping).

The contribution of local wind setup to the spatial distribution of water level inside the bay can be approximated following Wong and Moses-Hall (1998). The bay can be assumed to be a simple, long, and well-mixed embayment for which the cross-bay gradients and vertical stratification can be ignored. The linearized vertically integrated momentum and mass conservation equations are as follows:

$$\begin{array}{}\text{(8)}& {\displaystyle \frac{\partial U}{\partial t}}=-gh{\displaystyle \frac{\partial \mathit{\eta}}{\partial x}}+{\displaystyle \frac{\mathrm{1}}{{\mathit{\rho}}_{\mathrm{0}}}}\left({\mathit{\tau}}_{\mathrm{s}}-{\mathit{\tau}}_{\mathrm{b}}\right)=-gh{\displaystyle \frac{\partial \mathit{\eta}}{\partial x}}+{\mathit{\tau}}_{w}-r{\displaystyle \frac{U}{h}},\text{(9)}& {\displaystyle}\mathrm{and}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{\displaystyle \frac{\partial U}{\partial x}}=-{\displaystyle \frac{\partial \mathit{\eta}}{\partial t}},\end{array}$$

where *η* is the water level in the bay, *U* is the depth-integrated
along-bay velocity, *h* is the water depth, *τ*_{s} and *τ*_{b} are the
surface and bottom dynamic stresses, respectively, and *ρ*_{0} is the
water density. ${\mathit{\tau}}_{w}={\mathit{\tau}}_{\mathrm{s}}/{\mathit{\rho}}_{\mathrm{0}}$ is
the spatially invariant kinematic wind stress and *r* is linearized bottom
friction.

Under the assumption of $\mathit{\eta}=\stackrel{\mathrm{\u0303}}{\mathit{\eta}}{e}^{i\mathit{\omega}t}$and
$u=\stackrel{\mathrm{\u0303}}{u}{e}^{i\mathit{\omega}t}$, where *ω* is the cyclic frequency, the
resulting equation is as follows:

$$\begin{array}{}\text{(10)}& {\displaystyle \frac{{\partial}^{\mathrm{2}}\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}{\partial {x}^{\mathrm{2}}}}+\stackrel{\mathrm{\u0303}}{\mathit{\eta}}\left({\displaystyle \frac{{\mathit{\omega}}^{\mathrm{2}}}{gh}}-{\displaystyle \frac{i\mathit{\omega}r}{g{h}^{\mathrm{2}}}}\right)={\displaystyle \frac{{\partial}^{\mathrm{2}}\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}{\partial {x}^{\mathrm{2}}}}+\stackrel{\mathrm{\u0303}}{\mathit{\eta}}{k}^{\mathrm{2}}=\mathrm{0},\end{array}$$

with boundary conditions $\mathit{\eta}\left(x=\mathrm{0},\mathit{\omega}\right)=\mathrm{0}$ assuming no
offshore forcing at the entrance (this assumption will be revisited in the
next section) and $\frac{\partial \stackrel{\mathrm{\u0303}}{\mathit{\eta}}\left(x=L,\mathit{\omega}\right)}{\partial x}=\frac{\stackrel{\mathrm{\u0303}}{{\mathit{\tau}}_{w}}\left(\mathit{\omega}\right)}{gh}$
assuming no flux at the head (*x*=*L*). $\stackrel{\mathrm{\u0303}}{{\mathit{\tau}}_{w}}\left(\mathit{\omega}\right)$ represents the magnitude of the kinematic wind stress that
results in water level fluctuations at a specific frequency.

The solution is as follows:

$$\begin{array}{}\text{(11)}& {\displaystyle}\stackrel{\mathrm{\u0303}}{\mathit{\eta}}\left(\mathit{\omega}\right)={\displaystyle \frac{\stackrel{\mathrm{\u0303}}{{\mathit{\tau}}_{w}}\left(\mathit{\omega}\right)\mathrm{sin}\left(kx\right)}{ghk\phantom{\rule{0.125em}{0ex}}\mathrm{cos}\left(kL\right)}},\text{(12)}& {\displaystyle}\mathrm{where}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}k={\left({\displaystyle \frac{{\mathit{\omega}}^{\mathrm{2}}}{gh}}-{\displaystyle \frac{i\mathit{\omega}r}{g{h}^{\mathrm{2}}}}\right)}^{\mathrm{1}/\mathrm{2}}.\end{array}$$

The wavenumber *k* determines the spatial response of the transfer between the
wind stress and the bay water level. The imaginary wavenumber part leads to
exponential decay based on the frictional characteristics. The magnitude of
water level is obtained from Eq. (11), while the ratio of real to
imaginary parts provides information about the phase lag between wind stress
and water level.

The local and remote effects can be combined in following the approach by
Wong and Moses-Hall (1998). While bays can exhibit complex spatial responses
to wind forcing especially in terms of currents (Csanady, 1973; Hunter and
Hearn, 1987; Cioffi et al., 2005), the basic response can be summarized as
the sum of local (wind) and remote (surge) forcings. The boundary condition
for the local wind effect can be altered to account for the influence of
the offshore water level, *η*_{o}. The resulting model is a modification of the
wind effect model that considers the analytical offshore influence in
Sect. 4.1. In a system with a single inlet, the solution can be simply
stated, as in Wong and Moses-Hall (1998):

$$\begin{array}{}\text{(13)}& \stackrel{\mathrm{\u0303}}{\mathit{\eta}}\left(\mathit{\omega}\right)={\displaystyle \frac{\stackrel{\mathrm{\u0303}}{{\mathit{\tau}}_{w}}\left(\mathit{\omega}\right)\phantom{\rule{0.125em}{0ex}}\mathrm{sin}\left(kx\right)}{ghk\phantom{\rule{0.125em}{0ex}}\mathrm{cos}\left(kL\right)}}+{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}\left(\mathit{\omega}\right){\displaystyle \frac{\mathrm{cos}\left(k\left(L-x\right)\right)}{\mathrm{cos}\left(kL\right)}}.\end{array}$$

In a system with multiple connections with the offshore, the solution can be
more complex. One limitation of the approach is that it utilizes a linear
friction approximation. To produce a better approximation that takes into
account the complex frictional conditions of the bay (e.g., varying
geometry, diverse bottom conditions, enhanced attenuation over submerged
vegetation), we can take a numerical solution of the bay that resolves the
tidal and subtidal water level conditions under realistic friction and
adjust the spatial distribution of the transfer from offshore accordingly.
As most of the water level variability in the bay is associated with the
*M*_{2} semidiurnal tidal constituent (Fig. 3) and the distribution of the
tide has been properly validated in the numerical simulations (Defne and
Ganju, 2015), we can take the spatial distribution of the *M*_{2} tidal
amplitude as a proxy for the internal frictional effects in the bay. Bottom
friction caused by both wind-driven and tidal effects is considered in the
numerical simulations. By adjusting the water level based on the numerical
*M*_{2} spatial distribution, we are approximating the complete frictional
characteristics of the bay. The adjustment is applied to each of the
sub-embayments using the following expression:

$$\begin{array}{}\text{(14)}& {\displaystyle \frac{{\widehat{\mathit{\eta}}}_{j}}{{\mathit{\eta}}_{\mathrm{o}}}}\left(x,\mathit{\omega}\right)=\mathrm{1}-{\displaystyle \frac{\left(\mathrm{1}-\frac{{\mathit{\eta}}_{j}}{{\mathit{\eta}}_{\mathrm{o}}}\left(\mathit{\omega}\right)\right)\left(\mathrm{1}-\frac{\mathit{\eta}\left(x,\phantom{\rule{0.125em}{0ex}}{\mathit{\omega}}_{{M}_{\mathrm{2}}}\right)}{\mathit{\eta}\left(\mathrm{offshore},{\mathit{\omega}}_{{M}_{\mathrm{2}}}\right)}\right)}{\left(\mathrm{1}-\frac{{\mathit{\eta}}_{j}}{{\mathit{\eta}}_{\mathrm{o}}}\left({\mathit{\omega}}_{{M}_{\mathrm{2}}}\right)\right)}},\end{array}$$

where *η*_{j}∕*η*_{o}(ω) is the transfer coefficient of the
sub-embayment *i* (single value) at frequency *ω*, $\mathit{\eta}\left(x,{\mathit{\omega}}_{{M}_{\mathrm{2}}}\right)$ is the amplitude of the *M*_{2} tidal fluctuations from
the numerical model solution (spatially variable), and
${\widehat{\mathit{\eta}}}_{j}/{\mathit{\eta}}_{\mathrm{o}}\left(x,\mathit{\omega}\right)$ is the spatially variable adjusted transfer
coefficient for sub-embayment *j*. The resulting adjusted transfer coefficients
provide estimates of the spatial changes not only between adjacent
sub-embayments but also inside each of the sub-embayments. The local wind
effects on bay water level can be added to the impact from offshore
fluctuations to obtain a combined local and remote water level response
estimate. In cases with simultaneous presence of wind and offshore level
fluctuations, the system can respond in a weakly nonlinear manner and
departures from the presented basic addition of the process are expected.

5 Results

Back to toptop
The maximum energy in water level spectra (Fig. 3, Table 2) was associated
with the *M*_{2} semidiurnal tidal constituent for the offshore proxies (SH
and AC) and at the stations TUC and ETH in the southern part of the BBLEH
area. For the locations in Barnegat Bay (WAR, SEH, MAN), maximum energy was
in the low-frequency band. Large spectral energy also occurred in the other
semidiurnal tidal frequencies (*S*_{2} and *N*_{2}), the diurnal frequencies
(*O*_{1} and *K*_{1}), the storm band (periods between 2 and 5 d), and
the low-frequency band (Table 2). The energy in the remaining bands
exhibited average fluctuations less than 0.03 m in size offshore, while in
the bay fluctuations were less than 0.01 m.

Transfer functions between Atlantic City (AC) and the five stations inside the bay (Fig. 4) for the longest available length of record showed a north to south gradient. The transfer of the offshore fluctuations was 50 %–80 % at periods between 2 and 5 d (storm band), except at Tuckerton (TUC; over 95 %). The transfers at diurnal periods were about 35 % for the three Barnegat Bay stations (WAR, SEH, MAN), about 45 % in Little Egg Harbor (ETH), and 80 % in Great Bay (TUC). For frequencies associated with the semidiurnal tides, the transfers were even more attenuated, with values of about 15 % (between 14 % and 16 %) inside Barnegat, between 30 % and 35 % at ETH, and between 60 % and 70 % at TUC. As the numerical model solution was only available for the period March–December 2012, the long-term (2007–2018) transfers were compared with shorter-term observations. The transfers were similar within the uncertainty envelopes for each station (not shown) for both datasets at most frequencies, except at Mantoloking (MAN), which showed enhanced transfers for periods between 1 and 5 d in the 2012 record, and at Seaside Heights (SEH), where transfers in the storm band were slightly attenuated during 2012. Transfer estimates using Sandy Hook (SH) as the offshore proxy instead of Atlantic City produced similar results (not shown). The transfer between stations AC and SH on the open coast (proxies for offshore fluctuations) has been shown to be close to one (Wong and Wilson, 1984; Aretxabaleta et al., 2014), confirming that the offshore forcing at all three inlets is about the same.

Transfers estimated from the numerical model solution (Fig. 5) showed similar magnitudes to the observed transfers within uncertainty envelopes provided by the Bendat and Piersol (1986) formulation at most frequencies. The observed and modeled transfer at diurnal and semidiurnal transfers were similar (within a few percentage points) at all stations, except that the model overestimated the semidiurnal transfer at TUC. Differences between modeled and observed estimates at MAN and SEH were only significant at frequencies that contained minimal energy. The model reproduced the enhanced transfer in the storm band at Mantoloking during 2012, suggesting a physical mechanism for the change that the model was able to capture but remains unexplained. The likely explanation is that the location of the Azores–Bermuda high-pressure system over the Atlantic in 2012 (Mattingly et al., 2015), associated with the negative phase of the North Atlantic Oscillation, resulted in average winds that lined up with the axis of the bay and caused enhanced wind setup in the northern part of the bay. The model overestimated the transfer at ETH in the storm band and underestimated the low-frequency transfer at Waretown. The likely cause for some of the discrepancies, especially at low frequencies, is the relatively short length of the available record.

The analytical model of offshore impact that considered five sub-embayments
(Sect. 4.1) was fit to the observed transfers to obtain an estimate of
linear friction. The fit considered the unevenly distributed energy spectra
(Fig. 3) with adjusted weight to the semidiurnal and low-frequency
components. The resulting friction was *r*=0.021 m s^{−1}. The associated
frictional adjustment time, ${t}_{\mathrm{fr}}=h/r$,
was about 1–5 min depending on the depth of the inlet. The analytical
curves (Fig. 6) matched the observed transfer function shape at most
frequencies. The analytical model with five sub-embayment domains captured
the north–south spatial differences in transfer. The analytical model for
the central Barnegat Bay sub-embayment (*A*_{3}) approximated the transfer
estimates from observations at Waretown at most frequencies (less than 5 % in the storm band). The analytical model for Great Bay (*A*_{1})
adequately matched observed transfers at Tuckerton at diurnal and
semidiurnal frequencies (less than 5 % difference) but underestimated the
transfer in the storm band (model estimates about 90 %, while observations
were above 95 %). Meanwhile, the analytical model for Little Egg Harbor
(*A*_{2}) matched the observed transfers at ETH within the uncertainty envelope,
except for a slight underprediction at diurnal frequencies (less than
5 %). The observed transfers at the northern stations (MAN and SEH) were
reproduced by the analytical model (*A*_{4}, *A*_{5}, respectively) at
diurnal and semidiurnal tides but were underpredicted for the higher storm
band frequencies (5 %–10 % less transfer in frequencies close to 2 d
periods) and overpredicted at low frequencies (about 10 % differences).
The analytical model was used to explore the effect on transfer of the
breach (*U*_{7}) at Mantoloking that opened during Hurricane Sandy. The
transfers were so minimally affected that the curves are indistinguishable,
with only a negligible enhancement (*<* 0.2 %) in transfer in the
northernmost sub-embayment (*A*_{5}). The breach was too small and shallow
for any significant volume transport to occur that would affect the large
bay.

The spectrum of the along-bay component (rotated 20^{∘}) of the wind
(Fig. 7a) from the offshore buoy NDBC 44065 (10 years, 2008–2018) showed
high energies in the storm band (2–5 d) and in low frequencies. The
largest single peak of energy was associated with 24 h period
oscillations, likely associated with sea breeze and matched energy values at
5 d frequencies. There was a small peak at inertial frequencies.

The local wind contribution to water level setup inside the bay was approximated using the Wong and Moses-Hall (1998) approach (Sect. 4.2). The resulting formulation showed the largest setup magnitudes near the head of the bay (e.g., northern part with wind blowing from the south) with a decay as distance from the head increased (Fig. 7b, c). The magnitude of the setup depended on the magnitude of the linear friction, with a smaller setup under stronger friction (Fig. 7b, c). The setup responded exponentially to fetch (distance), except over long durations and under low friction conditions, which were predominantly linear (Fig. 7b). The frictional control was less important at higher frequencies (Fig. 7c). As frequency increased there was less wind energy (Fig. 7a), so the frictional control is mostly important for low frequency and storm band wind fluctuations.

The resulting effect of the wind setup (or set-down) was small (less than
0.1 m with an along-bay wind stress of 0.1 Pa) for most of the domain
(Fig. 8). The estimate assumed a linear friction of the same magnitude as
in Sect. 5.1 (*r*=0.021 m s^{−1}). Under persistent wind stress of 0.1 Pa (about
8 m s^{−1} wind speed) in the along-bay direction, the resulting setups varied
depending on the frequency considered. Setup magnitudes over 0.2 m were
estimated for the 5 d period wind (Fig. 8c), while under half of that
magnitude was achieved for the 2 d persistent wind (Fig. 8b), and a much
smaller water level setup (peak smaller than 0.1 m) was estimated for the
sea breeze (Fig. 8a). During extreme events like Hurricane Sandy, under
intense wind stress, two additional effects should be considered: the depth
of the bay increases via the transfer of offshore surge, resulting in altered
setup response (Sect. 4.2), and the frictional effect is enhanced (a
larger linear friction would be needed) by the presence of wave-induced
roughness.

6 Discussion

Back to toptop
Following the approach described in Sect. 4.3, estimates of spatially
variable water level impact from offshore can be calculated (Fig. 9). The
*M*_{2} tidal constituent transfer (Fig. 9a) showed a large north to south
gradient, with values going from around 10 % in the north to over 80 % in
the vicinity of Little Egg Inlet. The role of Barnegat Inlet in enhancing
tidal transfer was greatly reduced, as most of the tide was attenuated in the
inlet. The contribution of the Point Pleasant Canal was also small as expected
from the tidal amplitudes (Chant, 2001; Defne and Ganju, 2015). The transfer
in the storm band of 2 d fluctuations (Fig. 9b) also showed a strong
north–south gradient, with values of about 50 % in Barnegat Bay, around
70 %–80 % in Little Egg Harbor, and larger values in Great Bay. The 5 d
offshore fluctuations were transferred more efficiently into the bay (Fig. 9c) with values over 70 % in the entire bay, reaching 80 %–90 % in Little
Egg Harbor, and over 90 % in Great Bay. Both storm band transfer estimates
were controlled by the exchange through Little Egg Inlet, with very localized
transfer enhancements in the vicinity of Barnegat Inlet and the Point Pleasant
Canal. While the presented estimates used Atlantic City as the offshore proxy,
similar results were obtained when Sandy Hook was used as the offshore
reference (as expected from Aretxabaleta et al., 2014).

When the magnitude of the fluctuations associated with a specific storm is
available (*η*_{o}), then an estimate of the average water level in the
bay during the storm can be obtained. For instance, for Hurricane Sandy the
offshore surge associated with the storm was of the order of 2–3 m.
Considering that the storm lasted for over a day, the water level transfer
would have been above 50 % in Barnegat Bay and above 70 % in Little Egg
Harbor. The resulting surge estimate in the bay was between 1 and 2 m just
considering the exchange through the existing inlets. There was reported
overtopping of the barrier island during the storm (McKenna et al., 2016)
that might have further increased water level in the bay that the proposed
method does not consider.

The wind setup effect inside the bay due to local wind can also be estimated for Hurricane Sandy using the approach in Sect. 4.2. Maximum wind stress during the storm was about 1 Pa. To obtain a maximum effect (worst-case scenario) it was assumed that the wind was persistently in the along-bay direction and that maximum stress was maintained for the duration of the storm. The maximum resulting water level considering the Wong and Moses-Hall method is linear with regard to wind stress magnitude (Fig. 7b) and would have been 10 times larger than the setup in Fig. 8b. The maximum wind setup would have been between 1 and 2 m, which was of the same order of magnitude as the surge produced from offshore sources. The cross-bay contribution to the wind setup during Hurricane Sandy was comparatively small, as wind direction was predominantly along-bay. Surge estimates from simple analytical formulations (State Committee for the Zuiderzee, 1926; Pugh, 1987) that do not consider storm duration produce similar magnitude results and are also dependent on the frictional response of the bay. Nonlinear interactions between local and remote effects may alter the total bay response, but these effects are likely second order.

The approach thus far was based on the combination of observations,
analytical models, and numerical models. In many systems, long-term
observations that allow for the estimation of transfer coefficients might
not be available. Also, numerical solutions for back-barrier bay systems tend
to be computationally expensive and might not be available for the period of
interest. We propose a relatively simple approach for some of these systems
based on the availability of high-resolution tidal solutions for the system.
The EC2015 ADCIRC tidal database (Szpilka et al., 2016) showed sufficient
resolution (down to 13 m in some areas) in many bays along the east coast of
the United States to resolve the tidal conditions with skill when compared
to NOAA CO-OPS stations and historic International Hydrographic Organization
(IHO) data. The EC2015 tidal database provides estimates for 37 tidal
constituents. Based on those constituents and assuming that the totality of
the offshore fluctuations at zero frequency reach the interior of the bay,
an estimate can be provided for the storm band frequencies. A weighted least-squares interpolation in the frequency domain was performed based on the
*M*_{4}, *K*_{2}, *S*_{2}, *L*_{2}, *M*_{2}, *N*_{2}, *K*_{1}, *P*_{1},
*O*_{1}, *Q*_{1} tidal amplitudes ratios between each point of the ADCIRC
domain inside the bay and a point in the offshore. Higher weight was given
to zero frequency to nudge toward 100 % transfer at zero frequency.
Estimates were calculated based on multiple locations inside the bay and
averaged to achieve a more robust calculation and also obtain an
approximation of the uncertainty associated with the estimate.

The resulting transfer estimates (Fig. 10) exhibited the same general
spatial patterns shown in the previous estimates (Fig. 9) with slight
differences. Some of the smaller features present in the COAWST numerical
solution (Defne and Ganju, 2015) were not present in the ADCIRC EC2015
domain. The *M*_{2} transfer estimate based on the tidal database (Fig. 10a) presented approximately the same magnitudes in most areas (average
difference less than 3 %). The 5 d transfer (Fig. 10c) was also
comparable to the solution described in Sect. 6.1 with values over 70 %
in the entire domain and the southern areas exceeding 90 % transfer. The
2 d transfer from ADCIRC (Fig. 10b) was 5 %–10 % higher than the direct
estimates (Fig. 9b). One of the benefits of the ADCIRC approach was the
possibility of providing an approximation of the uncertainty (Fig. 11).
The uncertainty estimate of the *M*_{2} transfer was about 1 %–2 % (Fig. 11a) with higher values in the southern part of the domain. The 2 d
transfer uncertainty (Fig. 11b) was above 4 % in Barnegat Bay in areas
of larger discrepancy between the ADCIRC and complete approaches. The
uncertainty estimates in the 5 d offshore water level transfer (Fig. 11c) in the northern part of the domain did not exceed 2.5 %.

The magnitude of the difference between the ADCIRC tidal database approach and the complete method highlighted in Sect. 4.3 was of the same order of magnitude or even smaller than the difference between observations and the analytical model (Fig. 6) or between observed transfers and numerical solutions (Fig. 5). This result emphasizes the validity of using the tidal database to calculate offshore transfer estimates, especially when water level observations inside the bay or numerical solutions are not available.

The effect of local wind setup will also need to be added to the ADCIRC-based estimate, especially during severe storms. The approach discussed in Sect. 5.2 or an even simpler surge calculation (e.g., from the steady state vertically averaged momentum equations, as in Pugh, 1987, from the traditional report of the State Committee for the Zuiderzee, 1926, or the updated frequency domain equivalent from Reef et al., 2018) could be used, and the resulting elevation could be added to the offshore transfer estimate obtain based on the ADCIRC tides. Thus, the production of bay water level predictions will require accurate wind forecast products and the quantification of the nonlinear interaction between local and remote effects.

The method presented offers a new methodology for coastal hazards assessment and risk analysis. While many methodologies are being used for open-coast regions (Thieler and Hammar-Klose, 1999; Kunreuther et al., 2000), vulnerability evaluation to coastal hazards in back-barrier bays remains underdeveloped. Evaluating bay hazards usually requires expensive computational simulations at appropriately high resolutions to characterize the spatial and temporal effects. The method presented here, using existing ADCIRC results, provides a less expensive approach that is able to properly estimate the spatial differences in vulnerability in response to flooding at different timescales (e.g., perigean spring tides, storms of different duration). It provides guidance for planning in response to “nuisance” flooding at a relatively low cost. It can be expanded to all back-barriers without the need to simulate each storm in each embayment while applying a consistent methodology.

Careful consideration needs to be given to the estimation of coastal hazards, especially for the forecast of intense storm effects. The inclusion of meticulously validated methodologies that consider both offshore influences (e.g., using the transfer estimated from ADCIRC tides) and local wind setup (e.g., Wong and Moses-Hall, 1998; Reef et al., 2018) is necessary. Skill assessment of storm hazard estimates using adequate observations is critical to avoid producing underpredictions or overpredictions of flooding and inundation.

As part of the general needs for hazard assessment (Ludwig et al., 2018), the important hazard characteristics that decision makers require include spatial extent, duration, and magnitude. The proposed methodology provides an approximation to both the extent of the area and the magnitude and also variations based on storm duration. Additionally, the fact that uncertainty estimates accompany the vulnerability provided by the present method enhances the potential value to decision makers. The extension to other bays in the United States will be included as part of the U.S. Geological Survey Coastal Change Hazards portal (USGS, 2018).

7 Summary

Back to toptop
The results presented here demonstrate a strategy for estimating the impact of offshore sea level and local wind setup in back-barrier bay water levels. The transfer estimates from offshore to the bay water level used a combination of observations, analytical models based on appropriate simplifications of the bay system, and numerical simulations that provide the needed spatial distribution and more realistic frictional control.

The resulting maps of water level response to offshore forcing showed larger attenuation of the relatively higher-frequency fluctuations such as the semidiurnal tides. Smaller transfers were associated with shorter duration storms than longer duration storms, and transfer was most spatially uniform for storms of long duration. The description of the magnitude and spatial dependence of transfer on storm duration will assist planning for flooding in back-barrier bays.

In the specific case of the Barnegat Bay–Little Egg Harbor system, larger transfers were estimated for the southern embayments (Great Bay and Little Egg Harbor) when compared to Barnegat Bay. The reason for the difference was the dominant role of Little Egg Inlet (wider and deeper) in controlling the exchange between the offshore and bay systems. During relatively small storms, the contribution of local wind to the bay water level setup was smaller than the transfer from offshore fluctuations. During intense events, like hurricanes, local wind setup was of the same order of magnitude or even larger than offshore influences, depending on wind magnitude and especially on the relative angle of the wind with respect to bay orientation.

We introduced two approaches that depend on the availability of observations and numerical solutions. The approach requiring fewer data based on the ADCIRC tidal database provides spatial offshore transfer estimates and measures of uncertainty. In both cases, the inclusion of the local wind setup could be achieved based on simple surge analytical estimates. The approach that includes an analytical model allows for a simple tool to study the response of back-barrier bay systems to alternative conditions and forcing (e.g., geomorphic changes, changing duration of storms, sea level rise).

The proposed method represents an effective and inexpensive approach to flooding hazard evaluation in back-barrier bays and inland waters. The method provides detailed spatial estimates of vulnerabilities and uncertainties that could be an intuitive tool for coastal managers.

Data availability

Back to toptop
Data availability.

The model data used are available in the permanent repository from Defne and Ganju (2018) and through the catalog USGS (2019). The observational data used are listed in the references, tables and the repository at http://waterdata.usgs.gov/nwis/inventory/?site_no=xxx (last access: 15 August 2019), where xxx stands for the USGS station number in Table 1.

Appendix A

Back to toptop
The offshore impact on the water level of the bay can be approximated with an analytical model that solves the linearized depth-averaged momentum equations. The system of equations for an idealized simplification of Barnegat Bay (Fig. 2) that includes five sub-embayments (based on constrictions inside the bay system) consists of 13 equations and unknowns.

$$\begin{array}{}\text{(A1)}& \begin{array}{rl}{\displaystyle \frac{\partial}{\partial t}}\left(\begin{array}{l}{u}_{\mathrm{1}}\\ {u}_{\mathrm{2}}\\ {u}_{\mathrm{3}}\\ {u}_{\mathrm{4}}\\ {u}_{\mathrm{5}}\\ {u}_{\mathrm{6}}\\ {u}_{\mathrm{7}}\\ {u}_{\mathrm{8}}\end{array}\right)& =g\left(\begin{array}{l}\mathrm{1}/{L}_{\mathrm{1}}\\ \mathrm{1}/{L}_{\mathrm{2}}\\ \mathrm{1}/{L}_{\mathrm{3}}\\ \mathrm{1}/{L}_{\mathrm{4}}\\ \mathrm{1}/{L}_{\mathrm{5}}\\ \mathrm{1}/{L}_{\mathrm{6}}\\ \mathrm{1}/{L}_{\mathrm{7}}\\ \mathrm{1}/{L}_{\mathrm{8}}\end{array}\right)\left(\begin{array}{l}{\mathit{\eta}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{LEI}}-{\mathit{\eta}}_{\mathrm{1}}\\ {\mathit{\eta}}_{\mathrm{1}}-{\mathit{\eta}}_{\mathrm{2}}\\ {\mathit{\eta}}_{\mathrm{2}}-{\mathit{\eta}}_{\mathrm{3}}\\ {\mathit{\eta}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{BI}}-{\mathit{\eta}}_{\mathrm{3}}\\ {\mathit{\eta}}_{\mathrm{3}}-{\mathit{\eta}}_{\mathrm{4}}\\ {\mathit{\eta}}_{\mathrm{4}}-{\mathit{\eta}}_{\mathrm{5}}\\ {\mathit{\eta}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{breach}}-{\mathit{\eta}}_{\mathrm{5}}\\ {\mathit{\eta}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{PPC}}-{\mathit{\eta}}_{\mathrm{5}}\end{array}\right)\\ & -\left(\begin{array}{l}{u}_{\mathrm{1}}{r}_{\mathrm{1}}/{h}_{\mathrm{1}}\\ {u}_{\mathrm{2}}{r}_{\mathrm{2}}/{h}_{\mathrm{2}}\\ {u}_{\mathrm{3}}{r}_{\mathrm{3}}/{h}_{\mathrm{3}}\\ {u}_{\mathrm{4}}{r}_{\mathrm{4}}/{h}_{\mathrm{4}}\\ {u}_{\mathrm{5}}{r}_{\mathrm{5}}/{h}_{\mathrm{5}}\\ {u}_{\mathrm{6}}{r}_{\mathrm{6}}/{h}_{\mathrm{6}}\\ {u}_{\mathrm{7}}{r}_{\mathrm{7}}/{h}_{\mathrm{7}}\\ {u}_{\mathrm{8}}{r}_{\mathrm{8}}/{h}_{\mathrm{8}}\end{array}\right)\end{array}\end{array}$$

$$\begin{array}{}\text{(A2)}& \begin{array}{rl}& \left(\begin{array}{l}{A}_{\mathrm{1}}\\ {A}_{\mathrm{2}}\\ {A}_{\mathrm{3}}\\ {A}_{\mathrm{4}}\\ {A}_{\mathrm{5}}\end{array}\right){\displaystyle \frac{\partial}{\partial t}}\left(\begin{array}{l}{\mathit{\eta}}_{\mathrm{1}}\\ {\mathit{\eta}}_{\mathrm{2}}\\ {\mathit{\eta}}_{\mathrm{3}}\\ {\mathit{\eta}}_{\mathrm{4}}\\ {\mathit{\eta}}_{\mathrm{5}}\end{array}\right)=\\ & \left(\begin{array}{l}{h}_{\mathrm{1}}{W}_{\mathrm{1}}{u}_{\mathrm{1}}-{h}_{\mathrm{2}}{W}_{\mathrm{2}}{u}_{\mathrm{2}}\\ {h}_{\mathrm{2}}{W}_{\mathrm{2}}{u}_{\mathrm{2}}-{h}_{\mathrm{3}}{W}_{\mathrm{3}}{u}_{\mathrm{3}}\\ {h}_{\mathrm{3}}{W}_{\mathrm{3}}{u}_{\mathrm{3}}+{h}_{\mathrm{4}}{W}_{\mathrm{4}}{u}_{\mathrm{4}}-{h}_{\mathrm{5}}{W}_{\mathrm{5}}{u}_{\mathrm{5}}\\ {h}_{\mathrm{5}}{W}_{\mathrm{5}}{u}_{\mathrm{5}}-{h}_{\mathrm{6}}{W}_{\mathrm{6}}{u}_{\mathrm{6}}\\ {h}_{\mathrm{6}}{W}_{\mathrm{6}}{u}_{\mathrm{6}}+{h}_{\mathrm{7}}{W}_{\mathrm{7}}{u}_{\mathrm{7}}+{h}_{\mathrm{8}}{W}_{\mathrm{8}}{u}_{\mathrm{8}}\end{array}\right)\end{array}\end{array}$$

${\mathit{\varphi}}_{\mathrm{LEI}},{\mathit{\varphi}}_{\mathrm{BI}},{\mathit{\varphi}}_{\mathrm{breach}}$, and
*ϕ*_{PPC} are linear frequency-dependent relationships between the water
levels at offshore proxy stations (Sandy Hook or Atlantic City) and the
water level just offshore of Little Egg Inlet, Barnegat Inlet, the breach at
Mantoloking caused by Hurricane Sandy, and the Point Pleasant Canal.

By performing Fourier transforms on the momentum equations ($\mathit{\eta}=\stackrel{\mathrm{\u0303}}{\mathit{\eta}}{e}^{i\mathit{\omega}t}$ and $u=\stackrel{\mathrm{\u0303}}{u}{e}^{i\mathit{\omega}t}$, where $\stackrel{\mathrm{\u0303}}{\mathit{\eta}}$ and $\stackrel{\mathrm{\u0303}}{u}$ represent the magnitude of the water level and velocity oscillations, respectively), we obtain the following equations.

$$\begin{array}{}\text{(A3)}& \begin{array}{rl}i\mathit{\omega}\left(\begin{array}{l}\stackrel{\mathrm{\u0303}}{{u}_{\mathrm{1}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{2}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{3}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{4}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{5}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{6}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{7}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{8}}}\end{array}\right)& =g\left(\begin{array}{l}\mathrm{1}/{L}_{\mathrm{1}}\\ \mathrm{1}/{L}_{\mathrm{2}}\\ \mathrm{1}/{L}_{\mathrm{3}}\\ \mathrm{1}/{L}_{\mathrm{4}}\\ \mathrm{1}/{L}_{\mathrm{5}}\\ \mathrm{1}/{L}_{\mathrm{6}}\\ \mathrm{1}/{L}_{\mathrm{7}}\\ \mathrm{1}/{L}_{\mathrm{8}}\end{array}\right)\left(\begin{array}{l}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{LEI}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{BI}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{breach}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{PPC}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}\end{array}\right)\\ & -\left(\begin{array}{l}\stackrel{\mathrm{\u0303}}{{u}_{\mathrm{1}}}{r}_{\mathrm{1}}/{h}_{\mathrm{1}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{2}}}{r}_{\mathrm{2}}/{h}_{\mathrm{2}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{3}}}{r}_{\mathrm{3}}/{h}_{\mathrm{3}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{4}}}{r}_{\mathrm{4}}/{h}_{\mathrm{4}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{5}}}{r}_{\mathrm{5}}/{h}_{\mathrm{5}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{6}}}{r}_{\mathrm{6}}/{h}_{\mathrm{6}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{7}}}{r}_{\mathrm{7}}/{h}_{\mathrm{7}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{8}}}{r}_{\mathrm{8}}/{h}_{\mathrm{8}}\end{array}\right),\end{array}\end{array}$$

and then

$$\begin{array}{}\text{(A4)}& \begin{array}{rl}\left(\begin{array}{l}\stackrel{\mathrm{\u0303}}{{u}_{\mathrm{1}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{2}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{3}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{4}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{5}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{6}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{7}}}\\ \stackrel{\mathrm{\u0303}}{{u}_{\mathrm{8}}}\end{array}\right)& =g\left(\begin{array}{l}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{LEI}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{BI}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{breach}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}\\ {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{PPC}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}\end{array}\right)/\\ & \left[\left(\begin{array}{l}{L}_{\mathrm{1}}\\ {L}_{\mathrm{2}}\\ {L}_{\mathrm{3}}\\ {L}_{\mathrm{4}}\\ {L}_{\mathrm{5}}\\ {L}_{\mathrm{6}}\\ {L}_{\mathrm{7}}\\ {L}_{\mathrm{8}}\end{array}\right)\left(\begin{array}{l}i\mathit{\omega}+{r}_{\mathrm{1}}/{h}_{\mathrm{1}}\\ i\mathit{\omega}+{r}_{\mathrm{2}}/{h}_{\mathrm{2}}\\ i\mathit{\omega}+{r}_{\mathrm{3}}/{h}_{\mathrm{3}}\\ i\mathit{\omega}+{r}_{\mathrm{4}}/{h}_{\mathrm{4}}\\ i\mathit{\omega}+{r}_{\mathrm{5}}/{h}_{\mathrm{5}}\\ i\mathit{\omega}+{r}_{\mathrm{6}}/{h}_{\mathrm{6}}\\ i\mathit{\omega}+{r}_{\mathrm{7}}/{h}_{\mathrm{7}}\\ i\mathit{\omega}+{r}_{\mathrm{8}}{h}_{\mathrm{8}}\end{array}\right)\right].\end{array}\end{array}$$

Performing the Fourier transform on the continuity equations (Eq. A2) and substituting the velocity values from Eq. (A4) we obtain the following equation.

$$\begin{array}{}\text{(A5)}& \begin{array}{rl}& i\mathit{\omega}\left(\begin{array}{l}{A}_{\mathrm{1}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}\\ {A}_{\mathrm{2}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}\\ {A}_{\mathrm{3}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}\\ {A}_{\mathrm{4}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}\\ {A}_{\mathrm{5}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}\end{array}\right)=g\\ & \left(\begin{array}{l}{h}_{\mathrm{1}}{W}_{\mathrm{1}}\left({\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{LEI}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}\right)/{L}_{\mathrm{1}}\left(i\mathit{\omega}+{r}_{\mathrm{1}}/{h}_{\mathrm{1}}\right)-{h}_{\mathrm{2}}{W}_{\mathrm{2}}\left({\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}\right)/\\ {L}_{\mathrm{2}}\left(i\mathit{\omega}+{r}_{\mathrm{2}}/{h}_{\mathrm{2}}\right)\\ {h}_{\mathrm{2}}{W}_{\mathrm{2}}\left({\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}\right)/{L}_{\mathrm{2}}\left(i\mathit{\omega}+{r}_{\mathrm{2}}/{h}_{\mathrm{2}}\right)-{h}_{\mathrm{3}}{W}_{\mathrm{3}}\left({\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}\right)/\\ {L}_{\mathrm{3}}\left(i\mathit{\omega}+{r}_{\mathrm{3}}/{h}_{\mathrm{3}}\right)\\ {h}_{\mathrm{3}}{W}_{\mathrm{3}}\left({\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}\right)/{L}_{\mathrm{3}}\left(i\mathit{\omega}+{r}_{\mathrm{3}}/{h}_{\mathrm{3}}\right)+{h}_{\mathrm{4}}{W}_{\mathrm{4}}\left({\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{BI}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}\right)/\\ {L}_{\mathrm{4}}\left(i\mathit{\omega}+{r}_{\mathrm{4}}/{h}_{\mathrm{4}}\right)-{h}_{\mathrm{5}}{W}_{\mathrm{5}}\left({\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}\right)/{L}_{\mathrm{5}}\left(i\mathit{\omega}+{r}_{\mathrm{5}}/{h}_{\mathrm{5}}\right)\\ {h}_{\mathrm{5}}{W}_{\mathrm{5}}\left({\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}\right)/{L}_{\mathrm{5}}\left(i\mathit{\omega}+{r}_{\mathrm{5}}/{h}_{\mathrm{5}}\right)-{h}_{\mathrm{6}}{W}_{\mathrm{6}}\left({\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}\right)/\\ {L}_{\mathrm{6}}\left(i\mathit{\omega}+{r}_{\mathrm{6}}/{h}_{\mathrm{6}}\right)\\ {h}_{\mathrm{6}}{W}_{\mathrm{6}}\left({\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}\right)/{L}_{\mathrm{6}}\left(i\mathit{\omega}+{r}_{\mathrm{6}}/{h}_{\mathrm{6}}\right)+{h}_{\mathrm{7}}{W}_{\mathrm{7}}\left({\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{breach}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}\right)/\\ {L}_{\mathrm{7}}\left(i\mathit{\omega}+{r}_{\mathrm{7}}/{h}_{\mathrm{7}}\right)+{h}_{\mathrm{8}}{W}_{\mathrm{8}}\left({\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{0}}{\mathit{\varphi}}_{\mathrm{PPC}}-{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}\right)/{L}_{\mathrm{8}}\left(i\mathit{\omega}+{r}_{\mathrm{8}}/{h}_{\mathrm{8}}\right)\end{array}\right).\end{array}\end{array}$$

Then, we can define ${K}_{n}=\frac{{h}_{n}{W}_{n}g}{{L}_{n}\left(\frac{{r}_{n}}{{h}_{n}}+i\mathit{\omega}\right)}$ for each channel *n*=1, …,
8 and with the proper rearrangement we get the following equation:

$$\begin{array}{}\text{(A6)}& \left(\begin{array}{l}i\mathit{\omega}{A}_{\mathrm{1}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}={K}_{\mathrm{1}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{LEI}}-{K}_{\mathrm{1}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}-{K}_{\mathrm{2}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}+{K}_{\mathrm{2}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}},\\ i\mathit{\omega}{A}_{\mathrm{2}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}={K}_{\mathrm{2}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}-{K}_{\mathrm{2}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}-{K}_{\mathrm{3}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}+{K}_{\mathrm{3}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}},\\ i\mathit{\omega}{A}_{\mathrm{3}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}={K}_{\mathrm{3}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}-{K}_{\mathrm{3}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}+{K}_{\mathrm{4}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{BI}}-{K}_{\mathrm{4}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}},\\ -{K}_{\mathrm{5}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}+{K}_{\mathrm{5}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}},\\ i\mathit{\omega}{A}_{\mathrm{4}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}={K}_{\mathrm{5}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}-{K}_{\mathrm{5}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}-{K}_{\mathrm{6}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}+{K}_{\mathrm{6}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}},\\ i\mathit{\omega}{A}_{\mathrm{5}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}={K}_{\mathrm{6}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}-{K}_{\mathrm{6}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}+{K}_{\mathrm{7}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{breach}},\\ -{K}_{\mathrm{7}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}+{K}_{\mathrm{8}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{PPC}}-{K}_{\mathrm{8}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}\end{array}\right).\end{array}$$

The system of equations can be solved by substitution.

$$\begin{array}{}\text{(A7)}& {\displaystyle}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{5}}={\displaystyle \frac{\left({K}_{\mathrm{7}}{\mathit{\varphi}}_{\mathrm{breach}}+{K}_{\mathrm{8}}{\mathit{\varphi}}_{\mathrm{PPC}}\right){\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}+{K}_{\mathrm{6}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}}{i\mathit{\omega}{A}_{\mathrm{5}}+{K}_{\mathrm{6}}+{K}_{\mathrm{7}}+{K}_{\mathrm{8}}}},\text{(A8)}& {\displaystyle}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{4}}={\displaystyle \frac{{K}_{\mathrm{5}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}+\frac{{K}_{\mathrm{6}}\left({K}_{\mathrm{7}}{\mathit{\varphi}}_{\mathrm{breach}}+{K}_{\mathrm{8}}{\mathit{\varphi}}_{\mathrm{PPC}}\right){\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}}{i\mathit{\omega}{A}_{\mathrm{5}}+{K}_{\mathrm{6}}+{K}_{\mathrm{7}}+{K}_{\mathrm{8}}}}{i\mathit{\omega}{A}_{\mathrm{4}}+{K}_{\mathrm{5}}+{K}_{\mathrm{6}}-\frac{{K}_{\mathrm{6}}{K}_{\mathrm{6}}}{i\mathit{\omega}{A}_{\mathrm{5}}+{K}_{\mathrm{6}}+{K}_{\mathrm{7}}+{K}_{\mathrm{8}}}}},\text{(A9)}& {\displaystyle}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{1}}={\displaystyle \frac{{K}_{\mathrm{1}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{LEI}}+{K}_{\mathrm{2}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}}{i\mathit{\omega}{A}_{\mathrm{1}}+{K}_{\mathrm{1}}+{K}_{\mathrm{2}}}},\text{(A10)}& {\displaystyle}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{2}}={\displaystyle \frac{{K}_{\mathrm{3}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}+\frac{{K}_{\mathrm{2}}{K}_{\mathrm{1}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{LEI}}}{i\mathit{\omega}{A}_{\mathrm{1}}+{K}_{\mathrm{1}}+{K}_{\mathrm{2}}}}{i\mathit{\omega}{A}_{\mathrm{2}}+{K}_{\mathrm{2}}+{K}_{\mathrm{3}}-\frac{{K}_{\mathrm{2}}{K}_{\mathrm{2}}}{i\mathit{\omega}{A}_{\mathrm{1}}+{K}_{\mathrm{1}}+{K}_{\mathrm{2}}}}},\end{array}$$

and finally,

$$\begin{array}{}\text{(A11)}& {\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{3}}={\displaystyle \frac{{\displaystyle \begin{array}{c}{K}_{\mathrm{4}}{\mathit{\eta}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{BI}}+{\displaystyle \frac{\frac{{K}_{\mathrm{1}}{K}_{\mathrm{2}}{K}_{\mathrm{3}}{\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}{\mathit{\varphi}}_{\mathrm{LEI}}}{i\mathit{\omega}{A}_{\mathrm{1}}+{K}_{\mathrm{1}}+{K}_{\mathrm{2}}}}{i\mathit{\omega}{A}_{\mathrm{2}}+{K}_{\mathrm{2}}+{K}_{\mathrm{3}}-\frac{{K}_{\mathrm{2}}{K}_{\mathrm{2}}}{i\mathit{\omega}{A}_{\mathrm{1}}+{K}_{\mathrm{1}}+{K}_{\mathrm{2}}}}}\\ +\frac{\frac{{K}_{\mathrm{5}}{K}_{\mathrm{6}}\left({K}_{\mathrm{7}}{\mathit{\varphi}}_{\mathrm{breach}}+{K}_{\mathrm{8}}{\mathit{\varphi}}_{\mathrm{PPC}}\right){\stackrel{\mathrm{\u0303}}{\mathit{\eta}}}_{\mathrm{o}}}{i\mathit{\omega}{A}_{\mathrm{5}}+{K}_{\mathrm{6}}+{K}_{\mathrm{7}}+{K}_{\mathrm{8}}}}{i\mathit{\omega}{A}_{\mathrm{4}}+{K}_{\mathrm{5}}+{K}_{\mathrm{6}}-\frac{{K}_{\mathrm{6}}{K}_{\mathrm{6}}}{i\mathit{\omega}{A}_{\mathrm{5}}+{K}_{\mathrm{6}}+{K}_{\mathrm{7}}+{K}_{\mathrm{8}}}}\end{array}}}{{\scriptscriptstyle \begin{array}{c}i\mathit{\omega}{A}_{\mathrm{3}}+{K}_{\mathrm{3}}+{K}_{\mathrm{4}}+{K}_{\mathrm{5}}-\frac{{K}_{\mathrm{3}}{K}_{\mathrm{3}}}{i\mathit{\omega}{A}_{\mathrm{2}}+{K}_{\mathrm{2}}+{K}_{\mathrm{3}}-\frac{{K}_{\mathrm{2}}{K}_{\mathrm{2}}}{i\mathit{\omega}{A}_{\mathrm{1}}+{K}_{\mathrm{1}}+{K}_{\mathrm{2}}}}\\ -{\scriptscriptstyle \frac{{K}_{\mathrm{5}}{K}_{\mathrm{5}}}{i\mathit{\omega}{A}_{\mathrm{4}}+{K}_{\mathrm{5}}+{K}_{\mathrm{6}}-\frac{{K}_{\mathrm{6}}{K}_{\mathrm{6}}}{i\mathit{\omega}{A}_{\mathrm{5}}+{K}_{\mathrm{6}}+{K}_{\mathrm{7}}+{K}_{\mathrm{8}}}}}\end{array}}}}.\end{array}$$

As it only depends on the offshore water level, the solution for the water level of the central sub-embayment (Eq. A11) can be used to recursively calculate the solutions for the rest of the sub-embayments by substituting the water level estimates into Eqs. (A7), (A8), (A9), and (A10).

Author contributions

Back to toptop
Author contributions.

ALA developed the methodology and wrote most of the manuscript. NKG and RPS suggested improvements to the methodology. ZD contributed to the figure generation, the COAWST simulation analysis, and the discussion of results. All authors contributed to the final version.

Competing interests

Back to toptop
Competing interests.

The authors declare that they have no conflict of interest.

Disclaimer

Back to toptop
Disclaimer.

Use of firm and product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

Acknowledgements

Back to toptop
Acknowledgements.

The authors thank Brad Butman for his helpful comments on the manuscript. This work was supported by the U.S. Geological Survey, Coastal and Marine Hazards/Resources Program. The hydrodynamic model outputs used in this study are available from Defne and Ganju (2018). The numerical model is the open-source model COAWST (Warner et al., 2019) available from https://code.usgs.gov/coawstmodel/COAWST (last access: 15 August 2019) after registration with John C. Warner (jcwarner@usgs.gov) at the U.S. Geological Survey. The data used are listed in the references, tables, and the repository at http://waterdata.usgs.gov/nwis/inventory/?site_no=xxx (last access: 15 August 2019), where xxx stands for the USGS station number in Table 1. The paper benefited from conversations with Chris Sherwood, Tarandeep Kalra, and John Warner from USGS.

Financial support

Back to toptop
Financial support.

This work was supported by the US Geological Survey, Coastal and Marine Hazards/Resources Program.

Review statement

Back to toptop
Review statement.

This paper was edited by Ira Didenkulova and reviewed by three anonymous referees.

References

Back to toptop
Aretxabaleta, A. L., Butman, B., and Ganju, N. K.: Water level response in back-barrier bays unchanged following Hurricane Sandy, Geophys. Res. Lett., 41, 3163–3171, https://doi.org/10.1002/2014GL059957, 2014.

Aretxabaleta, A. L., Ganju, N. K., Butman, B., and Signell, R. P.: Observations and a linear model of water level in an interconnected inlet-bay system, J. Geophys. Res., 122, 2760–2780, https://doi.org/10.1002/2016JC012318, 2017.

Aretxabaleta, A. L., Doran, K. S., Long, J. W., Erikson, L., and Storlazzi, C. D.: Toward a national coastal hazard forecast of total water levels, Proceedings of the Coastal Sediments, 2019, 1373–1384, 2019.

Bendat, J. S. and Piersol, A. G.: Random data, Analysis and measurement procedures, John Wiley & Sons, New York, 566 pp., 1986.

Chant, R. J.: Tidal and subtidal motion in a shallow bar-built multiple inlet/bay system, J. Coast. Res., 32, 102–114, 2001.

Chuang, W.-S. and Swenson, E. M.: Subtidal Water level variations in Lake Pontchartrain, Louisiana, J. Geophys. Res., 86, 4198–4204, https://doi.org/10.1029/JC086iC05p04198, 1981.

Cioffi, F., Gallerano, F., and Napoli, E.: Three-dimensional numerical simulation of wind-driven flows in closed channels and basins, J. Hydraul. Res., 43, 290–301, 2005.

Csanady, G. T.: Wind-induced barotropic motions in long lakes, J. Phys. Oceanogr., 3, 429–438, 1973.

Defne, Z. and Ganju, N. K.: Quantifying the Residence Time and Flushing Characteristics of a Shallow, Back-Barrier Estuary: Application of Hydrodynamic and Particle Tracking Models, Estuar. Coast., 38, 1719–1734, https://doi.org/10.1007/s12237-014-9885-3, 2015.

Defne, Z. and Ganju, N. K.: USGS Barnegat Bay hydrodynamic model for March to September 2012: US Geological Survey data release, https://doi.org/10.5066/F7SB44QS, 2018.

Garvine, R. W.: A simple model of estuarine subtidal fluctuations forced by local and remote wind stress, J. Geophys. Res., 90, 11945–11948, 1985.

Gornitz, V. M., Daniels, R. C., White, T. W., and Birdwell, K. R.: The development of a coastal risk assessment database: vulnerability to sea-level rise in the US Southeast, J. Coast. Res., 327–338, 1994.

Hunter, J. R. and Hearn, C. J.: Lateral and vertical variations in the wind-driven circulation in long, shallow lakes, J. Geophys. Res., 92, 13106–13114, 1987.

Keulegan, G. H.: Tidal flow in entrances: water-level fluctuations of basins in communication with seas, Committee on Tidal Hydraulics, Technical Bulletin No. 14, U.S. Army Corps of Engineers, Vicksburg, Miss., 89 pp., 1967.

Klein, R. J. and Nicholls, R. J.: Assessment of coastal vulnerability to climate change, Ambio, 28, 182–187, 1999.

Kunreuther, H., Platt, R., Baruch, S., Bernknopf, R. L., Buckley, M., Burkett, V., Conrad, D., Davidson, T., Deutsch, K., Geis, D., and Jannereth, M.: The hidden costs of coastal hazards: Implications for risk assessment and mitigation, H. John Heinz III Center for Science, Economics, and the Environment, Panel on Risk, Vulnerability, and the True Costs of Coastal Hazards, Island Press, Washington, D.C., 220 pp., 2000.

Ludwig, K. A., Ramsey, D. W., Wood, N. J., Pennaz, A. B., Godt, J. W., Plant, N. G., Luco, N., Koenig, T. A., Hudnut, K. W., Davis, D. K., and Bright, P. R.: Science for a risky world – A U.S. Geological Survey plan for risk research and applications: U.S. Geological Survey Circular, 1444, 57 pp., https://doi.org/10.3133/cir1444, 2018.

Mattingly, K. S., McLeod, J. T., Knox, J. A., Shepherd, J. M., and Mote, T. L.: A climatological assessment of Greenland blocking conditions associated with the track of Hurricane Sandy and historical North Atlantic hurricanes, Int. J. Climatol., 35, 746–760, 2015.

McKenna, K. K., Farrell, S. C., and Gebert, J. A.: Hurricane Sandy: Beach-dune recovery at New Jersey Beach Profile Network (NJBPN) sites, Shore & Beach, 84, 5–17, 2016.

Moftakhari, H. R., AghaKouchak, A., Sanders, B. F., Allaire, M., and Matthew, R. A.: What is nuisance flooding? Defining and monitoring an emerging challenge, Water Resour. Res., 54, 4218–4227, 2018.

Neumann, B., Vafeidis, A. T., Zimmermann, J., and Nicholls, R. J.: Future coastal population growth and exposure to sea-level rise and coastal flooding-a global assessment, PloS One, 10, e0118571, https://doi.org/10.1371/journal.pone.0118571, 2015.

Nicholls, R. J., Wong, P. P., Burkett, V. R., Codignotto, J. O., Hay, J. E., McLean, R. F., Ragoonaden, S., and Woodroffe, C. D.: Coastal systems and low-lying areas, Climate Change 2007: Impacts, Adaptation and Vulnerability, in: Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Parry, M. L., Canziani, O. F., Palutikof, J. P., van der Linden, P. J., and Hanson, C. E., Cambridge University Press, Cambridge, UK, 315–356, 2007.

Nicholls, R. J., Hanson, S. E., Lowe, J. A., Warrick, R. A., Lu, X., and Long, A. J.: Sea-level scenarios for evaluating coastal impacts, WIRES-Clim. Change, 5, 129–150, https://doi.org/10.1002/wcc.253, 2014.

Pugh, D. T.: Tides, surges and mean sea-level, Wiley, Chichester, ISBN 0-471-91505-X, 1987.

Rahmstorf, S.: Rising hazard of storm-surge flooding, P. Natl. Acad. Sci. USA, 114, 11806–11808, https://doi.org/10.1073/pnas.1715895114, 2017.

Reef, K. R., Lipari, G., Roos, P. C., and Hulscher, S. J.: Time-varying storm surges on Lorentz's Wadden Sea networks, Ocean Dynam., 68, 1051–1065, 2018.

State Committee for the Zuiderzee: Verslag van de Staatscommissie Zuiderzee, Report, 1926 (in Dutch).

Szpilka, C., Dresback, K., Kolar, R., Feyen, J., and Wang, J.: Improvements for the Western North Atlantic, Caribbean and Gulf of Mexico ADCIRC Tidal Database (EC2015), J. Mar. Sci. Eng., 4, 72, https://doi.org/10.3390/jmse4040072, 2016.

Thieler, E. R. and Hammar-Klose, E. S.: National Assessment of Coastal Vulnerability to Sea-Level Rise: Preliminary Results for the US Atlantic Coast, Woods Hole, MA: United States Geological Survey (USGS), Open File Report 99-593, 1999.

USGS – US Geological Survey: Coastal Change Hazards Portal web page, available at: http://marine.usgs.gov/coastalchangehazardsportal/ (last access: 15 August 2019), 2018.

USGS – US Geological Survey: U.S. Geological Survey Woods Hole Coastal and Marine Science Center Hydrodynamic Model Simulations, available at: https://www.sciencebase.gov/catalog/item/5c62d00ae4b0fe48cb34c7ad, last access: 16 August 2019.

Vitousek, S., Barnard, P. L., Fletcher, C. H., Frazer, N., Erikson, L., and Storlazzi, C. D.: Doubling of coastal flooding frequency within decades due to sea-level rise, Sci. Rep.-UK, 7, 1399, https://doi.org/10.1038/s41598-017-01362-7, 2017.

Wahl, T., Haigh, I. D., Nicholls, R. J., Arns, A., Dangendorf, S., Hinkel, J., and Slangen, A. B.: Understanding extreme sea levels for broad-scale coastal impact and adaptation analysis, Nat. Commun., 8, 16075, https://doi.org/10.1038/ncomms16075, 2017.

Warner J. C., Armstrong, B., He, R., and Zambon, J. B.: Development of a Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) Modeling System, Ocean Model., 35, 230–244, https://doi.org/10.1016/j.ocemod.2010.07.010, 2010.

Warner, J. C., Ganju, N., Sherwood, C. R., Kalra, T., Aretxabaleta, A., Olabarrieta, M., He, R., Zambon, J., and Kumar, N.: A Coupled Ocean Atmosphere Wave Sediment Transport Numerical Modeling System (COAWST), USGS Github Code Repository, available at: https://code.usgs.gov/coawstmodel/COAWST, last access: 15 August 2019.

Wilkin, J. L. and Hunter, E. J.: An assessment of the skill of real-time models of Mid-Atlantic Bight continental shelf circulation, J. Geophys. Res., 118, 2919–2933, 2013.

Wong, K.-C. and DiLorenzo, J.: The response of Delaware's inland bays to ocean forcing, J. Geophys. Res., 93, 12525–12535, 1988.

Wong, K.-C. and Moses-Hall, J. E.: On the relative importance of the remote and local wind effects to the subtidal variability in a coastal plain estuary, J. Geophys. Res., 103, 18393–18404, 1998.

Wong, K.-C. and Wilson, R. E.: Observations of low-frequency variability in Great South Bay and relations to atmospheric forcing, J. Phys. Oceanogr., 14, 1893–1900, 1984.

Short summary

Water levels in bays are affected by open-ocean changes and wind. Tides are more dampened in the bays than storm surges and sea level rise. We compare observed and modeled levels with ocean conditions and combine them with analytical models. We consider the local setup, caused by wind along the bay. Expansion using the ADCIRC tidal database will allow coverage of other bay systems on the United States East Coast. Spatial estimates of water level can inform decisions about bay flooding hazards.

Water levels in bays are affected by open-ocean changes and wind. Tides are more dampened in the...

Natural Hazards and Earth System Sciences

An interactive open-access journal of the European Geosciences Union