Compound Inundation Impacts of Coastal Climate Change: Sea-Level Rise, Groundwater Rise, and Coastal Precipitation

Abstract

The importance of considering the compound effects of multiple hazards has increased in recent years due to their catastrophic impacts on human lives and property. Compound effects correspond to events with multiple concurrent or consecutive drivers, e.g., heavy storms, coastal flooding, high tides, and sea level rise (SLR). There is a recent evidence on inundation caused by SLR-driven groundwater rise, and there is a distinct knowledge gap in understanding the compound inundation effects of this phenomenon considering the important hydrologic and hydraulic considerations under compound events. To fill this knowledge gap, we developed a novel analytical framework to understand the movements of the surface flow under typical precipitation events considering their interaction with uprising groundwater and SLR in a coastal watershed located in Oakland Flatlands, CA, USA, home to several disadvantaged communities. This modelling approach simulates the dynamics of compound flooding in two dimensions of the earth surface in a fine resolution, which is critical for devising proper flood management strategies. The reason to focus on disadvantaged coastal communities is that such communities typically encounter disproportionate environmental injustices due to the lack of sufficient drainage capacity in their infrastructure. Our results show that by considering the compound effect of SLR, groundwater inundation and precipitation flooding, the drainage capacity of infrastructure will be substantially exceeded, such that over 700 acres of the built infrastructure could be flooded. This is a considerable increase compared to scenarios that do not consider compound effect, or scenarios that consider inappropriate combinations of driving factors. In sum, our results highlight the significance of considering compound effects in the coastal inundation analyses, with a particular emphasis on the role of groundwater rise.

1. Introduction

Built infrastructure of coastal regions have been recently found vulnerable to marine inundation as a result of sea-level rise (SLR), particularly in disadvantaged communities. For instance, the sea-level changes that have already occurred in the City of Imperial Beach, California—which is identified as a Disadvantaged Community according to the California Department of Water Resources (DWR) with the majority of its population as Hispanic or Latino according to the 2019 census estimates—have made 5 major pump stations, many of the key stormwater outlets, and other water infrastructure elements, vulnerable to marine inundation [1]. Even 0.5 m of additional SLR will make miles of pipes, additional pump stations, stormwater outlets and system elements vulnerable to inundation in Imperial Beach [1]. Figure 1 illustrates all California communities that are expected to face SLR inundation challenges in this century [2]. The global mean sea level may rise 0.18–0.48 m by mid-century, and 0.5–1.4 m by the end of the century [3,4]. The recent SLR has been primarily attributed to greenhouse gas (GHG) emissions, which have accelerated the melting of polar ice caps and thermal expansion of sea water [3,4,5,6]. In California, records indicate an observed SLR rate of 17–20 cm per century over the past few decades, which is analogous to the global rise [4]. However, the rate of global SLR has accelerated in recent years, and projections suggest the potential for substantially greater SLR over this century and next [4,6,7]. Hauer, Evans [8] considered population growth in their assessment of the effect of SLR on communities in United States, where they found that a 2100 SLR of 0.9 m places a land area projected to house 4.2 million people at risk of inundation, while 1.8 m of rise could affect 13.1 million people—approximately two times larger than indicated by current populations.

The impacts of SLR on communities and the environment along the California coast are of substantial concern because the majority of the State’s population lives in coastal counties (approximately 87% of the population as of 2014) [9,10,11]. The California coastline, with a length of over 1800 km consists of wetlands, estuaries, and coastal resources, and is known as a major asset in the State [9,12]. The unprecedented SLR in modern human history may cause catastrophic impacts on this asset [4]. Marine inundation as a result of SLR will shift the coastline landward, erode beaches, accelerate cliff failure, degrade sensitive coastal habitats, and damage coastal infrastructure and humans [2,12,13,14,15]. SLR can also contribute to the degradation of coastal aquifers that have already been under pressure in California over the last century due to excessive water extraction and persistent, severe drought [16,17,18,19]. Saltwater intrusion into fresh water aquifers in California’s coastal regions has been attributed predominantly to ground water overdraft, especially in southern California [16,17,20,21]. However, recent studies demonstrate that SLR could also contribute to saltwater intrusion in coastal regions by raising the interface between intruding saltwater and overlying fresh water [14,22]. Therefore, SLR is expected to impact shallow coastal aquifers resulting in groundwater shoaling and emergence, which is considered another potential source of inundation [14,23,24].

