Quantifying Aggravated Threats to Stormwater Management Ponds by Tropical Cyclone Storm Surge and Inundation under Climate Change Scenarios

: Stormwater management ponds (SMPs) protect coastal communities from ﬂooding caused by heavy rainfall and runoff. If the SMPs are submerged under seawater during a tropical cyclone (TC) and its storm surge, their function will be compromised. Under climate change scenarios, this threat is exacerbated by sea level rise (SLR) and more extreme tropical cyclones. This study quantiﬁes the impact of tropical cyclones and their storm surge and inundation on South Carolina SMPs under various SLR scenarios. A coupled hydrodynamic model calculates storm surge heights and their return periods using historical tropical cyclones. The surge decay coefﬁcient method is used to calculate inundation areas caused by different return period storm surges under various SLR scenarios. According to the ﬁndings, stormwater management ponds will be aggravated by sea level rise and extreme storm surge. In South Carolina, the number of SMPs at risk of being inundated by tides and storm surges increases almost linearly with SLR, by 10 SMPs for every inch of SLR for TC storm surges with all return periods. Long Bay, Charleston, and Beaufort were identiﬁed as high-risk coastal areas. The ﬁndings of this study indicate where current SMPs need to be redesigned and where more SMPs are required. The modeling and analysis system used in this study can be employed to evaluate the effects of SLR and other types of climate change on SMP facilities in other regions.


Introduction
Tropical cyclones (TCs) threaten coastal communities with flooding, including both freshwater flash floods caused by heavy rain and seawater inundation caused by storm surges. Stormwater management ponds (SMP), including both (wet) retention ponds which hold water year-round and (dry) detention ponds which only hold water for a short time before the water enters streams, are a critical type of infrastructure that protects coastal communities from flood hazards caused by heavy rainfalls and runoff. A study in southwestern Taiwan showed that detention ponds reduce the inundated potential, inundated volume, flood damage, and flood stage of peak flow [1]. Sahoo and Pekkat [2] found that detention ponds reduce the effect of urban flood in an urbanizing catchment, reducing flood depth by 46.5% and the inundated area by 43%. Similarly, Chrétien et al. [3] analyzed retention ponds in a small agricultural catchment of Canada and stated that retention ponds could reduce peak flows by 38%.
Retention ponds fail to function when the rainfall exceeds the pond's capacity, e.g., with the underprediction of rainfall intensity or unrealistic predictions of post-development We used the above modeled historical storm surge height data to estimate the storm surge heights for the South Carolina coast with 5-year, 10-year, 50-year, 100-year, and 1000-year return periods. The storm-surge-height return periods were calculated based on the classical Hazen method [16] 2 2 1 (1) where is the return period of a specific storm surge height, n is the total number of events (132 in this study), and m is the rank of a specific storm surge height. m = 1 when the maximum storm surge height is the 132-year high.
The resultant storm surge height return period is shown in Figure 2. The storm surge heights as a function of return periods fit into h(t) = −0.0394 + 1.38ln(t) − 0.107ln(t) 2 , where h is the storm surge height (in m) and t is the return period (in year). The storm surge heights for the return periods of 5, 10, 50, 100, and 1000 years are 1.9 m, 2.6 m, 3.7 m, 4.0 m, and 4.4 m, respectively. The 1000-year storm surge height has never occurred before and its value is extrapolated from past storm surge records. Note that these storm surge modeling results were based solely on changes in water level caused by TCs. Astronomical tide is another important factor affecting the actual storm surge height. We used the above modeled historical storm surge height data to estimate the storm surge heights for the South Carolina coast with 5-year, 10-year, 50-year, 100-year, and 1000-year return periods. The storm-surge-height return periods were calculated based on the classical Hazen method [16] T p = 2n 2m − 1 (1) where T p is the return period of a specific storm surge height, n is the total number of events (132 in this study), and m is the rank of a specific storm surge height. m = 1 when the maximum storm surge height is the 132-year high. The resultant storm surge height return period is shown in Figure 2. The storm surge heights as a function of return periods fit into h(t) = −0.0394 + 1.38ln(t) − 0.107ln(t) 2 , where h is the storm surge height (in m) and t is the return period (in year).
where is the return period of a specific storm surge height, n is the total number of events (132 in this study), and m is the rank of a specific storm surge height. m = 1 when the maximum storm surge height is the 132-year high.
The resultant storm surge height return period is shown in Figure 2. The storm surge heights as a function of return periods fit into h(t) = −0.0394 + 1.38ln(t) − 0.107ln(t) 2 , where h is the storm surge height (in m) and t is the return period (in year). The storm surge heights for the return periods of 5, 10, 50, 100, and 1000 years are 1.9 m, 2.6 m, 3.7 m, 4.0 m, and 4.4 m, respectively. The 1000-year storm surge height has never occurred before and its value is extrapolated from past storm surge records. Note that these storm surge modeling results were based solely on changes in water level caused by TCs. Astronomical tide is another important factor affecting the actual storm surge height. The storm surge heights for the return periods of 5, 10, 50, 100, and 1000 years are 1.9 m, 2.6 m, 3.7 m, 4.0 m, and 4.4 m, respectively. The 1000-year storm surge height has never occurred before and its value is extrapolated from past storm surge records. Note that these storm surge modeling results were based solely on changes in water level caused by TCs. Astronomical tide is another important factor affecting the actual storm surge height. For example, the once-in-100-years (return period = 100 years) storm surge height for the South Carolina coast is 4.0 m, which may coincide with a high tide of 1.5 m-the average high tide for the South Carolina coast-and become 5.5 m, or coincide with a low tide of −1.5 m and become 2.5 m, or coincide with a tide somewhere between the two extremes. Due to the uncertain possibilities of the tidal phase coinciding with storm surge height, the astronomical tide-caused water level changes are not included in the calculation of the storm-surge-height return periods. Instead, the impact of tides is considered as a worst-case scenario in the subsequent inundation modeling and the assessment of the risks of SMPs, in which case it is assumed that the storm surge height with a certain return period coincides with a high tide.

