Estuarine response to storm surge and sea-level rise associated with channel deepening: a flood vulnerability assessment of southwest Louisiana, USA

This study investigates the sensitivity of the Calcasieu Lake estuarine region to channel deepening in southwest Louisiana in the USA. We test the hypothesis that the depth increase in a navigational channel in an estuarine region results in the amplification of the inland penetration of storm surge, thereby increasing the flood vulnerability of the region. We run numerical experiments using the Delft3D modeling suite (validated with observational data) with different historic channel depth scenarios. Model results show that channel deepening facilitates increased water movement into the lake–estuary system during a storm surge event. The inland peak water level increases by 37% in the presence of the deepest channel. Moreover, the peak volumetric flow rate increases by 291.6% along the navigational channel. Furthermore, the tidal prism and the volume of surge prism passing through the channel inlet increase by 487% and 153.3%, respectively. In our study, the presence of the deepest channel results in extra 56.72 km2 of flooded area (approximately 12% increase) which is an indication that channel deepening over the years has rendered the region more vulnerable to hurricane-induced flooding. The study also analyzes the impact of channel deepening on storm surge in estuaries under different future sea-level rise (SLR) scenarios. Simulations suggest that even the most conservative scenario of SLR will cause an approximately 51% increase in flooded area in the presence of the deepest ship channel, thereby suggesting that rising sea level will cause increased surge penetration and increased flood risk.