Habel, Fletcher [25] described a modeling approach that simulates narrowing of the unsaturated space and groundwater inundation generated by SLR-induced lifting of coastal groundwater. The methodology combines terrain modeling, groundwater monitoring, estimation of tidal influence, and numerical groundwater-flow modeling to simulate future flood scenarios considering user-specified tide stages and magnitudes of SLR. They applied their method to the heavily urbanized and low-lying areas in Honolulu, Hawaii, USA. Results indicate that SLR of nearly 1 m generates groundwater inundation across 23% of the 13 km² study area, threatening $5 billion of real estate and 48 km of roadway. Befus, Barnard [26] assessed the increasing threat of coastal groundwater hazards from SLR in California. They modeled a range of SLR scenarios to assess the responses of water tables across the diverse topography and climates of the California coast. This research presented that with 1 m of SLR, flooded areas are predicted to expand 50–130 m inland, and low-lying coastal communities around San Francisco Bay are at a high risk. They illustrated that topography is a controlling factor; long-term rising water tables will intercept low-elevation drainage features, allowing for groundwater discharge that damps the extent of shoaling in 70% of California’s coastal water tables. Davtalab, Mirchi [27] also assessed the effect of sea level rise on groundwater inundation and its consequences on the capacity of retention ponds near Tampa Bay, Florida, USA. They found out that groundwater level rise can significantly reduce the soil storage capacity and cause infiltration at retention ponds, which reduces the ponds’ capacity in capturing upcoming storms. Partial groundwater inundation can affect retention ponds’ ability to reduce peak flow rates and keep the post-development outflow lower than or equal to pre-development conditions. Apart from causing SLR, climate change can shift the pattern and frequency of precipitation, causing unprecedented surface flooding and inundation [28,29,30]. In coastal regions, if rainfall over inland areas coincides with higher sea levels (e.g., high tide), impacts of flooding could be exacerbated due to the diminished ability of terrestrial drainage. In sum, emerging inundation sources in coastal regions under a changing climate are (Figure 2): (i) surface flooding as a result of precipitation, (ii) SLR-driven marine inundation, and (iii) SLR-driven groundwater rise.

Therefore, while a natural hazard is conventionally thought of as an infrequent (extreme) event with a short timescale, understanding the compound impacts of above-mentioned stressors occurring on longer timescales or with higher frequency is necessary for devising appropriate adaptation strategies to protect humans, infrastructure, and habitats against impacts of coastal climate change. Such impacts are undoubtedly understudied and have distinct knowledge gaps in identifying appropriate timescales in which these stressors are likely to impact the communities and the operation of infrastructure. Moftakhari, Salvadori [31] proposed a bivariate flood hazard assessment approach that accounts for compound flooding from river flow and coastal water level. They showed that a univariate approach may not appropriately characterize the flood hazard if there are compounding effects. Co-occurrence is of paramount importance as it may result in flood hazards with potential impacts larger than each extreme in isolation. Muñoz, Moftakhari [32] used a bivariate statistical analysis linked to hydrodynamic modeling, to quantify compounding effects of flood drivers and generate flood hazard maps near Savannah, Georgia, USA. Using statistical measures, they analyzed compounding effects of terrestrial/coastal flood drivers and wetland elevation correction on maximum floodwater height, and velocity for 50-year return period scenarios. Strauss, Ziemlinski [33] employed high-resolution edition of the National Elevation Dataset, and a tidal model covering the contiguous US, together with data from the 2010 Census, to quantify low-lying coastal land, housing and population relative to local mean high tide levels, which range from 0 to 3 m in elevation (North American Vertical Datum, a.k.a. NAVD of 1988). They used tidally adjusted approach, and estimated the contiguous US population living on land within 1 m of high tide to be 3.7 million. In 544 municipalities and 38 counties, they found that over 10% of the population lives below this line. Sukop, Rogers [34] modeled high temporal resolution of the impact of SLR, tides, and precipitation on water table flooding in the Arch Creek basin, Florida. Their study showed that rain-induced short-term groundwater table rise can be viewed as an important driver factor of flooding events under current conditions.