Inundation
The method for calculating inundation extent as a result of storm surge during a tropical cyclone is illustrated in Figure 3. While storm surge heights caused by past tropical cyclones were simulated using the coupled hydrodynamic model, inundation areas were calculated based on the geographic information system (GIS) technique using the QGIS platform. It is well established that storm surge height intrudes on land with decay [17]. The surge decay coefficient (SDC) is defined as follows: where H s is the storm surge height, H i is inundation elevation (the average elevation where the effect of the surge can reach), W i is inundation width (the longest distance of inundated area away from coast), and W c is the width of the constant surge from the coast. For example, the once-in-100-years (return period = 100 years) storm surge height for the South Carolina coast is 4.0 m, which may coincide with a high tide of 1.5 m-the average high tide for the South Carolina coast-and become 5.5 m, or coincide with a low tide of −1.5 m and become 2.5 m, or coincide with a tide somewhere between the two extremes. Due to the uncertain possibilities of the tidal phase coinciding with storm surge height, the astronomical tide-caused water level changes are not included in the calculation of the storm-surge-height return periods. Instead, the impact of tides is considered as a worstcase scenario in the subsequent inundation modeling and the assessment of the risks of SMPs, in which case it is assumed that the storm surge height with a certain return period coincides with a high tide.

Inundation
The method for calculating inundation extent as a result of storm surge during a tropical cyclone is illustrated in Figure 3. While storm surge heights caused by past tropical cyclones were simulated using the coupled hydrodynamic model, inundation areas were calculated based on the geographic information system (GIS) technique using the QGIS platform. It is well established that storm surge height intrudes on land with decay [17]. The surge decay coefficient (SDC) is defined as follows: (2) where is the storm surge height, is inundation elevation (the average elevation where the effect of the surge can reach), is inundation width (the longest distance of inundated area away from coast), and is the width of the constant surge from the coast. We assume that the inundation begins at the coastline so that = 0 m. We first run an idealized 2D hydrodynamic model simulation by setting sea level to the storm surge heights (the in Equation (2)) of 1.9 m, 2.6 m, 3.7 m, 4.0 m, and 4.4 m to determine the widths of the inundated area and the end inundation elevation . The widths of the inundated area and the end inundation elevation were obtained manually based on flooding extent maps produced by the hydrodynamic model Delft3D-FLOW.
Delft3D-FLOW is a 2D or 3D simulation program that calculates water's non-steady flow and transport caused by tidal and meteorological forcing on a rectilinear or a curvilinear fitted grid. The spatial resolution varies from 127 to 1268 m in the x direction and 55 to 557 m in the y direction. The time step is 0.4 s. Such grid spacing and time step allow stable integration in shallow water, e.g., when the water depth is 100 m, the Courant number is dt/(dx/sqrt(gH)) = 0.2, which is much smaller than required for computational stability under CFL conditions. The Holland model [18] is used to create the wind data based on hurricane track and intensity to drive Delft3D-FLOW. We assume that the inundation begins at the coastline so that W c = 0 m. We first run an idealized 2D hydrodynamic model simulation by setting sea level to the storm surge heights (the H s in Equation (2) Delft3D-FLOW is a 2D or 3D simulation program that calculates water's non-steady flow and transport caused by tidal and meteorological forcing on a rectilinear or a curvilinear fitted grid. The spatial resolution varies from 127 to 1268 m in the x direction and 55 to 557 m in the y direction. The time step is 0.4 s. Such grid spacing and time step allow stable integration in shallow water, e.g., when the water depth is 100 m, the Courant number is dt/(dx/sqrt(gH)) = 0.2, which is much smaller than required for computational stability under CFL conditions. The Holland model [18] is used to create the wind data based on hurricane track and intensity to drive Delft3D-FLOW.
The SDC values can then be determined based on the above equation. Additionally, later, these SDCs were applied to the inundation area calculation under different SLR scenarios using the geographic information system technique on the QGIS platform, which is more computationally efficient than numerically modeling each individual scenario. The QGIS Geospatial Data Abstraction Library (GDAL) was used to process the inundation maps. The steps for determining whether a giver inland location is inundated are as follows: 1.
Calculate the raster distance (RD) from the coastline to the inland location using the Proximity function. Note RD = 0 m at the coastline.