Introduction
Global climate change is expected to intensify the strength of hurricanes (Emanuel 2005(Emanuel , 2020; Knutson et al. 2010; references therein and many others). Despite some debate, the frequency of hurricanes is anticipated to escalate as well (Holland and Webster 2007;Emanuel 2021;Vecchi et al. 2021;Klotzbach et al. 2022; references therein and many others). Moreover, rising sea levels are expected to exacerbate the flooding caused due to the storm surges associated with these tropical cyclones, posing a dire threat to life and property of coastal communities (Karim and Mimura 2008;Shepard et al. 2012;Tasnim et al. 2015; references therein and many others). Coincident with climate change are anthropogenic modifications to coastal areas-channelization and dredging for shipping purposes, reclamation, diking, filling of intertidal and shallow subtidal waters, hardening of shorelines and regulating river discharges. These anthropogenic developments have resulted in some adverse consequences, such as (1) increased salinity due to saltwater intrusion that threatens drinking water supplies (Kim and Johnson 1998;Siles-Ajamil et al. 2019); (2) circulation and exchange flow alterations affecting salinity (Zhu et al. 2015;Chant et al. 2018; Wang and Shen 2020); (3) increased risk of erosion due to increase in the tidal prism (Cox et al. 2003); (4) tidal range, current and energy alterations that changes sediment transport patterns (Pareja-Roman et al. 2020); (5) increased sediment concentration (Kerner 2007;van Maren et al. 2015;Dijkstra et al. 2019) and upstream movement of estuarine turbidity maximum (de Jonge et al. 2014;Burchard et al. 2018) affecting ecology and water quality; (6) decrease in dissolved oxygen (Kaur et al. 2007;Sangita et al. 2014); (7) flow and discharge augmentation that causes sediment regime alteration (Rahman and Yunus 2016); and (8) storm surge and tidal amplification (Familkhalili and Talke 2016;Familkhalili et al. 2020;Talke et al. 2021; references therein and many others).
Channel deepening can have varying effects on the water-level dynamics of a given location during any given storm event. Although it usually increases the water level resulting from storm surge and tidal amplification, it can lead to attenuated mean water level (Jensen et al. 2003;Jay et al. 2011;) and tidal amplitude (Cai et al. 2012a) as well, owing to increased hydraulic conveyance. There may even be no relationship between channel deepening and storm surge and tidal dynamics (Cai et al. 2012a;Bilskie 2013). Thus, the same process of channel deepening can both increase and decrease the total water level. The geometric properties of an estuary play a vital role in this regard since tidal dynamics and the friction to water movement depend on depth, cross-sectional area, width, convergence and other bathymetric properties (Friedrichs and Aubrey 1994;Orton et al. 2015;Familkhalili and Talke 2016).
Moreover, the impact of sea-level rise (SLR) on coastal flooding is uncertain on its own. The flood area does not always increase because of SLR, particularly in tidal river estuaries, where streamflows and the resulting hydraulic head in the upriver location control storm-driven flood levels (Orton et al. 2020) and where water levels are primarily influenced by the freshwater discharges (Chen and Liu 2016). Thus, climate change and its associated uncertainties suggest the importance of incorporating climate change effects in a coastal hazard risk assessment. In spite of increasing demand for the integration of climate change adaptation and disaster risk reduction, effective integration of the two in the disaster management planning process remains limited due to logistical problems associated with funding; lack of quality governance, coordination, collaboration, implementation, mainstreaming and community involvement; geopolitical factors; and 1 3 information, knowledge and policy gaps (Mitchell et al. 2010;Gero et al. 2011;Serrao-Neumann et al. 2015;Islam et al. 2020).
In this paper, we investigate the sensitivity of the Calcasieu Lake estuarine region in southwest Louisiana (Fig. 1) to the common anthropogenic modification of channel deepening. We study the changes in storm surge and flood vulnerability in the estuarine region that can be tied to water depth changes of the Calcasieu Ship Channel. This is critical in the sense that the Calcasieu Lake estuarine region is one of Louisiana's most at-risk hurricane zones and is prone to widespread flooding (Rego et al. 2010;Manning 2021). Regardless of the hurricane and flood risk, the region is one of America's most active energy corridors, with the country's 6th largest refinery, two of the country's largest liquefied natural gas (LNG) facilities, a complex of energy pipelines and the Henry Hub, a major natural gas hub for the USA. The port also serves as a big driver of the country's GDP, which is an estimated $39 billion (Calcasieu Ship Channel 2021;Martin Associates 2021). To sustain such heavy port activities, channel deepening is essential for navigation. However, the impact of channel deepening on ecological and environmental integrity is a major concern.
Despite all the previous studies on the various adverse consequences of channel deepening, little focus has been given on how the water level responds or storm surge propagates inland during a hurricane owing to the presence of a deepened channel in an estuarine region. Our study aims to fill this knowledge gap by testing the hypothesis that the deepening of the Calcasieu Ship Channel has increased the flood vulnerability of the Calcasieu Lake estuarine region. We accomplish this using the Delft3D modeling suite based on measured 1874,1903,1941 and 2020 channel bathymetric conditions, thereby testing our hypothesis on a broad, shell-shaped, sluggish lake-estuary system using measured meteorological forcing data of Hurricane Laura (2020). We also incorporate SLR effects in our study by analyzing the impact of channel deepening on storm surge in estuaries under future sea-level rise scenarios.

Study domain
Calcasieu Lake is a large, shallow, brackish lake in southwestern Louisiana (Lee et al. 1990;). The lake is surrounded by brackish marsh and runs approximately 28 km from north to south. The broad southern shore extends approximately 19 km in the east-west direction. Calcasieu Lake is connected to the Gulf of Mexico through a narrow, 9 km channel called the Calcasieu Pass. Apart from the channel, the estuarine region is shallow with an average depth of 2 m. The region is part of the Calcasieu-Sabine Basin.
Like many estuaries of the world, this lake-estuary system is no longer in its natural configuration owing to the Calcasieu Ship Channel, which is maintained at 40 feet (12.19 m) deep by 400 feet (121.92 m) wide and extends from the Gulf of Mexico to Lake Charles (CWPPRA 1992;McGinnis et al. 2019). The U.S. Army Corps of Engineers (USACE) first constructed the 5-feet-deep × 80-feet-wide × 7500-feet-long (1.52 m × 24.38 m × 2286 m) channel in 1874. Since 1874, the navigational channel has been dredged and widened several times leading to alterations in the hydrology and water quality of the estuarine region (McGinnis et al. 2019). In 1903, the channel was deepened to 13 feet (3.96 m) and was completed in 1941 with a depth of 30 feet (9.14 m) and width of 250 feet (76.2 m). Finally, in 1968, it was dredged to its current dimensions.

Hurricane Laura
Hurricane Laura (2020) was a destructive, Category 4 hurricane that made landfall in this estuarine region near Cameron in southwest Louisiana in the early hours of August 27, 2020 (Jafari et al. 2020;Roueche et al. 2021;Fig. 2). It is tied with the Last Island Hurricane of 1856 as the strongest hurricane that made landfall in the state of Louisiana in terms of maximum sustained wind. It is also tied as the fifth strongest hurricane to make landfall in the continental USA and has set the record for the fastest intensification rate in the Gulf of Mexico (Roueche et al. 2021). Hurricane Laura formed from a large tropical wave that originated off the West African coast on August 16, 2020, then converted into a tropical depression on August 20, 2020, and finally made landfall on early August 27, 2020 (U.S. Census Bureau 2020; Roueche et al. 2021). Southwest Louisiana and southeast Texas were hit by winds up to 240 km/h and storm surge more than 3 m causing an estimated $19 billion in damages. Moreover, 33 people died in the state of Louisiana alone (NCEI 2021).

Data, model setup and model validation
Delft3D is an integrated modeling suite that can operate in different modes, such as 2D, Q3D or 3D (Elias et al. 2001). The 2D mode of Delft3D solves the depth-integrated, incompressible, Reynolds-averaged Navier-Stokes equations, or the nonlinear shallowwater equations (Symonds et al. 2017;Deltares 2020). It has been shown to be an effective tool to study coastal processes, including storm surge and hydrodynamic response to channel deepening (Hu et al. 2015;Zarzuelo et al. 2015;Wang et al. 2017;Hopkins et al. 2018;Liu et al. 2018;Johnson et al. 2021; references therein and many others). We initialized Delft3D using the bathymetric data from USGS ( Fig. 1). We forced the model using the observed wind and water-level data obtained from the NOAA station 8768094 at Calcasieu Pass, LA. The time series of the wind data (wind speed and wind direction) and the water-level data are shown in Fig. 3. We resolved the channel bathymetry using the NOAA nautical charts, combined with the coarser USGS map. These inputs were applied to a structured grid with uniform size of 20 m which spans 68.2 km in the horizontal and 54.66 km in the vertical direction ( Fig. 1). We validated the model using observed water-level data from the USGS stations-08017118 in Calcasieu River at Cameron, LA and 08017095 in North Calcasieu Lake near Hackberry, LA ( Fig. 1). From here on out, we will mention USGS 08017118 as USGS south and USGS 08017095 as USGS north for simplicity.

Channel depth scenarios
Four Calcasieu Ship Channel depth scenarios were simulated using Hurricane Laura wind forcing and observed water level to study how changing conditions affect the storm surge and the corresponding flooding situations. The channel depth scenarios are defined according to specific channel depth of specific years as follows: We did not consider the potential climate change effects on precipitation, evapotranspiration, subsidence, storm frequency and storm intensity for the purpose of this study, isolating the effects of sea-level rise.
One thing to be noted here is that the tests on the channel deepening are projections backward (2020-1874), while scenarios for SLR are projections forward. Some of these scenarios is never expected to happen. For example, the scenario combining the channel depth of 1874 and SLR will never happen for this specific site. Nonetheless, we run these combined scenarios to study how the same navigational channel with different historic water depths will respond differently to the same SLR implications. The goal is to examine how channel deepening over the years made the region more vulnerable to coastal flooding and how future SLR will make the already bad situation worse.

Results and discussion
We observe satisfactory model-data agreement (Fig. 4). The model skillfully simulates the fluctuations and the patterns of the observed water level at the two USGS station locations, although it is not perfect in capturing the minimum and the maximum. The maximum is under-predicted, and the minimum is over-predicted along with a slight temporal difference, north of Calcasieu Lake (Fig. 4b). Considering the order of magnitude of the flood regime being studied (kilometers), these inaccuracies are expected to affect the flood vulnerability analysis only slightly.
Model results from different channel depth scenarios suggest that shallower channel blocks the water from going into the lake, whereas deeper channel facilitates water movement into the lake (Fig. 5). We believe that two processes drive this phenomenon-reduction in frictional resistance to water movement and increased channel cross-sectional area. In the momentum balance equation, the frictional term is inversely proportional to the depth (Jay 1991;Friedrichs and Aubrey 1994;Pullen et al. 2007;van Rijn 2011;Arns et al. 2017). At the same time, increase in channel depth increases the cross-sectional area of the channel which in turn increases the volume of water passing through the channel (Yorke 1978;Hughes 2002). As seen from the model results in Fig. 5a, the '5ft' (1.52 m) experiment produces a higher peak water level at the USGS-south station location in the ship channel. This peak gradually decreases for the other scenarios-lower for '13ft' (3.96 m) scenario and even lower for '30ft' (9.14 m) and '40ft' (12.19 m) cases (Fig. 5a). The USGS-north station shows the opposite-maximum peak water level for '40ft' scenario and gradually decreasing with the minimum being that for '5ft' scenario (Fig. 5b). In comparison with the '5ft' scenario, the peak water level at this location has increased by 0.23 m (19%), 0.39 m (32.2%) and 0.45 m (37%) for '13ft, ' '30ft' and '40ft' scenarios, respectively. Peak water levels at different locations (and different times, given that the peak occurs earlier in the south than the north as the storm moves inland) in the study domain also suggest increased water movement through the ship channel with higher water depth during Hurricane Laura (Fig. 6). The figures infer that the higher peak water level for lower channel depth at the channel-lake mouth is caused due to 'piling-up' effect where the surge entering the channel from Gulf of Mexico is piling-up at the channel-lake mouth and not making it north. We observe substantial increase in the peak water level at different locations in and around Lake Calcasieu for '40ft' scenario as compared to the '5ft' one ( Fig. 6d). This suggests that the channel depth can control storm surge movement into the lake.
Examining the momentum balance along the channel also suggests reduction in frictional resistance with increased channel depth and piling-up of water in the ship channel (Fig. 7). Generally, in the momentum balance equation for channels, pressure term (g ∂η/∂x; g = acceleration due to gravity; η = water level) and friction term (c d u |u|/R; c d = drag coefficient; R = hydraulic radius ≈ h if b > > h; h = channel depth; b = channel width) are the most important. Local acceleration (∂u/∂t; u = cross-sectionally averaged velocity) and convective acceleration (u ∂u/∂x) terms are relatively small (Friedrichs and Aubrey 1994;van Rijn 2011). Model results from different channel depth scenarios Fig. 6 a Peak water level relative to NAVD88 in meter at each grid point of the developed model in the presence of the '5ft' (1.52 m) channel after Hurricane Laura made landfall. The difference in the peak water level at each grid point between that of b '13ft' (3.96 m), c '30ft' (9.14 m) and d '40ft' (12.19 m) scenario, and the '5ft' scenario, with the red color meaning higher peak for the '5ft' scenario and blue meaning higher peak for the other scenarios. Notice that the water level reaches a higher peak in the channel-lake mouth for the '5ft' scenario, and at different locations in and around the lake for the other scenarios also show that the local acceleration term and the convective acceleration term approximately differ by an order of magnitude from the pressure term and the friction term (Fig. 7). As seen from Fig. 7c, water surface slope between the gulf and the channel drives the momentum balance as the surge water flows through the channel. We observe higher water surface slope along the channel for the lower channel depth scenario as compared to the greater channel depth scenario. This suggests that for a shallower channel, the excess water gets piled-up more at the bay entrance and in the channel instead of flowing to the north. On the other hand, as seen from Fig. 7d, the friction term along the channel decreases with the increase in the channel depth, which facilitates water movement through the channel and into the lake. Thus, the channel depth controls the flow of the storm surge in the lake-estuary system.
Further, model results show that peak volumetric flow rates along Calcasieu Ship Channel increase with the increase in depth (Fig. 8). The '40ft' scenario produces the highest flow rate along the channel, whereas the '5ft' scenario produces the lowest (orange versus purple in Fig. 8). On average, as compared to '5ft' scenario, the percentage increase in peak flow rate along the channel is 73.7% for the '13ft' scenario, 215.8% for the '30ft' scenario and 291.6% for the '40ft' scenario. Moreover, model results also show that tidal prism of the lake and the volume of surge water passing through the inlet increase with the increase in depth (Fig. 9). Tidal prism is the volume of water passing through an estuarine cross section or inlet on each tidal cycle (Lopes et al. 2013). Empirical relationships between tidal prism and channel or inlet crosssectional area have been established in previous studies, and subsequent studies have confirmed their validity (O'Brien 1931(O'Brien , 1969Jarrett 1976;Marchi 1990;D'Alpaos et al. 2010;Fig. 8 Peak flow rate at various locations along the ship channel. The x-axis represents the distance in kilometers from mouth of Calcasieu Ship Channel and Gulf of Mexico, and the y-axis shows the peak flow rate. The plots are color-coded according to channel depth-purple for '5ft' (1.52 m), green for '13ft' (3.96 m), blue for '30ft' (9.14 m) and orange for '40ft' (12.19 m), respectively. The black dashed line marks the start of Calcasieu Lake Fig. 9 a Tidal prism at Calcasieu Ship Channel inlet computed from simulation output (x-axis) versus empirical relationship (y-axis). b Volume of surge water passing through the same inlet cross section during Hurricane Laura computed from simulation output (x-axis) versus empirical relationship (y-axis). The plots are color-coded according to channel depth-purple for '5ft' (1.52 m), green for '13ft' (3.96 m), blue for '30ft' (9.14 m) and orange for '40ft' (12.19 m), respectively Hoitink et al. 2017; references therein and many others). The relationship, termed as the O'Brien-Jarrett-Marchi (OBJM) law (D'Alpaos et al. 2009(D'Alpaos et al. , 2010, states that, Ω = k*P α . Here, Ω is the channel or inlet cross-sectional area and P is the tidal prism with an exponent α = 6/7 and k = 1.2 × 10 −3 m 2−3α , when Ω and P are expressed in m 2 and m 3 , respectively (Marchi 1990;D'Alpaos et al. 2010). Tidal prism output from the Delft3D model of this study and the above empirical relationship are comparable, further validating the OBJM law for astronomical tides. However, in case of the storm tide, volume of surge water passing though the channel inlet, obtained from our model output, deviates slightly from the OBJM law. Nonetheless, it is found from the analysis that the '40ft' scenario produces the maximum tidal prism at the channel inlet, whereas the '5ft' scenario produces the minimum (orange versus purple in Fig. 9a). The trend is similar for the volume of surge water passing though the inlet: maximum for the '40ft' scenario and minimum for the '5ft' scenario (orange versus purple in Fig. 9b). Compared to the '5ft' scenario, the percentage increase in tidal prism is 90% for the '13ft' scenario, 320.5% for the '30ft' scenario and 487% for the '40ft' scenario. The corresponding percentage increase in the volume of surge water passing through the channel inlet is 22.6% for the '13ft' scenario, 78.1% for the '30ft' scenario and 153.3% for the '40ft' scenario.
Finally, simulations show an increase in flooded area for the '40ft' scenario as compared to other scenarios (Figs. 10, 11a). Compared to '5ft' experiment, each scenario exhibits notable increase in flooded area (blue, cyan and black areas in Fig. 10); the maximum is being observed for the '40ft' experiment. Area of new induced flooding for each of '13ft,' '30ft' and '40ft' scenario as compared to the '5ft' scenario is 24.9 km 2 (5.4% increase in the flooded area), 45.57 km 2 (9.8% increase) and 56.72 km 2 (12.2% increase), respectively. Some induced drying is also observed for each scenario near the coast and southern shore Fig. 10 New areas of induced flooding and drying for '13ft' (3.96 m), '30ft' (9.14 m) and '40ft' (12.19 m) scenario as compared to '5ft' (1.52 m) scenario. Blue color shows new areas of induced flooding for '13ft' scenario as compared to '5ft' scenario. Blue and cyan colors together are areas of new flooding for '30ft' scenario. Blue, cyan and black areas together are new flooding areas for '40ft' scenario. Red color shows new areas of induced drying for all three scenarios together when compared with '5ft' scenario. Green contour lines indicate the land edge of Lake Calcasieu (red areas in Fig. 10), likely because deeper channel conveys more water flow beyond the channel and into the lake, whereas shallower channel blocks the flow causing water to pile up inside the channel and leads to flooding of the overlying areas. However, areas of induced drying are negligible compared to those of induced flooding.
The wetlands in the surrounding area are expected to act as a vanguard against hurricane-induced flooding by preventing the storm surge from propagating inland. The potential of wetlands to mitigate hurricane-induced storm surge and the associated flooding is widely recognized (Loder et al. 2009;Wamsley et al. 2009;Barbier et al. 2013;Haddad et al. 2016;Zhao and Chen 2016;Narayan et al. 2017; references therein and many others). However, our results suggest that channel deepening can lead to increased storm surge penetration even in the presence of wetland vegetation, with the flooded area inland increasing drastically from changing the channel depth alone (Fig. 10).
Globally, channel deepening has been seen to produce alterations of tidal and storm surge amplitudes. Previous studies show instances of tidal amplification owing to channel deepening in various locations around the world including the Scheldt, Elbe, Ems and Loire River estuaries in Europe , the Thames in the UK (Amin 1983), the Rhine-Meuse in Europe (Vellinga et al. 2014), the Weser estuary in Bremen, Germany , the Gironde estuary in France (Jalón-Rojas et al. 2018), the Cape Fear River Estuary in North Carolina, USA (Familkhalili and Talke 2016), the St. Johns River in Florida, USA (Bilskie 2013;Talke et al. 2021), the Chesapeake Bay estuary in the USA (Ross et al. 2017), the Pearl River estuary in China (Zhang et al. 2010, the Delaware River estuary on the East coast of the USA (DiLorenzo et al. 1993;Ross et al. 2017;Pareja-Roman et al. 2020), the Rotterdam waterway in the Netherlands , the Hudson River in the USA , the Modaomen waterway in China (Cai et al. 2012b), the Columbia River in North America (Jay et al. 2011;Helaire et al. 2019), the Bay of Cádiz in Spain (Zarzuelo et al. 2015), Tampa Bay estuary in Florida, USA (Zhu et al. 2015), Newark Bay in New Jersey, USA (Chant et al. 2018), the Guadalquivir River estuary in Spain (Siles-Ajamil et al. 2019) and James River in Virginia, USA (Wang and Shen 2020). Storm surge, being shallow-water long waves like tides, undergoes amplification  (Familkhalili and Talke 2016), Delaware Bay on the East Coast of the USA (Familkhalili et al. 2020) and Saint Johns River Estuary in Florida (Talke et al. 2021). Moreover, linkage between navigational channel deepening and flood risk has been investigated in some studies. Deepening of the Ria de Aveiro lagoon main channels in Portugal has been seen to induce an increase in the lagoon flood extent area (Lopes et al. 2013). In another study, it was suggested that coastal flooding in New York City can be mitigated by shallowing the channels around the city (Orton et al. 2015). A case study of the Pearl River Delta shows that channel deepening due to sand mining has increased peak water level in the tide-dominated coastal regions and the associated flood risk (Bao et al. 2022). An analytical model study also shows that deepening of an estuarine channel results in larger tidal amplitude and greater surge height causing large changes in the amplitude and spatial distribution of coastal flooding (Familkhalili et al. 2022). However, by contrast, channel deepening may reduce the flood level caused by river discharge, owing to increased hydraulic conveyance (Jensen et al. 2003;Jay et al. 2011;Helaire et al. 2019;. For example, in Albany, New York in the USA, located in the upper tidal river, although the tidal range has doubled and storm surge magnitude has increased over the past 150 years, the flood risk caused by large river discharge has significantly decreased . Our findings are consistent with other studies that find that channel deepening facilitates inland storm surge propagation and increases the associated flood risk (Cai et al. 2012b;Lopes et al. 2013;Orton et al. 2015; references therein and many others).
The situation worsens under future SLR implications (Figs. 11b, 12). In the presence of the 5 ft (1.52 m) deep ship channel, the area of new flooding increases approximately 21%, 34% and 49% for the low, medium and high SLR scenarios respectively, when compared to no SLR (blue, cyan and black areas in Fig. 12a). With a 40ft (12.19 m) deep channel, we see even more flooding; again, for low, medium and high: approximately 36%, 60% and 86% increases, respectively (blue, cyan and black areas in

Conclusions
This study investigates the flood risk associated with channel deepening in an estuarine system during an extreme hurricane. Results from our numerical experiments suggest that channel deepening over the years has rendered the Calcasieu Lake estuarine region more vulnerable to hurricane-induced flooding. The presence of a deeper navigational channel facilitates water movement through the channel, resulting in increased surge penetration and increased flood risk. Flooded area increases almost 12% in case of the deepest channel scenario. Simulations also suggest that rising sea level is expected to make the already bad situation worse. According to our results, even the most conservative scenario of SLR will cause an approximately 51% increase in flooded area in the presence of the deepest ship channel. Considering the economic driver and other societal interests, it may be impractical not to administer dredging projects across the country. But proper research and cost-benefit analysis are imperative before implementing such projects, taking into consideration the damages resulting from increased flood events during extreme weather conditions, especially because of the current global climate change and rising sea levels. The Calcasieu Lake estuarine region may not be populated enough to cause as much damage as other densely populated coastal communities during such coastal flooding events. In a densely populated and highly urbanized coastal city, such as Boston or New York City, the damages can be ghastly. Numerous impervious parcels including unsewered paved roads and parking lots in urbanized areas will only add to the problem (Park et al. 2020;Blum et al. 2020;Feng et al. 2021). Therefore, the importance of additional state-of-the-art research into the matter cannot be overstated and needs to be continued.
Funding Open access funding provided by Northeastern University Library. Funding for the study has been provided in part by the U.S. National Science Foundation (NSF Grants 2139882). Maqsood Mansur's work was partially supported by the Northeastern University Department of Civil and Environmental Engineering Fellowship.
Data availability Dataset containing the input data, model setup and simulation output is freely available in DesignSafe-CI (Rathje et al. 2017;Mansur et al. 2022).