Numerical models are critically important tools for understanding flooding extent, hazard assessment and flood management planning. Available hydrologic and hydraulic models (H& H) of coastal watersheds can solve the continuity and flow motion equations in either 1 dimension along the stream line (1D), 2 dimensions along the land surface (2D), or in 3 dimensions (3D) [35]. Although 1D modeling approaches could be efficient and useful in some contexts, mainly for manmade channels, they face several limitations for overflow analysis [36]. When water begins to overflow, it propagates in 2 dimensions across the Earth’s surface, and there is a possibility that the flow return back to the main channel further downstream, which cannot be modeled in 1D platforms. Quirogaa, Kurea [35] have successfully applied the new HEC-RAS version 5 for flood hazardous assessment of the Mamore River flood and showed the 2D capability of the HEC-RAS model to generate flood depth, flow velocity and flood duration. Recent literature shows that 2D numerical models have been successfully applied for flood modeling [37,38,39]. Rangari et al. (2019), also used HEC-RAS 2D to assess the inundation risk in urban flood for different extreme rainfall scenarios [40].

The focus of the present study is to recognize the SLR-driven groundwater rise and assess its compound effect with the surface flooding caused by precipitation in the San Leandro watershed located in the Oakland Flatland which contains several designated “disadvantage community” blocks. A 2D model is developed for the study area using HEC RAS 5.0.7 to generate flood inundation maps. The paper is organized in four sections. The sections are as follows: methodology, results, conclusion, and discussion.

2. Methods

2.1. Study Area

The study area is San Leandro watershed (Figure 3) located in Oakland Flatlands, which is home to several disadvantaged community tracts according to DWR. The Flatlands are low-lying areas predominantly comprised of dredged sediment and construction debris adjacent to San Francisco Bay, Alameda Estuary, and San Leandro Bay [41]. Census data discloses a disproportionate concentration of poverty in the Flatlands, affecting a population that is majority African American, and Hispanic or Latino. Most of Oakland’s White population lives in the more prosperous foothills and hills neighborhoods [41]. There is literature available on the high probability of groundwater rise due to SLR in San Leandro downstream [23].

The San Leandro Creek watershed (Figure 3) covers 49.4 square miles in the counties of Alameda and Contra Costa. The Alameda County Flood Control and Water Conservation District reports that the creek itself is 22 miles long and is in a natural uncovered state [42]. The City of San Leandro is largely composed of sandstone, shale, serpentine outcrops (found underneath Lake Chabot), rhyolite, and conglomerate from the Cretaceous period [43]. The Oakland Flatlands generally sit on a base of Franciscan Formation, a sedimentary rock formation composed of sandstone, black shale, red chert, grey limestone, and greenstone, a volcanic rock [44]. San Leandro’s rich alluvial soils and temperate climate support a wide variety of plants and animals. Wetlands in the southwest part of San Leandro provide habitat for the salt marsh harvest mouse and other special status species. San Leandro Creek remains one of the few waterways in the urbanized East Bay that retains its natural character along most of its course. Oakland and San Leandro, CA, share the Mediterranean climate that characterizes the immediate San Francisco Bay Area [45]. Mild, wet winters and dry, cool summers are typical of the weather of both cities [45]. According to the San Francisco Bay Regional Water Quality Control Board, it is one of the few East Bay Creeks that is entirely above ground. The Water Board also reports that the annual rainfall in the watershed ranges from 18–28 inches per year [46]. A study performed by the San Francisco Estuary Institute for the years 1998–2007 estimated the average rainfall for the San Leandro Creek below Lake Chabot is 21 mm. According to that study, flooding information was as follows: minimum first flood volume was 64,000 m³, max first flood volume was 414,000 m³, and average first flood volume was 186,000 m³. The minimum largest flood volume was 230,000 m³, maximum largest flood volume was 907,000 m³, and the average largest flood volume was 490,000 m³. According to the USGS Stream Stats, the two year peak flow in the San Leandro watershed is around 1460 cfs.

Although energetic waves caused by coastal storms can be disastrous events in many coastal communities, in areas located around bays, e.g., San Francisco Bay, their impacts are diminished, and thus they are not included in the compound effect analysis presented herein. This claim has been proven by the USGS HERA inundation analyses, which represent the trivial impact of even a 100-year coastal storm on the inundation in San Leandro Bay [2].

