Beyond the Tide: A Comprehensive Guide to Sea Level Rise Inundation Mapping using FOSS4G

1. Introduction

Climate change has become a severe issue on the planet, with an increase in the average global temperature leading to more frequent and severe weather events [1,2]. This rise in temperature causes ice sheets and glaciers to melt, contributing to sea level rise [3]. With record-breaking temperatures this year in southern states (e.g., Nevada, Arizona, Texas, and Florida [4,5]), coupled with the increasing trend of ocean temperature and unprecedented ice melting speed, an inevitable sea level rise is occurring. A few feet of sea level rise may not sound harmful, but it becomes lethal when combined with king tides and storm surge, especially considering the vast number of people living near the shoreline. According to Small et al. [6], approximately 400 million people live within 20 meters of sea level and within 20 km of the coastline worldwide. Sea level rise could result in a series of consequences, such as saltwater intrusion into inland rivers contaminating freshwater resources and harming agricultural lands [7,8]. Coastal ecosystems like mangroves, marshes, and coral reefs can be altered, resulting in the loss of habitats for certain species [9,10]. For human communities, shoreline infrastructure and properties are vulnerable to flooding, and homeowners may be forced to relocate, significantly disrupting the local economy [11]. Thus, to cope with the impact of sea level rise and prepare for future climate change, it is crucial to map the inundation in different sea level rise scenarios.

The remainder of the paper is structured as follows. Section 1.1 gives an overview of related literature, while Section 1.2 explains the novelty of our contribution. Section 2 describes the datasets and algorithms employed in the study. More specifically, Section 2.1 discussess elevation and water surface datasets and related vertical reference systems, while Section 2.2 explains the mapping process. Our FOSS4G (Free and Open-Source Software for Geospatial) implementation is described in Section 2.3, followed by the description of a case study in Section 2.4. Results of a case study demonstrating the approach are presented in Section 3, and finally Section 4 concludes the work and places it into context.

1.1. Previous work

Currently, various models have been used in sea level rise inundation modeling. Based on different spatial scales, complexity of computation, and data involved in the model, SLR models can be generally divided in to Hydrodynamic models and GIS-based approaches [12]. Hydrodynamic models, such as the SLOSH (Sea, Lake, and Overland Surges from Hurricanes) Model tailored to specific shoreline intricacies [13], Delft3D focusing on hydrodynamics and sediment transport [14], and ADCIRC (Advance Circulation Model for Coastal Ocean Hydrodynamics) modeling coastal water levels and currents, aim to replicate the dynamics of water inundation [15]. Although they are effective for predicting hurricane-driven storm surges, these models are constrained to specific areas due to their intricate demands for expertise, data, and computational power. In contrast, GIS-based models are emerging as a more accessible method to map SLR inundation. Requiring decent data quality and computational resources, they are also easily adaptable to various regions [16]. Medium-resolution GIS techniques sufficiently address the needs of scientific research, coastal management, and spatial planning in understanding coastal processes and assessing hazards [12]. Prominent GIS-based models include SLAMM (Sea Level Affecting Marshes Model), which relies on elevation data, BTELSS (Barataria-Terrebonne ecological landscape spatial simulation) addressing regional coastal dynamics [17,18], and the widely-used Bathtub model [19,20]. Williams and Lück-Vogel [12]’s enhanced version of the Bathtub model offers improved accuracy in storm surge contexts by considering additional parameters. Some GIS models, like DIVA, even delve into the broader impacts of SLR, including socioeconomic and biophysical implications, as well as the adaptive strategies’ costs and benefits [21,22]. Web applications, such as NOAA’s (National Oceonagraphic and Athmospheric Administration) Sea Level Rise Viewer1 and Climate Central’s Coastal Risk Screening Tool2 are just a few examples that popularize SLR inundation mapping by making SLR scenario maps available to the public.

