1. Introduction
Sea level rise (SLR) is an ongoing climatic phenomenon caused by increased rates of thermal expansion, glacier mass loss, groundwater pumping and ice discharge from both ice-sheets [1] - [3] that continues to put coastal populations at risk for increased flooding after storms and during king tide events. This is of critical concern since nearly half the US population lives within 50 miles of the coast, involving most major commercial, leisure and import/export enterprises; hence much of the economic activity of the nation is tied to coastal communities. Coastal communities in the rest of the world have similar populations at risk.
Gregory et al. [3] note that within the last two decades, the global rate of sea-level rise has been larger than the 20th-century time-mean. Various researchers have already noted impacts on coastal and island environments [1] [4] - [13] , and most coastal communities are cognizant of the ongoing discussion about SLR, the associated loss of soil storage capacity and more intense storms overwhelming the current stormwater [12] [14] - [18] . Prior studies have shown that SLR impacts will be felt globally, but the Southeast Florida region, with its low-lying coasts, subtropical climate, porous subsurface formations, and hydrology, is one of the world’s most vulnerable areas [18] - [20] .
In southeast Florida, sea level rise (SLR) threatens both the 6.6 million people living and working in the built environment as well as the natural system that recharges the aquifer and drives much of the climate [12] [18] [21] . The region constitutes one-third of the state’s total population, which makes it essential to continue improving flood management strategies [8] [13] [22] - [24] . The region has the highest rates of projected population growth [21] , which creates the potential to place important infrastructure that we rely on every day at risk. SLR is expected to rise up to three feet by 2100 [18] [25] - [29] . Fortunately, SLR is a slow, steady creep, allowing communities to take the appropriate time to evaluate and plan for SLR to insure that stranded infrastructure and failure to construct projects are exceptions to the norm [30] .
While coastal populations are particularly at risk due to SLR inundation and storm surge, interior populations are also susceptible to rising water tables and extended periods of inundation caused by the inability to drain inland areas. Higher groundwater levels equate to reduced soil storage capacity, which means less capacity for soil to absorb precipitation, thereby increasing the risk of groundwater flooding [31] . Chang et al. [32] describe an overall “lifting process” by which there is a 1:1 ratio in water table elevation that correlated to sea-level rise. Confounding the global trend is that South Florida has distinct wet and dry seasons and is very flat topographically. Seventy percent of the annual precipitation typically falls from June to September, just before the king tides in late September and October. To address the seasonal flooding, flood protection efforts were constructed to drain the land and surficial aquifer quickly to permit development. Soil storage is limited because the aquifer levels are often just below the surface in the wet season, leading to flooding, necessitating the extensive drainage works facilities that discharge large volumes of water during the wet season. As a result there is a nexus that drives the public demands for future development in South Florida, while protecting infrastructure and property from flood damage. Such impetus requires making long-term decisions because infrastructure and development is not temporal―it is expected to last 50 years or more. Hence it is in the community’s interests to develop a stormwater planning framework to adapt to SLR and protect vulnerable infrastructure through a long-term plan. While uncertainties in the scale, timing and location of climate change impacts can make decision-making difficult, response strategies can be effective if planning is initiated early.
Prior work by the authors in Broward and Miami-Dade Counties indicated significant risk to existing and future populations due to the low lying land (50% of the urban area is under 5 ft. NAVD88-see [19] [30] [33] - [35] ). Palm Beach County is located just north of Broward and Miami-Dade Counties. While portions of Palm Beach County is coastal, the perception is that since it is north of the more seriously impacted areas of Miami and Fort Lauderdale and at higher elevations (15 ft. vs 5 ft.), there is less risk to the community, yet flooding remains an ongoing problem that appears to be increasing in frequency.
The City of West Palm Beach is a city of over 100,000 people in central Palm Beach County, Florida. The City is a major economic engine for the County, and hosts both governmental as well as commercial activities. The downtown area is located on the Intracoastal Waterway, separated from the Atlantic Ocean by the Town of Palm Beach. The City has experienced flooding in a number of neighborhoods, which appears to be increasing in frequency with time. As a result the City is in the process of developing a stormwater master plan (SWMP) that includes a component to evaluate and include the potential impacts of sea level rise (SLR) on long-term stormwater management in the City, something not previously done. The objectives of the SWMP were to develop a method for planning for SLR, and providing a means to prioritize improvements at the appropriate time. In addition, the goals were to provide guidance in developing a means to prioritize infrastructure to maximize benefit to the community.
However unlike Broward and Miami-Dade Counties, Palm Beach County was found to have far fewer monitoring wells, and generally less long-term data collection. Since long-term monitoring and monitoring well density are thought to be critical to developing mapping for analysis of sea level rise, the paucity of monitoring wells created a confounding issue at the start. To accomplish this task, the City, through its consultants, engaged the investigators to perform a Sea Level Rise Vulnerability Assessment using available groundwater and tidal data in ArcGIS with the following elements:
1) Establish baseline sea level and groundwater elevations and develop maps based on the best available information.
2) Develop topographic maps specific to climate change analysis.
3) Identify vulnerable areas.
4) Develop costs for stormwater infrastructure to address the impact of sea level rise.
5) Develop a toolbox of infrastructure mitigation and policy solutions.
2. Methodology
Southeast Florida is fortunate in that there is readily available LiDAR data for the coastal communities. The LiDAR data format used was the American Standard Code for Information Interchange (ASCII). This data format is easily handled by ArcGIS software. The ASCII format is comprised of the raw LAS LiDAR data type format, translated into a geographically referenced X, Y, Z global coordinate plane system [36] . Of the different topographical data repository sources the NOAA offered the data natively in ASCII format.
Prior work by Weiss [17] used a bathtub model for inundation. This leads to the needs to define the critical or “acceptable” level of service (frequency of flooding). The bathtub approach assumes that groundwater levels will be flat, and match the elevation of the ocean (typically the average) to determine vulnerability. This method is used by many governmental organizations due to the ease of data acquisition and model creation. The main disadvantages of this type of model is that is does not consider urban water control infrastructure such as dikes and canals that lead to overestimation of inundation [37] , nor the fact that groundwater levels cannot be “flat”, and generally correlate with high tide, not mean tide. The results can lead to underestimation of inundation because they do not identify low-lying inland areas that might flood at an earlier time than areas along the coast as a result of higher groundwater tables.
A modified bathtub model is a model that considers more than just static elevation to determine SLR vulnerability. The modified bathtub approach assumes that groundwater levels increase as one moves away from the coast, an assumption that is easily justified with groundwater monitoring data in many communities. The importance of the groundwater table in the model is that it is responsible for determining the soil storage capacity [38] . This method was previously explained in [18] [22] [31] [33] - [35] [39] [40] , among others. In these cases, the difference between the topographic surface and the groundwater surface was used to determine the potential vulnerability of infrastructure and land areas due to sea level rise in Miami-Dade and Broward Counties due to their high level of preconceived vulnerability due to having a low elevation ground surface. Projecting groundwater levels indicated a greater risk for flooding, and more rapid failure of roadway bases and buried infrastructure. As a result, water, sewer, stormwater and transportation infrastructure in low-lying areas may be compromised. Once vulnerable areas are identified, solutions and costs can be identified.
Developing the results for measuring vulnerability in the City required a series of tasks. The first was to create a topographic surface layer. The second task was an acceptable groundwater layer needed to be developed. In prior efforts, the number of monitoring wells was significant, but this was not the case for the City of West Palm Beach. The third was to find the difference between these layers in as Geographic Information System (GIS) called ARCMAP10.2 to define vulnerability. Finally sea level rise vulnerability under different scenarios could be created.
2.1. Surface Topography
The means to assess what is vulnerable to SLR requires detailed topographic information. Topography is a key parameter that influences many of the processes involved in coastal vulnerability, and thus, up-to-date, high-resolution, high-accuracy elevation data are required to model the coastal environment. Previous approaches to modeling inundation from simulated sea-level rise have been limited by coarse-resolution elevation datasets (surveys, field spot elevations, USGS maps) versus electronic imagery [11] [16] [41] - [44] . However communicating the importance of sea level rise to local entities requires better data [10] [45] [46] . Low resolution LiDAR is available in many areas, but the coarse vertical definition (±2 feet) is not useful for coastal areas where inches matter.
Duke et al. [44] and Romah [31] showed that while higher-resolution elevation data represent a significant advance for modeling sea-level rise impacts, there can be a large variability in inundation estimates depending on the horizontal fit of the raster data [10] . High-resolution elevation data resolves the topographic complexity on landscape processes, including drainage canals [13] [44] [47] . For this project, elevation data beyond the City limits is required for the analysis. The remaining LiDAR elevation outside of the City limits was obtained from the high resolution dataset shot by the City of Wellington (in 2015) and NOAA (in 2007), but calibrated by FAU (see Romah [31] ). An issue arises in LiDAR elevation data for the western portions of Palm Beach County since NOAA did not shoot high resolution LiDAR west of the Turnpike due to lesser population densities. The lower resolution LiDAR (±2 feet) data set was used and kriged to the high resolution set. The new Wellington high resolution LiDAR elevation data was incorporated into FAU’s database in January 2015.
2.2. Groundwater Data
After downloading and developing the LiDAR maps of the land surface into GIS layers, the second step for creating the model was incorporating a groundwater surface elevation. A series of data sources had to be mined and related to past events. These include data derived from the South Florida Water Management District’s data storage and retrieval system called DBHYDRO. The database includes historical information for monitoring wells, canal stages and permit data for modeling water supply wells. Priority was placed on monitoring locations with 15 or more years of data, and in particular, more recent data since tidal information indicates that king tides have been increasing over the last 15 years. The goal was to find common dates for high groundwater levels, defined as occurring in the upper 2 percent of all observations. These were compared to tidal data.
Groundwater surface elevations for monitoring wells downloaded from the SFWMD DBHYDRO database for the period from 1995 through a range of end to dates between May 2014 and December 2014 and then compared using Excel® spreadsheets. Older wells are less representative of current groundwater levels since groundwater appeared to be on an upward trend during peak seasons. However, Figure 1 shows that the number of available wells was limited, which could lead to large errors during the kriging process in ArcMap10.2. Only two lie in the City, so a search for additional data points was undertaken. The results were then tabulated in ascending order and reviewed
Figure 1. Monitoring wells with September 2012 values shown in NGVD29.
to determine common dates within the 98-100th percentile of highest elevations. Since 99% reliability is the level of service (LOS) discussed by the City (4 days/yr), four dates occurred regularly, but September 22, 2012 was the one common to the most data sets and that date was used to create wells (note it also occurs in the top 2% of tidal data).
A second dataset involved the City’s wellfield. In the late 1990’s the City had their wellfield modeled for water use permitting purposes by JLA Geosciences. This model was developed for a dry season-drought condition that matched a time period from 1996-1997, although the month was not given. Because there are a series of wellfields in and around the City, a static groundwater table cannot be created to properly mimic current conditions without taking into account these stressors. These drawdown areas were not picked up in the limited monitoring well data. As a result it was assumed that if the month of the modeled contours could be learned, this map might improve the krig (it was a different season and set of conditions-drought-than the 2012 data). As a result, the 1996-97 contours were used as a means to create a series of faux monitoring well levels using an ARCMAP10.2 as a means to develop a kriged groundwater surface along with assumed average tidal data. The intent of this map was to derive an overall surface for groundwater given pumping and other conditions. However there was uncertainty in using a model that was 20 years old, and lacking the background information on timing, calibration and tides for the period in question, while keeping in mind that the model was run for groundwater withdrawal purposes during extreme dry events when the aquifer was low.
In a fortunate circumstance, the South Florida Water Management District had a third set of data-MODFLOW files for a JLA model run for permitting purposes. The groundwater model was re-run by the authors and the output layers for each month were created. Because there are a series of wellfields in and around the City, the end of cycle results for the first layer was calibrated to December 1997. The actual tidal elevation for December 1997 was used as the eastern boundary layer for the map which is a kriged (smoothed) version of the output files.
Several issues arose at this point-3 different data sources, different datums, and different tidal conditions, and different seasons. The datums were all translated to NAVD88 using an algorithm developed by the South Florida Water Management District and made publically available. Then water level data for the common date was compared to the well levels measured in December 1997 and September 2012, and differences between the layers were developed. Consistently, groundwater elevations were higher along the coastal ridge (FEC railroad) and east, while remaining much the same on the west when compared. This is due to the water management within the Grassy Waters Preserve at basically the same elevation year-round. The result allowed for the development of a “kriged” wedge that reflects the difference in elevation between the December 1997 modeled water levels calibrated to 1997 groundwater levels and the groundwater levels measured on September 22, 2012.
2.3. Tidal Elevations
The purpose of SLR vulnerability modeling is to explore future vulnerabilities of infrastructure and property due to the increase in sea level by predicting how areas with low elevations may be affected by inundation from the ocean directly, from rising groundwater levels, and inundation from the inability of inland areas to drain. The Lake Worth pier data, corrected to the NAVD88 datum, was used for this data. The tide fluctuates during the year; during the spring months, the mean high tide is less than 1 foot nearly all of the time and periodically below 0.0. In the fall, mean high tide can be as high as 2 feet. If the current condition is developed as a baseline, the reference elevation needs to be defined as mean high tide or a higher percentile. It was determined that the 95 - 100 percentile tides occurred primarily in the September and October timeframe (by listing all tides in ascending order), which comports with the September 22 date for highest groundwater.
2.4. Comparing Topography with Groundwater for Measuring Vulnerability
Groundwater and surface topography are used collectively to help quantify soil storage capacities. Topography identifies vulnerable land, but soil storage capacity is what really identifies flooding potential. The difference between the land surface and the groundwater level is the soil storage amount. The US Army Corps of Engineers’ vulnerability definitions were used based on the difference between ground surface elevation and groundwater level elevation to define vulnerable areas with a difference of 0 or less, potentially vulnerable areas with a difference of 0 - 2 ft, and not vulnerable area with a difference greater than 2 feet. The term “potentially vulnerable” is used for areas that need further investigation to deal with the uncertainty of timing, surface improvements, etc. that might affect the situation. For ease of understanding, the protocol established by Romah [31] for SLR vulnerability was used: vulnerable areas are red, potentially vulnerable areas are yellow, and non-vulnerable areas are green. Results were developed for the current SLR scenarios of, 1 feet (0.3 m), 2 feet (0.6 m), and 3 feet (0.9 m).
2.5. Solutions
The final piece of the project was to identify costs and types of project to correct short and long-term deficiencies. This was accomplished by determining where current infrastructure might be lacking, identifying potential solutions for stormwater problems and assessing costs for same. The latter, due to the size of the City and lack of detailed stormwater modeling, were magnitude of scale estimates in 2015 dollars, based on bids reviewed or received by the investigators in southeast Florida.
3. Results
Figure 2(a) shows September 2012 monitoring well locations and the groundwater elevations at each point. Figure 2(b) is the kriged map of the data developed from the permit map. What this map shows is that compared with Figure 2(a), the groundwater table in coastal ridge is underestimated. Since that is an important part of the City’s
(a) (b) (c)
Figure 2. Comparison of the kirged water levels for the 2012 groundwater levels, map from permit files and groundwater model.
commercial area, the lack of data collection points would create an unrealistic (lower risk) view of the groundwater levels in these areas. They are higher than initially estimated. Figure 2(c) is the results for the MODLFOW model which extend the results of Figure 2(b) further. Figure 3 shows the delta surface between the September 2012 krig and the May 1996-97 krig which illustrates the change between the drought season and the wet season with the wells. Without the groundwater modeling data, the large drawdown areas would not be seen. These artificially keep the groundwater table down, and reduce potential flooding in these areas. This information is needed to compare the differences in the elevations between methods and identify how far off the krig might be with limited data as in the case of the 14 wells. Table 1 demonstrates that both the monitoring well and published data have nearly 2/3 of points within ±1 ft. of the actual model results despite noting that the wellfield drawdown points are underestimated in the 2012 monitoring well krig. Overall, while more data is probably better on a small scale, for this model the two data sets were not dramatically different in accuracy. As a result using the modeled results, increased to match the actual data in 2012, Figures 4(a)-(d) show the vulnerable areas of the City under the current sea level rise scenarios at the 99th percentile, 1ft. SLR, 2 ft. SLR and 3 ft. SLR.
Based on Figure 4, Table 2 outlines the area within the City subject to vulnerability
Table 1. Modflow vs kriged output from 14 monitoring wells or 1997 map.
Table 2. Summary of areas at risk and estimated capital needs.
Figure 3. Elevation change between 1997 model and 2012 high water when adjusted for tides and datum-September 2012 for use in assessing vulnerability between initial assessment and using groundwater model.
(red) and potential vulnerability (yellow) in the City. Note that the total area remains the same, but the red area increases as sea level rises. Based on current construction projects and their costs, an estimate of $15 million/200 acres was used to predict likely costs and timeframes (construction data reviewed by the authors). A much finer evaluation of neighborhoods is needed to determine if a given level of service is met or there are local circumstances that increase vulnerability.
The next step is to analyze vulnerability spatially, by overlaying development priorities with expected climate change on GIS maps to identify hotspots where adaptation activities should be focused. This effort includes identification of the critical data gaps which, when filled, will enable more precise identification of at-risk infrastructure and predictions of impacts on physical infrastructure and on communities. Based on findings of the vulnerable areas, the next step was the development of scenarios whereby toolbox options are utilized to address flooding in the community and where they might be applicable (see Table 3 and Table 4). The goal is identifying successful flood mitigation strategies used by other cities facing similar drainage and construction problems. Not all strategies are hard infrastructure (Table 3 only). Which strategies to pursue are locally driven in consultation with staff, based on identified vulnerabilities and cost effectiveness. These two issues are then combined to develop a framework to evaluate the impacts of climate change on infrastructure and economic development (as they are intrinsically intertwined).
The strengths of this framework are the initial focus on location-specific science, the use of both economic and social evaluation criteria, and the notion that the plan is not a fixed document, but rather a process that evolves in harmony with a changing environment. The final two steps occur at regular intervals by the community with associated adjustments made.
Whatever projection is placed on sea level rise dates, they need to consider uncertainty in the rate of warming, deglaciation, and other factors. When planning 50 - 100 years out, other factors can come into play as well. As a result, from the perspective of the authors, to allow flexibility in the analysis due to the range of increases within the different time periods, an approach that uses incremental increases of 1, 2, and 3 feet of SLR was considered for the scenarios.
4. Conclusions
The purpose of SLR vulnerability modeling is to explore future vulnerabilities of infrastructure, buildings and facilities on public and private property due to the increase in sea level by predicting how areas with low ground surface elevations may be affected by inundation from the rising ocean directly, from rising groundwater levels, and inundation from the inability of inland areas to drain.
The objectives of this paper were to: 1) establish baseline sea level and groundwater elevations and develop maps based on the best available information; 2) develop topographic maps specific to climate change analysis; 3) identify vulnerable areas; 4) develop costs for stormwater infrastructure to address the impact of sea level rise; and 5)
Table 3. Hard infrastructure improvements.
develop a toolbox of infrastructure mitigation and policy solutions. The results for the models that incorporate the loss of soil storage capacity created by rising groundwater levels indicate that the persistent flooding of roadways will start in the western portion of the study region farthest from the coast.
It was determined that in areas with limited monitoring well data, the potential to mis-estimate vulnerability is high. However regional groundwater models can be substituted for actual groundwater data if they are adjusted to the level of service condition. Likewise, good topographic information is needed to predict vulnerability from sea level and groundwater level rise. High resolution LiDAR is a must when inches matter. A final issue not considered in many prior SLR modeling efforts is the impact of soil storage capacity and groundwater levels on the potential to flood or damage infrastructure as they primarily focus on coastal regions and average mean tides (or mean high tides).
For the City of West Palm Beach a number of potential options are available to deal with sea level rise. Their vulnerability is similar to that of their southern neighbors in Broward and Miami-Dade County, so many options work in all three locations. The
Table 4. Soft infrastructure improvements.
perception that they were less vulnerable was found to be incorrect given that the coastal ridge appears to control groundwater levels in the City. The data revealed that over $400 million in current dollars might be needed to address stormwater issues arising from sea level rise before 2100.