In this study, a 3 m DEM (digital elevation model) is used as an input of the HEC-RAS 2D model. The DEM is available on the USDA website and is a part of the National Elevation Dataset (NED) [47]. NED is a seamless mosaic of best-available elevation data. The vertical datum is NAVD88. Other parameters required to develop the model are surface roughness, boundary conditions, and rainfall data. The surface roughness can be input to the model as Manning’s roughness coefficient. Manning’s roughness coefficient (n) for urban catchments can be evaluated from land use classification and the properties of the soil. In the study area of this paper, we have defined 6 different land uses, as river, roads, open areas, buildings, asphalt/concrete, and marshes, and calibrated the value of Manning’s roughness coefficient for each land use. Regarding the boundary condition, a typical precipitation event, i.e., 2-year peak flow (1460 cfs), is used as an upstream boundary condition to show expectable situations, as opposed to extreme events (e.g., 100- or 500-year events). The value of the 2-year peak flow was obtained from the USGS Stream Stats application. The effects of SLR and groundwater rise will occur on longer timescales as compared with acute storms. Therefore, SLR and groundwater rise have been modeled as the model’s boundary condition via three different scenarios:

• In the first scenario, the value of MHHW (mean higher high water) which is 1.94 m is used as downstream boundary condition [48];
• In the second scenario 1 m SLR is added to MHHW value; and
• In the third scenario groundwater inundation is also considered as a downstream boundary condition.

2.2. Modelling

A. Framework and Data Processing: Figure 4 represents the flowchart of the proposed methodology:

The core part of the methodology is constructing a representative 2D flooding model which has three steps as it can be seen in Figure 4. The first step is to processing the input files, including DEM data, roughness layers, and the geometry of the watershed. As mentioned in the previous section, we have used 3 m DEM, and 6 different roughness value for the defined domain of the downstream watershed. The second step is to develop groundwater inundation depth values which will be used as one of boundary conditions in the third step. For groundwater rise projection, we have obtained and processed data from literature as follows. In shallow, unconfined aquifers along the coast, SLR will increase recharge from seawater, increasing groundwater levels over longer time scales [14,23]. Where unconfined fresh groundwater is in contact with underlying seawater, the fresh groundwater floats on the higher-density seawater, and the average elevation of the water table will be above mean sea-level [49]. Near the coast, both the sea water and saline or overlying fresh groundwater respond to forcing, with the magnitude of the response diminishing towards the inland from the coast [14,50,51]. SLR will cause the groundwater table to rise in these areas, and in low-lying areas the water table could approach and ultimately rise above the ground surface [14,24]. This study uses the findings of a previous study of Plane, Hill [23] to identify potential groundwater flooding hotspots as sea levels rises in coastal cities. In their study, they used empirical depth-to-water data and a DEM of the San Francisco Bay Area to construct an interpolated surface of estimated minimum depth to water. Then, by subtracting the surface from the DEM, the groundwater inundation map is generated.

Finally, in the third step, inundation maps will be generated as the results of HEC-RAS output. As mentioned before, this study develops maps for three different downstream boundary conditions in order to assess the compound effect of mentioned factors on the inundation of low-lying coastal area. In the next paragraph, the detailed calculation of HEC-RAS model is presented. The outputs of this model are the illustrative inundation map that are represented in the Result section.

B. Modeling Concepts:

A distributed model was developed using a 2-dimensional, unsteady modeling of inundation under a 2-year precipitation event. Groundwater inundation and SLR heights were used as the boundary conditions. This modelling approach simulates the dynamics of compound inundation in a fine resolution, which is critical for devising proper flood management strategies. To validate the Stream Stats estimates, the IDF (intensity duration frequency) curves are analyzed, using the frequency of interest (2-year) and duration of interest (24-h, as the basis for flood analysis) to achieve the precipitation intensity. Then, precipitation intensity is translated to pick discharge rates using Modified Rational Method [52] for the watersheds of interest. Next, to enable unsteady modeling, runoff hydrographs are developed using the soils conservation service (SCS) triangular hydrograph method [53]. The following equations are used to define the shape of the triangular hydrograph.

$${\mathrm{T}}_{\mathrm{L}}=\frac{{\mathrm{l}}^{0.8}{\left(\frac{1000}{\mathrm{CN}}-9\right)}^{0.7}}{1900{\mathrm{Y}}^{0.5}}$$

(1)