2.
Determine the water depth (WD) based on H s , SDC, RD, and DEM elevation (DEM) using the formula: WD > 0 indicates inundation. The Raster calculator function in QGIS is used in this step.
The storm surge height, inundation width, the inundation elevation, and surge decay coefficient for storm surges with different return periods are listed in Table 1.

Experiments-Impact of SLR to Stormwater Ponds
For each return period of storm surge height, the inundation is calculated under the following scenarios: (1) no tide and no SLR; (2) with tide and no SLR; (3) with tide and SLR = 0.5 m; (4) with tide and SLR = 1.0 m; (5) with tide and SLR = 1.5 m; and (6) with tide and SLR = 2.0 m. These SLR scenarios include the low-, intermediate-low-, intermediatehigh-, and high-risk scenarios proposed by Parris [12]. A total of 36 calculations (6 TC storm surge return periods each under 6 different SLR scenarios) were carried out. In each calculation, the number of submerged SMPs was counted. The impact of SLR scenarios on the number of SMPs affected during non-storm (return period of zero years) and storm conditions (return periods of 5, 10, 50, 100, and 1000 years) were assessed. The high-risk "hotspots" were identified where mitigation and adaptation efforts are required under future climate SLR scenarios.
Furthermore, the South Carolina coast was divided into three domains based on topography and the three SMP clusters centered in these three regions: Long Bay, Charleston, and Beaufort. In addition to the assessment of the entire South Carolina coast as a single domain, the impact of SLR on SMPs was also evaluated in detail for each of these three domains individually. Note that some SMPs are outside of these three regions.