While multiple examples of flood-related FOSS4G solutions exist ([23,24,25], most models and studies rely on proprietary software. For example, Lichter and Felsenstein [26] crafted an automated method utilizing ArcGIS geo-processing commands in Python to locally gauge the financial implications of SLR and intense flooding. Li et al. [27] and Perini et al. [28] both developed their coastal inundation models with the aid of ESRI’s Cost-Distance tool, while both models take into account hydrological connectivity to the coastline, and the latter also factors in the beach slope that sea-borne inundation water must traverse. Additionally, NOAA’s popular SLR Viewer tool employed ArcGIS software for SLR inundation mapping [29]. While the capabilities of proprietary software undeniably useful for modeling SLR scenarios, their proprietary nature limits broader adoption and utilization.

With the advancement of the open data policy issued by the Office of Management and Budget (OMB) in 2013, accessibility to geospatial data has significantly improved [30]. In the context of U.S. policy, the Federal Geographic Data Committee (FGDC), founded in 1990, aims to promote the coordinated development, utilization, sharing, and dissemination of geospatial data across all governmental levels, the private sector, non-profit entities, and academic institutions [31]. In 2001, the geospatial one-stop initiative introduced a portal for the discovery and retrieval of geospatial assets from various governmental bodies and it has been embraced by professionals, from urban planners and researchers to policymakers, all of whom depend on spatial data for informed decision-making [32]. The Geospatial Data Act of 2018 further emphasized the importance of public accessibility, endorsing data sharing across governmental tiers, all while heeding privacy and security concerns [33]. This act has underscored the principles of collaboration, transparency, and efficiency in geospatial data usage. With the progress made in open data policies, the development of Free and Open Source Software (FOSS) has also kept pace. The adoption of FOSS offers several advantages, including cost reduction, vendor independence, heightened customizability, assured interoperability, and robust community support. There’s been a notable momentum in the geospatial domain towards such open-source paradigms, namely Free and Open Source Software for Geospatial (FOSS4G). FOSS4G’s value, particularly in an economic context, has been highlighted in developing nations [34]. Addressing the issue of global coastal inundation, particularly in resource-limited settings, necessitates the creation of an open-source tool for sea level rise inundation mapping. Building on recent open data policies and advancements in free and open software, the groundwork is now laid for crafting a dedicated approach to mapping sea level rise inundation at local scales.

1.2. Contribution and significance of study

This research presents an open-source framework designed to map SLR inundation following NOAA’s passive coastal inundation mapping framework [35]. We leveraged various open-source GIS-based tools, including GDAL, QGIS, and GRASS GIS. Beyond the simple bathtub model, our methodology incorporates hydrological connectivity, but due to computing constraints, some influencing factors of SLR are omitted from the model (see limitations in Section 3.4). Harnessing a blend of high-resolution elevation models, tidal data, and projected SLR scenarios, we formulated an integrated workflow using FOSS4G tools to chart inundation extents across diverse SLR situations. Our approach highlight that FOSS4G tools present a viable and economic approach to map SLR inundation. This not only aids stakeholders in gauging coastal vulnerabilities but also helps in shaping resilient adaptation strategies. The adoption of an open-source strategy promotes widespread collaboration and ensures reproducibility, laying the foundation for a more exhaustive understanding of SLR effects on coastal systems.

2. Materials and Methods

To successfully map SLR inundation extents, a variety of data sources are needed throughout the process, which are detailed below.

2.1. Surface and Water Elevation Data

Digital elevation models (DEMs) are the representation of land elevation (topographic) or water depth (bathymetric). These datasets serve as the most valuable resources for coastal mapping applications [36]. The process of SLR inundation mapping deals with the interaction of land and water, and therefore relies on not only topographic and bathymetric elevation, but datasets representing water levels [37]. Since multiple elevation datasets have to be used in the process, the first challenge is maintaining consistency between them in terms of their spatial reference. While horizontal datum transformations are basic functionality in most GIS software, vertical datum transformations are less commonly used. However, to map coastal inundation, all datasets must reference elevations in the same system. These reference systems developed to measure measure heights and depths on Earth’s surface are called vertical datums [38]. They can be established using ellipsoidal, orthometric or tidal models. Figure 1a illustrates ellipsoidal model heights (h, N) that use a simplified, geometric surface to represent Earth’s shape and size, and orthometric model heights (H), which is a physical, gravimetric model that approximates mean sea level. In the United States, DEMs are created by various government agencies and other actors in varying detail. This includes local, state and federal organizations. An overview of this complex data ecoystem is provided by Juhász et al. [39]. While the nationally available 1 arc-second spatial resolution (30 m) National Elevation Dataset is not adequate for accurate coastal inundation mapping, high-resolution datasets are increasingly being created and made available publicly [40]. These provide enough detail that can facilitate coastal and marine spatial planning on localized scales.

Tidal datums reference water levels according to certain tidal stages [38]. While geodetic and orthometric datums are more precise and they can be used over a larger area, tidal levels vary by time and by location, therefore tidal datums always refer to one particular location. Figure 1b shows common tidal datums, while Table 1 lists the difference in elevation values of a specific tidal station in Virginia Key, FL. The Mean Higher-High Water (MHHW) is the average height of the highest tide recorded at a station over a period. It is more meaningful for example for land development purposes than lower tidal datums [41], and therefore it is a common reference datum to measure coastal inundation, since areas that are below the MHHW surface are likely to be flooded daily during a normal high tide. For other purposes, such as permitting, other tidal surfaces such as Mean High Water (MHW) are used. Table 1 shows that in Virginia Key, FL, the MHHW surface is 7 cm above the NAVD88 (North American Vertical Datum of 1988) orthometric datum, and 64 cm above the mean sea level at that location. In addition, oceans and seas are not flat surfaces and they deviate locally and regionally due to land movement, varying temperature and other factors [42,43]. This necessitates to integrate accurately representation of tidal surfaces in the inundation mapping process as opposed to using a constant elevation value in modeling inundation extent. A simple, commonly used approach to generate a digital representation of tidal surfaces is using data from tidal stations (see e.g. [44]), however, enough number of tidal gauges must be available to cover the study area. Another approach in the United States is to use VDatum, a database and software developed by NOAA, designed to vertically transform geospatial data among a variety of tidal, orthometric and ellipsoidal vertical datums [45]. However, VDatum covers only the continental United States, Puerto Rico, and the U.S. Virgin Islands, therefore, other areas need to use other approaches to general local tidal surfaces [46]. Software with similar functionality to VDAtum exist for other areas, such as the Australian Coastal Vertical Datum Transformation (AusCoastVDT) Software [47] and Vertical Datum Conversion Tool from the Canadian Geodetic Survey [48].

2.2. Bathtub Approach to Mapping SLR Inundation

Coastal areas become flooded as they submerge below surrounding water levels due to the rise in sea levels. This is typically simulated using a "bathtub" approach. This simple model assumes that areas that are below a water reference surface are filled with water like a bathtub. This approach relies on digital representations of the topography (DEM) and water surface, as described is Section 2.1. The process can be formalized as follows:

$$I_x{,}_y=\left\{\begin{array}{cc}1\hfill & \mathrm{if}{D}_x{,}_y\le W_x{,}_y\hfill \\ 0\hfill & \mathrm{if}{D}_x{,}_y>W_x{,}_y\hfill \end{array}\right.$$

(1)

where I is a binary variable indicating whether an area is inundated (1) or not inundated (0) at location $(x,y)$, while $D\mathrm{and}W$ are the DEM elevation value and projected water level respectively at that location. The bathtub can be "filled" with or without hydrological connectivity. In the latter case, only areas that are connected to the main water body (i.e. ocean or sea) are considered flooded. This can be achieved by multiplying I in Equation 1 by $C_x{,}_y$ which represents hydrological connectivity, i.e. whether a location is connected (1) or not connected (0) to the main body of water. Areas where $I_x{,}_y\ne (I_x{,}_y\times C_x{,}_y)$ can be considered low-lying.

Modeling SLR inundation in a GIS environment is a multi-step process that is usually computed using a combination of standard, cell-based raster operations. While most commonly used desktop GIS software do not have a dedicated tool to perform this task, the process is well-documented and has been described before (see e.g. [50]). A general flowchart is given in Figure 2. The process relies on three inputs, a raster DEM representing elevation values, a baseline raster surface representing the surface on which SLR is measured, and and a numerical SLR value that represents the rise in sea levels above that reference surface. MHHW is commonly used as a baseline surface (tidal surface), since it is a reliable reference point for understanding the upper limits of daily tidal fluctuations. We refer to SLR scenario as an x meter increase in sea levels above MHHW. The process starts by creating an increased water surface that is x m above the baseline surface (step 1). With reference to Equation 1, $W=t+x$, where t is the tidal surface. This new surface W represents the water surface in the bathtub modeling process to which the topography (DEM) is compared to. Next, the water surface is subtracted from the DEM (step 2). This resulting surface is negative if the water surface is below the DEM. In step 3, a cell-based classification is performed that classifies this surface into a binary map where cell values of 1 represent areas that fall below the water surface, in other words, "bathtubs" filled with water. Hydrological connectivity is evaluated in steps 4-6, where the binary map produced in step 3 is first recategorized by grouping cells that form physically discrete areas together. Then, categories that represent inundated areas connected with the main body of water are retained. This surface is identical to I in Equation 1. Subsequent steps are optional, but can be used for specific purposes. For example projected water depths can be generated by using I as a mask surface to overlay with the water depth surface generated in step 2. Additionally, a raster to vector conversion can be used to extract inundation polygons to use in vector-based spatial operations and visualization (step 8).

2.3. FOSS4G Implementation

While the implementation of the bathtub approach has been described before, it has been mainly done with proprietary tools, such as the ESRI ArcGIS suite of software [29,46,50]. In this section, we describe an alternative environment that utilizes FOSS4G software tools to implement bathtub modeling of SLR inundation.

2.3.1. Software and Tools

Our implementation uses three widely used FOSS4G software, namely the Geospatial Data Abstraction Library with the OGR Simple Features Library (GDAL/OGR), Geographic Resources Analysis Support System (GRASS GIS) and QGIS. These software frameworks represent an ensemble of geospatial tools and algorithms, each purpose-built to excel in specific geospatial tasks such as raster algebra, vectorization, data transformation, and more. It is the seamless integration and collaboration among these tools and algorithms that culminate in the creation of a comprehensive desktop GIS software solution.

In essence, GDAL/OGR3 provides abstract data models to represent raster and vector geospatial data. It also contains useful utility tools to translate between formats as well as to process data [51]. GDAL/OGR facilitates seamless integration into scientific workflows, enabling spatial data professionals to efficiently manage and extract meaningful insights from complex geospatial datasets.

GRASS GIS4 is a multi-purpose GIS software used for the production, analysis and presentation of geospatial data. It is one of the most comprehensive GIS software available, with its development work dating back to 1982 [52]. With an extensive array of tools for raster and vector data processing, terrain analysis, hydrological modeling, and more, GRASS GIS provides the means to conduct geospatial analyses, all within a flexible and customizable environment. Its integration of geostatistical, image analysis, and visualization capabilities renders it as an indispensable tool for spatial research.

QGIS5 is a desktop GIS software widely employed in scientific research due to its robust geospatial data handling capabilities and extensive analytical tools [53]. Offering a user-friendly interface, QGIS enables researchers to visualize, analyze, and manipulate spatial data. Its adaptability through plugins and support for numerous data formats further enhances its utility, facilitating complex geospatial analyses and fostering collaborative scientific exploration. In addition to its native algorithms, QGIS also provides an interface to run GDAL/OGR and GRASS GIS algorithms.

All three software are relased as open-source software and are official projects of the Open Source Geospatial Foundation (OSGeo6), which is an umbrella organization that fosters the development and promotion of open-source geospatial technologies and supports collaborative efforts within the geospatial community.

Table 2 relates the general workflow described in Figure 2 to specific software and algorithms. Specific implementation of these steps in the context of the case study presented in Section 2.4 and Section 3 are given in Appendix A and are also provided7.

2.3.2. Software Interfaces

These software and tools presented in the previous section can be interfaced in various ways. GDAL/OGR, QGIS and GRASS GIS all provide application programming interfaces (APIs) that allows to integrate their algorithms in custom-built software. However this can be considered advanced use-case, and for end uses, such as GIS analysts and coastal professionals, this is usually not necessary. More convenient ways to access functionality necessary for coastal inundation mapping are through graphical user interfaces (GUIs), such as desktop software environments, and through the command line.

2.4. Description of the Case Study

We demonstrate the implementation with a case study in Virginia Key, which is an approximately 3.5km² barrier island in Miami, FL, USA. The island constitutes of a mix of built-up and natural areas, and houses the Miami Seaquarium and a wastewater facility, as well as is home to large mangrove wetlands. A NOAA tidal station (#87232148) is located on southern part of the island. We adapt the extrapolated intermediate sea level scenario specific to this location for the year 2060. The amount of SLR is 0.44 m (1.44 ft). In the USA, this can be determined for each tidal station with sufficient observations using the Interagency Sea Level Rise Scenario Tool maintained by NASA9. Other scenarios and detailed methodological description of determining SLR are given by Sweet et al. [54]. The other main input is a LiDAR-based, bare-earth, hydro-enforced DEM provided by Miami-Dade County from 2021 [55]. Elevation values are referenced in meters in relation to NAVD88.

Since areas below the MHHW are susceptible to inundation, we use this surface to identify inundated areas. MHHW is determined at the Vriginia Key tidal station, and is 0.07m above the NAVD88 reference surface [49]. However, the ocean is not a flat surface, and there can be significant amount of variation in the height of tidal surfaces [56]. For this reason, we sampled MHHW at multiple locations and interpolate a continuous tidal surface using these distinct points. MHHW was also expressed in meters in reference to NAVD88. Chapter 3 demonstrates this process for Virginia Key. Step-by-step instructions are provided in Appendix A.

Figure 3. Satellite photo of the Virginia Key study area (red outline) (a); LiDAR-based DEM of the study area with NAVD88 elevations (b); Locations where MHHW was sampled (NAVD88) (c); and interpolated, continuous MHHW surface (d)

Figure 3. Satellite photo of the Virginia Key study area (red outline) (a); LiDAR-based DEM of the study area with NAVD88 elevations (b); Locations where MHHW was sampled (NAVD88) (c); and interpolated, continuous MHHW surface (d)

3. Results

The following section describes our FOSS4G approach for the Virginia Key Study areas. Throughout the section, references to specific steps in Figure 2 are also given.

3.1. Preparation of the Datasets

As described in Section 2.1, coastal inundation modeling requires high quality input elevation datasets that represent various surfaces. The DEM representing ground elevations was reprojected to the North American Datum 1983 (NAD83) / Universal Transverse Mercator (UTM) zone 17N coordinate system and resampled to 2.5 m pixel resolution. This resolution balances computational efficiency and results, and is considered adequate for local-scale modeling of inundation. The study area and the corresponding DEM is shown in Figure 1a-b. The difference in elevation within the study area is  17m between the lowest and highest points. The highest elevation area are located within the coastal dunes of the northern tip of the island.

While a common approach is using the MHHW surface value from the nearest tidal station as reference, we recognize that MHHW varies in space. We use NOAA’s VDatum software to sample MHHW at distinct locations shown in Figure 1c around the study area. VDatum is able to convert between various orthometric and tidal datums. We express MHHW in reference to NAVD88 so that the elevations can be directly compared with the DEM. Figure A3 shows the correct setup to transform MHHW to NAVD88 values in VDatum. It has to be noted that VDatum coverage area does not extend far beyond the shoreline. Unlike our study area, coastal mapping projects may aim to extend inland, which requires special steps to extend the tidal surface inland as well [29,57]. In the next step, we use the TIN (triangulated irregular network) approach to interpolate a continuous surface using the same spatial reference and pixel size as the DEM. TIN is an exact interpolation technique that retains the sampled values at those locations, and returns the interpolated elevation value based on linear interpolation between the sampled points within each triangle of the network. The interpolated surface is shown in Figure 1d. The difference between MHHW is 5cm between the Biscayne Bay and Atlantic sides of the study area. While this seems negligible, it already represents 11% of the projected SLR value in 2060. In addition, the difference over larger study areas can be more significant, which highlights the importance of using a more refined MHHW surface than using a constant value across the entire study area.

After preparing these datasets, all inputs of the modeling approach (see Figure 2 are available. The SLR value for the intermediate scenario in 2060 is 0.44m. The DEM and the MHHW surface are available as single-band rasters in the NAD83 / UTM Zone 17N coordinate system with 2.5 m pixel resolution, pixel values representing MHHW elevations above NAVD88 orthometric datum.

3.2. Simulating SLR Inundation

The SLR scenario value (0.44 m) can be added to the interpolated MHHW surface using cell-based raster algebra. The new water surface represents MHHW in 2060 according to the intermediate SLR scenario (step 1). Equation 1 is implemented in two steps, first by subtracting the new water surface (W) from the DEM (D) (step 2), then evaluating pixel values to determine inundation extent (step 3). Figure 4a shows inundation depth after step 2. Since areas where $D>W$ cannot be considered inundated, these values are retained as 0. In step 3, a binary map is created where negative values are reclassified and given the value of 1, representing areas that are below the SLR scenario (Figure 4b). The previous steps are executed using raster algebra which is usually implemented in both FOSS and GIS software. GDAL provides a command-line interface through the gdal_calc python script10. This is also accessible from the QGIS GUI in addition to QGIS’ own implementation.

Next, hydrological connectivity is evaluated. Pixels in the the binary water extent map generated in step 3 (Figure 4b) take the value of 1 if the pixel is below the SLR scenario (hence potentially inundated) or 0 if it lies above the SLR scenario (i.e. not inundated). Neighboring inundated cells can be grouped together to form physically distinct areas with the r.clump GRASS GIS algorithm (step 4). This tool finds all contiguous areas of the same cell value and assigns a unique category to these areas (also referred to as "clumps"). Figure 4c illustrates clumps within the study area. Unique categories are distinguished with different colors. In this figure, the category representing the main body of water is shown in purple with the category ID 538. Identifying this category is a manually process and can be done in the QGUS GUI through the Identify tool (step 5).

Once the category number of the water body is identified in the clumped map, this area can be used to extract hydrologically connected areas (step 6). This is done by reclassifying the clumped raster map into a binary map using gdal_calc. The result map is shown in Figure 4d. The difference between this output and the previously created map in Figure 4b is that inland areas below the SLR scenario are removed.

In step 7, the connected areas raster layer can be used as a mask to extract low-lying areas, that are below the SLR scenario but not hydrologically connected, as well as to exclude these areas from the inundation depth raster created earlier in step 2 (Figure 4a).

In step 8, raster-vector conversion is performed on the binary inundation map from step 3 and the connected areas map from step 6. This conversion is done using the r.to.vect function available in GRASS GIS. The conversion process transforms the low-lying areas and the hydrologically connected areas from binary rasters into vector polygons. These vector datasets are more suitable for certain spatial analyses and can be integrated with other vector-based datasets for a comprehensive assessment. Other advantages of the vector format in this context is the potentially smaller size of the dataset over larger areas as well as vector-based visualization techniques. Visualization options are explored in Section 3.3 below. its

3.3. Visualization of SLR and Low-Lying Areas

Figure 5a-b show the study area with and without displaying coastal inundation under the 2060 intermediate SLR scenario (+0.44 m). SLR inundation is usually visualized using pre-rendered outputs similar to Figure 4a and Figure 4e. There are two main approaches. First, similarly to Figure 4a and Figure 5c, inundation depth can be displayed with graduated colors. This approach is popularized by NOAA’s Sea Level Rise Viewer11 [58] that can display SLR for the contiguous United States. Furthermore, NOAA also allows to integrate these visuals in other application through serving the data through application programming interfaces (APIs), OGC (Open Geospatial Consortium) compliant services and ESRI-compatible maps. This allows other users to ingest SLR visualizations within their own GIS system [58]. A similar approach shown in Figure 5d displays pre-rendered, cell-based raster images using a single color to represent inundated areas. Low-lying areas can be displayed with both of these approaches.

As described in Section 3.2, inundated and low-lying areas can also be converted into vector polygons, which opens up the flexibility of visualizing vector datasets. An example is given in Figure 5e, where in addition to the semi-transparent inundation polygon, low-lying areas are symbolized with dashed polygons. Outside of desktop GIS environments, the flexibility of vector visualization can be exploited in interactive web applications that show the impact of SLR [39].

3.4. Limitations of the Bathtub Mapping Approach

While the passive inundation modeling approach provides a valuable tool for understanding potential flood risks in various scenarios, it has its inherent limitations. These are given below:

• The approach relies on static input data, however, over time, natural and artificial processes can significantly alter the landscape, potentially leading to inaccurate flood predictions.
• The approach doesn’t consider the dynamic interplay of water flow, wave action, and wind. This can lead to an oversimplification of flood scenarios, especially in areas prone to storm surges or rapidly changing water levels.
• We assume universal rise in water levels across the entire study domain, however factors such as tidal variations, river discharges, and localized rainfall can cause significant disparities in water level changes across a region. This is more prominent over large study areas
• The accuracy and resolution of the data inputs, such as DEMs directly impact the reliability of passive inundation models

It is essential for researchers and policymakers to recognize these limitations when interpreting results from passive inundation models. While they offer valuable insights, they should ideally be used in conjunction with other modeling techniques and real-world data for a comprehensive understanding of flood risks.

4. Discussion and Conclusions

The escalating threat of sea level rise due to climate change underscores the need for accurate and accessible tools to predict and visualize potential inundation scenarios. Historically, most studies and applications in this domain have relied heavily on proprietary tools, which often come with significant costs and limited accessibility. This creates barriers, especially for stakeholders in developing nations where resources are constrained, and the need for such tools is pressing.

Our study bridges this gap by presenting a comprehensive methodology for mapping SLR inundation using the bathtub approach, but with a distinct emphasis on leveraging the capabilities of free and open-source tools. This democratizes the inundation mapping process, making it accessible to a broader range of users, from researchers in resource-limited settings to local governments in developing nations. The case study of Virginia Key not only demonstrated the practical application of this approach but also emphasized the importance of accounting for spatial variations in tidal surfaces.

While the results derived from the bathtub model provide valuable insights, they are inherently constrained by the static nature of the input data and the absence of dynamic hydrological considerations. However, the adaptability and cost-effectiveness of FOSS4G tools ensure that the generated outputs can be tailored to diverse audiences, from policymakers to local communities.

In conclusion, our work offers a significant contribution by providing an alternative to proprietary tools for SLR inundation mapping. By harnessing the power of FOSS4G tools, we present a solution that is both accessible and customizable, catering especially to those in resource-limited settings. As the global community confronts the challenges of climate change, such tools will be pivotal in driving informed decision-making and fostering proactive adaptation strategies.

Appendix A.1. Step 1

The command above adds 0.44 (amount of SLR scenario in meters) to the input raster tidal_surface_navd88.tif using cell-based raster algebra. The -A flag specifies the input dataset, the --calc parameter contains the raster algebra expression, in this case a simple addition operation. The -co COMPRESS=DEFLATE flag applies the DEFLATE losless compression method to the output GeoTIFF. Figure A1 illustrates how to run the same tool within the QGIS GUI.

Figure A1. Running gdal_calc in QGIS

Figure A1. Running gdal_calc in QGIS

Appendix A.2. Step 2

The command above uses the DEM (project_dem_metric.tif) and the output of step 1 (water_surface.tif) labeled as A and B respectively. The expression evaluates whether the DEM <= MHHW increase with the SLR scenario, and if so, subtracts these values resulting in the water depth. All other values are retained as 0s since the logical multiplier A <= B evaluates to 0.

Appendix A.3. Step 3

Instead of retaining water depth, this command generates a binary water extent raster.

Appendix A.4. Step 4

Grouping connected cells is done by running the r.clump GRASS GIS algorithm. Although GRASS GIS can be run natively, this command above uses the QGIS Process Executor engine. The advantage of this approach is that data conversion from GeoTIFF to GRASS’ native format and other settings are handled automatically.

Appendix A.5. Step 5

The water surface can be visually identified in QGIS using the Identify tool (Figure A2. In the case study, this corresponds to the value of 538.

Figure A2. Identifying the category representing the water extent in QGIS.

Figure A2. Identifying the category representing the water extent in QGIS.

Appendix A.6. Step 6

The command above uses the category identified in the previous step. It retains pixels with the value of 1 if they belong to the hydrologically connected main body of water. All other pixels are classified as 0.

Appendix A.7. Step 7

The command above retains low-lying areas with the pixel value of 1 by comparing the hydrologically connected output from the previous step with areas that are below the SLR scenario (step 3).

Appendix A.8. Step 8

In a 2-step process, first the GRASS GIS algorithm r.to.vect is used to convert the hydrologically connected raster dataset into an ESRI shapefile. Then, ogr2ogr can be used to retain only the polygons that were created with inundated pixels.