$${\mathrm{T}}_{\mathrm{p}}= \frac{\Delta \mathrm{D}}{2}+{ \mathrm{T}}_{\mathrm{L}}$$

(2)

$${\mathrm{Q}}_{\mathrm{p}}= \frac{484\mathrm{A}}{{\mathrm{T}}_{\mathrm{p}}}$$

(3)

$${\mathrm{T}}_{\mathrm{b}}=1.67{\mathrm{T}}_{\mathrm{p}}$$

(4)

where T_L is lag time (h), l is flow length (ft), CN is curve number, Y is average watershed slope (%), T_p is time to peak (h), ΔD is duration of rainfall (h), Q_p is peak discharge (cfs), A is area (mi²), and T_b is time to base (h). Figure 5 represents the flow hydrograph of the catchment developed by applying SCS triangular hydrograph method. The duration of rainfall, and CN are considered as 24 h and 0.85, respectively. Excess rainfall is the rainfall that directly contributes to the surface runoff. It is calculated by subtracting the losses and/or abstractions from the total rainfall.

For simulating inundation under SLR, groundwater rise and extreme precipitation conditions, we have used the diffusion wave equation (Equation (5)) and unsteady differential form of the conservation of the mass equation (Equation (6)). The diffusion wave form of the momentum equation disregards the unsteady, advection, turbulence and Coriolis terms [35]:

$$-\mathrm{g}\nabla \mathrm{H}={\mathrm{c}}_{\mathrm{f}}\mathrm{V}$$

(5)

$$\frac{\partial \mathrm{H}}{\partial \mathrm{t}}+\nabla ·\mathrm{hV}+\mathrm{q}=0$$

(6)

where g is gravitational acceleration, $\nabla (\partial /\partial \mathrm{x},\partial /\partial \mathrm{y}$) is vector of the partial derivative in the x and y directions, H is water surface elevation, c_f is bottom friction coefficient, V (u, v) is velocity (where u and v are respectively the velocity components in the x and y directions), t is time, h is water depth, and q is a source or sink term. Using Manning’s equation to represent bottom friction, Equation (5) will turn into:

$$-\mathrm{g}\nabla \mathrm{H}=\left(\frac{{\mathrm{n}}^2\mathrm{g}\left|\mathrm{V}\right|}{{\mathrm{R}}^{\raisebox{1ex}{$4$}\!\left/ \!\raisebox{-1ex}{$3$}\right.}}\right)\mathrm{V}$$

(7)

where n is Manning’s roughness coefficient, and R is hydraulic radius. The velocity can be determined by assuming that flow movement is driven by the balance between the barotropic pressure gradient and bottom friction (8). Lastly, the diffusion wave is substituted into the continuity equation to derive a single equation model, the diffusion wave approximation, which is used to solve for H, the water surface elevation (WSE) (9):

$$\mathrm{V}=-\left({\left(\mathrm{R}\left(\mathrm{H}\right)\right)}^{\raisebox{1ex}{$2$}\!\left/ \!\raisebox{-1ex}{$3$}\right.}/\mathrm{n}\right)\frac{\left(\nabla \mathrm{H}\right)}{{\left|\nabla \mathrm{H}\right|}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}}$$

(8)

$$\frac{\partial \mathrm{H}}{\partial \mathrm{t}}-\nabla ·\frac{{\left(\mathrm{R}\left(\mathrm{H}\right)\right)}^{5/3}}{\mathrm{n}{\left|\nabla \mathrm{H}\right|}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}}\nabla \mathrm{H}+\mathrm{q}=0$$

(9)

For 2D calculations, the study area is discretized, encompassing the rivers and adjacent floodplain areas, in a collection of individual grid cells. Each grid cell contains elevation and roughness data to represent the terrain elevation and friction effects used to solve the diffusion wave equations. At each cell face, hydraulic properties (cross-sectional area, wetted perimeter, and hydraulic radius) are computed. Then WSE-volume relationships are computed for each cell and WSE-hydraulic properties relationships are computed for every cell face. Lastly, the diffuse wave equations were solved with an iterative scheme.

The unsteady flood inundation was simulated using the United States Army Corps of Engineer’s HEC–RAS 2D. Although the program can solve for either the 2D full Saint-Venant equations or the 2D diffusion wave equations, the 2D diffusion wave equations is preferred for analysis, because it provides around the same stable results as the Saint-Venant method for urban inundation problems, but in a much shorter computation time.