Storm Water Pond and Topography Data
The South Carolina storm management pond survey [19] based on aerial images, which identified approximately 21,594 SMPs, was used in this study for the SMP locations ( Figure 4). The aerial images used to identify SMPs are provided by the National Agriculture Imagery Program (NAIP). The spatial resolution of these aerial images is about 1 m. The SMPs are concentrated in three high-density clusters in the Long Bay, Charleston, and Beaufort areas. These topography data are an important input for producing the inundation area. The topography data used in this study were obtained from the United States Geological Survey 3D Elevation Program (3DEP), with a spatial resolution of 30 m. inundation area. The topography data used in this study were obtained from the United States Geological Survey 3D Elevation Program (3DEP), with a spatial resolution of 30 m.  Figure 5 shows the number of SMPs in the Long Bay region that are vulnerable to TC storm surges with 0-year (non-storm), 5-year, 10-year, 50-year, 100-year, and 1000-year return periods under 0 m, 0.5 m, 1.0 m, 1.5 m, and 2.0 m SLR scenarios. A non-storm high tide could inundate 29 SMPs under current climate conditions. For stronger hurricanes, the number could rise sharply, reaching about 300 SMPs for a 100-year TC storm surge. However, SMPs at risk of 100-year and 1000-year storm surges appear to be similar. When the tide is not included, the number is significantly lower, emphasizing the importance of taking tides into account. The worst-case scenario method is used, which involves adding an average high tide to the storm surge heights for all return periods.  However, SMPs at risk of 100-year and 1000-year storm surges appear to be similar. When the tide is not included, the number is significantly lower, emphasizing the importance of taking tides into account. The worst-case scenario method is used, which involves adding an average high tide to the storm surge heights for all return periods.