The 2D flow area in HEC RAS specifies the extent of area within which 2D flow calculations are performed. It is marked as a hatched polygon layer for the study area in Figure 6. A computational 2D mesh is established within the defined 2D flow area at a cell spacing of 10 m which generated 60,582 computational cells. The model delivers sufficient accuracy and stability with Courant numbers as high as 5.0 while using the diffusion wave equations. When running a model for a small urban area, Courant numbers less than 1.0 should be used for more output accuracy. Most importantly, the time step used must be sufficient to produce stable results which can be identified rather quickly. Generally, the time step or computation interval for a model simulation run should be small enough to allow water movement through computational cells [54]. The HEC RAS manual suggests 30 s or shorter as a reasonable time step for a 60 m cell spacing, though it may vary as long as reasonable results are produced without sacrificing accuracy. In this study, a time interval of 10 s is selected to produce the highest accuracy in the output without affecting model stability.

The upstream boundary condition is a unit hydrograph resulted from SCS triangular method explained earlier. For downstream boundary condition, three different scenarios are defined. In the first scenario, the value of MHHW which is 1.94 m is used as the downstream boundary condition [48]. In the second scenario, 1 m SLR is added to MHHW value, and in the third scenario groundwater inundation is also considered as a downstream boundary condition. Inundations caused by groundwater rise are modeled as boundary conditions for the model because it occurs over a much longer period of time. Inundation caused by sea-level rise in a high tide condition is also modeled as the downstream boundary condition for the model. Figure 6 represents the DEM of the watershed modeled in HEC-RAS. It is important to note that we have manipulated the DEM downloaded from USDA to show buildings on earth’s surface.

3. Results

As discussed earlier, there is a knowledge gap in understanding the compound effects of SLR-driven inundation and precipitation-based flooding using proper timescales. Flow motion in coastal watersheds during rainfall events can cause inundation in additional areas that a static analysis of compound effects is unable to identify. It is important to note that there has been an attempt in recent literature to understand the compound effects of flooding sources due to their catastrophic impacts on coastal communities, but those efforts are limited to extreme events, and do not include precipitation-based flooding, as our study did [32,55,56]. In previous literature, coastal storms are typically considered with SLR inundation without understanding the interconnections of SLR to groundwater rise [57].

Figure 7 illustrates the results for our study area. As mentioned earlier, San Leandro is home to several “disadvantaged community blocks”, according to CA Department of Water Resources (DWR). The figure, which is obtained using the specified methods in the previous section, show the inundated area under Scenario 1 shown in Figure 7a (high tide conditions and 2-year precipitation) and Scenario 2 shown in Figure 7b (high tide conditions, 1 m SLR, and 2-year precipitation) is not significant. The San Leandro creek is not expected to overflow and to affect the buildings under these scenarios. In other words, this projection shows that no part of the built infrastructure is expected to face inundation by such compound impacts. The inundated area in Scenario 2 increases by 9% compared to Scenario 1, i.e., when the SLR is taken into account.

In contrast, by considering groundwater inundation, a 1 m SLR, and a 2-year precipitation event “together” as the boundary condition (Scenario 3), a considerable amount of built infrastructure will be inundated (Figure 8). The inundated area in this scenario is 768 acres. As it is evident, elevated tailwater conditions caused by groundwater rise affects the capacity of the channel which causes overflow and will inundate the roads and residential areas, especially those in the downstream portion. Based on these results, we have found that groundwater rise attributed to SLR can both (i) significantly occupy the capacity of drainage infrastructure and (ii) cause elevated tailwater conditions resulting in stream overflow in upstream areas. The extent shown in Figure 8 is identified as a Disadvantaged Community Tract according to DWR.

In order to better represent the water elevation in different scenarios, two cross sections with significant inundation projections are selected to be analyzed. These two cross sections are identified in Figure 8. Figure 9 represents the side view of water elevation in these two cross sections. According to this figure, compound effects can cause a flooding depth of over 0.5 m in the built infrastructure, which may cause safety issues, especially for the transportation infrastructure.

Another variable that we analyzed here is the flow velocity. Our findings did not show a significant increase in the velocity of overland flow under the studied scenarios with the risk of impacting the safety of residents.

4. Discussion