Long Bay
Unsurprisingly, more SMPs would be at risk of TC storm surges under future SLR scenarios. According to Figure 5, under a non-storm condition, the number of affected SMPs would double from 29 to 54, with an SLR of 0.5 m. The number of affected SMPs continues to grow linearly reaching 157 with an SLR of 2 m. Under storm conditions, the number of affected SMPs also increases almost linearly with SLR. For the Long Bay area, each 0.5 m SLR would expose about 50 more SMPs to TC storm surge and inundation. Note that the number of affected SMPs listed above does not represent the total number of SMPs that will be inundated when a single hurricane hits the SC coast to cause storm surge and inundation at an arbitrary location. It is instead the number of SMPs that could be inundated by a local storm surge. Unsurprisingly, more SMPs would be at risk of TC storm surges under future SLR scenarios. According to Figure 5, under a non-storm condition, the number of affected SMPs would double from 29 to 54, with an SLR of 0.5 m. The number of affected SMPs continues to grow linearly reaching 157 with an SLR of 2 m. Under storm conditions, the number of affected SMPs also increases almost linearly with SLR. For the Long Bay area, each 0.5 m SLR would expose about 50 more SMPs to TC storm surge and inundation. Note that the number of affected SMPs listed above does not represent the total number of SMPs that will be inundated when a single hurricane hits the SC coast to cause storm surge and inundation at an arbitrary location. It is instead the number of SMPs that could be inundated by a local storm surge. Figure 6 depicts the inundation map for the Long Bay area caused by hurricanes with different return periods under an SLR = 2.0 m scenario. The green-yellow colors indicate topography above sea level, and the blue-red colors represent the water depth, which includes perennial wet areas (rivers and wetlands) as well as the flooded areas. More areas are inundated when the storm surge heights increase from non-storm to 5-, 10-, 50-, 100and 1000-year return periods, as indicated by the increase in the blue-red areas. Under non-storm conditions, SMPs in several locations along the coast of Long Bay immediately adjacent to the coastline are vulnerable to high tides. As the storm surge height increases, several high-risk hotspots become evident: the SMPs near Winyah Bay, particularly those near Georgetown City and Parleys Island, are at the greatest risk, as indicated by the clustering of the black-shaded areas. Note that Figure 6 has been rotated clockwise toward the northeast.  Figure 6 depicts the inundation map for the Long Bay area caused by hurricanes with different return periods under an SLR = 2.0 m scenario. The green-yellow colors indicate topography above sea level, and the blue-red colors represent the water depth, which includes perennial wet areas (rivers and wetlands) as well as the flooded areas. More areas are inundated when the storm surge heights increase from non-storm to 5-, 10-, 50-, 100and 1000-year return periods, as indicated by the increase in the blue-red areas. Under non-storm conditions, SMPs in several locations along the coast of Long Bay immediately adjacent to the coastline are vulnerable to high tides. As the storm surge height increases, several high-risk hotspots become evident: the SMPs near Winyah Bay, particularly those near Georgetown City and Parleys Island, are at the greatest risk, as indicated by the clustering of the black-shaded areas. Note that Figure 6 has been rotated clockwise toward the northeast. The number might increase sharply for stronger TCs, reaching more than 400 SMPs for a 100-year storm surge. However, the number of SMPs at risk of a 100-year storm surge and a 1000-year one appears to be similar. The worst-case scenario method is used for the Long Bay region, where an average high tide is added to the storm surge heights for all return periods. Overall, the SMPs in the Charleston area are more vulnerable to high tides and TC-caused storm surges than those in the Long Bay area, mainly due to the low-lying terrain of the Charleston Harbor area.  Figure 7 shows the number of SMPs in the Charleston region that are vulnerable to submersion from TC storm surge heights with 0-year (non-storm), 5-year, 10-year, 50year, 100-year, and 1000-year return periods under scenarios of 0 m, 0.5 m, 1.0 m, 1.5 m, and 2.0 m SLR. A non-storm high tide could threaten ~56 SMPs under current climate conditions. The number might increase sharply for stronger TCs, reaching more than 400 SMPs for a 100-year storm surge. However, the number of SMPs at risk of a 100-year storm surge and a 1000-year one appears to be similar. The worst-case scenario method is used for the Long Bay region, where an average high tide is added to the storm surge heights for all return periods. Overall, the SMPs in the Charleston area are more vulnerable to high tides and TC-caused storm surges than those in the Long Bay area, mainly due to the low-lying terrain of the Charleston Harbor area.  In Charleston Harbor, similar to Long Bay, more SMPs would be vulnerable to TC storm surges under future SLR scenarios. According to Figure 7, with an SLR of 0.5 m, the number of affected SMPs would double from 56 to 105 under non-storm high-tide conditions. The number would keep rising almost linearly, reaching 258 with an SLR of 2 m. The number of vulnerable SMPs also increases almost linearly with SLR during TC storm surges. In the Charleston Harbor area, there would be about 70 more SMPs at risk of storm surge and inundation for every 0.5 m SLR. Figure 8 depicts the inundation map for the Charleston Harbor area caused by TC storm surges under an SLR of 2.0 m. The green-yellow colors represent topography above sea level, and the blue-red colors represent the water depth, including perennially wet areas and flooded areas, i.e., areas that were previously dry but became wet during the  The green-yellow colors represent topography above sea level, and the blue-red colors represent the water depth, including perennially wet areas and flooded areas, i.e., areas that were previously dry but became wet during the TC storm surges. As shown by the increase in blue-red areas, when storm surge heights increase from non-storm to 5-, 10-, 50-, 100-, and 1000-year return periods, more areas and SMPs are inundated. The majority of SMPs in the Charleston area are located around the Charleston Harbor, and few are immediately adjacent to the coastline. SMPs in several locations around Charleston Harbor are at risk of non-storm high tides. As the TC storm surge height rises, one hotspot becomes evident: the SMPs near Mount Pleasant. Note that Figure 8 has also been rotated clockwise toward the northeast.  Figure 9 shows the number of SMPs in Beaufort County at risk of TC storm surges with 0-year (non-storm), 5-, 10-, 50-, 100-, and 1000-year return periods under 0 m, 0.5 m, 1.0 m, 1.5 m, and 2.0 m SLR scenarios. A non-storm high tide could inundate ~69 SMPs in the current climate. The number rises for stronger storm surges, reaching more than 700 SMPs for a 100-year storm surge. Notably, 100-year and 1000-year TC storm surges appear to expose similar numbers of SMPs. Again, the worst-case scenario method adds a high tide to the storm surge heights for all return periods. In the Beaufort area, there are significantly more vulnerable SMPs than in the Long Bay area and Charleston Harbor area, mainly due to its low-lying terrain.  Figure 9 shows the number of SMPs in Beaufort County at risk of TC storm surges with 0-year (non-storm), 5-, 10-, 50-, 100-, and 1000-year return periods under 0 m, 0.5 m, 1.0 m, 1.5 m, and 2.0 m SLR scenarios. A non-storm high tide could inundate~69 SMPs in the current climate. The number rises for stronger storm surges, reaching more than 700 SMPs for a 100-year storm surge. Notably, 100-year and 1000-year TC storm surges appear to expose similar numbers of SMPs. Again, the worst-case scenario method adds a high tide to the storm surge heights for all return periods. In the Beaufort area, there are significantly more vulnerable SMPs than in the Long Bay area and Charleston Harbor area, mainly due to its low-lying terrain. Similar to the other two regions, more SMPs would be at risk of storm surges under future SLR scenarios in Beaufort County. Figure 9 shows that an SLR of 0.5 m would double the number of SMPs that are affected by high tides from 69 to 131, and with an SLR of 2 m, the number would reach 357. Under TC storm surge conditions, the number of affected SMPs also increases almost linearly with SLR. For every 0.5 m SLR, 80 more SMPs in Beaufort would be at risk of TC storm surge and inundation. Figure 10 shows the inundation maps for Beaufort County area caused by TC storm surges under an SLR of 2.0 m. Unsurprisingly, when storm surge heights increase from non-storm to 5-, 10-, 50-, 100-, and 1000-year return periods, more areas are inundated. Beaufort County's SMPs, unlike those in Charleston Harbor, concentrate on the Hilton Head Island immediately adjacent to the coastline. The SMPs on Hilton Head Island along the coast are at risk during non-storm high tides, and they have higher risks as storm surge height rises: even a mild 5-year TC storm surge would submerge nearly half of Hilton Head Island's SMPs. Similar to the other two regions, more SMPs would be at risk of storm surges under future SLR scenarios in Beaufort County. Figure 9 shows that an SLR of 0.5 m would double the number of SMPs that are affected by high tides from 69 to 131, and with an SLR of 2 m, the number would reach 357. Under TC storm surge conditions, the number of affected SMPs also increases almost linearly with SLR. For every 0.5 m SLR, 80 more SMPs in Beaufort would be at risk of TC storm surge and inundation. Figure 10 shows the inundation maps for Beaufort County area caused by TC storm surges under an SLR of 2.0 m. Unsurprisingly, when storm surge heights increase from non-storm to 5-, 10-, 50-, 100-, and 1000-year return periods, more areas are inundated. Beaufort County's SMPs, unlike those in Charleston Harbor, concentrate on the Hilton Head Island immediately adjacent to the coastline. The SMPs on Hilton Head Island along the coast are at risk during non-storm high tides, and they have higher risks as storm surge height rises: even a mild 5-year TC storm surge would submerge nearly half of Hilton Head Island's SMPs. Figure 11 shows the number of South Carolina SMPs at risk of TC storm surges with various return periods under SLR scenarios. A non-storm high tide could inundatẽ 208 SMPs in the current climate. The number increases sharply for stronger hurricanes, reaching 1700 SMPs for a 100-year TC storm surge. The SMPs at risk of 100-and a 1000year TC storm surges appear similar. Under future SLR scenarios, more SMPs would be vulnerable to TC storm surges. Figure 11 shows that with an SLR of 0.5 m, the number of affected SMPs would double from 208 to 405 for non-storm high tides. The number continues to rise almost linearly, reaching 998 with a 2 m SLR. Under storm conditions, the number also increases linearly with SLR. For every 0.5 m SLR, 200 more South Carolina SMPs would be at risk from storm surge and inundation, or 10 more SMPs for every inch of SLR.  Figure 11 shows the number of South Carolina SMPs at risk of TC storm surges with various return periods under SLR scenarios. A non-storm high tide could inundate ~208 SMPs in the current climate. The number increases sharply for stronger hurricanes, reaching 1700 SMPs for a 100-year TC storm surge. The SMPs at risk of 100-and a 1000-year TC storm surges appear similar. Under future SLR scenarios, more SMPs would be vulnerable to TC storm surges. Figure 11 shows that with an SLR of 0.5 m, the number of affected SMPs would double from 208 to 405 for non-storm high tides. The number continues to rise almost linearly, reaching 998 with a 2 m SLR. Under storm conditions, the number also increases linearly with SLR. For every 0.5 m SLR, 200 more South Carolina SMPs would be at risk from storm surge and inundation, or 10 more SMPs for every inch of SLR.  Figure 12 shows SMPs at risk of TC storm surges with different return periods under an SLR = 2.0 m scenario. The danger zone hotspots identified in Figures 6, 8 and 10 are also evident in Figure 12, including Hilton Head Islands (Beaufort County), Mount Pleasant (Charleston area), and Pawleys Island (Long Bay). In addition to these hotspots, SMPs in Kiawah Island, between Charleston Harbor and Beaufort County, also have a high risk of being inundated by storm surges under SLR scenarios.  Figure 12 shows SMPs at risk of TC storm surges with different return periods under an SLR = 2.0 m scenario. The danger zone hotspots identified in Figures 6, 8 and 10 are also evident in Figure 12, including Hilton Head Islands (Beaufort County), Mount Pleasant (Charleston area), and Pawleys Island (Long Bay). In addition to these hotspots, SMPs in Kiawah Island, between Charleston Harbor and Beaufort County, also have a high risk of being inundated by storm surges under SLR scenarios.

Discussion
The result of this study inevitably contains uncertainties. The 30 m topography data, for example, may be too coarse for some SMPs to be resolved accurately. Furthermore, we used the high-tide water levels to represent the effect of astronomical tides during storm surges, which may overestimate the storm surge water level and the hence the damage caused by the storm surge to the SMPs. As a result of these uncertainties, the findings of

Discussion
The result of this study inevitably contains uncertainties. The 30 m topography data, for example, may be too coarse for some SMPs to be resolved accurately. Furthermore, we used the high-tide water levels to represent the effect of astronomical tides during storm surges, which may overestimate the storm surge water level and the hence the damage caused by the storm surge to the SMPs. As a result of these uncertainties, the findings of this study should be used with caution and interpreted as a general trend of aggravated threats to SMPs from climate change, such as sea level rise and more extreme tropical cyclones, rather than an accurate prediction of which storm pounds will be damaged in any given storm.

Conclusions
Tropical cyclone-induced storm surges cause inundation that could threaten the proper function of stormwater management ponds in the coastal South Carolina area, and the risks could be aggravated by future sea level rise.
A framework is developed to assess the impact of TC-induced storm surges and inundation on stormwater management ponds. A hydrodynamic model simulates storm surges caused by historical hurricanes that struck the coast of South Carolina between 1851 and 2018. Storm surge heights with various return periods are established. Inundation caused by tides and storm surges could threaten storm management ponds in coastal regions. The rise in sea level could exacerbate these risks. Under various SLR scenarios, the inundations caused by storm surge heights with return periods of 0 (non-storm high tides), 5, 10, 50, 100, and 1000 years are simulated, and the SMPs at risk of being affected by these events are counted. For the entire South Carolina region, the number of SMPs at risk of being inundated by tides and storm surges would increase almost linearly with SLR by around 10 SMPs for every inch of SLR, for TC storm surges with all return periods.
Three coastal regions-Long Bay, Charleston, and Beaufort-are analyzed in detail to identify the communities whose SMPs are at high risk in the face of future SLR scenarios. Hilton Head Islands, Charleston Harbor, Kiawah Island, George Town City, and Pawleys Island are identified as these high-risk hotspot areas. If submerged by storm surge and seawater inundation, current SMPs installed in South Carolina could be incapable of functioning to mitigate the increased surface runoff and flash flooding caused by extreme rainfalls. Therefore, the findings of this study provide necessary information regarding the locations where current SMPs must be re-designed and where more SMPs must be developed. The modeling and analysis system developed for this study can be applied to assess the impact of SLR and other climate change scenarios on SMP facilities in other regions.