There are studies available focusing on the compound effect of different driving factors of flooding in coastal areas. Santiago-Collazo, Bilskie [57] presented a comprehensive review of the most recent studies that quantify extreme flooding using variations of a compound inundation model. Their review represents that even most of the recent publications do not consider typical events and they are mostly focus on the “extreme hurricane and coastal surges”. Therefore, our method for considering the annual/typical precipitation coupled with other drivers like SLR and groundwater inundation is an innovative approach to assess the expectable inundation, especially in locations where coastal surges do not have paramount role in flooding. Xu, Xu [55] assessed the compound effect of precipitation and storm tides as flood drivers in coastal zones in Haikou, China. They figured out that the complex interplay between these two drivers can exacerbate the impacts of flooding. This study investigates the bivariate return period of compounding rainfall and storm tide events based on copula functions and the failure probability is used to assess the variation of bivariate flood risk during the entire project lifetime. Although, this study did not consider the effect of SLR and subsequent impacts, e.g., groundwater rise. Muñoz, Moftakhari [32] assessed the compound effect of high river discharge and extreme coastal water level in Savannah, Georgia, USA. They used a bivariate statistical analysis linked to hydrodynamic modeling to quantify compounding effects of flood drivers. Their results represent that compound flooding analysis in estuaries accounting for not only the likelihood of high river discharge and extreme coastal water levels but also the accuracy of maximum floodwater height and maximum floodwater velocity in coastal flood hazard maps. Furthermore, Liu, Cheng [56] developed a framework for exploring joint effects of conditional factors on compound floods in Texas, USA. These compound drivers consist of combinations of precipitation and surface runoff with the El Nino-Southern Oscillation and rising temperatures as underlying conditions. They mentioned that although the individual variable of precipitation and runoff may not itself be extreme, large exceedances can lead to flooding situations when combined. Sukop, Rogers [34], assessed the impact of SLR, rainfall and tide on groundwater inundation in Southeast Florida, USA. However, they did not work on the compound impact of rainfall, SLR, and groundwater inundation on flooding. By reviewing these recent articles that worked on compound flooding of coastal areas, it is evident that considering the compound effect of SLR, precipitation and groundwater inundation on coastal inundation is understudied. This paper proposed an applicable methodology which can address the flooding issues in coastal areas by providing 2D inundation projections, and can help decision makers to manage the flooding hazards in a practical manner.

5. Conclusions

In recent years, few studies addressed projecting the inundation caused by SLR-driven groundwater rise, and no study paid particular attention to the compound effects of these phenomena. The gap in literature addressing these particular concerns is substantial considering the important hydrologic and hydraulic effects. In the present study, a novel analytical framework was developed to understand the movements of the surface flow under a typical (as opposed to extreme) precipitation event considering the composite interaction with uprising groundwater and sea level rise. In the proposed approach, the dynamics of compound inundation is simulated in a fine resolution, which is critical for devising proper flood management strategies.

In the first step, a contour of the water table that was available in the literature was transformed to the model’s boundary conditions. The data required for running HEC-RAS 2D (including DEM, manning roughness coefficients, and boundary conditions) were gathered. In the next step, the HEC-RAS 2D model was developed for three different scenarios, (1) compound effect of 2-year precipitation and MHHW; (2) MHHW, 1m SLR, and 2-year precipitation; and (3) MHHW, 1m SLR, groundwater inundation, and 2-year precipitation. The results show that SLR itself is not expected to cause flooding problems in the built infrastructure. However, by incorporating groundwater rise, the inundation projection exceeded 768 acres affecting the urban areas including buildings and transportation infrastructure, especially those located in the downstream of the watershed. These results show the importance of considering SLR-driven groundwater rise when it coincides with precipitation flooding. The importance of combining these has not been addressed in previous studies. For these reasons, this paper aimed to propose an analytical framework to help decision makers understand the importance of these compound effects in coastal communities, specifically in disadvantaged communities.

The model can be applied to any region with minimum input data requirements such as precipitation events, the digital elevation model and catchment properties. The flood risk maps can then be developed by utilizing H&H models to help develop adaptation strategies. Additionally, the proposed model can be utilized to simulate real-time flood events to map inundated areas and to help local authorities decide on suitable adaptation strategies, such as implementing rainwater harvesting in the upper part of the watershed in order to reduce the amount of precipitation caused flooding.